Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase velocity”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Elucidating the Microscale Behavior and Phase Separation Kinetics of Thermally Responsive Ionic Liquid–Water Mixtures

Thermally responsive ionic liquids (ILs) exhibit liquid-liquid phase separation into a water-rich (WR) and ionic-liquid-rich (ILR) phase when heated above a lower critical solution temperature (LCST). This phase behavior has been leveraged for applications ranging from forward osmosis (FO) desalination, where the IL acts as a draw solute, to refrigeration and dehumidification cycles, where the IL acts as a liquid desiccant. While significant effort has been devoted to characterizing the thermodynamic and thermophysical properties of LCST ILs, their phase separation kinetics have not been investigated. In this work, we describe the macroscale phase separation kinetics (phase separation time) by gleaning insight into the microscale colloidal behavior of aqueous mixtures of four different materials, P 4444 TFA (tetrabutylphosphonium-2,4-trifluoroacetate), P 4444 DMBS (tetrabutylphosphonium-2,4-dimethyl-benzenesulfonate), N 4444 Sal (tetrabutylammonium salicylate), and P 4444 Sal (tetrabutylphosphonium salicylate) as a function of IL concentration at a separation temperature of 70 °C. We report the discontinuous microscale size distributions for each material and correlate their theoretical settling velocities to experimental phase separation times. The results indicate that a simple Stokes' law model can predict the phase separation time within reasonable accuracy. Overall, this work lays the foundation for understanding the micro- to macroscale phase separation behavior and kinetics of LCST ILs for various water-energy applications.

LCST↗
DOE OSTI2021

Sequentially prepared Mo-V-Based SCR catalyst for simultaneous Hg 0 oxidation and NO reduction

Molybdenum (Mo)-vanadium (V)-based selective catalytic reduction (SCR) catalyst synthesized by the sequential impregnation of Mo and W followed by V was investigated for simultaneous elemental mercury (Hg 0 ) oxidation and nitrogen oxide (NO) reduction in an existing SCR unit with respect to different TiO 2 phases, calcination temperatures, flue gas constituents, gas velocities, and reaction temperatures. Anatase phase TiO 2 and the calcination temperatures of 400 and 500 degrees C resulted in similar to 69% Hg 0 oxidation at 10 ppmv HCl and 350 degrees C. The high calcination temperature of 700 degrees C resulted in TiO 2 phase transformation from anatase to rutile and agglomeration. The modified SCR catalyst prepared with the impregnation sequence of Mo and W followed by V using anatase TiO 2 and calcination temperature 500 degrees C showed similar to 99% Hg 0 oxidation and 87% NO reduction conversions at an NH 3 /NO molar ratio of 0.9 under 350 degrees C and 5,000 hr -1 space velocity in typical sub-bituminous and lignite coal simulated flue gases. Finally, the effects of these parameters and conditions were further investigated using various characterization techniques including BET, TEM, XRD, NH 3 TPD and XAFS.

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

Resolving Anodic Current and Temperature Distributions in a Polymer Electrolyte Membrane Water Electrolysis Cell Using a Pseudo-Two-Phase Computational Fluid Dynamics Model

Expanding upon our prior experimental work, we constructed a three-dimensional model of a polymer electrolyte membrane water electrolyzer using computational fluid dynamics. We applied the assumption of pseudo-two-phase flow, the flow of two phases with equal velocity. Experimental data were used to obtain parameters and to determine the conditions under which this model was valid. Anodic distributions of current density, temperature, liquid saturation, and relative humidity were obtained at various flow rates. The overall current density and temperature difference from inlet to outlet at the anode agreed strongly with experimental measurements under most circumstances. This verification allowed us to further examine the apparent gas coverage calculated from experimental and model temperature data. Results suggested a low liquid saturation and low relative humidity at the anode due to the consumption of liquid water and water vapor. However, we questioned the accuracy of the pseudo-two-phase assumption at low water feed rates. We concluded that the model was applicable to systems with liquid water feed rates greater than 0.6 ml min -1 cm -2 . Therefore, it is a fair screening method that can advise which operating conditions lead to excessive temperatures or drying at the anode, thereby promoting the longevity of the membrane and catalyst.

