Search NASA⌕ Search

SEARCH · Search NASA

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 19 records

Flare loop radiative hydrodynamics. III - Nonlocal radiative transfer effects

The study has three goals. The first is to demonstrate that processes exist whose intrinsic nonlocal nature cannot be represented by local approximations. The second is to elucidate the physical nature and origins of these nonlocal processes. The third is to suggest that the methods and results described here may prove useful in constructing semiempirical models of the chromosphere by means more efficient than trial and error. Matrices are computed that describe the effect of a temperature perturbation at an arbitrary point in the loop on density, hydrogen ionized fraction, total radiative loss rate, and radiative loss rate of selected hydrogen lines and continua at all other points. It is found that the dominant nonlocal radiative transfer effects can be separated into flux divergence coefficient effects and upper level population effects. The former are most important when the perturbation takes place in a region of significant opacity. Upper level population effects arise in both optically thick and thin regions in response to nonlocal density, ionization, and interlocking effects.

Canfield, R. C.↗

Test problems for radiation and radiation-hydrodynamics codes

A number of test problems for radiation and radiation-hydrodynamics computer codes are described. These include evolution to radiative equilibrium, cooling from radiative equilibrium, subcritical and supercritical radiating shocks, and a radiating blast wave in a power-law density distribution. For each test problem, example input parameters and plots of the results are presented. Some test problems for pure hydrodynamics are also suggested. The radiation-hydrodynamics code used to perform the example test problems and the equations it solves are described in some detail.

Ensman, Lisa↗

Improved Finite-Volume Method for Radiative Hydrodynamics

Fully coupled simulations of hydrodynamics and radiative transfer are essential to a number of fields ranging from astrophysics to engineering applications. Of particular interest in this work are hypersonic atmospheric entries and associated experimental apparatus, e.g., shock tubes and high enthalpy testing facilities. The radiative transfer calculations must supply to the CFD a heating term in the energy equation in the form of the divergence of the radiative heat flux and the radiative heat fluxes to bounding surfaces. It is most efficient to solve the radiative transfer equation on the same grid as the CFD solution, and this work presents an algorithm with improved accuracy for such simulations on structured and unstructured grids compared to more conventional approaches. Results will be shown for shock radiation during hypersonic reentry. Issues of parallelization within a radiation sweep will also be discussed.

Wray, Alan↗

AGN Obscuration Through Dusty Infrared Dominated Flows: Radiation-Hydrodynamics Solution for the Wind - 1

We construct a radiation-hydrodynamics model for the obscuring toroidal structure in active galactic nuclei. In this model the obscuration is produced at parsec scale by a dense, dusty wind which is supported by infrared radiation pressure on dust grains. To find the distribution of radiation pressure, we numerically solve the 2D radiation transfer problem in a flux limited diffusion approximation. We iteratively couple the solution with calculations of stationary 1D models for the wind, and obtain the z-component of the velocity. Our results demonstrate that for AGN luminosities greater than 0.1 L(sub edd) external illumination can support a geometrically thick obscuration via outflows driven by infrared radiation pressure. The terminal velocity of marginally Compton-thin models (0.2 < tau(sub T) < 0.6), is comparable to or greater than the escape velocity. In Compton thick models the maximum value of the vertical component of the velocity is lower than the escape velocity, suggesting that a significant part of our torus is in the form of failed wind. The results demonstrate that obscuration via normal or failed infrared-driven winds is a viable option for the AGN torus problem and AGN unification models. Such winds can also provide an important channel for AGN feedback.

Dorodnitsyn, A.↗

Radiative-hydrodynamical simulations of X-ray burst-induced accretion disk coronae

A fully self-consistent radiative hydrodynamical algorithm developed earlier is used to study the 'steady' state coronae in order to simulate the response of the evaporated outflow to a time-dependent irradiation of the disk's surface. It is found that the burst definition is altered at both low and high inclination angles i. At i greater than or approximately equal to 45 deg, the detectable flux is enhanced (as much as 50 percent) by X-rays scattered into the line of sight from the radiation propagation through the corona, and (greater than or approximately equal to 30 percent) by X-rays reflected off the corona when i is less than or approximately equal to 45 deg.

