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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Assessing Melting and Solid–Solid Transition Properties of Choline Chloride via Molecular Dynamics Simulations

Choline chloride (ChCl) is used extensively as a hydrogen bond donor in deep eutectic solvents (DESs). However, determining its melting properties experimentally is challenging due to decomposition upon melting, leading to widely varying literature values. Accurate melting properties are crucial for understanding the solid–liquid phase behavior of ChCl-containing DESs. Here, we employ molecular dynamics simulations to compute the phase transitions of ChCl, testing a variety of atomistic force fields. We find that the results are sensitive to the choice of force field, but a melting temperature of 627 K and a melting enthalpy of 7.8 kJ/mol seem most reasonable, in good agreement with some literature values. Furthermore, we suggest these as the likely melting properties of ChCl, though the results are tentative due to limited experimental data for the liquid ChCl phase.

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

Electric-Field-Induced Spatially Dynamic Heterogeneity of Solvent Motion and Cation Transference in Electrolytes

While electric fields primarily result in migration of charged species in electrolytic solutions, the solutions are dynamically heterogeneous. Solvent molecules within the solvation shells of the cation will be dragged by the field while free solvent molecules will not. Here we combine electrophoretic NMR measurements of ion and solvent velocities under applied electric fields with molecular dynamics simulations to interrogate different solvation motifs in a model liquid electrolyte. Measured values of the cation transference number $(t^0_+)$ agree quantitatively with simulation-based predictions over a range of electrolyte concentrations. Solvent-cation interactions strongly influence the concentration-dependent behavior of $t^0_+$. We identify a critical concentration at which most of the solvent molecules lie within solvation shells of the cations. The dynamic heterogeneity of solvent molecules is minimized at this concentration where $t^0_+$ is approximately equal to zero.

42 ENGINEERING↗
DOE OSTI2021

Density of Fe-Ni-C Liquids at High Pressures and Implications for Liquid Cores of Earth and the Moon

The presence of light elements in the metallic cores of the Earth, the Moon, and other rocky planetary bodies has been widely proposed. Carbon is among the top candidates in light of its high cosmic abundance, siderophile nature, and ubiquity in iron meteorites. It is, however, still controversial whether carbon-rich core compositional models can account for the seismic velocity observations within the Earth and lunar cores. In this paper we report the density and elasticity of Fe 90 Ni 10 -3 wt.% C and Fe 90 Ni 10 -5 wt.% C liquid alloys using synchrotron-based X-ray absorption experiments and first-principles molecular dynamics simulations. Our results show that alloying of 3 wt.% and 5 wt.% C lowers the density of Fe 90 Ni 10 liquid by ~2.9-3.1% at 2 GPa, and ~3.4-3.6% at 9 GPa. More intriguingly, our experiments and simulations both demonstrate that the bulk moduli of the Fe-Ni-C liquids are similar to or slightly higher than those of Fe-Ni liquids. Thus, the calculated compressional velocity ($v_p$) of Fe-Ni-C liquids are higher than that of pure Fe-Ni alloy, promoting carbon as a possible candidate to explain the elevated v p in the Earth's outer core. However, the values and slopes of both density and v p of the studied two Fe-Ni-C liquids do not match the outer core seismic models, suggesting that carbon may not be the sole principal light element in Earth's outer core. The high $v_p$ of Fe-Ni-C liquids does not match the presumptive $v_p$ of the lunar outer core well, indicating that carbon is less likely to be its dominant light element.

58 GEOSCIENCES↗
DOE OSTI2020

Liquid-Vapor Coexistence and Critical Point of Mg 2 SiO 4 From Ab Initio Simulations

Hypervelocity impact-driven vaporization is characteristic of late-stage planet formation. Yet the behavior and properties of liquid-vapor mixtures of planetary materials of interest are typically unknown. Multiphase equations of state used in hydrodynamic simulations of planet impacts therefore lack reliable data for this important phenomenon. Here, we present the first constraints on the liquid-vapor critical point and coexistence phase boundary of Mg2SiO4 computed from ab initio molecular dynamics simulations. We found that the vapor is depleted in magnesium and enriched in silica and oxygen, while the coexisting liquid is enriched in magnesium and depleted in oxygen, from which we infer vaporization is incongruent. The critical point was estimated from an equation of state fit to the data. The results are in line with recent calculations of MgSiO3 and together confirm that extant multiphase equation of state (EOS) models used in planetary accretion modeling significantly underestimate the amount of supercritical material postimpact.

