Search NASA⌕ Search

SEARCH · Search NASA

Results for “Radiative transfer equation”

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 379 records · Page 21

Radiative transfer with partial frequency redistribution in inhomogeneous atmospheres - Application to the Jovian aurora

A direct finite difference numerical solution for the equation of radiative transfer by the Feautrier method is developed for use in planetary atmospheres. The procedure described here uses a plane-parallel atmosphere, and can treat partial frequency redistribution, inhomogeneity, external or internal sources, and various boundary conditions. Isotropic scattering is assumed, but in the case of no frequency redistribution, Rayleigh scattering can also be handled. A program utilizing this method is tested in a variety of situations against more powerful and elaborate methods. The case of the Lyman alpha aurora on Jupiter is then considered, where the effects of partial frequency redistribution are shown to be of great importance. New results for the detailed line profiles for Lyman alpha in the Jovian aurora are presented. The method is quite versatile, and should be especially useful in studying a wide range of problems related to auroral or dayglow emissions in planetary atmospheres.

Gladstone, G. R.↗

A Model of Radiative and Conductive Energy Transfer in Planetary Regoliths

The thermal regime in planetary regoliths involves three processes: propagation of visible radiation, propagation of thermal radiation, and thermal conduction. The equations of radiative transfer and heat conduction are formulated for particulate media composed of anisotropically scattering particles. Although the equations are time dependent, only steady state problems are considered in this paper. Using the two-stream approximation, solutions are obtained for two cases: a layer of powder heated from below and an infinitely thick regolith illuminated by visible radiation. Radiative conductivity, subsurface temperature gradients, and the solid state greenhouse effect all appear intrinsically in the solutions without ad hoc additions. Although the equations are nonlinear, approximate analytic solutions that are accurate to a few percent are obtained. Analytic expressions are given for the temperature distribution, the optical and thermal radiance distributions, the hemispherical albedo, the hemispherical emissivity, and the directional emissivity. Additional applications of the new model to three problems of interest in planetary regoliths are presented by Hapke.

Hapke, Bruce↗

Comparison of Polarized Radiative Transfer Codes Used by the EHT Collaboration

Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT Collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I , Q , U , and V , respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

Analytic approach to photoelectron transport.

The equation governing the transport of photoelectrons in the ionosphere is shown to be equivalent to the equation of radiative transfer. In the single-energy approximation this equation is solved in closed form by the method of discrete ordinates for isotropic scattering and for a single-constituent atmosphere. The results include prediction of the angular distribution of photoelectrons at all altitudes and, in particular, the angular distribution of the escape flux. The implications of these solutions in real atmosphere calculations are discussed.

Stolarski, R. S.↗

Influence of aerosol scattering on atmospheric blurring of surface features

The influence of the shape of atmospheric scattering phase functions on blurring of surface detail in images acquired from space is investigated. The effects are characterized by computing atmospheric transfer functions and by solving the multidimensional equation of radiative transfer using a Fourier transform method. It is predicted that increased forward scattering enhances the detectability of atmospheric blurring near reflectance boundaries. Results for off-nadir viewing are also presented.

Diner, D. J.↗

Radiative transfer in dusty nebulae

The effects of dust scattering on observable optical and infrared parameters, and the accuracy of approximate solutions were examined. The equation of radiative transfer in a static and homogeneous, but not necessarily uniform, distribution gas and dust around a central empty core with a point source of energy at its center was solved. The dust properties were characterized by a phenomenological extinction cross section, albedo and parameters describing the anisotropy of dust scattering. For ultraviolet photons, ionization equilibrium equations for the gas were solved, and for infrared photons a self-consistent dust temperature was calculated. Ray tracing was used to solve for the angular dependence of the intensity.

Dana, R. A.↗

A note on the variational method of Stokes and DeMarcus for radiative transfer in planetary atmospheres

A generalization of the variational principle of Stokes and DeMarcus (1971) by the construction of two bifunctionals is proposed. The equation of radiative transfer is obtained by imposing the stationary condition on such functionals, and the reflectivity and transmissivity at arbitrary angle are given by their extremum values. Further simplification is possible when trail functions are used that are linear in the variational parameters, in which case the source function derived from the method of Stokes and DeMarcus can give the reflectivity at arbitrary emergent angle to second order. The usefulness and limitations of these variational priinciples are explored with extensive numerical results.

Chow, K. S. K.↗

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↗

Microwave remote sensing of a two-layer random medium

A two-variable expansion technique is used to solve for the mean Green's. function from the Dyson equation under the nonlinear approximation. The Bethe-Salpeter equation then gives rise to a set of modified radiative transfer (MRT) equations which accommodate coherent effects essential to bounded media. It is found that the nonlinear approximation, instead of the more popular bilocal approximation, should be used for the case of bounded media. The two approximations yield identical results for unbounded media. The MRT equations are then solved for a two-layer random medium. The MRT equations give rise to simple and useful solutions which are applicable to both active and passive microwave remote sensing.

Tsang, L.↗

A Lattice Boltzmann Scheme for Polarized Radiative Transfer in Planetary Atmospheres

