Radiative transfer in absorbing, emitting and linearly anisotropic-scattering inhomogeneous cylindrical medium

Radiative integral transfer equations for one-dimensional solid cylinder with absorbing, emitting and linearly anisotropic-scattering inhomogeneous medium were derived by Abulwafa et al. (JQSRT 62 (1999) 755). The anisotropic terms in the integral equations and their results for anisotropic benchmark problems (JQSRT 66 (2000) 487) are inaccurate. In this study, the integral equations for absorbing, emitting and linearly anisotropic-scattering medium are rederived, and the integral equations for one-dimensional solid cylindrical medium are solved. The results are compared with selected cases using the discrete ordinates S₁₆ and the exact solutions available in the literature.     

Radiative transfer in cylindrical media

A numerical method for the solution of the radiative transfer equation in a circularly symmetric, cylindrical region is developed. The transfer equation is formulated as a second-order differential equation resulting in a set of tridiagonal difference equations. This form is particularly well suited to line formation and energy balance calculations using the complete linearization method. Several numerical examples are presented.     

Radiative transfer in finite participating atmospheric aerosol media

The properties of radiation transfer through a plane-parallel atmospheric aerosol medium has been studied. It has been done by employing Mie theory to calculate the radiation transfer scattering parameters of the medium in the form of extinction, scattering, and absorption efficiencies. Then, the equation of radiative transfer through a plane-parallel atmosphere of aerosol has been solved for partial heat fluxes using two different analytical techniques, namely, the Variational Pomraning -Eddington approximation and Galerkin technique. Average efficiencies over log-normal and modified gamma size distributions are calculated. Therefore, the radiative properties of Carbon, Anthracite, Bituminous, Lignite, and Fly ash have been calculated. The obtained numerical results show very good agreement with each other in addition to the previous published work.     

Radiative transfer in stochastically inhomogeneous media

New radiative transfer theory is developed for stochastically inhomogeneous scattering media. The three-dimensional shapes and large scale (compared to the mean free path) structures of the media are modeled by stochastic interfaces separating regions of different scattering properties. The small scale fluctuations are characterized by a pair-correlation function. The radiative transfer equation is extended to include individual scattering and propagation probabilities of a ray for each subregion as well as the probability for a ray to cross the interface between two subregions. The propagation probability is found to depend on the entire preceding path of the ray; the present formulation accounts for the two previous scatterings. A new adding/doubling algorithm is developed to solve this problem numerically. Transmission through a cloud layer and backward scattering seem to be particularly sensitive to inhomogeneities.     

A general integration method for radiative transfer in 3-D non-homogeneous cylindrical media with anisotropic scattering

A general set of integral equations is presented to solve 3-D radiative heat transfer problems in emitting, absorbing and linear anisotropic scattering finite hollow or solid cylinders with non-homogeneous media. By tracing a ray to compute the intensity,it is much easier to handle the spatial change properties including extinction coefficient. Both the continuous change property and step-change property are dealt with without difficulties. The solid angle integration in getting the incident radiation and heat fluxes is represented by the bounding surface integration. In order to avoid the singularity problem near the bounding surface, the surface integrations are transformed to new modified integral equations by mathematical methods. By doing so, we get more flexible general integral equations applicable to all cases (3-D solid cylinders, 3-D hollow cylinders, finite cylinders or infinite cylinders). This scheme has been verified by comparing the results with published data in the literature. It is believed that this method will be useful in combined radiation and convection heat transfer problems.     

Radiative transfer in finite inhomogeneous plane-parallel atmospheres

The F_N method is used to deduce accurate numerical results for the exit distributions of radiation relevant to a finite, plane-parallel atmosphere with an exponentially varying albedo for single scattering.     

Radiative transfer in finite volcanic eruption clouds

Radiation transfer through a volcanic aerosol medium has been studied. The radiation transfer properties of the medium as scattering, absorption and extinction coefficients are calculated using the Mie scattering theory. Average coefficients over the size parameter and the radiation wavelength are calculated. The radiation heat fluxes for volcanic eruption ash medium are calculated using the Variational Pomraning–Eddington approximation and compared with those obtained from the Galerkin method. The comparison showed very good agreement.     

Reduction of radiative transfer problems in semi-infinite media to linear Fredholm integral equations

