Pomraning-Eddington approximation for radiative transfer in a spherical turbid medium

Abstract

The Pomraning-Eddington approximation is used to solve the radiative transfer problem for anisotropic scattering in a spherical homogeneous turbid medium with diffuse and specular reflecting boundaries. This approximation replaces the radiative transfer integro-differential equation by a second-order differential equation which has an analytical solution in terms of the modified Bessel function. Here, we calculate the partial heat flux at the boundary of anisotropic scattering on a homogeneous solid sphere. The calculations are carried out for spherical media of radii 0.1, 1.0 and 10 mfp and for scattering albedos between 0.1 and 1.0. In addition, the calculations are given for media with transparent, diffuse reflecting and diffuse and specular reflecting boundaries. Two different weight functions are used to verify the boundary conditions. Our results are compared with those given by the Galerkin technique and show greater accuracy for thick and highly scattering media.     

Effect of reflecting modes on combined heat transfer within an anisotropic scattering slab

Under various interface reflecting modes, different transient thermal responses will occur in the media. Combined radiative-conductive heat transfer is investigated within a participating, anisotropic scattering gray planar slab. The two interfaces of the slab are considered to be diffuse and semitransparent. Using the ray tracing method, an anisotropic scattering radiative transfer model for diffuse reflection at boundaries is set up, and with the help of direct radiative transfer coefficients, corresponding radiative transfer coefficients (RTCs) are deduced. RTCs are used to calculate the radiative source term in energy equation. Transient energy equation is solved by the full implicit control-volume method under the external radiative-convective boundary conditions. The influences of two reflecting modes including both specular reflection and diffuse reflection on transient temperature fields and steady heat flux are examined. According to numerical results obtained in this paper, it is found that there exits great difference in thermal behavior between slabs with diffuse interfaces and that with specular interfaces for slabs with big refractive index.     

Heat transfer in a spherical turbid medium with conduction and radiation

The heat transfer through a spherical media with conduction and radiation is considered. The medium is considered to be turbid and anisotropically scattering with diffusely reflecting boundaries of constant temperatures. The radiative transfer problem is solved using the Galerkin method. An iterative method is used to solve the nonlinear relation between the radiative transfer equation and the conductive energy equation. Calculations are carried out and compared for a homogeneous, isotropically scattering medium with isothermal, transparent boundaries. The results show good agreement with previous work. Calculations are carried out for inhomogeneous media with isotropic, and forward and backward anisotropic scattering. The boundaries of the media are considered to be isothermal and may be transparent or diffusely reflecting boundaries. The calculations are used to study the effects of the single scattering albedo, the anisotropic scattering parameter, the conduction-radiation parameter and the heat source.     

Pomraning–Eddington approximation for radiative transfer in a homogeneous solid cylinder

Abstract

The radiative transfer in a solid cylinder containing a homogeneous turbid medium with anisotropic scattering is considered. The medium has a diffuse and specular reflecting boundary illuminated by an external incidence and contains an internal energy source. This general problem can be solved in terms of the solution of the corresponding source-free problem with a specular reflecting boundary and isotropic external incidence. The Pomraning–Eddington approximation is used to solve the source-free problem. Three different weight functions are used to verify the boundary condition to find the constants of the solution. The partial flux, the irradiance and the net flux at the boundary for the general problem are calculated.     

Stochastic radiative transfer in a finite plane-parallel medium with general boundary conditions

The time-independent radiative transfer problem in a scattering and absorbing planar random medium with general boundary conditions and internal energy source is considered. The medium is assumed to consist of two randomly mixed immiscible fluids, with the mixing statistics described as a two-state homogeneous Markov process. The problem is solved in terms of the solution of the corresponding free-source problem with simple boundary conditions which is solved using Pomraning-Eddington approximation in the deterministic case. A formalism, developed to treat radiative transfer in statistical mixtures, is used to obtain the ensemble-averaged solution. The average partial heat fluxes are calculated in terms of the albedoes of the source-free problem. Results are obtained for isotropic and anisotropic scattering for specular and diffused reflecting boundaries.     

Coupled conductive-radiative heat transfer problem for two-layer slab

The coupled conductive radiative transfer problem in two homogeneous layers slab of anisotropic scattering with specularly reflecting boundaries has been considered. A Galerkin-iterative technique is used to solve the coupled conductive radiative heat equations in integral forms for the two layers. Numerical results are obtained for the temperature, the conductive, radiative and the total heat fluxes for the two homogeneous layers with isotropic and anisotropic scattering. The calculations are also carried out for homogeneous plane parallel medium with anisotropic scattering which show good agreement with the published calculations.     

