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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

First EarthCARE CPR Observations of Dynamics–Microphysics Coupling in Marine Cold‐Air Outbreak Clouds

Marine cold-air outbreak (MCAO) clouds begin as shallow cloud streets near the ice edge and evolve into broken open-cell convection downstream. Their rapid transitions in cloud structure and microphysics are closely linked to vertical motions, but this interplay has not previously been observed from space. The EarthCARE satellite, carrying the first spaceborne Doppler radar (Cloud Profiling Radar, CPR), now measures hydrometeor sedimentation velocity and vertical air motion. Using CPR observations over the Norwegian and Barents Seas (December 2024–May 2025), we analyze MCAO and non-MCAO clouds. We find that strong updrafts in MCAO clouds coincide with larger supercooled liquid water (SLW) path (LWP). These conditions are accompanied by increases in sedimentation velocity and its vertical gradient, reflecting rapid riming-driven ice growth. In contrast, non-MCAO clouds show lower riming efficiency at similar LWP due to relatively weaker dynamical support. This emphasizes the importance of dynamics in shaping polar cloud structure and microphysics.

Kim, Jiseob [McGill Univ., Montreal, QC (Canada)] ↗
DOE OSTI2022

Mechanism of Solid Electrolyte Interphase Formation on Lithium Anodes Studied by Molecular Dynamics Simulations

We investigated the initial formation of the SEI on lithium-metal anodes in presence of several electrolyte compositions using ReaxFF, in absence of a bias potential. We observed the formation of well-defined regions with different density and oxidation states of Li. Dissolution of Li metal atoms results in an uneven distribution of Li in the liquid phase. The electrolyte decomposition mechanism was studied for various electrolyte mixtures. Overall, the reaction mechanisms observed with ReaxFF are in good agreement with previous reports using DFT. However, we found some discrepancies in specific predictions by the ReaxFF that should be addressed in future work.

25 ENERGY STORAGE↗
DOE OSTI2026

Reduced-Order CFD Modeling to Support Waste Loading Optimization in Hanford WTP Vitrification

The U.S. DOE Hanford Site stores over 56 million gallons of radioactive liquid tank waste that must be treated and immobilized for long-term disposal The Waste Treatment and Immobilization Plant (WTP) will vitrify this waste by feeding it into Joule-heated melters, where it is incorporated into a stable borosilicate glass Computational fluid dynamics (CFD) simulations of glass melters can provide insight into the maximum achievable waste loading under varying melter operating conditions Fully resolved VOF multiphase simulations were used as the reference model to capture bubble-driven convection in the melter, including bubble formation, rise behavior, and induced glass melt circulation Effective bubble column diameter and rise velocity were extracted from the resolved simulations, compared with empirical correlations, and refit across relevant viscosity and gas flow rate conditions Explicit gas–liquid interface tracking was replaced with a single-phase momentum source term model, enabling faster steady-state CFD simulations while preserving the dominant hydrodynamic effects of bubbling New empirical correlations were developed for effective bubble column diameter and bubble rise velocity by fitting resolved simulation data across expected melter viscosity and gas flow rate ranges, providing improved inputs for the momentum source term model compared with existing literature correlations The momentum source term model reduced fluid-domain mesh size by 89% and achieved an 8.4× computational speedup relative to resolved bubbling simulations The validated momentum source term approach enables prediction of process-relevant heat transfer behavior in the integrated melter model, including heat transfer from the molten glass to the cold cap, plenum, refractory walls, and surrounding structural regions under varying melter operating conditions

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2023

Synchrotron radiography of Richtmyer–Meshkov instability driven by exploding wire arrays