25 ENERGY STORAGE↗
DOE OSTI2024

Effect of the volume fraction gradient on the phase interaction force model for disperse two-phase flows

In this work, the effects of the particle volume fraction gradient on fluid-particle interactions are studied. The phase interaction force is decomposed into three terms. For the first term, namely the symmetrized force density, we present theoretical reasoning and numerical evidence to assume that it is independent of the particle volume fraction gradient. The second term is the particle volume fraction gradient times a newly introduced diffusion stress. The third term is the divergence of the particle-fluid-particle (PFP) stress. If this assumption of independence of the particle volume fraction gradient for the first term can be verified, to the first order of the ratio of the mean distance between particles to the macroscopic lengthscale, all three terms can be studied and modeled in flows with uniform particle distributions. Models thus obtained are applicable to statistically inhomogeneous flows, with the second and third terms accounting for statistical inhomogeneity. To verify this assumption, numerical simulations of flows passing fixed arrays of particles are performed. Both uniform and nonuniform particle volume fractions are studied and compared for disperse multiphase flows with the particle Reynolds numbers ranging from 1 to 100, and particle volume fraction ranging from 1% to 26% in statistically steady states. It is found that the symmetrized force (first) term can be well approximated by the drag force obtained from studies of uniform flows. The diffusion stress is positive along the flow direction and negative in the directions perpendicular to the flow. In the case of moving particles, this stress could potentially cause particle aggregation in the flow direction and dispersion in the directions perpendicular to the flow. Finally, the diffusion stress is only important when there is a volume fraction gradient, while the PFP stress can be important in inhomogeneous flows with either nonuniform particle concentrations or nonuniform average relative velocities between the phases.

42 ENGINEERING↗
DOE OSTI2021

CFD-DEM study of bubble properties in a cylindrical fluidized bed of Geldart Group D particles and comparison with prior MRI data

Numerical simulations of a freely bubbling cylindrical fluidized bed are carried out using a coupled computational fluid dynamics and discrete element method (CFD-DEM) model and compared to recent experimental data. The experiments were conducted using high-resolution and high-frequency magnetic resonance imaging providing high-fidelity data of the bubbling within a central 10 mm slice of the bed. Qualitatively, we nd more regular (i.e., less chaotic) structures in the simulated beds than observed experimentally. Quantitatively, however, the bubble diameter and number of bubbles as a function of height within the bed is predicted well by the base model. Unfortunately, the regularity in the simulations manifests as a considerable discrepancy in the speed of the (dense) emulsion phase. The simulated velocity probability distribution functions show an accumulation of low-speed regions and deficiency of high-speed regions. A simple parametric study of the base model is also carried out considering many of the most common CFD-DEM modeling parameters. It is found that the fluid grid size, geometry resolution, transfer kernel and drag law did not have a significant effect on bubble or particle dynamics.

42 ENGINEERING↗
DOE OSTI2022

Crust and Upper Mantle Structure Beneath the Eastern United States

Abstract The Eastern United States (EUS) has a complex geological history and hosts several seismic active regions. We investigate the subsurface structure beneath the broader EUS. To produce reliable images of the subsurface, we simultaneously invert smoothed P‐wave receiver functions, Rayleigh‐wave phase and group velocity measurements, and Bouguer gravity observations for the 3D shear‐wave speed. Using surface‐wave observations (3–250 s) and spatially smoothed receiver functions, our velocity models are robust, reliable, and rich in detail. The shear‐wave velocity models fit all three types of observations well. The resulting velocity model for the eastern U.S. shows thinner crust beneath New England, the east coast, and the Mississippi Embayment (ME). A relatively thicker crust was found beneath the stable North America craton. A relatively slower upper mantle was imaged beneath New England, the east coast, and western ME. A comparison of crust thickness derived from our model against four recent published models shows first‐order consistency. A relatively small upper mantle low‐speed region correlates with a published P‐wave analysis that has associated the anomaly with a 75 Ma kimberlite volcanic site in Kentucky. We also explored the relationship between the subsurface structure and seismicity in the eastern U.S. We found that earthquakes often locate near regions with seismic velocity variations, but not universally. Not all regions of significant subsurface wave speed changes are loci of seismicity. A weak correlation between upper mantle shear velocity and earthquake focal mechanism has been observed.

