Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamic codes”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Achieving Higher Order Accuracy in Space in Hydrodynamic Simulations of Self-Gravitating Gas

Modern astrophysical simulation codes employ a variety of numerical algorithms capable of achieving higher-order accuracy in both space and time. Albeit they succeed in achieving an effective higher spatial resolution and in suppressing the numerical damping of waves, to our knowledge, all current astrophysical simulations invoking self-gravity are limited to second-order accuracy in space. If we can devise an algorithm to evaluate self-gravity with a higher-order spatial accuracy, we can better the evaluation of the gravitational acceleration and gravitational energy release which dictate the evolution of many astrophysical systems. Herein, we present a numerical algorithm for self-gravitating hydrodynamics capable of achieving fourth order accuracy for a given density distribution on a Cartesian uniform grid. First, we derive the cell-averaged gravitational potential at fourth-order accuracy from the cell-averaged density by solving the Poisson equation. Next, we obtain the cell average of the product of the density and gravitational acceleration, which differs from the cell-averaged density multiplied by the cell-averaged gravitational acceleration. We then show the verification of the algorithm by applying it to critical test problems: (1) maintaining equilibria of self-gravitating slabs, even upon advection, (2) evolving a polytropic sphere with a massive power-law envelope, and (3) conservation of specific entropy during the propagation of a sound wave.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Physics of the implosion up until the time of ignition in a revolver (triple-shell) capsule

Here, a revolver is designed to have two dynamical fuel implosion stages to reach ignition: a shock phase in which the fuel is pre-heated up until the shock collapses at the capsule center, followed by an adiabatic compression phase to ignition. This picture assumes ideal hydrodynamics. Employing the state-of-the-art, hybrid (kinetic-ion/fluid electron), multi-ion Vlasov-Fokker-Planck code, iFP, along with semi-analytic predictions from ideal hydrodynamics in spherical geometry, we confirm this two phase picture—resolving some outstanding questions in Revolver design along the way. Although we find that shock kinetic effects and non-ideality are present in the course of the implosion, these effects do not change the overall dynamics (which is well described by ideal hydrodynamics theory). Additionally, we put the assumption of adiabaticity, for the post-collapse phase, to the test. Finally, we verify the ideal hydro-conjecture that the shock collapse (fuel) convergence is ~2.12 given a single shock, and we confirm that the convergence is about 2.5 in the case of two shocks—in agreement with previous radiation-hydrodynamics simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Three-dimensional Core-collapse Supernova Simulations with 160 Isotopic Species Evolved to Shock Breakout

We present three-dimensional simulations of core-collapse supernovae using the FLASH code that follow the progression of the explosion to the stellar surface, starting from neutrino radiation hydrodynamic simulations of the neutrino-driven phase performed with the Chimera code. We consider a 9.6 M ⊙ zero-metallicity progenitor starting from both 2D and 3D Chimera models and a 10 M ⊙ solar-metallicity progenitor starting from a 2D Chimera model, all simulated until shock breakout in 3D while tracking 160 nuclear species. The relative velocity difference between the supernova shock and the metal-rich Rayleigh–Taylor (R-T) "bullets" determines how the metal-rich ejecta evolves as it propagates through the density profile of the progenitor and dictates the final morphology of the explosion. We find maximum 56 Ni velocities of ~1950 and ~1750 km s -1 at shock breakout from 2D and 3D 9.6 M ⊙ Chimera models, respectively, due to the bullets' ability to penetrate the He/H shell. When mapping from 2D, we find that the development of higher-velocity structures is suppressed when the 2D Chimera model and 3D FLASH model meshes are aligned. The development of faster-growing spherical-bubble structures, as opposed to the slower-growing toroidal structure imposed by axisymmetry, allows for interaction of the bullets with the shock and seeds further R-T instabilities at the He/H interface. We see similar effects in the 10 M ⊙ model, which achieves maximum 56 Ni velocities of ~2500 km s -1 at shock breakout.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2020

