Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics simulations”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Shock physics in warm dense matter: A quantum hydrodynamics perspective

Abstract Warm dense matter (WDM), an exotic, highly compressed state of matter between solid and plasma phases, is of high current interest, in particular for astrophysics and inertial confinement fusion. For the latter, in particular the propagation of compression shocks is crucial. The main unknown in the shock propagation in WDM is the behaviour of the electrons since they are governed by correlations, quantum and spin effects that need to be accounted for simultaneously. Here we describe the shock dynamics of the warm dense electron gas using a quantum hydrodynamic model. From the numerical hydrodynamic simulations, we observe that the quantum Bohm pressure induces shear force that weakens the formation and strength of the shock. In addition, the Bohm pressure induces an electron density response that takes the form of oscillations. This is confirmed by the theoretical analysis of the early stage of the shock formation. Our theoretical and numerical analyses allow us to identify characteristic dimensionless shock propagation parameters at which the effect of the Bohm force is important.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Double cylinder implosion experiments at the National Ignition Facility

Here, cylindrical implosion experiments are used to directly measure instability growth in a convergent geometry, providing a wealth of data for model validation. Double cylinders are a natural extension of the platform and enable measurements at a classically unstable interface, the outer surface of the inner cylinder, which experiences no ablative stabilization from the laser drive. However, the utility of this platform relies upon maintaining adequate axial uniformity of the inner cylinder during the implosion. Although previous smaller-scale double cylinder experiments exhibited acceptable levels of axial uniformity, radiation-hydrodynamics simulations of larger-scale double cylinders predict more axial non-uniformity induced by the impedance mismatch as the shock wraps around the axial ends of the inner cylinder. A mechanism to reduce axial non-uniformity in these larger double cylinder implosions is presented, and preliminary experimental data confirms the efficacy of the selected mitigation approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Hydro-scaling of direct-drive cylindrical implosions at the OMEGA and the National Ignition Facility

Deceleration-phase Rayleigh–Taylor instability (RTI) growth during inertial confinement fusion capsule implosions significantly affects the performance as it mixes cold ablator material into the fuel. Precise measurements of such instability growth are essential for both validating the existing simulation codes and improving our predictive capability. RTI measurements on the inner surface of a spherical shell are limited and are often inferred indirectly at limited convergence. In contrast, cylindrical implosions allow for direct diagnostic access to the converging interface by imaging down the cylinder axis while retaining the effects of convergence. We have performed direct-drive cylindrical implosion experiments at both the OMEGA and the NIF laser facilities using scaled targets. RTI growth is demonstrated to be scale-invariant between the cylindrical targets at OMEGA and similar targets at the NIF, which are scaled up by a factor of three in the radial dimension. Single-mode (m = 0) instability growth factors of ~14 are measured at a convergence ratio (CR) ~ 2.5 with nearly identical mode growth at both scales. The measurements are in agreement with xRAGE radiation-hydrodynamics simulations. In addition, we have developed the Bayesian-inference-engine method to account for the variations in the target alignment, magnification, and the parallax effect in the measurement, allowing a more precise comparison between the experimental data and the simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Material flow behavior and microstructural refinement of AA6061 alloy during friction extrusion

In this work we used friction extrusion (FE), a solid phase processing technique, to produce dense, fully consolidated 5?mm rods of aluminum alloy 6061 (AA6061). The combination of large shear stresses and high temperatures at the tool-billet interface during extrusion produced equiaxed, dynamically recrystallized grains and finely distributed precipitates. Texture and microstructure evolution during extrusion was investigated in detail using scanning electron microscopy, electron backscatter diffraction, and transmission electron microscopy. Smoothed particle hydrodynamics simulation was performed to study complex material flow during extrusion. Simulation results suggest spiral material flow during extrusion which corroborated experimental results. Advantages of friction extruded microstructure over conventionally extruded counterparts are also explored using flash annealing for solution treatment followed by artificial aging. Mechanical properties of the as-friction extruded, and the artificially aged specimens were evaluated using tensile testing and compared with conventional extruded material properties.

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

Single-shot in-line x-ray phase-contrast imaging of void-shockwave interactions in fusion energy materials

