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,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Predicted Melt Curve and Liquid Shear Viscosity of RDX up to 30 GPa

Recent grain-scale simulations of HMX and TATB have shown that predictions for hot spot formation in high explosives are particularly sensitive to accurate determinations of the pressure-dependent melt curve and the shear viscosity of the liquid phase. These physics terms are poorly constrained beyond ambient pressure for the explosive RDX. Here we adopt an all-atom modeling approach using molecular dynamics (MD) simulations to predict the melt curve of RDX near to detonation conditions (30 GPa) and determine the shear viscosity of the liquid as a function of temperature and pressure above the melt curve. Phase-coexistence simulations were used to determine the melt curve, which is predicted to vary by almost 1100 K as the pressure increases from 0 GPa to 30 GPa. Equilibrium MD simulations and the Green-Kubo formalism were used to obtain the pressure-temperature-dependent shear viscosity. The shear viscosity of RDX is predicted to be of similar magnitude to the viscosity of TATB at low GPa-range pressures, and to be roughly an order of magnitude lower than the viscosity of HMX. The temperature dependence of the shear viscosity is Arrhenius at a given pressure, and the exponential pre-factor and activation term exhibit a strong, yet complicated, pressure dependence. An empirical pressure-temperature dependent function for RDX shear viscosity is developed that simultaneously captures a wide range of MD predictions while taking an analytic form that extrapolates smoothly beyond the fitted regime. The relative strength of the pressure and temperature dependencies of these two physics terms is found to be of similar magnitude for RDX, HMX, and TATB, which motivates incorporating these results in future RDX grain-scale modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Critical fluid dynamics in two and three dimensions

We describe a numerical method for simulating stochastic fluid dynamics near a critical point in the Ising universality class. This theory is known as model H, and is expected to govern the nonequilibrium dynamics of quantum chromodynamics (QCD) near a possible critical endpoint of the phase transition between a hadron liquid and the quark-gluon plasma. The numerical algorithm is based on a Metropolis scheme, and automatically ensures that the distribution function of the hydrodynamic variables in equilibrium is independent of the transport coefficients and only governed by the microscopic free energy. We verify dynamic scaling near the critical point of a two and three-dimensional fluid and extract the associated critical exponent z. Here, we find z≃3 in three dimensions, and z≃2 for a two-dimensional fluid. In a finite system, we observe a crossover between the mean field value z=4 and the true critical exponent z≃3 (z≃2 in d=2). This crossover is governed by the values of the correlation length and the renormalized shear viscosity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Characterizing the seismic response of a molten salt nuclear reactor

Abstract One Generation IV nuclear reactor, which uses a fluoride salt as a coolant, graphite reflector blocks as a moderator, and circulating buoyant TRISO pebbles as fuel is at an advanced stage of development. To characterize the seismic behavior of components of this reactor, validate numerical models for analysis, and develop recommendations for design, a set of earthquake‐simulator experiments on a scaled model of the reactor vessel and its internals was executed on a six‐degree‐of‐freedom earthquake simulator. The model was seismically isolated at its base using two types of spherical sliding bearings. The scaled model involved representations of the prototype reactor vessel, core barrel, reflector blocks, coolant, and spherical fuel pebbles. The material and geometric properties of different test components were selected based on a dynamic similitude scaling analysis and an approximate length scale of 0.4. Four sets of three‐component earthquake motions were used as inputs for testing. Instrumentation on the test specimen recorded the dynamic responses of the outer vessel, core barrel, and reflector‐block assembly, the hydrodynamic responses (sloshing and hydrodynamic pressure) of the liquid coolant, pebble consolidation under earthquake shaking, and the behavior of the isolation systems. This paper describes the design of the experiments and presents key results from the tests. The dynamic responses of the outer vessel, core barrel, and the reflector blocks revealed that the components responded as a unit for the intense shaking used in the experiments. The sloshing response of the fluid in a thin annulus near the perimeter of the vessel was heavily damped. The change in the packing fraction of the pebble bed under repeated, intense 3D earthquake shaking was less than 3%. Seismically isolating the vessel substantially reduced demands on its internal components.

Engineering↗
DOE OSTI2025

The Dynamical Role of Optical Phonons and Sublattice Screening in a Solid-State Ion Conductor

