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 847 records · Page 47

Aligning the Induced Anisotropy of Isotropic Nanoparticles with Liquid Crystals

Thermotropic liquid crystals (LCs) present opportunities for synergistic interactions with ligand-functionalized nanoparticles (NPs). Understanding the dynamics and structures of ligand-coated NPs in an anisotropic LC environment allows for better design and control of versatile LC-NP hybrid materials. Here, simulations and experiments yield direct evidence that cyanobiphenyl LCs induce anisotropy in the biphenylalkyl ligand shells of spherical NPs. The ellipsoidal NP aligns with the LC director. Magnetic fields can dictate the directional ordering of mesogens and consequently allow control of the orientation of the nonmagnetic NPs. Further, controlling the alignment and deformation advances the design and engineering of these hybrid nanomaterials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Identifying High Ionic Conductivity Compositions of Ionic Liquid Electrolytes Using Features of the Solvation Environment

Binary mixtures of ionic liquids with molecular solvents are gaining interest in electrochemical applications due to the improvement in their performance over neat ionic liquids. Dilution with suitable molecular solvents can reduce the viscosity and facilitate faster diffusion of ions, thereby yielding substantially higher ionic conductivity than that for a pure ionic liquid. Although viscosity and diffusion coefficients typically behave as monotonic functions of concentration, ionic conductivity often passes through a peak value at an optimum molar ratio of the molecular solvent to the ionic liquid. The ionic conductivity maximum is generally explained in terms of a balance between the ease of charge transport and the concentration of the charge carriers. In this work, fluctuation in the local environment surrounding an ion is invoked as a plausible explanation for the ionic conductivity mechanism with a binary mixture of 1-ethyl-3-methylimidazolium tetrafluoroborate and ethylene glycol as an example. The magnitude of the dynamism in the local environment is captured by measuring the spatial and temporal features of the solvation environment. Standard deviation in the number of ions in the solvation environment serves as a spatial feature, while the cage correlation lifetimes for oppositely charged ions within the first solvation shell serve as a temporal feature. Large standard deviations in the cluster ion population and short cage correlation lifetimes are indicators of highly dynamic ionic environment at the molecular level and consequently yield high ionic conductivity. Such compositions were found to be in good agreement with the optimum ionic liquid mole fractions obtained through experimental measurement. Short cage correlation lifetimes enable the identification of optimum mixture compositions using simulation trajectories significantly shorter than those required to implement the Nernst–Einstein or Einstein formalisms for calculating ionic conductivity. We validated the applicability of this approach across force fields and in six ionic liquid-molecular solvent electrolytes formed with combination of cations, anions, and solvents. We offer a computationally efficient approach of screening ionic liquid-molecular solvent binary mixture electrolytes to identify molar ratios that yield high ionic conductivity.

25 ENERGY STORAGE↗

Study of dielectric breakdown in liquid xenon with XeBrA: The xenon breakdown apparatus

Maintaining the electric fields necessary for the current generation of noble liquid time projection chambers (TPCs), with drift lengths exceeding 1 m, requires a large negative voltage applied to their cathode. Delivering such high voltage is associated with an elevated risk of electrostatic discharge and electroluminescence, which would be detrimental to the performance of the experiment. The Xenon Breakdown Apparatus (XeBrA) is a 5-l, high voltage test chamber built to investigate the contributing factors to electrical breakdown in noble liquids. Here, in this work, we present the main findings after conducting scans over stressed electrode areas, surface finish, pressure, and high voltage ramp speed in the medium of liquid xenon. Area scaling and surface finish were observed to be the dominant factors affecting breakdown, whereas no significant changes were observed with varying pressure or ramp speed. A general rise in both the anode current and photon rate was observed in the last 30 s, leading up to a breakdown, with a marked increase in the last couple of seconds. In addition, the position of breakdowns was reconstructed with a system of high-speed cameras and a moderate correlation with the Fowler–Nordheim field emission model was found. Tentative evidence for bubble nucleation being the originating mechanism of breakdown in the liquid was also observed. We deem the results presented in this work to be of particular interest for the design of future, large TPCs, and practical recommendations are provided.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