Recent breakthroughs in nuclear fusion, specifically the report of reactions exceeding scientific breakeven at the National Ignition Facility (NIF), highlight the potential of inertial fusion energy (IFE) as a sustainable and virtually limitless energy source. However, further progress in IFE requires characterization of defects in ablator materials and how they affect fuel capsule compression. Voids within the ablator can degrade energy yield, but their impact on the density distribution has primarily been studied through simulations, with limited high-resolution experimental validation. To address this, we used the x-ray free-electron laser (XFEL) at the matter in extreme conditions (MECs) instrument at the Linac coherent light source (LCLS) to capture 2D x-ray phase-contrast (XPC) images of a void-bearing sample with a composition similar to inertial confinement fusion (ICF) ablators. By driving a compressive shockwave through the sample using MEC's long-pulse laser system, we analyzed how voids influence shockwave propagation and density distribution during compression. To quantify this impact, we extracted phase information using two phase retrieval algorithms. First, we applied the contrast transfer function (CTF) method, paired with Tikhonov regularization and a fast optimization approach to generate an initial phase estimate. We then refined the result using a projected gradient descent (PGD) method that works directly with the sample's refractive index. Comparing these results with radiation adaptive grid Eulerian (xRAGE) radiation hydrodynamic simulations enables identification of model validation needs or improvements. By calculating phase maps in situ, it becomes possible to reconstruct areal density maps, improving understanding of laser-capsule interactions and advancing IFE research.