Synthetic Observations of Helioseismic and Magnetic Imager

Synthetic observations of the FeI line (6173A), which corresponds to the operational wavelength of Helioseismic and Magnetic Imager (HMI) on board of the NASA Solar Dynamics Observatory (SDO). Synthetic observations represent 4D maps for four angular distances from the solar disk center (0deg), 30deg, 45deg, and 60deg at 30degrees latitude in range +/- 0.75 Angstrom from the reference line of 6173.341 Angstrom, number of points across the spectral line is 256, spatial resolution of 100km, and temporal resolution of seconds. Synthetic data set obtained using SPINOR code (Frutiger et al. 2000). Synthetic observations computed from hydrodynamic simulations, which cover 25Mm of upper layers of the convection zone and 1Mm of atmosphere, the horizontal size of the domain of 80x80 Mm. Simulation results obtained the “StellarBox” code (Wray et al., 2015, 2018) with imposed rotation rate at 30deg latitude using F-plane approximation.

Irina Kitiashvili↗
DOE OSTI2021

Exascale models of stellar explosions: Quintessential multi-physics simulation

The ExaStar project aims to deliver an efficient, versatile, and portable software ecosystem for multi-physics astrophysics simulations run on exascale machines. The code suite is a component-based multi-physics toolkit, built on the capabilities of current simulation codes (in particular Flash-X and Castro), and based on the massively parallel adaptive mesh refinement framework AMReX. It includes modules for hydrodynamics, advanced radiation transport, thermonuclear kinetics, and nuclear microphysics. The code will reach exascale efficiency by building upon current multi- and many-core packages integrated into an orchestration system that uses a combination of configuration tools, code translators, and a domain-specific asynchronous runtime to manage performance across a range of platform architectures. The target science includes multi-physics simulations of astrophysical explosions (such as supernovae and neutron star mergers) to understand the cosmic origin of the elements and the fundamental physics of matter and neutrinos under extreme conditions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Dust-void evolution driven by turbulent dust flux can induce runaway migration of Earth-mass planets

Torques from asymmetric dust structures (so-called dust-void and filamentary structures) formed around low-mass planets embedded in a nonturbulent dust-gas disk can exceed the torques produced by the gas disk component and then go on to dominate the planet’s orbital dynamics. Here, we investigate how these structures (hence the dust torque) change when the effect of turbulent dust diffusion and dust feedback are included, along with the direct implications on the migration of Earth-like planets. Using the FARGO3D code, we performed 2D and 3D multifluid hydrodynamic simulations, focusing on a non-migrating planet with a mass of M p = 1.5 M ⊕ in 2D and on migrating planets with M p ∈ [1.5, 12] M ⊕ in 3D. We varied the δ-dimensionless diffusivity parameter in the range [0, 3 × 10 −3 ] and considered three different Stokes numbers, St = {0.04, 0.26, 0.55}, which are representative of the gas-dominated, the transitional, and the gravity-dominated regimes, respectively. In our 2D models, we find that turbulent diffusion of dust prevents the formation of the dust-void and filamentary structures when δ > 3 × 10 −4 . Otherwise, dust structures survive turbulent diffusion flow. However, dust and total torques become positive only in transitional and gravity-dominated regimes. In our 3D models, we find that the dust-void is drastically modified and the high-density ring-shaped barrier delineating the dust-void disappears if δ ≳ 10 −4 , due to the effect of dust turbulent diffusion along with the back-reaction of the dust. For all values of δ, the filament in front of the planet is replaced by a low-density trench. Remarkably, as we have allowed the planets to migrate, the evolving dust-void can drive either runaway migration or outward (inward) oscillatory-torque migration. Our study thus suggests that low-mass Earth-like planets can undergo runaway migration in dusty disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Tidal disruption discs formed and fed by stream–stream and stream–disc interactions in global GRHD simulations

