Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamics methods”

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

Nebular Models of Sub-Chandrasekhar Mass Type Ia Supernovae: Clues to the Origin of Ca-rich Transients

In this work we use non-local thermal equilibrium radiative transport modeling to examine observational signatures of sub-Chandrasekhar mass double detonation explosions in the nebular phase. Results range from spectra that look like typical and subluminous Type Ia supernovae (SNe) for higher mass progenitors to spectra that look like Ca-rich transients for lower mass progenitors. This ignition mechanism produces an inherent relationship between emission features and the progenitor mass as the ratio of the nebular [Ca II ]/[Fe III ] emission lines increases with decreasing white dwarf mass. Examining the [Ca II ]/[Fe III ] nebular line ratio in a sample of observed SNe we find further evidence for the two distinct classes of SNe Ia identified in Polin et al. by their relationship between Si ii velocity and B-band magnitude, both at time of peak brightness. This suggests that SNe Ia arise from more than one progenitor channel, and provides an empirical method for classifying events based on their physical origin. Furthermore, we provide insight to the mysterious origin of Ca-rich transients. Low-mass double detonation models with only a small mass fraction of Ca (1%) produce nebular spectra that cool primarily through forbidden [Ca II ] emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Peculiarities of planar shockwave interaction with air–water interface and solid target

We present the results of the research on underwater electrical explosion of planar copper wire arrays, accompanied by the generation of a planar shock, using the synchrotron-based phase-contrast radiography imaging capabilities of the ID19 beamline at the European Synchrotron Radiation Facility. It is shown that the interaction of a strong shock with an air–water interface leads to additional acceleration of the wire expansion due to a rarefaction wave along with the emergence of micro-jets. In the case of a target placed above the array, we observed the formation of a cavity between the array and the target due to the interaction of two rarefaction waves. Furthermore, the results of two-dimensional hydrodynamic simulations of the wire explosion and the interaction of the generated shock with the water–air and water-target interfaces showed good agreement with experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Recent progress on coarse graining simulations

We focus on coarse graining simulations based on the primary conservation equations, effectively codesigned physics and algorithms, and low-Mach-number corrected (LMC) hydrodynamics. Simulation methods involve LANL’s x-Radiation-Adaptive-Grid-Eulerian Large-Eddy Simulation, Besnard-Harlow-Rauenzahn (BHR) Reynolds-Averaged Navier-Stokes (RANS) approach, and Dynamic BHR – a paradigm bridging RANS and LES. A relevant question addressed relates to whether 3D RANS and RANS/LES hybrids – the industry standards for aerospace and automotive research, are presently relevant for practical variable-density applications involving shocked and accelerated interface instabilities. Furthermore, recent simulations of the GaTECH inclined mixing-layer shock-tube and NIF ICF-capsule experiments are used to demonstrate issues, challenges, and potential for 3D coarse grained LMC simulation strategies for robustly simulating complex transitional and coupled hydrodynamics-multiphysics with coarser resolution. Present LES readiness to provide accurate predictions at scale is demonstrated – whereas 3D RANS and RANS/LES bridging do not appear impactful in this context.

42 ENGINEERING↗
DOE OSTI2020

Dynamics of Laterally Propagating Flames in X-Ray Bursts. I. Burning Front Structure

We investigate the structure of laterally propagating flames through the highly stratified burning layer in an X-ray burst. Two-dimensional hydrodynamics simulations of flame propagation are performed through a rotating plane-parallel atmosphere, exploring the structure of the flame. We discuss the approximations needed to capture the length and timescales at play in an X-ray burst and describe the flame acceleration observed. Finally, our studies complement other multidimensional studies of burning in X-ray bursts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Elucidating size effects on the yield strength of single-crystal Cu via the Richtmyer–Meshkov instability

