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 451 records · Page 25

The Strong Impact of Weak Horizontal Convergence on Continental Shallow Convection

Idealized large-eddy simulation (LES) is a basic tool for studying three-dimensional turbulence in the planetary boundary layer. LES is capable of providing benchmark solutions for parameterization development efforts. However, real small-scale atmospheric flows develop in heterogeneous and transient environments with locally varying vertical motions inherent to open multiscale interactive dynamical systems. These variations are often too subtle to detect them by state-of-the-art remote and in situ measurements, and are typically excluded from idealized simulations. The present study addresses the impact of weak [i.e., O(10-6) s-1] short-lived low-level large-scale convergence/divergence perturbations on continental shallow convection. The results show a strong response of shallow nonprecipitating convection to the applied weak large-scale dynamical forcing. Evolutions of CAPE, mean liquid water path, and cloud-top heights are significantly affected by the imposed convergence/divergence. In contrast, evolving cloud-base properties, such as the area coverage and mass flux, are only weakly affected. To contrast those impacts with microphysical sensitivity, the baseline simulations are perturbed assuming different observationally based cloud droplet number concentrations and thus different rainfall. For the tested range of microphysical perturbations, the imposed convergence/divergence provides significantly larger impact than changes in the cloud microphysics. Simulation results presented here provide a stringent test for convection parameterizations, especially important for large-scale models progressing toward resolving some nonhydrostatic effects.

Meteorology & Atmospheric Sciences↗

Atomic-level mechanisms of short-circuit diffusion in materials

This paper reviews the recent progress in understanding the atomic mechanisms of short-circuit diffusion along materials interfaces, such as grain and interphase boundaries, as well as lattice and interfacial dislocations/disconnections. Recent atomistic computer simulations have shown that short-circuit diffusion is dominated by collective atomic rearrangements in the form of strings and rings of mobile atoms. The process is dynamically heterogeneous in space and time and has many features in common with atomic dynamics in supercooled glass-forming liquids. We discuss examples of grain boundary, interphase boundary, and dislocation diffusion in metals and alloys, including the solute effect on the diffusion rates and mechanisms. Interphase boundaries are exemplified by Al–Si interfaces with diverse orientation relationships and atomic structures. The hierarchy of short-circuit diffusion paths in materials is reviewed by comparing the rates of grain boundary, interphase boundary, and dislocation diffusion. Future directions in the field of short-circuit diffusion in defect core regions are discussed.

36 MATERIALS SCIENCE↗

Coupling Finite Element and Finite Volume Simulation Within MOOSE

This work focuses on adding capability to couple finite element (FE) variables into finite volume (FV) physics within INL’s Multiphysics Object-Oriented Simulation Environment (MOOSE). This coupling can allow for improvement of multiphysics simulations where one set of physics is best suited for an FE discretization, while another set is best suited for an FV discretization. Electrohydrodynamics, which govern liquid metal reactor concepts and plasma dynamics, is a fitting example where the electromagnetic field equations are solved with FE and the fluid flow is solved with FV. The new FE to FV coupling method can be summarized as taking the element or face average of the FE variable value or gradient and applying that quantity directly in the FV equation objects.

97 MATHEMATICS AND COMPUTING↗

Molecular simulation of the separation of toluene and p-xylene with the thermally-robust ionic liquid triphenyl-p-phenyl sulfonyl phenyl phosphonium

In this work, we use molecular dynamics simulations to evaluate the behavior of toluene (TOL) and p-xylene (XYL) within three different solvents: n -hexane (HEX), n -heptane (HEP), and a thermally-robust ionic liquid (triphenyl-p-phenyl sulfonyl phenyl phosphonium + bis(trifluoromethylsulfonyl)imide (TPSP + Tf 2 N)). Several different temperatures are explored, while the solvation structures, energetics, and dynamics of TOL and XYL are evaluated. The solvation free energy (ΔG solv ) and transfer free energy (ΔG transfer ) of TOL and XYL are calculated using thermodynamic integration (TI). The transfer free energy of XYL from HEP/HEX to TPSP + Tf 2 N is more favorable than that of TOL over the entire temperature range, while the transfer of both TOL and XYL becomes less favorable as the temperature increases. The diffusion rates of TOL and XYL are lower in the ionic liquid when compared to HEP and HEX, but the thermal stability of TPSP + Tf 2 N may allow for higher operating temperatures and accelerated mass transport rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reversible phase transformations between Pb nanocrystals and a viscous liquid-like phase