The Iterative-DRESOR method to solve radiative transfer in a plane-parallel, anisotropic scattering medium with specular-diffuse boundaries

Using the intensity with high directional resolution obtained by the Basic-DRESOR method as an initial guess, which is substituted into the integrated radiative transfer equation (IRTE), an iterative algorithm is proposed, called the Iterative-DRESOR method. This method can reduce the error levels of the intensity from several percent using the Basic-DRESOR method to a level of less than 1.0×10⁻⁶ with acceptable computation costs. The method is also validated against the exact heat flux in literature in some cases. It further clarifies some uncertain results for the reflectance in a pure, linearly anisotropic scattering medium with specular-diffuse boundaries. The directional distributions of intensity are obviously influenced by the reflecting modes of the boundary, especially in the zone near the boundary. The reflecting mode of an emitting boundary has little effect on the transmittance or reflectance. The reflecting mode of a non-emitting boundary also has little effect on the transmittance, but it obviously influences the reflectance. The difference between the reflectance for specular and diffuse boundaries increases at first, and then decreases, as the optical thickness of the medium increases. The difference will decrease as the scattering albedo of the medium increases, and it is negligible when the medium is pure scattering. The effect of the scattering phase function of the medium on the difference can also not be ignored. The Iterative-DRESOR method is expected to strengthen the capability of the Monte Carlo method to produce accurate results and to validate the results of other methods to solve RTE.     

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.     

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.     

Comparison between radiative transfer theory and the simplified spherical harmonics approximation for a semi-infinite geometry

In this study, the third-order simplified spherical harmonics equations (SP3), an approximation of the radiative transfer equation, are solved for a semi-infinite geometry considering the exact simplified spherical harmonics boundary conditions. The obtained Green's function is compared to radiative transfer calculations and the diffusion theory. In general, it is shown that the SP3 equations provide better results than the diffusion approximation in media with high absorption coefficient values but no improvement is found for small distances to the source.     

Radiation transport and internal reflection in a sphere

We construct an integral equation for the flux intensity in a scattering and absorbing medium using the integro-differential form of the radiative transfer equation in a sphere. The sphere is uniformly irradiated by an external source of arbitrary angular distribution. The Fresnel boundary conditions, which incorporate reflection and refraction, are used. For the special cases of a non-scattering medium, and in the limit of an optically transparent medium, we obtain exact solutions for specular and diffuse refection. Some numerical examples are given which give qualitative agreement with some recent work of Tian and Chiu (JQSRT, 2005).     

Attenuation of solar radiation by a water mist from the ultraviolet to the infrared range

A model is developed for the hemispherical transmittance of direct and scattered solar radiation from a cloudless atmosphere by a mist layer of water droplets in order to investigate the potential of water misting systems to serve as a protection from solar irradiation with particular emphasis on harmful UV radiation. The proposed model is based on published spectral experimental data for solar irradiation, Mie theory for interaction of the radiation with single spherical droplets, and radiative transfer theory. Known limiting solutions are employed to simplify the Mie calculations. The modified two-flux approximation is used to account for both direct and diffuse irradiation in lieu of a numerical solution for the full radiative transfer equation in anisotropically scattering media. The role of the governing parameters of a disperse water curtain of water droplets, water content, and droplet size for sample conditions is studied in some detail, particularly in the near-ultraviolet part of the spectrum where radiation can result in human tissue damage.     

A theoretical study of radiation between two concentric spheres using a modified discrete ordinates method associated with Legendre transform

A modified discrete ordinates method (DOM) is used in spherical participating media. The radiative intensity is broken up into two components. One component is traced back to the enclosure's source. It is called direct intensity. The other component is rather traced back to the contribution of the medium itself. It is called diffuse intensity. Thus, the radiative transfer equation (RTE) is transformed into two simultaneous equations: a direct RTE and a diffuse RTE. The direct RTE is solved analytically. The diffuse RTE is solved numerically using the DOM. The streaming angular derivative term appearing in spherical geometry is modeled by making use of the Finite Legendre Transform. We study a pure radiation transfer problem between two concentric spheres. The medium is assumed to be gray and isotropically scattering. The limiting spheres are considered to be opaque, gray, diffusely emitting and diffusely reflecting with uniform emissivity over each surface. The obtained results are compared with available cases reported in the literature. In particular, relative importance of the direct radiation in optically thin media is studied.     

Error analysis of the adding/doubling method for radiative scattering in inhomogeneous media

The doubling method is a fast and exact procedure for calculating radiative transfer in a homogeneous, scattering, plane-parallel medium. It can also be used in the adding mode for an inhomogeneous medium that is approximated by a finite number of homogeneous sublayers with different radiative properties. The errors caused by this approximation are analyzed in this paper through comparison with invariant imbedding calculations that are slow but exact for inhomogeneous media. A procedure is developed so errors can be estimated and controlled when using the faster adding/doubling calculations.     

