Meshless approach for coupled radiative and conductive heat transfer in one-dimensional graded index medium

A meshless local Petrov-Galerkin (MLPG) approach is employed for solving the coupled radiative and conductive heat transfer in a one-dimensional slab with graded index media. The angular distribution term in discrete ordinate equation of radiative transfer within a one-dimensional graded index slab is discretized by a step scheme, and the meshless approach for radiative transfer is based on the discrete ordinate equation. A moving least-squares approximation is used to construct the shape function. Two particular test cases for coupled radiative and conductive heat transfer within a one-dimensional graded index slab are examined to verify this new approximate method. The temperatures and the radiative heat fluxes are obtained. The results are compared with the other benchmark approximate solutions. By comparison, the results show that the MLPG approach has a good accuracy in solving the coupled radiative and conductive heat transfer in one-dimensional graded index media.     

A radiative transfer model for an idealized and non-scattering atmosphere and its application for ground-based remote sensing

Inversion of tropospheric profiles from ground-based microwave measurements requires a simple and accurate model for calculating the brightness temperatures as received by the radiometer. In the first part, an analytic solution of the radiative transfer equation is derived for an exponentially decaying absorption coefficient and a linear temperature gradient. Based on the obtained analytic expressions, a discretized radiative transfer scheme is developed in the second part. The new scheme incorporates the generic behavior of the atmosphere with the effect that brightness temperatures can be modeled more accurately and with fewer grid points compared to commonly used radiative transfer schemes. The brightness temperature modeling accuracy was improved by a factor of six. The results suggest that the model could be employed for the retrieval of temperature and humidity profiles.     

A discrete bidirectional reflectance model in remote sensing

A numerical scheme is devised to develop a discrete bidirectional reflectance model. A pseudospectral method is utilized with which the discrete solution of the radiative transfer equation is made part of this development. To produce discrete bidirectional reflectance values, a set of algebraic equations is solved. Illustrative examples are given to demonstrate the performance of the developed model.     

A numerical scheme for the solution of axisymmetric radiative transfer problems in irregular domains filled by media with opaque and transparent diffuse and specular boundaries

This paper presents a new numerical scheme of the discrete ordinates method for the solution of axisymmetric radiative transfer problems in irregular domains filled by media with opaque and transparent diffuse and specular (Fresnel) boundaries and interfaces. New test problems of radiative transfer, which describe radiative transfer in domains with Fresnel interfaces, are proposed in this paper. These problems admit analytic solutions and can be used as benchmark ones. The proposed scheme is applied to the solution of the problems. Numerical results show that the presence of Fresnel interfaces leads to an appreciably larger error in numerical solution. This is connected with the “discontinuity” of the Fresnel reflectivity, which, through numerical diffusion, leads to the distortion of numerical solution. Modification of the scheme allows to reduce the numerical error.     

空间差分格式对有限体积法数值散射的影响

数值散射是辐射传递方程近似算法中最常见的离散误差。本文主要讨论空间差分格式对有限体积法数值散射的影响。构造激光平行及倾斜入射的物理模型,验证和比较阶梯格式、中心差分格式及指数格式下温度场的计算精度及数值散射特性。计算结果表明,在激光平行入射与倾斜入射两种情况下,阶梯格式引起的的数值散射比菱形格式及指数格式要多,但其计算精度高于菱形及指数格式。不同激光入射条件下,各种差分格式表现出的数值散射分布有明显的差异。     

Spectral collocation method with a flexible angular discretization scheme for radiative transfer in multi-layer graded index medium

The spectral collocation method (SCM) is employed to solve the radiative transfer in multi-layer semitransparent medium with graded index. A new flexible angular discretization scheme is employed to discretize the solid angle domain freely to overcome the limit of the number of discrete radiative direction when adopting traditional S_N discrete ordinate scheme. Three radial basis function interpolation approaches, named as multi-quadric (MQ), inverse multi-quadric (IMQ) and inverse quadratic (IQ) interpolation, are employed to couple the radiative intensity at the interface between two adjacent layers and numerical experiments show that MQ interpolation has the highest accuracy and best stability. Variable radiative transfer problems in double-layer semitransparent media with different thermophysical properties are investigated and the influence of these thermophysical properties on the radiative transfer procedure in double-layer semitransparent media is also analyzed. All the simulated results show that the present SCM with the new angular discretization scheme can predict the radiative transfer in multi-layer semitransparent medium with graded index efficiently and accurately.     