Phase transformations have been a prominent topic of study for both fundamental and applied science. Solid-liquid reaction–induced phase transformations can be hard to characterize, and the transformation mechanisms are often not fully understood. Here, we report reversible phase transformations between a metal (Pb) nanocrystal and a viscous liquid-like phase unveiled by in situ liquid cell transmission electron microscopy. The reversible phase transformations are obtained by modulating the electron current density (between 1000 and 3000 electrons Å –2 s –1 ). The metal-organic viscous liquid-like phase exhibits short-range ordering with a preferred Pb-Pb distance of 0.5 nm. Assisted by density functional theory and molecular dynamics calculations, we show that the viscous liquid-like phase results from the reactions of Pb with the CH 3 O fragments from the triethylene glycol solution under electron beam irradiation. Such reversible phase transformations may find broad implementations.

36 MATERIALS SCIENCE↗

Role of phonon softening induced by anisotropic fluctuations in the enhanced mobility at free glassy surfaces

The surface of a glassy material exhibits enhanced mobility compared to the bulk counterpart, however the underlying mechanism for this remains elusive. Herein, we present studies of the dynamical properties of a prototypical glass-forming metallic liquid Zr 50 Cu 50 as a function of the distance from both the free surface and pinned surface using molecular dynamics simulations. We found that the surface mobility increases gradually on approaching the free surface, with a concomitant increase of the non-Gaussianity. The phonon density of states at the free surface exhibits lower characteristic frequencies than in the bulk and pinned surface. These findings suggest phonon softening caused by anisotropic fluctuations at free surfaces as an alternative physical mechanism leading to the enhanced dynamics at free glassy surfaces from the perspective of collective excitations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Bayesian calibration of bubble size dynamics applied to CO 2 gas fermenters

To accelerate the scale-up of gaseous CO 2 fermentation reactors, computational models need to predict gas-to-liquid mass transfer which requires capturing the bubble size dynamics, i.e. bubble breakup and coalescence. However, the applicability of existing models beyond air–water mixtures remains to be established. Here, an inverse modeling approach, accelerated with a neural network surrogate, calibrates the breakup and coalescence closure models, that are used in class methods for population balance modeling (PBM). The calibration is performed based on experimental results obtained in a CO 2 -air–water-coflowing bubble column reactor. Bayesian inference is used to account for noise in the experimental dataset and bias in the simulation results. To accurately capture gas holdup and interphase mass transfer, the results show that the breakage rate needs to be increased by one order of magnitude. In conclusion, the inferred model parameters are then used on a separate configuration and shown to also improve bubble size distribution predictions.

09 BIOMASS FUELS↗

Understanding Relationships between Free Volume and Oxygen Absorption in Ionic Liquids

In this report understanding the factors that govern gas absorption in ionic liquids is critical to the development of high-capacity solvents for catalysis, electrochemistry, and gas separations. Here, we report experimental probes of liquid structure that provide insights into how free volume impacts the O 2 absorption properties of ionic liquids. Specifically, we establish that isothermal compressibility - measured rapidly and accurately through small-angle X-ray scattering - reports on the size distribution of transient voids within a representative series of ionic liquids and is correlated with O 2 absorption capacity. Additionally, O 2 absorption capacities are correlated with thermal expansion coefficients, reflecting the beneficial effect of weak intermolecular interactions in ionic liquids on free volume and gas absorption capacity. Molecular dynamics simulations show that the void size distribution - in particular, the probability of forming larger voids within an ionic liquid - has a greater impact on O 2 absorption than the total free volume. These results establish relationships between the ionic liquid structure and gas absorption properties that offer design strategies for ionic liquids with high gas solubilities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Anharmonicity and Phase Diagram of Magnesium Oxide in the Megabar Regime

