Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermal Radiative Transfer”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A Family of Multi-Dimensional Thermal Radiative Transfer Test Problems

Many thermal radiative transfer (TRT) test problems have been introduced over the years. Here, we will combine the good features of several tests into one new problem. Our goals are to have a multi-dimensional, multi-group TRT problem that stresses the physics and numerics of codes in a realistic way, but the geometry and materials are idealized just enough to be simple to set up in any code. The proposed problem has no connection to any physical experiment; it is simply meant to have a combination of properties that stress the solvers with both the optically thick and thin limits. To achieve this, we are willing to give up any form of analytic solution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A Monte Carlo Thermal Radiative Transfer Solver with Nonlinear Elimination

Here in this paper, we present a new Monte Carlo method for solving the thermal radiative transfer (TRT) equations via the method of nonlinear elimination (NLEM). This method is inspired by the previous application of NLEM to thermal radiation diffusion. Our approach, called diffusion accelerated Implicit Monte Carlo (DAIMC), is a hybrid technique which combines a Monte Carlo method for solving a purely-absorbing transport equation and a diffusion solution that accounts for effective scattering, or absorption–reemission. The method aims to improve the implicitness of the traditional implicit Monte Carlo (IMC) method. We derive DAIMC generally for 3D Cartesian geometries, but in this paper, we present results and analysis in 1D slab geometry. These preliminary results indicate that DAIMC implementations may provide more accurate and robust TRT solutions than IMC in certain test problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Machine Learning for Memory Reduction in the Implicit Monte Carlo Simulations of Thermal Radiative Transfer [Slides]