Solid-state electrolytes (SSEs) require ionic conductivities that are competitive with liquid electrolytes to realize applications in all-solid-state batteries. Although candidate SSEs have been discovered, the underlying mechanisms enabling superionic conduction (>1 mS cm –1 ) remain elusive. In particular, the role of ultrafast lattice dynamics in mediating ion migration, which involves couplings between ions, phonons, and electrons, is rarely explored experimentally at their corresponding time scales. To investigate the complex contributions of coupled lattice dynamics on ion migration, we modulate the charge density occupations within the crystal framework and then measure the time-resolved change in impedance on picosecond time scales for a candidate SSE, Li 0.5 La 0.5 TiO 3 (LLTO). Upon perturbation, we observe enhanced ion migration at ultrafast time scales. The respective transients match the time scales of optical and acoustic phonon vibrations, suggesting their involvement in ion migration. We further computationally evaluate the effect of a charge transfer from the O 2p to the Ti 3d band on the electronic and physical structure of LLTO. We hypothesize that the charge-transfer excitation distorts the TiO 6 polyhedra by altering the local charge density occupancy of the hopping site at the migration pathway saddle point, thereby causing a reduction in the migration barrier for the Li + hop. We rule out the contribution of photogenerated electron carriers and laser heating. Overall, our investigation introduces a new spectroscopic tool to probe fundamental ion hopping mechanisms transiently at ultrafast time scales, which has previously only been achieved in a time-averaged manner or solely via computational methods.

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

Infrared absorption spectroscopy of dynamically compressed water

Streaked visible spectroscopy is well established in dynamic compression research. Infrared measurements remain problematic, however, due to the diminishing sensitivity of streak camera photocathodes beyond 800 nm. Time-stretch techniques offer an alternative method for probing infrared features during single-event experiments. This paper discusses the development of a time-stretch spectroscopy diagnostic using dispersed supercontinuum laser pulses. The technique is applied to near-infrared measurements of liquid water during multiple shock compression.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRSDate not supplied

Modeling Self-Pressurization and Spray Bar Pressure Control of a Cryogenic Storage Tank in Normal Gravity

This paper presents Computational Fluid Dynamics (CFD) models for simulating self-pressurization and spray bar cooling processes in a large-scale liquid hydrogen storage tank under normal gravity. The self-pressurization model utilizes the kinetics-based Schrage equation with the Volume-Of-Fluid (VOF) model to account for the evaporative and condensing interfacial mass fluxes. It uses a laminar approach to model natural convection and compares the pressures and temperatures predicted by the CFD model with experimental data from the Multipurpose Hydrogen Test Bed (MHTB) self-pressurization experiment. Additionally, the predicted interfacial mass flux is presented. A CFD model for simulating pressure control using a spray bar is also presented. An Eulerian-Lagrangian approach models the interactions between the discrete droplets and the continuous ullage phase. The spray model is coupled with the VOF model by tracking particles in the ullage, removing particles from the ullage when they reach the interface, and then adding their contributions to the liquid. The T-sat model is used to calculate the droplet-ullage heat and mass transfer. In this model, droplets warm up to the saturation temperature corresponding to the ullage vapor pressure thereafter evaporating while remaining at the saturation temperature. The tank pressure, vapor, and liquid temperature evolutions are validated against data from the MHTB spray bar mixing experiment.

Self-Pressurization↗
DOE OSTI2023

Dynamics of polymer electrolyte with LiTFSI via Quasi-Elastic Neutron Scattering

Most lithium batteries offer a wide range of applications. However, safety issues are still an unresolved issue for several applications. To solve the safety issue of Li-ion batteries, solid polymer electrolyte is a promising candidate to replace commercial liquid electrolyte. A 4-arm star poly(ethylene oxide) polymer with LiTFSI salt as an electrolyte was studied. The dynamics of this polymer were explored with the Quasi-Elastic Neutron Scattering technique. Furthermore, the influence of temperature and Li salt concentration on the polymer dynamics was investigated. The dynamics of the polymer ends of the arm show much higher flexibility than the core parts making those types of polymers attractive for further studies in battery research.

Wang, Hui↗
DOE OSTI2020

Machine-learned interatomic potentials by active learning: amorphous and liquid hafnium dioxide

