Search NASA⌕ Search

SEARCH · Search NASA

Results for “radiative transport”

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

Multigroup Scattering in Monte Carlo Radiation Transport Codes [Slides]

Legendre truncations to multigroup scattering distributions are not amenable to Monte Carlo sampling due to negative values. We have implemented two moment-preserving methods in MGMC that capture the shape of the truncation, are non-negative over [-1, 1], and can be efficiently sampled on CPUs and GPUs. MGMC can now simulate neutrons with anisotropic scattering mechanics and MGMC shows good agreement with LANL production codes PARTISN.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A finite element method for angular discretization of the radiation transport equation on spherical geodesic grids

Discrete ordinate (S N ) and filtered spherical harmonics (FP N ) based schemes have been proven to be robust and accurate in solving the Boltzmann transport equation but they have their own strengths and weaknesses in different physical scenarios. We present a new method based on a finite element approach in angle that combines the strengths of both methods and mitigates their disadvantages. The angular variables are specified on a spherical geodesic grid with functions on the sphere being represented using a finite element basis. A positivity-preserving limiting strategy is employed to prevent non-physical values from appearing in the solutions. Here, the resulting method is then compared with both S N and FP N schemes using four test problems and is found to perform well when one of the other methods fail.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2002

Space Radiation Transport Methods Development

Improved spacecraft shield design requires early entry of radiation constraints into the design process to maximize performance and minimize costs. As a result, we have been investigating high-speed computational procedures to allow shield analysis from the preliminary design concepts to the final design. In particular, we will discuss the progress towards a full three-dimensional and computationally efficient deterministic code for which the current HZETRN evaluates the lowest order asymptotic term. HZETRN is the first deterministic solution to the Boltzmann equation allowing field mapping within the International Space Station (ISS) in tens of minutes using standard Finite Element Method (FEM) geometry common to engineering design practice enabling development of integrated multidisciplinary design optimization methods. A single ray trace in ISS FEM geometry requires 14 milliseconds and severely limits application of Monte Carlo methods to such engineering models. A potential means of improving the Monte Carlo efficiency in coupling to spacecraft geometry is given in terms of reconfigurable computing and could be utilized in the final design as verification of the deterministic method optimized design.

Wilson, J. W.↗
NASA NTRS2004

A space radiation transport method development

Improved spacecraft shield design requires early entry of radiation constraints into the design process to maximize performance and minimize costs. As a result, we have been investigating high-speed computational procedures to allow shield analysis from the preliminary design concepts to the final design. In particular, we will discuss the progress towards a full three-dimensional and computationally efficient deterministic code for which the current HZETRN evaluates the lowest-order asymptotic term. HZETRN is the first deterministic solution to the Boltzmann equation allowing field mapping within the International Space Station (ISS) in tens of minutes using standard finite element method (FEM) geometry common to engineering design practice enabling development of integrated multidisciplinary design optimization methods. A single ray trace in ISS FEM geometry requires 14 ms and severely limits application of Monte Carlo methods to such engineering models. A potential means of improving the Monte Carlo efficiency in coupling to spacecraft geometry is given in terms of re-configurable computing and could be utilized in the final design as verification of the deterministic method optimized design. Published by Elsevier Ltd on behalf of COSPAR.

NASA Center JSC↗
DOE OSTI2024

DG-IMEX method for a two-moment model for radiation transport in the $\mathscr{O}$($v$/$c$) limit

