Numerical Simulation of Detonation in Condensed Phase Explosives

DTIC Science & Technology

1998-08-01

34Numerical modelling of shocks in solids with elastic-plastic conditions", Shock Waves, 3: 55-66. 22. Jones, D.A., Oran, E.S. and Guirguis , R. (1990). "A...China Lake, CA 93555-6001, preprint. 55. P.J. Miller , P.J. and G.T. Sutherland, G.T. (1996) Reaction Rate Modelling of PBXN- 110, Shock Compression...report describes the development of a two-dimensional multi-material Eulerian hydrocode to model the effects of detonating condensed phase explosives on

Parametric study of graphite foam fins and application in heat exchangers

NASA Astrophysics Data System (ADS)

Collins, Michael

This thesis focuses on the simulation and experimental studies of finned graphite foam extended surfaces to test their heat transfer characteristics and potential applications in condensers. Different fin designs were developed to conduct a parametric study on the thermal effectiveness with respect to thickness, spacing and fin offset angle. Each fin design was computationally simulated to estimate the heat transfer under specific conditions. The simulations showed that this optimal fin configuration could conduct more than 297% the amount of thermal energy as compared to straight aluminum fins. Graphite foam fins were then implemented into a simulation of the condenser system. The condenser was simulated with six different orientations of baffles to examine the incoming vapor and resulting two-phase flow patterns. The simulations showed that using both horizontal and vertical baffling provided the configuration with the highest heat transfer and minimized the bypass regions where the vapor would circumvent the graphite foam. This baffle configuration increased the amount of vapor flow through the inner graphite fins and cold water pipes, which gave this configuration the highest heat transfer. The results from experimental tests using the condenser system confirmed that using three baffles will increase performance consistent with the simulation results. The experimental data showed that the condenser using graphite foam had five times the heat transfer compared to the condenser using only aluminum fins. Incorporating baffles into the condenser using graphite foam enabled this system to conduct nearly ten times more heat transfer than the condenser system which only had aluminum fins without baffles. The results from this research indicate that graphite foam is a far superior material heat transfer enhancement material for heat transfer compared to aluminum used as an extended surface. The longitudinal and horizontal baffles incorporated into the condenser system greatly enhanced the heat transfer because of the increased interaction with the porous graphite foam fins.

The calculation of the phase equilibrium of the multicomponent hydrocarbon systems

NASA Astrophysics Data System (ADS)

Molchanov, D. A.

2018-01-01

Hydrocarbon mixtures filtration process simulation development has resulted in use of cubic equations of state of the van der Waals type to describe the thermodynamic properties of natural fluids under real thermobaric conditions. Binary hydrocarbon systems allow to simulate the fluids of different types of reservoirs qualitatively, what makes it possible to carry out the experimental study of their filtration features. Exploitation of gas-condensate reservoirs shows the possibility of existence of various two-phase filtration regimes, including self-oscillatory one, which occurs under certain values of mixture composition, temperature and pressure drop. Plotting of the phase diagram of the model mixture is required to determine these values. A software package to calculate the vapor-liquid equilibrium of binary systems using cubic equation of state of the van der Waals type has been created. Phase diagrams of gas-condensate model mixtures have been calculated.

Reactive simulation of the chemistry behind the condensed-phase ignition of RDX from hot spots.

PubMed

Joshi, Kaushik L; Chaudhuri, Santanu

2015-07-28

Chemical events that lead to thermal initiation and spontaneous ignition of the high-pressure phase of RDX are presented using reactive molecular dynamics simulations. In order to initiate the chemistry behind thermal ignition, approximately 5% of RDX crystal is subjected to a constant temperature thermal pulse for various time durations to create a hot spot. After application of the thermal pulse, the ensuing chemical evolution of the system is monitored using reactive molecular dynamics under adiabatic conditions. Thermal pulses lasting longer than certain time durations lead to the spontaneous ignition of RDX after an incubation period. For cases where the ignition is observed, the incubation period is dominated by intermolecular and intramolecular hydrogen transfer reactions. Contrary to the widely accepted unimolecular models of initiation chemistry, N-N bond dissociations that produce NO2 species are suppressed in the condensed phase. The gradual temperature and pressure increase in the incubation period is accompanied by the accumulation of short-lived, heavier polyradicals. The polyradicals contain intact triazine rings from the RDX molecules. At certain temperatures and pressures, the polyradicals undergo ring-opening reactions, which fuel a series of rapid exothermic chemical reactions leading to a thermal runaway regime with stable gas-products such as N2, H2O and CO2. The evolution of the RDX crystal throughout the thermal initiation, incubation and thermal runaway phases observed in the reactive simulations contains a rich diversity of condensed-phase chemistry of nitramines under high-temperature/pressure conditions.

CFD Analysis of Evaporation-Condensation Phenomenon In an Evaporation Chamber of Natural Vacuum Solar Desalination

NASA Astrophysics Data System (ADS)

Ambarita, H.; Ronowikarto, A. D.; Siregar, R. E. T.; Setyawan, E. Y.

2018-01-01

Desalination technologies is one of solutions for water scarcity. With using renewable energy, like solar energy, wind energy, and geothermal energy, expected will reduce the energy demand. This required study on the modeling and transport parameters determination of natural vacuum solar desalination by using computational fluid dynamics (CFD) method to simulate the model. A three-dimensional case, two-phase model was developed for evaporation-condensation phenomenon in natural vacuum solar desalination. The CFD simulation results were compared with the avalaible experimental data. The simulation results shows inthat there is a phenomenon of evaporation-condensation in an evaporation chamber. From the simulation, the fresh water productivity is 2.21 litre, and from the experimental is 2.1 litre. This study shows there’s an error of magnitude 0.4%. The CFD results also show that, vacuum pressure will degrade the saturation temperature of sea water.

Imaging photonic crystals using hemispherical digital condensers and phase-recovery techniques.

PubMed

Alotaibi, Maged; Skinner-Ramos, Sueli; Farooq, Hira; Alharbi, Nouf; Alghasham, Hawra; de Peralta, Luis Grave

2018-05-10

We describe experiments where Fourier ptychographic microscopy (FPM) and dual-space microscopy (DSM) are implemented for imaging photonic crystals using a hemispherical digital condenser (HDC). Phase-recovery imaging simulations show that both techniques should be able to image photonic crystals with a period below the Rayleigh resolution limit. However, after processing the experimental images using both phase-recovery algorithms, we found that DSM can, but FPM cannot, image periodic structures with a period below the diffraction limit. We studied the origin of this apparent contradiction between simulations and experiments, and we concluded that the occurrence of unwanted reflections in the HDC is the source of the apparent failure of FPM. We thereafter solved the problem of reflections by using a single-directional illumination source and showed that FPM can image photonic crystals with a period below the Rayleigh resolution limit.

Pressurization of a Flightweight, Liquid Hydrogen Tank: Evaporation & Condensation at a Liquid/Vapor Interface

NASA Technical Reports Server (NTRS)

Stewart, Mark E. M.

2017-01-01

This paper presents an analysis and simulation of evaporation and condensation at a motionless liquid/vapor interface. A 1-D model equation, emphasizing heat and mass transfer at the interface, is solved in two ways, and incorporated into a subgrid interface model within a CFD simulation. Simulation predictions are compared with experimental data from the CPST Engineering Design Unit tank, a cryogenic fluid management test tank in 1-g. The numerical challenge here is the physics of the liquid/vapor interface; pressurizing the ullage heats it by several degrees, and sets up an interfacial temperature gradient that transfers heat to the liquid phase-the rate limiting step of condensation is heat conducted through the liquid and vapor. This physics occurs in thin thermal layers O(1 mm) on either side of the interface which is resolved by the subgrid interface model. An accommodation coefficient of 1.0 is used in the simulations which is consistent with theory and measurements. This model is predictive of evaporation/condensation rates, that is, there is no parameter tuning.

Condensation and dissociation rates for gas phase metal clusters from molecular dynamics trajectory calculations

DOE PAGES

Yang, Huan; Goudeli, Eirini; Hogan, Christopher J.

2018-04-24

In gas phase synthesis systems, clusters form and grow via condensation, in which a monomer binds to an existing cluster. While a hard sphere equation is frequently used to predict the condensation rate coefficient, this equation neglects the influences of potential interactions and cluster internal energy on the condensation process. Here, we present a collision rate theory-Molecular Dynamics simulation approach to calculate condensation probabilities and condensation rate coefficients; we use this approach to examine atomic condensation onto 6-56 atom Au and Mg clusters. The probability of condensation depends upon the initial relative velocity ( v) between atom and cluster andmore » the initial impact parameter ( b). In all cases there is a well-defined region of b-v space where condensation is highly probable, and outside of which the condensation probability drops to zero. For Au clusters with more than 10 atoms, we find that at gas temperatures in the 300-1200 K range, the condensation rate coefficient exceeds the hard sphere rate coefficient by a factor of 1.5-2.0. Conversely, for Au clusters with 10 or fewer atoms, and for 14 atom and 28 atom Mg clusters, as cluster equilibration temperature increases the condensation rate coefficient drops to values below the hard sphere rate coefficient. Calculations also yield the self-dissociation rate coefficient, which is found to vary considerably with gas temperature. Finally, calculations results reveal that grazing (high b) atom-cluster collisions at elevated velocity (> 1000 m s -1) can result in the colliding atom rebounding (bounce) from the cluster surface or binding while another atom dissociates (replacement). In conclusion, the presented method can be applied in developing rate equations to predict material formation and growth rates in vapor phase systems.« less

Condensation and dissociation rates for gas phase metal clusters from molecular dynamics trajectory calculations.

PubMed

Yang, Huan; Goudeli, Eirini; Hogan, Christopher J

2018-04-28

In gas phase synthesis systems, clusters form and grow via condensation, in which a monomer binds to an existing cluster. While a hard-sphere equation is frequently used to predict the condensation rate coefficient, this equation neglects the influences of potential interactions and cluster internal energy on the condensation process. Here, we present a collision rate theory-molecular dynamics simulation approach to calculate condensation probabilities and condensation rate coefficients. We use this approach to examine atomic condensation onto 6-56-atom Au and Mg clusters. The probability of condensation depends upon the initial relative velocity (v) between atom and cluster and the initial impact parameter (b). In all cases, there is a well-defined region of b-v space where condensation is highly probable, and outside of which the condensation probability drops to zero. For Au clusters with more than 10 atoms, we find that at gas temperatures in the 300-1200 K range, the condensation rate coefficient exceeds the hard-sphere rate coefficient by a factor of 1.5-2.0. Conversely, for Au clusters with 10 or fewer atoms and for 14- and 28-atom Mg clusters, as cluster equilibration temperature increases, the condensation rate coefficient drops to values below the hard-sphere rate coefficient. Calculations also yield the self-dissociation rate coefficient, which is found to vary considerably with gas temperature. Finally, calculations results reveal that grazing (high b) atom-cluster collisions at elevated velocity (>1000 m s -1 ) can result in the colliding atom rebounding (bounce) from the cluster surface or binding while another atom dissociates (replacement). The presented method can be applied in developing rate equations to predict material formation and growth rates in vapor phase systems.

Condensation and dissociation rates for gas phase metal clusters from molecular dynamics trajectory calculations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yang, Huan; Goudeli, Eirini; Hogan, Christopher J.

In gas phase synthesis systems, clusters form and grow via condensation, in which a monomer binds to an existing cluster. While a hard sphere equation is frequently used to predict the condensation rate coefficient, this equation neglects the influences of potential interactions and cluster internal energy on the condensation process. Here, we present a collision rate theory-Molecular Dynamics simulation approach to calculate condensation probabilities and condensation rate coefficients; we use this approach to examine atomic condensation onto 6-56 atom Au and Mg clusters. The probability of condensation depends upon the initial relative velocity ( v) between atom and cluster andmore » the initial impact parameter ( b). In all cases there is a well-defined region of b-v space where condensation is highly probable, and outside of which the condensation probability drops to zero. For Au clusters with more than 10 atoms, we find that at gas temperatures in the 300-1200 K range, the condensation rate coefficient exceeds the hard sphere rate coefficient by a factor of 1.5-2.0. Conversely, for Au clusters with 10 or fewer atoms, and for 14 atom and 28 atom Mg clusters, as cluster equilibration temperature increases the condensation rate coefficient drops to values below the hard sphere rate coefficient. Calculations also yield the self-dissociation rate coefficient, which is found to vary considerably with gas temperature. Finally, calculations results reveal that grazing (high b) atom-cluster collisions at elevated velocity (> 1000 m s -1) can result in the colliding atom rebounding (bounce) from the cluster surface or binding while another atom dissociates (replacement). In conclusion, the presented method can be applied in developing rate equations to predict material formation and growth rates in vapor phase systems.« less

Condensation and dissociation rates for gas phase metal clusters from molecular dynamics trajectory calculations

NASA Astrophysics Data System (ADS)

Yang, Huan; Goudeli, Eirini; Hogan, Christopher J.

2018-04-01

In gas phase synthesis systems, clusters form and grow via condensation, in which a monomer binds to an existing cluster. While a hard-sphere equation is frequently used to predict the condensation rate coefficient, this equation neglects the influences of potential interactions and cluster internal energy on the condensation process. Here, we present a collision rate theory-molecular dynamics simulation approach to calculate condensation probabilities and condensation rate coefficients. We use this approach to examine atomic condensation onto 6-56-atom Au and Mg clusters. The probability of condensation depends upon the initial relative velocity (v) between atom and cluster and the initial impact parameter (b). In all cases, there is a well-defined region of b-v space where condensation is highly probable, and outside of which the condensation probability drops to zero. For Au clusters with more than 10 atoms, we find that at gas temperatures in the 300-1200 K range, the condensation rate coefficient exceeds the hard-sphere rate coefficient by a factor of 1.5-2.0. Conversely, for Au clusters with 10 or fewer atoms and for 14- and 28-atom Mg clusters, as cluster equilibration temperature increases, the condensation rate coefficient drops to values below the hard-sphere rate coefficient. Calculations also yield the self-dissociation rate coefficient, which is found to vary considerably with gas temperature. Finally, calculations results reveal that grazing (high b) atom-cluster collisions at elevated velocity (>1000 m s-1) can result in the colliding atom rebounding (bounce) from the cluster surface or binding while another atom dissociates (replacement). The presented method can be applied in developing rate equations to predict material formation and growth rates in vapor phase systems.

Numerical simulation of submicron particles formation by condensation at coals burning

NASA Astrophysics Data System (ADS)

Kortsenshteyn, N. M.; Petrov, L. V.

2017-11-01

The thermodynamic analysis of the composition of the combustion products of 15 types of coals was carried out with consideration for the formation of potassium and sodium aluminosilicates and solid and liquid slag removal. Based on the results of the analysis, the approximating temperature dependences of the concentrations of condensed components (potassium and sodium sulfates) were obtained for the cases of two-phase and single-phase equilibriums; conclusions on the comparative influence of solid and liquid slag removal on the probability of the formation of submicron particles on the combustion of coals were made. The found dependences was make it possible to perform a numerical simulation of the bulk condensation of potassium and sodium sulfate vapors upon the cooling of coal combustion products in a process flow. The number concentration and size distribution of the formed particles have been determined. Agreement with experimental data on the fraction composition of particles has been reached at a reasonable value of a free parameter of the model.

Large-Scale Hybrid Density Functional Theory Calculations in the Condensed-Phase: Ab Initio Molecular Dynamics in the Isobaric-Isothermal Ensemble

NASA Astrophysics Data System (ADS)

Ko, Hsin-Yu; Santra, Biswajit; Distasio, Robert A., Jr.; Wu, Xifan; Car, Roberto

Hybrid functionals are known to alleviate the self-interaction error in density functional theory (DFT) and provide a more accurate description of the electronic structure of molecules and materials. However, hybrid DFT in the condensed-phase has a prohibitively high associated computational cost which limits their applicability to large systems of interest. In this work, we present a general-purpose order(N) implementation of hybrid DFT in the condensed-phase using Maximally localized Wannier function; this implementation is optimized for massively parallel computing architectures. This algorithm is used to perform large-scale ab initio molecular dynamics simulations of liquid water, ice, and aqueous ionic solutions. We have performed simulations in the isothermal-isobaric ensemble to quantify the effects of exact exchange on the equilibrium density properties of water at different thermodynamic conditions. We find that the anomalous density difference between ice I h and liquid water at ambient conditions as well as the enthalpy differences between ice I h, II, and III phases at the experimental triple point (238 K and 20 Kbar) are significantly improved using hybrid DFT over previous estimates using the lower rungs of DFT This work has been supported by the Department of Energy under Grants No. DE-FG02-05ER46201 and DE-SC0008626.

Quantum Kibble-Zurek Mechanism in a Spin-1 Bose-Einstein Condensate

NASA Astrophysics Data System (ADS)

Anquez, M.; Robbins, B. A.; Bharath, H. M.; Boguslawski, M.; Hoang, T. M.; Chapman, M. S.

2016-04-01

The dynamics of a quantum phase transition are explored using slow quenches from the polar to the broken-axisymmetry phases in a small spin-1 ferromagnetic Bose-Einstein condensate. Measurements of the evolution of the spin populations reveal a power-law scaling of the temporal onset of excitations versus quench speed as predicted from quantum extensions of the Kibble-Zurek mechanism. The satisfactory agreement of the measured scaling exponent with the analytical theory and numerical simulations provides experimental confirmation of the quantum Kibble-Zurek model.

Simulator predicts transient flow for Malaysian subsea pipeline

DOE Office of Scientific and Technical Information (OSTI.GOV)

Inayat-Hussain, A.A.; Ayob, M.S.; Zain, A.B.M.

1996-04-15

