Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Radiation Scattering Equations”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Computational attributes of the integral form of the equation of transfer

Difficulties can arise in radiative and neutron transport calculations when a highly anisotropic scattering phase function is present. In the presence of anisotropy, currently used numerical solutions are based on the integro-differential form of the linearized Boltzmann transport equation. This paper, departs from classical thought and presents an alternative numerical approach based on application of the integral form of the transport equation. Use of the integral formalism facilitates the following steps: a reduction in dimensionality of the system prior to discretization, the use of symbolic manipulation to augment the computational procedure, and the direct determination of key physical quantities which are derivable through the various Legendre moments of the intensity. The approach is developed in the context of radiative heat transfer in a plane-parallel geometry, and results are presented and compared with existing benchmark solutions. Encouraging results are presented to illustrate the potential of the integral formalism for computation. The integral formalism appears to possess several computational attributes which are well-suited to radiative and neutron transport calculations.

Frankel, J. I.↗
NASA NTRS2016

Expansion of Tabulated Scattering Matrices in Generalized Spherical Functions

An efficient way to solve the vector radiative transfer equation for plane-parallel turbid media is to Fourier-decompose it in azimuth. This methodology is typically based on the analytical computation of the Fourier components of the phase matrix and is predicated on the knowledge of the coefficients appearing in the expansion of the normalized scattering matrix in generalized spherical functions. Quite often the expansion coefficients have to be determined from tabulated values of the scattering matrix obtained from measurements or calculated by solving the Maxwell equations. In such cases one needs an efficient and accurate computer procedure converting a tabulated scattering matrix into the corresponding set of expansion coefficients. This short communication summarizes the theoretical basis of this procedure and serves as the user guide to a simple public-domain FORTRAN program.

scattering matrix↗
NASA NTRS1983

A mathematical characterization of vegetation effect on microwave remote sensing from the Earth

In passive microwave remote sensing of the earth, a theoretical model that utilizes the radiative transfer equations was developed to account for the volume scattering effects of the vegetation canopy. Vegetation canopies such as alfalfa, sorghum, and corn are simulated by a layer of ellipsoidal scatterers and cylindrical structures. The ellipsoidal scatterers represent the leaves of vegetation and are randomly positioned and oriented. The orientation of ellipsoids is characterized by a probability density function of Eulerian angles of rotation. The cylindrical structures represent the stalks of vegetation and their radii are assumed to be much smaller than their lengths. The underlying soil is represented by a half-space medium with a homogeneous permittivity and uniform temperature profile. The radiative transfer quations are solved by a numerical method using a Gaussian quadrature formula to compute both the vertical and horizontal polarized brightness temperature as a function of observation angle. The theory was applied to the interpretation of experimental data obtained from sorghum covered fields near College Station, Texas.

Choe, Y.↗
NASA NTRS1995

On the Calculation of Anisotropic Extinction Coefficients for Rigid Fibrous Ceramic Insulations

The specific anisotropic extinction coefficient e* couples the effective Radiative properties of a fibrous insulation into the radiation diffusion equation. This coefficient can be calculated using Mie scattering theory if fiber diameters, refractive indices and fiber orientation distributions are known. In general, fiber orientation distributions are not readily accessible and past calculations have considered fibers as either randomly distributed or normal to the heat flow direction. In certain rigid fibrous ceramic insulations neither of these cases apply well and a simple procedure is described for approximating e* from values calculated for the random and normal orientation cases. The intrinsic error associated with this scaling procedure is investigated. Numerical computations for several test structures and fiber materials show the average error to be less than 5% for net heat flux and radiation conductivity calculations.

Marshall, Jochen↗
NASA NTRS1996

Two-Flux Green's Function Analysis for Transient Spectral Radiation in a Composite

An analysis is developed for obtaining transient temperatures in a two-layer semitransparent composite with spectrally dependent properties. Each external boundary of the composite is subjected to radiation and convection. The two-flux radiative transfer equations are solved by deriving a Green's function. This yields the local radiative heat source needed to numerically solve the transient energy equation. An advantage of the two-flux method is that isotropic scattering is included without added complexity. The layer refractive indices are larger than one. This produces internal reflections at the boundaries and the internal interface; the reflections are assumed diffuse. Spectral results using the Green's function method are verified by comparing with numerical solutions using the exact radiative transfer equations. Transient temperature distributions are given to illustrate the effect of radiative heating on one side of a composite with external convective cooling. The protection of a material from incident radiation is illustrated by adding scattering to the layer adjacent to the radiative source.

Siegel, Robert↗
NASA NTRS1993

Calculation of the angular radiance distribution for a coupled atmosphere and canopy

