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 OSTI2025

Reduced Model of Ionization Lag in Intense Laser-Produced Plasmas

A physics-based empirical formula is derived to predict the ionization lag in underdense plasmas generated by an intense laser. Time-dependent nonlocal thermodynamic equilibrium calculations demonstrate significantly delayed ionization, due to rapid changes in plasma conditions, which critically impacts plasma properties such as opacity, emissivity, and heat transport. The reduced model, based on these calculations, enables the estimation of ionization lag without requiring in-depth knowledge of nonlocal thermodynamic equilibrium modeling. Furthermore, modeling reveals that the two-step ionization process—collisional excitation followed by photoionization—plays a crucial role in this ionization delay, with collisional excitation setting the timescale for ionization. Simulations across a range of elements, from beryllium to germanium, demonstrate that ionization lag is a widespread phenomenon, underscoring the importance of incorporating such processes into ionization modeling in radiation hydrodynamic simulations for various laser-plasma experiments.

direct drive↗
NASA NTRS1988

Observations and simulations of recurrent novae: U Sco and V394 CrA

Observations and analysis of the Aug. 1987 outburst of the recurrent nova V394 CrA are presented. This nova is extremely fast and its outburst characteristics closely resemble those of the recurrent nova U Sco. Hydrodynamic simulations of the outbursts of recurrent novae were performed. Results as applied to the outbursts of V394 CrA and U Sco are summarized.

Starrfield, S.↗
DOE OSTI2020

Discovery of a Candidate Binary Supermassive Black Hole in a Periodic Quasar from Circumbinary Accretion Variability

Binary supermassive black holes (BSBHs) are expected to be a generic byproduct from hierarchical galaxy formation. The final coalescence of BSBHs is thought to be the loudest gravitational wave (GW) siren, yet no confirmed BSBH is known in the GW-dominated regime. While periodic quasars have been proposed as BSBH candidates, the physical origin of the periodicity has been largely uncertain. Here we report discovery of a periodicity (P=1607±7 days) at 99.95% significance (with a global p-value of ~10 –3 accounting for the look elsewhere effect) in the optical light curves of a redshift 1.53 quasar, SDSS J025214.67–002813.7. Combining archival Sloan Digital Sky Survey data with new, sensitive imaging from the Dark Energy Survey, the total ~20-yr time baseline spans ~4.6 cycles of the observed 4.4-yr (restframe 1.7-yr) periodicity. The light curves are best fit by a bursty model predicted by hydrodynamic simulations of circumbinary accretion disks. The periodicity is likely caused by accretion rate modulation by a milli-parsec BSBH emitting GWs, dynamically coupled to the circumbinary accretion disk. A bursty hydrodynamic variability model is statistically preferred over a smooth, sinusoidal model expected from relativistic Doppler boost, a kinematic effect proposed for PG1302–102. Furthermore, the frequency dependence of the variability amplitudes disfavors Doppler boost, lending independent support to the circumbinary accretion variability hypothesis. Given our detection rate of one BSBH candidate from circumbinary accretion variability out of 625 quasars, it suggests that future large, sensitive synoptic surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time may be able to detect hundreds to thousands of candidate BSBHs from circumbinary accretion with direct implications for Laser Interferometer Space Antenna.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Uncertainties in Density and Simulation Parameters for Radiographic Reconstructions Using Machine Learning

We will develop new and improved methods that will enable a fuller and more extensive use of experimental data, such as from radiographic imaging, towards better characterizing and reducing uncertainties in predictive modeling of weapons performance. We will achieve this by leveraging recent developments in the areas of computational imaging, statistical and machine learning, and reduced order modeling. Key new aspects of our approach include (a) a coupling of deep learning based density reconstruction with fast hydrodynamics simulators that will permit the development of certain guarantees on extrapolatability and bounds on certain other hydrodynamics-related para- metric uncertainties and (b) a dual approach towards better characterizing uncertainties in density reconstruction using deep learning based surrogates to accelerate Bayesian reconstructions on the one hand and using variational inference and Bayesian machine learning on the other hand.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2026

