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

Authors: Vasyl Chekurin; Yurij Boychuk
Journal: Advances in Mathematical Physics
 (2017)

---

## Abstract

The mathematical 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.


---

## Body

## 1. 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 Tr,t, dependent on spatial coordinates r∈V and time t, and spectral radiation intensity Iνr,ω,t being a function of IR radiation frequency ν, spatial coordinates r∈V, direction ω∈S2, and time t. The model also includes coupled relations expressing the conditions, in which the temperature Tr,t and intensity Iνr,ω,t should satisfy on the boundary ∂V of the body. Fundamentals of the theory for radiation transfer in participating heat conducting medium can be found in particular in the monographs [1–3].Nonlinear interconnected integrodifferential equations and nonlinear nonlocal boundary conditions make boundary-value problems, formulated within this model, too complicated for their solving in a closed form. Furthermore, the dimensionality of the problems creates additional computational complexity: in the case that 3D bodies intensityIνr,ω,t is dependent on seven scalar independent variables (frequency ν, three spatial coordinates x,y,z, zenith θ, azimuth φ angles determining direction ω, and time t), temperature Tr,t is dependent on four independent variables (x,y,z, and t). To overcome these complications, one can make some simplifying assumptions concerning the model, restrict the consideration by 1D or 2D geometry, and apply numerical methods.A significant simplification of the problem can be attained in the case, when the contribution of the radiative energy flux into the total medium’s heat flux is a negligible quantity. It can be true at sufficiently low temperatures. In this case the problem becomes uncoupled and a linear one. In fact in this case one can treat temperature field as a prescribed function of spatial coordinate and time and concentrate his effort on solving the radiation transfer problem. Such approximation was used by many authors to study 1D radiation transfer problems (see, for examples, monograph [1] and [4–6]). The principal difficulties in solving the problem in this case are the presence of the scattering integral in RTE and nonlocal boundary conditions, arising when scattering of radiation by the boundary is taken into consideration.There are several effective approaches to overcome this difficulty. Among them is the method of discrete ordinates, in agreement with which RTE is solved for a set of discrete directions, each of which is associated with a solid angle in which the intensity is assumed to be constant [7, 8]. Another method, known as spherical harmonic method, consists in a Legendre expansion of the radiation intensity and the phase function in the radiation transfer equation [4–6].The heat transfer problems in the coupled conductive-radiative formulation are fundamentally nonlinear. Usually HCE, used in “exact” coupled problems, is formulated proceeding from the energy conservation equation [1]. Whereupon HCE contains the divergence of the total energy flux q including two components, conductive heat flux qc and radiative energy flux qr: q=qc+qr. Hence, the boundary conditions for such HCE should be formulated in terms of temperature T and/or total energy flux q. The conductive flux can be presented, in agreement with Fourier law, as qr=-kr∇T or with the use of generalized form, as τ0∂qr/dt+qr=-kr∇T. Here kc is specific heat conductivity, ∇ denotes the gradient operator, and τ0 stands for relaxation time. These two presentations lead to parabolic and hyperbolic HCE correspondingly. To make these equations well defined the flux qr should be determined. In the “exact” formulation one can use the definition of radiation energy flux in terms of radiation intensity as [1]  qr=∫0∞∫4πωIνr,ω,tdωdν. Hence, to determine qr 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–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 radiative-conductive 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 h-=0.1, that is, for optically thin medium. In [8], where the influence of h- 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 thicknessh- of order unity or smaller. The reason is the presence of two spatial scales in the problem: the layer’s physical thickness h and the photon mean free path l0. As a consequence, the coefficient at the highest spatial derivative in HCE decreases with increasing of h-≡h/l0 as 1/h-2. This becomes apparent in occurrence of so-called boundary layer of thickness of order l0 at sufficiently large h-. The temperature and radiation intensity are sharply varied in the boundary layer. Whereupon the problem becomes singularly disturbed at large h-, 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 qr is also taken to be proportional to temperature gradient qr=-kr∇T, where kr is radiation conductivity, which is defined as kr=16/3σSBT3/κ [3], and σSB stands for the Stefan-Boltzmann constant σSB≈5.670367×10-8 W·m2·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 ask=kc+kr, which depends on the temperature. As a consequence, the second spatial scale l0 also vanishes, and the problem remains regular at arbitrary optical thickness h-.This approach was successfully used in various practical problems, in particular, for evaluation of effective heat resistance properties of low-conductivity coats [11–13]. But it is important to realize that the diffusive approximation is valid only at sufficiently low temperature gradients and only in inner medium’s domains, remote from the boundary ∂V on distance exceeding the mean free path l0. In a boundary layer with the thickness of order l0 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–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.


