Exact formulation of multiple scattering in a three-dimensional cylindrical geometry

Exact expressions for the source function, flux, and scattered intensity normal to the surface are developed in cylindrical coordinates for a three-dimensional, absorbing, emitting, isotropically scattering medium exposed to both diffuse and collimated radiation. Simplifications of these expressions for certain important geometries and uniform loading are presented. Also, superposition of these equations and radiative equilibrium are discussed. The generalized three-dimensional equations are shown to reduce to the familiar one-dimensional results. Also, the equations for a strongly anisotropic phase function which is made up of a spike in the forward direction superimposed on an otherwise isotropic phase function are expressed in terms of the isotropic expressions.     

Radiative transfer in a two-dimensional rectangular medium exposed to diffuse radiation

The integral equation for radiative transfer in a two-dimensional rectangular scattering medium exposed to diffuse radiation is solved numerically by removing the singularity. This method yielded accurate results except at very large optical thicknesses. Graphical and tabular results for the source function, flux, and intensity are presented. The source function is also calculated using the first term of a Taylor series expansion. The Taylor series is fairly accurate for small optical thicknesses and columnar geometries. A method is presented for extending these results to the problem of a strongly anisotropic scattering phase function which is made up of a spike in the forward direction superimposed on an isotropic phase function.     

Two-dimensional isotropic scattering in a semi-infinite cylindrical medium

Beginning with the integral equation for the source function, the solutions for the source function, flux and intensity at the boundary of a two-dimensional, isotropically scattering cylindrical medium are found. The incident radiation is collimated and normal to the surface of the medium and depends only on the radial coordinate. For a Bessel function boundary condition, separation of variables is used to reduce the source function integral equation to a one-dimensional equation. The resulting integral equation is shown to be the same as that for the two-dimensional planar case. Solutions for other boundary conditions are then shown to be superpositions of the Bessel function solution. Numerical results are presented for a Gaussian distribution of incident radiation which closely models a laser beam. These multiple scattering results are compared to the single scattering approximation. Also, the solution for a strongly anisotropic phase function which is made up of a spike in the forward direction superimposed on an otherwise isotropic phase function is expressed in terms of the isotropic results.     

Two-dimensional linearly anisotropic scattering in a semi-infinite cylindrical medium exposed to a laser beam

A modification of Ambarzumian's method is used to develop the integro-differential equations for the source function, flux, and intensity at the boundary of a two-dimensional, semi-infinite cylindrical medium which scatters linearly. The incident radiation is collimated, normal to the top surface of the medium, and is dependent only on the radial coordinate. The radial variation is assumed to be a Bessel function or a Gaussian distribution. The Gaussian boundary condition is used to simulate a laser beam. Numerical results are presented in graphical and tabular forms for both boundary conditions. Results for forward and backward scattering phase functions are compared with those for isotropic scattering. A method is presented for extending these results to the problem of a strongly anisotropic phase function which is made up of a spike in the forward direction superimposed on a linear phase function.     

Two-dimensional rayleigh scattering in a semi-infinite cylindrical medium exposed to a laser beam

A modification of Ambarzumian's method is used to develop the integro-differential equations for the source function, flux, and intensity at the boundary of a two-dimensional, semi-infinite cylindrical medium with second order Legendre phase function scattering. The incident radiation is collimated, normal to the top surface, and is dependent only on the radial coordinate. Boundary conditions which vary as a Bessel function and as a Gaussian distribution are investigated. The Gaussian distribution approximates a laser beam. Numerical results are presented in graphical and tabular forms for a Rayleigh scattering medium. The results are compared with those of isotropic scattering.     

Polarization of multiple scattered light in planetary atmospheres from solution of the coupled linear integral equations derived for mixed rayleigh-isotropic scattering