58 GEOSCIENCES↗
DOE OSTI2023

Flow field characteristics of a 127-pin rod bundle with hexagonal spacer grids

A determination of nominal flow phenomena in liquid metal fast reactor (LMFR) fuel assemblies is critical toward generation-IV reactor development. Axially positioned spacer grids are used to maintain the geometry of hexagonal rod bundles and simultaneously introduce perturbations in the flow. Three-dimensional (3D) printed asymmetric honeycomb spacer grids were installed in a prototypical 127-pin LMFR fuel assembly model to study complex fluid dynamics interactions induced by the spacer grid and rods. To characterize flow dynamics in this intricate geometry, time-resolved particle image velocimetry (TR-PIV) using the matched-index-of-refraction method was employed to obtain non-intrusive velocity measurements for three axial planes (one near-wall and two interior planes) at a Reynolds number of 6000. The statistical TR-PIV results compared sub-channel-dependent normalized time-averaged velocity, velocity fluctuations, Reynolds stress, vorticity, and turbulence kinetic energy distributions. TR-PIV line profiles characterized downstream spacer grid flow dynamics. Two-point spatial and spatial–temporal cross-correlation fields revealed local coherent structures and quantified convection velocities of traveling vortices. Spatial–temporal decomposition using dynamic mode decomposition (DMD) applied to the near-wall vorticity fields extracted turbulent structures and flow instabilities in the wake region of the spacer grid, along with their decay and frequency rates. Reduced-order velocity fields from DMD reconstructions identified the most energy-containing coherent structures persistent in the near-wall region. This research provides experimental data sets and analyses of flow behavior in rod bundles with hexagonal spacer grids. In conclusion, the results are critical toward LMFR design and geometry optimization, crucial for the validation of computational fluid dynamics and reduced-order flow models.

42 ENGINEERING↗
DOE OSTI2025

Low-energy edge signatures of the Kitaev spin liquid

Recent experimental work indicates that the Kitaev spin liquid may be realizable close to zero magnetic field in exfoliated 𝛼−RuCl 3 flakes, thus providing a more versatile setting for studying non-Abelian anyons. Here, we propose a robust nanoscale signature of the Kitaev spin liquid that results from its edge states and manifests in the low-energy spin dynamics. In particular, we highlight a singular peak in the dynamical spin structure factor of a zigzag edge whose energy scales linearly with only one component of the magnetic field. This sharp feature in the local spin dynamics directly reflects the bond-directional Kitaev spin interactions and, more generally, the projective symmetries of the Kitaev spin liquid. In conclusion, we demonstrate that our proposed signature survives in the presence of edge disorder as well as non-Kitaev interactions, and provide detailed guidelines for experimentally observing it using inelastic electron tunneling spectroscopy and color-center relaxometry.

Color centers↗
DOE OSTI2024

Impulse-driven release of gas-encapsulated drops

Gas-encapsulated drops, much like antibubbles, are drops enclosed in a bubble within a liquid. They show potential as payload carriers in fluid transport and mixing techniques where sound waves can be leveraged to induce the collapse of the gas core and the subsequent release of the drop. Here, the interaction of millimetre-sized gas-encapsulated drops with impulsive laser-induced shock waves is investigated to gain fundamental insights into the release process. Experimental synchrotron X-ray phase contrast imaging, which allows the drop dynamics to be visualised inside the encapsulating bubble, is complemented by numerical simulations to study the intricate physics at play. Three drop dynamical release regimes are discovered, namely the drop impact , partial deposition and jet impact regimes. The regime type is mainly dependent on the shape of the bubble interface impacting the drop and the associated Weber and Reynolds numbers. The drop dynamics of the drop impact and partial deposition regimes show similarities with the canonical configuration of drops impacting flat liquid surfaces, whereas the jet impact regime resembles binary drop collisions, which allows existing scaling laws to be applied to describe the underlying processes. The release of the drop is investigated numerically. The time evolution of the drop dissemination within the surrounding liquid discloses enhanced mixing for dynamics involving high Weber and Reynolds numbers such as the drop impact and jet impact regimes.