Radiative transfer theory with time delay for the effect of pulse entrapping in a resonant random medium: General transfer equation and point-like scatterer model

Abstract

A pulse propagation of a vector electromagnetic wave field in a discrete random medium under the condition of Mie resonant scattering is considered on the basis of the Bethe–Salpeter equation in the two-frequency domain in the form of an exact kinetic equation which takes into account the energy accumulation inside scatterers. The kinetic equation is simplified using the transverse field and far wave zone approximations which give a new general tensor radiative transfer equation with strong time delay by resonant scattering. This new general radiative transfer equation, being specified in terms of the low-density limit and the resonant point-like scatterer model, takes the form of a new tensor radiative transfer equation with three Lorentzian time-delay kernels by resonant scattering. In contrast to the known phenomenological scalar Sobolev equation with one Lorentzian time-delay kernel, the derived radiative transfer equation does take into account effects of (i) the radiation polarization, (ii) the energy accumulation inside scatterers, (iii) the time delay in three terms, namely in terms with the Rayleigh phase tensor, the extinction coefficient and a coefficient of the energy accumulation inside scatterers, respectively (i.e. not only in a term with the Rayleigh phase tensor). It is worth noting that the derived radiative transfer equation is coordinated with Poynting's theorem for non-stationary radiation, unlike the Sobolev equation. The derived radiative transfer equation is applied to study the Compton–Milne effect of a pulse entrapping by its diffuse reflection from the semi-infinite random medium when the pulse, while propagating in the medium, spends most of its time inside scatterers. This specific albedo problem for the derived radiative transfer equation is resolved in scalar approximation using a version of the time-dependent invariance principle. In fact, the scattering function of the diffusely reflected pulse is expressed in terms of a generalized time-dependent Chandrasekhar H-function which satisfies a governing nonlinear integral equation. Simple analytic asymptotics are obtained for the scattering function of the front and the back parts of the diffusely reflected Dirac delta function incident pulse, depending on time, the angle of reflection, the mean free time, the microscopic time delay and a parameter of the energy accumulation inside scatterers. These asymptotics show quantitatively how the rate of increase of the front part and the rate of decrease of the rear part of the diffusely reflected pulse become slower with transition from the regime of conventional radiative transfer to that of pulse entrapping in the resonant random medium.     

Matrix formulations of radiative transfer including the polarization effect in a coupled atmosphere-ocean system

A vector radiative transfer model has been developed for a coupled atmosphere-ocean system. The radiative transfer scheme is based on the discrete ordinate and matrix operator methods. The reflection/transmission matrices and source vectors are obtained for each atmospheric or oceanic layer through the discrete ordinate solution. The vertically inhomogeneous system is constructed using the matrix operator method, which combines the radiative interaction between the layers. This radiative transfer scheme is flexible for a vertically inhomogeneous system including the oceanic layers as well as the ocean surface. Compared with the benchmark results, the computational error attributable to the radiative transfer scheme has been less than 0.1% in the case of eight discrete ordinate directions. Furthermore, increasing the number of discrete ordinate directions has produced computations with higher accuracy. Based on our radiative transfer scheme, simulations of sun glint radiation have been presented for wavelengths of 670 nm and 1.6 μm. Results of simulations have shown reasonable characteristics of the sun glint radiation such as the strongly peaked, but slightly smoothed radiation by the rough ocean surface and depolarization through multiple scattering by the aerosol-loaded atmosphere. The radiative transfer scheme of this paper has been implemented to the numerical model named Pstar as one of the OpenCLASTR/STAR radiative transfer code systems, which are widely applied to many radiative transfer problems, including the polarization effect.     

Verification of numerical simulation method for entropy generation of radiation heat transfer in semitransparent medium

The numerical simulation method of radiative entropy generation in participating media presented by Caldas and Semiao [Entropy generation through radiative transfer in participating media: analysis and numerical computation. JQSRT 2005;96:423-37] is extended to analyze the radiative entropy generation in the enclosures filled with semitransparent media. A discrete ordinates method is used to solve radiative transfer equation and radiative entropy generation. Two different examples are employed to verify the numerical simulation method of radiative entropy generation in the enclosure. Numerical results of dimensionless radiative entropy generation of enclosure are identical to that of entire thermodynamics analysis for the enclosure system. This numerical simulation method can be used in the entropy generation analysis of high-temperature systems such as boilers and furnaces, in which radiation is the dominant mode of heat transfer.     

