Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Liquid dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Dynamics and diffusion mechanism of low-density liquid silicon

A first-order phase transition from a high-density liquid to a low-density liquid has been proposed to explain the various thermodynamic anomies of water. It also has been proposed that such liquid–liquid phase transition would exist in supercooled silicon. Computer simulation studies show that, across the transition, the diffusivity drops roughly 2 orders of magnitude, and the structures exhibit considerable tetrahedral ordering. The resulting phase is a highly viscous, low-density liquid silicon. Investigations on the atomic diffusion of such a novel form of liquid silicon are of high interest. Here we report such diffusion results from molecular dynamics simulations using the classical Stillinger–Weber (SW) potential of silicon. We show that the atomic diffusion of the low-density liquid is highly correlated with local tetrahedral geometries. We also show that atoms diffuse through hopping processes within short ranges, which gradually accumulate to an overall random motion for long ranges as in normal liquids. There is a close relationship between dynamical heterogeneity and hopping process. We point out that the above diffusion mechanism is closely related to the strong directional bonding nature of the distorted tetrahedral network. Here, our work offers new insights into the complex behavior of the highly viscous low density liquid silicon, suggesting similar diffusion behaviors in other tetrahedral coordinated liquids that exhibit liquid–liquid phase transition such as carbon and germanium.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Uniform hydrogen fuel layers for inertial fusion targets by microgravity

A critical concern in the fabrication of targets for inertial confinement fusion (ICF) is ensuring that the hydrogenic (D(sub 2) or DT) fuel layer maintains spherical symmetry. Solid layered targets have structural integrity, but lack the needed surface smoothness. Liquid targets are inherently smooth, but suffer from gravitationally induced sagging. One method to reduce the effective gravitational field environment is freefall insertion into the target chamber. Another method to counterbalance field gravitational force is to use an applied magnetic field combined with a gradient field to induce a magnetic dipole force on the liquid fuel layer. Based on time dependent calculations of the dynamics of the liquid fuel layer in microgravity environments, we show that it may be possible to produce a liquid layered ICF target that satisfies both smoothness and symmetry requirements.

Parks, P. B.↗
DOE OSTI2022

Regimes of evaporation and mixing behaviors of nanodroplets at transcritical conditions

The objective of this paper is to examine the fundamental mechanisms responsible for the transition between subcritical evaporation and supercritical dense-fluid-mixing in the absence of convection effects, specifically focusing on the liquid–vapor interfacial dynamics. To isolate the dynamics of this transition process, we characterize the different physical behaviors exhibited by an -dodecane nanoscale droplet placed in different nitrogen ambient conditions across the fuel’s critical point. We employ a continuum-based interface-resolving diffuse-interface method to explore the underlying phase-exchange mechanisms that bring about such distinct dynamics. Following the comparison against molecular dynamics simulations and experiments of evaporating droplets and experimental data for vapor–liquid equilibria, a parametric study at various ambient conditions and droplet sizing is performed to identify four regimes of evaporation/mixing behaviors: sub- and supercritical droplet evaporation, and sub- and supercritical dense-fluid-mixing. It is shown that the distinction in the phase-exchange mechanisms in these four regimes are brought about by the different thermodynamic phases the droplet center can exhibit during the evaporation/mixing process: subcritical liquid, supercritical liquid-like, subcritical gaseous, and supercritical gas-like, respectively. It is shown that the subcritical dense-fluid-mixing behavior is a direct result of nanoconfinement of the liquid–vapor interfacial structure and thus is not present for large droplet sizes. Finally, the present study also shows that the supercritical phase-exchange dynamics can follow two different pathways: supercritical droplet-like evaporation and supercritical dense-fluid-mixing. Furthermore, promoting the early transition to supercritical dense-fluid-mixing can significantly expedite the phase-exchange process through the disintegration of the liquid-like droplet core.

42 ENGINEERING↗
DOE OSTI2019

Pseudoplasticity and Dynamic Interfacial Tension Relaxation Effects on Nucleate Pool Boiling in Aqueous Polymeric Liquids