Here, we consider neutral particle systems described by moments of a phase-space density and propose a realizability-preserving numerical method to evolve a spectral two-moment model for particles interacting with a background fluid moving with nonrelativistic velocities. The system of nonlinear moment equations, with special relativistic corrections to $\mathscr{O}$($v$/$c$), expresses a balance between phase-space advection and collisions and includes velocity-dependent terms that account for spatial advection, Doppler shift, and angular aberration. The model is conservative for the correct $\mathscr{O}$($v$/$c$) Eulerian-frame number density and is consistent, to $\mathscr{O}$($v$/$c$), with Eulerian-frame energy and momentum conservation. This model is closely related to the one promoted by Lowrie et al. and similar to models currently used to study transport phenomena in large-scale simulations of astrophysical environments. The proposed numerical method is designed to preserve moment realizability, which guarantees that the moments correspond to a nonnegative phase-space density. The realizability-preserving scheme consists of the following key components: (i) a strong stability-preserving implicit-explicit (IMEX) time-integration method; (ii) a discontinuous Galerkin (DG) phase-space discretization with carefully constructed numerical uxes; (iii) a realizability-preserving implicit collision update; and(iv) a realizability-enforcing limiter. In time integration, nonlinearity of the moment model necessitates solution of nonlinear equations, which we formulate as fixed-point problems and solve with tailored iterative solvers that preserve moment realizability with guaranteed global convergence. We also analyze the simultaneous Eulerian-frame number and energy conservation properties of the semi-discrete DG scheme and propose a "spectral redistribution" scheme that promotes Eulerian-frame energy conservation. Through numerical experiments, we demonstrate the accuracy and robustness of this DG-IMEX method and investigate its Eulerian-frame energy conservation properties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1972

Evaluation of one-dimensional approximations for radiative transport in blunt body shock layers

One-dimensional approximations for the nongray radiative flux and flux divergence in radiating shock layers about a blunt entry body are compared with an exact three-dimensional treatment. A coupled radiative-gasdynamic analysis of the shock layer flow about the entire body provided the thermodynamic field used in these comparisons. In terms of calculating the total energy lost by the shock layer, the one-dimensional approximations to the radiative flux divergence are accurate to within a few percent. In terms of calculating the surface flux, the one-dimensional approximations introduce the largest errors of approximately 15 percent near the stagnation point. The source of these errors is the slab-like geometric representation of the shock layer inherent in all one-dimensional models. Finally, for both the radiative flux and its divergence, the tangent slab approximation provides more accurate results than differential approximation methods.

Wilson, K. H.↗
NASA NTRS2007

Comparison of Integrated Radiation Transport Models with TEPC Measurements for the Average Quality Factors in Spaceflights

The purpose of this work is to test our theoretical model for the interpretation of radiation data measured in space. During the space missions astronauts are exposed to the complex field of radiation type and kinetic energies from galactic cosmic rays (GCR), trapped protons, and sometimes solar particle events (SPEs). The tissue equivalent proportional counter (TEPC) is a simple time-dependent approach for radiation monitoring for astronauts on board the International Space Station. Another and a newer approach to Microdosimetry is the use of silicon-on-insulator (SOI) technology launched on the MidSTAR-1 mission in low Earth orbit (LEO). In the radiation protection practice, the average quality factor of a radiation field is defined as a function of linear energy transfer (LET), Q(sub ave)(LET). However, TEPC measures the average quality factor as a function of the lineal energy y, Q(sub ave)(y), defined as the average energy deposition in a volume divided by the average chord length of the volume. Lineal energy, y, deviates from LET due to energy straggling, delta-ray escape or entry, and nuclear fragments produced in the detector volume. Monte Carlo track structure simulation was employed to obtain the response of a TEPC irradiated with charged particle for an equivalent site diameter of 1 micron of wall-less counter. The calculated data of the energy absorption in the wall-less counter were compiled for various y values for several ion types at various discrete projectile energy levels. For the simulation of TEPC response from the mixed radiation environments inside a spacecraft, such as, Space Shuttle and International Space Station, the complete microdosimetric TEPC response, f( y, E, Z), were calculated with the Monte Carlo theoretical results by using the first order Lagrangian interpolation for a monovariate function at a given y value (y = 0.1 keV/micron 5000 keV/micron) at any projectile energy level (E = 0.01 MeV/u to 50,000 MeV/u) of each specific radiation type (Z = 1 to 28). Because the anomalous response has been observed at large event sizes in the experiment due to the escape of energy out of sensitive volume by delta-rays and the entry of delta-rays from the high-density wall into the low-density gas-volume cavity, Monte Carlo simulation was also made for the response of a walled-TEPC with wall thickness 2 mm and density 1 g/cm(exp 3). The radius of cavity was set to 6.35 mm and a gas density 7.874 x 10(exp -5) g/cm(exp 3). The response of the walled- and the wall-less counters were compared. The average quality factor Q(sub ave)(y) for trapped protons on STS-89 demonstrated the good agreement between the model calculations and flight TEPC data as shown. Using an integrated space radiation model (this includes the transport codes HZETRN and BRYNTRN, the quantum nuclear interaction model QMSFRG) and the resultant response distribution functions of walled-TEPC from Monte-Carlo track simulations, we compared model calculations with walled-TEPC measurements from NASA missions in LEO and made predictions for the lunar and the Mars missions. The Q(sub ave)(y) values for the trapped or the solar protons ranged from 1.9-2.5. This over-estimates the Qave(LET) values which ranged from 1.4-1.6. Both quantities increase with shield thickness due to nuclear fragmentation. The Q(sub ave)(LET) for the complete GCR spectra was found to be 3.5-4.5, while flight TEPCs measured 2.9-3.4 for Q(sub ave)(y). The GCR values are decreasing with the shield thickness. Our analysis for a proper interpretation of data supports the use of TEPCs for monitoring space radiation environment.