We present a new technique for the investigation of shock-driven hydrodynamic phenomena in gases, liquids, and solids in arbitrary geometries. The technique consists of a pulsed power-driven resistive wire array explosion in combination with multi-MHz synchrotron radiography. Compared to commonly used techniques, it offers multiple advantages: (1) the shockwave geometry can be shaped to the requirements of the experiment, (2) the pressure (P > 300 MPa) generated by the exploding wires enables the use of liquid and solid hydrodynamic targets with well-characterized initial conditions (ICs), (3) the multi-MHz radiography enables data acquisition to occur within a single experiment, eliminating uncertainties regarding repeatability of the ICs and subsequent dynamics, and (4) the radiographic measurements enable estimation of compression ratios from the x-ray attenuation. In addition, the use of a synchrotron x-ray source allows the hydrodynamic samples to be volumetrically characterized at a high spatial resolution with synchrotron-based microtomography. This experimental technique is demonstrated by performing a planar Richtmyer–Meshkov instability (RMI) experiment on an aerogel–water interface characterized by Atwood number A0∼−0.8 and Mach number M∼1.5. The qualitative and quantitative features of the experiment are discussed, including the energy deposition into the exploding wires, shockwave generation, compression of the interface, startup phase of the instability, and asymptotic growth consistent with Richtmyer's impulsive theory. Additional effects unique to liquids and solids—such as cavitation bubbles caused by rarefaction flows or initial jetting due to small perturbations—are observed. It is also demonstrated that the technique is not shape dependent by driving a cylindrically convergent RMI experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

3D Printing ofThermally Responsive Shape Memory LiquidCrystalline Epoxy Networks

Abstract A two-component liquid crystalline epoxy network (LCEN) with shape memory behavior was developed and evaluated as a candidate material for 3D printing. The cure kinetics of the uncured material and the shape memory properties of the cured LCEN were investigated by using parallel plate rheology and dynamic mechanical analysis, respectively. A commercially available fumed silica additive was introduced to the neat, uncured material to improve the rheological properties for 3D printing. The addition of fumed silica was found to increase the yield stress, shear-thinning behavior, and toughness of the uncured epoxy ink. Polarized light microscopy, differential scanning calorimetry, and wide-angle X-ray scattering measurements between the neat and additive-modified LCEN suggested a reduction in liquid crystalline alignment in the modified LCEN, owing to interactions between crystalline domains and fumed silica, which in turn influenced the mechanical behavior. Overall, the additive was found to be successful in preserving the shape memory properties of LCEN while improving its printability.

Chemistry↗
DOE OSTI2022

Constant-potential molecular dynamics simulations of molten salt double layers for FLiBe and FLiNaK

We report the results of constant-potential molecular dynamics simulations of the double layer interface between molten 2LiF–BeF 2 (FLiBe) and 23LiF–6NaF–21KF (FLiNaK) fluoride mixtures and idealized solid electrodes. Employing methods similar to those used in studies of chloride double layers, we compute the structure and differential capacitance of molten fluoride electric double layers as a function of applied voltage. The role of molten salt structure is probed through comparisons between FLiBe and FLiNaK, which serve as models for strong and weak associate-forming salts, respectively. In FLiBe, screening involves changes in Be–F–Be angles and alignment of the oligomers parallel to the electrode, while in FLiNaK, the electric field is screened mainly by rearrangement of individual ions, predominantly the polarizable potassium cation.

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

Boron nitride nanotube precursor formation during high-temperature synthesis: kinetic and thermodynamic modelling

In this work, we performed integrated modelling of the chemical pathways of formation for boron nitride nanotube (BNNT) precursors during high-temperature synthesis in a B/N 2 mixture. Integrated modelling includes quantum chemistry, Quantum–classical molecular dynamics, thermodynamic modelling, and kinetic approaches. We demonstrate that BN compounds are formed via the interaction of molecular nitrogen with small boron clusters, rather than through interactions with less reactive liquid boron. (This process can also be described as N 2 molecule fixation.) Liquid boron evaporates to produce these boron clusters (B m with m ≤ 5), which are subsequently converted into B m N n chains. The production of such chains is crucial to the growth of BNNTs because these chains form the building blocks of bigger and longer BN chains and rings, which are in turn the building blocks of fullborenes and BNNTs. Additionally, kinetic modelling revealed that B 4 N 4 and B 5 N 4 species in particular play a major role in the N 2 molecule fixation process. Furthermore, the formation of these species via reactions with B 4 and B 5 clusters is not adequately described under the assumption of thermodynamic equilibrium, as is demonstrated in our kinetic modelling. Thus, the accumulation of both B 4 N 4 and B 5 N 4 depends on the background gas pressure and the gas cooling rate. Long BN chains and rings, which are precursors of the fullborene and BNNT growth, form via self-assembly of components B 4 N 4 and B 5 N 4 . Our modelling results—particularly the increased densities of B 4 N 4 and B 5 N 4 species at higher gas pressures—explain the experimentally observed effect of gas pressure on the yield of high-quality BNNTs. The catalytic role of hydrogen was also studied; it is shown that HBNH molecules can be the main precursor of BNNT synthesis in the presence of hydrogen.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Extracted Water Induces Concentration Fluctuations in Model Ternary Liquid–Liquid Extraction System