In a step towards acquiring in-house capability in multiphase flow technology, Petronas Research and Scientific Services Sdn. Bhd., Kuala Lumpur, has developed two-phase flow simulation software for analyzing slow gas-condensate transient flow. Unlike its general-purpose contemporaries -- TACITE, OLGA, Traflow (OGJ, Jan. 3, 1994, p. 42; OGJ, Jan. 10, 1994, p. 52), and PLAC (AEA Technology, U.K.) -- ABASs is a dedicated software for slow transient flows generated during pigging operations in the Duyong network, offshore Malaysia. This network links the Duyong and Bekok fields to the onshore gas terminal (OGT) on the east coast of peninsular Malaysia. It predictsmore » the steady-state pressure drop vs. flow rates, condensate volume in the network, pigging dynamics including volume of produced slug, and the condensate build-up following pigging. The predictions of ABASs have been verified against field data obtained from the Duyong network. Presented here is an overview of the development, verification, and application of the ABASs software. Field data are presented for verification of the software, and several operational scenarios are simulated using the software. The field data and simulation study documented here will provide software users and developers with a further set of results on which to benchmark their own software and two-phase pipeline operating guidelines.« less

Modeling of two-phase flow instabilities during startup transients utilizing RAMONA-4B methodology

DOE Office of Scientific and Technical Information (OSTI.GOV)

Paniagua, J.; Rohatgi, U.S.; Prasad, V.

1996-10-01

RAMONA-4B code is currently under development for simulating thermal hydraulic instabilities that can occur in Boiling Water Reactors (BWRs) and the Simplified Boiling Water Reactor (SBWR). As one of the missions of RAMONA-4B is to simulate SBWR startup transients, where geysering or condensation-induced instability may be encountered, the code needs to be assessed for this application. This paper outlines the results of the assessments of the current version of RAMONA-4B and the modifications necessary for simulating the geysering or condensation-induced instability. The test selected for assessment are the geysering tests performed by Prof Aritomi (1993).

SIMULATING THE IN SITU CONDENSATION PROCESS OF SOLAR PROMINENCES

DOE Office of Scientific and Technical Information (OSTI.GOV)

Xia, C.; Keppens, R.; Antolin, P.

2014-09-10

Prominences in the solar corona are a hundredfold cooler and denser than their surroundings, with a total mass of 10{sup 13} up to 10{sup 15} g. Here, we report on the first comprehensive simulations of three-dimensional, thermally and gravitationally stratified magnetic flux ropes where in situ condensation to a prominence occurs due to radiative losses. After a gradual thermodynamic adjustment, we witness a phase where runaway cooling occurs while counter-streaming shearing flows drain off mass along helical field lines. After this drainage, a prominence-like condensation resides in concave upward field regions, and this prominence retains its overall characteristics for moremore » than two hours. While condensing, the prominence establishes a prominence-corona transition region where magnetic field-aligned thermal conduction is operative during the runaway cooling. The prominence structure represents a force-balanced state in a helical flux rope. The simulated condensation demonstrates a right-bearing barb, as a remnant of the drainage. Synthetic images at extreme ultraviolet wavelengths follow the onset of the condensation, and confirm the appearance of horns and a three-part structure for the stable prominence state, as often seen in erupting prominences. This naturally explains recent Solar Dynamics Observatory views with the Atmospheric Imaging Assembly on prominences in coronal cavities demonstrating horns.« less

The track structure in condensed matter

NASA Astrophysics Data System (ADS)

Kaplan, I. G.

1995-11-01

The physical stage of track formation in a condensed phase is discussed. For interaction of charged particles with condensed molecular media its most important specific features are: (a) the continuous oscillator strength distribution with the broak peak in the energy range 21-22 eV attributed to the collective plasmon-type state; (b) the lowering of ionization potential compared to a gas phase. These specific features must be taken into account for simulation of track structures. The great difference in mass and charge for a electron and heavy ions cause a qualitative difference in their track structures. We analyse the structure of heavy ion tracks and prove the impossibility to use the LET as a universal characteristic for the radiation action of different ions.

Simulating the control of molecular reactions via modulated light fields: from gas phase to solution

NASA Astrophysics Data System (ADS)

Thallmair, Sebastian; Keefer, Daniel; Rott, Florian; de Vivie-Riedle, Regina

2017-04-01

Over the past few years quantum control has proven to be very successful in steering molecular processes. By combining theory with experiment, even highly complex control aims were realized in the gas phase. In this topical review, we illustrate the past achievements on several examples in the molecular context. The next step for the quantum control of chemical processes is to translate the fruitful interplay between theory and experiment to the condensed phase and thus to the regime where chemical synthesis can be supported. On the theory side, increased efforts to include solvent effects in quantum control simulations were made recently. We discuss two major concepts, namely an implicit description of the environment via the density matrix algorithm and an explicit inclusion of solvent molecules. By application to chemical reactions, both concepts conclude that despite environmental perturbations leading to more complex control tasks, efficient quantum control in the condensed phase is still feasible.

Pressure-induced metallization of condensed phase β-HMX under shock loadings via molecular dynamics simulations in conjunction with multi-scale shock technique.

PubMed

Ge, Ni-Na; Wei, Yong-Kai; Zhao, Feng; Chen, Xiang-Rong; Ji, Guang-Fu

2014-07-01

The electronic structure and initial decomposition in high explosive HMX under conditions of shock loading are examined. The simulation is performed using quantum molecular dynamics in conjunction with multi-scale shock technique (MSST). A self-consistent charge density-functional tight-binding (SCC-DFTB) method is adapted. The results show that the N-N-C angle has a drastic change under shock wave compression along lattice vector b at shock velocity 11 km/s, which is the main reason that leads to an insulator-to-metal transition for the HMX system. The metallization pressure (about 130 GPa) of condensed-phase HMX is predicted firstly. We also detect the formation of several key products of condensed-phase HMX decomposition, such as NO2, NO, N2, N2O, H2O, CO, and CO2, and all of them have been observed in previous experimental studies. Moreover, the initial decomposition products include H2 due to the C-H bond breaking as a primary reaction pathway at extreme condition, which presents a new insight into the initial decomposition mechanism of HMX under shock loading at the atomistic level.

Analytical Approach in DeCoM

NASA Technical Reports Server (NTRS)

Patel, Deepak

2011-01-01

There are many papers on describing a LHP as an overall system, but few detail on the condenser section of a loop heat pipe. The DeCoM (Deepak Condenser Model) method utilizes user set initial parameters in-order to simulate a condenser by calculating the interactions between the fluid and the wall. Equations are derived for two sections of the condenser: a two-phase section and a subcooled (liquid) section. All Equations are based upon the conservation of energy theory, from which fluid temperature, and fluid quality values are solved. In order to solve for the heat transfer value, between fluid and the wall in two phase section, the Lockhart-Martinelli correlation method was implemented as a solution approach. For Liquid phase, the Reynolds number was used in-order to differentiate the flow state, from either turbulent or laminar, and Nusselt number was used to solve for the film coefficient. To represent these calculations for both sections a flow chart is presented in order to display the execution process of DeCoM. The benefit of DeCoM is that it is capable of performing preliminary analysis without requiring a license and without much of users knowledge on condensers.

Excess electron localization in solvated DNA bases.

PubMed

Smyth, Maeve; Kohanoff, Jorge

2011-06-10

We present a first-principles molecular dynamics study of an excess electron in condensed phase models of solvated DNA bases. Calculations on increasingly large microsolvated clusters taken from liquid phase simulations show that adiabatic electron affinities increase systematically upon solvation, as for optimized gas-phase geometries. Dynamical simulations after vertical attachment indicate that the excess electron, which is initially found delocalized, localizes around the nucleobases within a 15 fs time scale. This transition requires small rearrangements in the geometry of the bases.

Excess Electron Localization in Solvated DNA Bases

DOE Office of Scientific and Technical Information (OSTI.GOV)

Smyth, Maeve; Kohanoff, Jorge

2011-06-10

We present a first-principles molecular dynamics study of an excess electron in condensed phase models of solvated DNA bases. Calculations on increasingly large microsolvated clusters taken from liquid phase simulations show that adiabatic electron affinities increase systematically upon solvation, as for optimized gas-phase geometries. Dynamical simulations after vertical attachment indicate that the excess electron, which is initially found delocalized, localizes around the nucleobases within a 15 fs time scale. This transition requires small rearrangements in the geometry of the bases.

Solid H2 in the interstellar medium

NASA Astrophysics Data System (ADS)

Füglistaler, A.; Pfenniger, D.

2018-06-01

Context. Condensation of H2 in the interstellar medium (ISM) has long been seen as a possibility, either by deposition on dust grains or thanks to a phase transition combined with self-gravity. H2 condensation might explain the observed low efficiency of star formation and might help to hide baryons in spiral galaxies. Aims: Our aim is to quantify the solid fraction of H2 in the ISM due to a phase transition including self-gravity for different densities and temperatures in order to use the results in more complex simulations of the ISM as subgrid physics. Methods: We used molecular dynamics simulations of fluids at different temperatures and densities to study the formation of solids. Once the simulations reached a steady state, we calculated the solid mass fraction, energy increase, and timescales. By determining the power laws measured over several orders of magnitude, we extrapolated to lower densities the higher density fluids that can be simulated with current computers. Results: The solid fraction and energy increase of fluids in a phase transition are above 0.1 and do not follow a power law. Fluids out of a phase transition are still forming a small amount of solids due to chance encounters of molecules. The solid mass fraction and energy increase of these fluids are linearly dependent on density and can easily be extrapolated. The timescale is below one second, the condensation can be considered instantaneous. Conclusions: The presence of solid H2 grains has important dynamic implications on the ISM as they may be the building blocks for larger solid bodies when gravity is included. We provide the solid mass fraction, energy increase, and timescales for high density fluids and extrapolation laws for lower densities.

Numerical simulation of condensation on structured surfaces.

PubMed

Fu, Xiaowu; Yao, Zhaohui; Hao, Pengfei

2014-11-25

Condensation of liquid droplets on solid surfaces happens widely in nature and industrial processes. This phase-change phenomenon has great effect on the performance of some microfluidic devices. On the basis of micro- and nanotechnology, superhydrophobic structured surfaces can be well-fabricated. In this work, the nucleating and growth of droplets on different structured surfaces are investigated numerically. The dynamic behavior of droplets during the condensation is simulated by the multiphase lattice Boltzmann method (LBM), which has the ability to incorporate the microscopic interactions, including fluid-fluid interaction and fluid-surface interaction. The results by the LBM show that, besides the chemical properties of surfaces, the topography of structures on solid surfaces influences the condensation process. For superhydrophobic surfaces, the spacing and height of microridges have significant influence on the nucleation sites. This mechanism provides an effective way for prevention of wetting on surfaces in engineering applications. Moreover, it suggests a way to prevent ice formation on surfaces caused by the condensation of subcooled water. For hydrophilic surfaces, however, microstructures may be submerged by the liquid films adhering to the surfaces. In this case, microstructures will fail to control the condensation process. Our research provides an optimized way for designing surfaces for condensation in engineering systems.

Molecular simulation of steady-state evaporation and condensation in the presence of a non-condensable gas

NASA Astrophysics Data System (ADS)

Liang, Zhi; Keblinski, Pawel

2018-02-01

Using molecular dynamics simulations, we study evaporation and condensation of fluid Ar in the presence of a non-condensable Ne gas in a nanochannel. The evaporation and condensation are driven by the temperature difference, ΔTL, between the evaporating and condensing liquid surfaces. The steady-state evaporation and condensation fluxes (JMD) are also affected by the Ne concentration, ρNe, and the nanochannel length. We find that across a wide range of ΔTL and ρNe, JMD is in good agreement with the prediction from Stefan's law and from Schrage relationships. Furthermore, for ΔTL less than ˜20% of the absolute average temperature, we find that both steady-state heat and mass fluxes are proportional to ΔTL. This allows us to determine the interfacial resistance to the heat and mass transfer and compare it with the corresponding resistances in the gas phase. In this context, we derive an analytical expression for the effective thermal conductivity of the gas region in the nanochannel and the mass transport interfacial resistance equivalent length, i.e., the length of the nanochannel for which the resistance to the mass flow is the same as the interfacial resistance to the mass flow.

Condensation with two constraints and disorder

NASA Astrophysics Data System (ADS)

Barré, J.; Mangeolle, L.

2018-04-01

We consider a set of positive random variables obeying two additive constraints, a linear and a quadratic one; these constraints mimic the conservation laws of a dynamical system. In the simplest setting, without disorder, it is known that such a system may undergo a ‘condensation’ transition, whereby one random variable becomes much larger than the others; this transition has been related to the spontaneous appearance of non linear localized excitations in certain nonlinear chains, called breathers. Motivated by the study of breathers in a disordered discrete nonlinear Schrödinger equation, we study different instances of this problem in presence of a quenched disorder. Unless the disorder is too strong, the phase diagram looks like the one without disorder, with a transition separating a fluid phase, where all variables have the same order of magnitude, and a condensed phase, where one variable is much larger than the others. We then show that the condensed phase exhibits various degrees of ‘intermediate symmetry breaking’: the site hosting the condensate is chosen neither uniformly at random, nor is it fixed by the disorder realization. Throughout the article, our heuristic arguments are complemented with direct Monte Carlo simulations.

Overcoming time scale and finite size limitations to compute nucleation rates from small scale well tempered metadynamics simulations.

PubMed

Salvalaglio, Matteo; Tiwary, Pratyush; Maggioni, Giovanni Maria; Mazzotti, Marco; Parrinello, Michele

2016-12-07

Condensation of a liquid droplet from a supersaturated vapour phase is initiated by a prototypical nucleation event. As such it is challenging to compute its rate from atomistic molecular dynamics simulations. In fact at realistic supersaturation conditions condensation occurs on time scales that far exceed what can be reached with conventional molecular dynamics methods. Another known problem in this context is the distortion of the free energy profile associated to nucleation due to the small, finite size of typical simulation boxes. In this work the problem of time scale is addressed with a recently developed enhanced sampling method while contextually correcting for finite size effects. We demonstrate our approach by studying the condensation of argon, and showing that characteristic nucleation times of the order of magnitude of hours can be reliably calculated. Nucleation rates spanning a range of 10 orders of magnitude are computed at moderate supersaturation levels, thus bridging the gap between what standard molecular dynamics simulations can do and real physical systems.

Overcoming time scale and finite size limitations to compute nucleation rates from small scale well tempered metadynamics simulations

NASA Astrophysics Data System (ADS)

Salvalaglio, Matteo; Tiwary, Pratyush; Maggioni, Giovanni Maria; Mazzotti, Marco; Parrinello, Michele

2016-12-01

Condensation of a liquid droplet from a supersaturated vapour phase is initiated by a prototypical nucleation event. As such it is challenging to compute its rate from atomistic molecular dynamics simulations. In fact at realistic supersaturation conditions condensation occurs on time scales that far exceed what can be reached with conventional molecular dynamics methods. Another known problem in this context is the distortion of the free energy profile associated to nucleation due to the small, finite size of typical simulation boxes. In this work the problem of time scale is addressed with a recently developed enhanced sampling method while contextually correcting for finite size effects. We demonstrate our approach by studying the condensation of argon, and showing that characteristic nucleation times of the order of magnitude of hours can be reliably calculated. Nucleation rates spanning a range of 10 orders of magnitude are computed at moderate supersaturation levels, thus bridging the gap between what standard molecular dynamics simulations can do and real physical systems.

Simulation of Containment Atmosphere Mixing and Stratification Experiment in the ThAI Facility with a CFD Code

DOE Office of Scientific and Technical Information (OSTI.GOV)

Babic, Miroslav; Kljenak, Ivo; Mavko, Borut

2006-07-01

The CFD code CFX4.4 was used to simulate an experiment in the ThAI facility, which was designed for investigation of thermal-hydraulic processes during a severe accident inside a Light Water Reactor containment. In the considered experiment, air was initially present in the vessel, and helium and steam were injected during different phases of the experiment at various mass flow rates and at different locations. The main purpose of the proposed work was to assess the capabilities of the CFD code to reproduce the atmosphere structure with a three-dimensional model, coupled with condensation models proposed by the authors. A three-dimensional modelmore » of the ThAI vessel for the CFX4.4 code was developed. The flow in the simulation domain was modeled as single-phase. Steam condensation on vessel walls was modeled as a sink of mass and energy using a correlation that was originally developed for an integral approach. A simple model of bulk phase change was also included. Calculated time-dependent variables together with temperature and volume fraction distributions at the end of different experiment phases are compared to experimental results. (authors)« less

Molecular dynamics simulation of Bu4N+ in dimethylformamide: Solvation-induced volume changes

NASA Astrophysics Data System (ADS)

Kiselev, M. G.; Safonova, L. P.

2011-06-01

The structure of the Bu4N+-dimethylformamide system in the condensed and gas phases was studied by molecular dynamics simulation and quantum-chemical calculations. The calculation results were used to reveal the role played by steric effects in the volumetric characteristics of ion solvation.

Development and evaluation of the aerosol dynamics and gas phase chemistry model ADCHEM

NASA Astrophysics Data System (ADS)

Roldin, P.; Swietlicki, E.; Schurgers, G.; Arneth, A.; Lehtinen, K. E. J.; Boy, M.; Kulmala, M.

2011-06-01

The aim of this work was to develop a model suited for detailed studies of aerosol dynamics, gas and particle phase chemistry within urban plumes, from local scale (1 × 1 km2) to regional scale. This article describes and evaluates the trajectory model for Aerosol Dynamics, gas and particle phase CHEMistry and radiative transfer (ADCHEM). The model treats both vertical and horizontal dispersion perpendicular to an air mass trajectory (2-space dimensions). The Lagrangian approach enables a more detailed representation of the aerosol dynamics, gas and particle phase chemistry and a finer spatial and temporal resolution compared to that of available regional 3D-CTMs. These features make it among others well suited for urban plume studies. The aerosol dynamics model includes Brownian coagulation, dry deposition, wet deposition, in-cloud processing, condensation, evaporation, primary particle emissions and homogeneous nucleation. The organic mass partitioning was either modeled with a 2-dimensional volatility basis set (2D-VBS) or with the traditional two-product model approach. In ADCHEM these models consider the diffusion limited and particle size dependent condensation and evaporation of 110 and 40 different organic compounds respectively. The gas phase chemistry model calculates the gas phase concentrations of 61 different species, using 130 different chemical reactions. Daily isoprene and monoterpene emissions from European forests were simulated separately with the vegetation model LPJ-GUESS, and included as input to ADCHEM. ADCHEM was used to simulate the ageing of the urban plumes from the city of Malmö in southern Sweden (280 000 inhabitants). Several sensitivity tests were performed concerning the number of size bins, size structure method, aerosol dynamic processes, vertical and horizontal mixing, coupled or uncoupled condensation and the secondary organic aerosol formation. The simulations show that the full-stationary size structure gives accurate results with little numerical diffusion when more than 50 size bins are used between 1.5 and 2500 nm, while the moving-center method is preferable when only a few size bins are selected. The particle number size distribution in the center of the urban plume from Malmö was mainly affected by dry deposition, coagulation and vertical dilution. The modeled PM2.5 mass was dominated by organic material, nitrate, sulfate and ammonium. If the condensation of HNO3 and NH3 was treated as a coupled process (pH independent) the model gave lower nitrate PM2.5 mass than if considering uncoupled condensation. Although the time of ageing from that SOA precursors are emitted until condensable products are formed is substantially different with the 2D-VBS and two product model, the models gave similar total organic mass concentrations.