A set of coupled integral equations describing nonconservative multiple scattering for a mixed isotropic and Rayleigh single scattering phase function in inhomogeneous, plane-parallel planetary atmospheres is derived. The equations are applicable for the frequency redistributions MRE (monochromatic radiative equilibrium), CFR (complete frequency redistribution), or PFR (partial frequency redistribution). Solution of the equations permits one to calculate the intensity and degree of polarization in an arbitrary direction outside or inside the plane-parallel scattering medium. The equations are readily adaptable to more complicated geometries. Solutions for several cases are presented to demonstrate the versatility and validity of the method. These include a calculation of MRE pure Rayleigh scattering of sunlight in an optically-thick planetary atmosphere to demonstrate agreement with the results of Coulson⁽¹⁾et al. (1960), calculations of the altitude profile of the degree of polarization of the earth's Ly-α 1216A and helium 584A dayglow, and the center to limb variation of the degree of polarization of the sunlight diffusely reflected from a distant planet such that the scattering is pure Rayleigh and conservative at the top of the planet's atmosphere, varying smoothly to conditions of nonconservative and pure isotropic scattering deep in the planet's atmosphere. Tables of functions that one might utilize (without resort to a digital computer) to obtain solutions in the escape function approximation are also given.     

Two-dimensional isotropic scattering in a finite thick cylindrical medium exposed to a laser beam

Graphical and tabular results are presented for the back-scattered intensity from a finite two-dimensional cylindrical medium exposed to a Gaussian beam of radiation. Also, results for the source function and flux at the boundaries are presented. The influence of optical thickness and albedo are most pronounced at large optical radii. The semi-infinite results can be used to approximate the finite case for small optical radii. Ranges for single, double, and multiple scattering are discussed. For locations far from the incident beam, the results can be expressed in terms of universal functions independent of beam size. A method is presented for extending the isotropic results to the anisotropic case where the phase function is made up of a spike superimposed on an otherwise isotropic phase function.     

Solution of the kinetic equation for thin two-and three-dimensional metallic slabs with diffuse boundary scattering and arbitrary impurity scattering in the bulk

The integrodifferential equation in the phase space, describing the response to a homogeneous electric field of two-or three-dimensional isotropic metallic slab with diffuse scattering on the border and arbitrary scattering on impurities, is transformed into the system of inhomogeneous integral equations in the configuration space. It is shown that one of the functions fulfilling this system is simply connected with the current response of the slab. The general solution of this system is found in the form of asymptotical expansion provided that the mean-free path is much greater than the width of the slab. In the opposite limit, the asymptotic solution is found only at almost isotropic impurity scattering. In other words, the Fuchs problem (1938) has been solved for thin slabs at a real low-temperature collision integral.     

Spectral diffuse reflectance measurements of gadolinia, silica thin film systems using synchrotron radiation

Diffuse reflectance or optical scattering in thin films and multilayers can pose serious limiting factors to their desired or ultimate performances. Besides, such studies provide valuable information related to the buried microstructures and interfaces. Synchrotron radiation is the most appropriate source to record wavelength dependent polarized light scattering in thin films and multilayers. In the present experiment several gadolinia, silica thin films and multilayers were studied for their light scattering using the white light synchrotron beam. Various thin film layer geometries were selected to probe the results due to different types and combinations of interfaces. Due to phase coherent delay in certain optical non-wedged component used in the experimental setup very interesting spectral interference were noticed as the modulations in the diffuse reflectance signal. Appropriate modeling approach utilizing Gaussian function de-convolution technique is used to compute the pulse delay between the back reflected and forward propagating scattering signals that lead to such spectral interference. Alternatively inverse fast Fourier transform (IFFT) and analytical techniques were adopted to determine the group delays. The diffuse reflectance spectra were finally compared with their specular counter part and a shift both in the wavelength and phase were noticed. This may be explained on the basis of thin film roughness factors, different polarizations and incident geometries used in the measurements. However, all most all the spectral scattering signals exhibited the features resembling their specular reflection characteristic.     

Diffuse reflectance of TiO2 pigmented paints: Spectral dependence of the average pathlength parameter and the forward scattering ratio

Diffuse reflectance spectra of paint coatings with different pigment concentrations, normally illuminated with unpolarized radiation, have been measured. A four-flux radiative transfer approach is used to model the diffuse reflectance of TiO₂ (rutile) pigmented coatings through the solar spectral range. The spectral dependence of the average pathlength parameter and of the forward scattering ratio for diffuse radiation, are explicitly incorporated into this four-flux model from two novel approximations. The size distribution of the pigments has been taken into account to obtain the averages of the four-flux parameters: scattering and absorption cross sections, forward scattering ratios for collimated and isotropic diffuse radiation, and coefficients involved in the expansion of the single particle phase function in terms of Legendre polynomials.     