Capturing the dynamic response of a material under high strain-rate deformation often demands challenging and time consuming experimental effort. While shock hydrodynamic simulation methods can aid in this area, a priori characterizations of the material strength under shock loading and spall failure are needed in order to parameterize constitutive models needed for these computational tools. Moreover, parameterizations of strain-rate-dependent strength models are needed to capture the full suite of Richtmyer–Meshkov instability (RMI) behavior of shock compressed metals, creating an unrealistic demand for these training data solely on experiments. Herein, we sweep a large range of geometric, crystallographic, and shock conditions within molecular dynamics (MD) simulations and demonstrate the breadth of RMI in Cu that can be captured from the atomic scale. In this work, yield strength measurements from jetted and arrested material from a sinusoidal surface perturbation were quantified as Y RMI = 0.787 ± 0.374 GPa, higher than strain-rate-independent models used in experimentally matched hydrodynamic simulations. Defect-free, single-crystal Cu samples used in MD will overestimate Y RMI , but the drastic scale difference between experiment and MD is highlighted by high confidence neighborhood clustering predictions of RMI characterizations, yielding incorrect classifications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗
NASA NTRS1926

A resume of the advances in theoretical aeronautics made by Max M. Munk

In order to apply profitably the mathematical methods of hydrodynamics to aeronautical problems, it is necessary to make simplifications in the physical conditions of the latter. To begin with, it is allowable in many problems, as Prandtl has so successfully shown, to treat the air as having constant density and as free of viscosity. But this is not sufficient. It is also necessary to specify certain shapes for the solid bodies whose motion through the air is discussed, shapes suggested by the actual solids - airships or airfoils - it is true, but so chosen that they lead to solvable problems. In a valuable paper presented by Dr. Max M. Munk, of the National Advisory Committee for Aeronautics, Washington, to the Delft Conference in April, 1924, these necessary simplifying assumptions are discussed in detail. It is the purpose of the present paper to present in as simple a manner as possible some of the interesting results obtained by Dr. Munk's methods.

Ames, Joseph S↗
DOE OSTI2021

Practical Effects of Integrating Temperature with Strang Split Reactions

For astrophysical reacting flows, operator splitting is commonly used to couple hydrodynamics and reactions. Each process operates independent of one another, but by staggering the updates in a symmetric fashion (via Strang splitting) second order accuracy in time can be achieved. However, approximations are often made to the reacting system, including the choice of whether or not to integrate temperature with the species. Here we demonstrate through a simple convergence test that integrating an energy equation together with reactions achieves the best convergence when modeling reactive flows with Strang splitting. Additionally, second order convergence cannot be achieved without integrating an energy or temperature equation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Dynamics of Laterally Propagating Flames in X-Ray Bursts. II. Realistic Burning and Rotation

In this work, we continue to investigate two-dimensional laterally propagating flames in type I X-ray bursts using fully compressible hydrodynamics simulations. In the current study we relax previous approximations where we artificially boosted the flames. We now use more physically realistic reaction rates, thermal conductivities, and rotation rates, exploring the effects of neutron star rotation rate and thermal structure on the flame. We find that at lower rotation rates the flame becomes harder to ignite, whereas at higher rotation rates the nuclear burning is enhanced by increased confinement from the Coriolis force and the flame propagates steadily. At higher crustal temperatures, the flame moves more quickly and accelerates as it propagates through the atmosphere. If the temperature is too high, instead of a flame propagating across the surface the entire atmosphere burns uniformly. Our findings could have implications for the relationship between observed burst rise times and neutron star rotation and accretion rates. All of the software used for these simulations is freely available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Neural Networks for Nuclear Reactions in MAESTROeX

We demonstrate the use of neural networks to accelerate the reaction steps in the MAESTROeX stellar hydrodynamics code. A traditional MAESTROeX simulation uses a stiff ODE integrator for the reactions; here, we employ a ResNet architecture and describe details relating to the architecture, training, and validation of our networks. Our customized approach includes options for the form of the loss functions, a demonstration that the use of parallel neural networks leads to increased accuracy, and a description of a perturbational approach in the training step that robustifies the model. We test our approach on millimeter-scale flames using a single-step, 3-isotope network describing the first stages of carbon fusion occurring in Type Ia supernovae. We train the neural networks using simulation data from a standard MAESTROeX simulation, and show that the resulting model can be effectively applied to different flame configurations. This work lays the groundwork for more complex networks, and iterative time-integration strategies that can leverage the efficiency of the neural networks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Effect of causality constraints on Bayesian analyses of heavy-ion collisions