A direct method for the integration of the equation of radiative transfer in a turbid atmosphere

In this paper we apply a numerical method on the equation of radiative transfer in a turbid atmosphere. The solution is obtained by means of direct integration of the equation of radiative transfer without any circuitous series development.     

Adjoint problem of radiative transfer for a pseudo-spherical atmosphere and general viewing geometries

We formulate the adjoint radiative transfer for a pseudo-spherical atmosphere and various retrieval scenarios. The single scattering radiance is computed in a spherical atmosphere by using the source integration technique, while for the multiple scattering radiance we formulate an one-dimensional adjoint radiative transfer equation in a plane-parallel atmosphere. The adjoint solution of the radiative transfer equation is obtained by employing the discrete ordinate method with matrix exponential. We provide an abbreviated derivation of our formalism as well as a discussion of the numerical implementation of the theory.     

Half-Moment Closure for Radiative Transfer Equations

In this paper a moment method for radiative transfer equations is considered which has been developed and investigated using different approaches. Problems appearing for this moment system for boundary value problems using Maxwell-type boundary conditions are described. A new method based on the consideration of positive and negative half fluxes is developed and shown to overcome the above problems. Moreover, a numerical scheme and numerical results for the new moment system are presented.     

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.     

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.     

Numerical developments for short-pulsed Near Infra-Red laser spectroscopy. Part II: inverse treatment

Indirect optical spectroscopy or tomography, that is, mapping of optical properties in scattering and absorption inside a medium given a set of measurements at the boundaries, is highly dependent on the radiative transfer model used to track radiative energy propagation in semi-transparent materials. In the first part of this study, a numerical tool adapted for treating radiative transfer in the frame of short-pulsed laser beam interaction with non-homogeneous matter has been presented. In this paper, it is intended to show how such numerical tools can undergo inversion through adjoint treatment or reverse differentiation.Adjoint models, as well as reverse differentiation, are used in order to allow an efficient computation of the gradient, in the unknown optical parameters space, of an objective or cost function estimating the residual between data obtained at the boundary and predictions by numerical simulations. This gradient is a crucial indication as to update, through line minimization, the set of internal optical properties of the medium.First, the theoretical background of the inverse treatments, both reverse differentiation and adjoint model, for the transient radiative transfer equation model introduced in Part I is developed. Second, different reconstruction configurations are presented. Time-dependent sampling and time filtering effects of the measurements are addressed. Image reconstructions from simulated data are achieved for material phantoms of simple geometry.     

非规则形状介质内辐射-导热耦合传热的间断有限元求解

采用间断有限元法(discontinuous finite element method,DFEM)求解非规则形状介质内的辐射导热耦合传热问题,得到了典型非规则形状介质内辐射导热耦合传热问题的高精度数值结果.和传统连续型有限元方法不同,DFEM将计算区域划分成相互独立的离散单元,形函数的构造、未知量的加权近似以及控制方程的求解均在每一个离散单元上进行.通过在单元之间施加迎风格式的数值通量,DFEM保证了整个计算区域的连续性,因此这种方法兼具良好的几何灵活性和局部守恒性.推导了辐射传输方程和能量扩散方程的射导热耦合传热问题,得到了典型非规则形状介质内辐射导热耦合传热的高精度数值结果.     

Least-squares collocation meshless approach for radiative heat transfer in absorbing and scattering media

A least-squares collocation meshless method is employed for solving the radiative heat transfer in absorbing, emitting and scattering media. The least-squares collocation meshless method for radiative transfer is based on the discrete ordinates equation. A moving least-squares approximation is applied to construct the trial functions. Except for the collocation points which are used to construct the trial functions, a number of auxiliary points are also adopted to form the total residuals of the problem. The least-squares technique is used to obtain the solution of the problem by minimizing the summation of residuals of all collocation and auxiliary points. Three numerical examples are studied to illustrate the performance of this new solution method. The numerical results are compared with the other benchmark approximate solutions. By comparison, the results show that the least-squares collocation meshless method is efficient, accurate and stable, and can be used for solving the radiative heat transfer in absorbing, emitting and scattering media.     

圆筒形半透明介质内非稳态复合导热与辐射的研究

本文研究细长电加热体在圆筒形半透明介质中的温度响应，通过对控制非稳态导热与辐射复合传热过程的积分-微分方程直接进行数值求解，分析了热辐射对内部径向热流及温度变化的影响．模拟计算结果对热线法测量半透明介质导热系数的研究具有理论指导意义。