We propose an active learning scheme for automatically sampling a minimum number of uncorrelated configurations for fitting the Gaussian Approximation Potential (GAP). Our active learning scheme consists of an unsupervised machine learning (ML) scheme coupled with a Bayesian optimization technique that evaluates the GAP model. We apply this scheme to a Hafnium dioxide (HfO 2 ) dataset generated from a "melt-quench" ab initio molecular dynamics (AIMD) protocol. Our results show that the active learning scheme, with no prior knowledge of the dataset, is able to extract a configuration that reaches the required energy fit tolerance. Further, molecular dynamics (MD) simulations performed using this active learned GAP model on 6144 atom systems of amorphous and liquid state elucidate the structural properties of HfO 2 with near ab initio precision and quench rates (i.e., 1.0 K/ps) not accessible via AIMD. The melt and amorphous X-ray structural factors generated from our simulation are in good agreement with experiment. In addition, the calculated diffusion constants are in good agreement with previous ab initio studies.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Overscreening, Coion-Dominated Electroosmosis, and Electric Field Strength Mediated Flow Reversal in Polyelectrolyte Brush Functionalized Nanochannels

Controlling the direction and strength of nanofluidic electrohydrodyanmic transport in the presence of an externally applied electric field is extremely important in a number of nanotechnological applications. Here, we employ all-atom molecular dynamics simulations to discover the possibility of changing the direction of electroosmotic (EOS) liquid flows by merely changing the electric field strength in a nanochannel functionalized with polyelectrolyte (PE) brushes. In exploring this, we have uncovered three facets of nanoconfined PE brush behavior and resulting EOS transport. First, we identify the onset of an overscreening effect: such overscreening refers to the presence of more counterions (Na + ) within the brush layer than needed to neutralize the negative brush charges. Accordingly, as a consequence of the overscreening, in the bulk liquid outside the brush layer, there is a greater number of co-ions (Cl – ) than counterions in the presence of an added salt (NaCl). Second, this specific ion distribution ensures that the overall EOS flow is along the direction of motion of the co-ions. Such co-ion-dictated EOS transport directly contradicts the notion that EOS flow is always dictated by the motion of the counterions. Finally, for large-enough electric fields, the brush height reduces significantly, causing some of the excess overscreening-inducing counterions to squeeze out of the PE brush layer into the brush-free bulk. As a result, the overscreening effect disappears and the number of co-ions and counterions outside the PE brush layer become similar. Despite that there is an EOS transport, this EOS transport, unlike the standard EOS transport that occurs due to the imbalance of the co-ions and counterions, occurs since a larger residence time of the water molecules in the first solvation shell of the counterions (Na + ) ensures a water transport in the direction of motion of the counterions. In this work, the net effect is the reversal of the direction of the EOS transport by merely changing the strength of the electric field.

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

Direct Observation of Ultrafast Hydrogen Bond Strengthening in Liquid Water

Water is one of the most important, yet least understood, liquids in nature. Many anomalous properties of liquid water originate from its well-connected hydrogen bond network, including unusually efficient vibrational energy redistribution and relaxation. An accurate description of the ultrafast vibrational motion of water molecules is essential for understanding the nature of hydrogen bonds and many solution-phase chemical reactions. Most existing knowledge of vibrational relaxation in water is built upon ultrafast spectroscopy experiments. However, these experiments cannot directly resolve the motion of the atomic positions and require difficult translation of spectral dynamics into hydrogen bond dynamics. Here we measured the ultrafast structural response to the excitation of the OH stretching vibration in liquid water with femtosecond temporal and atomic spatial resolution using liquid ultrafast electron scattering. We observed a transient hydrogen bond contraction of roughly 0.04 Å on a timescale of 80 femtoseconds, followed by a thermalization on a timescale of ~1 picosecond. Molecular dynamics simulations reveal the need to treat the distribution of the shared proton in the hydrogen bond quantum mechanically to capture the structural dynamics on the femtosecond timescales. Our experiment and simulations unveiled the intermolecular character of the water vibration preceding the relaxation of the OH stretch.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The effect of blue and infrared laser melting frequency on oxide morphology in 304L