Nucleate pool boiling heat transfer and its ebullient dynamics in polymeric solutions at atmospheric pressure saturated conditions are experimentally investigated. Three grades of hydroxyethyl cellulose (HEC) are used, which have intrinsic viscosity in the range 5.29 ≤ [ η ] ≤ 10.31 [dl/g]. Their aqueous solutions in different concentrations, with zero-shear viscosity in the range 0.0021 ≤ η 0 ≤ 0.0118 [N∙s/m 2 ], exhibit shear-thinning rheology in varying degrees, as well as gas–liquid interfacial tension relaxation and wetting. Boiling heat transfer in solutions with constant molar concentrations of each additive, which are greater than their respective critical polymer concentration C * , is seen to have anomalous characteristics. There is degradation in the heat transfer at low heat fluxes, relative to that in the solvent, where the postnucleation bubble dynamics in the partial boiling regime is dominated by viscous resistance of the polymeric solutions. At higher heat fluxes, however, there is enhancement of boiling heat transfer due to a complex interplay of pseudoplasticity and dynamic surface tension effects. The higher frequency vapor bubbling train with high interfacial shear rates in this fully developed boiling regime tends to be influenced by increasing shear-thinning and time-dependent differential interfacial tension relaxation at the dynamic gas–liquid interfaces.

Engineering↗
DOE OSTI2020

Isotope effects in molecular structures and electronic properties of liquid water via deep potential molecular dynamics based on the SCAN functional

Here, Feynman path-integral deep potential molecular dynamics (PI-DPMD) calculations have been employed to study both light (H 2 O) and heavy water (D 2 O) within the isothermalisobaric ensemble. In particular, the deep neural network is trained based on ab initio data obtained from the strongly constrained and appropriately normed (SCAN) exchange-correlation functional. Because of the lighter mass of hydrogen than deuteron, the properties of light water are more influenced by nuclear quantum effect than those of heavy water. Clear isotope effects are observed and analyzed in terms of hydrogen-bond structure and electronic properties of water that are closely associated with experimental observables. The molecular structures of both liquid H 2 O and D 2 O agree well with the data extracted from scattering experiments. The delicate isotope effects on radial distribution functions and angular distribution functions are well reproduced as well. Our approach demonstrates that deep neural network combined with SCAN functional based ab initio molecular dynamics provides an accurate theoretical tool for modeling water and its isotope effects.

36 MATERIALS SCIENCE↗
DOE OSTI2014

The zero-multipole summation method for estimating electrostatic interactions in molecular dynamics: Analysis of the accuracy and application to liquid systems

In the preceding paper [I. Fukuda, J. Chem. Phys. 139, 174107 (2013)], the zero-multipole (ZM) summation method was proposed for efficiently evaluating the electrostatic Coulombic interactions of a classical point charge system. The summation takes a simple pairwise form, but prevents the electrically non-neutral multipole states that may artificially be generated by a simple cutoff truncation, which often causes large energetic noises and significant artifacts. The purpose of this paper is to judge the ability of the ZM method by investigating the accuracy, parameter dependencies, and stability in applications to liquid systems. To conduct this, first, the energy-functional error was divided into three terms and each term was analyzed by a theoretical error-bound estimation. This estimation gave us a clear basis of the discussions on the numerical investigations. It also gave a new viewpoint between the excess energy error and the damping effect by the damping parameter. Second, with the aid of these analyses, the ZM method was evaluated based on molecular dynamics (MD) simulations of two fundamental liquid systems, a molten sodium-chlorine ion system and a pure water molecule system. In the ion system, the energy accuracy, compared with the Ewald summation, was better for a larger value of multipole moment l currently induced until l ≲ 3 on average. This accuracy improvement with increasing l is due to the enhancement of the excess-energy accuracy. However, this improvement is wholly effective in the total accuracy if the theoretical moment l is smaller than or equal to a system intrinsic moment L. The simulation results thus indicate L ∼ 3 in this system, and we observed less accuracy in l = 4. We demonstrated the origins of parameter dependencies appearing in the crossing behavior and the oscillations of the energy error curves. With raising the moment l we observed, smaller values of the damping parameter provided more accurate results and smoother behaviors with respect to cutoff length were obtained. These features can be explained, on the basis of the theoretical error analyses, such that the excess energy accuracy is improved with increasing l and that the total accuracy improvement within l ⩽ L is facilitated by a small damping parameter. Although the accuracy was fundamentally similar to the ion system, the bulk water system exhibited distinguishable quantitative behaviors. A smaller damping parameter was effective in all the practical cutoff distance, and this fact can be interpreted by the reduction of the excess subset. A lower moment was advantageous in the energy accuracy, where l = 1 was slightly superior to l = 2 in this system. However, the method with l = 2 (viz., the zero-quadrupole sum) gave accurate results for the radial distribution function. We confirmed the stability in the numerical integration for MD simulations employing the ZM scheme. This result is supported by the sufficient smoothness of the energy function. Along with the smoothness, the pairwise feature and the allowance of the atom-based cutoff mode on the energy formula lead to the exact zero total-force, ensuring the total-momentum conservations for typical MD equations of motion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Role of Fast Dynamics in Conductivity of Polymerized Ionic Liquids