With density functional molecular dynamics simulations, we computed the phase diagram of MgO from 50 to 2000 GPa up to 20000 K. Via thermodynamic integration (TDI), we derive the Gibbs free energies of the B1, B2, and liquid phases and determine their phase boundaries. With TDI and a pseudo-quasi-harmonic approach, we show that anharmonic effects are important and stabilize the B1 phase in particular. As a result, the B1-B2 transition boundary in the pressure-temperature plane exhibits a steep slope. We predict the B1-B2-liquid triple point to occur at approximately T= 10000 K and P = 370 GPa, which is higher in pressure than was inferred with quasi-harmonic methods alone. We predict the principal shock Hugoniot curve to enter the B2 phase stability domain but only over a very small range of parameters. Here, this may render it difficult to observe this phase with shock experiments because of kinetic effects.

3-dimensional systems↗

Generalized Analysis of Dynamic Flow Fouling on Heat Transfer Surfaces

Several energy applications involve flow of untreated heat transfer fluids that contain trace amounts of mineral salts, such as calcium carbonate and calcium sulfate, which show reverse solubility with temperature. The gradual accretion of these salts onto the heat transfer surfaces causes fouling that degrades thermal performance over time. While fouling of bare surfaces has been extensively studied, fouling of nonwetting superhydrophobic or liquid-infused surfaces is less understood, especially under dynamic flow conditions. This study presents a systematic experimental analysis to elucidate dynamic flow fouling of calcium carbonate and calcium sulfate on smooth, superhydrophobic and liquid-infused surfaces in a generalized manner. Here, copper tubes with modified inner surface wettability are fabricated and subjected to forced convection at different flow rates and fouling salt concentrations. Fouling characteristics of the two salts on the different surfaces are examined through microstructural analysis, and quantified in terms of asymptotic fouling resistance that measures the total extent of fouling, and a threshold time, which correlates to the onset of significant fouling. Using Taguchi design of experiments, the study reports, for the first time, closed-form analytical relationships for the asymptotic fouling thermal resistance and threshold time on the Reynolds number, dimensionless concentration of fouling agent and dimensionless infused liquid viscosity that represents the different surface types in a unified manner. It is shown that the analytical model accurately predicts the fouling characteristics in up to about 97% of the data, based on which contour maps of optimum surface designs for minimizing fouling resistance are presented. The results provide information on designing surfaces for reduced fouling or for estimating the fouling characteristics for a given surface in applications.

36 MATERIALS SCIENCE↗

Rational Design of Weakly‐Solvating Molecules for Salt‐In‐Pre‐Ionic‐Liquid Electrolytes for Li Metal Batteries

Lithium metal batteries (LMBs) promise step‐changes in energy densities but suffer from poor cycle life due to unstable electrolyte‐lithium interfaces. Conventional carbonate electrolytes exhibit excessive lithium‐ion solvation and low oxidative stability, leading to rapid capacity loss. Herein, we report a rationally designed weakly‐solvating cyclic sulfonamide, 1‐trifluoromethanesulfonyl)amide pyrrolidine (TFMSPyr), which integrates an electron‐withdrawing trifluoromethanesulfonyl functional group at pyrrolidinic‐N. TFMSPyr acts as a pre‐ionic‐liquid solvent that forms intrinsically localized, anion‐dominated solvation, coupling molecular architecture, solvation topology, and transport dynamics. As a result, LiFSI based salt‐in‐pre‐ionic‐liquid (SIPIL) electrolytes exhibit high lithium‐ion transference number, oxidative stability > 5 V versus Li/Li + and anion‐derived solid electrolyte interphases (SEI). Li||Cu cells with SIPIL deliver a first cycle Coulombic efficiency (CE) of ≈ 99% with average CE of 99.2% for 100 cycles, and lithium half‐cells with lithium iron phosphate (LFP) cathode exhibit 82% capacity retention after 400 cycles with CE of 99.98%. In anode‐free full cells, 95% of initial capacity is retained after 63 cycles with an average CE of 99.5%. These results demonstrate that molecular engineering of solvents offers a powerful pathway to stabilize lithium metal interfaces and enable practical Anodeless LMBs.

