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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Dynamic compression effects of H 2 ⁡O in a dynamic diamond anvil cell: Origin of metastable ice VII and its crystal growth kinetics

We report on the structural verification of metastable ice VII solidifying in the phase space of ice VI at 1.80 GPa at room temperature. Using time-resolved (TR) x-ray diffraction and TR ruby luminescence paired with high-speed microphotography utilizing a dynamic diamond anvil cell, an initial compression rate range from 0.12 to 95.84 GPa/s was explored. The solidification pressure of metastable ice VII has a potential sigmoidal dependence upon compression rate with a turnover compression rate of ∼80 GPa/s. The preferred crystallization of ice VII in the stability field of ice VI is due to the increased nucleation rate of ice VII over ice VI at 1.77 GPa that is driven by the surface energy difference between the liquid and solid phases along with the change in Gibbs free energy of solidification. The dynamic pressure-volume–compression behaviors of ice phases (VI and VII) show a lattice stiffening in both phases, especially during the compression loading. It is also found that the compression rate greatly affects the solid-solid phase transition between ice VI and VII but does not affect the liquid-solid transition between water and ice VI as much. Lastly, a third phase transition was found to occur after metastable ice VII transforms into high-density amorphous (HDA) ice, which could be a disordered hydrogen-bonded network configuration of ice VII forming out of HDA ice facilitated by the decoupling of the oxygen movement and reorientation of the H 2⁡ O molecule. These results demonstrate the complexity of a seemingly simple molecule H 2⁡ O, how it can readily change its static properties with the modification of (de)compression rate, and highlight the need to use multiple TR structural and spectroscopic probes at higher time resolutions to realize the most comprehensive understanding.

Chemical bonding↗
DOE OSTI2022

In vitro investigation of protein assembly by combined microscopy and infrared spectroscopy at the nanometer scale

The nanoscale structure and dynamics of proteins on surfaces has been extensively studied using various imaging techniques, such as transmission electron microscopy and atomic force microscopy (AFM) in liquid environments. These powerful imaging techniques, however, can potentially damage or perturb delicate biological material and do not provide chemical information, which prevents a fundamental understanding of the dynamic processes underlying their evolution under physiological conditions. Here, we use a platform developed in our laboratory that enables acquisition of infrared (IR) spectroscopy and AFM images of biological material in physiological liquids with nanometer resolution in a cell closed by atomically thin graphene membranes transparent to IR photons. In this work, we studied the self-assembly process of S-layer proteins at the graphene-aqueous solution interface. The graphene acts also as the membrane separating the solution containing the proteins and Ca 2+ ions from the AFM tip, thus eliminating sample damage and contamination effects. The formation of S-layer protein lattices and their structural evolution was monitored by AFM and by recording the amide I and II IR absorption bands, which reveal the noncovalent interaction between proteins and their response to the environment, including ionic strength and solvation. Our measurement platform opens unique opportunities to study biological material and soft materials in general.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Cloud Process Coupling and Time Integration in the E3SM Atmosphere Model

Abstract In this study, we find significant sensitivity to the choice of time step for the Energy Exascale Earth System Model's atmospheric component, leading to large decreases in the magnitude of cloud forcing when the time step is reduced to 10 s. Reducing the time step size for the microphysics increases precipitation, leading to a drying of the atmosphere and an increase in surface evaporation. This effect is amplified when the microphysics is substepped together with other cloud physics processes. Coupling the model's dynamics and physics more frequently reduces cloud fraction at lower altitudes, while producing more cloud liquid at higher altitudes. Reducing the deep convection time step also reduces low cloud mass and cloud fraction. Together, these results suggest that cloud physics in a global circulation model can depend strongly on time step and, in particular, on the frequency with which cloud‐related processes are coupled with each other and with the model dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Exploring the Limits of Biological Complexity Amenable to Studies by Incoherent Neutron Spectroscopy

