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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Thermal management of a high temperature sodium sulphur battery stack

The sodium sulfur battery is an advanced secondary battery with high potential for grid-level storage due to their high energy density, low cost of the reactants, and high open-circuit voltage. However, as the operating temperature of the battery is high (about 300 °C), effective thermal management is required to prevent thermal runaway under high current density operation. Here, to develop efficient thermal management strategies, a detailed, thermo-electrochemical, non-isothermal, distributed dynamic model of a sodium-sulfur battery stack is developed. Models of three thermal management strategies are developed and analyzed in this work: active cooling, passive cooling, and hybrid cooling. The active cooling strategy uses air as the cooling medium whereas, the passive cooling strategy uses a phase change material. In the hybrid cooling strategy, both active and passive cooling strategies are used. Due to the high operating temperature of the sodium sulfur batteries, the rejected heat can be utilized effectively and consequently a high variability in the heat rejection rate may not be acceptable. Performance of these strategies is analyzed under high current density operation by evaluating the variability in the cell temperature during charge/discharge cycles, temperature difference between the cells based on their locations, and variability in the heat rejection rate. The phase change-based hybrid thermal management strategy with an embedded controller maintained the temperature variation within ±/- 2.8 °C. Considering the tradeoff between the capital cost for PCM, fan power requirement, and variability in the heat rejection rate, the optimal quantity of the PCM and the maximum air velocity for the maximum current density of 260 mA/cm 2 (0.55 C-rate) were found be 0.74 g/Wh and 3 m/s, respectively.

25 ENERGY STORAGE↗
DOE OSTI2021

Directional clogging and phase separation for disk flow through periodic and diluted obstacle arrays

We model collective disk flow though a square array of obstacles as the flow direction is changed relative to the symmetry directions of the array. At lower disk densities there is no clogging for any driving direction, but as the disk density increases, the average disk velocity decreases and develops a drive angle dependence. For certain driving angles, the flow is reduced or drops to zero when the system forms a heterogeneous clogged state consisting of high density clogged regions coexisting with empty regions. The clogged states are fragile and can be unclogged by changing the driving angle. For large obstacle sizes, we find a uniform clogged state that is distinct from the collective clogging regime. Within the clogged phases, depinning transitions can occur as a function of increasing driving force, with intermittent motion appearing just above the depinning threshold. The clogging is robust against the random removal or dilution of the obstacle sites, and the disks are able to form system-spanning clogged clusters even under increasing dilution. If the dilution becomes too large, however, the clogging behavior is lost.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Directed energy deposition of radially functionally graded nickel based superalloys IN625 and IN738: Microstructure and mechanical properties

This study investigates the microstructure and mechanical properties of radially functionally graded (RFG) materials fabricated from nickel-based superalloys IN625 and IN738 using directed energy deposition (DED). RFG cylindrical deposits with radial gradients of IN625-IN738 and IN738-IN625 were successfully fabricated and confirmed to be crack-free through tailored composition grading, scan strategy and process parameters that minimized thermal gradients and solidification velocities. Microstructural analysis revealed distinct variations across the gradient, characterized by differences in dendrite arm spacing, microsegregation, and the formation of distinct γ' precipitates and carbide/boride phases (MC, M 6 C, M 2 C, and MB 2 ) in as-deposited (AD) conditions, which evolved into M 23 C 6 and M 2 B upon subsequent heat treatment (HT). EBSD showed predominant columnar grains with a [100] texture, while microhardness profiles demonstrated consistent variation across the gradients. Mechanical testing highlighted the synergy of the graded structure, combining IN738's strength with IN625's ductility, resulting in enhanced properties such as a yield strength of 1131 MPa, ultimate tensile strength of 1382 MPa, and elongation of 12 % for HT IN625-IN738 grading. Numerical simulations employing CALPHAD thermodynamic calculations and finite element thermal modeling are used to establish phase evolution and microstructural feature correlation with thermal profiles, aiding in the optimization of process-structure-property relationship.

