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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Modeling of Non-Isothermal Cryogenic Fluid Sloshing

A computational fluid dynamic model was used to simulate the thermal destratification in an upright self-pressurized cryostat approximately half-filled with liquid nitrogen and subjected to forced sinusoidal lateral shaking. A full three-dimensional computational grid was used to model the tank dynamics, fluid flow and thermodynamics using the ANSYS Fluent code. A non-inertial grid was used which required the addition of momentum and energy source terms to account for the inertial forces, energy transfer and wall reaction forces produced by the shaken tank. The kinetics-based Schrage mass transfer model provided the interfacial mass transfer due to evaporation and condensation at the sloshing interface. The dynamic behavior of the sloshing interface, its amplitude and transition to different wave modes, provided insight into the fluid process at the interface. The tank pressure evolution and temperature profiles compared relatively well with the shaken cryostat experimental test data provided by the Centre National D'Etudes Spatiales.

Cryogenics↗
DOE OSTI2025

Multiphasic size-dependent growth dynamics of nanoparticle ensembles

Colloidal nanoparticles are of great interest in modern science and industry. However, the thermodynamic mechanism and dynamics of nanoparticle growth have yet to be understood. Addressing these issues, we tracked hundreds of in-situ growth trajectories of a nanoparticle ensemble using liquid-phase TEM and found that the nanoparticle growth, including coalescence, exhibits nanoparticle size-dependent multiphasic dynamics, unexplainable by current theories. Motivated by this finding, we developed a model and theory for an ensemble of growing nanoparticles, providing a unified, quantitative understanding of the time-dependent mean and fluctuation of nanoparticle size and size-dependent growth rate profiles across various nanoparticle systems and experimental conditions. Our work reveals that the chemical potential in a small nanoparticle strongly deviates from the Gibbs-Thomson equation, shedding light on how it governs the size-dependent growth dynamics of nanoparticles.

Kim, Ji-Hyun↗
NASA NTRS1996

Leakage Through a Channel Formed by a Gasket, a Sealing Surface, and a Filament Trapped Between Them

Plumbing for the transport of liquid Hydrogen or liquid Oxygen at the Kennedy Space Center (KSC) is very critical. Every piece of hardware for handling such a hazardous cryogen is subject to testing prior to installation and use. Safe, realistic testing of all such hardware is prohibitively expensive, which leads, perforce, to expidients such as: (1) lead testing with non-flammable tracer fluids (e.g, liquid nitrogen) and (2) leak testing with room temperature tracer fluids (e.g. liquid helium). Such expedients undermine the realism of the tests. If however, one could apply rational fluid dynamics methods to derive a general analytical expression with which one could relate the throughput of gaseous Helium through a given leak channel to the throughput of liquid Hydrogen through the same channel, then one could recover much of the information that one would otherwise forfeit through these expedients. These facts lead to the following questions: (1) What would be an example of a generic flaw in a gasket?; and (2) How can one calculate the flow of fluid in it? The report addresses these questions. It considers a particular leak geometry, namely one formed by a gasket, a sealing surface, and a filament trapped between them (so that the cross section of the leak channel is a flat bottomed curvilinear triangle, two sides of which are circular arcs and which has cusps on all three corners).

Russell, John↗
DOE OSTI2023

Quantum disordered ground state in the triangular-lattice magnet NaRuO 2

It has long been hoped that spin liquid states might be observed in materials that realize the triangular lattice Hubbard model. However, weak spin-orbit coupling and other small perturbations often induce conventional spin freezing or magnetic ordering. Sufficiently strong spin-orbit coupling, however, can renormalize the electronic wave function and induce anisotropic exchange interactions that promote magnetic frustration. Here we show that the cooperative interplay of spin-orbit coupling and correlation effects in the triangular lattice magnet NaRuO 2 produces an inherently fluctuating magnetic ground state. Despite the presence of a charge gap, we find that low-temperature spin excitations generate a metal-like term in the specific heat and a continuum of excitations in neutron scattering, reminiscent of spin liquid states previously found in triangular lattice organic magnets. Further cooling produces a crossover into a different, highly disordered, spin state whose dynamic spin autocorrelation function reflects persistent fluctuations. Furthermore, these findings establish NaRuO 2 as a cousin to organic, Heisenberg spin liquid compounds with a low-temperature crossover in quantum disorder.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Calculation of the detonation state of HN 3 with quantum accuracy