There have long been questions about the limits to the validity of relativistic fluid dynamics and whether it is being used outside its regime of validity in modern simulations of relativistic heavy-ion collisions. An important new tool for answering this question is a causality analysis in the nonlinear regime—if the solutions of the evolution equations do not respect relativistic causality, then they are not a faithful representation of the underlying relativistic theory (in this case, quantum chromodynamics). Using this nonlinear criterion, it has recently been shown that hydrodynamics is indeed being used outside its regime of validity in simulations, at least sometimes. Here we explore the phenomenological implications, particularly the quantitative effects of demanding limits on acausality in modern Bayesian parameter estimation. We find that, while typically only a small fraction of the system's energy is initially in an acausal regime, placing strict limits on the allowed energy fraction significantly changes the preferred properties of the initial condition, which in turn alters the extracted medium properties such as bulk viscosity, where large values are no longer favored. Furthermore, these findings highlight the importance of developing better theoretical descriptions of the early-time, out-of-equilibrium dynamics of relativistic heavy-ion collisions.

Bayesian methods↗
NASA NTRS2008

Multi-wavelength Observations and Modelling of a Solar Flare

Aims: The aim of this work is to investigate the dynamic behavior of a C-class solar flare through the evolution of temperature, emission measure, energy loss and velocity. In particular, the variation of these properties with time are studied using multi-wavelength observations in combination with a recently developed 0-D hydrodynamic model. Methods: The temperature and emission measure evolution were studied using several instruments covering a wide range of temperatures - the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI, >5 MK), GOES-12 (5- 30 MK), the Transition Region and Coronal Explorer (TRACE 171 A, 1 MK) and the Coronal Diagnostic Spectrometer (CDS, 0.03-8 MK). The temperature and emission measure were analysed through the systematic cooling of flare plasma through the response functions of these instruments. These parameters were then investigated using the Enthalpy Based Thermal Evolution of Loops model (EBTEL). The Doppler shifts at both flare footpoints were analysed using five emission lines seen by CDS. Results: The flare began with clear evidence for pre-flare heating. Upflows of approx.90 km/s and low level emission, both observed in Fe XIX before the main impulsive phase were explained by pre-flare gentle chromospheric evaporation. During the main impulsive phase, the flare plasma was heated to a temperature of >13 MK in approximately 10 minutes. Explosive chromospheric evaporation was observed, driving upflows of approx.80 km/s in Fe XIX and simultaneous downflows of approx.20 km/s in He I and O v. At the peak of the Rare, conduction modelled by EBTEL was found to be the dominant loss mechanism, working efficiently to both lower the temperatures and drive gentle chromospheric evaporation. As the temperature fell below approx.8 MK, radiation became the dominant loss mechanism. During the final stages of the decay phase, downflowing plasma was observed at the footpoints in He I, O v and Mg x at velocities of up to approx.40 km/s, suggesting loop draining occurred. Conclusions. This is the first extensive study of the evolution of flare plasma using both spectroscopic and broad-band instruments in conjunction with a comprehensive hydrodynamic model. The flare began with pre-flare heating and then evolved following the predictions of the standard flare model. Detailed analysis of the plasma heating mechanisms was carried out and the heating function most consistent with observations was found to be Gaussian in shape. The simulations suggested that both direct heating and heating by a non-thermal beam played significant roles in this event.

Raftery, Claire L.↗
DOE OSTI2022

Inferring the Thermal History of the Intergalactic Medium from the Properties of the Hydrogen and Helium Lyα Forest