Polymerized ionic liquids (PolyILs) are promising candidates for a broad range of technologies. However, the relatively low conductivity of PolyILs at room temperature has strongly limited their applications. In this work, we provide new insights into the roles of various microscopic parameters controlling ion transport in these polymers, which are crucial for their rational design and practical applications. Using broadband dielectric spectroscopy and neutron and light scattering techniques, we found a clear connection between the activation energy for conductivity, fast dynamics, and high-frequency shear modulus in PolyILs at their glass transition temperature (T g ). In particular, our analysis reveals a correlation between conductivity and the amplitude of fast picosecond fluctuations at T g , suggesting the possible involvement of fast dynamics in lowering the energy barrier for ion conductivity. We also demonstrate that both the activation energy for ion transport and the amplitude of the fast fluctuations depend on the high-frequency shear moduli of PolyILs, thus identifying a practically important parameter for tuning conductivity. Finally, the parameters recognized in this work and their connection to the ionic conductivity of PolyILs set the stage for a deeper understanding of the mechanism of ion transport in PolyILs in the glassy state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1989

Interaction of a sodium ion with the water liquid-vapor interface

Molecular dynamics results are presented for the density profile of a sodium ion near the water liquid-vapor interface at 320 K. These results are compared with the predictions of a simple dielectric model for the interaction of a monovalent ion with this interface. The interfacial region described by the model profile is too narrow and the profile decreases too abruptly near the solution interface. Thus, the simple model does not provide a satisfactory description of the molecular dynamics results for ion positions within two molecular diameters from the solution interface where appreciable ion concentrations are observed. These results suggest that surfaces associated with dielectric models of ionic processes at aqueous solution interfaces should be located at least two molecular diameters inside the liquid phase. A free energy expense of about 2 kcal/mol is required to move the ion within two molecular layers of the free water liquid-vapor interface.

NASA Center ARC↗
DOE OSTI2022

Investigation of γ - ( U , Z r ) structural properties and its interfacial properties with liquid sodium using ab initio molecular dynamics

In this study, the elastic properties, structural parameters, sound velocity, and Debye temperature of γ–(U,Zr) were computed using ab initio molecular dynamics (AIMD) at temperatures between 1000 K and 1400 K and for Zr content between 0 at. % and 100 at. %. UZr is used as a metallic fuel for Sodium Fast Reactors (SFRs). The study of the mechanical and thermal behavior of these alloys leads to a better data-informed fuel design. The bulk modulus, shear modulus, Young's modulus, and Poisson's ratio were calculated from the elastic constants and their dependence on Zr content and temperature was investigated, comparing the results with previous computational work and the available experimental data in the literature. Interfacial properties between UZr (up to 32 at. % which typically exists in nuclear fuel) and liquid sodium are also of interest due to the presence of a sodium bond between the fuel and the cladding in metallic nuclear fuel. The interfacial energy between γ–(U,Zr) and liquid sodium, the surface tension of liquid sodium, and the work of adhesion were computed at different temperatures and Zr concentrations. It was demonstrated that γ–(U,Zr) is completely wetted by liquid sodium at all the investigated temperatures and Zr concentrations. Finally, this work provides the basis for the determination of interfacial resistances in SFRs and their implementation into heat transfer fuel performance simulations, which will be the subject of future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

Understanding the nonlinear dynamics of driven particles in supercooled liquids in terms of an effective temperature

In active microrheology, the mechanical properties of a material are tested by adding probe particles which are pulled by an external force. In case of supercooled liquids, strong forcing leads to a thinning of the host material which becomes more pronounced as the system approaches the glass transition. In this work, we provide a quantitative theoretical description of this thinning behavior based on the properties of the Potential Energy Landscape (PEL) of a model glass-former. A key role plays the trap-like nature of the PEL. We find that the mechanical properties in the strongly driven system behave the same as in a quiescent system at an enhanced temperature, giving rise to a well-characterized effective temperature. Furthermore, this effective temperature turns out to be independent of the chosen observable and individually shows up in the thermodynamic and dynamic properties of the system. Based on this underlying theoretical understanding, we can estimate its dependence on temperature and force by the PEL-properties of the quiescent system. We furthermore critically discuss the relevance of effective temperatures obtained by scaling relations for the description of out-of-equilibrium situations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

System Dynamics of the Competition of Municipal Solid Waste to Landfill, Electricity, and Liquid Fuel in California