42 ENGINEERING↗
DOE OSTI2020

Infrared absorption spectroscopy of dynamically compressed water

Presentation describes development of a time-stretch spectroscopy diagnostic using dispersed supercontinuum laser pulses. The technique is applied to near-infrared measurements of liquid water during multiple shock compression.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Overcoming time and complexity limitations in molecular dynamics investigations of equilibrium melting

Abstract A hybrid Monte-Carlo molecular-dynamics method for determining solidus and liquidus compositions in multicomponent systems is presented that overcomes both the time limitations in conventional molecular dynamics that prevent the evolution of distinct solid and liquid compositions via diffusion and the complexity challenge that prevents use of thermodynamic assessment in systems of many components. This hybrid method is validated in the Cu–Ni system against an independent assessment of solidus and liquidus compositions based on the regular solution model. Strategies for efficient mapping of different phase diagrams, based on the thermodynamic parameter T 0 , the temperature at which two phases of the composition X 0 have equal free energies, are presented and then demonstrated for the copper-nickel fully miscible system and the gold–silicon eutectic system. A calculation of the solidus and liquidus sampled during the equilibrium melting of equiatomic CrMnFeCoNi is performed, indicating that this method has potential to be extended to the study of many component alloys.

Au-Si↗
DOE OSTI2023

A sui generis whipping-instability-based self-sequencing multi-monodisperse 2D spray from an anisotropic microfluidic liquid jet device

Well-defined aerosols pave the way for versatile basic and applied research. Here, we demonstrate a unique whipping instability that generates from a high-aspect-ratio microfluidic device resulting in a unique steady-state gas-focused whipping jet (WJ) without any need for electrification. This WJ device emanates a multi-monodisperse whipping spray jet with a two-dimensional (2D) profile. We demonstrate this phenomenon based on various fluidic parameters theoretically and experimentally. The 2D WJ’s unique behavior is derived using analytical fluid dynamics to explain jet diameter, whipping regime, and spreading angle. The phenomenon is further characterized experimentally by measuring the angle with respect to the flow rate, the distances between droplets, the droplet shapes, and the reproducibility of these parameters. We also explain the precise fabrication of such inexpensive devices. Lastly, we highlight these devices’ potential use as sample environments in versatile applications ranging from cryoelectron microscopy over mass spectrometry to drug formulation and structural studies at X-ray free-electron lasers.

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

Development, Verification, and Validation of an OpenFOAM-Based Solver for Modeling Inertial Fusion Energy Chambers

Our work seeks to introduce a computational tool tailored to the physics of inertial fusion energy chambers, in particular, those concepts based on thick liquid walls. In this approach, the structural materials are protected by several neutron mean-free-paths of renewable liquid and thus will be able to survive much longer than un-shielded walls, with virtually all structures lasting for the life of the plant and enabling the use of commercially available and qualified materials. The OpenFOAM-based solver named rhoCentralFoam has been used as a starting point. rhoCentralFoam belongs to the standard OpenFOAM solver toolset. It is a high-speed, explicit compressible flow solver with shock-capturing capability. While the main features have been retained, the solver had to be restructured to make use of tabular data for equations of states, a necessary addition to model the complex thermo-physical properties of ionized gasses. This entailed the need to change the independent state variables used by the solver, resulting in a new thermodynamic library and slightly different solution algorithm. Moreover, a radiation heat transfer model based on the P-1 approximation was added to the solver. The solver is verified against an analytical solution from the Sedov-Taylor-Neumann test problem to showcase the ability of the hydrodynamic solvers to handle strong shocks, whereas the P-1 model was verified using a simple one-dimensional problem with an analytical solution. Additionally, a validation case involving shock-wave propagation through jet array is presented, and the results are compared with experimental data from the open literature. Lastly, in order to showcase the utility of the solver for practical cases, we applied the refined solver to two representative scenarios: gas venting within the HYLIFE-II chamber and the compression of the gas following the partial ablation of the liquid wall.