The filamentary network of intergalactic medium (IGM) gas that gives origin to the Lyα forest in the spectra of distant quasars encodes information on the physics of structure formation and the early thermodynamics of diffuse baryonic material. Here we use a massive suite of more than 400 high-resolution cosmological hydrodynamical simulations run with the Graphics Processing Unit–accelerated code Cholla to study the IGM at high spatial resolution maintained over the entire computational volume. The simulations capture a wide range of possible IGM thermal histories by varying the photoheating and photoionizing background produced by star-forming galaxies and active galactic nuclei. A statistical comparison of synthetic spectra with the observed 1D flux power spectra of hydrogen at redshifts 2.2 ≤ z ≤ 5.0 and with the helium Lyα opacity at redshifts 2.4 < z < 2.9 tightly constrains the photoionization and photoheating history of the IGM. By leveraging the constraining power of the available Lyα forest data to break model degeneracies, we find that the IGM experienced two main reheating events over 1.2 Gyr of cosmic time. For our best-fit model, hydrogen reionization completes by z R ≈ 6.0 with a first IGM temperature peak of T 0 ≃ 1.3 × 10 4 K and is followed by the reionization of He ii that completes by z R ≈ 3.0 and yields a second temperature peak of T 0 ≃ 1.4 × 10 4 K. We discuss how our results can be used to obtain information on the timing and the sources of hydrogen and helium reionization.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Modeling Solids in Nuclear Astrophysics with Smoothed Particle Hydrodynamics

Abstract Smoothed particle hydrodynamics (SPH) is a frequently applied tool in computational astrophysics to solve the fluid dynamics equations governing the systems under study. For some problems, for example when involving asteroids and asteroid impacts, the additional inclusion of material strength is necessary in order to accurately describe the dynamics. In compact stars, that is white dwarfs and neutron stars, solid components are also present. Neutron stars have a solid crust, which is the strongest material known in nature. However, their dynamical evolution, when modeled via SPH or other computational fluid dynamics codes, is usually described as a purely fluid dynamics problem. Here, we present the first 3D simulations of neutron star crustal toroidal oscillations including material strength with the Los Alamos National Laboratory SPH code FleCSPH. In the first half of the paper, we present the numerical implementation of solid material modeling together with standard tests. The second half is on the simulation of crustal oscillations in the fundamental toroidal mode. Here, we dedicate a large fraction of the paper to approaches that can suppress numerical noise in the solid. If not minimized, the latter can dominate the crustal motion in the simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Betelgeuse as a Merger of a Massive Star with a Companion

We investigate the merger between a 16M ⊙ star, on its way to becoming a red supergiant (RSG), and a 4M ⊙ main-sequence companion. Our study employs three-dimensional hydrodynamic simulations using the state-of-the-art adaptive mesh refinement code Octo-Tiger. The initially corotating binary undergoes interaction and mass transfer, resulting in the accumulation of mass around the companion and its subsequent loss through the second Lagrangian point (L2). The companion eventually plunges into the envelope of the primary, leading to its spin-up and subsequent merger with the helium core. We examine the internal structural properties of the post-merger star, as well as the merger environment and the outflow driven by the merger. Our findings reveal the ejection of approximately ∼0.6 M ⊙ of material in an asymmetric and somewhat bipolar outflow. We import the post-merger stellar structure into the MESA stellar evolution code to model its long-term nuclear evolution. In certain cases, the post-merger star exhibits persistent rapid equatorial surface rotation as it evolves in the H–R diagram toward the observed location of Betelgeuse. These cases demonstrate surface rotation velocities of a similar magnitude to those observed in Betelgeuse, along with a chemical composition resembling that of Betelgeuse. In other cases, efficient rotationally induced mixing leads to slower surface rotation. This pioneering study aims to model stellar mergers across critical timescales, encompassing dynamical, thermal, and nuclear evolutionary stages.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Violent mergers revisited: The origin of the fastest stars in the Galaxy