X-ray and molecular dynamics study of the temperature-dependent structure of molten NaF-Zr⁢F 4

The local atomic structure of NaF-Zr⁢F 4 (53–47 mol%) molten system and its evolution with temperature are examined with x-ray scattering measurements which are then used to validate the quality of ab initio and neural network-based molecular dynamics (NNMD) calculations in the temperature range 515–700°⁢C. The machine-learning enhanced NNMD calculations offer improved efficiency while maintaining accuracy at higher distances compared to ab initio calculations. Looking at the evolution of the pair distribution function with increasing temperature, a fundamental change in the liquid structure within the selected temperature range, accompanied by a slight decrease in overall correlation is revealed. NNMD calculations indicate the coexistence of three different fluorozirconate complexes: [Zr⁢F 6 ] 2– , [Zr⁢F 7 ] 3– , and [Zr⁢F 8 ] 4– , with a shift in the dominant coordination state from the 7-coordinated Zr cation toward a 6-coordinated cation with increasing temperature. The study also highlights the metastability of different local coordination structures, with frequent interconversions between the 6- and 7-coordinate states. Analysis of the Zr-F-Zr angular distribution function reveals the presence of both “edge-sharing” and “corner-sharing” fluorozirconate complexes with specific bond angles and distances in accord with previous studies, while the next-nearest-neighbor cation-cation correlations demonstrate a clear preference for unlike cations as nearest-neighbor pairs, emphasizing nonrandom arrangement. Finally, these findings contribute to a comprehensive understanding of the complex local structure of the molten salt, providing insights into temperature-dependent preferences and correlations within the molten system.

36 MATERIALS SCIENCE↗

Free-spray characteristics and spray-wall interactions of methanol on a gasoline direct injector under flash-boiling and non-flash-boiling conditions

Methanol is considered a promising alternative fuel for internal combustion engines (ICEs) due to its high-octane number, fast laminar flame speed, and elevated latent heat of vaporization, all of which support higher compression ratios and improved thermal efficiency. However, its substantial latent heat of vaporization also poses cold-start challenges, such as misfire and fuel film deposition. This study aims to investigate methanol spray morphology and spray-wall interaction using the Spray M injector from the Engine Combustion Network within a constant-pressure flow vessel. A recently developed unified numerical framework capable of modeling both flash and non-flash boiling sprays is validated against experimental liquid volume fraction data acquired via 3-D computed tomography. Here, the results reveal that flash boiling significantly alters the spray morphology, leading to smaller droplets and spray collapse due to enhanced air-entrainment-induced turbulence. Quantitative agreement between experiments and simulations confirms this behavior. Coupled 0-D equilibrium and 3-D computational fluid dynamics analyses show that flash boiling accelerates evaporation and reduces fuel residence time, while non-flash conditions maintain a persistent liquid core more susceptible to wall wetting. Wall temperature diagnostics reveal that spray collapse alters heat transfer patterns by shifting cooling effects. Mixture fraction analysis indicates that evaporation is primarily governed by shear-layer turbulence, though deviations from adiabatic equilibrium mixing emerge under low-turbulence conditions. Finally, increasing fuel, ambient, and wall temperatures reduces wall wetting and film thickness, mitigating cold-start risks. These findings enhance the understanding of methanol sprays’ behavior and support its adoption as a viable, alternative fuel for ICEs.

Engine Combustion Network↗

Flexible cryogenic air separation unit—An application for low-carbon fossil-fuel plants