HN 3 is a unique liquid energetic material that exhibits ultrafast detonation chemistry and a transition to metallic states during detonation. We combine the Chebyshev interaction model for efficient simulation (ChIMES) many-body reactive force field and the extended-Lagrangian multiscale shock technique molecular dynamics method to calculate the detonation properties of HN 3 with the accuracy of Kohn–Sham density-functional theory. ChIMES is based on a Chebyshev polynomial expansion and can accurately reproduce density-functional theory molecular dynamics (DFT-MD) simulations for a wide range of unreactive and decomposition conditions of liquid HN 3 . We show that addition of random displacement configurations and the energies of gas-phase equilibrium products in the training set allows ChIMES to efficiently explore the complex potential energy surface. Schemes for selecting force field parameters and the inclusion of stress tensor and energy data in the training set are examined. Structural and dynamical properties and chemistry predictions for the resulting models are benchmarked against DFT-MD. We demonstrate that the inclusion of explicit four-body energy terms is necessary to capture the potential energy surface across a wide range of conditions. Our results generally retain the accuracy of DFT-MD while yielding a high degree of computational efficiency, allowing simulations to approach orders of magnitude larger time and spatial scales. The techniques and recipes for MD model creation we present allow for direct simulation of nanosecond shock compression experiments and calculation of the detonation properties of materials with the accuracy of Kohn–Sham density-functional theory.

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

Elevating density functional theory to chemical accuracy for water simulations through a density-corrected many-body formalism

Abstract Density functional theory (DFT) has been extensively used to model the properties of water. Albeit maintaining a good balance between accuracy and efficiency, no density functional has so far achieved the degree of accuracy necessary to correctly predict the properties of water across the entire phase diagram. Here, we present density-corrected SCAN (DC-SCAN) calculations for water which, minimizing density-driven errors, elevate the accuracy of the SCAN functional to that of “gold standard” coupled-cluster theory. Building upon the accuracy of DC-SCAN within a many-body formalism, we introduce a data-driven many-body potential energy function, MB-SCAN(DC), that quantitatively reproduces coupled cluster reference values for interaction, binding, and individual many-body energies of water clusters. Importantly, molecular dynamics simulations carried out with MB-SCAN(DC) also reproduce the properties of liquid water, which thus demonstrates that MB-SCAN(DC) is effectively the first DFT-based model that correctly describes water from the gas to the liquid phase.

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

To investigate the reciprocal interactions between dielectric barrier discharges and different biological targets in cell culture medium

The objective of this project was to investigate the dynamic interactions between non-thermal plasma (NTP) and biological targets suspended in a physiological liquid medium. The experiments under this collaborative study were designed to combine theoretical and experimental approaches to test the hypothesis that when living cells (and the growth medium surrounding them) encounter plasma, they reciprocally and dynamically modify each other.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1995

On the breakup of viscous liquid threads