## 2. Governing Equations

A plane isotropic and physically homogeneous dielectric layerB occupying a region V=x,y,z:x∈-∞,∞,y∈-b,b,z∈-∞,∞ is considered. Here b>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=S1∪S2. Here S1=x,y,z:x∈-∞,∞,y=b,z∈-∞,∞,  S2=x,y,z:x∈-∞,∞,y=-b,z∈-∞,∞.Interaction of the subsystems in the volumeV is defined by volumetric radiative parameters: the spectral refraction index nν, spectral absorption κνa, and emitting κνe and scattering σν coefficients (further we suppose κνe=κνa≡κν) and by the spectral phase function pνμ,μ′. 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 temperatureText and optical refraction index n=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 temperatureTy, depending on transversal coordinate y∈-b,b. The state of the layer’s radiation is identically determined by spectral intensity or radiation Iλy,μ, depending on y and direction of radiation’s propagation. Here μ∈-1,1, μ≡cos⁡θ, where θ is the angle between a direction and axes y.Under such consideration we can use the equation of energy balance in the layer’s medium to obtain HCE:(1)-ddykcdTydy=2π∫-11∫0∞κνIνy,μ-nν2IνbTydνdμ.Here kc stands for thermal conductivity of the layer’s medium.This equation, unlike HCE, used in [7–9] and others, do not contain the radiation heat flux qr. 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.Radiation intensityIνy,μ satisfies, in the region V, (2)μ∂Iνy,μ∂y+βνIνy,μ=κνnν2IνbTy+12σν∫-11pνμ′,μIνy,μ′dμ′,where βν=κν+σν stands for spectral extinction coefficient.It can be easily shown, that, proceeding from (1) and (2), HCE, used in [7–9], can be obtained. To do this, one should integrate (2) with respect to variables μ∈-1,1 and φ∈0,2π and sum the obtained equation with (1).The spectral phase functionpνμ,μ′ in (2) determines the fraction of radiation incident in direction μ′ and scattered into direction μ. It satisfies the normality condition [1]:(3)12∫-11pμ′,μdμ=1.Planck function, (4)IνbT≡2hν3c0exp⁡hν/kBT-1,in (1) and (2) determines the intensity of blackbody’s radiation in vacuum. Here h≃6.626075 J·s is Planck constant, kB≃1.380658 J/K is Boltzmann constant, and c0≃2.99792458 m/s is speed of light in vacuum.ConstituentsQeT and QaIν(5)QeT≡2π∫-11∫0∞κνnν2IνbTydνdμ,QaIν≡2π∫-11∫0∞κνIνy,μdνdμ,in right-hand side of (1) determine, respectively, the local rates of radiation energy emitted and absorbing by the medium in volume V.TermsSνsc and SνeT(6)SνscIν=12σν∫-11pμ′,μIνy,μ′dμ′,SνeT=nν2κνIνbT,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 intensityIνμy=b,μ∈0,1 and Iνextμy=b,μ∈-1,0, incident on the boundary plane S1 out of the volume V and from the external medium, is partly scattered back, partly absorbed by the boundary S1, and partly transmitted through this boundary. Let ρ1ν∈0,1 be a spectral hemispherical reflectivity of the boundary S1;  ε1ν∈0,1 is its spectral hemispherical absorptivity. Then ρ1ν and ε1ν1-ρ1ν are the fraction of radiation scattered back and absorbed by the boundary S1, whereas the part of transmitted radiation will be equal to τ1ν≡1-ρ1ν-ε1ν1-ρ1ν=1-ε1ν1-ρ1ν.Similarly the boundary planeS2, which is characterized by the radiative properties ρ2ν, ε2ν, and τ2ν=1-ε2ν1-ρ2ν, scatters, absorbs, and transmits the radiation of intensities Iνμy=-b,μ∈-1,0 and Iνextμy=-b,μ∈0,1 incident on it from the volume V and external medium correspondingly.The intensities of radiation emitting by the planesS1,  S2 into the volume V and in external medium are defined as I1νe,outμ=ε1IνBTy=b,μ∈0,1 and I1νe,inμ=ε1nν2IνBTy=b,μ∈-1,0,  I2νe,outμ=ε2IνBTy=-b,μ∈-1,0 and I2νe,inμ=ε2n2IνBTy=-b,μ∈0,1, correspondingly.Now, applying the energy conservation law, we can establish the conditions (7) and (8), which dependent variables Ty and Iνy,μ satisfy on the boundary planes S1 and S2:(7)-kcdTdyy=bk=-h1Ty=bk-Texty=bk-π∫0∞1+nν2εkνIνbTy=bkdν+2π∫μ∈mk∫0∞εkν1-ρkνIνμy=bkdνμdμ+2π∫μ∈ml∫0∞εkν1-ρkνIνextμy=bkdνμdμ,(8)Iνμy=bk=nν2εkνIνbTy=bk+ρkν∫μ′∈mkIνμ′y=bkdμ′+1-εkν1-ρkν∫μ′∈mlIνextμ′y=bkdμ′,μ∈ml.Here we used denotations: k,l=1,2;  l≠k,  b1=b,  b2=-b. m1=0,1,  m2=-1,0.Unlike [7–9], boundary conditions (7) for HCE in our case do not contain the radiation heat flux qr. 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 Sk. The third and fourth members take into account the energy of internal and external radiation incident on the boundary plane Sk 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 Sk into the volume V, the second ones are the spectral intensities of IR radiation, reflected by boundary Sk into the volume V, and the last one is intensities of the external radiation penetrating through the boundaries Sk into the volume V.The spectral intensities of radiation proceeding out from layerB through the planes Sk (layer’s spectral radiosity) can be determined as (9)Ikνoutμ=εkνIνbTy=bk+1-εkν1-ρkνIνμy=bk+ρkν∫μ′∈mlIνextμ′y=bkdμ′,μ∈mk.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 layerB through boundary Sk, can be determined as(10)qkout=hTy=bk-Texty=bk+2π∫μ∈mk∫0∞I1νoutdνμdμ.At the given values Texty=bk and Iνextμy=bk of external temperature and radiation intensity on the boundary planes, given thermal parameters kc,h, radiative properties of layer’s medium κν,σν,  pμ′,μ, and εkν,ρkν 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.In further consideration, we will apply the approximation of “gray” medium to the layer [1, 3]. This does not change the types of (1) and (2) and boundary conditions (7), (8) but permits substantially reducing the amount of data needed for definition of the radiative properties of layer B. Additionally we will restrict the consideration by the case of isotropic scattering, assuming the phase function pμ′,μ=1 [1] and consider thermal and radiative properties as independent of temperature T. In this approximation (1) and (2) and boundary conditions (7) and (8) take the forms (11) and (12) and (16)–(19), respectively:(11)kcd2Tydy2=2πκ∫-11n2IbTy-Iy,μdμ,(12)μ1β∂Iy,μ∂y+Iy,μ=κ-n2IbTy+12σ-∫-11Iy,μ′dμ′,(13)-kcdTdyy=bk=-hkTy=bk-Texty=bk-π1+n-2εkIbTy=bk+2πεk1-ρk∫μ∈mkIμy=bkμdμ+2πεk1-ρk∫μ∈mlIextμy=bμdμ,(14)Iμy=bk=εkIbTy=bk+ρk∫μ′∈mkIμ′y=bkdμ′+1-εk1-ρk∫μ′∈mlIextμ′y=bkdμ′,μ∈ml.Here IbT≡n2σSBT4/π, n is integral refraction index, κ and σ stand for integral absorption and scattering coefficients, β≡κ+σ is integral extinction coefficient, σ-≡σ/β stands for scattering albedo, κ-≡κ/β=1-σ-, and  ρk and εk are hemispherical total reflectivity and hemispherical total absorptivity of the boundary Sk;  Iy,μ and Iextμ are the integral intensities of internal and external radiation.Dimensionless parametersκ-∈0,1,  σ-∈0,1,  κ-+σ-=1. We will consider participating media for which κ-=O1,  σ-=O1.The integral radiosity of the layerB on the outer side of its boundary plane Sk takes the forms in this case:(15)Ikoutμ=εkn-2IbTy=bk+1-εk1-ρkIμ′y=bk+ρk∫μ′∈mlIextμ′y=bkdμ′,μ∈mk.