The radiative transfer equations for a coupled atmosphere and canopy are solved numerically by an improved Gauss-Seidel iteration algorithm. The radiation field is decomposed into three components: unscattered sunlight, single scattering, and multiple scattering radiance for which the corresponding equations and boundary conditions are set up and their analytical or iterational solutions are explicitly derived. The classic Gauss-Seidel algorithm has been widely applied in atmospheric research. This is its first application for calculating the multiple scattering radiance of a coupled atmosphere and canopy. This algorithm enables us to obtain the internal radiation field as well as radiances at boundaries. Any form of bidirectional reflectance distribution function (BRDF) as a boundary condition can be easily incorporated into the iteration procedure. The hotspot effect of the canopy is accommodated by means of the modification of the extinction coefficients of upward single scattering radiation and unscattered sunlight using the formulation of Nilson and Kuusk. To reduce the computation for the case of large optical thickness, an improved iteration formula is derived to speed convergence. The upwelling radiances have been evaluated for different atmospheric conditions, leaf area index (LAI), leaf angle distribution (LAD), leaf size and so on. The formulation presented in this paper is also well suited to analyze the relative magnitude of multiple scattering radiance and single scattering radiance in both the visible and near infrared regions.

Liang, Shunlin↗
NASA NTRS2013

Some New Lidar Equations for Laser Pulses Scattered Back from Optically Thick Media Such as Clouds, Dense Aerosol Plumes, Sea Ice, Snow, and Turbid Coastal Waters

I survey the theoretical foundations of the slowly-but-surely emerging field of multiple scattering lidar, which has already found applications in atmospheric and cryospheric optics that I also discuss. In multiple scattering lidar, returned pulses are stretched far beyond recognition, and there is no longer a one-to-one connection between range and return-trip timing. Moreover, one can exploit the radial profile of the diffuse radiance field excited by the laser source that, by its very nature, is highly concentrated in space and collimated in direction. One needs, however, a new class of lidar equations to explore this new phenomenology. A very useful set is derived from radiative diffusion theory, which is found at the opposite asymptotic limit of radiative transfer theory than the conventional (single-scattering) limit used to derive the standard lidar equation. In particular, one can use it to show that, even if the simple time-of-flight-to-range connection is irretrievably lost, multiply-scattered lidar light can be used to restore a unique profiling capability with coarser resolution but much deeper penetration into a wide variety of optical thick media in nature. Several new applications are proposed, including a laser bathymetry technique that should work for highly turbid coastal waters.

turbid water↗
NASA NTRS1977

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.↗
NASA NTRS1985

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.↗
NASA NTRS1999

Bidirectional Reflectance of a Macroscopically Flat, High-Albedo Particulate Surface: An Efficient Radiative Transfer Solution and Applications to Regoliths

Many remote sensing applications rely on accurate knowledge of the bidirectional reflection function (BRF) of surfaces composed of discrete, randomly positioned scattering particles. Theoretical computations of BRFs for plane-parallel particulate layers are usually reduced to solving the radiative transfer equation (RTE) using one of existing exact or approximate techniques. Since semi-empirical approximate approaches are notorious for their low accuracy, violation of the energy conservation law, and ability to produce unphysical results, the use of numerically exact solutions of RTE has gained justified popularity. For example, the computation of BRFs for macroscopically flat particulate surfaces in many geophysical publications is based on the adding-doubling (AD) and discrete ordinate (DO) methods. A further saving of computer resources can be achieved by using a more efficient technique to solve the plane-parallel RTE than the AD and DO methods. Since many natural particulate surfaces can be well represented by the model of an optically semi-infinite, homogeneous scattering layer, one can find the BRF directly by solving the Ambartsumian's nonlinear integral equation using a simple iterative technique. In this way, the computation of the internal radiation field is avoided and the computer code becomes highly efficient and very accurate and compact. Furthermore, the BRF thus obtained fully obeys the fundamental physical laws of energy conservation and reciprocity. In this paper, we discuss numerical aspects and the computer implementation of this technique, examine the applicability of the Henyey-Greenstein phase function and the sigma-Eddington approximation in BRF and flux calculations, and describe sample applications demonstrating the potential effect of particle shape on the bidirectional reflectance of flat regolith surfaces. Although the effects of packing density and coherent backscattering are currently neglected, they can also be incorporated. The FORTRAN implementation of the technique is available on the World Wide Web, and can be applied to a wide range of remote sensing problems. BRF computations for undulated (macroscopically rough) surfaces are more complicated and often rely on time consuming Monte Carlo procedures. This approach is especially inefficient for optically thick, weakly absorbing media (e.g., snow and desert surfaces at visible wavelengths since a photon may undergo many internal scattering events before it exists the medium or is absorbed. However, undulated surfaces can often be represented as collections of locally flat tilted facets characterized by the BRF found from the traditional plane parallel RTE. In this way the MOnte Carlo procedure could be used only to evaluate the effects of surface shadowing and multiple surface reflections, thereby bypassing the time-consuming ray tracing inside the medium and providing a great savings of CPU time.

Mishchenko, Michael I.↗
NASA NTRS2002

The Space-Time CE/SE Method for Solving Maxwell's Equations in Time-Domain

An innovative finite-volume-type numerical method named as the space-time conservation element and solution element (CE/SE) method is applied to solve time-dependent Maxwell's equations in this paper. Test problems of electromagnetics scattering and antenna radiation are solved for validations. Numerical results are presented and compared with the analytical solutions, showing very good agreements.

Wang, X. Y.↗
NASA NTRS2004