When a star passes close to a supermassive black hole (BH), the BH’s tidal forces rip it apart into a thin stream, leading to a tidal disruption event (TDE). In this work, we study the post-disruption phase of TDEs in general relativistic hydrodynamics (GRHD) using our GPU-accelerated code h-amr. We carry out the first grid-based simulation of a deep-penetration TDE (β = 7) with realistic system parameters: a black hole-to-star mass ratio of 10 6 , a parabolic stellar trajectory, and a non-zero BH spin. We also carry out a simulation of a tilted TDE whose stellar orbit is inclined relative to the BH midplane. We show that for our aligned TDE, an accretion disc forms due to the dissipation of orbital energy with ~20 percent of the infalling material reaching the BH. The dissipation is initially dominated by violent self-intersections and later by stream–disc interactions near the pericentre. The self-intersections completely disrupt the incoming stream, resulting in five distinct self-intersection events separated by approximately 12 h and a flaring in the accretion rate. We also find that the disc is eccentric with mean eccentricity e ≈ 0.88. For our tilted TDE, we find only partial self-intersections due to nodal precession near pericentre. Although these partial intersections eject gas out of the orbital plane, an accretion disc still forms with a similar accreted fraction of the material to the aligned case. These results have important implications for disc formation in realistic tidal disruptions. For instance, the periodicity in accretion rate induced by the complete stream disruption may explain the flaring events from Swift J1644+57.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Halted-pendulum Relaxation: Application to White Dwarf Binary Initial Data

Studying compact-star binaries and their mergers is integral to determining progenitors for observable transients. Today, compact-star mergers are typically studied via state-of-the-art computational fluid dynamics codes. One such numerical technique, smoothed particle hydrodynamics (SPH), is frequently chosen for its excellent mass, energy, and momentum conservation. The natural treatment of vacuum and the ability to represent highly irregular morphologies make SPH an excellent tool for the study of compact-star binaries and mergers. For many scenarios, including binary systems, the outcome of simulations is only as accurate as the initial conditions. For SPH, it is essential to ensure that the particles are distributed regularly, representing the initial density profile but without long-range correlations. Particle noise in the form of high-frequency local motion and low-frequency global dynamics must be damped out. Damping the latter can be as computationally intensive as the actual simulation. We discuss a new and straightforward relaxation method, halted-pendulum relaxation (HPR), to remove global oscillation modes of SPH particle configurations. In combination with effective external potentials representing gravitational and orbital forces, we show that HPR has an excellent performance in efficiently relaxing SPH particles to the desired density distribution and removing global oscillation modes. We compare the method to frequently used relaxation approaches and test it on a white dwarf binary model at its Roche-lobe overflow limit. We highlight the importance of our method in achieving accurate initial conditions and its effect on achieving circular orbits and realistic accretion rates when compared with other general relaxation methods.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1976

Impact of composite plates: Analysis of stresses and forces

The foreign object damage resistance of composite fan blades was studied. Edge impact stresses in an anisotropic plate were first calculated incorporating a constrained layer damping model. It is shown that a very thin damping layer can dramatically decrease the maximum normal impact stresses. A multilayer model of a composite plate is then presented which allows computation of the interlaminar normal and shear stresses. Results are presented for the stresses due to a line impact load normal to the plane of a composite plate. It is shown that significant interlaminar tensile stresses can develop during impact. A computer code was developed for this problem using the fast Fourier transform. A marker and cell computer code were also used to investigate the hydrodynamic impact of a fluid slug against a wall or turbine blade. Application of fluid modeling of bird impact is reviewed.

Moon, F. C.↗
DOE OSTI2023

Multidimensional Radiation Hydrodynamics Simulations of Pulsational Pair-instability Supernovae

