SEARCH · Search NASA
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.
Radiative transfer equation solution for spectral line formed in two dimensionally varying atmosphere extended to continuum radiation in inhomogeneous atmospheres
An understanding of radiative transfer in the earth's atmosphere is a necessity for the remote sensing of surface reflectivity from satellites and aircraft. The range of the adjacency effect, which represents the effect of bright areas on the radiance above dark areas, is the main parameter that distinguishes atmospheric radiative transfer over a nonuniform surface from that over a uniform one. A radiative transfer model which expresses this range correctly is, therefore, critical for developing remote sensing methods for the case of an atmosphere over a nonuniform surface. The present investigation is concerned with the development of a new approximate solution of the radiative transfer (RT) equation. The solution is not limited to nonabsorbing atmospheres, but it will still be limited to nadir observations. The results compare favorably with Monte Carlo simulations.
Recent studies indicate that a cloudy atmosphere absorbs more solar radiation than any current 1D or 3D radiation model can predict. The excess absorption is not large, perhaps 10-15 W/sq m or less, but any such systematic bias is of concern since radiative transfer models are assumed to be sufficiently accurate for remote sensing applications and climate modeling. The most natural explanation would be that models do not capture real 3D cloud structure and, as a consequence, their photon path lengths are too short. However, extensive calculations, using increasingly realistic 3D cloud structures, failed to produce photon paths long enough to explain the excess absorption. Other possible explanations have also been unsuccessful so, at this point, conventional models seem to offer no solution to this puzzle. The weakest link in conventional models is the way a size distribution of cloud particles is mathematically handled. Basically, real particles are replaced with a single average particle. This "ensemble assumption" assumes that all particle sizes are well represented in any given elementary volume. But the concentration of larger particles can be so low that this assumption is significantly violated. We show how a different mathematical route, using the concept of a cumulative distribution, avoids the ensemble assumption. The cumulative distribution has jumps, or steps, corresponding to the rarer sizes. These jumps result in an additional term, a kind of Green's function, in the solution of the radiative transfer equation. Solving the cloud radiative transfer equation with the measured particle distributions, described in a cumulative rather than an ensemble fashion, may lead to increased cloud absorption of the magnitude observed.
A discrete probability function (DPF) method for the equation of radiative transfer is derived. The DPF is defined as the integral of the probability density function (PDF) over a discrete interval. The derivation allows the evaluation of the PDF of intensities leaving desired radiation paths including turbulence-radiation interactions without the use of computer intensive stochastic methods. The DPF method has a distinct advantage over conventional PDF methods since the creation of a partial differential equation from the equation of transfer is avoided. Further, convergence of all moments of intensity is guaranteed at the basic level of simulation unlike the stochastic method where the number of realizations for convergence of higher order moments increases rapidly. The DPF method is described for a representative path with approximately integral-length scale-sized spatial discretization. The results show good agreement with measurements in a propylene/air flame except for the effects of intermittency resulting from highly correlated realizations. The method can be extended to the treatment of spatial correlations as described in the Appendix. However, information regarding spatial correlations in turbulent flames is needed prior to the execution of this extension.
The present volume on radar polarimetry for geoscience applications discusses wave properties and polarization, scattering matrix representation for simple targets, scattering models for point and distributed targets, polarimetric scatterometer systems and measurements, polarimetric radar system design, and polarimetric SAR applications. Attention is given to plane waves in a lossless homogeneous medium-wave polarization, polarization synthesis and response, and coordinate system transformations. Topics addressed include high- and low-frequency scattering, rough-surface scattering models, radiative transfer theory and deficiencies thereof, solutions for the radiative transfer equation, and a radiative transfer model for a forest canopy. Also discussed are network analyzer-based polarimetric scatterometers, calibration of polarimetric scatterometers, synthesized polarization response of distributed targets, and measurement of the propagation parameters of a forest canopy.
The derivation of the equations of statistical equilibrium are outlined, starting from the quantum density-matrix equations, drawing particular attention to the approximations and assumptions used in the development of tractable expressions. Then, using the quantum-fluctuation-regression theorem, emission and absorption coefficients are obtained for multilevel atomic systems which are nondegenerate except for m-substates. These coefficients are valid to first order in the incident intensity. Possible extensions to higher intensity broadband incoherent fields are suggested.
An approximation to the radiative transfer equation for solar radiation in relatively full, homogeneous plant canopies is presented and solved analytically for solar zenith angles less than 60 deg. The model predicts reflectance at any depth in the canopy and in any direction and may be inverted with bidirectional reflectance measurements. The model was fit to data at two sun angles and two wavebands (visible and NIR) to within the assumed errors on the reflectance data. The calculated albedos are insensitive to achievable measurement errors. Some of the parameter values themselves found by the inversion agree reasonably well with independent measurements, but the uncertainties introduced by the data noise are rather large. However, the agreement is good enough to demonstrate that the model is physically realistic.
Multiple scattering contributions to lidar returns from turbid atmospheres are derived by means of an analytical theory. It is assumed that scattering takes place mainly at small angles except for one event that scatters the light backward. The phase functions are approximated by the sum of Gaussian functions of the scattering angle in both the forward and backward directions. The three-dimensional radiative transfer equation is transformed to a one-dimensional problem by means of Fourier transforms. Neumann solutions to the transformed equation of radiative transfer are then found. A number of examples are presented for cloud, fog and haze models. The results are found to be in satisfactory agreement with results obtained from the Monte Carlo analysis of Kunkel (1974) and the theory of light pulses doubly scattered by turbid atmospheres which was developed by Eloranta (1972).
The recent focus of parameterization of the radiative properties of clouds has been to include the microphysical properties of the cloud. A variety of parameterization have been developed for both the shortwave and the longwave. In parameterizing the longwave properties of clouds, it is useful to consider the two stream solution of the radiative transfer equation appropriate for a thermal source. These radiative transfer equations are considered.
Numerical solutions for radiative transport in a class of anisotropically scattering materials are presented. Conditions for convergence and divergence of the iterative method are given and supported by computed results. The relation of two flux theories to the equation of radiative transfer for isotropic scattering is discussed. The adequacy of the two flux approach for the reflectance, radiative flux and radiative flux divergence of highly scattering media is evaluated with respect to solutions of the radiative transfer equation.
In this paper, three approximate physical microwave radiometry models have been used to calculate brightness temperatures from a forest canopy at L-band. These models are (1) tau-omega model (zero order scattering approximation to radiative transfer equations), (2) successive order of scattering model up to first order (first order scattering approximation to the radiative transfer equations), and (3) Peak technique utilizing the active solution obtained from the Distorted Born Approximation (DBA). These models are physically-based and treat vegetation as a layer of discrete scatterers over a rough surface. Vegetation components within the canopy are represented by canonical shapes such as dielectric discs and cylinders. The tau-omega model is based on a zero-order solution to the radiative transfer (RT) equations. The model ignores scattering except for the effect of the scatterers in the attenuation of the emission through the vegetation. Application of the tau-omega model to data acquired during airborne and ground-based campaigns over the years has solidified scientific understanding of microwave interactions with different landscapes. In particular, shrubland, grasslands, agricultural crops, and light to moderate vegetation have been investigated. Its applicability to areas with a significant tree fraction is unknown. The first order scattering model is based on an iterative solution of the RT equation up to the first order. The first order solution is obtained by substituting the zeroth-order solution into the scattering source term and then solving the resulting radiative transfer equations. This formulation adds a new scattering term to the tau-omega model. It represents emission by particles in the layer and emission by the ground that is scattered once by particles in the layer. The resulting model represents an improvement over the standard zero-order solution (the tau-omega model) since it accounts for the scattered vegetation and ground radiation that can have a pronounced effect on the observed brightness temperature. The third model is based on the Peake formulation in conjunction with the DBA. The procedure for calculation of forest emission is accomplished by first calculating the bistatic scattering cross section for each type of scatterer, then by using the DBA to calculate specular albedo of the ground and the diffused albedo of the layer. Once the albedos are determined, Peake s principle relating active and passive problems can be used to determine the effective emissivity of the forest layer.
Spectral reflectance of snow under diffuse illumination is studied using the two-stream approximation of the radiative transfer equation. The scattering and absorption parameters of the radiative transfer equation - the single scattering albedo, the optical depth, and the integrated phase function are obtained from the grain size and density of snow. Analytical expressions for the intensity within the snowpack, the reflectance, and the asymptotic flux extinction coefficient, are given. Good agreement is shown between the theory and available experimental data on visible and near-infrared reflectance, and the asymptotic flux extinction coefficient. The theory may also be used to explain the observed effect of aging on the snow reflectance.
In a prior effort, a thermal/fluid model of the interior of Penn State University's laboratory-scale Insulation Test Motor (ITM) was constructed to predict both the convective and radiative heat transfer to the interior walls of the ITM with a minimum of empiricism. These predictions were then compared to values of total and radiative heat flux measured in a previous series of ITM test firings to assess the capabilities and shortcomings of the chosen modeling approach. Though the calculated fluxes reasonably agreed with those measured during testing, this exercise revealed means of improving the fidelity of the model to, in the case of the thermal radiation, enable direct comparison of the measured and calculated fluxes and, for the total heat flux, compute a value indicative of the average measured condition. By replacing the P1-Approximation with the discrete ordinates (DO) model for the solution of the gray radiative transfer equation, the radiation intensity field in the optically thin region near the radiometer is accurately estimated, allowing the thermal radiation flux to be calculated on the heat-flux sensor itself, which was then compared directly to the measured values. Though the fully coupling the wall thermal response with the flow model was not attempted due to the excessive computational time required, a separate wall thermal response model was used to better estimate the average temperature of the graphite surfaces upstream of the heat flux gauges and improve the accuracy of both the total and radiative heat flux computations. The success of this modeling approach increases confidence in the ability of state-of-the-art thermal and fluid modeling to accurately predict SRM internal environments, offers corrections to older methods, and supplies a tool for further studies of the dynamics of SRM interiors.
In a prior effort, a thermal/fluid model of the interior of Penn State University's laboratory-scale Insulation Test Motor (ITM) was constructed to predict both the convective and radiative heat transfer to the interior walls of the ITM with a minimum of empiricism. These predictions were then compared to values of total and radiative heat flux measured in a previous series of ITM test firings to assess the capabilities and shortcomings of the chosen modeling approach. Though the calculated fluxes reasonably agreed with those measured during testing, this exercise revealed means of improving the fidelity of the model to, in the case of the thermal radiation, enable direct comparison of the measured and calculated fluxes and, for the total heat flux, compute a value indicative of the average measured condition. By replacing the P1-Approximation with the discrete ordinates (DO) model for the solution of the gray radiative transfer equation, the radiation intensity field in the optically thin region near the radiometer is accurately estimated, allowing the thermal radiation flux to be calculated on the heat-flux sensor itself, which was then compared directly to the measured values. Though the fully coupling the wall thermal response with the flow model was not attempted due to the excessive computational time required, a separate wall thermal response model was used to better estimate the average temperature of the graphite surfaces upstream of the heat flux gauges and improve the accuracy of both the total and radiative heat flux computations. The success of this modeling approach increases confidence in the ability of state-of-the-art thermal and fluid modeling to accurately predict SRM internal environments, offers corrections to older methods, and supplies a tool for further studies of the dynamics of SRM interiors.
An effective iterative technique is introduced to solve a nonlinear integral equation frequently associated with radiative transfer problems. The problem is formulated in such a way that each step of an iterative sequence requires the solution of a linear integral equation. The advantage of a previously introduced variational technique which utilizes a stepwise constant trial function is exploited to cope with the nonlinear problem. The method is simple and straightforward. Rapid convergence is obtained by employing a linear interpolation of the iterative solutions. Using absorption coefficients of the Milne-Eddington type, which are applicable to some planetary atmospheric radiation problems. Solutions are found in terms of temperature and radiative flux. These solutions are presented numerically and show excellent agreement with other numerical solutions.
My presentation will begin with the discussion of the Intercomparison of three-dimensional (3D) Radiative Codes (13RC) project that has been started in 1997. I will highlight the question of how well the atmospheric science community can solve the 3D radiative transfer equation. Initially I3RC was focused only on algorithm intercomparison; now it has acquired a broader identity providing new insights and creating new community resources for 3D radiative transfer calculations. Then I will switch to satellite remote sensing. Almost all radiative transfer calculations for satellite remote sensing are one-dimensional (1D) assuming (i) no variability inside a satellite pixel and (ii) no radiative interactions between pixels. The assumptions behind the 1D approach will be checked using cloud and aerosol data measured by the MODerate Resolution Imaging Spectroradiometer (MODIS) on board of two NASA satellites TERRA and AQUA. In the discussion, I will use both analysis technique: statistical analysis over large areas and time intervals, and single scene analysis to validate how well the 1D radiative transfer equation describes radiative regime in cloudy atmospheres.
A reconstruction technique for reducing computation burden in the 3D image processes, wherein the reconstruction procedure comprises an inverse and a forward model. The inverse model uses a hybrid dual Fourier algorithm that combines a 2D Fourier inversion with a 1D matrix inversion to thereby provide high-speed inverse computations. The inverse algorithm uses a hybrid transfer to provide fast Fourier inversion for data of multiple sources and multiple detectors. The forward model is based on an analytical cumulant solution of a radiative transfer equation. The accurate analytical form of the solution to the radiative transfer equation provides an efficient formalism for fast computation of the forward model.
Complete polarimetric signatures of a canopy of dielectric cylinders overlying a homogeneous half space are studied with the first and second order solutions of the vector radiative transfer theory. The vector radiative transfer equations contain a general nondiagonal extinction matrix and a phase matrix. The energy conservation issue is addressed by calculating the elements of the extinction matrix and the elements of the phase matrix in a manner that is consistent with energy conservation. Two methods are used. In the first method, the surface fields and the internal fields of the dielectric cylinder are calculated by using the fields of an infinite cylinder. The phase matrix is calculated and the extinction matrix is calculated by summing the absorption and scattering to ensure energy conservation. In the second method, the method of moments is used to calculate the elements of the extinction and phase matrices. The Mueller matrix based on the first order and second order multiple scattering solutions of the vector radiative transfer equation are calculated. Results from the two methods are compared. The vector radiative transfer equations, combined with the solution based on method of moments, obey both energy conservation and reciprocity. The polarimetric signatures, copolarized and depolarized return, degree of polarization, and phase differences are studied as a function of the orientation, sizes, and dielectric properties of the cylinders. It is shown that second order scattering is generally important for vegetation canopy at C band and can be important at L band for some cases.