In this study, the rapid integration of intermittent renewable sources into the electricity grid is driving the need for a flexible cryogenic air separation unit (ASU) coupled with a low-carbon fossil-fuel plant. However, the state-of-the-art ASU process is highly integrated and nonlinear, which can significantly restrict its ramping rate. In this work, we study the fundamental dynamics of a state-of-the-art double-column ASU, focusing on the dynamic characteristics of the highly integrated and nonlinear heat and mass transfers, and we propose basic control methods for achieving high process ramping rates. We found that the vapor–liquid countercurrent flows in the low-pressure column are critical to the cryogenic rectification of air, which governs the ramp rate of the ASU. This countercurrent flow structure is created through a complex heat integration process in the ASU. Here, this process is simplified as a countercurrent heat transfer to reduce the complexity of studying the ASU dynamics. To preserve this flow structure, the heat integration is maintained so that its heat duty follows the ASU load, using several basic controllers on the critical stream flowrates. A flow-driven dynamic ASU model is built in Aspen Plus Dynamics to capture the basic ramping dynamics. This model revealed a fundamental mismatch in the dynamics in the ASU column that causes a significant loss of O 2 product purity when ramping down the ASU and slightly increases the purity when ramping up. Based on these findings, we propose a basic control method for rapidly ramping down/up the ASU while maintaining O 2 product purity. Simulation results show the ASU basic dynamic process successfully ramps at a rate up to 10 %/min (40–100 % load) while maintaining the O 2 product purity at 95.2–95.6 mol%.

42 ENGINEERING↗

Studying Interfacial Dark Reactions of Glyoxal and Hydrogen Peroxide Using Vacuum Ultraviolet Single Photon Ionization Mass Spectrometry

Aqueous secondary organic aerosol (aqSOA) formation from volatile and semivolatile organic compounds at the air–liquid interface is considered as an important source of fine particles in the atmosphere. However, due to the lack of in situ detecting techniques, the detailed interfacial reaction mechanism and dynamics still remain uncertain. In this study, synchrotron-based vacuum ultraviolet single-photon ionization mass spectrometry (VUV SPI-MS) was coupled with the System for Analysis at the Liquid Vacuum Interface (SALVI) to investigate glyoxal dark oxidation products at the aqueous surface. Mass spectral analysis and determination of appearance energies (AEs) suggest that the main products of glyoxal dark interfacial aging are carboxylic acid related oligomers. Furthermore, the VUV SPI-MS results were compared and validated against those of in situ liquid time-of-flight secondary ion mass spectrometry (ToF-SIMS). The reaction mechanisms of the dark glyoxal interfacial oxidation, obtained using two different approaches, indicate that differences in ionization and instrument operation principles could contribute to their abilities to detect different oligomers. Therefore, the mechanistic differences revealed between the VUV SPI-MS and ToF-SIMS indicate that more in situ and real-time techniques are needed to investigate the contribution of the air–liquid interfacial reactions leading to aqSOA formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Topological orders of monopoles and hedgehogs: From electronic and magnetic spin-orbit coupling to quarks

Topological states of matter are, generally, quantum liquids of conserved topological defects. We establish this by constructing and analyzing topological field theories which introduce gauge fields to describe the dynamics of singularities in the original field configurations. In this work, homotopy groups are utilized to identify topologically protected singularities, and the conservation of their protected number is captured by a topological action term that unambiguously obtains from the given set of symmetries. Stable phases of these theories include quantum liquids with emergent massless Abelian and non-Abelian gauge fields, as well as topological orders with long-range quantum entanglement, fractional excitations, boundary modes, and unconventional responses to external perturbations. This paper focuses on the derivation of topological field theories and basic phenomenological characterization of topological orders associated with homotopy groups π n (S n ), n ≥ 1 . These homotopies govern monopole and hedgehog topological defects in d = n + 1 dimensions, and enable the generalization of both weakly interacting and fractional quantum Hall liquids of vortices to d > 2 . Hedgehogs have not been in the spotlight so far, but they are particularly important defects of magnetic moments because they can be stimulated in realistic systems with spin-orbit coupling, such as chiral magnets and d = 3 topological materials. We predict topological orders in systems with U(1) × Spin(d) symmetry in which fractional electric charge attaches to hedgehogs. Monopoles, the analogous defects of charge or generic U(1) currents, may bind to hedgehogs via Zeeman effect, or effectively emerge in purely magnetic systems. The latter can lead to spin liquids with different topological orders than that of the resonant valence bond spin liquid. Charge fractionalization of quarks in atomic nuclei is also seen as possibly arising from the charge-hedgehog attachment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Potentials of mean force fail to describe chemical bond-breaking in solution