Passive Freeze-Out of the Richtmyer-Meshkov Instability

The Richtmyer-Meshkov instability (RMI) poses a major challenge in inertial confinement fusion (ICF) due to its role in mixing and performance degradation. We report the first experimental observation of passive freeze-out of RMI in a low-pressure surrogate regime, an instability stagnation effect induced without modifying the driving pressure pulse or the target surface geometry. Using additively manufactured subsurface voids in a sinusoidal target, we convert a single shock into a sequence of weaker shocks that suppress instability growth upstream of the surface by over 70%. High-speed x-ray imaging and hydrodynamic simulations suggest that this suppression arises primarily from temporal shaping, with lesser contributions from spatial curvature and shock weakening. Our results demonstrate a driver-independent pathway for controlling shock-driven hydrodynamic instabilities relevant to ICF and other high energy density systems.

Materials science↗
DOE OSTI2023

CFD Benchmark of Pressure Drop in a 61-Pin Wire-Wrapped Assembly with Blocked Channels Using NekRS

Thermal hydraulic behaviors of wire-wrapped rod bundles continue to be a subject of intense research. One of the leading next-generation designs, the sodium fast reactor, relies on a hexagonal assembly of wire-wrapped fuel pins. The issue of potential channel blockage has been raised as a safety concern due to the tightly packed arrangement of the fuel pins. This has led to several recent experimental and computational studies working to quantify the potential impact on the fluid flow and heat transfer behaviors of such blockages. The objective of the present study is to benchmark the high-fidelity NekRS CFD solver in predicting pressure drop for large blockages against available experimental data. A 61-pin wire-wrapped fuel assembly with two flow blockage configurations has been simulated and investigated at various low to moderate Reynolds numbers. The NekRS solver has been shown to yield exponentially decreasing spatial discretization errors with increasing polynomial order. All simulated results agreed well with measured data, which indicates that the overall methodology is adequate and consistent. The results of this benchmark study demonstrate the accuracy of NekRS for sodium fast reactor hydrodynamic simulations, increasing the confidence in its use for design, licensing, and analysis activities.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1992

Statistics of the cosmic Mach number from numerical simulations of a cold dark matter universe

Results are presented of an analysis of the cosmic Mach number, M, the ratio of the streaming velocity, v, to the random velocity dispersion, sigma, of galaxies in a given patch of the universe, which was performed on the basis of hydrodynamical simulations of the cold dark matter scenario. Galaxy formation is modeled by application of detailed physical processes rather than by the ad hoc assumption of 'bias' between dark matter and galaxy fluctuations. The correlation between M and sigma is found to be very weak for both components. No evidence is found for a physical 'velocity bias' in the quantities which appear in the definition of M. Standard cold-dark-matter-dominated universes are in conflict, at a statistically significant level, with the available observation, in that they predict a Mach number considerably lower than is observed.

Suto, Yasushi↗
DOE OSTI2025

Simulation of a High-Confinement Tidal Energy Array (Abstract)

When the projected area of a current turbine array occupies a substantial fraction of a tidal or river channel cross-sectional area, the confinement associated with the channel boundaries substantially increases turbine efficiency. This effect can be exploited to substantially reduce the levelized cost of energy for current turbine arrays. However, the presence of such an array will affect inflow velocity, water depth, and circulation in the surrounding area. This has implications for array performance and environmental impacts. This project combines existing numerical modeling capability from PNNL with an understanding of confined array performance from UW. The work will 1) Simulate tidal currents and water surface elevation in a regional circulation model that includes the Knik Arm site at appropriate resolution; 2) Incorporate a varying thrust coefficient into this model that reflects the outcomes of tidal hydrodynamic simulations of confined turbine arrays; 3) Use the model to determine the effects of array operation on water levels and inflow conditions as a function of different blockage ratios and control strategies; and 4) Identify the “far field” environmental implications of operating an array of this type. UW will use the outputs from this model to evaluate the implications of “feedback effects” and variable inflow conditions on levelized cost of energy estimates for a confinement exploiting array of cross-flow turbines.