Radiative Transfer Modeling of a Large Pool Fire by Discrete Ordinates, Discrete Transfer, Ray Tracing, Monte Carlo and Moment Methods

Five computational methods for solution of the radiative transfer equation in an absorbing-emitting and non-scattering gray medium were compared on a 2 m JP-8 pool fire. The temperature and absorption coefficient fields were taken from a synthetic fire due to the lack of a complete set of experimental data for fires of this size. These quantities were generated by a code that has been shown to agree well with the limited quantity of relevant data in the literature. Reference solutions to the governing equation were determined using the Monte Carlo method and a ray tracing scheme with high angular resolution. Solutions using the discrete transfer method, the discrete ordinate method (DOM) with both S(sub 4) and LC(sub 11) quadratures, and moment model using the M(sub 1) closure were compared to the reference solutions in both isotropic and anisotropic regions of the computational domain. DOM LC(sub 11) is shown to be the more accurate than the commonly used S(sub 4) quadrature technique, especially in anisotropic regions of the fire domain. This represents the first study where the M(sub 1) method was applied to a combustion problem occurring in a complex three-dimensional geometry. The M(sub 1) results agree well with other solution techniques, which is encouraging for future applications to similar problems since it is computationally the least expensive solution technique. Moreover, M(sub 1) results are comparable to DOM S(sub 4).

Jensen, K. A.↗
NASA NTRS2013

Measurement of Electromagnetic Energy Flow Through a Sparse Particulate Medium: A Perspective

First-principle analysis of the functional design of a well-collimated radiometer (WCR) reveals that in general, this instrument does not record the instantaneous directional flow of electromagnetic energy. Only in special cases can a sequence of measurements with a WCR yield the magnitude and direction of the local time-averaged Poynting vector. Our analysis demonstrates that it is imperative to clearly formulate the physical nature of the actual measurement afforded by a directional radiometer rather than presume desirable measurement capabilities. Only then can the directional radiometer be considered a legitimate part of physically based remote sensing and radiation-budget applications. We also emphasize the need for a better understanding of the nature of measurements with panoramic radiometers.

directional radiometry↗
NASA NTRS1998

Combining data from lidar and in situ instruments to characterize the vertical structure of aerosol optical properties

Over the last decade, the quantification of tropospheric aerosol abundance, composition and radiative impacts has become an important research endeavor. For the most part, the interest in tropospheric aerosols is derived from questions related to the global and local (instantaneous) radiative forcing of climate due to these aerosols. One approach is to study local forcing under well-defined conditions, and to extrapolate such results to global scales. To estimate local aerosol forcing, appropriate radiative transfer models can be employed (e.g., the Fu-Liou radiative transfer code, [Fu and Liou, 1993]). In general, such models require information on derived aerosol properties [Toon, 1994]; namely the aerosol optical depth, single-scattering albedo, and asymmetry factor (phase function), all of which appear in the equations of radiative transfer. In this paper, we report on a method that utilizes lidar data and in situ aerosol size distribution measurements to deduce the vertical structure of the aerosol complex index of refraction in the near IR, thus identifying the aerosol type. Together with aerosol size distributions obtained in situ, the aerosol refractive index can be used to calculate the necessary derived aerosol properties. The data analyzed here were collected during NASA's PEM West-B (Pacific Exploratory Mission) experiment, which took place in February/March 1994. The platform for the measurements was the NASA DC-8 aircraft. The primary goal of the PEM West missions [Browell et al., 1996] was the assessment of potential anthropogenic perturbations of the chemistry in the Pacific Basin troposphere. For this purpose the timing of PEM West-B corresponded to the seasonal peak in transport from the Asian continent into the Pacific basin [Merrill et al., in press]. This period normally occurs during Northern Hemisphere spring, when the Japan jet is well developed.

Redemann, J.↗
NASA NTRS1996

Particle Scattering in the Resonance Regime: Models Based on Integral Equations and Hybrid Finite Elements - Integral Equation Techniques

Reliable descriptions of the optical properties of clouds and aerosols are essential for studies of radiative transfer in the terrestrial atmosphere...Here we explore the utility of two approaches for deriving the single scattering optical properties of particles with sharp corners and large axial ratios.

clouds aerosols axi-symmetric particles resonance ↗
NASA NTRS1983

A model function for ocean microwave brightness temperatures

A relatively simple, yet accurate, relationship between the microwave brightness temperature of the ocean and conventional oceanographic and meteorological parameters is derived. The equation for the brightness temperature upwelling from the sea surface through the intervening atmosphere is obtained, considering radiative emission and scattering by the sea surface along with radiative absorption and emission by the atmosphere. A number of approximations are applied to the integral brightness temperature equation and its supporting equations in order to obtain a simple equation for the brightness temperature that does not contain integrals. Values for a number of atmospheric parameters are determined, including temperature sensitivities, oxygen opacity, water vapor and liquid water normalized absorption coefficients, and effective columnar height. The sea surface emissivity model is then considered, modelling the sea surface as a composite of foam-free rough water and foam patches.

Wentz, F. J.↗