Search NASA⌕ Search

SEARCH · Search NASA

Results for “Liquid structure simulation”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Sockeye Theory Manual

Sockeye is an application that models heat pipe performance, based on the MOOSE framework. Its primary focus is on liquid-metal heat pipes with annular screen or porous wick structures, with the intended application being the simulation of heat pipes in microreactors. The purpose of this capability is to evaluate and understand heat pipe performance under different conditions. Furthermore, this heat pipe performance is intended to be coupled to multiphysics applications for modeling microreactors, so that the effects of heat pipe performance on the larger reactor system can be simulated.

42 ENGINEERING↗
DOE OSTI2025

FLAMES─Fast, Low-Storage, Accurate, and Memory-Efficient Adaptive Sampling─Approach to Resolve Spatially Dependent Dynamics of Molecular Liquids

Many critical phenomena in soft matter occur at large length scales, necessitating the resolution of their structure and dynamics at low wavenumbers. However, resolving wavenumber-dependent dynamics computationally via molecular dynamics simulations presents significant challenges, as these phenomena span several orders of magnitude in both time and length scales, resulting in high computational costs and memory demands. Here, this work highlights the computational and memory challenges associated with analyzing molecular trajectories in reciprocal space and demonstrates a method to address them. We introduce FLAMESFast, Low-storage, Accurate, and Memory-Efficient adaptive Sampling, which is a direct method for calculation of structure factors, allowing us to select only the required number of wavevectors for binning. We also use wavenumber-dependent time steps to extract dynamics. Our FLAMES approach effectively mitigates computational and memory/storage bottlenecks. We demonstrate the method using simulations of a model system, liquid octane, at various temperatures. Comparisons with experimental data and real space computation show that the FLAMES technique achieves high accuracy in resolving temperature- and spatially dependent dynamics while being significantly more computationally efficient and requiring less memory and storage than methods based on a uniform wavevector grid and fixed temporal spacing.