Thermodynamic integration based on classical atomistic simulations to determine the Gibbs energy of condensed phases: Calculation of the aluminum-zirconium system

NASA Astrophysics Data System (ADS)

Harvey, J.-P.; Gheribi, A. E.; Chartrand, P.

2012-12-01

In this work, an in silico procedure to generate a fully coherent set of thermodynamic properties obtained from classical molecular dynamics (MD) and Monte Carlo (MC) simulations is proposed. The procedure is applied to the Al-Zr system because of its importance in the development of high strength Al-Li alloys and of bulk metallic glasses. Cohesive energies of the studied condensed phases of the Al-Zr system (the liquid phase, the fcc solid solution, and various orthorhombic stoichiometric compounds) are calculated using the modified embedded atom model (MEAM) in the second-nearest-neighbor formalism (2NN). The Al-Zr MEAM-2NN potential is parameterized in this work using ab initio and experimental data found in the literature for the AlZr3-L12 structure, while its predictive ability is confirmed for several other solid structures and for the liquid phase. The thermodynamic integration (TI) method is implemented in a general MC algorithm in order to evaluate the absolute Gibbs energy of the liquid and the fcc solutions. The entropy of mixing calculated from the TI method, combined to the enthalpy of mixing and the heat capacity data generated from MD/MC simulations performed in the isobaric-isothermal/canonical (NPT/NVT) ensembles are used to parameterize the Gibbs energy function of all the condensed phases in the Al-rich side of the Al-Zr system in a CALculation of PHAse Diagrams (CALPHAD) approach. The modified quasichemical model in the pair approximation (MQMPA) and the cluster variation method (CVM) in the tetrahedron approximation are used to define the Gibbs energy of the liquid and the fcc solid solution respectively for their entire range of composition. Thermodynamic and structural data generated from our MD/MC simulations are used as input data to parameterize these thermodynamic models. A detailed analysis of the validity and transferability of the Al-Zr MEAM-2NN potential is presented throughout our work by comparing the predicted properties obtained from this formalism with available ab initio and experimental data for both liquid and solid phases.

The Exoplanet Cloud Atlas

NASA Astrophysics Data System (ADS)

Gao, Peter; Marley, Mark S.; Morley, Caroline; Fortney, Jonathan J.

2017-10-01

Clouds have been readily inferred from observations of exoplanet atmospheres, and there exists great variability in cloudiness between planets, such that no clear trend in exoplanet cloudiness has so far been discerned. Equilibrium condensation calculations suggest a myriad of species - salts, sulfides, silicates, and metals - could condense in exoplanet atmospheres, but how they behave as clouds is uncertain. The behavior of clouds - their formation, evolution, and equilibrium size distribution - is controlled by cloud microphysics, which includes processes such as nucleation, condensation, and evaporation. In this work, we explore the cloudy exoplanet phase space by using a cloud microphysics model to simulate a suite of cloud species ranging from cooler condensates such as KCl/ZnS, to hotter condensates like perovskite and corundum. We investigate how the cloudiness and cloud particle sizes of exoplanets change due to variations in temperature, metallicity, gravity, and cloud formation mechanisms, and how these changes may be reflected in current and future observations. In particular, we will evaluate where in phase space could cloud spectral features be observable using JWST MIRI at long wavelengths, which will be dependent on the cloud particle size distribution and cloud species.

Quantum fluids of light in acoustic lattices

NASA Astrophysics Data System (ADS)

Cerda-Méndez, E. A.; Krizhanovskii, D. N.; Skolnick, M. S.; Santos, P. V.

2018-01-01

In this topical review, we report on the recent advances on the manipulation of hybrid light-matter quasi-particles called exciton-polaritons and their quantum condensed phases by means of acoustic and static periodic potentials. Polaritons are a superposition of photons and excitons and form in optical microcavities with quantum wells embedded in it. They are low-mass bosons in the dilute limit and have strong inter-particle interactions inherited from the excitonic component. Their capability to form quantum-condensed phases at temperatures in the kelvin range and to behave like quantum fluids makes them very attractive for novel solid-state devices. Since their de Broglie wavelength is of the order of a few micrometers, polaritons can be manipulated using static or dynamic potentials with micrometer scales. We present here a summary of the techniques used to submit polaritons and their condensed phases to periodic potentials, with an emphasis in dynamic ones produced by surface acoustic waves. We discuss the interesting phenomena that occur under such a modulation, such as condensation in excited states of the Brillouin zone, fragmentation of a condensate, formation of self-localized wavepackets, and Dirac and massive polaritons in static hexagonal and kagome lattices, respectively. The different techniques explored open the way to implement polariton-based quantum simulators, nano-optomechanic resonators and polaritonic topological insulators.

A reduced theoretical model for estimating condensation effects in combustion-heated hypersonic tunnel

NASA Astrophysics Data System (ADS)

Lin, L.; Luo, X.; Qin, F.; Yang, J.

2018-03-01

As one of the combustion products of hydrocarbon fuels in a combustion-heated wind tunnel, water vapor may condense during the rapid expansion process, which will lead to a complex two-phase flow inside the wind tunnel and even change the design flow conditions at the nozzle exit. The coupling of the phase transition and the compressible flow makes the estimation of the condensation effects in such wind tunnels very difficult and time-consuming. In this work, a reduced theoretical model is developed to approximately compute the nozzle-exit conditions of a flow including real-gas and homogeneous condensation effects. Specifically, the conservation equations of the axisymmetric flow are first approximated in the quasi-one-dimensional way. Then, the complex process is split into two steps, i.e., a real-gas nozzle flow but excluding condensation, resulting in supersaturated nozzle-exit conditions, and a discontinuous jump at the end of the nozzle from the supersaturated state to a saturated state. Compared with two-dimensional numerical simulations implemented with a detailed condensation model, the reduced model predicts the flow parameters with good accuracy except for some deviations caused by the two-dimensional effect. Therefore, this reduced theoretical model can provide a fast, simple but also accurate estimation of the condensation effect in combustion-heated hypersonic tunnels.

Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces.

PubMed

Spezia, Riccardo; Martínez-Nuñez, Emilio; Vazquez, Saulo; Hase, William L

2017-04-28

In this Introduction, we show the basic problems of non-statistical and non-equilibrium phenomena related to the papers collected in this themed issue. Over the past few years, significant advances in both computing power and development of theories have allowed the study of larger systems, increasing the time length of simulations and improving the quality of potential energy surfaces. In particular, the possibility of using quantum chemistry to calculate energies and forces 'on the fly' has paved the way to directly study chemical reactions. This has provided a valuable tool to explore molecular mechanisms at given temperatures and energies and to see whether these reactive trajectories follow statistical laws and/or minimum energy pathways. This themed issue collects different aspects of the problem and gives an overview of recent works and developments in different contexts, from the gas phase to the condensed phase to excited states.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'. © 2017 The Author(s).

Two-phase/two-phase heat exchanger analysis

NASA Technical Reports Server (NTRS)

Kim, Rhyn H.

1992-01-01

A capillary pumped loop (CPL) system with a condenser linked to a double two-phase heat exchanger is analyzed numerically to simulate the performance of the system from different starting conditions to a steady state condition based on a simplified model. Results of the investigation are compared with those of similar apparatus available in the Space Station applications of the CPL system with a double two-phase heat exchanger.

Chiral phase transition of three flavor QCD with nonzero magnetic field using standard staggered fermions

NASA Astrophysics Data System (ADS)

Tomiya, Akio; Ding, Heng-Tong; Mukherjee, Swagato; Schmidt, Christian; Wang, Xiao-Dan

2018-03-01

Lattice simulations for (2+1)-flavor QCD with external magnetic field demon-strated that the quark mass is one of the important parameters responsible for the (inverse) magnetic catalysis. We discuss the dependences of chiral condensates and susceptibilities, the Polyakov loop on the magnetic field and quark mass in three degenerate flavor QCD. The lattice simulations are performed using standard staggered fermions and the plaquette action with spatial sizes Nσ = 16 and 24 and a fixed temporal size Nτ = 4. The value of the quark masses are chosen such that the system undergoes a first order chiral phase transition and crossover with zero magnetic field. We find that in light mass regime, the quark chiral condensate undergoes magnetic catalysis in the whole temperature region and the phase transition tend to become stronger as the magnetic field increases. In crossover regime, deconfinement transition temperature is shifted by the magnetic field when quark mass ma is less than 0:4. The lattice cutoff effects are also discussed.

Coarsening dynamics in condensing zero-range processes and size-biased birth death chains

NASA Astrophysics Data System (ADS)

Jatuviriyapornchai, Watthanan; Grosskinsky, Stefan

2016-05-01

Zero-range processes with decreasing jump rates are well known to exhibit a condensation transition under certain conditions on the jump rates, and the dynamics of this transition continues to be a subject of current research interest. Starting from homogeneous initial conditions, the time evolution of the condensed phase exhibits an interesting coarsening phenomenon of mass transport between cluster sites characterized by a power law. We revisit the approach in Godrèche (2003 J. Phys. A: Math. Gen. 36 6313) to derive effective single site dynamics which form a nonlinear birth death chain describing the coarsening behavior. We extend these results to a larger class of parameter values, and introduce a size-biased version of the single site process, which provides an effective tool to analyze the dynamics of the condensed phase without finite size effects and is the main novelty of this paper. Our results are based on a few heuristic assumptions and exact computations, and are corroborated by detailed simulation data.

Revisiting kinetic boundary conditions at the surface of fuel droplet hydrocarbons: An atomistic computational fluid dynamics simulation

PubMed Central

Nasiri, Rasoul

2016-01-01

The role of boundary conditions at the interface for both Boltzmann equation and the set of Navier-Stokes equations have been suggested to be important for studying of multiphase flows such as evaporation/condensation process which doesn’t always obey the equilibrium conditions. Here we present aspects of transition-state theory (TST) alongside with kinetic gas theory (KGT) relevant to the study of quasi-equilibrium interfacial phenomena and the equilibrium gas phase processes, respectively. A two-state mathematical model for long-chain hydrocarbons which have multi-structural specifications is introduced to clarify how kinetics and thermodynamics affect evaporation/condensation process at the surface of fuel droplet, liquid and gas phases and then show how experimental observations for a number of n-alkane may be reproduced using a hybrid framework TST and KGT with physically reasonable parameters controlling the interface, gas and liquid phases. The importance of internal activation dynamics at the surface of n-alkane droplets is established during the evaporation/condensation process. PMID:27215897

Analysis of nanoscale two-phase flow of argon using molecular dynamics

NASA Astrophysics Data System (ADS)

Verma, Abhishek Kumar; Kumar, Rakesh

2014-12-01

Two phase flows through micro and nanochannels have attracted a lot of attention because of their immense applicability to many advanced fields such as MEMS/NEMS, electronic cooling, bioengineering etc. In this work, a molecular dynamics simulation method is employed to study the condensation process of superheated argon vapor force driven flow through a nanochannel combining fluid flow and heat transfer. A simple and effective particle insertion method is proposed to model phase change of argon based on non-periodic boundary conditions in the simulation domain. Starting from a crystalline solid wall of channel, the condensation process evolves from a transient unsteady state where we study the influence of different wall temperatures and fluid wall interactions on interfacial and heat transport properties of two phase flows. Subsequently, we analyzed transient temperature, density and velocity fields across the channel and their dependency on varying wall temperature and fluid wall interaction, after a dynamic equilibrium is achieved in phase transition. Quasi-steady nonequilibrium temperature profile, heat flux and interfacial thermal resistance were analyzed. The results demonstrate that the molecular dynamics method, with the proposed particle insertion method, effectively solves unsteady nonequilibrium two phase flows at nanoscale resolutions whose interphase between liquid and vapor phase is typically of the order of a few molecular diameters.

Modeling Secondary Organic Aerosols over Europe: Impact of Activity Coefficients and Viscosity

NASA Astrophysics Data System (ADS)

Kim, Y.; Sartelet, K.; Couvidat, F.

2014-12-01

Semi-volatile organic species (SVOC) can condense on suspended particulate materials (PM) in the atmosphere. The modeling of condensation/evaporation of SVOC often assumes that gas-phase and particle-phase concentrations are at equilibrium. However, recent studies show that secondary organic aerosols (SOA) may not be accurately represented by an equilibrium approach between the gas and particle phases, because organic aerosols in the particle phase may be very viscous. The condensation in the viscous liquid phase is limited by the diffusion from the surface of PM to its core. Using a surrogate approach to represent SVOC, depending on the user's choice, the secondary organic aerosol processor (SOAP) may assume equilibrium or model dynamically the condensation/evaporation between the gas and particle phases to take into account the viscosity of organic aerosols. The model is implemented in the three-dimensional chemistry-transport model of POLYPHEMUS. In SOAP, activity coefficients for organic mixtures can be computed using UNIFAC for short-range interactions between molecules and AIOMFAC to also take into account the effect of inorganic species on activity coefficients. Simulations over Europe are performed and POLYPHEMUS/SOAP is compared to POLYPHEMUS/H2O, which was previously used to model SOA using the equilibrium approach with activity coefficients from UNIFAC. Impacts of the dynamic approach on modeling SOA over Europe are evaluated. The concentrations of SOA using the dynamic approach are compared with those using the equilibrium approach. The increase of computational cost is also evaluated.

Numerical simulation of superheated vapor bubble rising in stagnant liquid

NASA Astrophysics Data System (ADS)

Samkhaniani, N.; Ansari, M. R.

2017-09-01

In present study, the rising of superheated vapor bubble in saturated liquid is simulated using volume of fluid method in OpenFOAM cfd package. The surface tension between vapor-liquid phases is considered using continuous surface force method. In order to reduce spurious current near interface, Lafaurie smoothing filter is applied to improve curvature calculation. Phase change is considered using Tanasawa mass transfer model. The variation of saturation temperature in vapor bubble with local pressure is considered with simplified Clausius-Clapeyron relation. The couple velocity-pressure equation is solved using PISO algorithm. The numerical model is validated with: (1) isothermal bubble rising and (2) one-dimensional horizontal film condensation. Then, the shape and life time history of single superheated vapor bubble are investigated. The present numerical study shows vapor bubble in saturated liquid undergoes boiling and condensation. It indicates bubble life time is nearly linear proportional with bubble size and superheat temperature.

Design Rule for Colloidal Crystals of DNA-Functionalized Particles

NASA Astrophysics Data System (ADS)

Martinez-Veracoechea, Francisco J.; Mladek, Bianca M.; Tkachenko, Alexei V.; Frenkel, Daan

2011-07-01

We report a Monte Carlo simulation study of the phase behavior of colloids coated with long, flexible DNA chains. We find that an important change occurs in the phase diagram when the number of DNAs per colloid is decreased below a critical value. In this case, the triple point disappears and the condensed phase that coexists with the vapor is always liquid. Our simulations thus explain why, in the dilute solutions typically used in experiments, colloids coated with a small number of DNA strands cannot crystallize. We understand this behavior in terms of the discrete nature of DNA binding.

A Statistical Comparison of PSC Model Simulations and POAM Observations

NASA Technical Reports Server (NTRS)

Strawa, A. W.; Drdla, K.; Fromm, M.; Bokarius, K.; Gore, Warren J. (Technical Monitor)

2002-01-01

A better knowledge of PSC composition and formation mechanisms is important to better understand and predict stratospheric ozone depletion. Several past studies have attempted to compare modeling results with satellite observations. These comparisons have concentrated on case studies. In this paper we adopt a statistical approach. POAM PSC observations from several Arctic winters are categorized into Type Ia and Ib PSCs using a technique based on Strawa et al. The discrimination technique has been modified to employ the wavelengths dependence of the extinction signal at all wavelengths rather than only at 603 and 10 18 nm. Winter-long simulations for the 1999-2000 Arctic winter have been made using the IMPACT model. These simulations have been constrained by aircraft observations made during the SOLVE/THESEO 2000 campaign. A complete set of winter-long simulations was run for several different microphysical and PSC formation scenarios. The simulations give us perfect knowledge of PSC type (Ia, Ib, or II), composition, especially condensed phase HNO3 which is important for denitrification, and condensed phase H2O. Comparisons are made between the simulation and observation of PSC extinction at 1018 rim versus wavelength dependence, winter-long percentages of Ia and Ib occurrence, and temporal and altitude trends of the PSCs. These comparisons allow us to comment on how realistic some modeling scenarios are.

Heat and Mass Transfer with Condensation in Capillary Porous Bodies

PubMed Central

2014-01-01

The purpose of this present work is related to wetting process analysis caused by condensation phenomena in capillary porous material by using a numerical simulation. Special emphasis is given to the study of the mechanism involved and the evaluation of classical theoretical models used as a predictive tool. A further discussion will be given for the distribution of the liquid phase for both its pendular and its funicular state and its consequence on diffusion coefficients of the mathematical model used. Beyond the complexity of the interaction effects between vaporisation-condensation processes on the gas-liquid interfaces, the comparison between experimental and numerical simulations permits to identify the specific contribution and the relative part of mass and energy transport parameters. This analysis allows us to understand the contribution of each part of the mathematical model used and to simplify the study. PMID:24688366

Antifoam Degradation Products in Off Gas and Condensate of Sludge Batch 9 Simulant Nitric-Formic Flowsheet Testing for the Defense Waste Processing Facility