16 TIDAL AND WAVE POWER↗
NASA NTRS2014

Water Delivery and Giant Impacts in the 'Grand Tack' Scenario

A new model for terrestrial planet formation has explored accretion in a truncated protoplanetary disk, and found that such a configuration is able to reproduce the distribution of mass among the planets in the Solar System, especially the Earth/Mars mass ratio, which earlier simulations have generally not been able to match. Walsh et al. tested a possible mechanism to truncate the disk-a two-stage, inward-then-outward migration of Jupiter and Saturn, as found in numerous hydrodynamical simulations of giant planet formation. In addition to truncating the disk and producing a more realistic Earth/Mars mass ratio, the migration of the giant planets also populates the asteroid belt with two distinct populations of bodies-the inner belt is filled by bodies originating inside of 3 AU, and the outer belt is filled with bodies originating from between and beyond the giant planets (which are hereafter referred to as 'primitive' bodies). One implication of the truncation mechanism proposed in Walsh et al. is the scattering of primitive planetesimals onto planet-crossing orbits during the formation of the planets. We find here that the planets will accrete on order 1-2% of their total mass from these bodies. For an assumed value of 10% for the water mass fraction of the primitive planetesimals, this model delivers a total amount of water comparable to that estimated to be on the Earth today. The radial distribution of the planetary masses and the dynamical excitation of their orbits are a good match to the observed system. However, we find that a truncated disk leads to formation timescales more rapid than suggested by radiometric chronometers. In particular, the last giant impact is typically earlier than 20 Myr, and a substantial amount of mass is accreted after that event. This is at odds with the dating of the Moon-forming impact and the estimated amount of mass accreted by Earth following that event. However, 5 of the 27 planets larger than half an Earth mass formed in all simulations do experience large late impacts and subsequent accretion consistent with those constraints.

Planetary dynamics↗
DOE OSTI2020

Fully compressible simulations of waves and core convection in main-sequence stars

Context: Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and requires rather high artificial viscosity and thermal diffusivity for numerical stability. A direct comparison with observations is thus limited. Aims: In this study, we explore the capabilities of our compressible multidimensional Seven-League Hydro (SLH) code to simulate stellar oscillations. Methods: We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic 3 M ⊙ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results: Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with the observed spectra of massive stars. Conclusion: The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows for the coupling of gravity and pressure waves in the outer convective envelopes of evolved stars to be studied in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Turbulently Driven Detonation Initiation in Electron-degenerate Matter with Helium

Abstract Type Ia supernovae (SNe Ia) are standardizable cosmological candles that led to the discovery of the accelerating Universe. However, the physics of how white dwarfs (WDs) explode and lead to SNe Ia is still poorly understood. The initiation of the detonation front that rapidly disrupts the WD is a crucial element of the puzzle, and global 3D simulations of SNe Ia cannot resolve the requisite length scales to capture detonation initiation. In this work, we elucidate a theoretical criterion for detonation initiation in the distributed burning regime. We test this criterion against local 3D driven turbulent hydrodynamical simulations within electron-degenerate WD matter consisting initially of pure helium. We demonstrate a novel pathway for detonation, in which strong turbulent dissipation rapidly heats the helium, and forms carbon nuclei sufficient to lead to a detonation through accelerated burning via α captures. Simulations of strongly driven turbulent conditions lead to detonations at a mean density of 10 6 g cm −3 and mean temperature of 1.4–1.8 × 10 9 K, but fail to detonate at a lower density of 10 5 g cm −3 , in excellent agreement with theoretical predictions.

Astronomy & Astrophysics↗
DOE OSTI2024

Systematic Local Simulations of Fast Neutrino Flavor Conversions with Scattering Effects