Continuous Wave Laser Beam (LB) melting offers control over the localised heating and cooling of melt pools for welding, brazing, additive manufacturing and solidification. Research into laserliquid metal interactions has primarily focused on the heat and mass transport under large thermal gradients imposed by the localised melting conditions. Recent commercial availability of blue (450 nm) wavelength lasers has enabled the printing of material with poor IR (1064) nm attenuation like Cu. However, little research has been conducted into varying input laser frequency and scan rate at fixed absorption to understand the effects of laser light on surface oxide formation and the underlying liquid metal thermodynamic response. Here, this article conducts laser melting of 304L under Blue (450 nm) and Infrared (IR, 1064 nm) laser frequencies and examines their impact on laser oxide thickness, chemistry and coloration. We find that laser frequencies induce changes in the oxide layer thickness and chemistry that cannot be explained using conventional thermal absorption shifts and fluid dynamics. We suggest that light coupling may have thermodynamic implications for the chemical potential of the liquid metal, which may drive the observed phase behaviour. Identifying mechanisms for the tuning of the chemical potential via laser frequency control could enable the ability to control solidification microstructure and atomic ordering in laser welding and additive manufacturing processes.

304L↗
NASA NTRS1992

Liquid propulsion turbomachinery model testing

For the past few years an extensive experimental program to understand the fluid dynamics of the Space Shuttle Main Engine hot gas manifold has been in progress. This program includes models of the Phase II and II+ manifolds for each of the air and water flow facilities, as well as two different turbine flow paths and two simulated power levels for each manifold. All models are full-scale (geometric). The water models are constructed partially of acrylic to allow flow visualization. The intent of this paper is to discuss the concept, including the test objectives, facilities, and models, and to summarize the data for an example configuration, including static pressure data, flow visualization, and the solution of a specific flow problem.

Mcdaniels, David M.↗
NASA NTRS1993

Liquid propellant rocket engine combustion simulation with a time-accurate CFD method

Time-accurate computational fluid dynamics (CFD) algorithms are among the basic requirements as an engineering or research tool for realistic simulations of transient combustion phenomena, such as combustion instability, transient start-up, etc., inside the rocket engine combustion chamber. A time-accurate pressure based method is employed in the FDNS code for combustion model development. This is in connection with other program development activities such as spray combustion model development and efficient finite-rate chemistry solution method implementation. In the present study, a second-order time-accurate time-marching scheme is employed. For better spatial resolutions near discontinuities (e.g., shocks, contact discontinuities), a 3rd-order accurate TVD scheme for modeling the convection terms is implemented in the FDNS code. Necessary modification to the predictor/multi-corrector solution algorithm in order to maintain time-accurate wave propagation is also investigated. Benchmark 1-D and multidimensional test cases, which include the classical shock tube wave propagation problems, resonant pipe test case, unsteady flow development of a blast tube test case, and H2/O2 rocket engine chamber combustion start-up transient simulation, etc., are investigated to validate and demonstrate the accuracy and robustness of the present numerical scheme and solution algorithm.

Chen, Y. S.↗
NASA NTRS2009

Results of the Fluid Merging Viscosity Measurement International Space Station Experiment

The purpose of FMVM is to measure the rate of coalescence of two highly viscous liquid drops and correlate the results with the liquid viscosity and surface tension. The experiment takes advantage of the low gravitational free floating conditions in space to permit the unconstrained coalescence of two nearly spherical drops. The merging of the drops is accomplished by deploying them from a syringe and suspending them on Nomex threads followed by the astronaut s manipulation of one of the drops toward a stationary droplet till contact is achieved. Coalescence and merging occurs due to shape relaxation and reduction of surface energy, being resisted by the viscous drag within the liquid. Experiments were conducted onboard the International Space Station in July of 2004 and subsequently in May of 2005. The coalescence was recorded on video and down-linked near real-time. When the coefficient of surface tension for the liquid is known, the increase in contact radius can be used to determine the coefficient of viscosity for that liquid. The viscosity is determined by fitting the experimental speed to theoretically calculated contact radius speed for the same experimental parameters. Recent fluid dynamical numerical simulations of the coalescence process will be presented. The results are important for a better understanding of the coalescence process. The experiment is also relevant to liquid phase sintering, free form in-situ fabrication, and as a potential new method for measuring the viscosity of viscous glass formers at low shear rates.