Stars with masses of 80–130 M ⊙ can encounter pulsational pair-instability at the end of their lives, which triggers consecutive episodes of explosive burning that eject multiple massive shells. Collisions between these shells produce bright transients known as pulsational pair-instability supernovae (PPI SNe) that may explain some extreme supernovae. In this paper, we present the first 2D and 3D radiation hydrodynamics simulations of PPI SNe with the CASTRO code. Radiative cooling causes the collided shells to evolve into thin, dense structures with hot spots that can enhance the peak luminosity of the SN by factors of 2–3. The light curve peaks at 1.9–2.1 × 10 43 erg s -1 for 50 days and then plateaus at 2–3 × 10 42 erg s -1 for 200 days, depending on the viewing angle. The presence of 12 C and 16 O and the absence of 28 Si and 56 Fe in its spectra can uniquely identify this transient as a PPI SN in follow-up observations. Our models suggest that multidimensional radiation hydrodynamics is required to model the evolution and light curves of all shell-collision SNe, such as Type IIne, not just PPI SNe.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗
DOE OSTI2022

General Relativistic Implicit Monte Carlo Radiation-hydrodynamics

Abstract We report on a new capability added to our general relativistic radiation-magnetohydrodynamics code, Cosmos++ : an implicit Monte Carlo (IMC) treatment for radiation transport. The method is based on a Fleck-type implicit discretization of the radiation-hydrodynamics equations, but generalized for both Newtonian and relativistic regimes. A multiple reference frame approach is used to geodesically transport photon packets (and solve the hydrodynamics equations) in the coordinate frame, while radiation–matter interactions are handled either in the fluid or electron frames then communicated via Lorentz boosts and orthonormal tetrad bases attached to the fluid. We describe a method for constructing estimators of radiation moments using path-weighting that generalizes to arbitrary coordinate systems in flat or curved spacetime. Absorption, emission, scattering, and relativistic Comptonization are among the matter interactions considered in this report. We discuss our formulations and numerical methods, and validate our models against a suite of radiation and coupled radiation-hydrodynamics test problems in both flat and curved spacetimes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

OC6 Phase Ia - Nonlinear hydrodynamic loading validation dataset

Two validation campaigns were examined within the Offshore Code Comparison Collaboration, Continued, with Correlation and unCertainty (OC6) Phase 1 project to examine the modeling tools' underprediction of loads and motion of a floating wind semisubmersible (semi) at their surge and pitch natural frequencies. These campaigns were performed at the Maritime Research Institute Netherlands (MARIN) in 2017 and 2018. The load cases (LC) considered include: LC1 – Load measurements across semi under current loading; LC2 - Load measurements across semi under forced surge oscillation; LC3 – Load measurements across semi under wave loading, while held fixed; LC4 – Free-decay motion measurements in surge, pitch, and heave; and LC5 – Motion measurements under wave loading. Details on the results from the OC6 Phase Ia project can be found in the reference, “OC6 Phase 1: Investigating the underprediction of low-frequency hydrodynamic loads and responses of floating wind turbines”, J Phys: Conf Series 1618 032033.

17 WIND ENERGY↗
DOE OSTI2020

Theoretical Review of the MFIX Fluid and Two-Fluid Models

MFIX (Multiphase Flow with Interphase eXchanges) is a general purpose code that can be used for describing the hydrodynamics, chemical reactions, and heat transfer of dense or dilute multiphase fluid flows. MFIX calculations give detailed information on pressure, temperature, composition, and velocity distributions in the defined system. This report provides an overview of MFIX’s expressions for single phase flow and continuous dispersed multiphase flow (i.e., the two-fluid model). In addition to the two fluid model (MFIX-TFM), MFIX contains a discrete element model (MFIX-DEM) and a multi-phase particle in cell model (MFIX-PIC) that are not reviewed within this document.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1983

Analytical calculation of a single jet in cross-flow and comparison with experiment

A series of progressively finer grid systems to calculate the single jet in cross-flow experimentally measured were employed. The experimental measurements provide a collection of velocities, turbulence intensities, and jet concentration profiles with measurements of the inlet field. The use of a series of progressively finer grid systems allows a differentiation between numerical errors and the hydrodynamic modeling assumptions embodied in the 3D combustor code.