Binary systems of two carbon-oxygen white dwarfs are one of the most promising candidates for the progenitor systems of Type Ia supernovae. Violent mergers, where the primary white dwarf ignites when the secondary white dwarf smashes into it while being disrupted on its last orbit, were the first double degenerate merger scenario proposed that ignites dynamically. However, violent mergers likely contribute only a few percent to the total Type Ia supernova rate and do not yield normal Type Ia supernova light curves. Here we revisit the scenario, simulating a violent merger with better methods and, in particular, a more accurate treatment of the detonation. We find good agreement with previous simulations but with one critical difference: The secondary white dwarf being disrupted and accelerated towards the primary white dwarf and impacted by its explosion does not fully burn, and its core survives as a bound object. The explosion leaves behind a 0.16 M ⊙ star travelling 2800 km/s, making it an excellent (and so far the only) candidate to explain the origin of the fastest observed hypervelocity stars. We also show that before the explosion, 5 × 10 −3 M ⊙ of material predominantly consisting of helium, carbon, and oxygen had already been ejected at velocities above 1000 km/s. Finally, we argue that if a violent merger made the hypervelocity stars D6-1 and D6-3 and violent mergers require the most massive primary white dwarfs in binaries of two carbon-oxygen white dwarfs, there has to be a much larger population of white dwarf mergers with slightly lower mass primary white dwarfs. Because this population likely represents ≫10% of the Type Ia supernovae rate, it can essentially only give rise to normal Type Ia supernovae.

Astronomy and AstroPhysics↗
DOE OSTI2021

The discrete Green's function paradigm for two-way coupled Euler–Lagrange simulation

We outline a methodology for the simulation of two-way coupled particle-laden flows. The drag force that couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. In addition to hydrodynamics, the method can be extended to other physics present in particle-laden flows such as heat transfer and electromagnetism. The discrete Green's functions for the Navier–Stokes equations are obtained at low particle Reynolds number in a two-plane channel geometry. We perform verification at different Reynolds numbers for a particle settling under gravity parallel to a plane wall, for different wall-normal separations. Compared with other point-particle schemes, the Stokesian discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations. To account for degradation in accuracy away from the wall at finite Reynolds number, we extend the present methodology to an Oseen-like discrete Green's function. The extended discrete Green's function method is found to be accurate within 6% at all wall-normal separations for particle Reynolds numbers up to 24. Furthermore, the discrete Green's function approach is well suited to dilute systems with significant mass loading and this is highlighted by comparison against other Euler–Lagrange as well as particle-resolved simulations of gas–solid turbulent channel flow. Strong particle–turbulence coupling is observed in the form of turbulence modification and turbophoresis suppression, and these observations are placed in context of the different methods.

42 ENGINEERING↗
DOE OSTI2021

Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star

A realistic parametrization of convection and convective boundary mixing in conventional stellar evolution codes is still the subject of ongoing research. Furthermore, to improve the current situation, multidimensional hydrodynamic simulations are used to study convection in stellar interiors. Such simulations are numerically challenging, especially for flows at low Mach numbers which are typical for convection during early evolutionary stages. We explore the benefits of using a low-Mach hydrodynamic flux solver and demonstrate its usability for simulations in the astrophysical context. Simulations of convection for a realistic stellar profile are analyzed regarding the properties of convective boundary mixing. The time-implicit Seven-League Hydro (SLH) code was used to perform multidimensional simulations of convective helium shell burning based on a 25 M ⊙ star model. The results obtained with the low-Mach AUSM + -up solver were compared to results when using its non low-Mach variant AUSM B + -up. We applied well-balancing of the gravitational source term to maintain the initial hydrostatic background stratification. The computational grids have resolutions ranging from 180 × 90 2 to 810 × 540 2 cells and the nuclear energy release was boosted by factors of 3 × 10 3 , 1 × 10 4 , and 3 × 10 4 to study the dependence of the results on these parameters.

79 ASTRONOMY AND ASTROPHYSICS↗