Additive manufacturing↗
DOE OSTI2021

Modeling Subsurface Explosions Recorded on a Distributed Fiber Optic Sensor

Fiber optic distributed acoustic sensors (DAS) are becoming a widely used tool for seismic sensing. Here we examine recordings of two subsurface chemical explosions, DAG-1 and DAG-3, each of which was about one metric ton (TNT equivalent), that were recorded from a helical fiber installed in two boreholes 80 m away from the source location. Several clear phases including the initial P wave, a weak S wave, and a surface reflected P wave are observed on the helical DAS data. We estimate a velocity model using arrival times measured from the fiber. The DAS waveform data were compared with colocated accelerometers at specific depths in both frequency and time domains. Furthermore, the spectra of the DAS data matched spectra estimated from the accelerometer records. Comparisons of observed waveform shape between the accelerometer records and the fiber measurements (strain-rate) show reasonable agreement except for the data near the event depth. The DAS data and the accelerometer agreed in relative amplitudes but we had difficulties in matching absolute amplitudes, possibly due to errors in metadata. Synthetic strain-rate waveforms were calculated using a 2D wavenumber algorithm and matched the waveform shape and relative amplitudes. In general, DAS is effective at recording strong ground motions at high spatial density. Comparison of the synthetic seismograms with observed data indicate that the waveforms are not consistent with a pure isotropic explosion source and that the observed S waves originate from very near the source region.

58 GEOSCIENCES↗
DOE OSTI2021

Polymorphism of gold under laser-based ramp compression to 690 GPa

The high-pressure behavior and polymorphism of gold has attracted strong recent interest from theory as well as from both static and dynamic experiments. In this study, gold was dynamically ramp-compressed to 690 GPa via laser ablation using the Omega Laser Facility (U. Rochester) in order to investigate the stability regions of this noble metal's polymorphs. Here, stresses were determined from measurements of surface or interface velocities by laser interferometry. In situ X-ray diffraction was carried out to constrain the structure and density. The ambient face-centered-cubic (fcc) phase was observed at 162 GPa, followed by a mixed phase region consisting of the fcc and body-centeredcubic (bcc). The presence of only the bcc phase was observed between 377 and 690 GPa. Our results are consistent with recent shock compression experiments and demonstrate that gold transforms to the bcc phase under conditions of both ramp and shock loading.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Nonlinear dynamics, avalanches, and noise for driven Wigner crystals

Here we consider the driven dynamics of Wigner crystals interacting with random disorder. Using numerical simulations, we find a rich variety of transport phenomena as a function of charge density, drive, and pinning strength. For weak pinning, the system forms a defect-free crystal that depins elastically. When the pinning is stronger, a pinned glass phase appears that depins into a filamentary flow state, transitions at higher drives into a disordered flow phase, and finally forms a moving smectic. Within the filamentary flow phase, the conduction curves can show switching dynamics as well as negative differential conductivity in which switching events cause some flow channels to be blocked. The velocity noise in the filamentary flow regime exhibits narrow band characteristics due to the one-dimensional nature of the motion, while the moving smectic has narrow band velocity noise with a washboard frequency. In the disordered flow state, the noise power reaches a peak value and the noise has a 1/$\mathcal{f}$ character. Our transport results are consistent with recent experimental transport studies in systems where Wigner crystal states are believed to be occurring. Below the conduction threshold, we find that avalanches with a power-law size distribution appear when there are sudden local rearrangements of charges from one pinned configuration to another.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

The Magnitude and Waveform of Shock Waves Induced by X-ray Lasers in Water