DOE Office of Scientific and Technical Information (OSTI.GOV)

Smith, T.

Ten chemical processing cell (CPC) experiments were performed using simulant to evaluate Sludge Batch 9 for sludge-only and coupled processing using the nitric-formic flowsheet in the Defense Waste Processing Facility (DWPF). Sludge Receipt and Adjustment Tank (SRAT) and Slurry Mix Evaporator (SME) cycles were performed on eight of the ten. The other two were SRAT cycles only. Samples of the condensate, sludge, and off gas were taken to monitor the chemistry of the CPC experiments. The Savannah River National Laboratory (SRNL) has previously shown antifoam decomposes to form flammable organic products, (hexamethyldisiloxane (HMDSO), trimethylsilanol (TMS), and propanal), that are presentmore » in the vapor phase and condensate of the CPC vessels. To minimize antifoam degradation product formation, a new antifoam addition strategy was implemented at SRNL and DWPF to add antifoam undiluted.« less

Transferability of polarizable models for ion-water electrostatic interaction

NASA Astrophysics Data System (ADS)

Masia, Marco

2009-06-01

Studies of ion-water systems at condensed phase and at interfaces have pointed out that molecular and ionic polarization plays an important role for many phenomena ranging from hydrogen bond dynamics to water interfaces' structure. Classical and ab initio Molecular Dynamics simulations reveal that induced dipole moments at interfaces (e.g. air-water and water-protein) are usually high, hinting that polarizable models to be implemented in classical force fields should be very accurate in reproducing the electrostatic properties of the system. In this paper the electrostatic properties of three classical polarizable models for ion-water interaction are compared with ab initio results both at gas and condensed phase. For Li+- water and Cl--water dimers the reproducibility of total dipole moments obtained with high level quantum chemical calculations is studied; for the same ions in liquid water, Car-Parrinello Molecular Dynamics simulations are used to compute the time evolution of ionic and molecular dipole moments, which are compared with the classical models. The PD2-H2O model developed by the author and coworkers [Masia et al. J. Chem. Phys. 2004, 121, 7362] together with the gaussian intermolecular damping for ion-water interaction [Masia et al. J. Chem. Phys. 2005, 123, 164505] showed to be the fittest in reproducing the ab initio results from gas to condensed phase, allowing for force field transferability.

Cluster growth mechanisms in Lennard-Jones fluids: A comparison between molecular dynamics and Brownian dynamics simulations

NASA Astrophysics Data System (ADS)

Jung, Jiyun; Lee, Jumin; Kim, Jun Soo

2015-03-01

We present a simulation study on the mechanisms of a phase separation in dilute fluids of Lennard-Jones (LJ) particles as a model of self-interacting molecules. Molecular dynamics (MD) and Brownian dynamics (BD) simulations of the LJ fluids are employed to model the condensation of a liquid droplet in the vapor phase and the mesoscopic aggregation in the solution phase, respectively. With emphasis on the cluster growth at late times well beyond the nucleation stage, we find that the growth mechanisms can be qualitatively different: cluster diffusion and coalescence in the MD simulations and Ostwald ripening in the BD simulations. We also show that the rates of the cluster growth have distinct scaling behaviors during cluster growth. This work suggests that in the solution phase the random Brownian nature of the solute dynamics may lead to the Ostwald ripening that is qualitatively different from the cluster coalescence in the vapor phase.

Development and evaluation of the aerosol dynamic and gas phase chemistry model ADCHEM

NASA Astrophysics Data System (ADS)

Roldin, P.; Swietlicki, E.; Schurgers, G.; Arneth, A.; Lehtinen, K. E. J.; Boy, M.; Kulmala, M.

2010-08-01

The aim of this work was to develop a model ideally suited for detailed studies on aerosol dynamics, gas and particle phase chemistry within urban plumes, from local scale (1×1 km2) to regional or global scale. This article describes and evaluates the trajectory model for Aerosol Dynamics, gas and particle phase CHEMistry and radiative transfer (ADCHEM), which has been developed and used at Lund University since 2007. The model treats both vertical and horizontal dispersion perpendicular to an air mass trajectory (2-space dimensions), which is not treated in Lagrangian box-models (0-space dimensions). The Lagrangian approach enables a more detailed representation of the aerosol dynamics, gas and particle phase chemistry and a finer spatial and temporal resolution compared to that of available regional 3D-CTMs. These features make it among others ideally suited for urban plume studies. The aerosol dynamics model includes Brownian coagulation, dry deposition, wet deposition, in-cloud processing, condensation, evaporation, primary particle emissions and homogeneous nucleation. The gas phase chemistry model calculates the gas phase concentrations of 63 different species, using 119 different chemical reactions. Daily isoprene and monoterpene emissions from European forests were simulated separately with the vegetation model LPJ-GUESS, and included as input to ADCHEM. ADCHEM was used to simulate the ageing of the urban plumes from the city of Malmö in Southern Sweden (280 000 inhabitants). Several sensitivity tests were performed concerning the number of size bins, size structure method, coupled or uncoupled condensation, the volatility basis set (VBS) or traditional 2-product model for secondary organic aerosol formation, different aerosol dynamic processes and vertical and horizontal mixing. The simulations show that the full-stationary size structure gives accurate results with little numerical diffusion when more than 50 size bins are used between 1.5 and 2500 nm, while the moving-center method is preferable when only a few size bins are selected. The particle number size distribution in the center of the urban plume from Malmö is mainly affected by dry deposition, coagulation and condensation, and is relatively insensitive to moderate variations in the vertical and horizontal mixing, as long as the mixing height is relatively uniform. The modeled PM2.5 was dominated by organics, nitrate, sulfate and ammonium. If treating the condensation of HNO3 and NH3 as a coupled process (pH independent) the model gave lower nitrate PM2.5 values than if considering uncoupled condensation. However, both methods gave similar and significant temporal variation in the particulate nitrate content, primarily due to fluctuation in the relative humidity.

Condensed phase conversion and growth of nanorods and other materials instead of from vapor

DOEpatents

Geohegan, David B.; Seals, Roland D.; Puretzky, Alex A.; Fan, Xudong

2010-10-19

Compositions, systems and methods are described for condensed phase conversion and growth of nanorods and other materials. A method includes providing a condensed phase matrix material; and activating the condensed phase matrix material to produce a plurality of nanorods by condensed phase conversion and growth from the condensed phase matrix material instead of from vapor. The compositions are very strong. The compositions and methods provide advantages because they allow (1) formation rates of nanostructures necessary for reasonable production rates, and (2) the near net shaped production of component structures.

Combustion Of Metals In Reduced Gravity And Extraterrestrial Environments

NASA Technical Reports Server (NTRS)

Abbud-Madrid, A.; Modak, A.; Branch, M. C.

2003-01-01

The recent focus of this research project has been to model the combustion of isolated metal droplets and, in particular, to couple the existing theories and formulations of phenomena such as condensation, reaction kinetics, radiation, and surface reactions to formulate a more complete combustion model. A fully transient, one-dimensional (spherical symmetry) numerical model that uses detailed chemical kinetics, multi-component molecular transport mechanisms, condensation kinetics, and gas phase radiation heat transfer was developed. A coagulation model was used to simulate the particulate formation of MgO. The model was used to simulate the combustion of an Mg droplet in pure O2 and CO2. Methanol droplet combustion is considered as a test case for the solution method for both quasi-steady and fully transient simulations. Although some important processes unique to methanol combustion, such as water absorption at the surface, are not included in the model, the results are in sufficient agreement with the published data. Since the major part of the heat released in combustion of Mg, and in combustion of metals in general, is due to the condensation of the metal oxide, it is very important to capture the condensation processes correctly. Using the modified nucleation theory, an Arrhenius type rate expression is derived to calculate the condensation rate of MgO. This expression can be easily included in the CHEMKIN reaction mechanism format. Although very little property data is available for MgO, the condensation rate expression derived using the existing data is able to capture the condensation of MgO. An appropriate choice of the reference temperature to calculate the rate coefficients allows the model to correctly predict the subsequent heat release and hence the flame temperature.

Flavor symmetry breaking in lattice QCD with a mixed action

DOE Office of Scientific and Technical Information (OSTI.GOV)

Baer, Oliver; Golterman, Maarten; Shamir, Yigal

2011-03-01

We study the phase structure of mixed-action QCD with two Wilson sea quarks and any number of chiral valence quarks (and ghosts), starting from the chiral Lagrangian. A priori the effective theory allows for a rich phase structure, including a phase with a condensate made of sea and valence quarks. In such a phase, mass eigenstates would become admixtures of sea and valence fields, and pure-sea correlation functions would depend on the parameters of the valence sector, in contradiction with the actual setup of mixed-action simulations. Using that the spectrum of the chiral Dirac operator has a gap for nonzeromore » quark mass we prove that spontaneous symmetry breaking of the flavor symmetries can only occur within the sea sector. This rules out a mixed condensate and implies restrictions on the low-energy constants of the effective theory.« less

The influence of liquid/vapor phase change onto the Nusselt number

NASA Astrophysics Data System (ADS)

Popescu, Elena-Roxana; Colin, Catherine; Tanguy, Sebastien

2017-11-01

In spite of its significant interest in various fields, there is currently a very few information on how an external flow will modify the evaporation or the condensation of a liquid surface. Although most applications involve turbulent flows, the simpler configuration where a laminar superheated or subcooled vapor flow is shearing a saturated liquid interface has still never been solved. Based on a numerical approach, we propose to characterize the interaction between a laminar boundary layer of a superheated or subcooled vapor flow and a static liquid pool at saturation temperature. By performing a full set of simulations sweeping the parameters space, correlations are proposed for the first time on the Nusselt number depending on the dimensionless numbers that characterize both vaporization and condensation. As attended, the Nusselt number decreases or increases in the configurations involving respectively vaporization or condensation. More unexpected is the behaviour of the friction of the vapor flow on the liquid pool, for which we report that it is weakly affected by the phase change, despite the important variation of the local flow structure due to evaporation or condensation.

Multi-scale study of condensation in water jets using ellipsoidal-statistical Bhatnagar-Gross-Krook and molecular dynamics modeling

NASA Astrophysics Data System (ADS)

Li, Zheng; Borner, Arnaud; Levin, Deborah A.

2014-06-01

Homogeneous water condensation and ice formation in supersonic expansions to vacuum for stagnation pressures from 12 to 1000 mbar are studied using the particle-based Ellipsoidal-Statistical Bhatnagar-Gross-Krook (ES-BGK) method. We find that when condensation starts to occur, at a stagnation pressure of 96 mbar, the increase in the degree of condensation causes an increase in the rotational temperature due to the latent heat of vaporization. The simulated rotational temperature profiles along the plume expansion agree well with measurements confirming the kinetic homogeneous condensation models and the method of simulation. Comparisons of the simulated gas and cluster number densities, cluster size for different stagnation pressures along the plume centerline were made and it is found that the cluster size increase linearly with respect to stagnation pressure, consistent with classical nucleation theory. The sensitivity of our results to cluster nucleation model and latent heat values based on bulk water, specific cluster size, or bulk ice are examined. In particular, the ES-BGK simulations are found to be too coarse-grained to provide information on the phase or structure of the clusters formed. For this reason, molecular dynamics simulations of water condensation in a one-dimensional free expansion to simulate the conditions in the core of a plume are performed. We find that the internal structure of the clusters formed depends on the stagnation temperature. A larger cluster of average size 21 was tracked down the expansion, and a calculation of its average internal temperature as well as a comparison of its radial distribution functions (RDFs) with values measured for solid amorphous ice clusters lead us to conclude that this cluster is in a solid-like rather than liquid form. In another molecular-dynamics simulation at a much lower stagnation temperature, a larger cluster of size 324 and internal temperature 200 K was extracted from an expansion plume and equilibrated to determine its RDF and self-diffusion coefficient. The value of the latter shows that this cluster is formed in a supercooled liquid state rather than in an amorphous solid state.

Water vapour condensation in a partly closed structure. Comparison between results obtained with an inside wet or dry bottom wall

NASA Astrophysics Data System (ADS)

Batina, Jean; Peyrous, René

2018-04-01

We are interested in the determination of the more significant parameters acting on the water vapour condensation in a partly closed structure, submitted to external constraints (temperature and humidity), in view to recover the generated droplets as an additional source of potable water. External temperature variations, by inducing temperature differences between outside and inside of the structure, lead to convective movements and thermal variations inside this structure. Through an orifice, these movements permit a renewing of the humid inner air and can lead to the condensation of the water vapour initially contained in the inner air volume and/or on the walls. With the above hypotheses, and by using a numerical simulation [1] based on the ambient air characteristics and a finite volumes method, it appears that condensed water quantities are mainly depending on the boundary conditions imposed. These conditions are: 1) dimensions of the structure; 2) external temperature and relative hygrometry; 3) the phase φ (T/RH) linking thermal and hydrometric conditions; 4) the air renewing and its hygrometry for each phase; and 5) for each case, the fact that the inside bottom wall can be wet or dry. The resulting condensed water vapour quantities obtained, for the width section, point out clearly that they are very depending on this phase φ (T/RH) which appears as the more significant parameter and can be modified by the presence or not of a thin layer of water vapour on the inside bottom wall. Condensation phenomenon could be increased if φ could be optimized.

Initial decomposition of the condensed-phase β-HMX under shock waves: molecular dynamics simulations.

PubMed

Ge, Ni-Na; Wei, Yong-Kai; Ji, Guang-Fu; Chen, Xiang-Rong; Zhao, Feng; Wei, Dong-Qing

2012-11-26

We have performed quantum-based multiscale simulations to study the initial chemical processes of condensed-phase octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) under shock wave loading. A self-consistent charge density-functional tight-binding (SCC-DFTB) method was employed. The results show that the initial decomposition of shocked HMX is triggered by the N-NO(2) bond breaking under the low velocity impact (8 km/s). As the shock velocity increases (11 km/s), the homolytic cleavage of the N-NO(2) bond is suppressed under high pressure, the C-H bond dissociation becomes the primary pathway for HMX decomposition in its early stages. It is accompanied by a five-membered ring formation and hydrogen transfer from the CH(2) group to the -NO(2) group. Our simulations suggest that the initial chemical processes of shocked HMX are dependent on the impact velocity, which gain new insights into the initial decomposition mechanism of HMX upon shock loading at the atomistic level, and have important implications for understanding and development of energetic materials.

Cluster model studies of anion and molecular specificities via electrospray ionization photoelectron spectroscopy

DOE PAGES

Wang, Xue -Bin

2017-01-06

Ion specificity, a widely observed macroscopic phenomenon in condensed phases and at interfaces, is essentially a fundamental chemical physical issue. We have been investigating such effects using cluster models in an “atom-by-atom” and “molecule-by-molecule” fashion not possible with condensed-phase methods. We use electrospray ionization (ESI) to generate molecular and ionic clusters to simulate key molecular entities involved in local binding regions, and characterize them employing negative ion photoelectron spectroscopy (NIPES). Inter- and intramolecular interactions and binding configurations are directly obtained as functions of cluster size and composition, providing insightful molecular-level description and characterization over the local active sites that playmore » crucial roles in determining solution chemistry and condensed phase phenomena. Finally, the topics covered in this article are relevant to a wide scope of research fields ranging from ion specific effects in electrolyte solutions, ion selectivity/recognition in normal functioning of life, to molecular specificity in aerosol particle formation, as well as in rational material design and synthesis.« less

One-dimensional thermohydraulic code THESEUS and its application to chilldown process simulation in two-phase hydrogen flows

NASA Astrophysics Data System (ADS)

Papadimitriou, P.; Skorek, T.

THESUS is a thermohydraulic code for the calculation of steady state and transient processes of two-phase cryogenic flows. The physical model is based on four conservation equations with separate liquid and gas phase mass conservation equations. The thermohydraulic non-equilibrium is calculated by means of evaporation and condensation models. The mechanical non-equilibrium is modeled by a full-range drift-flux model. Also heat conduction in solid structures and heat exchange for the full spectrum of heat transfer regimes can be simulated. Test analyses of two-channel chilldown experiments and comparisons with the measured data have been performed.

Condensed phase conversion and growth of nanorods instead of from vapor

DOEpatents

Geohegan, David B.; Seals, Roland D.; Puretzky, Alex A.; Fan, Xudong

2005-08-02

Compositions, systems and methods are described for condensed phase conversion and growth of nanorods and other materials. A method includes providing a condensed phase matrix material; and activating the condensed phase matrix material to produce a plurality of nanorods by condensed phase conversion and growth from the condensed chase matrix material instead of from vacor. The compositions are very strong. The compositions and methods provide advantages because they allow (1) formation rates of nanostructures necessary for reasonable production rates, and (2) the near net shaped production of component structures.

A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations.

PubMed

Duan, Yong; Wu, Chun; Chowdhury, Shibasish; Lee, Mathew C; Xiong, Guoming; Zhang, Wei; Yang, Rong; Cieplak, Piotr; Luo, Ray; Lee, Taisung; Caldwell, James; Wang, Junmei; Kollman, Peter

2003-12-01

Molecular mechanics models have been applied extensively to study the dynamics of proteins and nucleic acids. Here we report the development of a third-generation point-charge all-atom force field for proteins. Following the earlier approach of Cornell et al., the charge set was obtained by fitting to the electrostatic potentials of dipeptides calculated using B3LYP/cc-pVTZ//HF/6-31G** quantum mechanical methods. The main-chain torsion parameters were obtained by fitting to the energy profiles of Ace-Ala-Nme and Ace-Gly-Nme di-peptides calculated using MP2/cc-pVTZ//HF/6-31G** quantum mechanical methods. All other parameters were taken from the existing AMBER data base. The major departure from previous force fields is that all quantum mechanical calculations were done in the condensed phase with continuum solvent models and an effective dielectric constant of epsilon = 4. We anticipate that this force field parameter set will address certain critical short comings of previous force fields in condensed-phase simulations of proteins. Initial tests on peptides demonstrated a high-degree of similarity between the calculated and the statistically measured Ramanchandran maps for both Ace-Gly-Nme and Ace-Ala-Nme di-peptides. Some highlights of our results include (1) well-preserved balance between the extended and helical region distributions, and (2) favorable type-II poly-proline helical region in agreement with recent experiments. Backward compatibility between the new and Cornell et al. charge sets, as judged by overall agreement between dipole moments, allows a smooth transition to the new force field in the area of ligand-binding calculations. Test simulations on a large set of proteins are also discussed. Copyright 2003 Wiley Periodicals, Inc. J Comput Chem 24: 1999-2012, 2003