Quantitative calculation of reflected polarized light from planets is of great significance to interpret the polarimetric observations of planetary atmospheres. In this work, we propose a lattice Boltzmann (LB) scheme for polarized radiative transfer in planetary atmospheres. The LB scheme obtains the Stokes vector through the vector LB equation performing simple collision and streaming processes. The vector radiative transfer equation is rigorously derived from the vector LB equation via the Maxwell iteration technique. Polarized radiative transfer of a single point on the planetary surface with Rayleigh and Mie scattering atmospheres are first accurately solved by our LB scheme. Afterward, we systematically investigate the disk-integrated polarization of finite and semi-infinite, conservative and nonconservative Rayleigh scattering planetary atmospheres. The disk-integrated phase curves of the Stokes parameters and degree of polarization and spherical and geometric albedos are produced, which are in good agreement with the benchmark results for different cases. The numerical results indicate that our LB scheme is efficient and accurate for polarized radiative transfer in planetary atmospheres. Our LB scheme is expected to provide a competitive numerical tool to interpret the polarimetric observations of planetary atmospheres.

47 OTHER INSTRUMENTATION↗

Unsteady, combined radiation and conduction in an absorbing, scattering, and emitting medium.

The transient cooldown of a gray, absorbing, isotropic scattering, emitting, and conducting medium bounded by gray, diffusely emitting and reflecting parallel plates is considered. Numerical solutions are obtained for the initial boundary-value problem with a discontinuous decrease in temperature at one boundary. The quasi-steady equation of radiative transfer is solved using Gaussian quadrature and a matrix eigenvector technique together with explicit numerical solution of the unsteady energy equation. Temperature and energy flux distributions are presented for variations of optical thickness, boundary emissivity, albedo, and conduction-radiation parameter.

Weston, K. C.↗

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↗

Elastic wave propagation in a highly scattering medium - A diffusion approach

The propagation of elastic waves in the moon, where the first seismograms were characterized by the presence of a long coda attributed to strongly scattered waves, is modeled with the aid of the time-dependent equation of radiative transfer. The average energy density as a function of time and space is described by the diffusion equation with linear dissipation on the assumption that all the energy present has been scattered many times and the time and distance scales of the problem are long compared to the scales of the scattering process. Ultrasonic experiments in the laboratory confirm the applicability of the formalism.

Dainty, A. M.↗

On time-dependent radiative transfer

An integral equation is developed for application to time-dependent laboratory experiments in which partial redistribution effects are important. The equation of transport with the Heasly-Kneer (1976) emission coefficient and the equation of statistical equilibrium lead to a time-dependent redistribution function containing an absorption - reemission term which decays exponentially in time and a scattering term which is instantaneous. This integral equation does not agree with an equation written by Payne et al. (1974) that has been used to compare theory with experiments. The difference between the Payne equation and the equation developed here needs to be examined in detail, since it might under some circumstances be on the same order as the difference between partial and complete redistribution.

Streater, A.↗

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↗

Solution of an integral equation encountered in studies on radiative transfer in completely absorbing leaf canopies

The leaf normal distribution function is related to the extinction coefficient of direct solar radiation by a Fredholm integral equation of the first kind. The integral equation is solved by the method of constrained least squares. The solutions are stable and correct only in extreme cases and, in general, it is not possible to recover an arbitrary leaf normal distribution function.

Myneni, Ranga B.↗

A Unified Interpretation of Variability in Precipitation Isotope Ratios

Abstract Several mechanisms have been proposed to explain why the isotope ratios of precipitation vary in space and time and why they correlate with other climate variables like temperature and precipitation. Here, we argue that this behavior is best understood through the lens of radiative transfer, which treats the depletion of atmospheric vapor transport by precipitation as analogous to the attenuation of light by absorption or scattering. Building on earlier work by Siler et al., we introduce a simple model that uses the equations of radiative transfer to approximate the two-dimensional pattern of the oxygen isotope composition of precipitation ( δ p ) from monthly mean hydrologic variables. The model accurately simulates the spatial and seasonal variability in δ p within a state-of-the-art climate model and permits a simple decomposition of δ p variability into contributions from gradients in evaporation and the length scale of vapor transport. Outside the tropics, δ p is mostly controlled by gradients in evaporation, whose dependence on temperature explains the positive correlation between δ p and temperature (i.e., the temperature effect). At low latitudes, δ p is mostly controlled by gradients in the transport length scale, whose inverse relationship with precipitation explains the negative correlation between δ p and precipitation (i.e., the amount effect). This suggests that the temperature and amount effects are both mostly explained by the variability in upstream rainout, but they reflect distinct mechanisms governing rainout at different latitudes. Significance Statement The isotopic composition of precipitation has long been used to make inferences about past climates based on its observed relationship with precipitation in the tropics and with temperature at higher latitudes. These relationships—known as the “amount effect” and “temperature effect,” respectively—have been attributed to many different mechanisms, most of which are thought to operate at either high or low latitudes but not both. Here, we present a unified framework for interpreting the isotope variability that can explain the latitude dependence of the temperature and amount effects despite making no distinction between high and low latitudes. Although our results are generally consistent with certain interpretations of the amount effect, they suggest that the temperature effect is widely misunderstood.

54 ENVIRONMENTAL SCIENCES↗