We investigate the dynamics of fast neutrino flavor conversions (FFCs) in the one-dimensional (1D) and zero-dimensional (0D) models, in which spatial advection is considered and ignored, respectively. In this study, we employ snapshots obtained by our self-consistent, realistic Boltzmann-neutrino-radiation-hydrodynamics simulations. We show that the FFC growth rate is considerably larger in the 1D model than in the 0D model, as expected from the previous linear analysis results. We find that the momentum space dimension does not significantly influence the neutrino transition probability in 1D models. On the other hand, in the 0D model without collisions, the FFC depends on the momentum space, and the azimuthal angle dependence breaks the periodicity of the FFC. Our study demonstrates that collisional instability can lead to further flavor conversions on a long timescale in 1D models after the asymptotic state of FFC has been reached. Such an effect should be taken into consideration when the fast and collisional flavor instabilities coexist.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2016

Observations and Simulations of the Na I D-1 Line profiles in an M-Class Solar Flare

We study the temporal evolution of the Na I D1 line profiles in the M3.9 flare SOL2014-06-11T21:03UT, using observations at high spectral resolution obtained with the Interferometric Bidimensional Spectrometer instrumentation the Dunn Solar Telescope combined with radiative hydrodynamic simulations. Our results show a significant increase in the intensities of the line core and wings during the flare. The analysis of the line profiles from the flare ribbons reveals that the Na I D1 line has a central reversal with excess emission in the blue wing (blue asymmetry).We combine RADYN and RH simulations to synthesize Na I D1 line profiles of the flaring atmosphere and find good agreement with the observations. Heating with a beam of electrons modifies the radiation field in the flaring atmosphere and excites electrons from the ground state 3s 2S to the first excited state 3p 2P, which in turn modifies the relative population of the two states. The change in temperature and the population density of the energy states make the sodium line profile revert from absorption into emission. Furthermore, the rapid changes in temperature break the pressure balance between the different layers of the lower atmosphere, generating upflow/downflow patterns. Analysis of the simulated spectra reveals that the asymmetries of the Na I D1 flare profile are produced by the velocity gradients in the lower solar atmosphere.

methods: numerical↗
DOE OSTI2023

The Concentration–Mass relation of massive, dynamically relaxed galaxy clusters: agreement between observations and ΛCDM simulations

ABSTRACT The relationship linking a galaxy cluster’s total mass with the concentration of its mass profile and its redshift is a fundamental prediction of the Cold Dark Matter (CDM) paradigm of cosmic structure formation. However, confronting those predictions with observations is complicated by the fact that simulated clusters are not representative of observed samples where detailed mass profile constraints are possible. In this work, we calculate the symmetry-peakiness-alignment (SPA) morphology metrics for maps of X-ray emissivity from the three hundred project hydrodynamical simulations of galaxy clusters at four redshifts, and thereby select a sample of morphologically relaxed, simulated clusters, using observational criteria. These clusters have on average earlier formation times than the full sample, confirming that they are both morphologically and dynamically more relaxed than typical. We constrain the concentration–mass–redshift relation of both the relaxed and complete sample of simulated clusters, assuming power-law dependences on mass (κm) and 1 + z (κζ), finding κm = −0.12 ± 0.07 and κζ = −0.27 ± 0.19 for the relaxed subsample. From an equivalently selected sample of massive, relaxed clusters observed with Chandra, we find κm = −0.12 ± 0.08 and κζ = −0.48 ± 0.19, in good agreement with the simulation predictions. The simulated and observed samples also agree well on the average concentration at a pivot mass and redshift providing further validation of the ΛCDM paradigm in the properties of the largest gravitationally collapsed structures observed. This also represents the first clear detection of decreasing concentration with redshift, a longstanding prediction of simulations in data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Unraveling electronic correlations in warm dense quantum plasmas

The study of matter at extreme densities and temperatures has emerged as a highly active frontier at the interface of plasma physics, material science and quantum chemistry with relevance for planetary modeling and inertial confinement fusion. A particular feature of such warm dense matter is the complex interplay of Coulomb interactions, quantum effects, and thermal excitations, making its rigorous theoretical description challenging. Here, we demonstrate how ab initio path integral Monte Carlo simulations allow us to unravel this intricate interplay for the example of strongly compressed beryllium, focusing on two X-ray Thomson scattering data sets obtained at the National Ignition Facility. We find excellent agreement between simulation and experiment with a very high level of consistency between independent observations without the need for any empirical input parameters. Our results call into question previously used chemical models, with important implications for the interpretation of scattering experiments and radiation hydrodynamics simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Impact of Self-shielding Minihalos on the Lyα Forest at High Redshift