Many liquid phase studies assume that the potential energy surfaces of reacting molecules are the same as in the gas phase, neglecting complex solvent dynamics that can completely alter the nature of chemical reactivity. Even studies that include solvent effects typically only consider them in an average, equilibrium way as part of a potential of mean force (PMF). In this work, we use mixed quantum/classical simulations to compare how equilibrium and non-equilibrium solvent motions affect the photodissociation of a simple diatomic molecule, NaK + , in liquid tetrahydrofuran. A PMF analysis shows that as the excited-state molecule dissociates with the solvent at equilibrium, the bonding electron remains associated with K + at short bond distances but eventually localizes on Na + at the end of dissociation. When we examine non-equilibrium dynamical photodissociation trajectories, however, we find that they fall into three distinct categories: about a quarter of them have the bonding electron mainly associated with Na + , another quarter stay mainly associated with K + , and about half have the bonding electron shared roughly equally between the two ions. The results show that equilibrium PMFs cannot accurately describe the dynamics of bond-breaking chemical reactions in solution because there is insufficient time for the solvent to reach equilibrium on the time scale over which bond dissociation occurs. Furthermore, our analysis shows that the solvent coupling between the electronic energy surfaces is similar at and away from equilibrium, suggesting that other factors, such as solute velocity-driven solvent memory effects, play a more important role in explaining the failure of the equilibrium PMF to predict the non-equilibrium dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

From bulk to surface: Structure and dynamics of amorphous alumina from deep potential molecular dynamics

Understanding the atomic-scale structure and dynamics of amorphous oxide surfaces is essential for interpreting their chemical reactivity, mechanical stability, and interfacial behavior, yet direct experimental characterization remains challenging. We employ Deep Potential (DP) molecular dynamics to generate large-scale, ab initio -quality models of amorphous Al 2 O 3 bulk glasses and melt-quenched free surfaces, enabling a quantitative analysis of both structure and relaxation dynamics with statistical confidence inaccessible to direct ab initio simulation. The trained DP model reproduces experimental liquid and glass structure, captures the cooling-rate dependence of the bulk glass transition, and corrects systematic biases in the polyhedral populations predicted by widely used classical force fields. At the free surface, mass density recovers to bulk values over ~10 Å, while local coordination requires a slightly wider subsurface region to fully converge. The outermost layer is oxygen-enriched, exhibits altered polyhedral connectivity with contracted Al–O bonds, and hosts a broad population of under-coordinated motifs (notably AlO 3 and OAl 2 ) whose abundances are governed by glass stability. These under-coordinated surface motifs exhibit distinct vibrational signatures and occur as locally paired Lewis acid and Brønsted base sites consistent with bond-valence compensation, yet remain spatially dispersed rather than aggregating into extended clusters. Despite this pronounced structural heterogeneity, surface relaxation and the glass-transition temperature remain comparable to their bulk counterparts, suggesting that the disordered surface is kinetically stable once formed. Together, these results establish a molecular-level picture of amorphous alumina surfaces and demonstrate the capability of machine-learned potentials to resolve structure–property relationships in disordered oxide interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bayesian force fields from active learning for simulation of inter-dimensional transformation of stanene

Abstract We present a way to dramatically accelerate Gaussian process models for interatomic force fields based on many-body kernels by mapping both forces and uncertainties onto functions of low-dimensional features. This allows for automated active learning of models combining near-quantum accuracy, built-in uncertainty, and constant cost of evaluation that is comparable to classical analytical models, capable of simulating millions of atoms. Using this approach, we perform large-scale molecular dynamics simulations of the stability of the stanene monolayer. We discover an unusual phase transformation mechanism of 2D stanene, where ripples lead to nucleation of bilayer defects, densification into a disordered multilayer structure, followed by formation of bulk liquid at high temperature or nucleation and growth of the 3D bcc crystal at low temperature. The presented method opens possibilities for rapid development of fast accurate uncertainty-aware models for simulating long-time large-scale dynamics of complex materials.