Mella, Fulvio↗

AGN Obscuration Through Dusty Infrared Dominated Flows. II. Multidimensional, Radiation-Hydrodynamics Modeling

We explore a detailed model in which the active galactic nucleus (AGN) obscuration results from the extinction of AGN radiation in a global ow driven by the pressure of infrared radiation on dust grains. We assume that external illumination by UV and soft X-rays of the dusty gas located at approximately 1pc away from the supermassive black hole is followed by a conversion of such radiation into IR. Using 2.5D, time-dependent radiation hydrodynamics simulations in a ux-limited di usion approximation we nd that the external illumination can support a geometrically thick obscuration via out ows driven by infrared radiation pressure in AGN with luminosities greater than 0:05 L(sub edd) and Compton optical depth, Tau(sub T) approx > & 1.

Dorodnitsyn, Anton↗

Radiative-hydrodynamical simulations of accretion disk coronae

An algorithm has been developed that will carry out detailed, two-dimensional, fully self-consistent radiative-hydrodynamical simulations of accretion disk coronae in X-ray-luminous compact sources. The calculation reported here, for an accreting neutron star radiating at 0.5 times the Eddington luminosity, reveals several striking features. (1) The corona is comprised of two main regions - an inner highly dynamic portion whose vertical structure varies cyclically on a dynamical time scale, and an outer, more stable zone in which the evaporated plasma rises to form a 'sheath' that gradually merges into a wind at large radii. (2) The coronal structure shows a density inversion that contrasts sharply with the Gaussian profiles of earlier hydrostatic models. (3) Interestingly, flow velocities as high as a few billion cm/s are not uncommon in portions of the corona. The relevance of this point to the large measured width of the emission features seen in low-mass X-ray binaries and the Galactic black hole candidates is discussed.

Melia, Fulvio↗

3D Radiative Hydrodynamics Modeling of Convection of Stars to Probe Their Interiors and Photospheric Properties

The dramatic flow of data from the Kepler and K2 missions opens the opportunity to significantly improve our knowledge of stellar interiors, surface dynamics, and structure. However, interpretation of these observations is a challenging task because it depends on tiny effects that can be studied only with advanced first-principles modeling. We present results of 3D time-dependent radiative hydrodynamic simulations of stellar outer convection zones and atmospheres taking into account chemical composition, radiative transfer, turbulence effects, and a realistic equation of state for main sequence stars. We will discuss properties of convective structure and dynamics, convective overshoot, effects of magnetic fields and rotation, as well as the potential influence of turbulent surface dynamics on high-precision RV measurements.

Photospheric↗

Data-Driven Radiative Hydrodynamic Modeling of the 2014 March 29 X1.0 Solar Flare

Spectroscopic observations of solar flares provide critical diagnostics of the physical conditions in the flaring atmosphere. Some key features in observed spectra have not yet been accounted for in existing flare models. Here we report a data-driven simulation of the well-observed X1.0 flare on 2014 March 29 that can reconcile some well-known spectral discrepancies. We analyzed spectra of the flaring region from the Interface Region Imaging Spectrograph (IRIS) in Mg II hk, the Interferometric BIdimensional Spectropolarimeter at the Dunn Solar Telescope (DSTIBIS) in H(alpha) 6563A and Ca II 8542A, and the Reuven Ramaty High Energy Solar Spectroscope Imager (RHESSI) in hard X-rays. We constructed a multithreaded flare loop model and used the electron flux inferred from RHESSI data as the input to the radiative hydrodynamic code RADYN to simulate the atmospheric response. We then synthesized various chromospheric emission lines and compared them with the IRIS and IBIS observations. In general, the synthetic intensities agree with the observed ones, especially near the northern footpoint of the flare. The simulated Mg II line profile has narrower wings than the observed one. This discrepancy can be reduced by using a higher microturbulent velocity (27 km/s) in a narrow chromospheric layer. In addition, we found that an increase of electron density in the upper chromosphere within a narrow height range of approx. 800 km below the transition region can turn the simulated Mg II line core into emission and thus reproduce the single peaked profile, which is a common feature in all IRIS flares.