A quantitative system dynamics model was created to evaluate the economic and environmental tradeoffs between biomass to electricity and to liquid fuel using MSW biomass in the state of California as a case study. From an environmental perspective, landfilling represents the worst use of MSW over time, generating more greenhouse gas (GHG) emissions compared to converting MSW to liquid fuel or to electricity. MSW to ethanol results in the greatest displacement of GHG emissions per dollar spent compared to MSW to electricity. MSW to ethanol could save the state of California approximately $60 billion in energy costs by 2050 compared to landfilling, while also reducing GHG emissions state-wide by approximately 140 million metric tons during that timeframe. MSW conversion to electricity creates a significant cost within the state's electricity sector, although some conversion technologies are cost competitive with existing renewable generation.

09 BIOMASS FUELS↗
DOE OSTI2020

Fast predictions of liquid-phase acid-catalyzed reaction rates using molecular dynamics simulations and convolutional neural networks

The rates of liquid-phase, acid-catalyzed reactions relevant to the upgrading of biomass into high-value chemicals are highly sensitive to solvent composition and identifying suitable solvent mixtures is theoretically and experimentally challenging. We show that the complex atomistic configurations of reactant–solvent environments generated by classical molecular dynamics simulations can be exploited by 3D convolutional neural networks to enable accurate predictions of Brønsted acid-catalyzed reaction rates for model biomass compounds. We develop a 3D convolutional neural network, which we call SolventNet, and train it to predict acid-catalyzed reaction rates using experimental reaction data and corresponding molecular dynamics simulation data for seven biomass-derived oxygenates in water–cosolvent mixtures. We show that SolventNet can predict reaction rates for additional reactants and solvent systems an order of magnitude faster than prior simulation methods. This combination of machine learning with molecular dynamics enables the rapid, high-throughput screening of solvent systems and identification of improved biomass conversion conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Structural Properties of [N1888][TFSI] Ionic Liquid: A Small Angle Neutron Scattering and Polarizable Molecular Dynamics Study

In this study, we investigate the quaternary ammonium-based ionic liquid (QAIL), methyltrioctylammonium bis(trifluoromethylsulfonyl)imide, [N 1888 ][TFSI], utilizing small angle neutron scattering (SANS) measurements and polarizable molecular dynamics (MD) simulations to characterize the shortand long-range liquid structure. Scattering structure factors show signatures of three length scales in reciprocal space indicative of alternating polarity (k ~ 0.44 Å –1 ), charge (k ~ 0.75 Å –1 ), and neighboring or adjacent (k ~ 1.46 Å –1 ) domains. Excellent agreement between simulation and experimental scattering structure factors validates various simulation analyses that provide detailed atomistic characterization of the different length scale correlations. The first solvation shell structure is illustrated by obtaining radial, angular, dihedral, and combined distribution functions, where two dominant spatial motifs, N + ···N – and N + ···O – , compete for optimal packing around the polar head of the [N 1888 ] + cation. Intermediate and long-range structures are governed by the balance between local electroneutrality and octyl chain networking, respectively. By computing the charge-correlation structure factor, S ZZ , and the spatial extent of the octyl chain network using graph theory, the bulk-phase structure of [N 1888 ][TFSI] is characterized in terms of electrostatic screening and apolar domain formation length scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Maximizing ion dynamics and electrochemical performance of ionic liquid-acetonitrile electrolyte in Ti 3 C 2 T x MXene

Modification of the structure and morphology of MXene electrodes and the formulation of the electrolytes used in their supercapacitor configurations are significant factors affecting the performance of electrochemical devices. In this study, we investigated the electrochemical performance and ion dynamics of 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide, [EmimTFSI], ionic liquid in the presence of acetonitrile (ACN) at different concentrations in Ti 3 C 2 T x MXene supercapacitor. We found an optimum concentration of ACN, at which more cations from the ionic liquid attach to the MXene electrode surface, providing higher electrochemical performance. This higher capacitance is also associated with increased microscopic dynamics of the cation away from the pore wall. These findings give a guideline to optimize the performance of MXene-based supercapacitors using organic solvents-ionic liquid-based electrolyte systems.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Aspects of model-based rocket engine condition monitoring and control