Radiative transfer with anisotropic scattering in an isothermal slab

Two types of anisotropic scattering, linear anisotropic scattering and Rayleigh anisotropic scattering, are considered in the analysis of radiative transfer for an isothermal, plane-parallel medium confined between gray, diffuse walls. The problem is formulated in terms of a coupled pair of integral equations containing the intensity-moments as the unknown variables. These intensity-moments are shown to be the components of the source function. The set of equations is then solved both numerically and in closed form. The results reveal clearly the effects of anisotropic scattering on important characteristics such as heat flux directional emittance and incident radiant energy per unit area. These effects are well predicted by the approximate closed-form solution.     

Application of the mutual-interaction method to a class of two-scatterer systems: I. Two discrete scatterers

In a recent paper we developed a formalism that fully accommodates the mutual interactions among scatterers separable by parallel planes. The total fields propagating away from these planes are the unknowns of a system of difference equations. Each scatterer is characterized by a scattering function that expresses the scattered wave amplitude as a function of the incident and scattered wavevectors for a unit-amplitude plane wave scattered from the object in isolation. This function can be derived completely from the scattered far field with the help of analytic continuation. For a two-scatterer system the mutual-interaction equations reduce to a single Fredholm integral equation of the second kind. It turns out that analytic solutions are tractable for those scattering functions that are Dirac deltas or a sum of products of separable functions of the incident and scattered wavevectors. Scattering functions for planes and isotropic scatterers, as well as electric and magnetic dipoles all possess this property and are considered. The exact scattering functions agree with results obtained by analytic continuation. This paper consists of two parts. Part I derives analytic solutions for two discrete scatterers (isotropic scatterers. electric dipoles, magnetic dipoles). Part II is devoted to scattering from an object (isotropic or dipole scatterer) near an interface separating two semi-infinite uniforn-media. Because the results in this paper are exact, the effects of near-field interactions can be assessed. The forms of the scattering solutions can be adapted to objects that are both radiating and scattering.     

Application of the mutual-interaction method to a class of two-scatterer systems: I. Two discrete scatterers

Abstract

In a recent paper we developed a formalism that fully accommodates the mutual interactions among scatterers separable by parallel planes. The total fields propagating away from these planes are the unknowns of a system of difference equations. Each scatterer is characterized by a scattering function that expresses the scattered wave amplitude as a function of the incident and scattered wavevectors for a unit-amplitude plane wave scattered from the object in isolation. This function can be derived completely from the scattered far field with the help of analytic continuation. For a two-scatterer system the mutual-interaction equations reduce to a single Fredholm integral equation of the second kind. It turns out that analytic solutions are tractable for those scattering functions that are Dirac deltas or a sum of products of separable functions of the incident and scattered wavevectors. Scattering functions for planes and isotropic scatterers, as well as electric and magnetic dipoles all possess this property and are considered. The exact scattering functions agree with results obtained by analytic continuation. This paper consists of two parts. Part I derives analytic solutions for two discrete scatterers (isotropic scatterers. electric dipoles, magnetic dipoles). Part II is devoted to scattering from an object (isotropic or dipole scatterer) near an interface separating two semi-infinite uniforn-media. Because the results in this paper are exact, the effects of near-field interactions can be assessed. The forms of the scattering solutions can be adapted to objects that are both radiating and scattering.     

Two-dimensional radiative transfer in a finite scattering planar medium

The source function, radiative flux, and intensity at the boundaries are calculated for a two-dimensional, scattering, finite medium subjected to collimated radiation. The scattering phase function is composed of a spike in the forward direction super-imposed on an isotropic background. Exact radiative transfer theory is used to formulate the problem and Ambarzumian's method is used to obtain results. Using the principle of superposition, the results for any step variation in incident radiation are expressed in terms of universal functions for the semi-infinite step case. Two-dimensional effects are most pronounced at large optical thicknesses and albedos.     

Equivalence of MTF of a turbid medium and radiative transfer field

The equivalence of the modulation transfer function (MTF) of a turbid medium and the transmitted radiance from the medium under isotropic diffuse illumination is demonstrated. MTF of a turbid medium can be fully evaluated by numerically solving a radiative transfer problem in a plane parallel medium. MTF for a homogenous single layer turbid medium is investigated as illustration. General features of the MTF in the low and high spatial frequency domains are provided through their dependence on optical thickness, single scattering albedo, asymmetrical factor, and phase function type.     