Chemistry↗

Global Phase Diagram of a Spin-Orbital Kondo Impurity Model and the Suppression of Fermi-Liquid Scale

Many correlated metallic materials are described by Landau Fermi-liquid theory at low energies, but for Hund metals the Fermi-liquid coherence scale $T_{FL}$ is found to be surprisingly small. In this Letter, we study the simplest impurity model relevant for Hund metals, the three-channel spin-orbital Kondo model, using the numerical renormalization group (NRG) method and compute its global phase diagram. In this framework, $T_{FL}$ becomes arbitrarily small close to two new quantum critical points that we identify by tuning the spin or spin-orbital Kondo couplings into the ferromagnetic regimes. We find quantum phase transitions to a singular Fermi-liquid or a novel non-Fermi-liquid phase. The new non-Fermi-liquid phase shows frustrated behavior involving alternating overscreenings in spin and orbital sectors, with universal power laws in the spin (ω -1/5 ), orbital (ω 1/ 5 ) and spin-orbital (ω 1 ) dynamical susceptibilities. These power laws, and the NRG eigenlevel spectra, can be fully understood using conformal field theory arguments, which also clarify the nature of the non-Fermi-liquid phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Numerical investigation of liquid wall ablation in inertial fusion energy chambers

This paper presents a novel approach for modeling liquid wall ablation in liquid wall-protected inertial fusion energy (IFE) chambers. These systems are promising candidates for the implementation of fusion technology, yet significant gaps remain in understanding the underlying physical processes and their implications for design. Following target ignition, a portion of the fusion energy is released as x-rays, which deposit their energy into an array of liquid jets, leading to partial vaporization. Accurately modeling this heat deposition and vaporization process remains challenging due to the complex geometries typical of (pre-conceptual) IFE chamber designs. Furthermore, the subsequent expansion of vaporized material into the chamber’s vacuum environment poses difficulties for conventional CFD methods based on continuum assumptions, which can lead to significant inaccuracies. To address some aspects of these challenges, this work introduces a ray-tracing-based methodology to map the spatial distribution of ablated material in liquid wall-protected systems. In addition, a vacuum-tracking scheme is developed to extend the applicability of an OpenFOAM-based solver to gas dynamics in rarefied environments. The proposed approach has been verified through numerical benchmarks and applied to a practical case involving the HYLIFE-II (High Yield Lithium Injection Fusion Energy) chamber. The methodology advances the modeling capabilities for liquid wall-protected IFE systems and provides valuable tools to support their design and optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Polymers in Deep Eutectic Solvents

The project investigated the behavior of polymers in ionic liquids such as deep eutectic solvents using an array of techniques. These included the development of atomistic force fields, coarse graining these force fields to the united atom level, large scale molecular dynamics (MD) simulations using new thermostat algorithms, and machine learning (ML) methods for the phase behavior. The project demonstrated the feasibility and accuracy of first principles force fields for ionic liquids, deep-eutectic solvents, and urea-water mixtures. A novel hierarchical coarse graining method was then used to develop accurate and efficient united-atom models, using which microsecond simulations were performed for polymers in ionic liquids. These simulations were in quantitative agreement with experiment, thus resolving previous controversies. Methods were also developed to obtain the potential of mean force between complex ions in solution. Finally, supervised ML methods were developed for the phase behavior of polymers in ionic liquids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Partition between supercooled liquid droplets and ice crystals in mixed-phase clouds based on airborne in situ observations