The high energy densities deposited in materials by focused X-ray laser pulses generate shock waves which travel away from the irradiated region, and can generate complex wave patterns or induce phase changes. We determined the time-pressure histories of shocks induced by X-ray laser pulses in liquid water microdrops, by measuring the surface velocity of the microdrops from images recorded during the reflection of the shock at the surface. Measurements were made with ~30 µm diameter droplets using 10 keV X-rays, for X-ray pulse energies that deposited linear energy densities from 3.5 to 120 mJ/m; measurements were also made with ~60 µm diameter drops for a narrower energy range. At a distance of 15 µm from the X-ray beam, the peak shock pressures ranged from 44 to 472 MPa, and the corresponding time-pressure histories of the shocks had a fast quasi-exponential decay with positive pressure durations estimated to range from 2 to 5 ns. Knowledge of the amplitude and waveform of the shock waves enables accurate modeling of shock propagation and experiment designs that either maximize or minimize the effect of shocks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Novel challenges in tracking self-interacting dark matter subhalos

Cosmological N-body simulations are among the primary tools for studying structure formation in the Universe. Analyses of these simulations critically depend on accurately identifying and tracking dark matter subhalos over time. In recent years, several new algorithms have been developed to improve the accuracy and consistency of subhalo tracking in cold dark matter (CDM) simulations. These algorithms should be revisited in the context of new physics beyond gravity, which can modify the evolution and final properties of subhalo populations. In this work, we apply the particle-tracking-based subhalo finder Symfind to velocity-dependent self-interacting dark matter (SIDM) simulations with large cross section amplitudes to assess the performance of particle-tracking methods beyond the CDM paradigm. We find that the core-particle-tracking technique, which is key to the success of these algorithms in CDM, does not always yield accurate results in SIDM. In particular, the interplay between dark matter self-interactions and tidal stripping can cause the diffusion of core particles to larger radii, leading particle-tracking-based algorithms to prematurely lose track of SIDM subhalos. For massive core-expansion subhalos and core-collapse subhalos that experience close or repeated pericentric passages, a significant fraction of core particles can be lost, and particle-tracking-based finders such as Symfind offer no clear advantage over traditional methods that rely on identifying phase-space overdensities. On the other hand, for subhalos with large pericentric distances or fewer, more distant passages, Symfind tends to outperform. These differences depend sensitively on the cross section amplitude and turnover velocity of the underlying SIDM model. We therefore recommend a hybrid approach that leverages the strengths of both techniques to produce complete and robust catalogs of core-expansion and core-collapse SIDM subhalos.

cosmological simulations↗
DOE OSTI2023

The Variable Continuous Bimaterial Interface in the San Jacinto Fault Zone Revealed by Dense Seismic Array Analysis of Fault Zone Head Waves

Key factors controlling earthquake ruptures include fault geometry, continuity, and seismic velocity structure around the fault. We present a novel tool that better informs deep bimaterial fault geometry embedded in distributed damage and seismicity, associated velocity contrasts across the fault, and their correlations with surface complexities. The method employs fault zone head and direct body waves and is applied to recordings from five spatiotemporally different seismic arrays along the complex San Jacinto fault zone (SJFZ) in southern California. We detect and distinguish these signals based on instantaneous phase coherence and relative energy in a cascading manner from one scale array to another. The analysis reveals a >70-km long continuous bimaterial interface within the SJFZ with several deep northeast dipping fault segments. The northern SJFZ, for instance, locates ~7 km northeast of its surface expression at 18-km depth. P-wave velocity contrasts range from near 0% to >15%, consistent with other bimaterial faults, and differ by a few % depending on fault-array azimuth, implying directional-dependent velocity contrasts. S-wave head waves and velocity contrasts are also imaged for the first time at the southern SJFZ, averaging to 2.9% in agreement with tomography results. The imaged geometry and continuity suggest the SJFZ initiated along remnant tectonic structures and translates to a rupture potential of M > 7.2, i.e., the sizes of its largest paleo-earthquakes. The P and S contrasts, and their ratios, have important implications for earthquake rupture speed, mode, directivity, and frictional heating along the SJFZ and other major faults globally.

58 GEOSCIENCES↗
DOE OSTI2026

Development and Application of High-Fidelity Models for Heterogeneous CO2 Frost Formation