Aerosol-cloud interactions in mixed-phase convective clouds - Part 1: Aerosol perturbations

NASA Astrophysics Data System (ADS)

Miltenberger, Annette K.; Field, Paul R.; Hill, Adrian A.; Rosenberg, Phil; Shipway, Ben J.; Wilkinson, Jonathan M.; Scovell, Robert; Blyth, Alan M.

2018-03-01

Changes induced by perturbed aerosol conditions in moderately deep mixed-phase convective clouds (cloud top height ˜ 5 km) developing along sea-breeze convergence lines are investigated with high-resolution numerical model simulations. The simulations utilise the newly developed Cloud-AeroSol Interacting Microphysics (CASIM) module for the Unified Model (UM), which allows for the representation of the two-way interaction between cloud and aerosol fields. Simulations are evaluated against observations collected during the COnvective Precipitation Experiment (COPE) field campaign over the southwestern peninsula of the UK in 2013. The simulations compare favourably with observed thermodynamic profiles, cloud base cloud droplet number concentrations (CDNC), cloud depth, and radar reflectivity statistics. Including the modification of aerosol fields by cloud microphysical processes improves the correspondence with observed CDNC values and spatial variability, but reduces the agreement with observations for average cloud size and cloud top height. Accumulated precipitation is suppressed for higher-aerosol conditions before clouds become organised along the sea-breeze convergence lines. Changes in precipitation are smaller in simulations with aerosol processing. The precipitation suppression is due to less efficient precipitation production by warm-phase microphysics, consistent with parcel model predictions. In contrast, after convective cells organise along the sea-breeze convergence zone, accumulated precipitation increases with aerosol concentrations. Condensate production increases with the aerosol concentrations due to higher vertical velocities in the convective cores and higher cloud top heights. However, for the highest-aerosol scenarios, no further increase in the condensate production occurs, as clouds grow into an upper-level stable layer. In these cases, the reduced precipitation efficiency (PE) dominates the precipitation response and no further precipitation enhancement occurs. Previous studies of deep convective clouds have related larger vertical velocities under high-aerosol conditions to enhanced latent heating from freezing. In the presented simulations changes in latent heating above the 0°C are negligible, but latent heating from condensation increases with aerosol concentrations. It is hypothesised that this increase is related to changes in the cloud field structure reducing the mixing of environmental air into the convective core. The precipitation response of the deeper mixed-phase clouds along well-established convergence lines can be the opposite of predictions from parcel models. This occurs when clouds interact with a pre-existing thermodynamic environment and cloud field structural changes occur that are not captured by simple parcel model approaches.

Sensitivity of Aerosol Mass and Microphysics to varying treatments of Condensational Growth of Secondary Organic Compounds in a regional model

NASA Astrophysics Data System (ADS)

Lowe, Douglas; Topping, David; McFiggans, Gordon

2017-04-01

Gas to particle partitioning of atmospheric compounds occurs through disequilibrium mass transfer rather than through instantaneous equilibrium. However, it is common to treat only the inorganic compounds as partitioning dynamically whilst organic compounds, represented by the Volatility Basis Set (VBS), are partitioned instantaneously. In this study we implement a more realistic dynamic partitioning of organic compounds in a regional framework and assess impact on aerosol mass and microphysics. It is also common to assume condensed phase water is only associated with inorganic components. We thus also assess sensitivity to assuming all organics are hygroscopic according to their prescribed molecular weight. For this study we use WRF-Chem v3.4.1, focusing on anthropogenic dominated North-Western Europe. Gas-phase chemistry is represented using CBM-Z whilst aerosol dynamics are simulated using the 8-section MOSAIC scheme, including a 9-bin VBS treatment of organic aerosol. Results indicate that predicted mass loadings can vary significantly. Without gas phase ageing of higher volatility compounds, dynamic partitioning always results in lower mass loadings downwind of emission sources. The inclusion of condensed phase water in both partitioning models increases the predicted PM mass, resulting from a larger contribution from higher volatility organics, if present. If gas phase ageing of VBS compounds is allowed to occur in a dynamic model, this can often lead to higher predicted mass loadings, contrary to expected behaviour from a simple non-reactive gas phase box model. As descriptions of aerosol phase processes improve within regional models, the baseline descriptions of partitioning should retain the ability to treat dynamic partitioning of organics compounds. Using our simulations, we discuss whether derived sensitivities to aerosol processes in existing models may be inherently biased. This work was supported by the Natural Environment Research Council within the RONOCO (NE/F004656/1) and CCN-Vol (NE/L007827/1) projects.

Sensitivity of Aerosol Mass and Microphysics to Treatments of Condensational Growth of Secondary Organic Compounds in a Regional Model

NASA Astrophysics Data System (ADS)

Topping, D. O.; Lowe, D.; McFiggans, G.; Zaveri, R. A.

2016-12-01

Gas to particle partitioning of atmospheric compounds occurs through disequilibrium mass transfer rather than through instantaneous equilibrium. However, it is common to treat only the inorganic compounds as partitioning dynamically whilst organic compounds, represented by the Volatility Basis Set (VBS), are partitioned instantaneously. In this study we implement a more realistic dynamic partitioning of organic compounds in a regional framework and assess impact on aerosol mass and microphysics. It is also common to assume condensed phase water is only associated with inorganic components. We thus also assess sensitivity to assuming all organics are hygroscopic according to their prescribed molecular weight.For this study we use WRF-Chem v3.4.1, focusing on anthropogenic dominated North-Western Europe. Gas-phase chemistry is represented using CBM-Z whilst aerosol dynamics are simulated using the 8-section MOSAIC scheme, including a 9-bin volatility basis set (VBS) treatment of organic aerosol. Results indicate that predicted mass loadings can vary significantly. Without gas phase ageing of higher volatility compounds, dynamic partitioning always results in lower mass loadings downwind of emission sources. The inclusion of condensed phase water in both partitioning models increases the predicted PM mass, resulting from a larger contribution from higher volatility organics, if present. If gas phase ageing of VBS compounds is allowed to occur in a dynamic model, this can often lead to higher predicted mass loadings, contrary to expected behaviour from a simple non-reactive gas phase box model. As descriptions of aerosol phase processes improve within regional models, the baseline descriptions of partitioning should retain the ability to treat dynamic partitioning of organic compounds. Using our simulations, we discuss whether derived sensitivities to aerosol processes in existing models may be inherently biased.This work was supported by the Nature Environment Research Council within the RONOCO (NE/F004656/1) and CCN-Vol (NE/L007827/1) projects.

Development of a wet vapor homogeneous liquid metal MHD power system

NASA Astrophysics Data System (ADS)

1989-04-01

During the period covered by this report (October 1988 to March 1989), the following work was done: the mixing stream condensation process was analyzed, and a theoretical model for simulating this process was modified. A parametric study is being conducted at the present time; the separation processes were analyzed; and the experimental system was specified and its design is at present in an advanced stage. The mixing stream condensation process was analyzed. For the parameters defined in the SOW of this project the process was found to be a mist flow direct contact condensation, where the hot gas mixture consisting of inert gas and vapor is the continuous phase, and the subcooled liquid on which the vapor is condensed if the droplets dispersed phase. Two possibilities of creating the mist flow were considered. The first, injecting the cold Liquid Metal (LM) into the Mixing Streams Condenser (MSC) entrance as a jet and breaking it into LM fragments and the fragments into droplets by momentum transfer breakup mechanism. The second, atomizing the cooled LM stream into little droplets (approximately 100 micrometers in diameter) and accelerating them by the gas. The second possibility was preferred due to its much higher heat and mass transfer surface and coefficients relative to the first one.

Numerical simulation of steady cavitating flow of viscous fluid in a Francis hydroturbine

NASA Astrophysics Data System (ADS)

Panov, L. V.; Chirkov, D. V.; Cherny, S. G.; Pylev, I. M.; Sotnikov, A. A.

2012-09-01

Numerical technique was developed for simulation of cavitating flows through the flow passage of a hydraulic turbine. The technique is based on solution of steady 3D Navier—Stokes equations with a liquid phase transfer equation. The approch for setting boundary conditions meeting the requirements of cavitation testing standard was suggested. Four different models of evaporation and condensation were compared. Numerical simulations for turbines of different specific speed were compared with experiment.

Numerical investigation of transient behaviour of the recuperative heat exchanger in a MR J-T cryocooler using different heat transfer correlations

NASA Astrophysics Data System (ADS)

Damle, R. M.; Ardhapurkar, P. M.; Atrey, M. D.

2016-12-01

In J-T cryocoolers operating with mixed refrigerants (nitrogen-hydrocarbons), the recuperative heat exchange takes place under two-phase conditions. Simultaneous boiling of the low pressure stream and condensation of the high pressure stream results in higher heat transfer coefficients. The mixture composition, operating conditions and the heat exchanger design are crucial for obtaining the required cryogenic temperature. In this work, a one-dimensional transient algorithm is developed for the simulation of the two-phase heat transfer in the recuperative heat exchanger of a mixed refrigerant J-T cryocooler. Modified correlation is used for flow boiling of the high pressure fluid while different condensation correlations are employed with and without the correction for the low pressure fluid. Simulations are carried out for different mixture compositions and numerical predictions are compared with the experimental data. The overall heat transfer is predicted reasonably well and the qualitative trends of the temperature profiles are also captured by the developed numerical model.

Effects of defects on thermal decomposition of HMX via ReaxFF molecular dynamics simulations.

PubMed

Zhou, Ting-Ting; Huang, Feng-Lei

2011-01-20

Effects of molecular vacancies on the decomposition mechanisms and reaction dynamics of condensed-phase β-HMX at various temperatures were studied using ReaxFF molecular dynamics simulations. Results show that three primary initial decomposition mechanisms, namely, N-NO(2) bond dissociation, HONO elimination, and concerted ring fission, exist at both high and lower temperatures. The contribution of the three mechanisms to the initial decomposition of HMX is influenced by molecular vacancies, and the effects vary with temperature. At high temperature (2500 K), molecular vacancies remarkably promote N-N bond cleavage and concerted ring breaking but hinder HONO formation. N-N bond dissociation and HONO elimination are two primary competing reaction mechanisms, and the former is dominant in the initial decomposition. Concerted ring breaking of condensed-phase HMX is not favored at high temperature. At lower temperature (1500 K), the most preferential initial decomposition pathway is N-N bond dissociation followed by the formation of NO(3) (O migration), although all three mechanisms are promoted by molecular vacancies. The promotion effect on concerted ring breaking is considerable at lower temperature. Products resulting from concerted ring breaking appear in the defective system but not in the perfect crystal. The mechanism of HONO elimination is less important at lower temperature. We also estimated the reaction rate constant and activation barriers of initial decomposition with different vacancy concentrations. Molecular vacancies accelerate the decomposition of condensed-phase HMX by increasing the reaction rate constant and reducing activation barriers.

Particle growth kinetics over the Amazon rainforest

NASA Astrophysics Data System (ADS)

Pinterich, T.; Andreae, M. O.; Artaxo, P.; Kuang, C.; Longo, K.; Machado, L.; Manzi, A. O.; Martin, S. T.; Mei, F.; Pöhlker, C.; Pöhlker, M. L.; Poeschl, U.; Shilling, J. E.; Shiraiwa, M.; Tomlinson, J. M.; Zaveri, R. A.; Wang, J.

2016-12-01

Aerosol particles larger than 100 nm play a key role in global climate by acting as cloud condensation nuclei (CCN). Most of these particles, originated from new particle formation or directly emitted into the atmospheric, are initially too small to serve as CCN. These small particles grow to CCN size mainly through condensation of secondary species. In one extreme, the growth is dictated by kinetic condensation of very low-volatility compounds, favoring the growth of the smallest particles; in the other extreme, the process is driven by Raoult's law-based equilibrium partitioning of semi-volatile organic compound, favoring the growth of larger particles. These two mechanisms can lead to very different production rates of CCN. The growth of particles depends on a number of parameters, including the volatility of condensing species, particle phase, and diffusivity inside the particles, and this process is not well understood in part due to lack of ambient data. Here we examine atmospheric particle growth using high-resolution size distributions measured onboard the DOE G-1 aircraft during GoAmazon campaign, which took place from January 2014 to December 2015 near Manaus, Brazil, a city surrounded by natural forest for over 1000 km in every direction. City plumes are clearly identified by the strong enhancement of nucleation and Aitken mode particle concentrations over the clean background. As the plume traveled downwind, particle growth was observed, and is attributed to condensation of secondary species and coagulation (Fig.1). Observed aerosol growth is modeled using MOSAIC (Model for Simulating Aerosol Interactions and Chemistry), which dynamically partitions multiple compounds to all particle size bins by taking into account compound volatility, gas-phase diffusion, interfacial mass accommodation, particle-phase diffusion, and particle-phase reaction. The results from both wet and dry seasons will be discussed.

Black holes as critical point of quantum phase transition.

PubMed

Dvali, Gia; Gomez, Cesar

We reformulate the quantum black hole portrait in the language of modern condensed matter physics. We show that black holes can be understood as a graviton Bose-Einstein condensate at the critical point of a quantum phase transition, identical to what has been observed in systems of cold atoms. The Bogoliubov modes that become degenerate and nearly gapless at this point are the holographic quantum degrees of freedom responsible for the black hole entropy and the information storage. They have no (semi)classical counterparts and become inaccessible in this limit. These findings indicate a deep connection between the seemingly remote systems and suggest a new quantum foundation of holography. They also open an intriguing possibility of simulating black hole information processing in table-top labs.

Quantum Monte Carlo Simulation of condensed van der Waals Systems

NASA Astrophysics Data System (ADS)

Benali, Anouar; Shulenburger, Luke; Romero, Nichols A.; Kim, Jeongnim; Anatole von Lilienfeld, O.

2012-02-01

Van der Waals forces are as ubiquitous as infamous. While post-Hartree-Fock methods enable accurate estimates of these forces in molecules and clusters, they remain elusive for dealing with many-electron condensed phase systems. We present Quantum Monte Carlo [1,2] results for condensed van der Waals systems. Interatomic many-body contributions to cohesive energies and bulk modulus will be discussed. Numerical evidence is presented for crystals of rare gas atoms, and compared to experiments and methods [3]. Sandia National Laboratories is a multiprogram laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. DoE's National Nuclear Security Administration under Contract No. DE-AC04-94AL85000.[4pt] [1] J. Kim, K. Esler, J. McMinis and D. Ceperley, SciDAC 2010, J. of Physics: Conference series, Chattanooga, Tennessee, July 11 2011 [0pt] [2] QMCPACK simulation suite, http://qmcpack.cmscc.org (unpublished)[0pt] [3] O. A. von Lillienfeld and A. Tkatchenko, J. Chem. Phys. 132 234109 (2010)

Synthesis of condensed phases containing polycyclic aromatic hydrocarbons fullerenes and nanotubes

DOEpatents

Reilly, Peter T. A.

2004-10-19

The invention relates to methods for producing polycyclic aromatic hydrocarbons, fullerenes, and nanotubes, comprising: a. heating at least one carbon-containing material to form a condensed phase comprising at least one polycyclic aromatic hydrocarbon; b. collecting at least some of the condensed phase; c. reacting the condensed phase to form fullerenes and/or nanotubes.

Zero point energy leakage in condensed phase dynamics: An assessment of quantum simulation methods for liquid water

NASA Astrophysics Data System (ADS)

Habershon, Scott; Manolopoulos, David E.

2009-12-01

The approximate quantum mechanical ring polymer molecular dynamics (RPMD) and linearized semiclassical initial value representation (LSC-IVR) methods are compared and contrasted in a study of the dynamics of the flexible q-TIP4P/F water model at room temperature. For this water model, a RPMD simulation gives a diffusion coefficient that is only a few percent larger than the classical diffusion coefficient, whereas a LSC-IVR simulation gives a diffusion coefficient that is three times larger. We attribute this discrepancy to the unphysical leakage of initially quantized zero point energy (ZPE) from the intramolecular to the intermolecular modes of the liquid as the LSC-IVR simulation progresses. In spite of this problem, which is avoided by construction in RPMD, the LSC-IVR may still provide a useful approximation to certain short-time dynamical properties which are not so strongly affected by the ZPE leakage. We illustrate this with an application to the liquid water dipole absorption spectrum, for which the RPMD approximation breaks down at frequencies in the O-H stretching region owing to contamination from the internal modes of the ring polymer. The LSC-IVR does not suffer from this difficulty and it appears to provide quite a promising way to calculate condensed phase vibrational spectra.

Zero point energy leakage in condensed phase dynamics: an assessment of quantum simulation methods for liquid water.

PubMed

Habershon, Scott; Manolopoulos, David E

2009-12-28

The approximate quantum mechanical ring polymer molecular dynamics (RPMD) and linearized semiclassical initial value representation (LSC-IVR) methods are compared and contrasted in a study of the dynamics of the flexible q-TIP4P/F water model at room temperature. For this water model, a RPMD simulation gives a diffusion coefficient that is only a few percent larger than the classical diffusion coefficient, whereas a LSC-IVR simulation gives a diffusion coefficient that is three times larger. We attribute this discrepancy to the unphysical leakage of initially quantized zero point energy (ZPE) from the intramolecular to the intermolecular modes of the liquid as the LSC-IVR simulation progresses. In spite of this problem, which is avoided by construction in RPMD, the LSC-IVR may still provide a useful approximation to certain short-time dynamical properties which are not so strongly affected by the ZPE leakage. We illustrate this with an application to the liquid water dipole absorption spectrum, for which the RPMD approximation breaks down at frequencies in the O-H stretching region owing to contamination from the internal modes of the ring polymer. The LSC-IVR does not suffer from this difficulty and it appears to provide quite a promising way to calculate condensed phase vibrational spectra.