25 ENERGY STORAGE↗

Growth of Hexagonal Boron Nitride from Molten Nickel Solutions: A Reactive Molecular Dynamics Study

Metal flux methods are excellent for synthesizing high-quality hexagonal boron nitride (hBN) crystals, but the atomic mechanisms of hBN nucleation and growth in these systems are poorly understood and difficult to probe experimentally. Here, we harness classical reactive molecular dynamics (ReaxFF) to unravel the mechanisms of hBN synthesis from liquid nickel solvent over time scales up to 30 ns. These simulations mimic experimental conditions by including relatively large liquid nickel slabs containing dissolved boron and a molecular nitrogen gas phase. Overall, the reaction takes place almost exclusively on the surface of the liquid nickel, owing to the low solubility of nitrogen in bulk nickel and the intermediate species’ preference for the metal–gas interface. The formation of hBN invariably begins by reaction of dinitrogen with nickel-solvated boron atoms at the surface, forming intermediate N–N–B species, which typically evolve into B–N–B units through a short-lived intermediate where a single nitrogen atom is coordinated by one nitrogen and two boron atoms. The resulting B–N–B units, in turn, coalesce with growing hBN nuclei and carry nitrogen between hBN nanocrystals in an Ostwald ripening process. The amount of hBN produced on the tens of nanosecond time scale depends critically on the boron concentration, while having a much weaker dependence on the N 2 pressure for the regime considered (N 2 pressures of 2.5–10 MPa, Ni–B solutions with 6–12% boron by atom fraction). The highest rate of hBN formation occurs at the lowest temperature considered (1750 K, just above the melting point of nickel), while no hBN sheets are formed at 2000 K or above. An analysis of the transition pathways for nitrogen atoms shows that the final step, incorporation of small B–N motifs into larger hBN sheets, is the rate-limiting step in the regimes considered. While raising the temperature from 1750 to 2000 K has little effect on the formation of intermediates (N–N–B, B–N–B, etc.), the lack of large hBN sheets at temperatures >1900 K is explained by decreased probability of the final step and increased probability of breakup of hBN into B–N motifs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Defect Spirograph: Dynamical Behavior of Defects in Spatially Patterned Active Nematics

We report that topological defects in active liquid crystals can be confined by introducing gradients of activity. Here, we examine the dynamical behavior of two defects confined by a sharp gradient of activity that separates an active circular region and a surrounding passive nematic material. Continuum simulations are used to explain how the interplay among energy injection into the system, hydrodynamic interactions, and frictional forces governs the dynamics of topologically required self-propelling + 1/2 defects. Our findings are rationalized in terms of a phase diagram for the dynamical response of defects in terms of activity and frictional damping strength. Different regions of the underlying phase diagram correspond to distinct dynamical modes, namely immobile defects, steady rotation of defects, bouncing defects, bouncing-cruising defects, dancing defects, and multiple defects with irregular dynamics. These dynamic states raise the prospect of generating synchronized defect arrays for microfluidic applications.

36 MATERIALS SCIENCE↗

Dynamics of Emim + in [Emim][TFSI]/LiTFSI Solutions as Bulk and under Confinement in a Quasi-liquid Solid Electrolyte

Quasi-liquid solid electrolytes are a promising alternative for next-generation Li batteries. Furthermore, these systems combine the safety of solid electrolytes with the desired properties of liquids and are typically formed by solutions of Li salts in ionic liquids incorporated into solid matrices. Here, we present a fundamental understanding of the transport properties in solutions of lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) in 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([Emim][TFSI]), either in bulk form or incorporated in a boron nitride (BN) matrix. We performed a series of quasi-elastic neutron scattering experiments that, given the high incoherent neutron scattering cross section of hydrogen, allowed us to focus on the Emim + dynamics. First, [Emim][TFSI]/LiTFSI solutions (0.5 and 2.5 mol·kg –1 ) were investigated and we show how the increase in the concentration reduces the Emim+ mobility and increases the activation energy of their long-range motions. Then, the 0.5 mol·kg –1 solution was incorporated into the BN matrix and we report that the diffusivities of the Emim + cations that remain mobile under confinement are highly accelerated in comparison with the bulk sample and the activation energy of these motions is drastically reduced. We present the experimental evidence that this effect is related to the content of the Emim+ cations immobilized near the surfaces of the BN pores.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modelling and analysis of nuclear reactor system coupled with a liquid metal battery