Abstract Dense gas in minihalos with masses of 10 6 −10 8 M ⊙ can shield themselves from reionization for ∼100 Myr after being exposed to the UV background. These self-shielded systems, often unresolved in cosmological simulations, can introduce strong absorption in quasar spectra. This paper is the first systematic study on the impact of these systems on the Ly α forest. We first derive the H i column density profile of photoevaporating minihalos by conducting 1D radiation–hydrodynamics simulations. We utilize these results to estimate the Ly α opacity from minihalos in a large-scale simulation that cannot resolve self-shielding. When the ionization rate of the background radiation is 0.03 × 10 −12 s −1 , as expected near the end of reionization at z ∼ 5.5, we find that the incidence rate of damped Ly α absorbers increases by a factor of ∼2−4 compared to at z = 4.5. The Ly α flux is, on average, suppressed by ∼3% of its mean due to minihalos. The absorption features enhance the 1D power spectrum up to ∼5% at k ∼ 0.1 h Mpc −1 (or 10 −3 km −1 s), which is comparable to the enhancement caused by inhomogeneous reionization. The flux is particularly suppressed in the vicinity of large halos along the line-of-sight direction at separations of up to 10 h −1 Mpc at r ⊥ ≲ 2 h −1 Mpc. However, these effects become much smaller for higher ionizing rates (≳0.3 × 10 −12 s −1 ) expected in the post-reionization Universe. Our findings highlight the need to consider minihalo absorption when interpreting the Ly α forest at z ≳ 5.5. Moreover, the sensitivity of these quantities to the ionizing background intensity can be exploited to constrain the intensity itself.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRSDate not supplied

Impact Of Rotation on the Sun’s Subsurface Structure

It is well known that the dynamical evolution of subsurface flows of the Sun and magnetic fields is strongly linked to rotation. Observational investigation of sub surface structuring and dynamics can be performed only with relatively low resolution using helioseismic inversion techniques. To confidently model the coupling of turbulent convection, magnetic fields, and rotation in the solar interior, we perform advanced radiative hydrodynamic and MHD modeling of local scales. To explore the coupling effects in the upper layers of the convection zone, we present results of analysis from a series of hydrodynamic simulations performed for a 25-Mm-deep computational domain at different latitudes obtained for 1) the Carrington rotation rate at all latitudes and 2) the realistic rotation rate from the global helioseismic inversions during the solar minimum. Using additional 3D radiative hydrodynamic models, we reveal how physical properties of the solar convection are affected by rotation.

Heliophysics↗
NASA NTRS2014

Kinematic Evolution of Simulated Star-Forming Galaxies

Recent observations have shown that star-forming galaxies like our own Milky Way evolve kinematically into ordered thin disks over the last approximately 8 billion years since z = 1.2, undergoing a process of "disk settling." For the first time, we study the kinematic evolution of a suite of four state of the art "zoom in" hydrodynamic simulations of galaxy formation and evolution in a fully cosmological context and compare with these observations. Until now, robust measurements of the internal kinematics of simulated galaxies were lacking as the simulations suffered from low resolution, overproduction of stars, and overly massive bulges. The current generation of simulations has made great progress in overcoming these difficulties and is ready for a kinematic analysis. We show that simulated galaxies follow the same kinematic trends as real galaxies: they progressively decrease in disordered motions (sigma(sub g)) and increase in ordered rotation (V(sub rot)) with time. The slopes of the relations between both sigma(sub g) and V(sub rot) with redshift are consistent between the simulations and the observations. In addition, the morphologies of the simulated galaxies become less disturbed with time, also consistent with observations. This match between the simulated and observed trends is a significant success for the current generation of simulations, and a first step in determining the physical processes behind disk settling.

Evolution↗