Contrail Formation in Aircraft Wakes Using Large-Eddy Simulations

NASA Technical Reports Server (NTRS)

Paoli, R.; Helie, J.; Poinsot, T. J.; Ghosal, S.

2002-01-01

In this work we analyze the issue of the formation of condensation trails ("contrails") in the near-field of an aircraft wake. The basic configuration consists in an exhaust engine jet interacting with a wing-tip training vortex. The procedure adopted relies on a mixed Eulerian/Lagrangian two-phase flow approach; a simple micro-physics model for ice growth has been used to couple ice and vapor phases. Large eddy simulations have carried out at a realistic flight Reynolds number to evaluate the effects of turbulent mixing and wake vortex dynamics on ice-growth characteristics and vapor thermodynamic properties.

Dynamics in multiple-well Bose-Einstein condensates

NASA Astrophysics Data System (ADS)

Nigro, M.; Capuzzi, P.; Cataldo, H. M.; Jezek, D. M.

2018-01-01

We study the dynamics of three-dimensional weakly linked Bose-Einstein condensates using a multimode model with an effective interaction parameter. The system is confined by a ring-shaped four-well trapping potential. By constructing a two-mode Hamiltonian in a reduced highly symmetric phase space, we examine the periodic orbits and calculate their time periods both in the self-trapping and Josephson regimes. The dynamics in the vicinity of the reduced phase space is investigated by means of a Floquet multiplier analysis, finding regions of different linear stability and analyzing their implications on the exact dynamics. The numerical exploration in an extended region of the phase space demonstrates that two-mode tools can also be useful for performing a partition of the space in different regimes. Comparisons with Gross-Pitaevskii simulations confirm these findings and emphasize the importance of properly determining the effective on-site interaction parameter governing the multimode dynamics.

Numerical simulation of cloud and precipitation structure during GALE IOP-2

NASA Technical Reports Server (NTRS)

Robertson, F. R.; Perkey, D. J.; Seablom, M. S.

1988-01-01

A regional scale model, LAMPS (Limited Area Mesoscale Prediction System), is used to investigate cloud and precipitation structure that accompanied a short wave system during a portion of GALE IOP-2. A comparison of satellite imagery and model fields indicates that much of the large mesoscale organization of condensation has been captured by the simulation. In addition to reproducing a realistic phasing of two baroclinic zones associated with a split cold front, a reasonable simulation of the gross mesoscale cloud distribution has been achieved.

Simulating a topological transition in a superconducting phase qubit by fast adiabatic trajectories

NASA Astrophysics Data System (ADS)

Wang, Tenghui; Zhang, Zhenxing; Xiang, Liang; Gong, Zhihao; Wu, Jianlan; Yin, Yi

2018-04-01

The significance of topological phases has been widely recognized in the community of condensed matter physics. The well controllable quantum systems provide an artificial platform to probe and engineer various topological phases. The adiabatic trajectory of a quantum state describes the change of the bulk Bloch eigenstates with the momentum, and this adiabatic simulation method is however practically limited due to quantum dissipation. Here we apply the "shortcut to adiabaticity" (STA) protocol to realize fast adiabatic evolutions in the system of a superconducting phase qubit. The resulting fast adiabatic trajectories illustrate the change of the bulk Bloch eigenstates in the Su-Schrieffer-Heeger (SSH) model. A sharp transition is experimentally determined for the topological invariant of a winding number. Our experiment helps identify the topological Chern number of a two-dimensional toy model, suggesting the applicability of the fast adiabatic simulation method for topological systems.

Antiferromagnetic spinor condensates in a bichromatic superlattice

NASA Astrophysics Data System (ADS)

Tang, Tao; Zhao, Lichao; Chen, Zihe; Liu, Yingmei

2017-04-01

A spinor Bose-Einstein condensate in an optical supelattice has been considered as a good quantum simulator for understanding mesoscopic magnetism. We report an experimental study on an antiferromagnetic spinor condensate in a bichromatic superlattice constructed by a cubic red-detuned optical lattice and a one-dimensional blue-detuned optical lattice. Our data demonstrate a few advantages of this bichromatic superlattice over a monochromatic lattice. One distinct advantage is that the bichromatic superlattice enables realizing the first-order superfluid to Mott-insulator phase transitions within a much wider range of magnetic fields. In addition, we discuss an apparent discrepancy between our data and the mean-field theory. We thank the National Science Foundation and the Oklahoma Center for the Advancement of Science and Technology for financial support.

Investigation of Kibble-Zurek Quench Dynamics in a Spin-1 Ferromagnetic BEC

NASA Astrophysics Data System (ADS)

Anquez, Martin; Robbins, Bryce; Hoang, Thai; Yang, Xiaoyun; Land, Benjamin; Hamley, Christopher; Chapman, Michael

2014-05-01

We study the temporal evolution of spin populations in small spin-1 87Rb condensates following a slow quench. A ferromagnetic spin-1 BEC exhibits a second-order gapless (quantum) phase transition due to a competition between the magnetic and collisional spin interaction energies. The dynamics of slow quenches through the critical point are predicted to exhibit universal power-law scaling as a function of quench speed. In spatially extended condensates, these excitations are revealed as spatial spin domains. In small condensates, the excitations are manifest in the temporal evolution of the spin populations, illustrating a Kibble-Zurek type scaling. We will present the results of our investigation and compare them to full quantum simulations of the system.

Analytical solutions for coagulation and condensation kinetics of composite particles

NASA Astrophysics Data System (ADS)

Piskunov, Vladimir N.

2013-04-01

The processes of composite particles formation consisting of a mixture of different materials are essential for many practical problems: for analysis of the consequences of accidental releases in atmosphere; for simulation of precipitation formation in clouds; for description of multi-phase processes in chemical reactors and industrial facilities. Computer codes developed for numerical simulation of these processes require optimization of computational methods and verification of numerical programs. Kinetic equations of composite particle formation are given in this work in a concise form (impurity integrated). Coagulation, condensation and external sources associated with nucleation are taken into account. Analytical solutions were obtained in a number of model cases. The general laws for fraction redistribution of impurities were defined. The results can be applied to develop numerical algorithms considerably reducing the simulation effort, as well as to verify the numerical programs for calculation of the formation kinetics of composite particles in the problems of practical importance.

Floquet engineering of Haldane Chern insulators and chiral bosonic phase transitions

NASA Astrophysics Data System (ADS)

Plekhanov, Kirill; Roux, Guillaume; Le Hur, Karyn

2017-01-01

The realization of synthetic gauge fields has attracted a lot of attention recently in relation to periodically driven systems and the Floquet theory. In ultracold atom systems in optical lattices and photonic networks, this allows one to simulate exotic phases of matter such as quantum Hall phases, anomalous quantum Hall phases, and analogs of topological insulators. In this paper, we apply the Floquet theory to engineer anisotropic Haldane models on the honeycomb lattice and two-leg ladder systems. We show that these anisotropic Haldane models still possess a topologically nontrivial band structure associated with chiral edge modes. Focusing on (interacting) boson systems in s -wave bands of the lattice, we show how to engineer through the Floquet theory, a quantum phase transition (QPT) between a uniform superfluid and a Bose-Einstein condensate analog of Fulde-Ferrell-Larkin-Ovchinnikov states, where bosons condense at nonzero wave vectors. We perform a Ginzburg-Landau analysis of the QPT on the graphene lattice, and compute observables such as chiral currents and the momentum distribution. The results are supported by exact diagonalization calculations and compared with those of the isotropic situation. The validity of high-frequency expansion in the Floquet theory is also tested using time-dependent simulations for various parameters of the model. Last, we show that the anisotropic choice for the effective vector potential allows a bosonization approach in equivalent ladder (strip) geometries.

Water interactions with condensed organic phases: a combined experimental and theoretical study of molecular-level processes

NASA Astrophysics Data System (ADS)

Johansson, Sofia M.; Kong, Xiangrui; Thomson, Erik S.; Papagiannakopoulos, Panos; Pettersson, Jan B. C.; Lovrić, Josip; Toubin, Céline

2016-04-01

Water uptake on aerosol particles modifies their chemistry and microphysics with important implications for air quality and climate. A large fraction of the atmospheric aerosol consists of organic aerosol particles or inorganic particles with condensed organic components. Here, we combine laboratory studies using the environmental molecular beam (EMB) method1 with molecular dynamics (MD) simulations to characterize water interactions with organic surfaces in detail. The over-arching aim is to characterize the mechanisms that govern water uptake, in order to guide the development of physics-based models to be used in atmospheric modelling. The EMB method enables molecular level studies of interactions between gases and volatile surfaces at near ambient pressure,1 and the technique may provide information about collision dynamics, surface and bulk accommodation, desorption and diffusion kinetics. Molecular dynamics simulations provide complementary information about the collision dynamics and initial interactions between gas molecules and the condensed phase. Here, we focus on water interactions with condensed alcohol phases that serve as highly simplified proxies for systems in the environment. Gas-surface collisions are in general found to be highly inelastic and result in efficient surface accommodation of water molecules. As a consequence, surface accommodation of water can be safely assumed to be close to unity under typical ambient conditions. Bulk accommodation is inefficient on solid alcohol and the condensed materials appear to produce hydrophobic surface structures, with limited opportunities for adsorbed water to form hydrogen bonds with surface molecules. Accommodation is significantly more efficient on the dynamic liquid alcohol surfaces. The results for n-butanol (BuOH) are particularly intriguing where substantial changes in water accommodation taking place over a 10 K interval below and above the BuOH melting point.2 The governing mechanisms for the observed water accommodation are discussed based on the combined EMB and MD results. The studies illustrate that the detailed surface properties of the condensed organic phase may substantially modify water uptake, with potential implications for the properties and action of aerosols and clouds in the Earth system. References: 1. X.R. Kong, E. S. Thomson, P. Papagiannakopoulos, S.M. Johansson, and J.B.C. Pettersson, Water Accommodation on Ice and Organic Surfaces: Insights from Environmental Molecular Beam Experiments. J. Phys. Chem. B 118 (2014) 13378-13386. 2. P. Papagiannakopoulos, X. Kong, E. S. Thomson, N. Marković, and J. B. C. Pettersson, Surface Transformations and Water Uptake on Liquid and Solid Butanol near the Melting Temperature. J. Phys. Chem. C 117 (2013) 6678-6685.

Atomistic Simulations of Detonation Instabilities in Condensed Phase Systems

NASA Astrophysics Data System (ADS)

Kober, Edward; Heim, Andrew; Germann, Timothy; Jensen, Niels

2007-06-01

We report the results of simulations of condensed phase detonation phenomena using a model diatomic system: 2AB -> A2 + B2. The initial set of parameters for this system corresponded to the Model 0 set of C. White et al, which exhibits a steady, Chapman-Jouget (CJ) detonation structure with a reaction zone length of 30-100 å. This has a highly compressed CJ state (V/V0˜0.5) that does not consist of discrete molecular species. The potential form was modified so that a more molecular CJ state resulted, consistent with the models for conventional organic explosives. The new system has a less dense CJ state (V/V0˜0.8), and the reaction zone was substantially extended. The reaction rate fits Arrhenius-type kinetics with an activation energy of ˜2 eV, with a minor density dependence. In contrast, the original Model 0 system had a lower activation energy (˜1 eV) with a stronger density dependence. The new system exhibits quite marked two dimensional instability structures with well-defined wavelengths similar to what has been observed for gas-phase detonations and for nitromethane. Depending on the exothermicity and the width of the periodic simulations, these instabilities can result in either detonation failure or quasi-steady propagation. The observed propagation velocities are several per cent higher than CJ values derived from thermodynamic analyses.

Electron emission from condensed phase material induced by fast protons.

PubMed

Shinpaugh, J L; McLawhorn, R A; McLawhorn, S L; Carnes, K D; Dingfelder, M; Travia, A; Toburen, L H

2011-02-01

Monte Carlo track simulation has become an important tool in radiobiology. Monte Carlo transport codes commonly rely on elastic and inelastic electron scattering cross sections determined using theoretical methods supplemented with gas-phase data; experimental condensed phase data are often unavailable or infeasible. The largest uncertainties in the theoretical methods exist for low-energy electrons, which are important for simulating electron track ends. To test the reliability of these codes to deal with low-energy electron transport, yields of low-energy secondary electrons ejected from thin foils have been measured following passage of fast protons. Fast ions, where interaction cross sections are well known, provide the initial spectrum of low-energy electrons that subsequently undergo elastic and inelastic scattering in the material before exiting the foil surface and being detected. These data, measured as a function of the energy and angle of the emerging electrons, can provide tests of the physics of electron transport. Initial measurements from amorphous solid water frozen to a copper substrate indicated substantial disagreement with MC simulation, although questions remained because of target charging. More recent studies, using different freezing techniques, do not exhibit charging, but confirm the disagreement seen earlier between theory and experiment. One now has additional data on the absolute differential electron yields from copper, aluminum and gold, as well as for thin films of frozen hydrocarbons. Representative data are presented.

Electron emission from condensed phase material induced by fast protons†

PubMed Central

Shinpaugh, J. L.; McLawhorn, R. A.; McLawhorn, S. L.; Carnes, K. D.; Dingfelder, M.; Travia, A.; Toburen, L. H.

2011-01-01

Monte Carlo track simulation has become an important tool in radiobiology. Monte Carlo transport codes commonly rely on elastic and inelastic electron scattering cross sections determined using theoretical methods supplemented with gas-phase data; experimental condensed phase data are often unavailable or infeasible. The largest uncertainties in the theoretical methods exist for low-energy electrons, which are important for simulating electron track ends. To test the reliability of these codes to deal with low-energy electron transport, yields of low-energy secondary electrons ejected from thin foils have been measured following passage of fast protons. Fast ions, where interaction cross sections are well known, provide the initial spectrum of low-energy electrons that subsequently undergo elastic and inelastic scattering in the material before exiting the foil surface and being detected. These data, measured as a function of the energy and angle of the emerging electrons, can provide tests of the physics of electron transport. Initial measurements from amorphous solid water frozen to a copper substrate indicated substantial disagreement with MC simulation, although questions remained because of target charging. More recent studies, using different freezing techniques, do not exhibit charging, but confirm the disagreement seen earlier between theory and experiment. One now has additional data on the absolute differential electron yields from copper, aluminum and gold, as well as for thin films of frozen hydrocarbons. Representative data are presented. PMID:21183539

Mass accommodation of water: bridging the gap between molecular dynamics simulations and kinetic condensation models.

PubMed

Julin, Jan; Shiraiwa, Manabu; Miles, Rachael E H; Reid, Jonathan P; Pöschl, Ulrich; Riipinen, Ilona

2013-01-17

The condensational growth of submicrometer aerosol particles to climate relevant sizes is sensitive to their ability to accommodate vapor molecules, which is described by the mass accommodation coefficient. However, the underlying processes are not yet fully understood. We have simulated the mass accommodation and evaporation processes of water using molecular dynamics, and the results are compared to the condensation equations derived from the kinetic gas theory to shed light on the compatibility of the two. Molecular dynamics simulations were performed for a planar TIP4P-Ew water surface at four temperatures in the range 268-300 K as well as two droplets, with radii of 1.92 and 4.14 nm at T = 273.15 K. The evaporation flux from molecular dynamics was found to be in good qualitative agreement with that predicted by the simple kinetic condensation equations. Water droplet growth was also modeled with the kinetic multilayer model KM-GAP of Shiraiwa et al. [Atmos. Chem. Phys. 2012, 12, 2777]. It was found that, due to the fast transport across the interface, the growth of a pure water droplet is controlled by gas phase diffusion. These facts indicate that the simple kinetic treatment is sufficient in describing pure water condensation and evaporation. The droplet size was found to have minimal effect on the value of the mass accommodation coefficient. The mass accommodation coefficient was found to be unity (within 0.004) for all studied surfaces, which is in agreement with previous simulation work. Additionally, the simulated evaporation fluxes imply that the evaporation coefficient is also unity. Comparing the evaporation rates of the mass accommodation and evaporation simulations indicated that the high collision flux, corresponding to high supersaturation, present in typical molecular dynamics mass accommodation simulations can under certain conditions lead to an increase in the evaporation rate. Consequently, in such situations the mass accommodation coefficient can be overestimated, but in the present cases the corrected values were still close to unity with the lowest value at ≈0.99.

Mass Accommodation of Water: Bridging the Gap Between Molecular Dynamics Simulations and Kinetic Condensation Models

PubMed Central

2012-01-01

The condensational growth of submicrometer aerosol particles to climate relevant sizes is sensitive to their ability to accommodate vapor molecules, which is described by the mass accommodation coefficient. However, the underlying processes are not yet fully understood. We have simulated the mass accommodation and evaporation processes of water using molecular dynamics, and the results are compared to the condensation equations derived from the kinetic gas theory to shed light on the compatibility of the two. Molecular dynamics simulations were performed for a planar TIP4P-Ew water surface at four temperatures in the range 268–300 K as well as two droplets, with radii of 1.92 and 4.14 nm at T = 273.15 K. The evaporation flux from molecular dynamics was found to be in good qualitative agreement with that predicted by the simple kinetic condensation equations. Water droplet growth was also modeled with the kinetic multilayer model KM-GAP of Shiraiwa et al. [Atmos. Chem. Phys.2012, 117, 2777]. It was found that, due to the fast transport across the interface, the growth of a pure water droplet is controlled by gas phase diffusion. These facts indicate that the simple kinetic treatment is sufficient in describing pure water condensation and evaporation. The droplet size was found to have minimal effect on the value of the mass accommodation coefficient. The mass accommodation coefficient was found to be unity (within 0.004) for all studied surfaces, which is in agreement with previous simulation work. Additionally, the simulated evaporation fluxes imply that the evaporation coefficient is also unity. Comparing the evaporation rates of the mass accommodation and evaporation simulations indicated that the high collision flux, corresponding to high supersaturation, present in typical molecular dynamics mass accommodation simulations can under certain conditions lead to an increase in the evaporation rate. Consequently, in such situations the mass accommodation coefficient can be overestimated, but in the present cases the corrected values were still close to unity with the lowest value at ≈0.99. PMID:23253100

