Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “radiation hydrodynamics”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

The F-CHROMA Grid of 1D RADYN Flare Models

Context. Solar flares are the result of the sudden release of magnetic energy in the corona. Much of this energy goes into accelerating charged particles to high velocity. These particles travel along the magnetic field and the energy is dissipated when the density gets high enough, primarily in the solar chromosphere. Modelling this region is difficult because the radiation energy balance is dominated by strong, optically thick spectral lines. Aims. Our aim is to provide the community with realistic simulations of a flaring loop with an emphasis on the detailed treatment of the chromospheric energy balance. This will enable a detailed comparison of existing and upcoming observations with synthetic observables from the simulations, thereby elucidating the complex interactions in a flaring chromosphere. Methods. We used the 1D radiation hydrodynamics code RADYN to perform simulations of the effect of a beam of electrons injected at the apex of a solar coronal loop. A grid of models was produced, varying the total energy input, the steepness, and low-energy cutoff of the beam energy spectrum. Results. The full simulation results for a grid of models are made available online. Some general properties of the simulations are discussed.

hydrodynamics↗
NASA NTRSDate not supplied

Comparison of Ring Diagrams Based on the Doppler Shifts of Synthetic Data Obtained with Bisector Method and SDO/HMI Pipeline

Ring diagrams are cross-sections of three-dimensional spatiotemporal power spectra of solar oscillations. The rings reveal information about sub-surface flows and represent an important tool for helioseismology. Ring diagrams can be constructed using Doppler velocity or intensity maps of spectral lines. How the velocities are computed is an important factor for accuracy of information we can retrieve on subsurface flows. In our work, the ring diagrams are generated from Doppler shift data of synthesized Fe I 6173 Å line. We compare ring diagrams computed by two methods–HMI line-of-sight pipeline and the bisector of Fe line. Fe I line is synthesized for StellarBox 3D Radiative hydrodynamic simulations under LTE assumption. We aim to answer the following questions: 1.How do power spectra obtained from velocities computed with the HMI pipeline and bisector of the Fe I 6173Å compare? 2.How do the power spectral density retrieved with each method vary with heliocentric angle? 3.What is the effect of changing resolution on the power spectral density in ring diagrams obtained with the two methods?

SMD↗
DOE OSTI2020

Revisiting Theoretical Limits for One-Degree-of-Freedom Wave Energy Converters: Preprint

This work revisits the theoretical limits of one-degree of freedom wave energy converters (WECs). This paper considers the floating sphere used in the Ocean Energy Systems (OES) Task 10 Wave Energy Converter modeling and verification effort for analysis. Analytical equations are derived to determine bounds on the motion amplitude, time-averaged power, and power take-off (PTO) force. A unique result was found that shows the time-averaged power absorbed by a WEC can be defined solely by the inertial properties and the radiation hydrodynamic coefficients. In addition, a unique expression for the PTO force amplitude was derived which has provided upper and lower bounds when resistive control is used to maximize power generation. For complex conjugate control, this same expression can only provide a lower bound as there is theoretically no upper bound. These bounds are used to compare the performance of a floating sphere if it were to extract energy using surge or heave motion. The analysis shows that because of the differences in hydrodynamic coefficients of each oscillating mode there will be different frequency ranges that provide better power capture efficiency. The influence of a motion constraint on power absorption while also utilizing a nonideal power take-off is examined and found to reduce the losses associated with bidirectional energy flow. The expression to calculate the time-averaged power with a nonideal PTO is modified by the mechanical-to-electrical efficiency and the ratio of the PTO spring and damping coefficients. These were separated in the expression which allows for limits to be set on the possible values of PTO coefficients to ensure a net flow of power to the grid.