Chamber dynamics↗
DOE OSTI2020

Imogolite in water: Simulating the effects of nanotube curvature on structure and dynamics

Imogolite is a fascinating inorganic nanotube that is found in nature or synthesized in a laboratory. The synthesis process is carried out in liquid media, and leads to the formation of almost monodisperse diameter nanotubes. Here we investigate, employing classical molecular dynamics simulations, the interaction of water and imogolite for nanotubes of several radii. We established that water penetrates the pores of N = 9 and larger nanotubes, and adopts a coaxial arrangement in it. Also, while water molecules can diffuse along the center of the nanotube, the molecules next to the inner imogolite walls have very low mobility. At the outer nanotube wall, an increase of water density is observed, this effect extends up to 1 nm, beyond which water properties are bulk-like. Both phenomena are affected by the imogolite curvature.

58 GEOSCIENCES↗
DOE OSTI2023

Structure and dynamics of ILs-based gel polymer electrolytes and its enhanced conductive properties with the incorporation of Al 2 O 3 nanofibers

Here, this work reports the enhanced mobility of ions in ionic liquid (IL)-based gel polymer electrolytes (GPEs) with the incorporation of Al 2 O 3 nanofibers. A combination of PVDF-HFP, EMIMTFSI and LiTFSI with 3 wt% Al 2 O 3 nanofibers has been prepared through solution casting technique. The room temperature ionic conductivity of PVDF-HFP: ILs electrolyte (45:55, weight ratio of 0.82) (GPE) is found to be 2.7 × 10 –5 S cm –1 , which increases up to 7.8 × 10 –5 S cm –1 in Al 2 O 3 containing GPE (Al-GPE). Pulsed field gradient (PFG) NMR results validate the increased ionic conductivity observed in Al-GPE. We found that the diffusivity of Li + , TFSI – and EMIM + increases when Al 2 O 3 nanofibers are well-distributed in the GPE matrix. The surface morphology and the amorphicity of GPEs are examined through SEM and XRD analyses. Lastly, the local structure of Al 2 O 3 fibers and the molecular-level interactions of ions with polymer, and their effect on the diffusivity of ions are established through solid-state NMR detecting 27 Al, 1 H, 13 C, 19 F nuclei including 2D 13 C{ 1 H} HETCOR NMR experiments. The 13 C DPMAS and CPMAS experiments highlight the dynamic heterogeneity associated with the ions that are embedded in the rigid and the mobile phase of GPEs. While some of the ionic species strongly interact with polymer chains in the rigid environment, the majority of them reside in the mobile phase and contribute to the overall increased conductivity. Most importantly, Al 2 O 3 nanofibers significantly affect the dynamics of ionic species that are present in the mobile phase between the polymer chains.

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

Exploring the formation of gold/silver nanoalloys with gas-phase synthesis and machine-learning assisted simulations

While nanoalloys are of paramount scientific and practical interest, the main processes leading to their formation are still poorly understood. Key structural features in the alloy systems, including the crystal phase, chemical ordering, and morphology, are challenging to control at the nanoscale, making it difficult to extend their use to industrial applications. Here, in this contribution, we focus on the gold/silver system that has two of the most prevalent noble metals and combine experiments with simulations to uncover the formation mechanisms at the atomic level. Nanoparticles were produced using a state-of-the-art inert-gas aggregation source and analyzed using transmission electron microscopy and energy-dispersive X-ray spectroscopy. Machine-learning-assisted molecular dynamics simulations were employed to model the crystallization process from liquid droplets to nanocrystals. Our study finds a preponderance of nanoparticles with five-fold symmetric morphology, including icosahedra and decahedra which is consistent with previous results on mono-metallic nanoparticles. However, we observed that gold atoms, rather than silver atoms, segregate at the surface of the obtained nanoparticles for all the considered alloy compositions. These segregation tendencies are in contrast to previous studies and have consequences on the crystallization dynamics and the subsequent crystal ordering. We finally showed that the underpinning of this surprising segregation dynamics is due to charge transfer and electrostatic interactions rather than surface energy considerations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Generalized approach for rapid entropy calculation of liquids and solids