Linear Fredholm integral equations are derived for the Stokes vector of polarized radiation, emergent from a scattering plane parallel semi-infinite medium, by means of the full range orthogonality and completeness properties of Case's eigensolutions. A renormalization concerning the eigenmode with the greatest discrete eigenvalue is applied, which permits us to obtain a new integral equation for the zeroth Fourier component of the radiation field. The kernel of the integral equations is given in terms of Case's eigenfunctions or of the Green's function matrix for an infinite medium. For isotropic scattering, it is shown that the integral equation can be solved by means of a very rapidly convergent Neumann series. Physical arguments lead to the conclusion that the renormalized Fredholm integral equations are well suited also for arbitrary phase matrices.     

Radiative transfer in plane media without azimuthal symmetry

A rigorous integral theory is presented in this paper for the solution of radiative heat transfer problems in stratified media, when dependence on the azimuth of the propagating radiation must be taken into account. Anisotropy of scattering and specular and diffuse reflection from the bounding walls are incorporated in the final system of linear integral equations of Fredholm's type. A simple case of physical interest is considered in more detail, and solved explicitly by a constructive technique. Numerical results are reported and briefly discussed.     

On the finite volume method and the discrete ordinates method regarding radiative heat transfer in acute forward anisotropic scattering media

The ability of the finite volume method (FVM) and the discrete ordinates method (DOM) to model radiative heat transfer in acute forward anisotropic scattering media has been investigated. The test case involves a purely scattering medium in a cubic enclosure, irradiated by one boundary with diffuse emission. Four phase functions have been considered: three of the Henyey-Greenstein type with respective asymmetry factors of 0.2, 0.8 and 0.93, and a Mie phase function with a strong forward scattering peak (computed for a size parameter of 245 and corresponding to an asymmetry factor of 0.93). Results obtained with the FVM are in good agreement with Monte Carlo reference solutions, whatever the level of acute anisotropic scattering (for asymmetry factors up to 0.93). The DOM combined with the renormalization procedures of the phase function proposed by Kim and Lee (Effect of anisotropic scattering on radiative heat transfer in two-dimensional rectangular enclosures. Int J Heat Mass Transfer 1988;31:1711-21. [1]) and Wiscombe (On initialization error and flux conservation in the doubling method. JQSRT 1976;18:637-58. [2]) provides accurate results only for the smallest asymmetry factor. As the asymmetry factor increases, the renormalization procedures induce strong modifications in the values of the discretized phase function resulting in an underestimation of the effective attenuation by scattering. This error has been found to increase with optical thickness. In fact, when using the DOM, results would be more accurate combining this method with a Delta-Eddington approximation of the phase function, instead of using the actual phase function which is altered too much by renormalization.     

Wavelet method for solving integral equations of simple liquids

A new algorithm is developed to solve integral equations for simple liquids. The algorithm is based on the discrete wavelet transform of radial distribution functions. The Coifman 2 basis set is employed for the wavelet treatment. Using the algorithm, we have calculated structural and thermodynamic properties of a Lennard–Jones fluid in a wide range of energy and size parameters of the fluid.     

Depth-controlled reconstruction of 3D integral image using synthesized intermediate sub-images

In this paper, we propose a method that controls the depth of the three-dimensional (3D) object existing over the depth-of-focus in integral imaging. The depth control method is performed only in a computer by synthesizing the intermediate sub-images between original sub-images obtained by transforming the captured elemental images. In the reconstruction process, we can obtain reconstructed 3D images with the better image quality within depth-of-focus than that reconstructed over the depth-of-focus. To demonstrate the feasibility of our method, optical and computational experiments are carried out and its results are presented.     

Optical display of true 3D objects in depth-priority integral imaging using an active sensor

In this paper, we propose a system combining the pickup process using an active sensor and the display process using depth-priority integral imaging (DPII) system to display true three-dimensional (3D) objects within large depth through real and virtual image fields. The active sensor provides depth map and color images of 3D objects. Using captured depth map and original color images, elemental images are computationally synthesized and displayed optically in DPII system. Proposed system provides scaling of 3D scenes for true 3D object. To show the usefulness of proposed system, we carry out the experiment for true 3D objects of three character patterns and present the experimental results.     