Carbon America has developed a cryogenic carbon capture technology ("FrostCC") that separates CO2 from point source emissions by solidifying it at cold temperatures through preferential desublimation. Cooling is achieved through a series of interlinked compression, heat exchange, and expansion operations. In the current system, frosting of CO2 happens in heat exchangers, followed by CO2 recovery in a separate extraction step. In this work, multiphysics computational fluid dynamics (CFD) models are developed and validated for compressible and low Mach flows to simulate the formation of solid CO2 in flue gas flowing in a heat exchanger geometry. The models track the mass transfer rate of CO2 from gas phase to solid phase, heat released from desublimation, and the evolution of the solid CO2 layer. Simulations are used to answer scientific questions related to the angle of heat exchanger pipes, where buoyancy effects from flow velocity and pipe orientation influence CO2 frosting. Results show that upwardly angled pipes produce notably different flow structures compared to horizontal or vertical configurations, and that carbon capture efficiency correlates with buoyancy effects for pipe angles within plus or minus 23 degrees of horizontal.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Flow pattern and void fraction characterization in nitrogen/water flows through a diamond-type triply periodic minimal surface lattice

Here, this study presents, to the authors knowledge, the first experimental investigation on the void fraction and flow patterns in two-phase flows through a diamond-type Triply Periodic Minimal Surface (TPMS) lattice. An additively manufactured TPMS structure was tested under upward co-current flow of a water/nitrogen mixture. Superficial velocities varied for a total of 42 test conditions (gas: 0.01–2.4 m/s; liquid: 0.01–2.4 m/s; mass flux: 20–2370 kg/m 2 ·s). High-speed video and X-ray imaging enabled time-averaged void fraction measurements and identified six distinct flow regimes which were used to develop a flow pattern map. Comparison of the void fraction data with correlations from literature demonstrated the Rouhani and Axelsson (1970) [49] correlation modified by Steiner (1993) [52] provided the best agreement, which was improved with empirically fit coefficients. This approach predicted the void fraction with errors < ±20% for 64% of the data, with a mean absolute percent deviation (MAPD) of 23%. The measured frictional pressure drop was compared to correlations from literature which captured the observed trends but did not provide good accuracy. The best agreement was found after optimizing the empirical coefficients of the Muller-Steinhagen & Heck (1986) [56] correlation; this approach captured 33% of the data within ±20% with a MAPD of 46% over the full range and captured 71% of the data within ±20% a MAPD of 13.3% at mass fluxes >1100 kg/m 2 -s. These results provide foundational insight into TPMS two-phase flow behavior and inform modeling and design of advanced heat exchange components incorporating TPMS geometries.

42 - ENGINEERING↗
DOE OSTI2024

Blade Designs For Improved Multi-Phase Performance In sCO 2 Compressors; Part II - Optical Diagnostics In sCO 2 And Experimental Evaluation With Particle Image Velocimetry

This paper presents the second part of a study in which the leading-edge and suction surface of a compressor blade was modified to delay onset of phase change for sCO 2 compressors operating near the critical point. Using a first-of-its-kind apparatus for the measurement of sCO 2 flow fields, Particle Image Velocimetry (PIV) is used for local flow field measurements of two compressor blade geometries: the modified “biased-wedge,” and a conventional constant thickness blade. Utilizing the developed hardware, the feasibility of a simple, laser-based diagnostic for qualitatively measuring liquid phase regions, is also presented. The design of the optical diagnostics rig, a discussion of numerous challenges, and necessary considerations involved in performing optical-based measurements like PIV, in sCO 2 , are discussed. Velocity field measurements for the modified compressor profile show a much lower suction peak compared to a conventional blade. Furthermore, these results validate numerical results at the tested conditions, where the suction side profile of the biased wedge works to minimize the local pressure gradient.

14 SOLAR ENERGY↗
DOE OSTI2025