## 3. Analysis of the Problem

At smallμ (12) is singularly perturbed as at μ=0 its type is being changed. So we will consider function Iy,μ as piecewise homogeneous with respect to variable μ:(16)Iy,μ=I1y,μ,μ∈0,1I2y,μ,μ∈-1,0.Due to this we will consider instead of (12) the coupled system (17)μ1β∂I1y,μ∂y+I1y,μ=κ-IbTy+12σ-∫01I1y,μ′dμ′+∫-10I2y,μ′dμ′,(18)μ1β∂I2y,μ∂y+I2y,μ=κ-IbTy+12σ-∫01I1y,μ′dμ′+∫-10I2y,μ′dμ′.In the new denotations boundary conditions (14) for (17), (18) take the following form: (19)Ikμy=bl=εlIbTy=bl+ρl∫μ′∈mlIlμ′y=bldμ′+1-εl1-ρl∫μ∈mkIlextμ′dμ′.Here Ikextμ≡Iextμy=bkWe also rewrite (11) in equivalent form:(20)kcd2Tydy2=2πκ2IbTy-∫01I1y,μdμ+∫-10I2y,μdμ.The nonlinear second-order differential (20) should be subordinated to both boundary conditions (13). With the use of the new dependent variables these conditions take the following forms:(21)-kcdTdyy=bk-hkTy=bk-Tkext-π1+n-2εkIbTy=bk+2πεk1-ρk∫μ∈mkIkμy=bkμdμ+2πεk1-ρk∫μ∈mlIkextμy=bkμdμ.We introduce dimensionless variables: coordinatey-=y/b, temperature T-y-=Tby-/T0, and intensity Jy-,μ=Iby-,μ/IbT0, where T0 is a characteristic temperature. With these variables equations (20), (17), and (18) and boundary conditions (21), (19) can be rewritten as(22)ε02d2T-y-dy-2=2κ-2k-2T-y-4-∫01J1y-,μdμ+∫-10J2y-,μdμ(23)με0∂Jky-,μ∂y-+Jky-,μ=κ-T-4y-+12σ-∫01J1y-,μ′dμ′+∫-10J2y-,μ′dμ′(24)-k-ε0dT-dy-y-=±1=-h-kκ-T-y-=±1-T-exty-=±1-1+n-2κ-εkT-4y-=±1+2εk1-ρkκ-∫μ∈mkJkμy-=±1μdμ+2εk1-ρkκ-∫μ∈mlJkextμy-=±1μdμ,(25)Jkμy-=∓1=εlT-4y-=∓1+ρl∫μ′∈mlJlμ′y-=∓1dμ′+1-εl1-ρl∫μ′∈mkJlextμ′dμ′.Here ε0≡1/βb=l0/b,  k-≡kcκ/n2σSBT03,  h-1≡h1/n2σSBT03. In relations (24) and (25), containing the expression y=±1 or y=∓1, the upper sign should be taken at k=1, whereas the lower one – should be taken at k=2.Sincen2=O1,  σSBT03/κ~krT0, where krT is the radiative conductivity [3], k-=Okc/krT0, so dimensionless parameter k- determines the contributions of conductive and radiative heat transfers in total heat flux: when k-≫1, the radiation practically does not affect temperature field in the layer. But in the case of k-=O1, the contributions of both mechanisms of heat transfer in the layer are comparable. Further we suppose the case k-=O1.Parameterε0 determines the ratio of photon’s mean free path l0≡1/β to layer’s half-thickness b. Hence for layers with b≫l0, ε0 is a small parameter. Value of extinction coefficients β for solid dielectrics varies in the bounds (500⋯10000) m−1. Hence the relation b≫l0 becomes valid for layers with b~1⋯10 mm already. So, we can conclude, that, at b≫l0 (the case of optically thick layer), (22)–(23) contain the small parameter at highest derivative; hence the problem (22)–(25) is singularly perturbed [10].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 l0.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 l0. In this case ε0=O1 and the problem becomes regular (unperturbed). It is convenient to use now the length l0 as the spatial scale. Introducing dimensionless spatial coordinate ξ=βy=y/l0, we rewrite the problem (30)–(36) in the following form:(26)d2T-ξdξ2=2κ-λ-2T-ξ4-∫01J1ξ,μdμ+∫-10J2ξ,μdμ(27)μ∂Jkξ,μ∂ξ+Jkξ,μ=κ-T-ξ4+12σ-∫01J1ξ,μ′dμ′+∫-10J2ξ,μ′dμ′(28)-dT-dξξ=b-k=-γ-kT-ξ=b-k-T-extξ=b-k-1+n-2λ--1εkT-4ξ=b-k+2εk1-ρkλ--1∫μ∈mkJkμξ=b-kμdμ+2εk1-ρkλ--1∫μ∈mlJkextμξ=b-kμdμ,(29)Jkμξ=b-l=εlT-4ξ=b-l+ρl∫μ′∈mlJlμ′ξ=b-ldμ′+1-εl1-ρl∫μ′∈mkJlextμ′dμ′,where λ-≡kcβ/n2σSBT03,  b-1=b-,  b-2=-b-,  b-≡βb=b/l0,  γ1≡h1/kcβ,γ2≡h2/kcβ.