Project Goal: Use parametric Machine Learning methods in order to reduce memory requirements at checkpointing & restarting in the IMC simulations of Thermal Radiative Transfer using: Expectation Maximization and Weighted Gaussian Mixture Model-based approach for `particle-data compression', introduced in Plasma Physics to model Maxwellian particle distributions by Luis Chacon and Guangye Chen; Expectation Maximization with Weighted Hyper-Erlang Model in order to compress isotropic IMC particle data in the frequency domain; and Expectation Maximization and von Mises-Fisher Mixture Model for compression of anisotropic IMC particle data in the angular domain (work-in-progress).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

A Variable Eddington Factor Model for Thermal Radiative Transfer with Closure Based on Data-Driven Shape Function

Here, a new variable Eddington factor (VEF) model is presented for nonlinear problems of thermal radiative transfer (TRT). The VEF model is data-driven and acts on known (a-priori) radiation-diffusion solutions for material temperatures in the TRT problem. A linear auxiliary problem is constructed for the radiative transfer equation (RTE) whose emission source and opacities are evaluated at these known material temperatures. The solution to this RTE approximates the specific intensity distribution in phase-space and time. It is applied as a shape function to define the Eddington tensor for the presented VEF model. The shape function computed via the auxiliary RTE problem will capture some degree of transport effects within the TRT problem. The VEF moment equations closed with this approximate Eddington tensor will thus carry with them these captured transport effects. In this study, the temperature data comes from multigroup P 1 , P 1/3 , and flux-limited diffusion radiative transfer models. The proposed VEF model can be interpreted as a transport-corrected diffusion reduced-order model. Numerical results are presented on the Fleck-Cummings test problem which models a supersonic wavefront of radiation. The VEF model is shown to improve accuracy by 1–2 orders of magnitude compared to the considered radiation-diffusion model solutions to the TRT problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Reduced order models for thermal radiative transfer problems based on moment equations and data-driven approximations of the Eddington tensor

Here a new group of structure and asymptotic preserving reduced-order models (ROMs) for multidimensional nonlinear thermal radiative transfer (TRT) problems is presented. They are formulated by means of the nonlinear projective approach and data compression techniques. The nonlinear projection is applied to the Boltzmann transport equation (BTE) to derive a hierarchy of low-order moment equations. Approximation of the Eddington tensor that provides exact closure for the system of moment equations is found with projection-based data-driven methodologies. These include the (i) proper orthogonal decomposition (POD), (ii) dynamic mode decomposition (DMD) and (iii) a variant of the DMD. A parameterization is derived for this ROM for the temperature of radiation incoming to the problem domain (the radiation drive temperature). This parameterization is informed from results of a dimensionless study of the TRT problem. Analysis of the ROMs is performed on the classical Fleck-Cummings TRT multigroup test problem in 2D geometry with a radiation-driven Marshak wave. Numerical results are presented to demonstrate the performance of these ROMs for the simulation of evolving radiation and heat waves. Results show these models to be sufficiently accurate for practical computations with rather low-rank representations of the Eddington tensor. As the rank of the approximation is increased, the errors of solutions generated by the ROMs gradually decreases.

42 ENGINEERING↗
DOE OSTI2024

A Survey of Multimaterial Treatments for Thermal Radiative Transfer

Arbitrary Lagrangian-Eulerian methods are a popular choice for hydrodynamic modeling in radiation (rad-hydro) simulations. Because these methods involve a relaxation step that moves the mesh relative to material boundaries, multimaterial spatial zones are generally present. Accurate treatments of these zones are needed to resolve various physical phenomena of interest for inertial confinement fusion applications. However, these codes are often paired with single-material, deterministic thermal radiative transfer (TRT) codes that are oblivious to the material compositions of each zone. These single-material TRT codes can only accept homogenized material properties (opacities, specific heats, etc.) from the hydrodynamic code and output homogenized solutions. After each TRT time step, the multimaterial hydrodynamic code must dehomogenize the quantities computed by the TRT package in order to update subzonal material temperatures. The process by which hydrodynamic codes perform this dehomogenization has not been well documented in previous literature, and the methods can vary significantly from code to code. The purpose of this paper is to document, study, and compare existing techniques used for rad-hydro simulations as well as present a new method with potentially promising results. We summarize several methods and give comparisons on infinite-medium problems as well a finite-medium problem for two of the methods.

42 ENGINEERING↗
NASA NTRS2019

Cloud Scattering Impact on Thermal Radiative Transfer and Global Longwave Radiation

The potential importance of longwave (LW) cloud scattering has been recognized but the actual estimate of this effect on thermal radiation varies greatly among different studies. General circulation models (GCMs) generally neglect or simplify the multiple scattering in the LW. In this study, we use a rigorous radiative transfer algorithm to explicitly consider LW multiple-scattering and apply the GCM to quantify the impact of cloud LW scattering on thermal radiation fluxes. Our study shows that the cloud scattering effect on downward thermal radiation at the surface is concentrated in the infrared atmospheric window spectrum (800–1250 cm9exp −1)). The scattering effect on the outgoing longwave radiation (OLR) is also present in the window region over low clouds but it is mainly in the far-infrared spectrum (300–600 cm(exp −1)) over high clouds. For clouds with small to moderate optical depth (τ < 10), the scattering effect on thermal fluxes shows large variation with the cloud τ and has a maximum at an optical depth of ∼3. For opaque clouds, the scattering effect approaches an asymptote and is smaller and less important. The 2-stream radiative transfer scheme could have an error over 10% with an RMS error around 3.5%–4.0% in the calculated LW flux. This algorithm error of the 2-stream approximation could readily exceed the no-scattering error in the LW, and thus it is worthless to include the time-consuming computation of multiple scattering in a 2-stream radiative transfer scheme. However, the calculation error rapidly decreases as stream number increases and the RMS error in LW flux using the 4-stream scheme is under 0.3%, an accuracy sufficient for most climate studies. We implement the 4-stream discrete-ordinate algorithm in the GISS GCM and run the GCM for 20 years with and without the LW scattering effect, respectively. When cloud LW scattering is included, we find that the global annual mean OLR is reduced by 2.7 W/m(exp 2), and the downward surface flux and the net atmospheric absorption are increased by 1.6 W/m2 and 1.8 W/m(exp 2), respectively. Using one year of ISCCP clouds and running the standalone radiative transfer offline, the global annual mean non-scattering errors in OLR, surface LW downward flux and net atmospheric absorption are 3.6 W/m(exp 2), −1.1 W/m(exp 2), and −2.5 W/m(exp 2), respectively. The global scattering impact of 2.7 W/m(exp 2) on the OLR is small when compared to the typical global OLR value of 240 W/m2, but it is significant when compared to cloud LW radiative forcing (30 W/m2) and net cloud forcing (−14 W/m(exp 2)). Overall, the effect of neglecting scattering on the thermal fluxes is comparable to the reported clear sky radiative effect of doubling CO2.

longwave cloud scattering↗
NASA NTRS2017

A Fast Hybrid (3-D/1-D) Model for Thermal Radiative Transfer in Cirrus via Successive Orders of Scattering

We investigate the impact of cirrus cloud heterogeneity on the direct emission by cloud or surface and on the scattering by ice particles in the thermal infrared (TIR). Realistic 3-D cirri are modeled with the 3DCLOUD code, and top-of-atmosphere radiances are simulated by the 3-D Monte Carlo radiative transfer (RT) algorithm 3DMCPOL for two (8.65 micrometers and 12.05 micrometers) channels of the Imaging Infrared Radiometer on CALIPSO. At nadir, comparisons of 1-D and 3-D RT show that 3-D radiances are larger than their 1-D counterparts for direct emission but smaller for scattered radiation. For our cirrus cases, 99% of the 3-D total radiance is computed by the third scattering order, which corresponds to 90% of the total computational effort, but larger optical thicknesses need more scattering orders. To radically accelerate the 3-D RT computations (using only few percent of 3-D RT time with a Monte Carlo code), even in the presence of large optical depths, we develop a hybrid model based on exact 3-D direct emission, the first scattering order from 1-D in each homogenized column, and an empirical adjustment linearly dependent on the optical thickness to account for higher scattering orders. Good agreement is found between the hybrid model and the exact 3-D radiances for two very different cirrus models without changing the empirical parameters. We anticipate that a future deterministic implementation of the hybrid model will be fast enough to process multiangle thermal imagery in a practical tomographic reconstruction of 3-D cirrus fields.

radiative transfer↗
DOE OSTI2020

Efficient smoothed particle radiation hydrodynamics I: Thermal radiative transfer

This work presents efficient solution techniques for radiative transfer in the smoothed particle hydrodynamics discretization. Two choices that impact efficiency are how the material and radiation energy are coupled, which determines the number of iterations needed to converge the emission source, and how the radiation diffusion equation is solved, which must be done in each iteration. The coupled material and radiation energy equations are solved using an inexact Newton iteration scheme based on nonlinear elimination, which reduces the number of Newton iterations needed to converge within each time step. During each Newton iteration, the radiation diffusion equation is solved using Krylov iterative methods with a multigrid preconditioner, which abstracts and optimizes much of the communication when running in parallel. The code is verified for an infinite medium problem, a one-dimensional Marshak wave, and a two and three-dimensional manufactured problem, and exhibits first-order convergence in time and second-order convergence in space. For these problems, the number of iterations needed to converge the inexact Newton scheme and the diffusion equation is independent of the number of spatial points and the number of processors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Data-driven acceleration of thermal radiation transfer calculations with the dynamic mode decomposition and a sequential singular value decomposition

In this work, we present a method for accelerating discrete ordinates radiative transfer calculations for radiative transfer. Our method works with nonlinear positivity fixes, in contrast to most acceleration schemes. The method is based on the dynamic mode decomposition (DMD) and using a sequence of rank-one updates to compute the singular value decomposition needed for DMD. Using a sequential method allows us to automatically determine the number of solution vectors to include in the DMD acceleration. We present results for slab geometry discrete ordinates calculations with the standard temperature linearization. Compared with positive source iteration, our results demonstrate that our acceleration method reduces the number of transport sweeps required to solve the problem by a factor of about 3 on a standard diffusive Marshak wave problem, a factor of several thousand on a cooling problem where the effective scattering ratio approaches unity, and a factor of 20 improvement in a realistic, multimaterial radiating shock problem.

97 MATHEMATICS AND COMPUTING↗