Discovery of Powerful Multivelocity Ultrafast Outflows in the Starburst Merger Galaxy IRAS 05189–2524 with XRISM

We observed the X-ray-bright ultraluminous infrared galaxy IRAS 05189−2524 with XRISM during its performance verification phase. The unprecedented energy resolution of the onboard X-ray microcalorimeter revealed complex spectral features at ∼7–9 keV, which can be interpreted as blueshifted Fe XXV/XXVI absorption lines with various velocity dispersions, originating from ultrafast outflow (UFO) components with multiple bulk velocities of ∼0.076c, ∼0.101c, and ∼0.143c. In addition, a broad Fe–K emission line was detected around ∼7 keV, forming a P Cygni profile together with the absorption lines. The onboard X-ray CCD camera revealed a 0.4–12 keV broadband spectrum characterized by a neutrally absorbed power-law continuum with a photon index of ∼2.3 and intrinsic flare-like variability on timescales of ∼10 ks, both of which are likely associated with near-Eddington accretion. We also found potential variability of the UFO parameters on a timescale of ∼140 ks. Using these properties, we propose new constraints on the outflow structure and suggest the presence of multiple outflowing regions on scales of about tens to 100 Schwarzschild radii, located within roughly 2000 Schwarzschild radii. Since both the estimated momentum and energy outflow rates of the UFOs exceed those of galactic molecular outflows, our results indicate that powerful, multivelocity UFOs are already well developed during a short-lived evolutionary phase following a major galaxy merger, characterized by intense starburst activity and likely preceding the quasar phase. This system is expected to evolve into a quasar, sustaining strong UFO activity and suppressing star formation in the host galaxy.

Astronomy and AstroPhysics↗
DOE OSTI2021

Quantifying particle-scale 3D granular dynamics during rapid compaction from time-resolved in situ 2D x-ray images

Understanding the particle-scale dynamics of granular materials during rapid compaction and flow is of fundamental importance for manufacturing, planetary science, geology, and defense applications. Time-resolved 2D radiography and static 3D x-ray tomography are powerful in situ tools for studying particle-scale dynamics but provide detail only in 2D or with significant time-scale limitations, respectively. Here, we introduce a new method that uses 2D in situ x-ray imaging for determining time-resolved 3D particle-scale dynamics in rapidly compressed granular materials. The method employs initial particle packing structures obtained from x-ray tomography, a 2D x-ray image generation algorithm, and an optimization algorithm. We first describe and validate the method using finite element simulations. Furthermore, we then apply the technique to x-ray phase-contrast images obtained during rapid compaction of granular materials with varying particle morphology and sample thickness. The depth-resolved particle-scale dynamics reveal complex velocity and porosity fields evolving heterogeneously along and perpendicular to the compaction direction. We characterize these features, their fluctuations near the compaction front, and the compaction front thickness. Our technique can be applied to understanding granular dynamics during rapid compaction events, and rearrangements during slower, but non-quasi-static, flows.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Designing Monte Carlo Simulation and an Optimal Machine Learning to Optimize and Model Space Missions

This paper investigates applying artificial intelligence (AI) algorithms to attitude control system of satellites to optimally tune the controller using high performance computing. This methodology is applied to the Virtual Telescope for X-ray Observation mission, which is a precise formation of two separate spacecraft observing multiple objects in the space in the X-ray domain. The mission is divided into phases based on the instrumentation and the mission goal. To reach an stable precise formation robust to stochastic slew and slew rate (i.e., Euler angles and angular velocities) in a minimal constrained time T , consumed energy of the attitude control system, denoted as E, and root-mean-square state error of attitude control system, denoted as e, are minimized. Monte-Carlo simulation is used for the sensitivity analysis of optimization and designing a controller. Deep neural networks (DNN), Gaussian processes (GP), and support vector regression (SVR) learn this optimization as a surrogate model, while their hyperparameters are optimized in a novel approach. THETA supercomputer at Argonne Leadership Computing Facility (ALCF) is used for optimizing the hyperparameters of DNN. The surrogate model meets the requirements of the mission, and it shows a better performance over the optimization and Monte-Carlo. The optimal DNN can satisfy the mission requirements e and T while reducing E for 90% compared to the other given methods.