A one-dimensional model evolution equation is used to describe the nonlinear dynamics that can lead to the breakup of a cylindrical thread of Newtonian fluid when capillary forces drive the motion. The model is derived from the Stokes equations by use of rational asymptotic expansions and under a slender jet approximation. The equations are solved numerically and the jet radius is found to vanish after a finite time yielding breakup. The slender jet approximation is valid throughout the evolution leading to pinching. The model admits self-similar pinching solutions which yield symmetric shapes at breakup. These solutions are shown to be the ones selected by the initial boundary value problem, for general initial conditions. Further more, the terminal state of the model equation is shown to be identical to that predicted by a theory which looks for singular pinching solutions directly from the Stokes equations without invoking the slender jet approximation throughout the evolution. It is shown quantitatively, therefore, that the one-dimensional model gives a consistent terminal state with the jet shape being locally symmetric at breakup. The asymptotic expansion scheme is also extended to include unsteady and inerticial forces in the momentum equations to derive an evolution system modelling the breakup of Navier-Stokes jets. The model is employed in extensive simulations to compute breakup times for different initial conditions; satellite drop formation is also supported by the model and the dependence of satellite drop volumes on initial conditions is studied.

Papageorgiou, Demetrios T.↗
NASA NTRSDate not supplied

Validation of Accurate Cryogenic Fluid Vapor-Liquid Boundary Conditions Via Molecular Simulations

Vapor-liquid interfaces drive many important phenomena in cryogenic fluid management, including heat transfer, evaporation, and capillary flow. Design of cryogenic fluid systems, such as propellant storage, requires accurate predictions of fluid behavior, including evaporation rates. Many models have been proposed for heat and mass transfer at vapor liquid interfaces, but the accuracy of these models in the context of cryogenic fluids has not been performed. We use molecular dynamics simulations, which allow for nanometer scale resolution of fluid phenomena, to evaluate the accuracy of a variety of vapor-liquid boundary conditions at evaporating and condensing interfaces. We find that an anisotropic temperature distribution is a critical ingredient for accurate prediction of intensive evaporation and condensation.

Daniel Vigil↗
DOE OSTI2025

Consequences of the failure of equipartition for the p – V behavior of liquid water and the hydration free energy components of a small protein

Earlier we showed that in the molecular dynamics simulation of a rigid model of water it is necessary to use an integration time-step δ t ≤ 0.5 fs to ensure equipartition between translational and rotational modes. Here we extend that study in the NVT ensemble to NpT conditions and to an aqueous protein. We study neat liquid water with the rigid, SPC/E model and the protein BBA (PDB ID: 1FME) solvated in the rigid, TIP3P model. We examine integration time-steps ranging from 0.5 fs to 4.0 fs for various thermostat plus barostat combinations. We find that a small δ t is necessary to ensure consistent prediction of the simulation volume. Hydrogen mass repartitioning alleviates the problem somewhat, but is ineffective for the typical time-step used with this approach. The compressibility, a measure of volume fluctuations, and the dielectric constant, a measure of dipole moment fluctuations, are also seen to be sensitive to δ t . Using the mean volume estimated from the NpT simulation, we examine the electrostatic and van der Waals contribution to the hydration free energy of the protein in the NVT ensemble. These contributions are also sensitive to δ t . In going from δ t = 2 fs to δ t = 0.5 fs, the change in the net electrostatic plus van der Waals contribution to the hydration of BBA is already in excess of the folding free energy reported for this protein.

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

Molecular Dynamics Investigation of Oil Wetting on Synthetic Polymer Substrates

Understanding the interaction of polymer surfaces with nonpolar, low surface tension liquids, or whether a substrate is oleophobic versus oleophilic, is critical for applications ranging from antifingerprint coatings to oil–water separation membranes and oil spill remediation. Despite its technological importance, the molecular mechanisms underpinning polymer–oil wetting are not well characterized. Here, we employ molecular dynamics simulations to investigate the behavior of n-hexadecane in contact with chemically distinct polymer surfaces, spanning eight constitutional unit chemistries as well as amorphous and crystalline morphologies. This permits a critical examination of both thermodynamic and dynamic descriptors of oleophobicity, including oil contact angle, dewetting free energy, interfacial diffusivity, and a proposed “ghost probe energy,” which does not require explicit simulation of oil–polymer interactions. We find that the oil contact angle does not reliably distinguish oleophobic behavior across polymer chemistries, whereas other metrics provide clearer and more consistent differentiation. Analysis of the results reveals that polymer–oil wetting behavior is primarily governed by interfacial van der Waals interactions and modulated by surface flexibility and morphology. Collectively, this work establishes a computational framework for characterizing oil wetting and provides additional insight into what molecular-level factors dictate trends in oleophobicity.