Temperature response in participating media with anisotropic scattering caused by pulsed lasers

It is not by isotropic scattering but by anisotropic scattering that radiant energy is redistributed in some materials containing real particles, fibers, or impurities. In some instances, great difference can be caused in transient thermal behavior between isotropic scattering and anisotropic scattering media. Ray tracing method combined with Hottel's zonal method is introduced to deduce thermal radiative source term for various optical boundary conditions induced by collimated incidence passing through translucent boundary. Temperature response caused by laser pulse at non-incident side of participating and anisotropic scattering media is examined. We investigate effects of scattering albedo, scattering phase function, initial temperature of media and thickness of media on temperature response. Results obtained for anisotropic scattering media are compared with those for isotropic scattering one and show that anisotropic scattering must be considered in the simulating measurement of thermophysical properties by the laser flash method for some materials with big scattering albedo which behave anisotropically, or big error will be introduced; forward scattering can increase excess temperature and backward scattering can decrease it at non-incident side of the considered sample irradiated by laser pulse.     

Analysis of collimated radiation in participating media using the discrete transfer method

A general formulation of the discrete transfer method is provided to analyze radiative heat transfer problems in a participating medium subjected to collimated radiation. The formulation is validated by considering 1-D planar absorbing, emitting and anisotropically scattering gray medium in radiative equilibrium. Anisotropy of the medium is approximated by linear anisotropic phase function. For the purpose of comparison, the problem is also solved analytically. Results are obtained for different angles of incidence of the collimated radiation. At a given angle of incidence, results are obtained for forward, isotropic and backward scattering situations. Heat flux results are compared over a wide range of values of the extinction coefficient. Emissive power distributions in the medium are also obtained for some cases. The discrete transfer method results are found to compare very well with the analytic results.     

Pure-triplet scattering for radiative transfer in binary discrete random finite media with reflective boundaries and internal source of radiation

The stochastic solution of the monoenergetic radiative transfer equation in a finite slab random medium with pure-triplet anisotropic scattering is considered. The random medium is assumed to consist of two randomly mixed immiscible fluids labelled by 1 and 2. The extinction function, the scattering kernel, and the internal source of radiation are treated as discrete random variables, which obey the same statistics. The theoretical model used here for stochastic media transport assumes Markovian processes and exponential chord length statistics. The boundaries of the medium under consideration are considered to have specular and diffuse reflectivities with an internal source of radiation inside the medium. The ensemble-average partial heat fluxes are obtained in terms of the average albedos of the corresponding source-free problem, whose solution is obtained by using the Pomraning-Eddington approximation. Numerical results are calculated for the average forward and backward partial heat fluxes for different values of the single scattering albedo with variation of the parameters that characterize the random medium. Compared to the results obtained by Adams et al. in the case of isotropic scattering based on the Monte Carlo technique, it can be demonstrated that we have good comparable data.     

A comparison of variational and discrete ordinate methods for solving radiative transfer problems

In this paper, we compare the variational and discrete ordinate methods applied to their simplest representative cases: the 2-step function variational method and the 2-stream approach. It is shown that both methods require minimal computational time. Numerical results of the source functions, intensities and fluxes for the problem of diffuse reflection by different plane-parallel isotropic atmospheres are compared between the two methods. Limitations of both approaches for different single scattering albedos and optical depths are discussed. It is also shown that the variational technique is a more efficient algorithm in handling multi-layer cloud conditions for radiative transfer problems.     

Two-dimensional linearly anisotropic scattering in a finite-thick cylindrical medium exposed to a laser beam

Integro-differential equations are developed for the source function, flux and intensity at the boundaries of a two-dimensional finite-thick medium which scatters in a linear fashion. The incident radiation is collimated, normal to the upper surface of the medium and dependent only on the radial coordinate. Two radial distributions are investigated: (1) a Bessel function and (2) a Gaussian laser beam. The solution for the Gaussian beam is constructed from the Bessel solution. Numerical results are presented in graphical and tabular forms for both boundary conditions. Comparisons are made between forward and isotropic scattering and between the finite and semi-infinite cases.