The wavelengths of neutrons available at neutron scattering facilities are comparable with intra- and inter-molecular distances, while their energies are comparable with molecular vibrational energies, making such neutrons highly suitable for studies of molecular-level dynamics. The unmistakable trend in neutron spectroscopy has been towards measurements of systems of greater complexity. Several decades of studies of dynamics using neutron scattering have witnessed a progression from measurements of solids to liquids to protein complexes and biomembranes, which may exhibit properties characteristic of both solids and liquids. Over the last two decades, the frontier of complexity amenable to neutron spectroscopy studies has reached the level of cells. Considering this a baseline for neutron spectroscopy of systems of the utmost biological complexity, we briefly review what has been learned to date from neutron scattering studies at the cellular level and then discuss in more detail the recent strides into neutron spectroscopy of tissues and whole multicellular organisms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Gas-Liquid Flow Modeling for Renewable Fuels Production

Aerobic/anaerobic and gas fermentation pathways have emerged as promising new technologies for the generation of renewable fuels/chemicals from biomass derived sugars, and mixtures of greenhouse/energy rich gas streams (CO2/CH4/H2/CO) via microbial action. Example pathways include sugars-to-ethanol conversion, biomethanation (CO2/H2 to CH4), biogas upgrading, CO fermentation and wet-waste conversion. Gas and liquid phase transport, mass-transfer, and mixing physics at large length scales can significantly affect microbial conversion rates, particularly when the microbial reaction requires a narrow set of conditions. These phenomena are difficult to study in small-scale bench-top reactors that are typically well-mixed. Predictive computational fluid dynamics (CFD) based simulations can therefore aid in the scale-up, design and optimization of these reactors. This work presents multiphase Euler-Euler CFD simulations of at-scale (~500 m3) bioreactors. Our mathematical model treats the gas and liquid as interpenetrating phases. This approach reduces the computational complexity of tracking individual gas bubbles that are several orders of magnitude smaller than reactor dimensions. We solve the Reynolds averaged Navier-Stokes (RANS) multiphase equations that account for phase and chemical species transport, interphase mass and momentum transfer and uses a phenomenological model for gas uptake by microbes. We use a customized solver derived from open-source CFD toolbox, OpenFOAM [1], to perform these simulations, which has been validated against small-scale reactors in our previous work [2]. There is currently a knowledge-gap regarding bubble-size distributions when using gas mixtures with vastly different properties, which can have a significant impact overall mass-transfer. For example, hydrogen bubbles are more buoyant compared to other relatively heavier gases (CO2/CH4/CO), resulting in a large distribution of residence times and bubble sizes. This work therefore develops a deeper understanding of bubble dynamics and interphase mass transfer in such heterogenous gas mixtures through well-resolved computational models. We use a population balance model (PBM) for bubble-size-distribution modeling that is validated against small-scale experiments in our solver with an uncertainty quantification study for bubble coalescence and break-up model parameters. Results pertaining to multiple simulations of gas-fermentation reactors are presented where gas mixtures with varying compositions of CO2/CH4/CO/H2 are imposed at the sparger boundaries. The spatio-temporal variations in bubble-size distribution and mass transfer coefficient are analyzed for varying superficial velocities and gas-compositions for varying sizes of bubble-column and airlift reactors. This work will also examine the performance of different reactor designs, viz. bubble column reactor, airlift reactor with an internal draft tube, and a stirred-tank reactor with Rushton impellers. Reactor mass-transfer coefficient, gas hold-up, and dissolved gas distribution are critically analyzed among reactors, and sensitivity studies pertaining to gas flow rates and reactor geometry will be presented. [1] Weller, H., Tabor, G., Jasak, H. and Fureby, C., A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12, 6, 620--631, 1998. [2] Rahimi, M., Sitaraman, H., Humbird, D. and Stickel, J., Computational fluid dynamics study of full-scale aerobic bioreactors: Evaluation of gas-liquid mass transfer, oxygen uptake, and dynamic oxygen distribution, Chemical Engineering Research and Design, 139: 283-295.

BIOMASS FUELS↗
DOE OSTI2021

Dynamics of local magnetic moments induced by itinerant Weyl electrons