58 GEOSCIENCES↗
DOE OSTI2022

Spin-piston problem for a ferromagnetic thin film: Shock waves and solitons

In this report the unsteady, nonlinear magnetization dynamics induced by spin injection in an easy-plane ferromagnetic channel subject to an external magnetic field are studied analytically. Leveraging a dispersive hydrodynamic description, the Landau-Lifshitz equation is recast in terms of hydrodynamic-like variables for the magnetization's perpendicular component (spin density) and azimuthal phase gradient (fluid velocity). Spin injection acts as a moving piston that generates nonlinear, dynamical spin textures in the ferromagnetic channel with downstream quiescent spin density set by the external field. In contrast to the classical problem of a piston accelerating a compressible gas, here, variable spin injection and field lead to a rich variety of nonlinear wave phenomena from oscillatory spin shocks to solitons and rarefaction (expansion) waves. A full classification of solutions is provided using nonlinear wave modulation theory by identifying two key aspects of the fluid-like dynamics: subsonic/supersonic conditions and convex/nonconvex hydrodynamic flux. Familiar waveforms from the classical piston problem such as rarefaction waves and shocks manifest in their spin-based counterparts as smooth and highly oscillatory transitions, respectively, between two distinct magnetic states. The spin shock is an example of a dispersive shock wave, which arises in many physical systems. New features without a gas dynamics counterpart include composite wave complexes with “contact” spin shocks and rarefactions. Magnetic supersonic conditions lead to two pronounced piston edge behaviors including a stationary soliton and an oscillatory wave train. These coherent wave structures have physical implications for the generation of high frequency spin waves from pulsed injection and persistent, stable stationary and/or propagating solitons in the presence of magnetic damping. The analytical results are favorably compared with numerical simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Anharmonic Eigenvectors and Acoustic Phonon Disappearance in Quantum Paraelectric SrTiO 3

Pronounced anomalies in the SrTiO 3 dynamical structure factor, S(Q,E), including the disappearance of acoustic phonon branches at low temperatures, were uncovered with inelastic neutron scattering (INS) and simulations. The striking effect reflects anharmonic couplings between acoustic and optic phonons and the incipient ferroelectric instability near the quantum critical point. It is rationalized using a first-principles renormalized anharmonic phonon approach, pointing to nonlinear Ti-O hybridization causing unusual changes in real-space phonon eigenvectors, frequencies, group velocities, and scattering phase space. Lastly, our method is general and establishes how T dependences beyond the harmonic regime, assessed by INS mapping of large reciprocal-space volumes, provide real-space insights into anharmonic atomic dynamics near phase transitions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

CSAPR2 Optimized Convective Cell Tracking Data during TRACER

One of the challenges in analyzing convective cell properties is to observe the quick evolution of individual convective cells. While the operational radar data provide a volumetric data set to analyze radar observables of convective precipitation clouds, previous studies also suggested the quick evolution of cell life cycle that might not be captured by conventional radar volume scan strategies that take ~5-7 minutes. Aiming at enhancing our understanding of the links between convective cloud kinematic and microphysical processes as well as life cycles, the Tracking Aerosol Convection Interactions ExpeRiment (TRACER; Jensen et al. 2019) was conducted at Houston, Texas, in 2022. The TRACER campaign deployed the 2nd generation C-band Scanning ARM Precipitation Radar (CSAPR2), which performed frequent updates of range height indicator (RHI) and sector plan position indicator (PPI) scans to track individual convective cells every < 2 minutes, guided by a new cell tracking framework, Multisensor Agile Adaptive Sampling (MAAS; Kollias et al. 2020). This allows for capturing fast-evolving radar observables. We provide the processed CSAPR2 cell tracking data in CfRadial format collected during the TRACER field campaign from June to September 2022. The data files include processed radar variables: noise-masked reflectivity and differential reflectivity corrected for rain attenuation and systematic biases, noise-masked dealiased radial velocity, specific differential phase, locations of target cells (latitude, longitude, radar range), and radar-echo classification. Figure 1 provides an example of a 3D image of CSAPR2 reflectivity from the lowest PPI scan and an RHI scan after data processing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Synthesizing Stellar Populations in South Pole Telescope Galaxy Clusters. I. Ages of Quiescent Member Galaxies at 0.3 < $\textit{z}$ < 1.4

