A solution of a time-dependent equation of radiative transfer
Time dependent radiative transfer equation for plane-parallel isotropic scattering medium in first Gaussian approximation
Code, A. D.↗
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.
Modified radiative transfer theory for a two-layer anisotropic random medium
The modified radiative transfer (MRT) equations which describe propagation and scattering of the electromagnetic field intensity in a layered anisotropic random medium are derived from the Bethe-Salpeter equation with the ladder approximation and the Dyson equation with the nonlinear approximation. The extinction matrices and scattering phase matrices for both mean and incoherent field intensities are obtained. Backscattering enhancement is observed due to the spatial coherent effect between upward and downward propagating waves. It also occurs for the half-space case, because of coupling between ordinary and extraordinary waves. The depolarization effect is predicted in the first-order renormalization approximation to the MRT equation.
Lee, J. K.↗
Radiation and scattering from loaded microstrip antennas over a wide bandwidth
The integral equation and moment method solution is developed for two different antennas in the presence of an infinite grounded dielectric substrate. The first antenna is a rectangular microstrip patch antenna. This antenna is analyzed for excitation by an incident plane wave in free space and a vertical filament of uniform current in the dielectric. This antenna can be loaded by a lumped impedance in a vertical filament of uniform current extending from the patch through the dielectric to the ground plane. The radar cross section of the microstrip antenna is found from the plane wave excitation and shows good agreement to measurement for both an unloaded and loaded antenna. The input impedance is found from the current filament excitation. This is compared to the measured input impedance of a coaxially fed microstrip antenna and shows good agreement for both unloaded and loaded antennas when the dielectric substrate is much less than a wavelength. The second antenna is a vertical thin wire extending from the ground plane into or through the dielectric substrate. The mutual impedance between two imbedded monopoles is compared to a previous calculation.
Forrai, D. P.↗
Errors induced by the neglect of polarization in radiance calculations for Rayleigh-scattering atmospheres
Although neglecting polarization and replacing the rigorous vector radiative transfer equation by its approximate scalar counterpart has no physical background, it is a widely used simplification when the incident light is unpolarized and only the intensity of the reflected light is to be computed. We employ accurate vector and scalar multiple-scattering calculations to perform a systematic study of the errors induced by the neglect of polarization in radiance calculations for a homogeneous, plane-parallel Rayleigh-scattering atmosphere (with and without depolarization) above a Lambertian surface. Specifically, we calculate percent errors in the reflected intensity for various directions of light incidence and reflection, optical thicknesses of the atmosphere, single-scattering albedos, depolarization factors, and surface albedos. The numerical data displayed can be used to decide whether or not the scalar approximation may be employed depending on the parameters of the problem. We show that the errors decrease with increasing depolarization factor and/or increasing surface albedo. For conservative or nearly conservative scattering and small surface albedos, the errors are maximum at optical thicknesses of about 1. The calculated errors may be too large for some practical applications, and, therefore, rigorous vector calculations should be employed whenever possible. However, if approximate scalar calculations are used, we recommend to avoid geometries involving phase angles equal or close to 0 deg and 90 deg, where the errors are especially significant. We propose a theoretical explanation of the large vector/scalar differences in the case of Rayleigh scattering. According to this explanation, the differences are caused by the particular structure of the Rayleigh scattering matrix and come from lower-order (except first-order) light scattering paths involving right scattering angles and right-angle rotations of the scattering plane.
Mishchenko, M. I.↗
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.↗
Modifications Of Discrete Ordinate Method For Computations With High Scattering Anisotropy: Comparative Analysis
A numerical accuracy analysis of the radiative transfer equation (RTE) solution based on separation of the diffuse light field into anisotropic and smooth parts is presented. The analysis uses three different algorithms based on the discrete ordinate method (DOM). Two methods, DOMAS and DOM2+, that do not use the truncation of the phase function, are compared against the TMS-method. DOMAS and DOM2+ use the Small-Angle Modification of RTE and the single scattering term, respectively, as an anisotropic part. The TMS method uses Delta-M method for truncation of the phase function along with the single scattering correction. For reference, a standard discrete ordinate method, DOM, is also included in analysis. The obtained results for cases with high scattering anisotropy show that at low number of streams (16, 32) only DOMAS provides an accurate solution in the aureole area. Outside of the aureole, the convergence and accuracy of DOMAS, and TMS is found to be approximately similar: DOMAS was found more accurate in cases with coarse aerosol and liquid water cloud models, except low optical depth, while the TMS showed better results in case of ice cloud.
Korkin, Sergey V.↗
Angular Broadening: Effects of Nonzero, Spatially Varying Plasma Frequency Between the Source and Observer
Angular broadening of radiation due to scattering by density irregularities is usually described using geometric optics (GO) or the parabolic wave equation (PWE) with the assumptions that the radiation frequency f greatly exceeds the local plasma frequency f(sub p0) or that f(sub p0)/f is constant along the path. These assumptions are inappropriate for many solar system radio phenomena. Here the PWE and GO formalisms are extended to treat angular broadening in plasmas with nonzero, spatially varying ratios, f(sub p0)(z)/f < 1. The new PWE results show that the correlation function, scattered angular spectrum, and other quantities are modified by inclusion of a denominator factor [1 - f(sup 2, sub p0)(z prime)/f prime] inside the path integral over z prime, while the mean-square scattering angle depends on both f(sub p0)(z)/f at the observer and the foregoing factor inside the path integral. The PWE and GO predictions for are identical and involve equivalent assumptions. Previous GO and PWE results are recovered in the limits that f(sub p0)(z prime)/f is constant or zero. The new PWE and GO results will permit more accurate calculation of angular broadening for solar system and astrophysical sources. Moreover and importantly, due to the PWE and GO results for being identical, previous GO analyses of in solar system contexts are essentially correct, except for the neglect of or minor deficiencies in the treatment of nonzero, spatially varying f(sub p0)/f effects. The identical GO and PWE results for and the form of the PWE equation for the correlation function raise questions as to whether diffraction is unimportant for angular broadening (under the usual PWE conditions). Future direct comparisons of the PWE predictions with angular spectra calculated using existing GO ray-tracing codes should answer these questions. Diffraction effects are probably important when the medium and turbulence are not sufficiently homogeneous transverse to the central ray path. The implications are also discussed for studies of the 2-3.5 kHz radiation observed in the outer heliosphere.
Cairns, Iver H.↗
Two-Flux and Green's Function Method for Transient Radiative Transfer in a Semi-Transparent Layer
A method using a Green's function is developed for computing transient temperatures in a semitransparent layer by using the two-flux method coupled with the transient energy equation. Each boundary of the layer is exposed to a hot or cold radiative environment, and is heated or cooled by convection. The layer refractive index is larger than one, and the effect of internal reflections is included with the boundaries assumed diffuse. The analysis accounts for internal emission, absorption, heat conduction, and isotropic scattering. Spectrally dependent radiative properties are included, and transient results are given to illustrate two-band spectral behavior with optically thin and thick bands. Transient results using the present Green's function method are verified for a gray layer by comparison with a finite difference solution of the exact radiative transfer equations; excellent agreement is obtained. The present method requires only moderate computing times and incorporates isotropic scattering without additional complexity. Typical temperature distributions are given to illustrate application of the method by examining the effect of strong radiative heating on one side of a layer with convective cooling on the other side, and the interaction of strong convective heating with radiative cooling from the layer interior.
Siegel, Robert↗
Comparison of the Radiative Two-Flux and Diffusion Approximations
Approximate solutions are sometimes used to determine the heat transfer and temperatures in a semitransparent material in which conduction and thermal radiation are acting. A comparison of the Milne-Eddington two-flux approximation and the diffusion approximation for combined conduction and radiation heat transfer in a ceramic material was preformed to determine the accuracy of the diffusion solution. A plane gray semitransparent layer without a substrate and a non-gray semitransparent plane layer on an opaque substrate were considered. For the plane gray layer the material is semitransparent for all wavelengths and the scattering and absorption coefficients do not vary with wavelength. For the non-gray plane layer the material is semitransparent with constant absorption and scattering coefficients up to a specified wavelength. At higher wavelengths the non-gray plane layer is assumed to be opaque. The layers are heated on one side and cooled on the other by diffuse radiation and convection. The scattering and absorption coefficients were varied. The error in the diffusion approximation compared to the Milne-Eddington two flux approximation was obtained as a function of scattering coefficient and absorption coefficient. The percent difference in interface temperatures and heat flux through the layer obtained using the Milne-Eddington two-flux and diffusion approximations are presented as a function of scattering coefficient and absorption coefficient. The largest errors occur for high scattering and low absorption except for the back surface temperature of the plane gray layer where the error is also larger at low scattering and low absorption. It is shown that the accuracy of the diffusion approximation can be improved for some scattering and absorption conditions if a reflectance obtained from a Kubelka-Munk type two flux theory is used instead of a reflection obtained from the Fresnel equation. The Kubelka-Munk reflectance accounts for surface reflection and radiation scattered back by internal scattering sites while the Fresnel reflection only accounts for surface reflections.
Spuckler, Charles M.↗
Passive and active remote sensing of atmospheric precipitation
Both passive and active remote sensing of atmospheric precipitation are studied with the vector radiative transfer equations by making use of the Mie scattering phase functions and incorporating the raindrop-size distributions. For passive remote sensing, the Gaussian quadrature method is employed to solve for the brightness temperatures. For active remote sensing, an iterative approach carrying out to the second order in albedo is used to calculate for the bistatic coefficients, the backscattering cross sections/unit volume, and the interchannel cross talks. The calculated results are plotted as a function of rainfall rates and compared to various available experimental data. The theoretical model is easily applied to the remote sensing of aerosol particles, smoke, fog, and haze at infrared and visible frequencies.
Jin, Y. Q.↗
Methods for the solution of radiative transfer equation
To obtain an exact solution of the radiative-transfer equation in media where both absorption and scattering are significant, the usual approach is to use a numerical method. Three methods are known in the literature: invariant imbedding, eigenvalue-eigenfunction, and matrix doubling. This paper examines the practical application of these methods to the problem of emission from an inhomogeneous (Rayleigh) layer, the effects of layer parameters on the stability. It is found that invariant imbedding is most suitable for computing emission from an inhomogeneous layer with a temperature profile but tends to be unstable as the optical thickness of the layer increases beyond 0.5. On the other hand, the matrix-doubling method is stable for arbitrary optical thickness but is not suitable for handling multilayers. The eigenvalue-eigenfunction method is more stable than the invariant imbedding as optical thickness increases up to 2.0. It also permits temperature profile in the layer, but the computation is much more complicated. It is less stable than the matrix-doubling method when optical thickness is larger than 2.0. In general, the choice of a method is dependent on the nature of the problem.
Chen, M. F.↗
Sinusoidal reaction formulation for radiation and scattering from conducting surfaces
A piecewise-sinusoidal reaction technique is proposed for the study of scattering and radiation from perfectly conducting bodies of arbitrary shape. The reaction technique described is a moment method based on the integral equation of an electric field. A theory of sinusoidal reaction formulation for radiation and scattering problems is outlined, in which the scattered field is obtained by integrating the surface currents. Numerical results are presented for the scattering patterns of rectangular plates and corner reflectors and the radiation patterns of corner-reflector antennas.
Wang, N. N.↗
Derivation of phase functions from multiply scattered sunlight transmitted through a hazy atmosphere
The intensity of sunlight multiply scattered in model atmospheres is derived from the equation of radiative transfer by an analytical small-angle approximation. The approximate analytical solutions are compared to rigorous numerical solutions of the same problem. Results obtained from an aerosol-laden model atmosphere are presented. Agreement between the rigorous and the approximate solutions is found to be within a few per cent. The analytical solution to the problem which considers an aerosol-laden atmosphere is then inverted to yield a phase function which describes a single scattering event at small angles. The effect of noisy data on the derived phase function is discussed.
Weinman, J. A.↗
Scattering and absorption from poly-dispersed aerosols
Results from numerical computations of the diffuse sky radiances and degree of polarization are given. These results were obtained by solving the equation of radiative transfer in the formulation of Eschelbach. In these computations, multiple scattering and absorption by aerosol particles were considered. In addition, results from experimental measurements of the complex refractive index of aerosol samples at different relative humidities are presented. From the radiation flux divergences which were computed based on the determined properties of atmospheric aerosol particles, atmospheric heating rates were derived which were found to be comparable to the heating rates by water vapor.
Bullrich, K.↗
Three-dimensional radiative transfer using a Fourier-transform matrix-operator method
The three-dimensional equation of transfer for a scattering medium with planar geometry is solved by using a spatial Fourier transform and extending matrix-operator techniques developed previously for the one-dimensional equation. Doubling and adding algorithms were derived by means of an interaction principle for computing the Fourier-transformed radiation field. The resulting expressions fully describe the radiative transfer process in a scattering medium, inhomogeneous in the x-, y- and z-directions, illuminated from above by an arbitrarily general intensity field and bounded from below by a surface with completely general reflection properties.
Martonchik, J. V.↗
A Fast All-Sky Radiation Model for Solar Applications with Narrowband Irradiances on Tilted Surfaces (FARMS-NIT): Part II. The Cloudy-Sky Model
The Fast All-sky Radiation Model for Solar applications with Narrowband Irradiances on Tilted surfaces (FARMSNIT) reported in Part I of this study is enhanced to include the requirements for cloudy-sky conditions. Surface radiances in 2002 narrow-wavelength bands from 0.28 to 4.0 μm are analytically computed by solving the radiative transfer equation for five independent photon paths accounting for clear-sky absorption, Rayleigh scattering, and cloud absorption and scattering. The Simple Model of the Atmospheric Radiative Transfer of Sunshine (SMARTS) is used to provide the optical thickness of the clear-sky atmosphere. Unlike Part I, which approximates the computation of aerosol scattering using the single-scattering phase function, the cloud transmittance and reflectance are efficiently retrieved from a comprehensive look-up table pre-computed by a 32-stream DIScrete Ordinates Radiative Transfer (DISORT) model for possible cloud conditions as well as solar and viewing geometries. A resolution analysis is performed to assess the optimal balance between the computational efficiency and accuracy in the development of the look-up table. Model simulations by DISORT and TMYSPEC are used to evaluate the performance of FARMS-NIT under cloudy-sky conditions. Compared to DISORT, FARMS-NIT yields 2–3% uncertainties on average, but it substantially reduces the computational time because of the independent computation of cloud properties and the implementation of the look-up table. In contrast to TMYSPEC, which uses successive steps to empirically compute plane-of-array (POA) irradiances and spectral irradiances, FARMS-NIT directly solves spectral radiances from the radiative transfer equation, which profoundly increases the accuracy in surface irradiances, especially over inclined photovoltaics (PV) panels.
Xie, Yu↗
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.↗
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.↗