## 4. An Iterative Method for Solving of the Unperturbed Problem

We will develop an iterative method for solving the problem (26)–(29). We assume Jext≡0 and modify the boundary conditions (28), (29) for k=2. We suppose boundary S2 as a blackbody with prescribed fixed temperature T0, what enables us to rewrite them in the following form:(30)T-ξ=-b-=1,(31)J1μξ=-b-=1.The boundary conditions (28) and (29) for k=1 remain in the following form:(32)-dT-dξξ=b-=-γ-1T-ξ=b--T-extξ=b--1+n-2λ--1ε1T-4ξ=b-+2ε11-ρ1λ--1∫01J1μξ=b-μdμ.(33)J2μξ=b-=ε1T-4ξ=b-+ρ1∫01J1μ′ξ=b-dμ′,We expect to come to the solution of the nonlinear problem (26), (27), (30)–(33), considering sequences of functions T-mξ,J1mξ,μ,J2mξ,μ,m=0,1,2,…, which satisfy the next iterative sequence of linear problems:(34)d2T-m+1dξ2=Θmξ,(35)μ∂J1m+1∂ξ+J1m+1=Φmξ,μ,(36)μ∂J2m+1∂ξ+J2m+1=Φmξ,μ(37)dT-m+1dξξ=b-=ϑm,T-m+1ξ=-b-=1,(38)J1m+1ξ=-b-=1,J2m+1ξ=b-=ε1T-k4ξ=b-+2ρ1∫01J1kξ=b-μdμ.Here we used the denotations:(39)Θmξ=2κ-λ-2T-mξ4-∫01J1mξ,μdμ+∫-10J2mξ,μdμ,(40)Φmξ=κ-T-mξ4+12σ-∫01J1mξ,μ′dμ′+∫-10J2mξ,μ′dμ′.(41)ϑm=γ-1T-mξ=b--T-extξ=b-+1+n-2λ--1ε1T-mξ=b-4-2ε11-ρ1λ--1∫01J1mμξ=b-μdμIntegrating ordinary differential equations (34)–(36) and subordinated their solutions to boundary conditions (37) and (38), we obtain(42)T-m+1ξ=∫-b-ξ∫-b-zΘmζdζdz+ϑm-∫-b-bΘmξdξξ+b-+1(43)J1m+1ξ,μ=∫-b-ξΦmze-ξ-z/μμdz+e-ξ+b-/μ(44)J2m+1ξ,μ=∫b-ξΦmze-ξ-z/μμdz+e-ξ-b-/με1T-my-=b-4+2ρ1∫01J1kμ′dμ′ξ=b-.To avoid iterative integration, we approximate the solutions by linear finite elements. To do this we chose on the segment -b-,b-N+1 regular node points ξi=-b-+i·2b-/N,i=0,N¯, and represent the functions T-mξ,  J1(m)ξ,μ and J2(m)ξ,μ by their node values: T-im≡Tmξi,  J1i(m)μ≡J1(m)ξi,μ,  J2i(m)μ≡J2(m)ξi,μ. Approximating T-mξ,  J1(m)ξ,μ and J2(m)ξ,μ on the elements ξi,ξi+1,i=0,N-1¯, by linear functions,(45)T-mη=T-imφ1η+T-i+1mφ2η,J1mη,μ=J1imμφ1η+J1i+1mμφ2η,J2mη,μ=J2imμφ1η+J2i+1mμφ2η,where (46)φ1η=121+η,φ2η=121-η.We replace the integrals in the right-hand sides of relations (42)–(44) with corresponding sums and obtain in the issue(47)T-im+1=Δ24∑j=0i-12∑p=0jΘpm+Θp+1m-Θjm-Θj+1m+ϑim-Δ2∑j=0i-1Θjm+Θj+1mξi+b-+1(48)J1im+1μ=e-ξi/μμ∑j=1i-1Φjmvjμ+wj-1μ+Φ0mv0μ+Φimwi-1μ+μe-b-/μ(49)J2im+1μ=-e-ξi/μμ∑j=i+1N-1Φjmvjμ+wj-1μ+Φimviμ+ΦNmwN-1μ+μeb-/με1T-Nm4+2ρ1∫01J1Nmμ′μ′dμ′.Here,(50)vjμ=-μΔΔ+μeξj/μ-μeξj/μ,wjμ=μΔΔ-μeξj+1/μ+μeξj/μ,Δ≡2bNΘjm≡Θmξi,Φjm≡Φmξi.