Using stellar population synthesis models to infer star formation histories (SFHs), we analyze photometry and spectroscopy of a large sample of quiescent galaxies that are members of Sunyaev–Zel'dovich (SZ)-selected galaxy clusters across a wide range of redshifts. We calculate stellar masses and mass-weighted ages for 837 quiescent cluster members at 0.3 < $\textit{z}$ < 1.4 using rest-frame optical spectra and the Python-based Prospector framework, from 61 clusters in the SPT-GMOS Spectroscopic Survey (0.3 < $\textit{z}$ < 0.9) and three clusters in the SPT Hi-z cluster sample (1.25 < $\textit{z}$ < 1.4). We analyze spectra of subpopulations divided into bins of redshift, stellar mass, cluster mass, and velocity-radius phase-space location, as well as by creating composite spectra of quiescent member galaxies. We find that quiescent galaxies in our data set sample a diversity of SFHs, with a median formation redshift (corresponding to the lookback time from the redshift of observation to when a galaxy forms 50% of its mass, $t_{50}$) of $\textit{z}$ = 2.8 ± 0.5, which is similar to or marginally higher than that of massive quiescent field and cluster galaxy studies. We also report median age–stellar mass relations for the full sample (age of the universe at $t_{50}$ (Gyr) = 2.52 (±0.04)–1.66 (±0.12) log$_{10}(M/10^{11}M_⊙$)) and recover downsizing trends across stellar mass; we find that massive galaxies in our cluster sample form on aggregate ~0.75 Gyr earlier than lower-mass galaxies. We also find marginally steeper age–mass relations at high redshifts, and report a bigger difference in formation redshifts across stellar mass for fixed environment, relative to formation redshifts across environment for fixed stellar mass.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Dephasingless Laser Wakefield Acceleration

Laser wakefield accelerators (LWFAs) produce significantly high gradients enabling compact accelerators and radiation sources, but face design limitations, such as dephasing, occurring when trapped electrons outrun the accelerating phase of the wakefield. In this report we combine spherical aberration with a novel cylindrically symmetric echelon optic to spatiotemporally structure an ultra-short, high-intensity laser pulse that can overcome dephasing by propagating at any velocity over any distance. The ponderomotive force of the spatiotemporally shaped pulse can drive a wakefield with a phase velocity equal to the speed of light in vacuum, preventing trapped electrons from outrunning the wake. Simulations in the linear regime and scaling laws in the bubble regime illustrate that this dephasingless LWFA can accelerate electrons to high energies in much shorter distances than a traditional LWFA—a single 4.5 m stage can accelerate electrons to TeV energies without the need for guiding structures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Imaging the Shallow Structure of the Yucca Flat at the Source Physics Experiment Phase II Site with Horizontal-to-Vertical Spectral Ratio Inversion and a Large- N Seismic Array