Hodge, D. S. [Colorado State Univ., Fort Collins, ↗
DOE OSTI2021

Dark energy survey year 1 results: Constraining baryonic physics in the Universe

ABSTRACT Measurements of large-scale structure are interpreted using theoretical predictions for the matter distribution, including potential impacts of baryonic physics. We constrain the feedback strength of baryons jointly with cosmology using weak lensing and galaxy clustering observables (3 × 2pt) of Dark Energy Survey (DES) Year 1 data in combination with external information from baryon acoustic oscillations (BAO) and Planck cosmic microwave background polarization. Our baryon modelling is informed by a set of hydrodynamical simulations that span a variety of baryon scenarios; we span this space via a Principal Component (PC) analysis of the summary statistics extracted from these simulations. We show that at the level of DES Y1 constraining power, one PC is sufficient to describe the variation of baryonic effects in the observables, and the first PC amplitude (Q1) generally reflects the strength of baryon feedback. With the upper limit of Q1 prior being bound by the Illustris feedback scenarios, we reach $\sim 20{{\ \rm per\ cent}}$ improvement in the constraint of $S_8=\sigma _8(\Omega _{\rm m}/0.3)^{0.5}=0.788^{+0.018}_{-0.021}$ compared to the original DES 3 × 2pt analysis. This gain is driven by the inclusion of small-scale cosmic shear information down to 2.5 arcmin, which was excluded in previous DES analyses that did not model baryonic physics. We obtain $S_8=0.781^{+0.014}_{-0.015}$ for the combined DES Y1+Planck EE+BAO analysis with a non-informative Q1 prior. In terms of the baryon constraints, we measure $Q_1=1.14^{+2.20}_{-2.80}$ for DES Y1 only and $Q_1=1.42^{+1.63}_{-1.48}$ for DESY1+Planck EE+BAO, allowing us to exclude one of the most extreme AGN feedback hydrodynamical scenario at more than 2σ.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Cooling and instabilities in colliding flows

ABSTRACT Collisional self-interactions occurring in protostellar jets give rise to strong shocks, the structure of which can be affected by radiative cooling within the flow. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in both one and three dimensions with a power-law cooling function. The characteristic length and time-scales for cooling are temperature dependent and thus may vary as shocked gas cools. When the cooling length decreases sufficiently and rapidly, the system becomes unstable to the radiative shock instability, which produces oscillations in the position of the shock front; these oscillations can be seen in both the one- and three-dimensional cases. Our simulations show no evidence of the density clumping characteristic of a thermal instability, even when the cooling function meets the expected criteria. In the three-dimensional case, the nonlinear thin shell instability (NTSI) is found to dominate when the cooling length is sufficiently small. When the flows are subjected to the radiative shock instability, oscillations in the size of the cooling region allow NTSI to occur at larger cooling lengths, though larger cooling lengths delay the onset of NTSI by increasing the oscillation period.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Alignment of the magnetic field in star-forming regions and why it might be difficult to observe

ABSTRACT Magnetic fields are an important component of the interstellar medium (ISM) and exhibit strongly varying field strengths and a non-trivial correlation with the gas density. Its dynamical impact varies between individual regions of the ISM and correlates with the orientation of the field with respect to the gas structures. Using high-resolution magneto-hydrodynamical simulations of the ISM, we explore the connection between the orientation of the field and the dynamical state of the gas. We find that the onset of gravitational instability in molecular gas above a density of $\rho \sim 10^{-21}\, \mathrm{g\, cm}^{-3}$$(n\sim 400\, \mathrm{cm}^{-3})$ coincides with an alignment of the magnetic field lines and the gas flow. At this transition, the gradient of the density changes from mainly perpendicular to preferentially parallel to the field lines. A connection between the three-dimensional alignment and projected two-dimensional observables is non-trivial, because of a large dispersion of the magnetic field orientation along the line of sight. The turbulent correlation lengths can be small compared to the typical integration lengths. As a consequence, the small-scale signal of the orientation can sensitively depend on the line of sight or the dynamical state of the cloud can fluctuate stochastically or be completely averaged out. With higher spatial resolution more small-scale structures are resolved, which aggravates the link between magneto-hydrodynamical quantities and projected observables.

Girichidis, Philipp (ORCID:0000000293009914)↗
DOE OSTI2023

Turbulent dynamo action and its effects on the mixing at the convective boundary of an idealized oxygen-burning shell

Convection is one of the most important mixing processes in stellar interiors. Hydrodynamic mass entrainment can bring fresh fuel from neighboring stable layers into a convection zone, modifying the structure and evolution of the star. Because flows in stellar convection zones are highly turbulent, multidimensional hydrodynamic simulations are fundamental to accurately capture the physics of mixing processes. Under some conditions, strong magnetic fields can be sustained by the action of a turbulent dynamo, adding another layer of complexity and possibly altering the dynamics in the convection zone and at its boundaries. In this study, we used our fully compressible S EVEN -L EAGUE H YDRO code to run detailed and highly resolved three-dimensional magnetohydrodynamic simulations of turbulent convection, dynamo amplification, and convective boundary mixing in a simplified setup whose stratification is similar to that of an oxygen-burning shell in a star with an initial mass of 25 M ⊙ . We find that the random stretching of magnetic field lines by fluid motions in the inertial range of the turbulent spectrum (i.e., a small-scale dynamo) naturally amplifies the seed field by several orders of magnitude in a few convective turnover timescales. During the subsequent saturated regime, the magnetic-to-kinetic energy ratio inside the convective shell reaches values as high as 0.33, and the average magnetic field strength is ∼10 10 G. Such strong fields efficiently suppress shear instabilities, which feed the turbulent cascade of kinetic energy, on a wide range of spatial scales. The resulting convective flows are characterized by thread-like structures that extend over a large fraction of the convective shell. The reduced flow speeds and the presence of magnetic fields with strengths up to 60% of the equipartition value at the upper convective boundary diminish the rate of mass entrainment from the stable layer by ≈20% as compared to the purely hydrodynamic case.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Observational signatures of circumbinary discs – I. Kinematics

We present five morphological and kinematic criteria to aid in asserting the binary nature of a protoplanetary disc, based on 3D hydrodynamical simulations of circumbinary discs post-processed with Monte Carlo radiative transfer. We find that circumbinary discs may be identified by (i) a central cavity, (ii) spiral arms both in and outside of their central cavities, (iii) non-localized perturbations in their iso-velocity curves, (iv) asymmetry between the lines of maximum speed of the blue shifted and redshifted wings, and (v) asymmetry between the area of the blue shifted and redshifted wings. In this study, we provide quantitative metrics for the last two criteria that can be used, in conjunction with the morphological criteria, to signal whether a protoplanetary disc is likely to be a circumbinary disc.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Production efficiencies of sound waves in the intracluster medium driven by AGN jets

Feedback from active galactic nuclei (AGNs) is believed to be the most promising solution to the cooling flow problem in cool-core clusters. Dissipation of sound waves is considered as one of the possible heating mechanisms; however, its relative contribution to heating remains unclear. To estimate the energy budget for heating, we perform 3D hydrodynamic simulations of AGN jet injections in a Perseus-like cluster and quantify the amount of energy stored in the forms of weak shocks and waves. For this work, e find that, for a single jet injection with typical parameters in cool-core clusters, ${\sim}9{{\ \rm per\ cent}}$ of the total jet energy is stored in compressional waves (including both shocks and waves). However, due to the destructive effects among randomly phased waves as well as the dissipation of shock energies, in our simulations including self-regulated AGN feedback, no more than $3{{\ \rm per\ cent}}$ of the total injected energy goes into compressional waves. We further separate the energy contribution from shocks and waves and find that, for a single outburst, the shocks can only contribute to ${\sim}20{-}30{{\ \rm per\ cent}}$ of the total compressional energy in the inner radii and quickly dissipate away. In the self-regulated case where shocks are repeatedly generated, shocks completely dominate over sound waves in the inner region and can still provide ${\sim}40{-}50{{\ \rm per\ cent}}$ of the total compressional energy at outer radii. Our results suggest that the production of sound waves is not as efficient as what was previously found, and thus sound wave dissipation may be a subdominant source of heating in cool-core clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The landscape of disc outflows from black hole–neutron star mergers

ABSTRACT We investigate mass ejection from accretion discs formed in mergers of black holes (BHs) and neutron stars (NSs). The third observing run of the LIGO/Virgo interferometers provided BH–NS candidate events that yielded no electromagnetic (EM) counterparts. The broad range of disc configurations expected from BH–NS mergers motivates a thorough exploration of parameter space to improve EM signal predictions. Here we conduct 27 high-resolution, axisymmetric, long-term hydrodynamic simulations of the viscous evolution of BH accretion discs that include neutrino emission/absorption effects and post-processing with a nuclear reaction network. In the absence of magnetic fields, these simulations provide a lower limit to the fraction of the initial disc mass ejected. We find a nearly linear inverse dependence of this fraction on disc compactness (BH mass over initial disc radius). The dependence is related to the fraction of the disc mass accreted before the ouflow is launched, which depends on the disc position relative to the innermost stable circular orbit. We also characterize a trend of decreasing ejected fraction and decreasing lanthanide/actinide content with increasing disc mass at fixed BH mass. This trend results from a longer time to reach weak freezout and an increasingly dominant role of neutrino absorption at higher disc masses. We estimate the radioactive luminosity from the disc outflow alone available to power kilonovae over the range of configurations studied, finding a spread of two orders of magnitude. For most of the BH–NS parameter space, the disc outflow contribution is well below the kilonova mass upper limits for GW190814.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Universality of scaled particle spectra in ultrarelativistic heavy-ion collisions

We study the transverse momentum spectra of identified particles in ultrarelativistic collisions of large and small collision systems. In order to isolate information contained in the momentum dependence, we propose to scale the spectra by the total particle number and mean transverse momentum, global quantities which are already well studied. We observe an interesting, nearly universal, centrality-independent shape in the scaled spectra, similar to scalings that have been studied previously. This scaling behavior breaks down at large transverse momentum and for very small systems, such as those produced in p − p collisions. We perform hybrid hydrodynamic simulations and show that, in these simulations, a centrality-independent shape is a consequence of an event-by-event independence. Here, our results motivate further theoretical and experimental investigations of the regime of validity of this scaling phenomenon and their physical interpretation at different collision energies and systems.

Collective flow↗
DOE OSTI2025

A Wind-driven Origin for the Firework Morphology of the Supernova Remnant Pa 30

Pa 30—the likely remnant of the Galactic Type Iax supernova of 1181 AD—displays an unusual, firework-like morphology, consisting of radial filaments extending from a common center, where a white dwarf (WD) currently drives a very fast wind (speed ≳10 4 km s −1 ). We propose the filaments arose from the Rayleigh–Taylor-unstable nature of the interface between the circumstellar medium (CSM) and the shocked wind launched by the natal WD; the filaments then elongated intact due to the Kelvin–Helmholtz-stable nature of the large initial density contrast between the wind and CSM, supplemented by the slowly declining wind density profile (relative to homologously expanding ejecta). To support this interpretation, we present two-dimensional hydrodynamical simulations and derive the filament properties, including their speed, density, and temperature, all of which are consistent with observations. We suggest the filaments elongate until the wind and CSM densities become comparable at the contact discontinuity, which occurs within 1–10 yr, and then truncate because the Rayleigh–Taylor instability (RTI) halts. The subsequent Kelvin–Helmholtz instability growth timescale across the current width of the filaments is longer than the age of Pa 30, so they remain intact. The filament-less central region in Pa 30 is therefore more likely a consequence of the finite timescale over which the RTI operates, rather than a wind termination shock. In general, firework-like filaments may form in other systems, provided there is a sufficiently large density contrast between the ejecta and its surroundings.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Capture and migration of Jupiter and Saturn in mean motion resonance in a gaseous protoplanetary disc

Here, we study the dynamical evolution of Jupiter and Saturn embedded in a gaseous, solar nebula-type disc by means of hydrodynamics simulations with the FARGO2D1D code. We study the evolution for different initial separations of the planets’ orbits, Δa SJ , to investigate whether they become captured in mean motion resonance (MMR) and the direction of the subsequent migration of the planet (inwards or outwards). We also provide an assessment of the planet’s orbital dynamics at different epochs of Saturn’s growth. We find that the evolution of initially compact orbital configurations is dependent on the value of Δa SJ . This implies that an evolution as that proposed in the Grand Tack model depends on the precise initial orbits of Jupiter and Saturn and on the time-scales for their formation. Capture in the 1:2 MMR and inward or (nearly) stalled migration are highly favoured. Within its limits, our work suggests that the reversed migration, associated with the resonance capture of Jupiter and Saturn, may be a low-probability evolutionary scenario, so that other planetary systems with giant planets are not expected to have experienced a Grand Tack-like evolutionary path.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Heavy Elements and Electromagnetic Transients from Neutron Star Mergers

Compact binary mergers involving neutron stars can eject a fraction of their mass to space. Being extremely neutron rich, this material undergoes rapid neutron capture nucleosynthesis, and the resulting radioactivity powers fast, short-lived electromagnetic transients known as kilonova or macronova. Such transients are exciting probes of the most extreme physical conditions and their observation signals the enrichment of the Universe with heavy elements. Here the current understanding of the mass ejection mechanisms, the properties of the ejecta, and the resulting radioactive transients are reviewed. The first well-observed event in the aftermath of GW170817 delivered a wealth of insights, but much of today's picture of such events is still based on a patchwork of theoretical studies. Apart from summarizing the current understanding, questions where no consensus has been reached yet are also pointed out, and possible directions for the future research are sketched. In an appendix, a publicly available heating rate library based on the WinNet nuclear reaction network is described, and a simple fit formula to alleviate the implementation in hydrodynamic simulations is provided.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Revisiting the Kinematics of the Cylinder Test

The cylinder expansion experiment is a well-established performance test for condensed explosives that is utilized routinely to determine pressure-energy-volume relationships for detonation products. The modern cylinder test employs optical interferometric techniques to measure velocity, as opposed to older realizations where streak cameras were commonplace. Despite their widespread use, questions sometimes remain as to what kind of data the velocity diagnostics in a cylinder test provide along with how to interpret such measurements and relate them to physical quantities of interest. Here in this study, equations are derived that fully describe the kinematics of the cylinder wall during expansion. These equations are then applied to experimental as well as synthetic data generated via a hydrodynamic simulation in order to verify the mathematical framework developed.

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

Sediment transport and morphodynamic changes due to upstream variation and sea level rise on the Western Mekong Delta, Vietnam

The hydrological and morphological dynamics across the Vietnam Mekong Delta (VMD) coastline are increasingly shaped by human interventions and climate change, affecting the sustainability of socio-economic environments. This study employs the TELEMAC-MASCARET modeling system (TELEMAC-2D, TOMAWAC, and SISYPHE modules) to simulate hydrodynamics, sediment transport, and morphological changes for a 10-year period under baseline (SC1) and future climate (SC2) scenarios, considering wave climatology of Northeast (NEM) and Southwest monsoons (SWM). The results of sedimentation and erosion were analyzed across five key coastal zones from Ca Mau Cape (KV1) to Ha Tien (KV5). Model results highlight distinct seasonal patterns of sediment redistribution shaped by monsoonal flows, wave action, and seabed morphology. During the NEM, sediment deposition is concentrated near estuarine zones, while the SWM favors offshore sediment accumulation due to perpendicular wave actions, and shallow, gently sloping bathymetry that dissipates wave energy. Coastal erosion predominantly occurs at convex shoreline features subjected to stronger coastal currents, particularly in fine-grained sediment areas. Under future sea level rise and land subsidence scenarios, reduced flow velocities lead to lower bottom shear stress, promoting deposition and mitigating erosion problems despite increased wave energy. Besides that, sediment transport analyses across 14 cross-sections (C-S) reveal a southward trend during the NEM and a northward reversal during the SWM, with upstream sediment reductions primarily influencing near-estuarine zones and exerting minimal impact beyond the mid-coastal regions. Furthermore, these findings advance the understanding of sediment redistribution processes under climate change and seasonal change, and provide a scientific basis for adaptive coastal management and sustainable development.

Climate change↗