Abstract. The onset of ice nucleation in mixed-phase clouds determines the lifetime and microphysical properties of ice clouds. In this work, we develop a novel method that differentiates between various phases of mixed-phase clouds, such as clouds dominated by pure liquid or pure ice segments, compared with those having ice crystals surrounded by supercooled liquid water droplets or vice versa. Using this method, we examine the relationship between the macrophysical and microphysical properties of Southern Ocean mixed-phase clouds at −40 to 0 °C (e.g. stratiform and cumuliform clouds) based on the in situ aircraft-based observations during the US National Science Foundation Southern Ocean Clouds, Radiation, Aerosol Transport Experimental Study (SOCRATES) flight campaign. The results show that the exchange between supercooled liquid water and ice crystals from a macrophysical perspective, represented by the increasing spatial ratio of regions containing ice crystals relative to the total in-cloud region (defined as ice spatial ratio), is positively correlated with the phase exchange from a microphysical perspective, represented by the increasing ice water content (IWC), decreasing liquid water content (LWC), increasing ice mass fraction, and increasing ice particle number fraction (IPNF). The mass exchange between liquid and ice becomes more significant during phase 3 when pure ice cloud regions (ICRs) start to appear. Occurrence frequencies of cloud thermodynamic phases show a significant phase change from liquid to ice at a similar temperature (i.e. −17.5 °C) among three types of definitions of mixed-phase clouds based on ice spatial ratio, ice mass fraction, or IPNF. Aerosol indirect effects are quantified for different phases using number concentrations of aerosols greater than 100 or 500 nm (N>100 and N>500, respectively). N>500 shows stronger positive correlations with ice spatial ratios compared with N>100. This result indicates that larger aerosols potentially contain ice-nucleating particles (INPs), which facilitate the formation of ice crystals in mixed-phase clouds. The impact of N>500 is also more significant in phase 2 when ice crystals just start to appear in the mixed phase compared with phase 3 when pure ICRs have formed, possibly due to the competing aerosol indirect effects on primary and secondary ice production in phase 3. The thermodynamic and dynamic conditions are quantified for each phase. The results show stronger in-cloud turbulence and higher updraughts in phases 2 and 3 when liquid and ice coexist compared with pure liquid or ice (phases 1 and 4, respectively). The highest updraughts and turbulence are seen in phase 3 when supercooled liquid droplets are surrounded by ice crystals. These results indicate both updraughts and turbulence support the maintenance of supercooled liquid water amongst ice crystals. Overall, these results illustrate the varying effects of aerosols, thermodynamics, and dynamics through various stages of mixed-phase cloud evolution based on this new method that categorizes cloud phases.

54 ENVIRONMENTAL SCIENCES↗

Dynamics of fractionalized mean-field theories: Consequences for Kitaev materials

There have been substantial recent efforts, both experimentally and theoretically, to find a material realization of the Kitaev spin liquid—the ground state of the exactly solvable Kitaev model on the honeycomb lattice. Candidate materials are now plentiful, but the presence of non-Kitaev terms makes comparison between theory and experiment challenging. Here, we rederive time-dependent Majorana mean-field theory and extend it to include quantum phase information, allowing the direct computation of the experimentally relevant dynamical spin-spin correlator, which reproduces exact results for the unperturbed model. In contrast to previous work, we find that small perturbations do not substantially alter the exact result, implying that α-RuCl 3 is perhaps farther from the Kitaev phase than originally thought. Our approach generalizes to any correlator and to any model where Majorana mean-field theory is a valid starting point.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Fingerprints of hydrogen bonding in the terahertz dynamics of ethanol and water: An inelastic x-ray scattering study

Here, we used inelastic x-ray scattering methods to measure the terahertz spectrum of density fluctuations of ethanol in both liquid and solid phases. The results of a Bayesian inference-based lineshape analysis with a multiple excitation model and the comparison with a previous similar analysis on water indicate that the different structures induced by hydrogen bonds in ethanol and water have a profound influence on the respective dynamic responses, the latter being characterized by longer living and better resolved high-frequency acoustic excitations. In addition, we compare these findings with those obtained with an alternative approach based on the exponential expansion theory and ensuring sum rules fulfillment, demonstrating that the model’s choice directly impacts the number of spectral modes detected.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