The Source Physics Experiment (SPE) is a series of chemical explosions at the Nevada National Security Site (NNSS) with the goal of understanding seismic-wave generation and propagation of underground explosions. To understand explosion source physics, accurate geophysical models of the SPE site are needed. Here, we utilize a large-N seismic array deployed at the SPE phase II site to generate a shallow subsurface model of shear-wave velocity. The deployment consists of 500 geophones and covers an area of, approximately, 2.5 × 2 km. The array is located in the Yucca Flat in the northeast corner of the NNSS, Nye County, Nevada. Using ambient-noise recordings throughout the large-N seismic array, we calculate horizontal-to-vertical spectral ratios (HVSRs) across the array. We obtain 2D seismic images of shear-wave velocities across the SPE phase II site for the shallow structure of the basin. In this work the results clearly image two significant seismic impedance interfaces at ~150–500 and ~350–600 m depth. The shallower interface relates to the contrast between Quaternary alluvium and Tertiary volcanic rocks. The deeper interface relates to the contrast between Tertiary volcanic rocks and the Paleozoic bedrock. The 2D subsurface models support and extend previous understanding of the structure of the SPE phase II site. This study shows that the HVSR method in conjunction with a large-N seismic array is a quick and effective method for investigating shallow structures.

58 GEOSCIENCES↗
DOE OSTI2020

A conservative phase-space moving-grid strategy for a 1D-2V Vlasov–Fokker–Planck Solver

In this work, we develop a conservative configuration- and velocity-space (i.e., phase-space) moving-grid strategy for the Vlasov–Fokker–Planck (VFP) equation in a planar geometry. The velocity-space grid is normalized and shifted in terms of the thermal speed and the bulk-fluid velocity, respectively. The configuration-space grid is moved according to a mesh-motion-partial-differential equation (MMPDE), which equidistributes a monitor function that is inversely proportional to the gradient-length scales of the macroscopic plasma quantities. The resulting inertial terms in the transformed VFP equations are discretized to ensure the discrete conservation of mass, momentum, and energy. To satisfy the discrete conservation theorems in the presence of phase-space mesh motion, we employ the method of discrete nonlinear constraints – explored in previous studies – but the underlying symmetries are determined in a much more efficient manner than before. The conservative grid-adaptivity strategy provides an efficient scheme that resolves important physical structures in the phase-space while controlling the computational complexity at all times. We demonstrate the favorable features of the proposed algorithm through a set of test cases of increasing complexity. The problems test independent components of the algorithms, as well as the integrated capability on settings relevant to inertial confinement fusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A numerical study of particle jetting in a dense particle bed driven by an air-blast

In this work, the particle jetting behavior in a blast-driven dense particle bed is studied at early times. Four-way coupled Euler–Lagrange simulations are performed using a high-order discontinuous Galerkin spectral element solver coupled with a high-order Lagrangian particle solver, wherein the inter-particle collisions are resolved using a discrete element method collision model. Following the experiments of Rodriguez et al. [“Formation of particle jetting in a cylindrical shock tube,” Shock Waves 23(6), 619–634 (2013)] and the simulations of Osnes et al. [“Numerical simulation of particle jet formation induced by shock wave acceleration in a Hele-Shaw cell,” Shock Waves 28(3), 451–461 (2018)], the simulations are performed in a quasi-two-dimensional cylindrical geometry (Hele-Shaw cell). Parametric studies are carried out to assess the impact of the coefficient of restitution and the strength of the incident shock on the particle jetting behavior. The deposition of vorticity through a multiphase (gas–particle) analog of Richtmyer–Meshkov instability is observed to play a crucial role in channeling the particles into well-defined jets at the outer edge of the particle bed. This is confirmed by the presence of vortex pairs around the outer jets. Furthermore, the effect of the relaxation of the relative velocity between the two phases on the vorticity generation is explored by analyzing the correlation between the radial velocity of particles and the radial velocity of the gas at the particle location.

42 ENGINEERING↗
DOE OSTI2020

Extended anharmonic collapse of phonon dispersions in SnS and SnSe