lipids↗
DOE OSTI2025

Crystal nucleation rates in one-component Yukawa systems

Nucleation in the supercooled Yukawa system is relevant for addressing current challenges in understanding a range of crystallizing systems including white dwarf (WD) stars. We use both brute force and seeded molecular dynamics simulations to study homogeneous nucleation of crystals from supercooled Yukawa liquids. With our improved approach to seeded simulations, we obtain quantitative predictions of the crystal nucleation rate and cluster size distributions as a function of temperature and screening length. These quantitative results show trends towards fast nucleation with short-ranged potentials. They also indicate that for temperatures T > 0.9⁢T m , where T m is the melt temperature, classical homogeneous nucleation is too slow to initiate crystallization but transient clusters of ~100 particles should be common. As a result, we apply these general results to a typical WD model and obtain a delay of ~0.6 Gyr in the onset of crystallization that may be observable.

36 MATERIALS SCIENCE↗
NASA NTRS2016

Li-Doped Ionic Liquid Electrolytes: From Bulk Phase to Interfacial Behavior

Ionic liquids have been proposed as candidate electrolytes for high-energy density, rechargeable batteries. We present an extensive computational analysis supported by experimental comparisons of the bulk and interfacial properties of a representative set of these electrolytes as a function of Li-salt doping. We begin by investigating the bulk electrolyte using quantum chemistry and ab initio molecular dynamics to elucidate the solvation structure of Li(+). MD simulations using the polarizable force field of Borodin and coworkers were then performed, from which we obtain an array of thermodynamic and transport properties. Excellent agreement is found with experiments for diffusion, ionic conductivity, and viscosity. Combining MD simulations with electronic structure computations, we computed the electrochemical window of the electrolytes across a range of Li(+)-doping levels and comment on the role of the liquid environment. Finally, we performed a suite of simulations of these Li-doped electrolytes at ideal electrified interfaces to evaluate the differential capacitance and the equilibrium Li(+) distribution in the double layer. The magnitude of differential capacitance is in good agreement with our experiments and exhibits the characteristic camel-shaped profile. In addition, the simulations reveal Li(+) to be highly localized to the second molecular layer of the double layer, which is supported by additional computations that find this layer to be a free energy minimum with respect to Li(+) translation.

molecular dynamics↗
DOE OSTI2025

Development of a conduction-based model for analyzing frozen startup of alkali-metal heat pipes

One key area of interest in heat pipe modeling/simulation is to analyze the startup behavior of the liquid-metal heat pipes (LMHPs) from a frozen state. This so-called ‘frozen startup’ process involves a complex set of nonlinear mass and heat transport phenomena, including phase transitions from solid to liquid and vapor, multiphase interactions, microporous wick flow, and compressible vapor dynamics. The complexity of these processes makes it challenging to simulate LMHP’s frozen startup using conventional numerical methods or commercial computational fluid dynamics (CFD) software. This paper presents a simplified conduction-based modeling approach that can provide practical insights into the entire LMHP frozen startup process, while alleviating the challenges of modeling its complex physics. The theoretical foundation and physical assumptions of the proposed model are based solely on heat-conduction equation, allowing for a more tractable simulation without sacrificing essential physical accuracy. The proposed model was implemented in a commercial CFD software, and its prediction was compared with the experimental data obtained from sodium heat-pipe startup experiments. The comparison highlights the proposed model's ability to capture the transient thermal behavior of LMHP during frozen startup. This study not only validates the conduction-based frozen startup modeling method but also shows its potential as a practical and efficient tool for understanding the startup performance of the LMHP systems.

Microreactor↗