An Iterative Method for Solving of Coupled Equations for Conductive-Radiative Heat Transfer in Dielectric Layers

Themathematical model for describing combined conductive-radiative heat transfer in a dielectric layer, which emits, absorbs, and scatters IR radiation both in its volume and on the boundary, has been considered. A nonlinear stationary boundary-value problem for coupled heat and radiation transfer equations for the layer, which exchanges by energy with external medium by convection and radiation, has been formulated. In the case of optically thick layer, when its thickness is much more of photon-free path, the problem becomes a singularly perturbed one. In the inverse case of optically thin layer, the problem is regularly perturbed, and it becomes a regular (unperturbed) one, when the layer’s thickness is of order of several photon-free paths. An iterative method for solving of the unperturbed problem has been developed and its convergence has been tested numerically. With the use of the method, the temperature field and radiation fluxes have been studied. The model and method can be used for development of noncontact methods for temperature testing in dielectrics and for nondestructive determination of its radiation properties on the base of the data obtained by remote measuring of IR radiation emitted by the layer.


Introduction
Necessity for studying of coupled conductive and radiative heat exchange in continuous mediums arises in various theoretical and applied disciplines: in thermal physics, in plasma physics, in solid state physics, in space technologies, in technologies of semiconductor, glass, and polymer materials, and in methodologies for thermal measurement and nondestructive testing.Various methods are used for analytical and/or numerical solving of boundary-value problems, formulated with the use of certain mathematical models describing the coupled processes of conductive and radiative heat transfer.
In the case of opaque bodies (metals), possessing high thermal conductivity, one can neglect radiative heat transfer in the body's volume and take into account only the radiative heat flux emitting by body's boundary outward.Outer IR radiation, which is incident on the boundary of such body, is totally absorbed by the boundary.In agreement with the Stefan-Boltzmann law the radiative flux is proportional to the fourth degree of surface temperature [1].Hence we obtain in this case a nonlinear boundary-value problem for linear heat conduction equation.
Dielectric objects (nonconductors), such as glass, plastics, heat insulations, semiconductors, and biological tissues, are translucent for IR radiation.Radiation emitted by inner points of such object can spread on considerable distances in the body's volume and can stream away from the volume through the boundary.These two radiation energy fluxes can be comparable with volumetric conductive flux and surface convective flux correspondingly.Hence they should be accounted in the total energy balance.
To describe energy transfer in such mediums, the "exact" macroscopic model can be used.It includes coupled integrodifferential equations, heat conduction equation (HCE) and radiation transfer equation (RTE) [1].They involve two unknown functions defined in the body's volume V: temperature (r, ), dependent on spatial coordinates r ∈ V and time , and spectral radiation intensity  ] (r, , ) being a function of IR radiation frequency ], spatial coordinates of radiation intensity as [1]   = ∫ ∞ 0 ∫ 4  ] (r, , ) ].Hence, to determine q  one should anyway solve RTE, which, in turn, is depending on the temperature field.
With the use of an appropriate iterative process HCE and RTE can be solved consecutively.
Examples of solutions of such coupled problems for a plane layer with the use of numerical methods can be found in [6][7][8][9].Parabolic HCE was used in [6,9], in which temperature field in semitransparent medium subjected to a pulse irradiation was studied numerically.In [9], in particular, implicit central-difference scheme was employed for solving HCE, while a discrete curved ray-tracing method was used to solve RTE.In [7,8], the hyperbolic HCE was used in studying of non-Fourier effects on transient coupled radiativeconductive heat transfer in one-dimensional semitransparent medium subjected to a periodic irradiation.Here hyperbolic HCE was solved by the flux-splitting method, and RTE was solved by the discrete ordinate method.
We would notice here that the results of numerical study presented, in particular, in [9], were obtained for the layer of optical thickness ℎ = 0.1, that is, for optically thin medium.In [8], where the influence of ℎ on temperature field in the layer was studied, the maximum used value of optical thickness was taken to be equal to 1.0.
As we will show further, the problem in coupled "exact" formulation can be solved with the use of conventional method only for layers with optical thickness ℎ of order unity or smaller.The reason is the presence of two spatial scales in the problem: the layer's physical thickness ℎ and the photon mean free path  0 .As a consequence, the coefficient at the highest spatial derivative in HCE decreases with increasing of ℎ ≡ ℎ/ 0 as 1/ℎ 2 .This becomes apparent in occurrence of so-called boundary layer of thickness of order  0 at sufficiently large ℎ.The temperature and radiation intensity are sharply varied in the boundary layer.Whereupon the problem becomes singularly disturbed at large ℎ, special approaches should be used to analyze it numerically [10].
But some problems for optically thick layers can be solved numerically with the use of an approach, which known as Rosseland's or diffusion approximation [1,3].In this approximation, radiation flux q  is also taken to be proportional to temperature gradient q  = −  ∇, where   is radiation conductivity, which is defined as   = 16/3 SB  3 / [3], and  SB stands for the Stefan-Boltzmann constant  SB ≈ 5.670367 × 10 −8 W⋅m 2 ⋅K −4 .
Hence, applying the diffusion approximation, we exclude RTE from the problems of conductive-radiative heat exchange.Now we can deal with HCE for the medium with the modified heat conductivity defined as  =   +   , which depends on the temperature.As a consequence, the second spatial scale  0 also vanishes, and the problem remains regular at arbitrary optical thickness ℎ.
This approach was successfully used in various practical problems, in particular, for evaluation of effective heat resistance properties of low-conductivity coats [11][12][13].But it is important to realize that the diffusive approximation is valid only at sufficiently low temperature gradients and only in Advances in Mathematical Physics 3 inner medium's domains, remote from the boundary V on distance exceeding the mean free path  0 .In a boundary layer with the thickness of order  0 this approach does not work.
According to this approximation the IR radiation streaming out the object into external medium is determined solely by its surface temperature, whereas the radiation of the boundary layer's points is not take into account.
But there are problems in which the exact value of IR energy, emitted by the object, is significant.In this connection, we can mention, as an example, IR thermography and its application to determination of temperature fields of solid objects [14,15].In the case of opaque bodies (metals), the IR radiation streaming away from the body is totally emitted by its surface.So, measuring radiosity of the target object, one can unambiguously judge about its surface temperature.But in the case of dielectric objects it is not so.Outgoing radiation energy flux in this case contains two constituents: the radiation emitted by the surface and the flux, emitted in the boundary layer, which passes through the surface outward.Now, to interpret the measurement data, one should separate the constituents in the total measuring flux.Obviously the needed calculation cannot be performed to be sufficiently exact with the use of diffusion model.
This also applies to the model considered in [7][8][9], in which HCE is obtained starting from the total energy conservation law and the assumption about black boundary was applied.This assumption enabled satisfying simultaneously boundary conditions both for total energy flux q and for radiation intensity and to simplify the computational procedure.But many of real world objects possess the boundary, which is not absolutely black.It not only absorbs and emits IR radiation but also reflects and transmits it.Hence, this model, as well as the diffusion approximation, cannot be applied to such objects and it cannot be used for the inverse problems mentioned before.
The objectives of the paper are mathematical formulation of coupled problem for combined conductive-radiative heat transfer in dielectric layer of participating medium, which emits, absorbs, and scatters IR radiation in the volume and on the boundary; qualitative study of formulated problem depending on thermophysical and radiative properties of medium and layer's thickness; developing an iterative method for solving the coupled problem and studying of its convergence.
Here  > 0 is the layer's half-thickness.The body B is considered as a thermodynamic system consisting of two interacting subsystems, the participating medium and IR radiation, which occupy the same region V.The interaction of the subsystems is realized by absorbing, emitting, and coherent scattering of IR radiation by the layer's medium both in the volume V and on its boundary V = S 1 ∪ S 2 .Here Interaction of the subsystems in the volume V is defined by volumetric radiative parameters: the spectral refraction index  ] , spectral absorption   ] , and emitting   ] and scattering  ] coefficients (further we suppose   ] =   ] ≡  ] ) and by the spectral phase function  ] (,   ).Interaction of the subsystems on the boundary V is defined by its surface radiative parameters: the spectral absorptivity/emittivity  ] and reflectivity  ] .
The layer is situated within nonparticipating external medium with prescribed temperature  ext and optical refraction index  = 1.The subsystems exchange by energy with external medium by convection, emitting, and absorbing of the radiation through the boundary V.The energy exchange of the layer's medium with the external medium through the boundary V is defined by material parameters, coefficient of convective heat exchange, and by spectral absorptivity/emittivity  ] .The energy exchange of the layer's radiation with the external medium through the boundary V is defined by boundary's spectral transmissivity  ] .
1D stationary thermal state of the layer's medium is identically determined by temperature (), depending on transversal coordinate  ∈ [−, ].The state of the layer's radiation is identically determined by spectral intensity or radiation   (, ), depending on  and direction of radiation's propagation.Here  ∈ [−1, 1],  ≡ cos , where  is the angle between a direction and axes .
Under such consideration we can use the equation of energy balance in the layer's medium to obtain HCE: Here   stands for thermal conductivity of the layer's medium.This equation, unlike HCE, used in [7][8][9] and others, do not contain the radiation heat flux   .But instead it contains the source term in its right-hand side which takes into consideration the interaction of the layer's medium with the medium's radiation.
The spectral phase function  ] (,   ) in (2) determines the fraction of radiation incident in direction   and scattered into direction .It satisfies the normality condition [1]: Planck function, in ( 1) and ( 2) determines the intensity of blackbody's radiation in vacuum.Here ℎ ≃ 6.626075 J⋅s is Planck constant,   ≃ 1.380658 J/K is Boltzmann constant, and  0 ≃ 2.99792458 m/s is speed of light in vacuum.
Constituents   () and   ( ] ) in right-hand side of (1) determine, respectively, the local rates of radiation energy emitted and absorbing by the medium in volume V.
in right-hand side of (2) are the scattering integral and radiation inflow caused by its emitting by the medium in the volume V.
We consider the boundary V as translucent one, which diffusely emits, scatters, and absorbs IR radiation.
Radiation with intensity  ] ()| = ,  ∈ (0, 1) and  ext ] ()| = ,  ∈ (−1,0), incident on the boundary plane S 1 out of the volume V and from the external medium, is partly scattered back, partly absorbed by the boundary S 1 , and partly transmitted through this boundary.Let  1] ∈ [0, 1] be a spectral hemispherical reflectivity of the boundary ) are the fraction of radiation scattered back and absorbed by the boundary S 1 , whereas the part of transmitted radiation will be equal to Similarly the boundary plane S 2 , which is characterized by the radiative properties  2] ,  2] , and ), scatters, absorbs, and transmits the radiation of intensities  ] ()| =− ,  ∈ (−1,0) and  ext ] ()| =− ,  ∈ (0,1) incident on it from the volume V and external medium correspondingly.
Unlike [7][8][9], boundary conditions (7) for HCE in our case do not contain the radiation heat flux q  .The first and second terms in the right-hand parts of this relation are the surface energy loss caused correspondingly by the convection and the radiative emission of the boundary plane   .The third and fourth members take into account the energy of internal and external radiation incident on the boundary plane   and absorbed by it.
The first terms in the right-hand side of the relations (8) determine spectral intensities of the radiation emitting by the boundaries S  into the volume V, the second ones are the spectral intensities of IR radiation, reflected by boundary S  into the volume V, and the last one is intensities of the external radiation penetrating through the boundaries S  into the volume V.
The spectral intensities of radiation proceeding out from layer B through the planes S  (layer's spectral radiosity) can be determined as The second term in the right-hand part of the relations (9) takes into account the internal radiation, passing through the boundary into the external medium, whereas the last term in this part determines the spectral intensities of the external radiation, reflecting back into external medium by the boundary.Energy flux, flowing-off layer B through boundary S  , can be determined as At the given values  ext | =  and  ext ] ()| =  of external temperature and radiation intensity on the boundary planes, given thermal parameters   , ℎ, radiative properties of layer's medium  ] ,  ] , (  , ), and  ] ,  ] of its boundary V, (1) and ( 2) together with boundary conditions ( 7) and (8) determine correct formulation of direct boundary-value problem for combined conductive-radiative heat exchange in layer B, which exchanges by energy with external medium by convection and radiation.
Solution of this problem enables determining, with the use of relations ( 9) and ( 10), the angle distribution of radiation intensity and total energy flux emitting by the layer into the external medium.These parameters are important when inverse problems for noncontact determination of temperature field in the layer and/or radiative parameters of the layer's medium are considered.
The integral radiosity of the layer B on the outer side of its boundary plane S  takes the forms in this case:

Analysis of the Problem
At small  (12) is singularly perturbed as at  = 0 its type is being changed.So we will consider function (, ) as piecewise homogeneous with respect to variable : Due to this we will consider instead of (12) the coupled system
Since  2 = (1),  SB  3 0 / ∼   ( 0 ), where   () is the radiative conductivity [3],  = (  /  ( 0 )), so dimensionless parameter  determines the contributions of conductive and radiative heat transfers in total heat flux: when  ≫ 1, the radiation practically does not affect temperature field in the layer.But in the case of  = (1), the contributions of both mechanisms of heat transfer in the layer are comparable.Further we suppose the case  =  (1).
In the case of optically thin layer, when  0 ≫ 1, the problem becomes regularly perturbed, and it can be solved by expansion in small parameter 1/ 0 [10].
The case  0 ≪ 1 can be effectively analyzed quantitatively with the use of Rosseland's approximation [1].But, as it was mentioned before, this approximation does not describe a temperature distribution and radiation fluxes within the boundary layers of thickness about  0 .
There are various methods for quantitative analysis of singularly perturbed problems [10].One of the approaches to regularize the problems is considering them in the interior of domain and in the boundary layer separately with consequent joining of the two obtained solutions.
To implement such approach for problem (30)-(36) it is necessary to develop a method for its numerical solving for layers with thickness about photon mean path  0 .In this case  0 = (1) and the problem becomes regular (unperturbed).It is convenient to use now the length  0 as the spatial scale.