We derive the effective interactions between local magnetic moments which are mediated by Weyl electrons in magnetic topological semimetals. The resulting spin dynamics is governed by the induced Heisenberg, Kitaev, and Dzyaloshinskii-Moriya (DM) interactions with extended range and oscillatory dependence on the distance between the spins. These interactions are realized in multiple competing channels shaped by the multitude of Weyl nodes in the electron spectrum. Microscopic spins need to be spatially modulated with a channel-dependent wave vector in order to take advantage of the interactions. The DM vector is parallel to the displacement between the two interacting spins, and requires the presence of Weyl electron Fermi surfaces. Here, we also derive the Weyl-induced chiral three-spin interaction in the presence of an external magnetic field. This interaction has an extended range as well, and acts upon the spatially modulated spins in various channels. Its tendency is to produce a skyrmion lattice or a chiral spin liquid which exhibits topological Hall effect. Ultimately, the theory developed here addresses magnetic dynamics in relativistic metals even when chiral magnetism is microscopically precluded. We discuss insights into the ordered state of the magnetic Weyl semimetal NdAlSi.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Dataset for manuscript "Rotational Memory Function of SPC/E water"

Memory effect are essential for dynamics of condensed materials and are responsible for non-exponential relaxation of correlation functions of dynamic variables through the memory function entering the memory equation. Memory functions of dipole rotations for polar liquids have never been calculated. We present here calculations of memory functions and single-dipole rotations and of the overall system dipole moment for SPC/E water measured by dielectric spectroscopy. The memory functions for single-particle and collective dynamics turn out to be nearly identical. This result validates theories of dielectric spectroscopy in terms of single-particle time correlation function and the connection between the collective and single-particle relaxation times in terms of the Kirkwood factor. The dataset includes single particle and system dipole moments, including their time-dependence.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

High-speed characterization of two-phase flow and bubble dynamics in titanium felt porous media for hydrogen production

Bubble dynamics and two-phase flow phenomena are closely related to the performance of proton exchange membrane electrolyzer cells (PEMECs). This paper reports an in-situ study of the oxygen bubble behavior and associated multiphase evolutions in the anode side of PEMECs with titanium (Ti) felt liquid gas diffusion layers (LGDLs) by a high-speed visualization system. The micro oxygen bubble dynamics was captured and analyzed at different locations and virous operating conditions. The results show that the bubble detachment frequency and detachment diameter greatly increase with the operating current density. Additionally, they are significantly impacted by the local pore structure and morphology of Ti felt LGDLs. In the flow channels, there exist only several discrete micro bubbles at a low current density (0.04 A/cm 2 ) and a large flow velocity (133 mm/s). At a current density (0.2 A/cm 2 ) and a flow velocity (67 mm/s), a number of gas slugs are formed in the follow channels, in addition to discrete micro bubbles. At a high current density (1 A/cm 2 ) and a flow velocity (67 mm/s), more bubbles appear in the channel, and the flow field is dominated by slug or annular flows. These investigations can help to better understand the two-phase flow and bubble detachment mechanism, and provide a foundation for electrochemical reaction, multiphase flow studies and optimize the design of gas diffusion layers and flow fields for PEMECs in the future.

08 HYDROGEN↗
DOE OSTI2024

Elevated Temperature Effects on Protein Turnover Dynamics in Arabidopsis thaliana Seedlings Revealed by 15 N-Stable Isotope Labeling and ProteinTurnover Algorithm

Global warming poses a threat to plant survival, impacting growth and agricultural yield. Protein turnover, a critical regulatory mechanism balancing protein synthesis and degradation, is crucial for the cellular response to environmental changes. We investigated the effects of elevated temperature on proteome dynamics in Arabidopsis thaliana seedlings using 15 N-stable isotope labeling and ultra-performance liquid chromatography-high resolution mass spectrometry, coupled with the ProteinTurnover algorithm. Analyzing different cellular fractions from plants grown under 22 °C and 30 °C growth conditions, we found significant changes in the turnover rates of 571 proteins, with a median 1.4-fold increase, indicating accelerated protein dynamics under thermal stress. Notably, soluble root fraction proteins exhibited smaller turnover changes, suggesting tissue-specific adaptations. Significant turnover alterations occurred with redox signaling, stress response, protein folding, secondary metabolism, and photorespiration, indicating complex responses enhancing plant thermal resilience. Conversely, proteins involved in carbohydrate metabolism and mitochondrial ATP synthesis showed minimal changes, highlighting their stability. This analysis highlights the intricate balance between proteome stability and adaptability, advancing our understanding of plant responses to heat stress and supporting the development of improved thermotolerant crops.

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