Organic phase aggregation and phase splitting are two important (and connected) phenomena in liquid–liquid extraction of metal ions. Previously, we demonstrated how structure in binary mixtures of extractant and diluent originates from concentration fluctuations associated with the so-called “third-phase formation” liquid–liquid phase instability. For more complex systems, we sought to understand how extracted aqueous components affect these concentration fluctuations. Here, in this study, we use small-angle X-ray scattering and molecular dynamics simulations to investigate organic phase structural changes upon contact with water in a model ternary system of water/alkane/extractant. Critical exponents were determined from scattering with temperature variation and were consistent with mean field values. Then, by controlling the water content in the organic phase through direct dissolution of fixed quantities at constant temperature, we relate the observed increase in the fluctuation length scale to changes to the spinodal temperature of the third-phase formation phase boundary. We find a roughly linear relationship between water concentration and the spinodal temperature and a similar slope in that relationship for different linear alkane diluents. This suggests that the mechanism by which water extraction impacts the organic phase structure is to enhance concentration fluctuations through the introduction of additional polar extractant-water complexes, which increases the spinodal temperature and thereby decreases the reduced temperature, i.e., the distance to the instability. This explanation for how solutes induce aggregation in liquid–liquid extraction organic phases may extend to more complex systems, explaining the large increases in correlation lengths reported under high acid and metal loading.