Electromagnetic integral equations requiring small numbers of Krylov-subspace iterations

We present a new class of integral equations for the solution of problems of scattering of electromagnetic fields by perfectly conducting bodies. Like the classical Combined Field Integral Equation (CFIE), our formulation results from a representation of the scattered field as a combination of magnetic- and electric-dipole distributions on the surface of the scatterer. In contrast with the classical equations, however, the electric-dipole operator we use contains a regularizing operator; we call the resulting equations Regularized Combined Field Integral Equations (CFIE-R). Unlike the CFIE, the CFIE-R are Fredholm equations which, we show, are uniquely solvable; our selection of coupling parameters, further, yields CFIE-R operators with excellent spectral distributions—with closely clustered eigenvalues—so that small numbers of iterations suffice to solve the corresponding equations by means of Krylov subspace iterative solvers such as GMRES. The regularizing operators are constructed on the basis of the single layer operator, and can thus be incorporated easily within any existing surface integral equation implementation for the solution of the classical CFIE. We present one such methodology: a high-order Nyström approach based on use of partitions of unity and trapezoidal-rule integration. A variety of numerical results demonstrate very significant gains in computational costs that can result from the new formulations, for a given accuracy, over those arising from previous approaches.     

Three-dimensional image correlator using fast computational integral imaging reconstruction method based on pixel-to-pixel mapping

In this paper, a three-dimensional (3D) image correlator using a fast computational integral imaging reconstruction (CIIR) method based on a pixel-to-pixel mapping is proposed. In order to implement the fast CIIR method, we replace the magnification process in the conventional CIIR by a pixel-to-pixel mapping. The proposed fast CIIR method reconstructs two sorts of plane images; a plane image whose quality is sufficient, and a dot pattern plane image insufficient to view. This property is very useful to enhance the performance of a CIIR-based image correlator. Thus, we apply the fast CIIR method to a CIIR-based image correlator. To show the feasibility of the proposed method, some preliminary experiments on both pattern correlation and computational cost are carried out, and the results are presented. Our experimental results indicate that the proposed image correlator is superior to the previous method in terms of both correlation performance and complexity.     

On “Radiative transfer in three-dimensional rectangular enclosures containing inhomogeneous, anisotropically scattering media”

Our 1985 paper (JQSRT 1985; 33: 533-549) reported the result of the research we conducted back then to better understand heat transfer processes in large-scale combustion chambers, especially in pulverized coal-fired furnaces. It was one of the first works exploring radiative transfer in three-dimensional enclosures where absorption and scattering coefficients due to combustion particles and gases were allowed to vary within the medium. This flexibility of the mathematical model made it useful for applications to realistic furnaces and different types of high-temperature systems. This note briefly discusses the motivation behind the paper and the immediate extension of the idea to different systems.     

Radiative transfer in spatially heterogeneous,two-dimensional,anisotropically scattering media

A method is presented for solving the radiative transfer equation for a general anisotropically scattering and emitting medium exposed to arbitrary boundary radiation conditions. The method allows, in principle, for quite arbitrary spatial variability in the scattering and extinction properties of the medium. We formulate the method in the context of 2-dimensional radiative transfer and describe general solution procedures, based on the principles of invariant imbedding, which are applied in the form of doubling algorithms to obtain solutions for optically thick media. Some selected results are shown to demonstrate the versatility of the approach.     

Radiative transfer equation for graded index medium in cylindrical and spherical coordinate systems

In graded index medium, the ray goes along a curved path determined by Fermat principle, and the curved ray-tracing is very difficult and complex. To avoid the complicated and time-consuming computation of curved ray trajectory, the methods not based on ray-tracing technique need to be developed for the solution of radiative transfer in graded index medium. For this purpose, in this paper the streaming operator along a curved ray trajectory in original radiative transfer equation for graded index medium is transformed and expressed in spatial and angular ordinates and the radiative transfer equation for graded index medium in cylindrical and spherical coordinate systems are derived. The conservative and the non-conservative forms of radiative transfer equation for three-dimensional graded index medium are given, which can be used as base equations to develop the numerical simulation methods, such as finite volume method, discrete ordinates method, and finite element method, for radiative transfer in graded index medium in cylindrical and spherical coordinate systems.