Catalytic disproportionation on carbon superstructures enables long-life, high-loading Li–S batteries

Electrocatalysis has been widely explored as an effective strategy to accelerate polysulfide (PS) conversion and suppress the shuttle effect in lithium–sulfur (Li–S) batteries. However, the underlying mechanisms remain elusive, and electrocatalytic reactions are inactive during cell resting. In this work, we reveal and quantitatively analyze a previously unrecognized sulfur reduction route (SRR) driven by catalytic disproportionation at the carbon cathode surface—fundamentally distinct from conventional electrocatalysis. Unlike conventional stepwise pathways, this SRR enables high-order polysulfides (Sₓ²⁻, x = 5–8) to directly convert into S₈ and Li₂S₂, bypassing low-order intermediates. This sulfur-reduction shortcut is systematically elucidated through high-performance liquid chromatography, revealing the intrinsic catalytic contribution of carbon frameworks and the dynamic evolution of PS species. We demonstrate that carbon superstructures (CSS-0.5), assembled from nanosheet subunits with abundant N/O functionalities and interconnected charge-migration channels, synergistically promote this catalytic process. Benefiting from these features, CSS-0.5 delivers superior electrochemical performance under practical conditions, enabling high sulfur loading (6.0 mg cm⁻²) pouch cells with 80.5% capacity retention over 210 cycles. This study provides the first quantitative evidence of electrocatalytic disproportionation in Li–S batteries, offering mechanistic insights and design principles for advanced sulfur cathodes.

25 ENERGY STORAGE↗
DOE OSTI2021

A novel semi-structured kinetic model of methanotroph-photoautotroph cocultures for biogas conversion

Through metabolic coupling of methane oxidation and oxygenic photosynthesis, methanotroph-photoautotroph (M-P) cocultures offer a highly promising technology platform for biogas conversion. However, there has not been any quantitative modeling of the coculture growth kinetics. This is mainly due to the inherent difficulty associated with real time characterization of the M-P cocultures and the complex interactions such as the cross-feeding mechanism within the coculture. To address this challenge, we recently developed a novel experimental-computational (E-C) protocol to accurately characterize the M-P coculture in real-time, and validated its accuracy through cell counting. Enabled by the E-C protocol, this work presents the very first kinetic model for M-P cocultures. By explicitly modeling the exchange of in situ produced O 2 /CO 2 within the M-P coculture and coupling the individual biomass growth with mass transfer between the gas and liquid phases, the semi-structured kinetic model accurately predicts the growth dynamics of the M-P coculture under a wide range of growth conditions. The proposed model is validated by a series of wet-lab experiments using Methylomicrobium buryatense 5GB1 - Arthrospira platensis as the model coculture. Although it has been speculated that there may exist other emergent metabolic interactions within the M-P coculture, in addition to the exchange of in situ produced O 2 /CO 2 , there has not been any experimental validation prior to this study. By integrating designed experiments with the semi-structured kinetic model, this study is the first to confirm the existence of the additional emergent metabolic exchanges within the coculture. Furthermore, this study further quantifies the effect of these unknown metabolic interactions on the growth of both species in the model coculture, supporting further research to identify these exchanged metabolites for metabolic engineering.

42 ENGINEERING↗
DOE OSTI2022

Insights into Sorption and Molecular Transport of Aqueous Glucose into Zeolite Nanopores