Bonnett, Brittany L. [Argonne National Laboratory ↗
DOE OSTI2025

Flow instabilities in helical-coil steam generators for small modular reactors: A review

Here, this study covers the research and discoveries in two-phase flow-boiling instabilities available in the literature—specifically for a helical-coil steam generator (HCSG), including experimental findings, theoretical research, computational models, and system code analyses—supporting research and development of representative small modular reactors (SMRs). Like other new and advanced reactor systems, water-cooled SMRs require experimental data from both integral and separate thermal-hydraulics test facilities for the verification and validation (V&V) of the computational models and computer codes in order to design and obtain regulatory approval. The complex dynamics of two-phase flow-boiling instabilities includes flow regimes physics phenomena, flow-channel geometries, heat-transfer behavior, and interactions among the solid–liquid-gas within the system boundary, all of which are pivotal for understanding the design and operational challenges of SMRs. This study focuses on identifying the relevant knowledge gaps on boiling instabilities—specifically for a HCSG—and provides insights about future research direction optimizing the transport of thermal energy, mass-flow rates, and boundary conditions that ensure the adequate heat-transfer performance, operational stability, and safety associated with SMR systems.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2020

Temperature Dependence of Short and Intermediate Range Order in Molten MgCl 2 and Its Mixture with KCl

This article addresses the non-Debye-Waller temperature behavior in the intermediate range order for molten MgCl 2 and its mixtures with KCl from a theory, Molecular Dynamics and experimental X-ray scattering perspective and puts these findings in the context of discussions and controversies extending at least four decades. We find that these liquids are defined by two structural motifs. The first motif is associated with chains of positive-negative charge alternation; the second motif, which results in a prepeak in the structure function S(q), is associated with the interaction of Mg 2+ and Cl - ions that do not belong to the same charge alternation chain or aggregate. Our complementary X-ray scattering and computational results provide a quantitative explanation for the increase in intensity of the prepeak with temperature as opposed to the behavior of other peaks following normal Debye-Waller behavior. We note temperature has opposite effects on the prevalence of each of the two structural motifs, and the enhancement of one pattern appears to be at the detriment of the other. Whereas the intensity in S(q) associated with the charge alternation motif is diminished at higher temperature the opposite is true for the prepeak associated with the intermediate range order due to the second structural motif.

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

Direct observation of kink evolution due to Hund’s coupling on approach to metal-insulator transition in NiS 2- x Se x

Understanding characteristic energy scales is a fundamentally important issue in the study of strongly correlated systems. In multiband systems, an energy scale is affected not only by the effective Coulomb interaction but also by the Hund’s coupling. Direct observation of such energy scale has been elusive so far in spite of extensive studies. Here, we report the observation of a kink structure in the low energy dispersion of NiS 2- x Se x and its characteristic evolution with x , by using angle resolved photoemission spectroscopy. Dynamical mean field theory calculation combined with density functional theory confirms that this kink originates from Hund’s coupling. We find that the abrupt deviation from the Fermi liquid behavior in the electron self-energy results in the kink feature at low energy scale and that the kink is directly related to the coherence-incoherence crossover temperature scale. Our results mark the direct observation of the evolution of the characteristic temperature scale via kink features in the spectral function, which is the hallmark of Hund’s physics in the multiorbital system.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Triangular lattice quantum dimer model with variable dimer density

Abstract Quantum dimer models are known to host topological quantum spin liquid phases, and it has recently become possible to simulate such models with Rydberg atoms trapped in arrays of optical tweezers. Here, we present large-scale quantum Monte Carlo simulation results on an extension of the triangular lattice quantum dimer model with terms in the Hamiltonian annihilating and creating single dimers. We find distinct odd and even $${{\mathbb{Z}}}_{2}$$ Z 2 spin liquids, along with several phases with no topological order: a staggered crystal, a nematic phase, and a trivial symmetric phase with no obvious broken symmetry. We also present dynamic spectra of the phases, and note implications for experiments on Rydberg atoms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A blueprint for biomolecular condensation driven by bacterial microcompartment encapsulation peptides

Bacterial microcompartments are protein organelles with diverse metabolic capabilities. Their functional diversity is determined by an enzymatic core that is sequestered within a structurally conserved protein shell architecture. Segregation of protein cargo into the bacterial microcompartment is enabled by encapsulation peptides, which are short helical domains fused to core proteins through a disordered linker. Here, we investigate how encapsulation peptides drive multicomponent cargo assembly into biomolecular condensates. In vitro experiments supported by molecular dynamics simulations demonstrate the importance of both conserved hydrophobic packing and electrostatic interactions in stabilizing trimeric encapsulation peptide bundles. Topological rearrangements of encapsulation peptide domains can drive programmable liquid- or gel-like partitioning in vitro and in vivo. This partitioning is found to be encapsulation peptide-specific, modular, and can co-assemble at least three fluorescent reporters. In summary, we describe the molecular features necessary to drive biomolecular condensation using a widespread peptide tag. This work can serve as a blueprint for implementing encapsulation peptide biotechnology across diverse applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Universal fluctuating regime in triangular chromate antiferromagnets

In this work, we report x-ray diffraction, magnetic susceptibility, heat capacity, 1 H nuclear magnetic resonance (NMR), and muon spin relaxation (μSR) measurements, as well as density-functional band-structure calculations for the frustrated S=3/2 triangular lattice Heisenberg antiferromagnet (TLHAF) α-HCrO 2 (trigonal, space group: $R\bar{3}m$). This compound undergoes a clear magnetic transition at T N ≃22.5 K, as seen from the drop in the muon paramagnetic fraction and concurrent anomalies in the magnetic susceptibility and specific heat capacity. Local probes (NMR and μSR) reveal a broad regime with slow fluctuations down to 0.7T N , this temperature corresponding to the maximum in the μSR relaxation rate and in the NMR wipe-out. From the comparison with NaCrO 2 and α-KCrO 2 , the fluctuating regime and slow dynamics below T N appear to be hallmarks of the TLHAF with ABC stacking. We discuss the role of interlayer frustration, which may have impacted recent spin-liquid candidates with triangular geometry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Wettability Effects on Falling Film Flow and Heat Transfer Over Horizontal Tubes in Jet Flow Mode

The performance of a falling-film heat exchanger is strongly linked to the surface characteristics and the heat transfer processes that take place over the tubes. The primary aim of this numerical study is to characterize the influence of surface wettability on the film flow behavior and its associated surface heat transfer in the jet-flow mode. Volume of fluid (VOF) based simulations are carried out for horizontal tubes with different surface wettabilities. The wettability of the tube surfaces is represented using the Kistler's dynamic contact angle model. Surface wettability effects ranging from superhydrophilic to superhydrophobic are studied by varying the equilibrium contact angle from 2 deg to 175 deg. Two different liquid mass flow rates of 0.06 and 0.18 kg/m-s corresponding to the inline and staggered jet flow modes are studied. Results are presented in terms of the liquid film thickness, the contact areas between the different phases (solid–liquid and liquid–air), and the heat transfer coefficient or Nusselt number. The resistance imposed by the increasing contact angles inhibits the extent of the liquid spreading over the tube surface, and this, in turn, influences the liquid film thickness, and the wetted area of the tube surface. A significant decrement in the heat transfer rate from the tube surfaces was observed as the equilibrium contact angle increased from 2 deg to 175 deg. The local distributions of the Nusselt number over the tube surface are strongly influenced by the flow recirculation in the liquid bulk.

Engineering↗
DOE OSTI2024

Assessing Critical Conditions for Scour Near Obstructions using Bed Shear, Particle Onset of Motion Balances, and CFD-DEM Modeling of Granular Beds

Computational Fluid Dynamics combined with a Discrete Element Method is one of the computational methods that can be used to model multiphase flows. In this method various phases, gas and liquid or solid, are present in the same computational domain. The local averaged Navier–Stokes equations determine the flow of the continuous phase fluid and are solved using the traditional CFD finite volume approach. DEM is based on a Lagrangian formulation, which solves the equations of motion, expressed in ordinary differential equations, for representative particles as they move in space and time. The interactions between the continuous fluid phase and discrete solid phase are modeled with the use of Newton’s laws of motion via drag force. The particles interact with each other and with the boundaries of the fluid continuum, and the resulting contact forces are included in the equations of motion. The properties of solid particles and boundaries are treated as elastic bodies, with specified density, elastic modulus, and Poisson’s ratio. Particle shapes may vary from single spherical particles to more complex-shaped composite particles. The particles may be introduced into the domain by random or structured injection at a point, surface, or volume, depending on the application. More details on the formulation can be found in the Simcenter STAR-CCM+ User’s Manual and OpenFOAM website.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Filament Formation and Melt Spinning of Coal-Based Mesophase Pitch for Carbon Fiber Production

High-performance carbon fibers excel as high specific strength and modulus materials, and are utilized in composites applications ranging across aerospace, automotive, energy, infrastructure, and sports equipment sectors. Compared to PAN-based carbon fiber, mesophase pitch-based carbon fiber has lower tensile strength but provides a higher modulus, higher thermal conductivity, and a lower cost potential. However, challenges in stable, continuous melt spinning processing are a serious limiting factor. Unlike typical melt spinning of long, linear chain polymers, mesophase pitch is a comprised of relatively shorter polycyclic aromatic hydrocarbons which form a liquid crystal. In its ‘green’ or ‘as-spun’ state, the fibers are very fragile. Moreover, its temperature of processing is quite high, often approaching 400ºC. Complex flow dynamics combined with short length and time scale heat transfer render its melt spinning a formidable processing challenge. Improved understanding of the root causes for nascent filament breakage and spinning instabilities are needed. This work aims to determine fundamental phenomena that govern mesophase pitch filament formation by comparison of required draw force and uninterrupted spinning minutes with applied draw down ratio. The effect of perturbations on spinning stability are analyzed using capillary rheometry, microscopy of elongating filaments, and filtration.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Phonons reveal coupled cholesterol-lipid dynamics in ternary membranes

Experimental studies of collective dynamics in lipid bilayers have been challenging due to the energy resolution required to observe these low-energy phonon-like modes. However, inelastic X-ray scattering (IXS) measurements — a technique for probing vibrations in soft and biological materials — are now possible with sub-meV resolution, permitting direct observation of low energy, phonon-like modes in lipid membranes. Here, in this study, IXS measurements with sub-meV energy resolution reveal a low-energy optic-like phonon mode at roughly 3 meV in the liquid-ordered (L o ) and liquid-disordered (L d ) phases of a ternary lipid mixture. This mode is only observed experimentally at momentum transfers greater than 5 nm -1 in the L o system. A similar gapped mode is also observed in all-atom molecular dynamics (MD) simulations of the same mixture, indicating that the simulations accurately represnt the fast, collective dynamics in the L o phase. It’s optical nature and the Q range of the gap together suggest that the observed mode is due to the coupled motion of cholesterol-lipid pairs, separated by several hydrocarbon chains within the membrane plane. Analysis of the simulations provide molecular insight into the origin of the mode in transient, nanoscale substructures of hexagonally packed hydrocarbon chains. This nanoscale hexagonal packing was previously reported based on molecular dynamics simulations and later by NMR measurements. Here, however, the integration of IXS and MD simulations identifies a new signature of the L o substructure in the collective lipid dynamics, thanks to the recent confluence of IXS sensitivity and MD simulation capabilities.

36 MATERIALS SCIENCE↗