A rigorous propulsion system modelling method suitable for control and condition monitoring purposes is developed. Previously developed control oriented methods yielding nominal models for gaseous medium propulsion systems are extended to include both nominal and anomalous models for liquid mediums in the following two ways. First, thermodynamic and fluid dynamic properties for liquids such as liquid hydrogen are incorporated into the governing equations. Second, anomalous conditions are captured in ways compatible with existing system theoretic design tools so that anomalous models can be constructed. Control and condition monitoring based methods are seen as an improvement over some existing modelling methods because such methods typically do not rigorously lead to low order models nor do they provide a means for capturing anomalous conditions. Applications to the nominal SSME HPFP and degraded HPFP serve to illustrate the approach.

Karr, Gerald R.↗
NASA NTRS2016

Combustion Dynamics and Stability Modeling of a Liquid Oxygen/RP-2 Oxygen-Rich Staged Combustion Preburner and Thrust Chamber Assembly with Gas-Centered Swirl Coaxial Injector Elements

The Combustion Stability Tool Development (CSTD) project, funded by the Air Force Space and Missile Systems Center, began in March 2015 supporting a renewed interest in the development of a liquid oxygen/hydrocarbon, oxygen-rich combustion engine. The project encompasses the design, assembly, and hot-fire testing of the NASA Marshall Space Flight Center 40-klbf Integrated Test Rig (MITR). The test rig models a staged-combustion configuration by combining an oxygen-rich preburner (ORPB), to generate hot gas, with a thrust chamber assembly (TCA) using gas-centered swirl coaxial injector elements. There are five separately designed interchangeable injectors in the TCA that each contain 19- or 27- injector elements. A companion paper in this JANNAF conference describes the design characteristics, rationale, and fabrication issues for all the injectors. The data acquired from a heavily instrumented rig encompasses several injectors, several operating points, and stability bomb tests. Another companion paper in this JANNAF conference describes this test program in detail. In this paper, dynamic data from the hot-fire testing is characterized and used to identify the responses in the ORPB and TCA. A brief review of damping metrics are discussed and applied as a measure of stability margin for damped acoustic modes. Chug and longitudinal combustion stability models and predictions are described which includes new dynamic models for compressible flow through an orifice and a modification to incorporate a third feed line for inclusion of the fuel-film coolant. Flow-acoustics finite element modeling is used to investigate the anticipated TCA acoustics, the effects of injector element length on stability margin, and the potential use of an ORPB orifice trip ring for improving longitudinal stability margin.

Casiano, M. J.↗
DOE OSTI2015

Rock-salt structure lithium deuteride formation in liquid lithium with high-concentrations of deuterium: a first-principles molecular dynamics study

Because of lithium's possible use as a first wall material in a fusion reactor, a fundamental understanding of the interactions between liquid lithium (Li) and deuterium (D) is important. Here, we predict structural and dynamical properties of liquid Li samples with high concentrations of D, as derived from first-principles molecular dynamics simulations. Liquid Li samples with four concentrations of inserted D atoms (LiD$_{\beta}$ , $\beta =0.25$ , 0.50, 0.75, and 1.00) are studied at temperatures ranging from 470 to 1143 K. Densities, diffusivities, pair distribution functions, bond angle distribution functions, geometries, and charge transfer between Li and D atoms are calculated and analyzed. The analysis suggests liquid–solid phase transitions can occur at some concentrations and temperatures, forming rock-salt LiD within liquid Li. Finally, we observed the formation of some D 2 molecules at high D concentrations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Anisotropic charge dynamics in the quantum spin-liquid candidate κ - ( BEDT-TTF ) 2 Cu 2 ( CN ) 3

We have in detail characterized the anisotropic charge response of the dimer Mott insulator kappa-(BEDT-TTF)(2)Cu-2(CN)(3) by dc conductivity, Hall effect, and dielectric spectroscopy. At room temperature, the Hall coefficient is positive and close to the value expected from stoichiometry; the temperature behavior follows the dc resistivity rho(T). Within the planes the dc conductivity is well described by variable-range hopping in two dimensions; this model, however, fails for the out-of-plane direction. An unusually broad in-plane dielectric relaxation is detected below about 60 K; it slows down much faster than the dc conductivity following an Arrhenius law. At around 17 K, we can identify a pronounced dielectric anomaly concomitantly with anomalous features in the mean relaxation time and spectral broadening. The out-of-plane relaxation, on the other hand, shows a much weaker dielectric anomaly; it closely follows the temperature behavior of the respective dc resistivity. At lower temperatures, the dielectric constant becomes smaller both within and perpendicular to the planes; also, the relaxation levels off. The observed behavior bears features of relaxorlike ferroelectricity. Because heterogeneities impede its long-range development, only a weak-tunneling-like dynamics persists at low temperatures. We suggest that the random potential and domain structure gradually emerge due to the coupling to the anion network.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