Liquid-phase heterogeneous catalysis using zeolites is important for biomass conversion to fuels and chemicals. There is a substantial body of work on gas-phase sorption in zeolites with different topologies; however, studies investigating the diffusion of complex molecules in liquid medium into zeolitic nanopores are scarce. Here, we present a molecular dynamics study to understand the sorption and diffusion of aqueous β-d-glucose into β-zeolite silicate at T = 395 K and P = 1 bar. Through 2-μs-long molecular dynamics trajectories, we reveal the role of the solvent, the kinetics of the pore filling, and the effect of the water model on these properties. We find that the glucose and water loading is a function of the initial glucose concentration. Although the glucose concentration increases monotonically with the initial glucose concentration, the water loading exhibits a nonmonotonic behavior. At the highest initial concentration (~20 wt %), we find that the equilibrium loading of glucose is approximately five molecules per unit cell and displays a weak dependence on the water model. Glucose molecules follow a single-file diffusion in the nanopores due to confinement. Here, the dynamics of glucose and water molecules slows significantly at the interface. The average residence time for glucose molecules is an order of magnitude larger than that in the bulk solution, while it is about twice as large for the water molecules. Our simulations reveal critical molecular details of the glucose molecule’s local environment inside the zeolite pore relevant to catalytic conversion of biomass to valuable chemicals.

09 BIOMASS FUELS↗
DOE OSTI2023

Relevance of structural defects to the mechanism of mechanical deformation in metallic glasses

Abstract It is known that deformation in disordered materials such as metallic glasses and supercooled liquids occurs via the cooperative rearrangement of atoms or constituent particles at dynamical heterogeneities, commonly regarded as point-like defects. We show via molecular-dynamics simulations that there is no apparent relationship between atomic rearrangements and the local atomic environment as measured by the atomic-level stresses, kinetic and potential energies, and the per-atom Voronoi volume. In addition, there is only a weak correlation between atomic rearrangements and the largest and smallest eigenvalues of the dynamical matrix. Our results confirm the transient nature of dynamical heterogeneities and suggest that the notion of defects may be less relevant than that of a propensity for rearrangement.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Operando probing dynamic migration of copper carbonyl during electrocatalytic CO2 reduction

Single crystals and shape-controlled nanocrystals are well known to exhibit facet-dependent catalytic properties. However, few studies have investigated how those nanocrystals evolve and (de)activate during reactions, calling for the development of nanoscale time-resolved operando methods. In this context, we have designed Cu nanocubes as a model system to elucidate the underlying driving force of dynamic nanocatalyst reconstruction during the CO2 reduction reaction (CO2RR). Operando electrochemical liquid-cell scanning transmission electron microscopy (EC-STEM) and synchrotron-based X-ray spectroscopy reveal the size- and potential-dependent complete transformation from (100)-oriented Cu@Cu2O nanocubes to polycrystalline metallic Cu nanograins under CO2RR conditions. In addition, machine learning-assisted operando four-dimensional STEM reveals that large Cu nanograins derived from nanocubes form mainly crystalline domains, while their smaller counterparts are more amorphous due to faster evolution kinetics. In situ Raman spectroscopy and density functional theory calculations suggest that CO drives the ejection of single Cu atoms, resulting in few-nanometre Cu clusters and the surface migration of highly mobile copper carbonyl (Cu–CO) species. Combined, these multimodal operando methods and theoretical approaches pave the way for understanding the complex structural evolution of energy-related nanocatalysts under electrochemical conditions.

Yang, Yao↗
DOE OSTI2025

Influence of Pore Water on the Fracture of Silica Sand at Particle Scale