## 5. Numerical Study

We tested a convergence of the iterative process (42)–(44) numerically, depending on chosen values of dimensionless parameters b-, λ-, κ-, σ-, ε1, ρ1, γ1. The calculations were stopping when the criterion max⁡εT,εJ≤ε was meted, where ε is a given positive constant and(51)εT≡∑i=0NT-im+1-T-im2,εJ≡∑i=0N∫01J1im+1μ-J1imμdμ+∫-10J2im+1μ-J2imμdμ.In Figure 1, we can see how the dimensionless intensity J 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 S1 are presented in the figure.Figure 1
Dependencies of dimensionless intensity of radiation on spatial coordinateξ for different direction of its incidence at two values of hemispherical emissivity ε1: (а) – ε=0.2; (b) ε=0.5.
(a)
(b)In Figure2 the distributions of dimensionless temperature T-ξ and dimensionless radiation temperature T-rξ in the layer are shown at the same two values of hemispherical emissivity ε1. Radiation temperature is defined as T-r=G-r/21/4, where G-r is dimensionless spatial density of radiation G-rξ=∫01J1ξ,μdμ+∫-10J2ξ,μdμ.Figure 2
Dependencies of dimensionless temperatureT- and radiation temperature T-r on spatial coordinate for two different values of hemispherical emissivity ε1: (а) – ε=0.2; (b) ε=0.5.
(a)
(b)In Figure3 the coordinate dependencies of dimensionless conduction heat flux q-cξ=-λ-∂T-ξ/∂ξ, radiation energy flux q-rξ=∫01J1ξ,μμdμ+∫-10J2ξ,μμdμ, and total flux q-ξ=q-cξ+q-rξ are presented.Figure 3
Dependencies of dimensionless conductiveq-r, radiative q-r, and total q- energy fluxes on coordinate ξ for two different values of hemispherical emissivity ε1: (а) – ε1=0.2; (b) – ε1=0.5.
(a)
(b)Figure4 illustrates how the dimensionless energy flux q-out, inflowing into external medium from the layer through boundary S1, depends on boundary’s emissivity ε1. In considering the case the flux q-out contains three constituents: convective heat flux q-convout=γ1T-ξ=b--T-extξ=b-, energy flux q-S1out=ε1T-4ξ=b- emitted by the boundary S1, and energy flux q-Vout=1-ε21-ρ2∫01J1μξ=b-dμ coming out of the volume V. The fluxes q-S1out and q-Vout are also shown in Figure 4.Figure 4
Dependencies of the dimensionless energy fluxq-out outgoing from the layer through boundary S1 and its constituents q-S1out, q-Vout on emissivity ε1.Figure5 illustrates the influence of boundary’s S1 emissivity ε1 on directional dependence of dimensionless intensity JVoutμ=1-ε21-ρ2J1μξ=b- of radiation going out of the volume V through the boundary S1, where θ=arccos⁡μ is a zenith angle of radiation propagation.Figure 5
Directional dependence ofJVout for different value of boundary’s S1 emissivity ε1.The data presented on the figures were obtained for the next values of dimensionless parameters:b-=1.2,  λ-=7,  κ-=0.1,  σ-=0.9,  ρ1=0.2,  n=1.65,  γ1=0.01,  T-ext=0.3, T0=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 (b-=1.5).Analyzing Figures1 and 2, we can conclude that there are substantial differences between values of conductive T- and radiative T-r 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 Figures1–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.


## 6. 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 fluxqr 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 qr, 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].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.



---

*Source: 9139135-2017-08-16.xml*