Numerical Study
We tested a convergence of the iterative process (42)-(44) numerically, depending on chosen values of dimensionless parameters , , , ,  1 ,  1 ,  1 .The calculations were stopping when the criterion max{  ,   } ≤  was meted, where  is a given positive constant and In Figure 1, we can see how the dimensionless intensity  of radiation changes with coordinate , depending on the direction of its incidence.The data for two different values of hemispherical emissivity  1 of the boundary S 1 are presented in the figure.
In Figure 2 the distributions of dimensionless temperature () and dimensionless radiation temperature   () in the layer are shown at the same two values of hemispherical emissivity  1 .Radiation temperature is defined as   = (  /2) 1/4 , where   is dimensionless spatial density of radiation   () = ∫ The fluxes  out S 1 and  out V are also shown in Figure 4. Figure 5 illustrates the influence of boundary's S 1 emissivity  1 on directional dependence of dimensionless intensity The data presented on the figures were obtained for the next values of dimensionless parameters:  = 1.2,  = 7,  = 0.1,  = 0.9,  1 = 0.2,  = 1.65,  1 = 0.01,  ext = 0.3,  0 = 1000 K. Conducted numerical studies show that, for the chosen values of the parameters, the iterative process is convergent up to optical thickness of the layer that equals 3 ( = 1.5).
Analyzing Figures 1 and 2, we can conclude that there are substantial differences between values of conductive  and radiative   temperatures in the layer and between values of intensity of radiation incident in directions of positive  ∈ (0, 1) and negative  ∈ (−1, 0) hemispheres.It means that, at chosen values of the dimensionless parameters, the thermal state of the layer is far from a radiative equilibrium.
The data presented in Figures 1-3 show the substantial effect of surface radiative properties, in particular the boundary's hemispherical emissivity, on temperature field and radiation distribution in the layer.Figures 4 and 5 exhibit that these properties considerably affect the radiation emitted by the layer into external medium.It means that angular distribution of the radiation intensity emitting by the layer can be used as informative parameters for nondestructive determination of the volumetric and surface radiative properties of the medium.