The fracture of particles has a significant influence on the engineering behavior of granular materials. There are no reported conclusive answers in the literature to explain the influence of pore water on the fracture of sand, gravel, railroad ballast rock, aggregates, rockfill, and so on. Here, in this paper, two novel techniques were adopted to investigate the influence of pore water on the fracture of natural silica sand at the particle scale. The three-dimensional (3D) synchrotron microcomputed tomography (SMT) imaging technique was used to acquire and analyze multiple 3D images of specimens composed of wetted silica sand that were subjected to confined one-dimensional (1D) compression loading up to the fracture stage. A few representative particles were identified and digitally separated from the liquid water and gas phases. The 3D SMT images offer clear experimental evidence of the phase transition of water from liquid to gas within the opening cracks. In addition, a computational fluid dynamics numerical simulation framework for water flow into an opening crack was adopted and the results show that the fast opening of the crack induced cavitation, water phase transition, and water hammer effects. The results of the numerical simulations agree with the confined 1D compression experiments where the onset of the bubbles (gas phase) within the crack is mainly generated during the secondary cavitation stage, and the generated bubble near the crack mouth opening resulted from the primary cavitation, which was directly caused by the crack expansion. The results reported in this paper present a new cavitation phenomenon that occurs during particle fracture, which offers a physics-based explanation of why water-wetted sand fractures at smaller stresses than dry and can pave the way for more in-depth future studies on the fracture of water-wetted sand.

cavitation↗
DOE OSTI2020

Vortex guidance and transport in channeled pinning arrays

Here, we numerically examine vortices interacting with pinning arrays where a portion of the pinning sites have been removed in order to create coexisting regions of strong and weak pinning. The region without pinning sites acts as an easy-flow channel. For driving in different directions with respect to the channel, we observe distinct types of vortex flow. When the drive is parallel to the channel, the flow first occurs in the pin free region followed by a secondary depinning transition in the pinned region. At high vortex densities there is also an intermediate plastic flow phase due to the coupling between the weak and strong pinning regions. For driving applied perpendicular to the channel, we observe a jammed phase in which vortices accumulate on the boundary of the pinned region due to the vortex-vortex repulsion, while at higher drives the vortices begin to flow through the pinning array. For driving at an angle to the channel, depending on the filling we observe a drive-induced reentrant pinning effect as well as negative differential mobility which occurs when vortices move from the unpinned to the pinned portion of the sample.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A compartmentalized model of multiphase chemical kinetics

There are significant challenges in predicting multiphase chemical kinetics due to the complex coupling of reaction and mass transport across a phase boundary (i.e., interface). Here, we describe a framework for predicting multiphase kinetics that embeds the elementary kinetic steps of reaction, solvation, and diffusion into a coarse grain spatial description of two phases. The model is constructed to bridge the short-timescale interfacial dynamics observed in molecular simulations with the longer timescales observed in kinetic experiments. A simple set of governing differential equations is derived, which, when solved numerically or analytically, yield accurate predictions of multiphase kinetics in microdroplets. Although the equations are formulated for gas-liquid reactions, the underlying conceptual framework is general and can be applied to transformations in other two-phase systems (solid-liquid, liquid-liquid, etc.).

Chemical kinetics and dynamics↗
DOE OSTI2026

Filament extension atomizers: A study of integration potential into spray dryers for high-viscosity and non-Newtonian fluids

Filament extension atomizers (FEAs) are an emerging class of spray nozzles designed to atomize high-viscosity and non-Newtonian fluids that are challenging for conventional pressure or two-fluid nozzles. In this study, we investigate the influence of roller geometry, surface velocity, roller material, and fluid rheology on the atomization performance of FEA systems using concentrated whey protein suspensions (50–70 wt. %). Extensional and shear rheology experiments, along with high-speed imaging and particle image velocimetry, reveal that filament breakup dynamics are governed by competition between inertial, capillary, and viscoelastic stresses. High roller rotational velocity leads to narrower spray cones, contradicting rheology experiments and suggests a significant inertial contribution to filament breakup. Smaller rollers operating at the same rotational velocity led to broader spray cones consistent with expectations. An FEA nozzle was integrated into a conventional dryer producing particles in the 100 μm range. Results suggest that FEA technology enables atomization of highly viscous fluids at industrially relevant spray cone angles similar to those generated by pressure nozzles, offering a pathway to improve energy efficiency in spray drying by enabling higher solids loading feedstocks. Furthermore, these insights provide critical guidance for optimizing FEA nozzle designs and process parameters across a range of applications.

Energy efficiency↗