Chen, Guang [Argonne National Laboratory (ANL), Ar↗
DOE OSTI2024

Transferable Water Potentials Using Equivariant Neural Networks

Machine learning interatomic potentials (MLIPs) have emerged as a technique that promises quantum theory accuracy for reduced cost. It has been proposed [J. Chem. Phys. 2023, 158, 084111] that MLIPs trained on solely liquid water data cannot accurately transfer to the vapor–liquid equilibrium while recovering the many-body decomposition (MBD) analysis of gas-phase water clusters. This suggests that MLIPs do not directly learn the physically correct interactions of water molecules, limiting transferability. In this work, we show that MLIPs using equivariant architecture and trained on 3200 liquid water structures reproduces liquid-phase water properties (e.g., density within 0.003 g/cm 3 between 230 and 365 K), vapor–liquid equilibrium properties up to 550 K, the MBD analysis of gas-phase water cluster up to six-body interactions, and the relative energy and the vibrational density of states of ice phases. We show that potentials developed using equivariant MLIPs allow transferability for arbitrary phases of water that remain stable in nanosecond long simulations.

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

Three-dimensional atomic packing in amorphous solids with liquid-like structure

Liquids and solids are two fundamental states of matter. However, our understanding of their three-dimensional atomic structure is mostly based on physical models. Here we use atomic electron tomography to experimentally determine the three-dimensional atomic positions of monatomic amorphous solids, namely a Ta thin film and two Pd nanoparticles. We observe that pentagonal bipyramids are the most abundant atomic motifs in these amorphous materials. Instead of forming icosahedra, the majority of pentagonal bipyramids arrange into pentagonal bipyramid networks with medium-range order. Molecular dynamics simulations further reveal that pentagonal bipyramid networks are prevalent in monatomic metallic liquids, which rapidly grow in size and form more icosahedra during the quench from the liquid to the glass state. Furthermore, these results expand our understanding of the atomic structures of amorphous solids and will encourage future studies on amorphous–crystalline phase and glass transitions in non-crystalline materials with three-dimensional atomic resolution.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Path integral Monte Carlo approach to the structural properties and collective excitations of liquid $$^3{\text {He}}$$ without fixed nodes

Abstract Due to its nature as a strongly correlated quantum liquid, ultracold helium is characterized by the nontrivial interplay of different physical effects. Bosonic $$^4{\text {He}}$$ 4 He exhibits superfluidity and Bose-Einstein condensation. Its physical properties have been accurately determined on the basis of ab initio path integral Monte Carlo (PIMC) simulations. In contrast, the corresponding theoretical description of fermionic $$^3{\text {He}}$$ 3 He is severely hampered by the notorious fermion sign problem, and previous PIMC results have been derived by introducing the uncontrolled fixed-node approximation. In this work, we present extensive new PIMC simulations of normal liquid $$^3{\text {He}}$$ 3 He without any nodal constraints. This allows us to to unambiguously quantify the impact of Fermi statistics and to study the effects of temperature on different physical properties like the static structure factor $$S({\mathbf {q}})$$ S ( q ) , the momentum distribution $$n({\mathbf {q}})$$ n ( q ) , and the static density response function $$\chi ({\mathbf {q}})$$ χ ( q ) . In addition, the dynamic structure factor $$S({\mathbf {q}},\omega )$$ S ( q , ω ) is rigorously reconstructed from imaginary-time PIMC data. From simulations of $$^3{\text {He}}$$ 3 He , we derived the familiar phonon–maxon–roton dispersion function that is well-known for $$^4{\text {He}}$$ 4 He and has been reported previously for two-dimensional $$^3{\text {He}}$$ 3 He films (Nature 483:576–579 (2012)). The comparison of our new results for both $$S({\mathbf {q}})$$ S ( q ) and $$S({\mathbf {q}},\omega )$$ S ( q , ω ) with neutron scattering measurements reveals an excellent agreement between theory and experiment.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Machine learning accelerated turbulence modeling of transient flashing jets

We report modeling the sudden depressurization of superheated liquids through nozzles is a challenge because the pressure drop causes rapid flash boiling of the liquid. The resulting jet usually demonstrates a wide range of structures, including ligaments and droplets, due to both mechanical and thermodynamic effects. As the simulation comprises increasingly numerous phenomena, the computational cost begins to increase. One way to moderate the additional cost is to use machine learning surrogacy for specific elements of the calculation. This study presents a machine learning-assisted computational fluid dynamics approach for simulating the atomization of flashing liquids accounting for distinct stages, from primary atomization to secondary breakup to small droplets using the Σ - Y model coupled with the homogeneous relaxation model. Notably, the models for thermodynamic non-equilibrium (HRM) and Σ - Y are coupled, for the first time, with a deep neural network that simulates the turbulence quantities, which are then used in the prediction of superheated liquid jet atomization. The data-driven component of this method is used for turbulence modeling, avoiding the solution of the two-equation turbulence model typically used for Reynolds-averaged Navier-Stokes simulations for these problems. Both the accuracy and speed of the hybrid approach are evaluated, demonstrating adequate accuracy and at least 25% faster computational fluid dynamics simulations than the traditional approach. This acceleration suggests that perhaps additional components of the calculation could be replaced for even further benefit. Published under an exclusive license by AIP Publishing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Using Radiative Convective Equilibrium to Explore Clouds and Climate in the Community Atmosphere Model

Abstract Characteristics of, and fundamental differences between, the radiative‐convective equilibrium (RCE) climate states following the Radiative‐Convective Equilibrium Model Intercomparison Project (RCEMIP) protocols in the Community Atmosphere Model version 5 (CAM5) and version 6 (CAM6) are presented. This paper explores the characteristics of clouds, moisture, precipitation and circulation in the RCE state, as well as the tropical response to surface warming, in CAM5 and CAM6 with different parameterizations. Overall, CAM5 simulates higher precipitation rates that result in larger global average precipitation, despite lower outgoing longwave radiation compared to CAM6. Differences in the structure of clouds, particularly the amount and vertical location of cloud liquid, exist between the CAM versions and can, in part, be related to distinct representations of shallow convection and boundary layer processes. Both CAM5 and CAM6 simulate similar peaks in cloud fraction, relative humidity, and cloud ice, linked to the usage of a similar deep convection parameterization. These anvil clouds rise and decrease in extent in response to surface warming. More generally, extreme precipitation, aggregation of convection, and climate sensitivity increase with warming in both CAM5 and CAM6. This analysis provides a benchmark for future studies that explore clouds, convection, and climate in CAM with the RCEMIP protocols now available in the Community Earth System Model. These results are discussed within the context of realistic climate simulations using CAM5 and CAM6, highlighting the usefulness of a hierarchical modeling approach to understanding model and parameterization sensitivities to inform model development efforts.

Reed, Kevin A.↗
DOE OSTI2020

Low Mach number fluctuating hydrodynamics model for ionic liquids

We present a new mesoscale model for ionic liquids based on a low Mach number fluctuating hydrodynamics formulation for multicomponent charged species. The low Mach number approach eliminates sound waves from the fully compressible equations leading to a computationally efficient incompressible formulation. The model uses a Gibbs free energy functional that includes enthalpy of mixing, interfacial energy, and electrostatic contributions. These lead to a new fourth-order term in the mass equations and a reversible stress in the momentum equations. We calibrate our model using parameters for [DMPI+][F6P-], an extensively-studied room temperature ionic liquid (RTIL), and numerically demonstrate the formation of mesoscopic structuring at equilibrium in two and three dimensions. In simulations with electrode boundaries the measured double layer capacitance decreases with voltage, in agreement with theoretical predictions and experimental measurements for RTILs. Finally, we present a shear electroosmosis example to demonstrate that the methodology can be used to model electrokinetic flows.

42 ENGINEERING↗
DOE OSTI2019

Scale-Adaptive Simulation of Transient Two-Phase Flow in Continuous-Casting Mold

Scale-adaptive simulation (SAS) of the transient gas–liquid two-phase flow in a laboratory-scale continuous-casting mold is presented. The main objective is to investigate the applicability of the scale-adaptive unsteady Reynolds-averaged Navier–Stokes turbulent model (URANS SAS) for predicting the transient multiscale turbulent structures in a two-phase flow. Good quantitative agreements with the experimental data and the large eddy simulation (LES) results are obtained both for the time-averaged velocity field and for the transient turbulent characteristics. The introduction of the von Karman length-scale into the turbulence-scale equation allows the SAS model to dynamically adjust to the resolved turbulent structures. The LES-like pulsating behavior of the air gas and the large-scale liquid eddy magnitudes in the unsteady regions of flow field are captured by the SAS model. The classical − 5/3 law of power spectrum density (PSD) of the axial velocity is kept properly for the single-phase turbulent flow. For two-phase flow, the decay of PSD is too steep at the high-frequency region; the predicted PSD obtained with SAS is damped stronger than that estimated by LES. The SAS model offers an attractive alternative to the existing LES approach or to the other hybrid RANS/LES models for strongly unsteady flows.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Probing the solidification of quasicrystals via joint experiment and simulation (Final Report)

Quasicrystals (QCs) possess long-range positional order but non-crystallographic orientational order. Their classically ‘forbidden’ symmetry has long attracted the interest of scientists worldwide. Despite their frequent observation in both metal alloys and soft matter structures in the 35 years since their discovery, little is known about how QCs evolve from a liquid, amorphous, or crystalline precursor. In this project, we sought to resolve the enigma of QC self-assembly through a combined experimental and computational program. The Shahani group employed in situ electron and synchrotron X ray imaging to peer into the growth dynamics of QCs in a liquid, covering a broad range of length scales and solidification pathways. Glotzer's team developed new simulation models incorporating phasonic defects — which are unique to quasicrystals — and used this model to investigate two grains growing together, closely mimicking the experimental conditions. While QCs remain exceptional structures, most compounds in the realm of intermetallics adopt non-trivial geometries. In fact, only around 6% of phases are comprised of the simplest sphere packings that researchers so readily associate with metals. State-of-the-art knowledge of phase transformations at the outset of this collaborative project did not encompass the remaining 94% of intermetallic compounds that possess complex and aperiodic structure types. Therefore, we expect that outcomes from this project will have immediate and profound impact on synthesis and processing science: our efforts will be used to not only explain the growth dynamics of QCs but also complex intermetallics more broadly.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Spatially resolved structural order in low-temperature liquid electrolyte

Determining the degree and the spatial extent of structural order in liquids is a grand challenge. Here, we are able to resolve the structural order in a model organic electrolyte of 1 M lithium hexafluorophosphate (LiP 6 ) dissolved in 1:1 (v/v) ethylene carbonate:diethylcarbonate by developing an integrated method of liquid-phase transmission electron microscopy (TEM), cryo-TEM operated at -30°C, four-dimensional scanning TEM, and data analysis based on deep learning. This study reveals the presence of short-range order (SRO) in the high–salt concentration domains of the liquid electrolyte from liquid phase separation at the low temperature. Molecular dynamics simulations suggest the SRO originates from the Li + - (PF 6 - ) n (n > 2) local structural order induced by high LiP 6 salt concentration.

42 ENGINEERING↗
NASA NTRS2014

Investigation of Structure and Transport in Li-Doped Ionic Liquid Electrolytes: [pyr14][TFSI], [pyr13][FSI], [EMIM][BF4]

Ionic liquid electrolytes have been proposed as a means of improving the safety and cycling behavior of advanced lithium batteries; however, the properties of these electrolytes under high lithium doping are poorly understood. Here, we employ both polarizable molecular dynamics simulation and experiment to investigate the structure, thermodynamics and transport of three potential electrolytes, N-methyl-N-butylpyrrolidinium bis(trifluoromethylsufonyl)imide ([pyr14][TFSI]), N- methyl-N-propylpyrrolidinium bis(fluorosufonyl)imide ([pyr13][FSI]), and 1-ethyl-3-- methylimidazolium boron tetrafluoride ([EMIM][BF4]), as a function of Li (-) salt concentration and temperature. Structurally, Li(+) is shown to be solvated by three anion neighbors in [pyr14][TFSI] and four anion neighbors in both [pyr13][FSI] and [EMIM][BF4], and at all levels of xLi we find the presence of lithium aggregates. Furthermore, the computed density, diffusion, viscosity, and ionic conductivity show excellent agreement with experimental data. While the diffusion and viscosity exhibit a systematic decrease and increase, respectively, with increasing xLi, the contribution of Li(+) to ionic conductivity increases until reaching a saturation doping level of xLi 0.10. Comparatively, the Li(+) conductivity of [pyr14][TFSI] is an order of magnitude lower than that of the other liquids, which range between 0.1-0.3 mScm. The differences in Li(+) transport are reflected in the residence times of Li(+) with the anions, which are revealed to be much larger for [pyr14][TFSI] (up to 100 ns at the highest doping levels) than in either [EMIM][BF4] or [pyr13][FSI]. Finally, we comment on the relative kinetics of Li(+) transport in each liquid and we present strong evidence for transport through anion exchange (hopping) as opposed to the net motion of Li(+) with its solvation shell (vehicular).

ionic liquids↗
NASA NTRS2014

Investigation of Structure and Transport in Li-Doped Ionic Liquid Electrolytes: [pyr14][TFSI], [pyr13][FSI], [EMIM][BF4]

Ionic liquid electrolytes have been proposed as a means of improving the safety and cycling behavior of advanced lithium batteries; however, the properties of these electrolytes under high lithium doping are poorly understood. Here, we employ both polarizable molecular dynamics simulation and experiment to investigate the structure, thermodynamics and transport of three potential electrolytes, N-methyl-N-butylpyrrolidinium bis(trifluoromethylsufonyl)imide ([pyr14][TFSI]), N- methyl-N-propylpyrrolidinium bis(fluorosufonyl)imide ([pyr13][FSI]), and 1-ethyl-3-- methylimidazolium boron tetrafluoride ([EMIM][BF4]), as a function of Li (-) salt concentration and temperature. Structurally, Li(+) is shown to be solvated by three anion neighbors in [pyr14][TFSI] and four anion neighbors in both [pyr13][FSI] and [EMIM][BF4], and at all levels of xLi we find the presence of lithium aggregates. Furthermore, the computed density, diffusion, viscosity, and ionic conductivity show excellent agreement with experimental data. While the diffusion and viscosity exhibit a systematic decrease and increase, respectively, with increasing xLi, the contribution of Li(+) to ionic conductivity increases until reaching a saturation doping level of xLi 0.10. Comparatively, the Li(+) conductivity of [pyr14][TFSI] is an order of magnitude lower than that of the other liquids, which range between 0.1-0.3 mScm. The differences in Li(+) transport are reflected in the residence times of Li(+) with the anions, which are revealed to be much larger for [pyr14][TFSI] (up to 100 ns at the highest doping levels) than in either [EMIM][BF4] or [pyr13][FSI]. Finally, we comment on the relative kinetics of Li(+) transport in each liquid and we present strong evidence for transport through anion exchange (hopping) as opposed to the net motion of Li(+) with its solvation shell (vehicular).

molecular dynamics↗
DOE OSTI2025

Numerical assessment of triply periodic minimal surfaces for direct air capture of carbon dioxide

Direct air capture (DAC) systems often consist of packing material wetted by a capture fluid that reacts with CO 2 in the airstream. The efficiency of the contactor is determined by a complex relationship of fluid dynamics, heat and mass transfer, contactor geometry, and chemical properties. The efficiency of the contactor must be balanced with other factors, primarily pressure drop through the system. Triply periodic minimal surfaces (TPMS) are a class of differential surfaces that have been explored in multiple engineering applications and have been shown to exhibit excellent performance when used in heat exchangers. Their tortuous path provides a high surface-to-volume ratio and favorable trade-off between contact area and pressure drop. In this work, a gyroid-type TPMS contactor was evaluated using computational fluid dynamics for a variety of geometric parameters to explore the potential benefit of TPMS shapes for DAC applications. A thin-film model was employed to model the flow and distribution of the capture solvent, allowing efficient simulations of TPMS structures at scale by eliminating the need for a computationally intensive interface capturing method. A liquid-gas mass transfer model was implemented in the commercial software STAR-CCM+ and used to predict the CO 2 capture efficiency and study the trade-off between capture performance and pressure drop through analysis of capture rates, mass transfer coefficients, and other relevant variables. TPMS contactors with a variety of geometric parameters and two capture solvent options were investigated to determine the effect of design choices on the operational performance of DAC systems. In conclusion, results showed that while contactor geometry is the dominant factor in efficiency and pressure drop, the physiochemical properties of the solvent are an important secondary influence on the contactor performance.

CFD↗
DOE OSTI2019

Decoding Liquid Crystal Oligomer Phase Transitions: Toward Molecularly Engineered Shape Changing Materials

This work details an integrated investigation of liquid crystal (LC) oligomers that combines experiments and molecular dynamics simulations to obtain a detailed understanding of the molecular structure of LC oligomers and the mechanism underlying their phase transition temperatures. We synthesized and characterized a series of LC oligomers prepared from different lengths of methylene spacers in the reactive LC monomers and n-alkylamine chain extenders via the aza -Michael addition reaction. In parallel, we performed isothermal–isobaric ( NPT ) ensemble coarse-grained molecular dynamics (CG-MD) simulation of analogue mesogens that are connected to flexible spacers and extenders at varying temperatures, spacer lengths, and extender lengths. This approach allowed the effect of length in the flexible spacer as well as in the chain extender on the nematic–isotropic transition temperature ( T ni ) to be determined. The results showed that increasing the length of the extender decreases T ni for LC oligomers and amplifies the decrease of T ni in LC oligomers when the spacer length is short. We infer that the combination of spacer and extender changes the shape anisotropy of LC oligomers, changing the packing behavior of constituent mesogens, thus affecting their ability to transition from the isotropic to the nematic phase. Here, the detailed molecular structure–property relationships formulated enable prescribing design rules for LC oligomers geared toward molecularly engineered shape changing materials.

36 MATERIALS SCIENCE↗
NASA NTRS2014

Investigation of Structure and Transport in Li-Doped Ionic Liquid Electrolytes: [pyr14][TFSI], [pyr13][FSI] and [EMIM][BF4]

Ionic liquid electrolytes have been proposed as a means of improving the safety and cycling behavior of advanced lithium batteries; however, the properties of these electrolytes under high lithium doping are poorly understood. Here, we employ both polarizable molecular dynamics simulation and experiment to investigate the structure, thermodynamics and transport of three potential electrolytes, N-methyl-Nbutylpyrrolidinium bis(trifluoromethylsufonyl)imide ([pyr14][TFSI]), N- methyl-Npropylpyrrolidinium bis(fluorosufonyl)imide ([pyr13][FSI]), and 1-ethyl-3-- methylimidazolium boron tetrafluoride ([EMIM][BF4]), as a function of Li-salt concentration and temperature. Structurally, Li(+) is shown to be solvated by three anion neighbors in [pyr14][TFSI] and four anion neighbors in both [pyr13][FSI] and [EMIM][BF4], and at all levels of x(sub Li) we find the presence of lithium aggregates. Furthermore, the computed density, diffusion, viscosity, and ionic conductivity show excellent agreement with experimental data. While the diffusion and viscosity exhibit a systematic decrease and increase, respectively, with increasing x(sub Li), the contribution of Li(+) to ionic conductivity increases until reaching a saturation doping level of x(sub Li) is approximately 0.10. Comparatively, the Li(+) conductivity of [pyr14][TFSI] is an order of magnitude lower than that of the other liquids, which range between 0.1 - 0.3 mS/cm. The differences in Li(+) transport are reflected in the residence times of Li(+) with the anions, which are revealed to be much larger for [pyr14][TFSI] (up to 100 ns at the highest doping levels) than in either [EMIM][BF4] or [pyr13][FSI]. Finally, we comment on the relative kinetics of Li(+) transport in each liquid and we present strong evidence for transport through anion exchange (hopping) as opposed to the net motion of Li(+) with its solvation shell (vehicular).

lithium batteries↗
DOE OSTI2019

Study of non-isothermal liquid evaporation in synthetic micro-pore structures with hybrid lattice Boltzmann model

Non-isothermal liquid evaporation in micro-pore structures is studied experimentally and numerically using the lattice Boltzmann method. A hybrid thermal entropic multiple-relaxation-time multiphase lattice Boltzmann model (T-EMRT-MP LBM) is implemented and validated with experiments of droplet evaporation on a heated hydrophobic substrate. Then liquid evaporation is investigated in two specific pore structures, i.e. spiral-shaped and gradient-shaped micro-pillar cavities, referred to as SMS and GMS, respectively. In SMS, the liquid receding front follows the spiral pattern; while in GMS, the receding front moves layer by layer from the pillar rows with large pitch to the rows with small one. Both simulations agree well with experiments. Moreover, evaporative cooling effects in liquid and vapour are observed and explained with simulation results. Quantitatively, in both SMS and GMS, the change of liquid mass with time coincides with experimental measurements. The evaporation rate generally decreases slightly with time mainly because of the reduction of liquid–vapour interface. Isolated liquid films in SMS increase the evaporation rate temporarily resulting in local peaks in evaporation rate. Reynolds and capillary numbers show that the liquid internal flow is laminar and that the capillary forces are dominant resulting in menisci pinned to the pillars. Similar Péclet number is found in simulations and experiments, indicating a diffusive type of heat, liquid and vapour transport. Our numerical and experimental studies indicate a method for controlling liquid evaporation paths in micro-pore structures and maintaining high evaporation rate by specific geometry designs.

36 MATERIALS SCIENCE↗