Ethridge, Edwin C.↗
DOE OSTI2023

Kinematics and dynamics of disclination lines in three-dimensional nematics

An exact kinematic law for the motion of disclination lines in nematic liquid crystals as a function of the tensor order parameter Q is derived. Unlike other order parameter fields that become singular at their respective defect cores, the tensor order parameter remains regular. Following earlier experimental and theoretical work, the disclination core is defined to be the line where the uniaxial and biaxial order parameters are equal, or equivalently, where the two largest eigenvalues of Q cross. This allows an exact expression relating the velocity of the line to spatial and temporal derivatives of Q on the line, to be specified by a dynamical model for the evolution of the nematic. By introducing a linear core approximation for Q, analytical results are given for several prototypical configurations, including line interactions and motion, loop annihilation, and the response to external fields and shear flows. Behaviour that follows from topological constraints or defect geometry is highlighted. Finally, the analytic results are shown to be in agreement with three-dimensional numerical calculations based on a singular Maier–Saupe free energy that allows for anisotropic elasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Changes in short- and medium-range order in metallic liquids during undercooling

It has been widely speculated that dominant motifs, such as short-range icosahedral order, can influence glass formation and the properties of glasses. Experimental data on both fragile and strong undercooled liquids demonstrate corresponding changes in their thermophysical properties consistent with increasing development of a network of interconnect motifs based on molecular dynamics. Describing these regions of local order, how they connect, and how they are related to property changes have been challenging issues, both computationally and experimentally. Yet the consensus is that metallic liquids develop interconnected medium-range order consisting of some regions with lower mobility with deeper undercooling. Less well understood is how these motifs (or “crystal genes”) in the liquid can inhibit nucleation in the deeply undercooled liquid or influence phase selection upon devitrification. These motifs tend to have local packing unlike stable compounds with icosahedral order tending to dominate the best glass formers. The underlying kinetic and thermodynamic forces that guide the formation of these motifs and how they interconnect during undercooling remain open questions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The CP-PAW Code Package for First-Principles Calculations from a User’s Perspective

CP-PAW is a combined electronic structure and ab initio molecular dynamics code to perform mixed quantum and classical simulations of atomistic condensed phase systems, such as solids, liquids, and molecular systems. As the name suggests, the CP-PAW code unifies the all-electron projector augmented-wave (PAW) method with the Car–Parrinello (CP) approach to determine not only the electronic and nuclear ground states of condensed matter but also to study their properties and dynamics. In addition to briefly outlining the underlying theory, the focus will be on the unique aspects of CP-PAW and how to correctly employ them as a user. How to install CP-PAW using the new build system will also be briefly mentioned.

Blöchl, Peter E [Institute for Theoretical Physic↗
DOE OSTI2020

Reversible Electrochemical Phase Change in Monolayer to Bulk-like MoTe 2 by Ionic Liquid Gating

Transition-metal dichalcogenides (TMDs) exist in various crystal structures with semiconducting, semi-metallic, and metallic properties. The dynamic control of these phases is of immediate interest for next-generation electronics such as phase change memories. Of the binary Mo and W-based TMDs, MoTe 2 is attractive for electronic applications because it has the lowest energy difference (40 meV) between the semiconducting (2H) and semi-metallic (1T') phases, allowing for MoTe 2 phase change by electrostatic doping. Here, we report phase change between the 2H and 1T' polymorphs of MoTe 2 in thicknesses ranging from the monolayer to bulk-like case (73 nm) using an ionic liquid electrolyte at room temperature and in air. We find consistent evidence of a partially reversible 2H-1T' transition using in situ Raman spectroscopy where the phase change occurs in the topmost layers of the MoTe 2 flake. We find a thickness-dependent transition voltage where higher voltages are necessary to drive the phase change for thicker flakes. We also show evidence of electrochemical activity during the gating process by observation of Te metal formation. This finding suggests the formation of Te vacancies which have been reported to lower the energy difference between the 2H and 1T' phases, potentially aiding the phase change process. Furthermore, our discovery that the phase change can be achieved on the surface layer of bulk-like materials reveals that this electrochemical mechanism does not require isolation of a single layer and the effect may be more broadly applicable than previously thought.

2d materials↗