50 EE - Wind and Water Power Program - Water (EE-4↗
DOE OSTI2020

Challenges in simulating ground interacting nuclear explosions

This paper summarizes recent above-ground nuclear explosion simulations as part of a broader effort to better characterize conditions within a fireball that may influence the chemical evolution of bomb materials and other materials entrained from the local explosion environment. A critical component of this work is validation against historic footage of atmospheric testing, requiring that we understand how the frequency-dependent sensitivity of the utilized film footage influences data captured in such images. We focus first on the early physics of a nuclear explosion in the atmosphere before discussing some of the technical challenges we seek to capture in late-time models that include more complex emplacement conditions and subsurface features. We discuss required physics packages (compressible hydrodynamics, radiation transport, as well as necessary ancillary tables such as equations of state (EOS) and opacities). Additionally, we note reasonable “shortcuts” one may make and their limitations, e.g., using ideal gas EOS, replacing spectrally resolved radiation with spectrally averaged radiation, and exchanging deterministic transport with diffusion. We then discuss an approach to achieving an equilibrated initial stress state for problems where buoyancy and subsurface lithostatic stress are important. Our methodology is presented in the context of LLNL’s ALE3D multiphysics code but may readily be implemented in other codes. In this paper, we start with a description of the challenges of NUDET simulations, followed by a presentation of the simulated intensity (flux) as it would appear on an analysis of the Dixie test. We then progressively introduce additional complexity in subsequent sections (near-surface burst and gravity initialization) before concluding.

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

thornado-transport: Anderson- and GPU-accelerated nonlinear solvers for neutrino-matter coupling

Algorithms for neutrino-matter coupling in core-collapse supernovae (CCSNe) are investigated in the context of a spectral two-moment model, which is discretized in space with the discontinuous Galerkin method, integrated in time with implicit-explicit (IMEX) methods, and implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado). The model considers electron neutrinos and antineutrinos and tabulated opacities from Bruenn (1985), which includes neutrino-electron scattering and pair processes. The nonlinear system arising from implicit time discretization of the equations governing neutrino-matter coupling is iterated to convergence using Anderson-accelerated fixed-point methods, which avoid formation of Jacobians and inversion of dense linear systems. Numerical experiments show that, for a given tolerance, a nested iteration scheme which aims to reduce opacity evaluations can lower the computational cost. Our initial port to GPUs, using both OpenMP and OpenACC, shows an overall speedup of up to ~ 100× when compared to results using a single CPU core. These results indicate that the algorithms implemented in thornado are well-suited to GPU acceleration.

Laiu, Paul↗
DOE OSTI2026

The Challenge of Observing Patchy Reionization with CMB Optical-Depth Fluctuations

Spatial fluctuations in the Thomson optical depth encode information about the inhomogeneous nature of cosmic reionization. We compute the optical-depth angular power spectrum, $C_\ell^{ττ}$, using past lightcones constructed from five Cosmic Reionization on Computers (CROC) radiation-hydrodynamical simulations. By decomposing the electron-density field into patchy and density components, we quantify the separate contributions of ionization-fraction and baryon-density fluctuations to the optical-depth anisotropy. Because the simulations end at $z\approx5$, we supplement the reionization-era signal with an analytic estimate of the fully ionized low-redshift contribution. We find that baryon-density fluctuations dominate the high-redshift signal over most angular scales, while the accumulated low-redshift contribution exceeds the high-redshift signal across the full multipole range considered. Our results demonstrate that a significant fraction of the optical-depth power is not uniquely associated with reionization morphology, implying that future interpretations of $C_\ell^{ττ}$ must account for the density contribution in addition to patchy ionization.