Kim, Myung-Hee Y.↗
NASA NTRS2006

Comparisons of Integrated Radiation Transport Models with Microdosimetry Data in Spaceflight

Astronauts are exposed to galactic cosmic rays (GCR), trapped protons, and possible solar particle events (SPE) during spaceflight. For such complicated mixtures of radiation types and kinetic energies, tissue equivalent proportional counters (TEPC's) represent a simple time-dependent approach for radiation monitoring. Of interest in radiation protection is the average quality factor of a radiation field defined as a function of linear energy transfer, LET, Q(sub ave)(LET). However TEPC's measure the average quality factors as a function of lineal energy (y), Q(sub ave)(y) defined as the average energy deposition in a volume divided by the average chord length of the volume. Lineal energy, y deviates from LET due to energy straggling, delta-ray escape or entry, and nuclear fragments produced in the detector. Using integrated space radiation models that includes the transport code HZETRN/BRYNTRN, the quantum nuclear interaction model, QMSFRG, and results from Monte-Carlo track simulations of TEPC's response to ions, we consider comparisons of model calculations to TEPC results from NASA missions in low Earth orbit and make predictions for lunar and Mars missions. Good agreement between the model and measured spectra from past NASA missions is found. A finding of this work is that TEPC's values for trapped or solar protons of Q(sub ave)(y) range from 1.9-2.5, overestimating Q(sub ave)(LET), which ranges from 1.4-1.6 with both quantities increasing with shielding depth due to nuclear secondaries Comparisons for the complete GCR spectra show that Q(sub ave)(LET) for GCR is approximately 3.5-4.5, while TEPC's measure 2.9-3.4 for Q(sub ave)(y) with the GCR values decreasing with depth as heavy ions are absorbed in shielding material. Our results support the use of TEPC's for space radiation environmental monitoring when computational analysis is used for proper data interpretation.

Cucinotta, Francis A.↗
DOE OSTI2025

Global Simulations of Gravitational Instability in Protostellar Disks with Full Radiation Transport. II. Locality of Gravitoturbulence, Clumpy Spirals, and Implications for Observable Substructure

Spiral perturbations in a gravitationally unstable accretion disk regulate disk evolution through angular momentum transport and heating and provide an observational signature of gravitational instability (GI). We use global 3D simulations to systematically characterize and understand these spiral perturbations. The spiral perturbations and the resulting transport are overall insensitive to the cooling type, with the exception that radiative cooling, especially in the optically thick regime, reduces the amplitude of temperature perturbations. Spiral perturbations are localized around corotation, allowing transport to be approximated by a local α viscosity to zeroth order in aspect ratio (H/R), but only after averaging over multiple orbits in time and/or multiple scale heights in space. Meanwhile, large-amplitude perturbations from strong gravitoturbulence can cause $\mathcal{O}$(α 1/2 ) deviation in the cooling rate of the disk. We develop empirical prescriptions for the angular momentum transport, heating, and cooling in a gravitoturbulent disk that capture the deviation from a viscous, unperturbed disk to first order in H/R and α 1/2 . The spiral perturbations in saturated gravitoturbulence are clumpy, with dense clumps forming through the nonlinear coupling between multiple modes at different m. Observationally, the clumpy gravitoturbulence produced by saturated GI can be mistaken with observational noise or embedded companions, especially under finite resolution. Meanwhile, grand-design spirals with m-fold symmetry may be uncommon among disks in saturated gravitoturbulence, and we speculate that they may instead be a signature of recently triggered or decaying GI.

Gravitational instability↗
NASA NTRS1999

Diffusive and Radiative Transport in Fires (DARTFire): Opposed-Flow Flame Spread in Low-Velocity Flows

For flames spreading into a low-velocity flow that can only be obtained in microgravity, we have observed behavior that is different from that which is obtained at higher velocities where radiative effects are unimportant and species transport is relatively fast. Unfortunately, lack of a large body of low-gravity flame spread experimental data inhibits progress in developing an understanding of the physics of low-velocity, opposed-flow flame spread phenomena. Recent DARTFire sounding rocket experimental studies though, coupled with developing theory and modelling, have allowed some strides in understanding to be made, on which we report here. Four launches to date have resulted in six experiments for opposed-flow flame spread over a thick PMMA sample. During the 6 min microgravity period, the PMMA samples were ignited, and steady flame spread was studied under varied flow velocity, oxidizer atmospheric conditions, and, because radiative heat transfer is so important in these slowly spreading flames, external radiant flux. These were the first attempts at such experimental control and measurement in microgravity. A recent reflight of the Solid Surface Combustion Experiment (SSCE) has demonstrated, as modelling predicts, that for the thick fuel of the DARTFire experiment, flame spread in a quiescent environment is a transient process evolving from ignition to extinction on the order of 600 s (Altenkirch et al., 1999). Further study then of the effects of radiation in the very low-velocity opposing flows is of interest in understanding the transition from steady, sustained spread to the unsteady evolution to extinction as the opposing flow is reduced further and eventually removed.

Altenkirch, R. A.↗
NASA NTRS2008

Radiation Transport of Heliospheric Lyman-alpha from Combined Cassini and Voyager Data Sets

Heliospheric neutral hydrogen scatters solar Lyman-alpha radiation from the Sun with '27-day' intensity modulations observed near Earth due to the Sun's rotation combined with Earth's orbital motion. These modulations are increasingly damped in amplitude at larger distances from the Sun due to multiple scattering in the heliosphere, providing a diagnostic of the interplanetary neutral hydrogen density independent of instrument calibration. This paper presents Cassini data from 2003-2004 obtained downwind near Saturn at approximately 10 AU that at times show undamped '27-day' waves in good agreement with the single-scattering models of Pryor et al., 1992. Simultaneous Voyager 1 data from 2003- 2004 obtained upwind at a distance of 88.8-92.6 AU from the Sun show waves damped by a factor of -0.21. The observed degree of damping is interpreted in terms of Monte Carlo multiple-scattering calculations (e.g., Keller et al., 1981) applied to two heliospheric hydrogen two-shock density distributions (discussed in Gangopadhyay et al., 2006) calculated in the frame of the Baranov-Malama model of the solar wind interaction with the two-component (neutral hydrogen and plasma) interstellar wind (Baranov and Malama 1993, Izmodenov et al., 2001, Baranov and Izmodenov, 2006). We conclude that multiple scattering is definitely occurring in the outer heliosphere. Both models compare favorably to the data, using heliospheric neutral H densities at the termination shock of 0.085 cm(exp -3) and 0.095 cm(exp -3). This work generally agrees with earlier discussions of Voyager data in Quemerais et al., 1996 showing the importance of multiple scattering but is based on Voyager data obtained at larger distances from the Sun (with larger damping) simultaneously with Cassini data obtained closer to the Sun.

Pryor, W.↗
DOE OSTI2023

Variance Reduction within Implicit Monte Carlo Thermal Radiation Transport using the Local Importance Function Transform [Slides]

System of equations derived by Fleck and Cummings to model time-dependent thermal radiative transfer (TRT) problems: Solution of radiation specific intensity (I), non-linearly coupled to material internal energy (U); Absorption and re-emission approximated by effective scattering events. Useful for high-energy density physics simulations: i.e., astrophysics, inertial confinement fusion (ICF). Whereas deterministic methods are fully discretized and free of statistical noise, Monte Carlo methods allow for dynamic sampling of the phase space at the price of statistical noise; Monte Carlo methods are also characterized by a slow $\frac{1}{√Ν}$ convergence rate.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