Conclusion
The mathematical model for describing of combined conductive-radiative heat transfer in a dielectric layer which emits, absorbs, and scatters IR radiation both in the layer's volume and on its boundary has been considered.The layer is considered as thermodynamic system consisting of two interacting subsystems, the layer's medium and the layer's radiation, in correspondence to the fact that the heat transfer in the layer's volume is described by two interdependent processes, the heat conduction in the layer's medium and radiation transfer.Following this conception the heat conduction equation has been obtained proceeding from the equation of energy balance for layer's medium.In such form, HCE does not involve the radiation heat flux q  but it contains the source-type terms taking into account the interaction of the layer's medium with the radiation subsystem.The boundary conditions for HCE also do not contain the radiation heat flux q  , but they contain the terms describing the interaction of the layer's medium and the radiation on the layer boundary.This interaction is caused by the surface radiation properties of layer's medium, being hemispherical, emissivity/absorptivity, and reflectivity.
Such formulation of the model enables determining the spatial-and-angular distributions of parameters of IR radiation, emitting by semitransparent mediums.This radiation usually contains the energy fluxes emitting both by the boundary and by interior of the layer.This feature of the model makes it perspective for utilization in inverse problems intended on noncontact determination of the temperature fields and/or volumetric and surface radiation properties of such mediums proceeding from the data, obtained by remote measuring of parameters of IR radiation, emitted by the object.Such data can be gathered with the use of modern high-precision measurement techniques, such as IR thermography.
A nonlinear stationary boundary-value problem for coupled integrodifferential heat conduction equation and radiation transfer equation for the layer, which exchanges by energy with external medium by convection and radiation, has been formulated.Features of the problem, depending on the layer's thickness, have been studied with the use of dimensional analysis.It was established that, in the case of optically thick layer, when the layer's thickness is much more of photon-free path, the problem becomes singularly perturbed one.In the inverse case of optically thin layer, the problem is regularly perturbed, and it becomes a regular (unperturbed) one, when the layer's thickness is of order of several photon's free paths.
An iterative method for solving of the unperturbed problem has been developed.The method reduces the initial nonlinear problem to the sequence of linear problems.Due to this the solutions for iterations were obtained in closed form.The conducted numerical studies show that the iterative processes is convergent up to value of optical thickness that equals 3.
The method, in aggregate with the known diffusion approximation, can be used to implement an approach for solving of singularly perturbed problems and for quantitative analysis of the problem for conductive-radiative heat transfer in optically thick layer.Such approaches, consisting in considering of the singularly perturbed problem in the interior of   domain and in the boundary layer separately with consequent joining of the two obtained solutions, are known [10].

Advances in Mathematical Physics
With the use of the developed method the temperature field and distribution of radiation intensity in the layer have been studied.The obtained numerical results exhibit substantial dependence of the IR radiation, emitted by the layer into the external medium, on the temperature field in the layer and on its radiative properties.This makes possible to expect that developed model and method can be used to formulate and solve inverse problems for noncontact temperature testing in dielectrics and for nondestructive determination of its radiation properties (volumetric and surface) on the base of the data obtained by measuring of IR radiation emitted by the layer.

Figure 4
illustrates how the dimensionless energy flux  out , inflowing into external medium from the layer through boundary S 1 , depends on boundary's emissivity  1 .In considering the case the flux  out contains three constituents: convective heat flux  out conv =  1 (| = −  ext | = ), energy flux  out S 1 =  1  4 | = emitted by the boundary S 1 , and energy flux radiation going out of the volume V through the boundary S 1 , where  = arccos() is a zenith angle of radiation propagation.

Figure 4 :Figure 5 :
Figure 4: Dependencies of the dimensionless energy flux  out outgoing from the layer through boundary S 1 and its constituents  out S 1 ,  out V on emissivity  1 .