We build a comprehensive methodology for the fast computation of entropy across both solid and liquid phases. The proposed method utilizes a single trajectory of molecular dynamics (MD) to facilitate the calculation of entropy, which is composed of three components. The electronic entropy is determined through the temporal average acquired from density functional theory MD simulations. The vibrational entropy, typically the predominant contributor to the total entropy, even within the liquid state, is evaluated by computing the phonon density of states via the velocity autocorrelation function. The most arduous component to quantify, the configurational entropy, is assessed by probability analysis of the local structural arrangement and atomic distribution. We illustrate, through a variety of examples, that this method is both a versatile and valid technique for characterizing the thermodynamic states of both solids and liquids. Furthermore, this method is employed to expedite the calculation of melting temperatures, demonstrating its practical utility in computational thermodynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A collective elastic fluctuation mechanism for decoupling and stretched relaxation in glassy colloidal and molecular liquids

We propose a microscopic theory for the decoupling of self-diffusion and structural relaxation in glass-forming liquids within the Elastically Collective Nonlinear Langevin Equation (ECNLE) activated dynamics framework. Our central hypothesis is that the heterogeneity relevant to this problem is static fluctuations of local density on the scale of 3–4 particle diameters and how this changes local packing correlations. These fluctuations modify the degree of dynamical cage expansion that mechanistically couples intracage large amplitude hopping and longer range collective elasticity in ECNLE theory. Decoupling only emerges in the deeply supercooled regime where the strongly temperature dependent elastic barrier becomes non-negligible relative to its noncooperative local analog. The theory makes predictions for various aspects of the decoupling phenomenon, including apparent fractional power law Stokes-Einstein behavior, that appear to be consistent with experiments and simulations on hard sphere fluids and molecular liquids. Of central importance is a microscopic connection between the barrier fluctuation variance and most probable barrier height. Sensible results are also obtained for the nonexponential stretching of a generic relaxation time correlation function and its temperature evolution. Nonuniversality can arise from the relative importance of the local and collective barriers (related to fragility) and the precise magnitude of the length scale that defines the transition from local cage to elastic physics. Comparison is made with a traplike model based on a Gaussian distribution of barriers.

Xie, Shi-Jie (ORCID:0000000333174720)↗
DOE OSTI2022

A high-dimensional neural network potential for molecular dynamics simulations of condensed phase nickel and phase transitions

A high-dimensional neural network interatomic potential was developed and used in molecular dynamics simulations of condensed phase Ni and Ni systems with liquid–solid phase coexistence. The reference data set was generated by sampling the potential energy surface over a broad temperature-pressure domain using ab initio MD simulations to train a unified potential. Excellent agreement was achieved between bulk face-centered cubic nickel thermal expansion simulations and relevant experimental data. The same potential also yields accurate structures and diffusivities in the liquid state. The phase transition between liquid and solid phases was simulated using the two-phase interface method. The predicted melting point temperature is within a few kelvins of the literature value. Here, the general methodology could be applied to describe crystals with much more complex phase behaviors.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

pH‐ and Redox‐Responsive Pickering Emulsions Based on Cellulose Nanocrystal Surfactants

Abstract Nanoparticle surfactants (NPSs), formed by using dynamic interactions between nanoparticles and complementary ligands at the liquid‐liquid interface, have emerged as “smart emulsifiers” with attributes of high emulsification efficiency, long‐term stability, and on‐demand emulsification/demulsification capabilities. However, only pH‐responsiveness can be adopted for the assembly of reported NPSs formed by electrostatic interactions. Here, we propose an alternative design strategy, by taking advantage of the ferrocenium (Fc + ) sulfate ion pair, to develop a new type of cellulose nanocrystal (CNC) surfactant. The Fc + groups are sensitive to pH, redox reagents and voltage, imparting the CNC surfactants and derived Pickering emulsions with multi‐stimuli‐responsiveness, and showing promising applications in controllable delivery, release, and biphasic biocatalysis.

Yang, Yang↗