Claus, R. W.↗
NASA NTRS2012

X-Ray Spectra from MHD Simulations of Accreting Black Holes

We present the results of a new global radiation transport code coupled to a general relativistic magneto-hydrodynamic simulation of an accreting, nonrotating black hole. For the first time, we are able to explain from first principles in a self-consistent way the X-ray spectra observed from stellar-mass black holes, including a thermal peak, Compton reflection hump, power-law tail, and broad iron line. Varying only the mass accretion rate, we are able to reproduce the low/hard, steep power-law, and thermal-dominant states seen in most galactic black hole sources. The temperature in the corona is T(sub e) ~ 10 keV in a boundary layer near the disk and rises smoothly to T(sub e) greater than or approximately 100 keV in low-density regions far above the disk. Even as the disk's reflection edge varies from the horizon out to approximately equal to 6M as the accretion rate decreases, we find that the shape of the Fe Ka line is remarkably constant. This is because photons emitted from the plunging region are strongly beamed into the horizon and never reach the observer. We have also carried out a basic timing analysis of the spectra and find that the fractional variability increases with photon energy and viewer inclination angle, consistent with the coronal hot spot model for X-ray fluctuations.

Schnittman, Jeremy D.↗
NASA NTRS2019

SI IV Resonance Line Emission during Solar Flares: Non-LTE, Non-equilibrium, Radiation Transfer Simulations

The Interface Region Imaging Spectrograph routinely observes the Si IV resonance lines. When analyzing quiescent observations of these lines, it has typically been assumed that they form under optically thin conditions. This is likely valid for the quiescent Sun, but this assumption has also been applied to the more extreme flaring scenario. We used 36 electron-beam-driven radiation hydrodynamic solar flare simulations, computed using the RADYN code, to probe the validity of this assumption. Using these simulated atmospheres, we solved the radiation transfer equations to obtain the non-LTE, nonequilibrium populations, line profiles, and opacities for a model silicon atom, including charge exchange processes. This was achieved using the "minority species" version of RADYN. The inclusion of charge exchange resulted in a substantial fraction of Si IV at cooler temperatures than those predicted by ionization equilibrium. All simulations with an injected energy flux F> 5 x 10(exp 10) erg cm(exp -2) s(exp -1) resulted in optical depth effects on the Si IV emission, with differences in both intensity and line shape compared to the optically thin calculation. Weaker flares (down to F ~ 10(exp 9) erg cm(exp -2) s(exp -1)) also resulted in Si IV emission forming under optically thick conditions, depending on the other beam parameters. When opacity was significant, the atmospheres generally had column masses in excess of 5 x 10-6 g cm-2 over the temperature range 40-100 kK, and the Si IV formation temperatures were between 30 and 60 kK. We urge caution when analyzing Si IV flare observations, or when computing synthetic emission without performing a full radiation transfer calculation.

Kerr, Graham S.↗
DOE OSTI2024

A modern concept of Lagrangian hydrodynamics

Here, we offer a modern interpretation of Lagrangian hydrodynamics as employed in Lagrangian simulations of compressible fluid flow. Our main result is to show that artificial viscosity, traditionally viewed as a numerical artifice to control unphysical oscillations in flows with shocks, actually represents a physical process and is necessary to derive accurate simulations in any compressible flow. We begin by reviewing the origins of two numerical devices, artificial viscosity and finite-volume methods. We proceed to construct a mathematical (PDE) model that incorporates those numerics and in which a new length scale, the observer, arises representing the discretization. Associated with that length scale, there are new inviscid fluxes that are the artificial viscosity as first formulated by Richtmyer and an artificial heat flux postulated by Noh but typically not included in Lagrangian codes. We discuss the connection of our results to bivelocity hydrodynamics. We conclude with some speculation as to the direction of future developments in multidimensional Lagrangian codes as computers get faster and have larger memories.

97 MATHEMATICS AND COMPUTING↗