Traditionally, nuclear power plants in the U.S. provide baseload power to the power grid because they have less flexibility for ramping their output power than natural gas peaking plants. However, achieving climate goals to reduce the consumption of fossil‐based natural gas places pressure on nuclear power plants and other power generators to ramp up their power output to balance grid generation with demand. This paper presents the modelling and performance analysis of a nuclear reactor system (NRS) coupled to a liquid‐metal battery (LMB) to improve its dynamic response and enable its black start capability. The NRS and LMB thermal behaviour are modelled in Dymola, while the electrical dynamics of the LMB and power grid are modelled in RTDS‐RSCAD. Both simulation platforms are coupled and share their thermal and electrical data using a Transmission Control Protocol/Internet Protocol (TCP/IP) communication protocol. The dynamic performance of the NRS‐LMB integration is tested on the IEEE 9 bus, which demonstrates its ability to respond and provide frequency and voltage regulation. The black start capability of the NRS‐LMB is also evaluated by simulating a grid outage and using the LMB to supply the auxiliary loads required to bring the NRS back online as soon as possible. The results show that coupling an NRS to an LMB improves the system dynamic performance and enables it to black start after being disconnected from the grid for several days.

25 ENERGY STORAGE↗

Sound attenuation in the hyperhoneycomb Kitaev spin liquid

In recent years, it has been shown that the phonon dynamics may serve as an indirect probe of fractionalization of spin degrees of freedom. Furthermore we propose that the sound attenuation measurements allows for the characterization and identification of the Kitaev quantum spin liquid on the hyperhoneycomb lattice, which is particularly interesting since the strong Kitaev interaction was observed in the the hyperhoneycomb magnet β–Li 2 IrO 3 . To this end we consider the low-temperature scattering between acoustic phonons and gapless Majorana fermions with nodal-line band structure. Here, we find that the sound attenuation has a characteristic angular dependence, which is explicitly shown for the high-symmetry planes at temperatures below the flux energy gap.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electron and Electron-Nuclear Dynamics of Open Quantum Subsystems

Leveraging an open-subsystem formulation of Density Functional Theory (DFT) we aim at describing periodic and molecular systems alike, including their electronic and nuclear dynamics. Subsystem DFT enables first principles simulations to approach realistic time- and length-scales, and most importantly sheds light on the dynamical behavior of complex systems. Taking subsystem DFT to the time domain allows us to inspect the electron dynamics of condensed-phase systems in real time. In liquids and interfaces, we observe all the relevant regimes proper of non-Markovian open quantum system dynamics, such as electronic energy transfer, and screening. In addition, the ab-initio modeling of system-bath interactions brought us to observe and justify the holographic time-dependent electron density theorem. Contrary to interactions between molecular (finite) systems, when molecules interact with metal or semiconductor surfaces the electron dynamics is strongly non-Markovian with dramatic repercussions to the molecule’s response to external perturbations. Metals and semiconductors typically have large polarizabilities, and even in a regime of low coupling their effect on impinging molecular species is significant–line broadening, peak shift, and intensity borrowing are observed, characterized, and explained in terms of inter-subsystem dynamical interactions and a many-body decomposition of the system’s density-density response function in a way that transcends the canons of Fermi Golden Rule.

33 ADVANCED PROPULSION SYSTEMS↗

Holographic non-Fermi liquids at large d

Motivated by the results of Dynamical Mean Field Theory, we study the two-point function of fermions moving in a charged black brane background in AdS d +1 in the limit of large d . We observe the emergence of a locally critical form of the fermion self-energy, with a strongly constrained range of possible scaling behaviors at large d . Novelties compared to the analysis in d = 3, 4 include an enlarged regime of temperatures where the results apply, and the analytical tractability of the calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