The lattice dynamics and high-temperature structural transition in SnS and SnSe are investigated via inelastic neutron scattering, high-resolution Raman spectroscopy and anharmonic first-principles simulations. We uncover a spectacular, extreme softening and reconstruction of an entire manifold of low-energy acoustic and optic branches across a structural transition, reflecting strong directionality in bonding strength and anharmonicity. Further, our results solve a prior controversy by revealing the soft-mode mechanism of the phase transition that impacts thermal transport and thermoelectric efficiency. Our simulations of anharmonic phonon renormalization go beyond low-order perturbation theory and capture these striking effects, showing that the large phonon shifts directly affect the thermal conductivity by altering both the phonon scattering phase space and the group velocities. These results provide a detailed microscopic understanding of phase stability and thermal transport in technologically important materials, providing further insights on ways to control phonon propagation in thermoelectrics, photovoltaics, and other materials requiring thermal management.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Velocity structure functions in multiphase turbulence: interpreting kinematics of Hα filaments in cool-core clusters

ABSTRACT The central regions of cool-core galaxy clusters harbour multiphase gas, with gas temperatures ranging from $10$ to $10^7\, \mathrm{K}$. Feedback from active galactic nuclei jets prevents the gas from undergoing a catastrophic cooling flow. However, the exact mechanism of this feedback energy input is unknown, mainly due to the lack of velocity measurements of the hot-phase gas. However, recent observations have measured the velocity structure functions (VSFs) of the cooler molecular (${\sim} 10\, \mathrm{K}$) and Hα filaments (${\sim} 10^4\, \mathrm{K}$) and used them to indirectly estimate the motions of the hot phase. In the first part of this study, we conduct high-resolution (3843–15363 resolution elements) simulations of homogeneous isotropic subsonic turbulence, without radiative cooling. We analyse the second-order velocity structure functions (VSF2) in these simulations and study the effects of varying spatial resolution, the introduction of magnetic fields, and the effect of projection along the line of sight (LOS) on it. In the second part of the study, we analyse high-resolution (7683 resolution elements) idealized simulations of multiphase turbulence in the intracluster medium from the companion study Mohapatra et al. We compare the VSF2 for both the hot ($T\sim 10^7\, \mathrm{K}$) and cold ($T\sim 10^4\, \mathrm{K}$) phases and find that their amplitude depends on the density contrast between the phases. They have similar scaling with separation, but introducing magnetic fields steepens the VSF2 of only the cold phase. We also find that projection along the LOS steepens the VSF2 for the hot phase and mostly flattens it for the cold phase.

Mohapatra, Rajsekhar (ORCID:0000000216007552)↗
DOE OSTI2022

Seismic Waveform Inversion Capability on Resource-Constrained Edge Devices

Seismic full wave inversion (FWI) is a widely used non-linear seismic imaging method used to reconstruct subsurface velocity images, however it is time consuming, has high computational cost and depend heavily on human interaction. Recently, deep learning has accelerated it’s use in several data-driven techniques, however most deep learning techniques suffer from overfitting and stability issues. In this work, we propose an edge computing-based data-driven inversion technique based on supervised deep convolutional neural network to accurately reconstruct the subsurface velocities. Deep learning based data-driven technique depends mostly on bulk data training. In this work, we train our deep convolutional neural network (DCN) (UNet and InversionNet) on the raw seismic data and their corresponding velocity models during the training phase to learn the non-linear mapping between the seismic data and velocity models. The trained network is then used to estimate the velocity models from new input seismic data during the prediction phase. The prediction phase is performed on a resource-constrained edge device such as Raspberry Pi. Raspberry Pi provides real-time and on-device computational power to execute the inference process. In addition, we demonstrate robustness of our models to perform inversion in the presence on noise by performing both noise-aware and no-noise training and feeding the resulting trained models with noise at different signal-to-noise (SNR) ratio values. We make great efforts to achieve very feasible inference times on the Raspberry Pi for both models. Specifically, the inference times per prediction for UNet and InversionNet models on Raspberry Pi were 22 and 4 s respectively whilst inference times for both models on the GPU were 2 and 18 s which are very comparable. Finally, we have designed a user-friendly interactive graphical user interface (GUI) to automate the model execution and inversion process on the Raspberry Pi.

Manu, Daniel (ORCID:0000000154982677)↗