Takoudes, Nick [Chicago U., Astron. Astrophys. Ctr↗
DOE OSTI2023

Modeling Hydrodynamic Instabilities, Shocks, and Radiation Waves in High Energy Density Experiments [Dissertation]

This thesis presents the computational design, modeling, and analysis of three experiments in high-energy density physics (HEDP), all of them concerning fundamental radiation flows. The first experiment is a laboratory astrophysics experiment to investigate the role of the Kelvin-Helmholtz instability (KHI) in the process of galactic filaments supplying gas to galactic halos. The achieved goal was to provide a first study in which the role of the instability is maximal and predict behavior in future iterations of an experiment accessing a more radiative regime where the role of the instability is stifled. This experiment would help answer how certain galaxies are able to grow so rapidly and produce many stars as the KHI limits this process. The second experiment, COAX, is a radiation flow experiment with a novel spectroscopy diagnostic configuration, designed to spatially measure the temperature of a radiation wave as it travels down a doped foam. A key result of this work was the development of a synthetic spectroscopy application and application of modern spectroscopy comparison techniques to provide our first temperature reconstructions from the experimental data. This experimental platform serves as the launching ground for a number of new experiments that vary the basic premise and thus is foundational to our ongoing research. The final experiment is a full integration of modeling, design, and theoretical development for the Radishock experiment. This experimental platform studies the head-on collision of a radiation wave with a counter-propagating shock, and like COAX, uses spectroscopy to diagnose and detect the interaction. My research analyzes the successful shots, indicating aspects of successful detections and suggests improvements to future iterations of the design.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

The Radiative Magneto-Hydrodynamics of Colliding Flow: Instabilities, Reconnection and Exoplanet Atmosphere Connections

Hypersonic flows of magnetized plasmas naturally produce strong shocks when obstructions or self‐interactions (flow collisions) occur. Such complex flows occur in a wide variety of HEDP settings such as Z‐pinches and laser driven implosion experiments that include magnetic fields via explosive Helmholtz coils (MIFEDS) When thermal energy behind the shocks can be lost to optically thin radiation, the flow dynamics becomes even more complex, admitting a range of new instabilities such as the Nonlinear Thin Shell Instability (NTSI). Because of their inherent non‐ linearity, multi‐dimensionality and time‐dependence (i.e. turbulence), the ability to articulate the underlying physics of colliding hypersonic magnetized flows represents a significant frontier in HEDP studies. The addition of radiative cooling to these flows propels the study even further into the domain of extreme plasma physics. Exploration of these forms of HED plasmas is of interest to basic plasma physics as well as of interest for their application to Inertial Confinement Fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Effects of hydrodynamics and thermal radiation in the atmosphere after comet impacts

Radiation phenomena in the atmosphere after impacts of cosmic bodies have special features in comparison with the surface nuclear explosions. First, initial concentration of energy after the impact is lower, and second, a wake after the passage of the meteoroid through the atmosphere has a dramatic effect on the atmospheric flow and radiation transfer. Consequently, scaling laws can not be employed for prediction of the flow in the atmosphere and the light flux on the Earth's surface. If a density of high-velocity impactor is low relative to the ground, as in a case of a comet impact on rocks, a major part of the kinetic energy is converted to internal energy of dense hot vapors. But radiation effects can be essential even for fairly low velocities of the impactor. To clarify this issue we have undertaken calculations of 100-Mt explosions at the Earth's surface caused by small comets with velocities from 10 to 70 km/sec. That is, the initial concentration of energy has been varied. The calculations have shown that for velocities of the comet greater or about 20 km/sec a portion of energy emitted from the fireball exceeds 20% of the total energy of the explosion and this quantity does not change very much with the velocity. Other aspects of this investigation are discussed.

Nemchinov, I. V.↗
NASA NTRS2023

Optimizing the Resolution of Hydrodynamic Simulations for MCRaT Radiative Transfer Calculations

Despite their discovery about half a century ago, the Gamma-ray burst (GRB) prompt emission mechanism is still not well understood. Theoretical modeling of the prompt emission has advanced considerably due to new computational tools and techniques. One such tool is the PLUTO hydro dynamics code, which is used to numerically simulate GRB outflows. PLUTO uses Adaptive Mesh Refinement to focus computational efforts on the portion of the grid that contains the simulated jet. Another tool is the Monte Carlo Radiation Transfer (MCRaT) code, which predicts electromagnetic signatures of GRBs by conducting photon scatterings within a jet using PLUTO. The effects of the underlying resolution of a PLUTO simulation with respect to MCRaT post-processing radiative transfer results have not yet been quantified. We analyze an analytic spherical outflow and a hydrodynamically simulated GRB jet with MCRaT at varying spatial and temporal resolutions and quantify how decreasing both resolutions affect the resulting mock observations. We find that changing the spatial resolution changes the hydrodynamic properties of the jet, which directly affect the MCRaT mock observable peak energies. We also find that decreasing the temporal resolution artificially decreases the high energy slope of the mock observed spectrum, which increases both the spectral peak energy and the luminosity. We show that the effects are additive when both spatial and temporal resolutions are modified. Our results allow us to understand how decreased hydrodynamic temporal and spatial resolutions affect the results of post-processing radiative transfer calculations, allowing for the optimization of hydrodynamic simulations for radiative transfer codes.

Jose Arita-Escalante↗
NASA NTRS1992

Flux-limited diffusion in hydrodynamics

The radiation transport equation, correct to order v/c, is used to derive a flux limiter valid in high-velocity hydrodynamic flows. In the limit v = 0 the widely used result of Levermore and Pomraning (LP) (1982) is recovered. In the case of homologous flow (v varies as R) the LP result remains valid, but in the presence of shock waves and other nonhomologous flows differences may be large. We calculate the Eddington factors, which are significantly altered, and discuss the importance of having them possess the correct form in the presence of a shock wave, e.g., on the neutrino transport during the formation of a supernova.

Cooperstein, J.↗
DOE OSTI2020

Discretization Writeup for Grey Flux-Limited Radiation Diffusion

This report documents the time and space discretizations for grey flux-limited diffusion applied to the thermal radiative transfer (TRT) equations. We begin with a description of the physics being solved before moving into the diffusion approximation. Once we have the TRT system, we show a finite-volume-inspired discretization from Jim Morel (Texas A&M University, NUEN 627 class notes, lecture 8). As systems become hotter they emit more photons in the form of blackbody radiation. Because average photon energy of the blackbody source is proportional to the temperature of the system, we call these thermal photons or thermal radiation. As material temperatures increase, increasing fractions of the total energy in the system go into the radiation field. In addition, radiation can deposit energy and momentum non-locally, making it an important phenomenon for heating and impulse. In order to accurately study systems at high temperatures, we wish to add the physics of thermal radiation to our hydrodynamic system. In practice, coupling radiation and hydrodynamics is often done by operator-splitting each timestep into two consecutive, non-overlapping phases: (1) update the hydrodynamics for a fixed radiation state (2) update the radiation and internal energy for an otherwise fixed hydrodynamic state. Because of this clean separation of physics updates, in this report we show only the latter phase, which involves solely the TRT equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

MOCMC: Method of Characteristics Moment Closure, a Numerical Method for Covariant Radiation Magnetohydrodynamics

In this work, we present a conservative numerical method for radiation magnetohydrodynamics with frequency-dependent full transport in stationary spacetimes. This method is stable and accurate for both large and small optical depths and radiation pressures. The radiation stress–energy tensor is evolved in flux-conservative form, and closed with a swarm of samples that each transport a multigroup representation of the invariant specific intensity along a null geodesic. In each zone, the enclosed samples are used to efficiently construct a Delaunay triangulation of the unit sphere in the comoving frame, which in turn is used to calculate the Eddington tensor, average source terms, and adaptively refine the sample swarm. Furthermore, radiation four-forces are evaluated in the moment sector in a semi-implicit fashion. The radiative transfer equation is solved in invariant form deterministically for each sample. Since each sample carries a discrete representation of the full spectrum, the cost of evaluating the transport operator is independent of the number of frequency groups, representing a significant reduction of algorithmic complexity for transport in frequency-dependent problems. The major approximation we make in this work is performing scattering in an angle-averaged way. Local adaptivity in samples also makes this scheme more amenable to nonuniform meshes than a traditional Monte Carlo method. We describe the method and present results on a suite of test problems. We find that Method of Characteristics Moment Closure converges at least as ~N -1 , rather than the canonical Monte Carlo N -1/2 , where N is the number of samples per zone.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1989

Continuous infrared emission of proto- and young-planetary nebulae

Evolutionary sequences of planetary nebulae (PNe) were calculated beginning from the moment of shell ejection to its dissipation and IR spectra of the outgoing radiation were obtained for different stages of the evolution. To solve the coupled hydrodynamical and radiation transfer problem a computer code was used in which gas and dust are treated as two separate hydrodynamical components. The models, which are calculated for two grain materials (graphite and silicate) have a size distribution of particles based on that found for the interstellar medium. In the computations, the recent optical properties of graphite and astronomical silicate grains tabulated by Draine (1987) were employed. It was assumed that silicate of graphite grains exist below some critical temperature equal to 1300 or to 1700 K, respectively. Infrared model spectra are compared with the IR radiation emitted by some proto- and some young-planetary nebulae. The observed IR continua can be quite well matched with our models with grains having a reasonable size range.

Szczerba, Ryszard↗
DOE OSTI2023

Intricate structure of the plasma Rayleigh–Taylor instability in shock tubes

Spikes and bubbles grow on unstable interfaces that are accelerated in high-energy-density conditions. If a shock propagates ahead of the interface, the plasma can be heated to extreme conditions where conduction and radiation fluxes influence the hydrodynamics. For example, a National Ignition Facility experiment found reduced single-mode nonlinear mixed-width growth in conditions scaled from a supernova explosion. We present high-resolution two-dimensional radiation hydrodynamic simulations with the Flash code that quantitatively reproduce the experiment. Radiative fluxes are primarily responsible for ablating the spike and removing the mushroom caps. The ablated plasma increases the mixed mass and forms a low-density halo with spikes forming in both directions. This is considerably more complex than classical instability. The halo is sensitive to ablative physics, so radiographing it may aid in the verification of energy transport modeling. Although ablation changes the spike shape, it has little effect on the overall mixed width for these parameters. This is because ablation enhances the bubble velocity but it has the opposite effect on the spike. The radiation transport instead suppresses the growth via increasing the shocked foam density, thus decreasing the Atwood number. Furthermore, a terminal velocity model including the rarefaction expansion agrees with the experimental mixed-width growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