Radiation transport and internal reflection in a two region, turbid sphere

We construct an integral equation for the flux intensity in a scattering and absorbing two-region turbid spherical medium using the integro-differential form of the radiative transfer equation. The sphere is uniformly irradiated by an external source of arbitrary angular distribution and contains a distributed volume source. Anisotropic scattering is accounted for by the transport approximation. The Fresnel boundary conditions, which incorporate reflection and refraction, are used at the outer surface and at the interface between the two regions. In this respect, some new interfacial boundary conditions are introduced. For the special case of a non-scattering medium, we obtain exact solutions for specular reflection. Some numerical examples are given which show qualitative agreement with some recent work of other authors. Of particular interest are the emergent angular distribution and the albedo of the surface as a function of the refractive index and the radii of the two regions. We also draw attention to the fact that the boundary conditions at the interface differ according to the relative values of the refractive indices in the two regions. The interfacial boundary conditions for use in diffusion theory are derived and compared with those of Aronson [Boundary conditions for diffusion of light. J opt Soc Am 1995;12:2532]. In appendix B, we show how diffusion theory may be used to include scattering into the problem in a simple way.     

Finite element approximation of the Fokker-Planck equation for diffuse optical tomography

In diffuse optical tomography, light transport theory is used to describe photon propagation inside turbid medium. A commonly used simplification for the radiative transport equation is the diffusion approximation due to computational feasibility. However, it is known that the diffusion approximation is not valid close to the sources and boundary and in low-scattering regions. Fokker-Planck equation describes light propagation when scattering is forward-peaked. In this article a numerical solution of the Fokker-Planck equation using finite element method is developed. Approach is validated against Monte Carlo simulation and compared with the diffusion approximation. The results show that the Fokker-Planck equation gives equal or better results than the diffusion approximation on the boundary of a homogeneous medium and in turbid medium containing a low-scattering region when scattering is forward-peaked.     

Light scattering in disperse layers with a high concentration of optically soft particles

A model describing the propagation of radiation in disperse close-packed layers of optically soft particles is proposed. The model is based on the radiative transfer equation used for describing a light field in the multiple scattering regime and on the interference approximation for describing collective scattering effects; it also takes into account the Fresnel reflection at the layer boundaries. To describe the single-scattering characteristics, the Mie formulas are used. The model allows one to consider scattering in layers of any optical thickness. Calculations of the angular structure of the radiation leaving the layer and estimates of the errors of the calculations were performed. The possibility of forming a maximum at small scattering angles in media with correlated particles is shown. The theoretical results are in good agreement with the available experimental data.     

Evaluation method for radiative heat transfer in polydisperse water droplets

Simplifications of the model for nongray radiative heat transfer analysis in participating media comprised of polydisperse water droplets are presented. Databases of the radiative properties for a water droplet over a wide range of wavelengths and diameters are constructed using rigorous Mie theory. The accuracy of the radiative properties obtained from the database interpolation is validated by comparing them with those obtained from the Mie calculations. The radiative properties of polydisperse water droplets are compared with those of monodisperse water droplets with equivalent mean diameters. Nongray radiative heat transfer in the anisotropic scattering fog layer, including direct and diffuse solar irradiations and infrared sky flux, is analyzed using REM². The radiative heat fluxes within the fog layer containing polydisperse water droplets are compared with those in the layer containing monodisperse water droplets. Through numerical simulation of the radiative heat transfer, polydisperse water droplets can be approximated by using the Sauter diameter, a technique that can be useful in several research fields, such as engineering and atmospheric science. Although this approximation is valid in the case of pure radiative transfer problems, the Sauter diameter is reconfirmed to be the appropriate diameter for approximating problems in radiative heat transfer, although volume-length mean diameter shows better accordance in some cases. The CPU time for nongray radiative heat transfer analysis with a fog model is evaluated. It is proved that the CPU time is decreased by using the databases and the approximation method for polydisperse particulate media.     

漫反射各向异性散射介质内的温度热流场

本文采用射线踪迹结合节点分析法和谱带模型,研究了漫反射不透明边界下吸收、发射、各向异性散射介质内的热辐射传递过程。考虑介质辐射能的入射和散射方向,导出漫反射、不透明边界、各向异性散射介质的辐射传递系数。在辐射平衡的情况下,考察了表面发射率和散射反照率对介质内辐射热流和温度场的影响。研究表明,介质不透明边界处存在温度跃迁现象,而且,内界面发射率越大,相应界面温度跃迁越小。