Phase transition kinetics for a Bose Einstein condensate in a periodically driven band system

NASA Astrophysics Data System (ADS)

Michon, E.; Cabrera-Gutiérrez, C.; Fortun, A.; Berger, M.; Arnal, M.; Brunaud, V.; Billy, J.; Petitjean, C.; Schlagheck, P.; Guéry-Odelin, D.

2018-05-01

The dynamical transition of an atomic Bose–Einstein condensate from a spatially periodic state to a staggered state with alternating sign in its wavefunction is experimentally studied using a one-dimensional phase modulated optical lattice. We observe the crossover from quantum to thermal fluctuations as the triggering mechanism for the nucleation of staggered states. In good quantitative agreement with numerical simulations based on the truncated Wigner method, we experimentally investigate how the nucleation time varies with the renormalized tunneling rate, the atomic density, and the driving frequency. The effective inverted energy band in the driven lattice is identified as the key ingredient which explains the emergence of gap solitons as observed in numerics and the possibility to nucleate staggered states from interband excitations as reported experimentally.

Electron scattering in large water clusters from photoelectron imaging with high harmonic radiation.

PubMed

Gartmann, Thomas E; Hartweg, Sebastian; Ban, Loren; Chasovskikh, Egor; Yoder, Bruce L; Signorell, Ruth

2018-06-06

Low-energy electron scattering in water clusters (H2O)n with average cluster sizes of n < 700 is investigated by angle-resolved photoelectron spectroscopy using high harmonic radiation at photon energies of 14.0, 20.3, and 26.5 eV for ionization from the three outermost valence orbitals. The measurements probe the evolution of the photoelectron anisotropy parameter β as a function of cluster size. A remarkably steep decrease of β with increasing cluster size is observed, which for the largest clusters reaches liquid bulk values. Detailed electron scattering calculations reveal that neither gas nor condensed phase scattering can explain the cluster data. Qualitative agreement between experiment and simulations is obtained with scattering calculations that treat cluster scattering as an intermediate case between gas and condensed phase scattering.

Shaken Snow Globes: Kinematic Tracers of the Multiphase Condensation Cascade in Massive Galaxies, Groups, and Clusters

NASA Astrophysics Data System (ADS)

Gaspari, M.; McDonald, M.; Hamer, S. L.; Brighenti, F.; Temi, P.; Gendron-Marsolais, M.; Hlavacek-Larrondo, J.; Edge, A. C.; Werner, N.; Tozzi, P.; Sun, M.; Stone, J. M.; Tremblay, G. R.; Hogan, M. T.; Eckert, D.; Ettori, S.; Yu, H.; Biffi, V.; Planelles, S.

2018-02-01

We propose a novel method to constrain turbulence and bulk motions in massive galaxies, galaxy groups, and clusters, exploring both simulations and observations. As emerged in the recent picture of top-down multiphase condensation, hot gaseous halos are tightly linked to all other phases in terms of cospatiality and thermodynamics. While hot halos (∼107 K) are perturbed by subsonic turbulence, warm (∼104 K) ionized and neutral filaments condense out of the turbulent eddies. The peaks condense into cold molecular clouds (<100 K) raining in the core via chaotic cold accretion (CCA). We show that all phases are tightly linked in terms of the ensemble (wide-aperture) velocity dispersion along the line of sight. The correlation arises in complementary long-term AGN feedback simulations and high-resolution CCA runs, and is corroborated by the combined Hitomi and new Integral Field Unit measurements in the Perseus cluster. The ensemble multiphase gas distributions (from the UV to the radio band) are characterized by substantial spectral line broadening (σ v,los ≈ 100–200 {km} {{{s}}}-1) with a mild line shift. On the other hand, pencil-beam detections (as H I absorption against the AGN backlight) sample the small-scale clouds displaying smaller broadening and significant line shifts of up to several 100 {km} {{{s}}}-1 (for those falling toward the AGN), with increased scatter due to the turbulence intermittency. We present new ensemble σ v,los of the warm Hα+[N II] gas in 72 observed cluster/group cores: the constraints are consistent with the simulations and can be used as robust proxies for the turbulent velocities, in particular for the challenging hot plasma (otherwise requiring extremely long X-ray exposures). Finally, we show that the physically motivated criterion C ≡ t cool/t eddy ≈ 1 best traces the condensation extent region and the presence of multiphase gas in observed clusters and groups. The ensemble method can be applied to many available spectroscopic data sets and can substantially advance our understanding of multiphase halos in light of the next-generation multiwavelength missions.

Universal Themes of Bose-Einstein Condensation

NASA Astrophysics Data System (ADS)

Proukakis, Nick P.; Snoke, David W.; Littlewood, Peter B.

2017-04-01

Foreword; List of contributors; Preface; Part I. Introduction: 1. Universality and Bose-Einstein condensation: perspectives on recent work D. W. Snoke, N. P. Proukakis, T. Giamarchi and P. B. Littlewood; 2. A history of Bose-Einstein condensation of atomic hydrogen T. Greytak and D. Kleppner; 3. Twenty years of atomic quantum gases: 1995-2015 W. Ketterle; 4. Introduction to polariton condensation P. B. Littlewood and A. Edelman; Part II. General Topics: Editorial notes; 5. The question of spontaneous symmetry breaking in condensates D. W. Snoke and A. J. Daley; 6. Effects of interactions on Bose-Einstein condensation R. P. Smith; 7. Formation of Bose-Einstein condensates M. J. Davis, T. M. Wright, T. Gasenzer, S. A. Gardiner and N. P. Proukakis; 8. Quenches, relaxation and pre-thermalization in an isolated quantum system T. Langen and J. Schmiedmayer; 9. Ultracold gases with intrinsic scale invariance C. Chin; 10. Berezinskii-Kosterlitz-Thouless phase of a driven-dissipative condensate N. Y. Kim, W. H. Nitsche and Y. Yamamoto; 11. Superfluidity and phase correlations of driven dissipative condensates J. Keeling, L. M. Sieberer, E. Altman, L. Chen, S. Diehl and J. Toner; 12. BEC to BCS crossover from superconductors to polaritons A. Edelman and P. B. Littlewood; Part III. Condensates in Atomic Physics: Editorial notes; 13. Probing and controlling strongly correlated quantum many-body systems using ultracold quantum gases I. Bloch; 14. Preparing and probing chern bands with cold atoms N. Goldman, N. R. Cooper and J. Dalibard; 15. Bose-Einstein condensates in artificial gauge fields L. J. LeBlanc and I. B. Spielman; 16. Second sound in ultracold atomic gases L. Pitaevskii and S. Stringari; 17. Quantum turbulence in atomic Bose-Einstein condensates N. G. Parker, A. J. Allen, C. F. Barenghi and N. P. Proukakis; 18. Spinor-dipolar aspects of Bose-Einstein condensation M. Ueda; Part IV. Condensates in Condensed Matter Physics: Editorial notes; 19. Bose-Einstein condensation of photons and grand-canonical condensate fluctuations J. Klaers and M. Weitz; 20. Laser operation and Bose-Einstein condensation: analogies and differences A. Chiocchetta, A. Gambassi and I. Carusotto; 21. Vortices in resonant polariton condensates in semiconductor microcavities D. N. Krizhanovskii, K. Guda, M. Sich, M. S. Skolnick, L. Dominici and D. Sanvitto; 22. Optical control of polariton condensates G. Christmann, P. G. Savvidis and J. J. Baumberg; 23. Disorder, synchronization and phase-locking in non-equilibrium Bose-Einstein condensates P. R. Eastham and B. Rosenow; 24. Collective topological excitations in 1D polariton quantum fluids H. Terças, D. D. Solnyshkov and G. Malpuech; 25. Microscopic theory of Bose-Einstein condensation of magnons at room temperature H. Salman, N. G. Berloff and S. O. Demokritov; 26. Spintronics and magnon Bose-Einstein condensation R. A. Duine, A. Brataas, S. A. Bender and Y. Tserkovnyak; 27. Spin-superfluidity and spin-current mediated non-local transport H. Chen and A. H. MacDonald; 28. Bose-Einstein condensation in quantum magnets C. Kollath, T. Giamarchi and C. Rüegg; Part V. Condensates in Astrophysics and Cosmology: Editorial notes; 29. Bose-Einstein condensates in neutron stars C. J. Pethick, T. Schäfer and A. Schwenk; 30. A simulated cosmological metric: the superfluid 3He condensate G. R. Pickett; 31. Cosmic axion Bose-Einstein condensation N. Banik and P. Sikivie; 32. Graviton BECs: a new approach to quantum gravity G. Dvali and C. Gomez; Universal Bose-Einstein condensation workshop; Index.

Formation of double front detonations of a condensed-phase explosive with powdered aluminium

NASA Astrophysics Data System (ADS)

Kim, Wuhyun; Gwak, Min-cheol; Yoh, Jack J.

2018-03-01

The performance characteristics of aluminised high explosive are considered by varying the aluminium (Al) mass fraction in a hybrid non-ideal detonation model. Since the time scales of the characteristic induction and combustion of high explosives and Al particles differ, the process of energy release behind the leading detonation wave front occurs over an extended period of time. Two cardinal observations are reported: a decrease in detonation velocity with an increase in Al mass fraction and a double front detonation (DFD) feature when anaerobic Al reaction occurs behind the front. In order to simulate the performance characteristics due to the varying Al mass fraction, the tetrahexamine tetranitramine (HMX) is considered as a base high explosive when formulating the multiphase conservation laws of mass, momentum, and energy exchanges between particles and HMX product gases. While experimental studies have been reported on the effect of Al mass fraction on both gas-phase and solid-phase detonations, the numerical investigations have been limited to only gas-phase detonation for the varying Al particles in the mixture. In the current study, a two-phase model is utilised for understanding the volumetric effects of Al mass fraction in condensed phase detonations. A series of unconfined and confined rate sticks are considered for characterising the performance of aluminised HMX with a maximum Al mass fraction of 50%. The simulated results are compared with the experimental data for 5-25% mass fractions, and the higher mass fraction behaviours are consistent with the experimental observations.

Phase 2 testing results of immobilization of WTP effluent management facility vaporator bottoms simulant

DOE Office of Scientific and Technical Information (OSTI.GOV)

Reigel, M.; Cozzi, A.; McCabe, D.

2017-09-08

The Hanford Waste Treatment and Immobilization Plant (WTP) Low Activity Waste (LAW) vitrification facility will generate an aqueous condensate recycle stream (LAW Melter Off-Gas Condensate) from the primary off-gas system. This work examined three waste form formulations based on previous testing with related simulants: 8 wt% ordinary portland cement (OPC), 47 wt% blast furnace slag (BFS), 45 wt% fly ash (FA) known as Cast Stone formulation; 20 wt% Aquaset® II-GH and 80 wt% BFS; 20 wt% OPC and 80 wt% BFS. These tests successfully produced one waste form that set within five days (Cast Stone formulation); however the other twomore » formulations, Aquaset® II-GH/BFS and OPC/BFS, took approximately eight and fourteen days to set, respectively.« less

Diffuse CO2 degassing at Vesuvio, Italy

NASA Astrophysics Data System (ADS)

Frondini, Francesco; Chiodini, Giovanni; Caliro, Stefano; Cardellini, Carlo; Granieri, Domenico; Ventura, Guido

2004-10-01

At Vesuvio, a significant fraction of the rising hydrothermal-volcanic fluids is subjected to a condensation and separation process producing a CO2-rich gas phase, mainly expulsed through soil diffuse degassing from well defined areas called diffuse degassing structures (DDS), and a liquid phase that flows towards the outer part of the volcanic cone. A large amount of thermal energy is associated with the steam condensation process and subsequent cooling of the liquid phase. The total amount of volcanic-hydrothermal CO2 discharged through diffuse degassing has been computed through a sequential Gaussian simulation (sGs) approach based on several hundred accumulation chamber measurements and, at the time of the survey, amounted to 151 t d-1. The steam associated with the CO2 output, computed assuming that the original H2O/CO2 ratio of hydrothermal fluids is preserved in fumarolic effluents, is 553 t d-1, and the energy produced by the steam condensation and cooling of the liquid phase is 1.47×1012 J d-1 (17 MW). The location of the CO2 and temperature anomalies show that most of the gas is discharged from the inner part of the crater and suggests that crater morphology and local stratigraphy exert strong control on CO2 degassing and subsurface steam condensation. The amounts of gas and energy released by Vesuvio are comparable to those released by other volcanic degassing areas of the world and their estimates, through periodic surveys of soil CO2 flux, can constitute a useful tool to monitor volcanic activity.

Condensation heat transfer and flow friction in silicon microchannels

NASA Astrophysics Data System (ADS)

Wu, Huiying; Wu, Xinyu; Qu, Jian; Yu, Mengmeng

2008-11-01

An experimental investigation was performed on heat transfer and flow friction characteristics during steam condensation flow in silicon microchannels. Three sets of trapezoidal silicon microchannels, with hydraulic diameters of 77.5 µm, 93.0 µm and 128.5 µm respectively, were tested under different flow and cooling conditions. It was found that both the condensation heat transfer Nusselt number (Nu) and the condensation two-phase frictional multiplier (phi2Lo) were dependent on the steam Reynolds number (Rev), condensation number (Co) and dimensionless hydraulic diameter (Dh/L). With the increase in the steam Reynolds number, condensation number and dimensionless hydraulic diameter, the condensation Nusselt number increased. However, different variations were observed for the condensation two-phase frictional multiplier. With the increase in the steam Reynolds number and dimensionless hydraulic diameter, the condensation two-phase frictional multiplier decreased, while with the increase in the condensation number, the condensation two-phase frictional multiplier increased. Based on the experimental results, dimensionless correlations for condensation heat transfer and flow friction in silicon microchannels were proposed for the first time. These correlations can be used to determine the condensation heat transfer coefficient and pressure drop in silicon microchannels if the steam mass flow rate, cooling rate and geometric parameters are fixed. It was also found that the condensation heat transfer and flow friction have relations to the injection flow (a transition flow pattern from the annular flow to the slug/bubbly flow), and with injection flow moving toward the outlet, both the condensation heat transfer coefficient and the condensation two-phase frictional multiplier increased.

Simulation of homogeneous condensation of small polyatomic systems in high pressure supersonic nozzle flows using Bhatnagar-Gross-Krook model

NASA Astrophysics Data System (ADS)

Kumar, Rakesh; Levin, Deborah A.

2011-03-01

In the present work, we have simulated the homogeneous condensation of carbon dioxide and ethanol using the Bhatnagar-Gross-Krook based approach. In an earlier work of Gallagher-Rogers et al. [J. Thermophys. Heat Transfer 22, 695 (2008)], it was found that it was not possible to simulate condensation experiments of Wegener et al. [Phys. Fluids 15, 1869 (1972)] using the direct simulation Monte Carlo method. Therefore, in this work, we have used the statistical Bhatnagar-Gross-Krook approach, which was found to be numerically more efficient than direct simulation Monte Carlo method in our previous studies [Kumar et al., AIAA J. 48, 1531 (2010)], to model homogeneous condensation of two small polyatomic systems, carbon dioxide and ethanol. A new weighting scheme is developed in the Bhatnagar-Gross-Krook framework to reduce the computational load associated with the study of homogeneous condensation flows. The solutions obtained by the use of the new scheme are compared with those obtained by the baseline Bhatnagar-Gross-Krook condensation model (without the species weighting scheme) for the condensing flow of carbon dioxide in the stagnation pressure range of 1-5 bars. Use of the new weighting scheme in the present work makes the simulation of homogeneous condensation of ethanol possible. We obtain good agreement between our simulated predictions for homogeneous condensation of ethanol and experiments in terms of the point of condensation onset and the distribution of mass fraction of ethanol condensed along the nozzle centerline.

Dark-soliton dynamics in Bose-Einstein condensates at finite temperature

DOE Office of Scientific and Technical Information (OSTI.GOV)

Jackson, B.; Proukakis, N. P.; Barenghi, C. F.

2007-05-15

The dynamics of a dark soliton in an elongated Bose-Einstein condensate is studied at finite temperatures. In addition to accurately reproducing all stages of the decay of the soliton observed in the experiment of Burger et al. [Phys. Rev. Lett. 83, 5198 (1999)], our numerical simulations reveal the existence of an experimentally accessible parameter regime for which phase-imprinted dark solitons can execute at least one full axial oscillation prior to their decay. The dependence of the decay time scale on temperature and initial soliton depth is analyzed and the role of interatomic collisions quantified.

A statistical mechanical model of economics

NASA Astrophysics Data System (ADS)

Lubbers, Nicholas Edward Williams

Statistical mechanics pursues low-dimensional descriptions of systems with a very large number of degrees of freedom. I explore this theme in two contexts. The main body of this dissertation explores and extends the Yard Sale Model (YSM) of economic transactions using a combination of simulations and theory. The YSM is a simple interacting model for wealth distributions which has the potential to explain the empirical observation of Pareto distributions of wealth. I develop the link between wealth condensation and the breakdown of ergodicity due to nonlinear diffusion effects which are analogous to the geometric random walk. Using this, I develop a deterministic effective theory of wealth transfer in the YSM that is useful for explaining many quantitative results. I introduce various forms of growth to the model, paying attention to the effect of growth on wealth condensation, inequality, and ergodicity. Arithmetic growth is found to partially break condensation, and geometric growth is found to completely break condensation. Further generalizations of geometric growth with growth in- equality show that the system is divided into two phases by a tipping point in the inequality parameter. The tipping point marks the line between systems which are ergodic and systems which exhibit wealth condensation. I explore generalizations of the YSM transaction scheme to arbitrary betting functions to develop notions of universality in YSM-like models. I find that wealth vi condensation is universal to a large class of models which can be divided into two phases. The first exhibits slow, power-law condensation dynamics, and the second exhibits fast, finite-time condensation dynamics. I find that the YSM, which exhibits exponential dynamics, is the critical, self-similar model which marks the dividing line between the two phases. The final chapter develops a low-dimensional approach to materials microstructure quantification. Modern materials design harnesses complex microstructure effects to develop high-performance materials, but general microstructure quantification is an unsolved problem. Motivated by statistical physics, I envision microstructure as a low-dimensional manifold, and construct this manifold by leveraging multiple machine learning approaches including transfer learning, dimensionality reduction, and computer vision breakthroughs with convolutional neural networks.