42 ENGINEERING↗
DOE OSTI2020

KR Persei, a Mid-F Eclipsing Binary with a One-day Period

KR Per is a partially eclipsing binary with an orbital period of 0.9960798 days, very close to one sidereal day, making it difficult to obtain extensive phase coverage in a reasonable amount of time. We used the Wilson–Devinney program to determine its orbital elements and stellar absolute dimensions from recently acquired radial velocities and differential BVRI observations that were supplemented with previous differential UBV measurements and published times of minima. The two components are each F5 V stars with masses of 1.466±0.015 M{sub ⊙} and 1.458±0.015 M{sub ⊙}. The radii are 1.855±0.021 R{sub ⊙} and 1.824±0.022 R{sub ⊙}. The orbital period of the eclipsing system is variable and more times of minima observations are needed. A comparison with evolutionary tracks indicates that the system has an age of 2.1 ± 0.1 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Extreme Heating of Minor Ions in Imbalanced Solar-wind Turbulence

Minor ions in the solar corona are heated to extreme temperatures, far in excess of those of the electrons and protons that comprise the bulk of the plasma. These highly nonthermal distributions make minor ions sensitive probes of the collisionless processes that heat the corona and power the solar wind. The recent discovery of the "helicity barrier" offers a mechanism in which imbalanced Alfvénic turbulence in low-β plasmas preferentially heats protons over electrons, generating high-frequency, proton-cyclotron-resonant fluctuations. We use the hybrid-kinetic particle-in-cell code Pegasus++ to drive imbalanced Alfvénic turbulence in a 3D low-β plasma with additional passive ion species, He 2+ and O 5+ . A helicity barrier naturally develops, followed by clear phase-space signatures of oblique proton-cyclotron-wave heating and Landau-resonant heating from the imbalanced Alfvénic fluctuations. The former results in characteristically arced ion velocity distribution functions, whose non-bi-Maxwellian features are shown by linear ALPS calculations to be critical to the heating process. Additional features include a steep transition-range electromagnetic spectrum, proton-cyclotron waves propagating in the direction of the imbalance, significantly enhanced proton-to-electron heating ratios, ion temperatures that are considerably more perpendicular with respect to magnetic field, and extreme heating of heavier species in a manner consistent with mass scalings inferred from spacecraft measurements. None of these features are realized in an otherwise equivalent simulation of balanced turbulence. If seen simultaneously in the fast solar wind, these signatures of the helicity barrier would testify to the necessity of incorporating turbulence imbalance in a complete theory for the evolution of the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Utilizing TiO 2 amorphous precursors for polymorph selection: An in situ TEM study of phase formation and kinetics

Selective synthesis of metastable polymorphs requires a fundamental understanding of the complex energy landscapes in which these phases form. Recently, the development of in situ high temperature and controlled atmosphere transmission electron microscopy has enabled the direct observation of nucleation, growth, and phase transformations with near atomic resolution. In this work, we directly observe the crystallization behavior of amorphous TiO 2 thin films grown under different pulsed laser deposition conditions and quantify the mechanisms behind metastable crystalline polymorph stabilization. Films deposited at 10 mTorr chamber oxygen pressure crystallize into nanocrystalline Anatase at 325 degrees C, whereas films deposited at 2 mTorr crystallize into significantly larger needle-like grains of Brookite and Anatase at 270 degrees C. Increasing film deposition rate by a factor of 4 results in a 10x increase in the crystalline growth front velocity as well as a decrease in crystallization temperature from 270 degrees C to 225 degrees C. Engineering the amorphous precursor state through deposition conditions therefore provides routes to microstructure control and the accessibility of higher energy metastable phases.

36 MATERIALS SCIENCE↗