Da Costa, Fatima Rubio↗

Flare loop radiative hydrodynamics. I - Basic methods

The study presented here has two goals. The first is to examine the role of energy transport mechanisms in flare dynamics, for instance, thermal cone conduction in the corona and radiative transfer in the chromosphere. The second is to provide diagnostics of flare energization processes by predicting the evolutionary behavior of concurrent emissions from all regions of the flaring atmosphere. To attain these ends, a numerical method is developed for the simultaneous solution of thy continuity, momentum, and energy equations; the time-dependent atomic rate equations describing ionization and excitation; and the radiative transfer equations. With the physical model and computational methods used here, all components of the plasma move as a single fluid and in the initial studies are assumed to have a common kinetic temperature. The plasma is constrained to move in one dimension along the axis of a loop of nonuniform cross section, the geometry of which is defined by a sufficiently strong magnetic field.

Mcclymont, A. N.↗

Flare loop radiative hydrodynamics. IV - Dynamic evolution of unstable semiempirical loop models

The evolution of the unstable solar atmosphere into the nonlinear phase, in response to various perturbations, is followed. The initial dynamic evolution of the atmosphere follows the predictions of linear stability analysis. In the nonlinear phase, rapid changes are confined to the transition region; these changes are manifested as a propagation of the transition region through the plasma, i.e., chromospheric evaporation or condensation. Global evolution therefore proceeds on the coronal conductive time scale. The rate of propagation of the transition region is determined by the imbalance between the energy supplied by thermal conduction from the corona and radiative cooling within the transition region itself. Flow velocities in the lower corona during evaporation or condensation are, in the cases studied, of order 3 km/s. The observed dynamic evolution is consistent with the existence of relatively long-lived coronal loops whose brightnesses vary on the evaporative time scale.

An, C.-H.↗

Towards Surrogate Modeling of Subgrid Turbulent Transport for 3D Radiative Hydrodynamic Simulations of the Quiet Sun

In this work, we investigate the use of deep learn-ing techniques as surrogate models, to enhance the estimationof effects of subgrid turbulent transport for 3D radiatuve hy-drodynamic simulations of the quiet Sun. We develop two dis-tinct 3D Convolutional Neural Networks (3DCNNs) to capturespatio-temporal dependencies in 3D velocity fields, leveragingdifferent activation functions and architectural designs. Thesemodels integrate both averaged velocity vector components andscalar features such as plasma density to enhance predictionaccuracy. Additionally, a Multilayer Perceptron (MLP) modelis employed to approximate complex nonlinear relationships,offering a comparison in performance between convolutionaland fully connected architectures. Logarithmic transformationis applied to the targets to handle heavily skewed data, im-proving model performance. All models are compared againsta physics-based Gradient Model. Results show that the 3DCNNmodels excel at approximating Reynolds stress tensors, makingthem a candidate for assisting in producing reduced resolutionsimulations, and thereby reducing computational overheadwhile maintaining higher accuracy than the baseline. Thesefindings demonstrate the potential of deep learning, particu-larly CNNs, to advance scalable and accurate simulations ofsolar dynamics, offering a promising alternative to traditionalturbulence models.

Heliophysics↗

Surrogate Modeling of Subgrid Turbulent Transport Based on 3D Radiative Hydrodynamic Simulations of the Quiet Sun

Turbulent Transport: Plays a critical role in astrophysical plasmas, such as the solar interior, spanning multiple scales and challenging traditional modeling approaches. Objective: Develop machine learning (ML) models—MLP and CNN—to predict subgrid Reynolds stress tensors from StellarBox 3D simulations of the solar atmosphere. Benchmarking: Compare ML-driven models against physics-based Gradient and Smagorinsky approaches.

SMD↗