Disequilibrium condensation environments in space - A frontier in thermodynamics

NASA Technical Reports Server (NTRS)

De, B. R.

1979-01-01

The thermal-disequilibrium aspect of the problem of dust-particle formation from a gas phase in an open space environment is discussed in an effort to draw attention to the space condensation environment as an interesting arena for application and extension of the ideas and formalisms of nonequilibrium thermodynamics. It is shown that quasi-steady states with a disequilibrium between the gas-phase kinetic temperature and the condensed-phase internal temperature appear to be the norm of condensation environments in space. Consideration of the case of condensation onto a bulk condensed phase indicates that these quasi-steady states may constitute Prigogine dissipative structures. It is suggested that a proper study of the process of condensation in a space environment should include any effects arising from thermal disequilibrium.

Multiscale Reactive Molecular Dynamics

DTIC Science & Technology

2012-08-15

biology cannot be described without considering electronic and nuclear-level dynamics and their coupling to slower, cooperative motions of the system ...coupling to slower, cooperative motions of the system . These inherently multiscale problems require computationally efficient and accurate methods to...condensed phase systems with computational efficiency orders of magnitudes greater than currently possible with ab initio simulation methods, thus

Velocity and void distribution in a counter-current two-phase flow

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gabriel, S.; Schulenberg, T.; Laurien, E.

2012-07-01

Different flow regimes were investigated in a horizontal channel. Simulating a hot leg injection in case of a loss of coolant accident or flow conditions in reflux condenser mode, the hydraulic jump and partially reversed flow were identified as major constraints for a high amount of entrained water. Trying to simulate the reflux condenser mode, the test section now includes an inclined section connected to a horizontal channel. The channel is 90 mm high and 110 mm wide. Tests were carried out for water and air at ambient pressure and temperature. High speed video-metry was applied to obtain velocities frommore » flow pattern maps of the rising and falling fluid. In the horizontal part of the channel with partially reversed flow the fluid velocities were measured by planar particle image velocimetry. To obtain reliable results for the gaseous phase, this analysis was extended by endoscope measurements. Additionally, a new method based on the optical refraction at the interface between air and water in a back-light was used to obtain time-averaged void fraction. (authors)« less

CFD simulation of water vapour condensation in the presence of non-condensable gas in vertical cylindrical condensers.

PubMed

Li, Jun-De

2013-02-01

This paper presents the simulation of the condensation of water vapour in the presence of non-condensable gas using computational fluid dynamics (CFD) for turbulent flows in a vertical cylindrical condenser tube. The simulation accounts for the turbulent flow of the gas mixture, the condenser wall and the turbulent flow of the coolant in the annular channel with no assumptions of constant wall temperature or heat flux. The condensate film is assumed to occupy a negligible volume and its effect on the condensation of the water vapour has been taken into account by imposing a set of boundary conditions. A new strategy is used to overcome the limitation of the currently available commercial CFD package to solve the simultaneous simulation of flows involving multispecies and fluids of gas and liquid in separate channels. The results from the CFD simulations are compared with the experimental results from the literature for the condensation of water vapour with air as the non-condensable gas and for inlet mass fraction of the water vapour from 0.66 to 0.98. The CFD simulation results in general agree well with the directly measured quantities and it is found that the variation of heat flux in the condenser tube is more complex than a simple polynomial curve fit. The CFD results also show that, at least for flows involving high water vapour content, the axial velocity of the gas mixture at the interface between the gas mixture and the condensate film is in general not small and cannot be neglected.

CFD simulation of water vapour condensation in the presence of non-condensable gas in vertical cylindrical condensers

PubMed Central

Li, Jun-De

2013-01-01

This paper presents the simulation of the condensation of water vapour in the presence of non-condensable gas using computational fluid dynamics (CFD) for turbulent flows in a vertical cylindrical condenser tube. The simulation accounts for the turbulent flow of the gas mixture, the condenser wall and the turbulent flow of the coolant in the annular channel with no assumptions of constant wall temperature or heat flux. The condensate film is assumed to occupy a negligible volume and its effect on the condensation of the water vapour has been taken into account by imposing a set of boundary conditions. A new strategy is used to overcome the limitation of the currently available commercial CFD package to solve the simultaneous simulation of flows involving multispecies and fluids of gas and liquid in separate channels. The results from the CFD simulations are compared with the experimental results from the literature for the condensation of water vapour with air as the non-condensable gas and for inlet mass fraction of the water vapour from 0.66 to 0.98. The CFD simulation results in general agree well with the directly measured quantities and it is found that the variation of heat flux in the condenser tube is more complex than a simple polynomial curve fit. The CFD results also show that, at least for flows involving high water vapour content, the axial velocity of the gas mixture at the interface between the gas mixture and the condensate film is in general not small and cannot be neglected. PMID:24850953

Spin Order and Phase Transitions in Chains of Polariton Condensates.

PubMed

Ohadi, H; Ramsay, A J; Sigurdsson, H; Del Valle-Inclan Redondo, Y; Tsintzos, S I; Hatzopoulos, Z; Liew, T C H; Shelykh, I A; Rubo, Y G; Savvidis, P G; Baumberg, J J

2017-08-11

We demonstrate that multiply coupled spinor polariton condensates can be optically tuned through a sequence of spin-ordered phases by changing the coupling strength between nearest neighbors. For closed four-condensate chains these phases span from ferromagnetic (FM) to antiferromagnetic (AFM), separated by an unexpected crossover phase. This crossover phase is composed of alternating FM-AFM bonds. For larger eight-condensate chains, we show the critical role of spatial inhomogeneities and demonstrate a scheme to overcome them and prepare any desired spin state. Our observations thus demonstrate a fully controllable nonequilibrium spin lattice.

Phase collapse and revival of a 1-mode Bose-Einstein condensate induced by an off-resonant optical probe field and superselection rules

NASA Astrophysics Data System (ADS)

Arruda, L. G. E.; Prataviera, G. A.; de Oliveira, M. C.

2018-02-01

Phase collapse and revival for Bose-Einstein condensates are nonlinear phenomena appearing due to atomic collisions. While it has been observed in a general setting involving many modes, for one-mode condensates its occurrence is forbidden by the particle number superselection rule (SSR), which arises because there is no phase reference available. We consider a single mode atomic Bose-Einstein condensate interacting with an off-resonant optical probe field. We show that the condensate phase revival time is dependent on the atom-light interaction, allowing optical control on the atomic collapse and revival dynamics. Incoherent effects over the condensate phase are included by considering a continuous photo-detection over the probe field. We consider conditioned and unconditioned photo-counting events and verify that no extra control upon the condensate is achieved by the probe photo-detection, while further inference of the atomic system statistics is allowed leading to a useful test of the SSR on particle number and its imposition on the kind of physical condensate state.

Comparison of Chain Conformation of Poly(vinyl alcohol) in Solutions and Melts from Quantum Chemistry Based Molecular Dynamics Simulations

NASA Technical Reports Server (NTRS)

Jaffe, Richard; Han, Jie; Matsuda, Tsunetoshi; Yoon, Do; Langhoff, Stephen R. (Technical Monitor)

1997-01-01

Confirmations of 2,4-dihydroxypentane (DHP), a model molecule for poly(vinyl alcohol), have been studied by quantum chemistry (QC) calculations and molecular dynamics (MD) simulations. QC calculations at the 6-311G MP2 level show the meso tt conformer to be lowest in energy followed by the racemic tg, due to intramolecular hydrogen bond between the hydroxy groups. The Dreiding force field has been modified to reproduce the QC conformer energies for DHP. MD simulations using this force field have been carried out for DHP molecules in the gas phase, melt, and CHCl3 and water solutions. Extensive intramolecular hydrogen bonding is observed for the gas phase and CHCl3 solution, but not for the melt or aqueous solution, Such a condensed phase effect due to intermolecular interactions results in a drastic change in chain conformations, in agreement with experiments.

Dynamical diagnostics of the SST annual cycle in the eastern equatorial Pacific: Part II analysis of CMIP5 simulations

NASA Astrophysics Data System (ADS)

Chen, Ying-Ying; Jin, Fei-Fei

2017-12-01

In this study, a simple coupled framework established in Part I is utilized to investigate inter-model diversity in simulating the equatorial Pacific SST annual cycle (SSTAC). It demonstrates that the simulated amplitude and phase characteristics of SSTAC in models are controlled by two internal dynamical factors (the damping rate and phase speed) and two external forcing factors (the strength of the annual and semi-annual harmonic forcing). These four diagnostic factors are further condensed into a dynamical response factor and a forcing factor to derive theoretical solutions of amplitude and phase of SSTAC. The theoretical solutions are in remarkable agreement with observations and CMIP5 simulations. The great diversity in the simulated SSTACs is related to the spreads in these dynamic and forcing factors. Most models tend to simulate a weak SSTAC, due to their weak damping rate and annual harmonic forcing. The latter is due to bias in the meridional asymmetry of the annual mean state of the tropical Pacific, represented by the weak cross-equatorial winds in the cold tongue region.

A scale-bridging modeling approach for anisotropic organic molecules at patterned semiconductor surfaces

NASA Astrophysics Data System (ADS)

Kleppmann, Nicola; Klapp, Sabine H. L.

2015-02-01

Hybrid systems consisting of organic molecules at inorganic semiconductor surfaces are gaining increasing importance as thin film devices for optoelectronics. The efficiency of such devices strongly depends on the collective behavior of the adsorbed molecules. In the present paper, we propose a novel, coarse-grained model addressing the condensed phases of a representative hybrid system, that is, para-sexiphenyl (6P) at zinc-oxide (ZnO). Within our model, intermolecular interactions are represented via a Gay-Berne potential (describing steric and van-der-Waals interactions) combined with the electrostatic potential between two linear quadrupoles. Similarly, the molecule-substrate interactions include a coupling between a linear molecular quadrupole to the electric field generated by the line charges characterizing ZnO(10-10). To validate our approach, we perform equilibrium Monte Carlo simulations, where the lateral positions are fixed to a 2D lattice, while the rotational degrees of freedom are continuous. We use these simulations to investigate orientational ordering in the condensed state. We reproduce various experimentally observed features such as the alignment of individual molecules with the line charges on the surface, the formation of a standing uniaxial phase with a herringbone structure, as well as the formation of a lying nematic phase.

Gas network model allows full reservoir coupling

DOE Office of Scientific and Technical Information (OSTI.GOV)

Methnani, M.M.

The gas-network flow model (Gasnet) developed for and added to an existing Qatar General Petroleum Corp. (OGPC) in-house reservoir simulator, allows improved modeling of the interaction among the reservoir, wells, and pipeline networks. Gasnet is a three-phase model that is modified to handle gas-condensate systems. The numerical solution is based on a control volume scheme that uses the concept of cells and junctions, whereby pressure and phase densities are defined in cells, while phase flows are defined at junction links. The model features common numerical equations for the reservoir, the well, and the pipeline components and an efficient state-variable solutionmore » method in which all primary variables including phase flows are solved directly. Both steady-state and transient flow events can be simulated with the same tool. Three test cases show how the model runs. One case simulates flow redistribution in a simple two-branch gas network. The second simulates a horizontal gas well in a waterflooded gas reservoir. The third involves an export gas pipeline coupled to a producing reservoir.« less

Modeling the Phase Composition of Gas Condensate in Pipelines

NASA Astrophysics Data System (ADS)

Dudin, S. M.; Zemenkov, Yu D.; Shabarov, A. B.

2016-10-01

Gas condensate fields demonstrate a number of thermodynamic characteristics to be considered when they are developed, as well as when gas condensate is transported and processed. A complicated phase behavior of the gas condensate system, as well as the dependence of the extracted raw materials on the phase state of the deposit other conditions being equal, is a key aspect. Therefore, when designing gas condensate lines the crucial task is to select the most appropriate methods of calculating thermophysical properties and phase equilibrium of the transported gas condensate. The paper describes a physical-mathematical model of a gas-liquid flow in the gas condensate line. It was developed based on balance equations of conservation of mass, impulse and energy of the transported medium within the framework of a quasi-1D approach. Constitutive relationships are given separately, and practical recommendations on how to apply the research results are provided as well.

Pore diameter effects on phase behavior of a gas condensate in graphitic one-and two-dimensional nanopores.

PubMed

Welch, William R W; Piri, Mohammad

2016-01-01

Molecular dynamics (MD) simulations were performed on a hydrocarbon mixture representing a typical gas condensate composed mostly of methane and other small molecules with small fractions of heavier hydrocarbons, representative of mixtures found in tight shale reservoirs. The fluid was examined both in bulk and confined to graphitic nano-scale slits and pores. Numerous widths and diameters of slits and pores respectively were examined under variable pressures at 300 K in order to find conditions in which the fluid at the center of the apertures would not be affected by capillary condensation due to the oil-wet walls. For the bulk fluid, retrograde phase behavior was verified by liquid volumes obtained from Voronoi tessellations. In cases of both one and two-dimensional confinement, for the smallest apertures, heavy molecules aggregated inside the pore space and compression of the gas outside the solid structure lead to decreases in density of the confined fluid. Normal density/pressure relationships were observed for slits having gaps of above 3 nm and pores having diameters above 6 nm. At 70 bar, the minimum gap width at which the fluid could pass through the center of slits without condensation effects was predicted to be 6 nm and the corresponding diameter in pores was predicted to be 8 nm. The models suggest that in nanoscale networks involving pores smaller than these limiting dimensions, capillary condensation should significantly impede transmission of natural gases with similar composition.

Communication: On the consistency of approximate quantum dynamics simulation methods for vibrational spectra in the condensed phase.

PubMed

Rossi, Mariana; Liu, Hanchao; Paesani, Francesco; Bowman, Joel; Ceriotti, Michele

2014-11-14

Including quantum mechanical effects on the dynamics of nuclei in the condensed phase is challenging, because the complexity of exact methods grows exponentially with the number of quantum degrees of freedom. Efforts to circumvent these limitations can be traced down to two approaches: methods that treat a small subset of the degrees of freedom with rigorous quantum mechanics, considering the rest of the system as a static or classical environment, and methods that treat the whole system quantum mechanically, but using approximate dynamics. Here, we perform a systematic comparison between these two philosophies for the description of quantum effects in vibrational spectroscopy, taking the Embedded Local Monomer model and a mixed quantum-classical model as representatives of the first family of methods, and centroid molecular dynamics and thermostatted ring polymer molecular dynamics as examples of the latter. We use as benchmarks D2O doped with HOD and pure H2O at three distinct thermodynamic state points (ice Ih at 150 K, and the liquid at 300 K and 600 K), modeled with the simple q-TIP4P/F potential energy and dipole moment surfaces. With few exceptions the different techniques yield IR absorption frequencies that are consistent with one another within a few tens of cm(-1). Comparison with classical molecular dynamics demonstrates the importance of nuclear quantum effects up to the highest temperature, and a detailed discussion of the discrepancies between the various methods let us draw some (circumstantial) conclusions about the impact of the very different approximations that underlie them. Such cross validation between radically different approaches could indicate a way forward to further improve the state of the art in simulations of condensed-phase quantum dynamics.

Conversion of fogwater and aerosol organic nitrogen to ammonium, nitrate, and NOx during exposure to simulated sunlight and ozone.

PubMed

Zhang, Qi; Anastasio, Cort

2003-08-15

Although organic nitrogen (ON) compounds are apparently ubiquitous in the troposphere, very little is known about their fate and transformations. As one step in addressing this issue, we have studied the transformations of bulk (uncharacterized) organic nitrogen in fogwaters and aerosol aqueous extracts during exposure to simulated sunlight and O3. Our results show that over the course of several hours of exposure a significant portion of condensed-phase organic nitrogen is transformed into ammonium, nitrite, nitrate, and NOx. For nitrite, there was both photochemical formation and destruction, resulting in a slow net loss. Ammonium and nitrate were formed at initial rates on the order of a few micromolar per hour in the bulk fogwaters, corresponding to formation rates of approximately 10 and 40 ng m(-3) h(-1), respectively, in ambient fog. The average initial formation rate (expressed as ng (m of air)(-3) h(-1)) of NH4+ in the aqueous extracts of fine particles (PM2.5) was approximately one-half of the corresponding fogwater value. Initial formation rates of NOx (i.e., NO + NO2) were equivalent to approximately 2-11 pptv h(-1) in the three fogwaters tested. Although the formation rates of ammonium and nitrate were relatively small as compared to their initial concentrations in fogwaters (approximately 200-2000 microM) and aerosol particles (approximately 400-1500 ng m(-3)), this photochemical mineralization and "renoxification" from condensed-phase organic N is a previously uncharacterized source of inorganic N in the atmosphere. This conversion also represents a new component in the biogeochemical cycle of nitrogen that might have significant influences on atmospheric composition, condensed-phase properties, and the ecological impacts of N deposition.

Production and condensation of organic gases in the atmosphere of Titan

NASA Technical Reports Server (NTRS)

Sagan, C.; Thompson, W. R.

1982-01-01

The rates and altitudes for the dissociation of atmospheric constituents on Titan are calculated for solar ultraviolet radiation, the solar wind, Saturn magnetospheric particles, the Saturn co-rotating plasma, and cosmic rays. Laboratory experiments show that a variety of simple gas phase organic molecules and more complex organic solids called tholins are produced by such irradiations of simulated Titanian atmospheres. Except for ultraviolet wavelengths longward of the methane photodissociation continuum, most dissociation events occur between about 3100 and 3600 km altitude, corresponding well to the region of EUV opacity detected by Voyager. For a wide variety of simple to moderately complex organic gases in the Titanian atmosphere, condensation occurs below the top of the main cloud deck at about 2825 km. It is proposed that such condensates, beginning with CH4 at about 2615 km, comprise the principal mass of the Titan clouds. There is a distinct tendency for the atmosphere of Titan to act as a fractional distillation device, molecules of greater complexity condensing out at higher altitudes.

Mechanics of Sister Chromatids studied with a Polymer Model