CN105675646A - Intrinsic photothermal information based method for simultaneous measurement of absorption coefficient and thermal conductivity of high-temperature translucent medium - Google Patents

Abstract

The invention belongs to the technical field of physical property measurement of translucent media and discloses an intrinsic photothermal information based method for simultaneous measurement of absorption coefficient and thermal conductivity of a high-temperature translucent medium. The method includes: in a measurement process, using a heater for heating the translucent medium to a certain high temperature; using a detector for measuring intrinsic temperature response and directional radiation intensity of a to-be-detected sample; adopting an inverse problem solution technique for indirectly obtaining thermal conductivity and band absorption coefficient of the to-be-detected sample along with temperature variations. By establishment of conduction-radiation coupled heat transfer forward and inverse problem models of the thermal conductivity and band absorption coefficient of the translucent medium along with temperature variations, under the premise that other parameters of the medium are known, a method of adopting a particle swarm optimization algorithm for simultaneous inversion of the band absorption coefficient and temperature related thermal conductivity of the high-temperature translucent medium is provided, and precision in measurement of thermophysical properties of the translucent medium is improved. The problem that existing measurement methods are large in result error and require more measurement information for inversion of temperature related thermophysical properties is solved.

Description

Based on the method that intrinsic light and heat information measures high temperature translucent medium thermal conductivity and absorptance simultaneously

Technical field

The present invention relates to the method simultaneously obtaining high temperature translucent medium temperature correlation thermal conductivity and band absorption coefficient, belong to translucent medium physical measurement technical field.

Background technology

High temperature translucent medium radiation physical property and thermal physical property parameter be translucent medium is analyzed in its application process, important parameter needed for Design and optimization. In recent years, along with developing rapidly of the modern high technologies such as the infrared characteristic of Aero-Space, infrared acquisition, target and environment, laser, electronic device, biomedicine, the translucent medium physical parameter that varies with temperature when high temperature, multidimensional becomes particularly important. The research carrying out participating medium heat radiation physical property and related discipline is respectively provided with significance for dual-use field.

The research of phase-change heat-storage material has great importance for thermal energy storage. Wherein the research for its thermal conductivity and absorptance is particularly important. Deeply understand this thermal physical property parameter and it is carried out experiment measuring and theory analysis in the field such as material science and solar energy thermal-power-generating, also there is important using value. And under normal circumstances, absorptance and spectral correlation, and thermal conductivity can change with material temperature. Therefore, for the measurement of the thermal conductivity varied with temperature and band absorption coefficient in actual application by significant.

Due in actual measurement process, there is certain measurement error in experimental facilities, be used alone light in some cases or thermal information can not to complete the resultant error measured or obtain of radiant heat physical property relatively big, and the inverting for the hot physical property of temperature correlation needs more metrical information.

Summary of the invention

A kind of method simultaneously measuring high temperature translucent medium temperature Thermal Conductivity Varying With Temperature and band absorption coefficient based on intrinsic light and heat information of offer is provided, to solve existing measuring method, to there is resultant error relatively big, and for problem that the inverting of the hot physical property of temperature correlation needs more metrical information.

The present invention solves that above-mentioned technical problem adopts the technical scheme that:

First the translucent medium will measured heats to temperature T₀, then use the radiant intensity of detector measurement sample upper surface different directions. And use the Temperature Distribution of diverse location in thermocouple measurement medium. Using the measured result initial condition as inverse problem, then utilize the method that Particle Swarm Optimization obtains high temperature translucent medium temperature correlation thermal conductivity and band absorption coefficient simultaneously.

A kind of method simultaneously measuring high temperature translucent medium temperature Thermal Conductivity Varying With Temperature and band absorption coefficient of the present invention, concretely comprises the following steps:

Step one, by testing sample heating to a certain fixed temperature T₀;

Step 2, stops heating, uses detector to measure the radiant intensity q (θ of its different directions at the upper surface of testing sample_j), and use the temperature T within thermocouple measurement testing sample_i, wherein, θ_jDiverse location residing for the radiation direction respectively different with i and thermocouple; As shown in Figure 1.

Step 3, utilizes reverse temperature intensity method, it is assumed that three key band absorptances of testing sample are κ_λ1, κ_λ2And κ_λ3, assume that the thermal conductivity that testing sample varies with temperature is λ (T)=a simultaneously₁+a₂·T+a₃·T²; The unit of λ (T) is W/ (m K); A in formula₁、a₂、a₃Represent three warm variable coefficients of thermal conductivity respectively;

Step 4, solves radiation transfer equation and Heat Conduction Differential Equations, it is thus achieved that testing sample interior temperature distribution T_i,est, lower footnote est represents value of calculation;

Step 5, utilizes the testing sample internal temperature T that step 2 obtains_iValue of calculation T corresponding with step 4_i,est, in conjunction with formula:

$F_{1,obj}=\frac{1}{2}\underset{i=1}{\overset{n}{\Σ}}{\[T_{i,est}/T_i-1.0\]}^2---\left(1\right)$

Obtain the object function F in reverse temperature intensity algorithm_1,obj, wherein n is total thermocouple number;

Step 6, it is judged that whether the object function in step 5 is less than setting threshold epsilon₁, thermal conductivity λ (the Τ)=a of testing sample that if so, then will assume in step 3₁+a₂·T+a₃·T²W/ (m K) exports as a result and (namely exports a₁、a₂、a₃Value), otherwise adopt Particle Swarm Optimization correction sample thermal conductivity and three key band absorptances, return step 4;

Step 7, warm Thermal Conductivity Varying With Temperature step 6 exported is as final sample thermal conductivity, and three key band absorptances step 6 exported are κ_λ1, κ_λ2And κ_λ3Initial value as Particle Swarm Optimization;

Step 8, solves radiation transfer equation and Heat Conduction Differential Equations, it is thus achieved that calculate the radiant intensity q of testing sample upper surface different directions_est(θ_j), lower footnote est represents value of calculation;

Step 9, utilizes the different directions radiant intensity q (θ measuring gained in step 2_j) with step in corresponding value of calculation q_est(θ_j), in conjunction with formula:

$F_{2,obj}=\frac{1}{2}\underset{j=1}{\overset{m}{\Σ}}{\[q_{est}\left({\mathrm{\θ}}_j\right)/q\left({\mathrm{\θ}}_j\right)-1.0\]}^2---\left(2\right)$

Obtain the object function F in Particle Swarm Optimization_2,obj; Wherein, m is the direction number measuring the radiant intensity obtained;

Step 10, it is judged that whether the object function in step 8 is less than setting threshold epsilon₂, if so, then by the band absorption coefficient κ of the testing sample of acquisition_λ1, κ_λ2And κ_λ3Export as a result, the method completing simultaneously to obtain high temperature translucent sample temperature Thermal Conductivity Varying With Temperature and band absorption coefficient, otherwise adopt three key band absorptances of Particle Swarm Optimization correction, and warm Thermal Conductivity Varying With Temperature step 6 exported is as sample thermal conductivity, returns step 8.

Step 4 and step 8 obtain the method in the temperature field in computational fields: utilize Heat Conduction Differential Equations

$\mathrm{\λ}\left(T\right)\frac{{\∂}^{2}T}{\∂z^2}-\frac{\∂q^r}{\∂z}=0---\left(3\right)$

$q_{w1}^r-\mathrm{\λ}\frac{\∂T}{\∂z}{|}_{z=0}=h_{w1}(T_{\mathrm{\∞}}-T_{w1})---\left(4\right)$

$-\mathrm{\λ}\frac{\∂T}{\∂z}{|}_{z=L}=h_{w2}(T_{\mathrm{\∞}}-T_{w2})---\left(5\right)$

Realize, wherein ρ and c_pRepresenting density and the specific heat capacity of testing sample respectively, T and h represents sample temperature and convection transfer rate respectively;Q^rRepresenting radiant heat flux density, wherein footnote w1 and w2 represents coboundary and the lower boundary of testing sample respectively; T_∞Represent the temperature of surrounding fluid.

Step 4 and step 8 obtain the method for the radiation field intensity in computational fields: utilize radiation transfer equation

$\frac{{\mathrm{dI}}_k(z,\mathrm{\θ})}{dz}=-{\mathrm{\κ}}_{\mathrm{\λ}k}{I}_k\left(z\right)+{\mathrm{\κ}}_{\mathrm{\λ}k}{I}_{b,\mathrm{\λ}k}\left(z\right)---\left(6\right)$

Realize, κ in formula_λkRepresenting the kth band absorption coefficient of testing sample, k represents different bands of a spectrum, and I represents radiant intensity, I_bRepresenting the radiant intensity of black matrix at identical temperature, z represents thickness of sample; I_b,λkRepresent the radiant intensity of the kth bands of a spectrum of black matrix;

The method obtaining the heat flow density in Heat Conduction Differential Equations is: utilize equation

$q_{w1}^r=\mathrm{\ϵ}\[{\mathrm{\σT}}_{w1}^4-\underset{k=1}{\overset{3}{\Σ}}{\∫}_{\cos \mathrm{\θ}>0}2{\mathrm{\πI}}_k(L,\mathrm{\θ})\left|\cos \mathrm{\θ}\right|\sin \mathrm{\θ}d\mathrm{\θ}\]---\left(7\right)$

$\frac{\∂q^r}{\∂z}=\underset{k=1}{\overset{3}{\Σ}}4\mathrm{\π}\·{\mathrm{\κ}}_{\mathrm{\λ}k}\·\[I_{b,\mathrm{\λ}k}\left(z\right)-\frac{1}{2}{\∫}_0^{\mathrm{\π}}{I}_k(z,\mathrm{\θ})d\mathrm{\θ}\]---\left(8\right)$

Realizing, in formula, ε represents the emissivity of testing sample upper surface, and σ represents this special fence-Boltzmann constant, I_b,λkRepresent the radiant intensity of the kth bands of a spectrum of black matrix.

The invention has the beneficial effects as follows:

The present invention provides a kind of method that experiment obtains high temperature translucent medium temperature correlation thermal conductivity and band absorption coefficient in conjunction with inversion algorithm simultaneously. The inventive method is based on intrinsic light and heat information and realizes. The measuring method that the present invention proposes introduces intrinsic light and heat information integration technology on the basis of reverse temperature intensity, it is possible to be greatly improved the precision for translucent medium thermophysical property measurement.

Measurement process use heater by translucent medium heating to a certain high temperature, by temperature-responsive (thermal signal) and own direction radiant intensity (optical signal) of detector measurement testing sample, indirectly obtain, finally by reverse temperature intensity technology, thermal conductivity and the spectral absorptance that testing sample varies with temperature. The present invention is by setting up the forward and inverse problem model of the translucent medium heat conduction radiation coupled and heat-exchange of thermal conductivity and the band absorption coefficient varied with temperature, under the premise that other parameters of medium are known, it is proposed that the method adopting Particle Swarm Optimization Simultaneous Inversion high temperature translucent medium temperature correlation thermal conductivity and band absorption coefficient. Basic ideas are the temperature-responsive and the direction radiant intensity own that record testing sample by experiment, the method simultaneously obtaining high temperature translucent medium temperature correlation thermal conductivity and band absorption coefficient then in conjunction with Particle Swarm Optimization.

The present invention is by setting up band absorption coefficient and the direct problem of translucent medium heat conduction radiation coupled and heat-exchange that thermal conductivity varies with temperature and reverse temperature intensity model, solve high temperature translucent medium temperature correlation thermal conductivity and many bands of a spectrum absorptance can not directly measure problem inaccurate with measurement result, it is proposed that a kind of method simultaneously obtaining high temperature translucent medium temperature correlation thermal conductivity and band absorption coefficient. Advantage is in that: model is simple, it is simple to theory solves; Adopt Particle Swarm Optimization, during this Algorithm for Solving optimization problem, have simple, efficient and sensitivity advantages of higher. This invention provides one method fast and accurately for researching high-temperature translucent medium temperature correlation thermal conductivity and many bands of a spectrum absorptance, and space flight, defense and commercial industry tool are of great significance.

Therefore say that the invention solves existing measuring method exists resultant error relatively greatly, and the problem of the inverting more metrical informations of needs for the hot physical property of temperature correlation. Measurement process use heater by translucent medium heating to a certain high temperature, intrinsic temperature by detector measurement testing sample responds and direction radiant intensity, indirectly obtains, finally by reverse temperature intensity technology, thermal conductivity and the band absorption coefficient that testing sample varies with temperature. By setting up the forward and inverse problem model of the translucent medium heat conduction radiation coupled and heat-exchange of thermal conductivity and the band absorption coefficient varied with temperature, under the premise that other parameters of medium are known, propose the method adopting Particle Swarm Optimization Simultaneous Inversion high temperature translucent medium temperature correlation thermal conductivity and band absorption coefficient, improve the precision for translucent medium thermophysical property measurement.

Accompanying drawing explanation

Fig. 1 is the experimental provision schematic diagram simultaneously obtaining high temperature translucent medium temperature correlation thermal conductivity and band absorption coefficient described in detailed description of the invention one.

Detailed description of the invention

The method simultaneously obtaining high temperature translucent medium temperature correlation thermal conductivity and band absorption coefficient described in detailed description of the invention one, present embodiment, the concrete operation step of the method is:

Step one, by testing sample heating to a certain fixed temperature T₀;

Step 2, stops heating, uses detector to measure the radiant intensity q (θ of its different directions at the upper surface of testing sample_j), and use the temperature T within thermocouple measurement testing sample_i, wherein, θ_jDiverse location residing for the radiation direction respectively different with i and thermocouple; As shown in Figure 1.

Step 3, utilizes reverse temperature intensity method, it is assumed that three key band absorptances of testing sample are κ_λ1, κ_λ2And κ_λ3, assume that the thermal conductivity that testing sample varies with temperature is λ (T)=a simultaneously₁+a₂·T+a₃·T²; The unit of λ (T) is W/ (m K);

Step 4, solves radiation transfer equation and Heat Conduction Differential Equations, it is thus achieved that testing sample interior temperature distribution T_i,est, lower footnote est represents value of calculation;

Step 5, utilizes the testing sample internal temperature T that step 2 obtains_iValue of calculation T corresponding with step 4_i,est, in conjunction with formula:

$F_{1,obj}=\frac{1}{2}\underset{i=1}{\overset{n}{\Σ}}{\[T_{i,est}/T_i-1.0\]}^2---\left(1\right)$

Obtain the object function F in reverse temperature intensity algorithm_1,obj, wherein n is total thermocouple number;

Step 6, it is judged that whether the object function in step 5 is less than setting threshold epsilon₁, thermal conductivity λ (the Τ)=a of testing sample that if so, then will assume in step 3₁+a₂·T+a₃·T²W/ (m K) exports as a result, otherwise adopts Particle Swarm Optimization correction sample thermal conductivity and three key band absorptances, returns step 4;

Step 7, warm Thermal Conductivity Varying With Temperature step 6 exported is as final sample thermal conductivity, and three key band absorptances step 6 exported are κ_λ1, κ_λ2And κ_λ3Initial value as Particle Swarm Optimization;

Step 8, solves radiation transfer equation and Heat Conduction Differential Equations, it is thus achieved that calculate the radiant intensity q of testing sample upper surface different directions_est(θ_j), lower footnote est represents value of calculation;

Step 9, utilizes the different directions radiant intensity q (θ measuring gained in step 2_j) with step in corresponding value of calculation q_est(θ_j), in conjunction with formula:

$F_{2,obj}=\frac{1}{2}\underset{j=1}{\overset{m}{\Σ}}{\[q_{est}\left({\mathrm{\θ}}_j\right)/q\left({\mathrm{\θ}}_j\right)-1.0\]}^2---\left(2\right)$

Obtain the object function F in Particle Swarm Optimization_2,obj; Wherein, m is the direction number measuring the radiant intensity obtained;

Step 10, it is judged that whether the object function in step 8 is less than setting threshold epsilon₂, if so, then by the band absorption coefficient κ of the testing sample of acquisition_λ1, κ_λ2And κ_λ3Export as a result, the method completing simultaneously to obtain high temperature translucent sample temperature Thermal Conductivity Varying With Temperature and band absorption coefficient, otherwise adopt three key band absorptances of Particle Swarm Optimization correction, and warm Thermal Conductivity Varying With Temperature step 6 exported is as sample thermal conductivity, returns step 8.

The high temperature translucent medium transient radiation heat-transfer couple physical model that first present embodiment designs three key band absorptances and thermal conductivity varies with temperature, then set up corresponding mathematical model and method for solving, by measuring Temperature Distribution and the direction radiant intensity own obtaining testing sample, reverse temperature intensity reconstruction is utilized to go out temperature correlation thermal conductivity and the band absorption coefficient of high temperature translucent medium.

Detailed description of the invention two, present embodiment are further illustrating method described in detailed description of the invention one, and step 4 and step 8 obtain the method in the temperature field in computational fields and be: utilize Heat Conduction Differential Equations

$\mathrm{\λ}\left(T\right)\frac{{\∂}^{2}T}{\∂z^2}-\frac{\∂q^r}{\∂z}=0---\left(3\right)$

$q_{w1}^r-\mathrm{\λ}\frac{\∂T}{\∂z}{|}_{z=0}=h_{w1}(T_{\mathrm{\∞}}-T_{w1})---\left(4\right)$

$-\mathrm{\λ}\frac{\∂T}{\∂z}{|}_{z=L}=h_{w2}(T_{\mathrm{\∞}}-T_{w2})---\left(5\right)$

Realize, wherein ρ and c_pRepresenting density and the specific heat capacity of testing sample respectively, T and h represents temperature and convection transfer rate respectively;Q^rRepresenting radiant heat flux density, wherein footnote w1 and w2 represents coboundary and the lower boundary of testing sample respectively.

Detailed description of the invention three, present embodiment are further illustrating method described in detailed description of the invention one or two, and step 4 and step 8 obtain the method for the radiation field intensity in computational fields and be: utilize radiation transfer equation

$\frac{{\mathrm{dI}}_k(z,\mathrm{\θ})}{dz}=-{\mathrm{\κ}}_{\mathrm{\λ}k}{I}_k\left(z\right)+{\mathrm{\κ}}_{\mathrm{\λ}k}{I}_{b,\mathrm{\λ}k}\left(z\right)---\left(6\right)$

Realize, κ in formula_λkRepresenting the kth band absorption coefficient of testing sample, k represents different bands of a spectrum, and I represents radiant intensity, I_bRepresenting the radiant intensity of black matrix at identical temperature, z represents thickness of sample.

Detailed description of the invention four, present embodiment are further illustrating method described in detailed description of the invention one, two or three, and the method obtaining the heat flow density in Heat Conduction Differential Equations is: utilize equation

$q_{w1}^r=\mathrm{\ϵ}\[{\mathrm{\σT}}_{w1}^4-\underset{k=1}{\overset{3}{\Σ}}{\∫}_{\cos \mathrm{\θ}>0}2{\mathrm{\πI}}_k(L,\mathrm{\θ})\left|\cos \mathrm{\θ}\right|\sin \mathrm{\θ}d\mathrm{\θ}\]---\left(7\right)$

$\frac{\∂q^r}{\∂z}=\underset{k=1}{\overset{3}{\Σ}}4\mathrm{\π}\·{\mathrm{\κ}}_{\mathrm{\λ}k}\·\[I_{b,\mathrm{\λ}k}\left(z\right)-\frac{1}{2}{\∫}_0^{\mathrm{\π}}{I}_k(z,\mathrm{\θ})d\mathrm{\θ}\]---\left(8\right)$

Realizing, in formula, ε represents the emissivity of testing sample upper surface, and σ represents this special fence-Boltzmann constant.

Claims (4)

1. the method simultaneously measuring high temperature translucent medium thermal conductivity and absorptance based on intrinsic light and heat information, it is characterised in that concretely comprising the following steps of the method:

Step one, by testing sample heating to a certain fixed temperature T₀;

Step 2, stops heating, uses detector to measure the radiant intensity q (θ of its different directions at the upper surface of testing sample_j), and use the temperature T within thermocouple measurement testing sample_i, wherein, θ_jDiverse location residing for the radiation direction respectively different with i and thermocouple;

Step 3, utilizes reverse temperature intensity method, it is assumed that three key band absorptances of testing sample areWithAssume that the thermal conductivity that testing sample varies with temperature is λ (T)=a simultaneously₁+a₂·T+a₃·T²; The unit of λ (T) is W/ (m K); A in formula₁、a₂、a₃Represent three warm variable coefficients of thermal conductivity respectively;

Step 4, solves radiation transfer equation and Heat Conduction Differential Equations, it is thus achieved that testing sample interior temperature distribution T_i,est, lower footnote est represents value of calculation;

Step 5, utilizes the testing sample internal temperature T that step 2 obtains_iValue of calculation T corresponding with step 4_i,est, in conjunction with formula:

$F_{1,obj}=\frac{1}{2}\underset{i=1}{\overset{n}{\Σ}}{\[T_{i,est}/T_i-1.0\]}^2---\left(1\right)$

Obtain the object function F in reverse temperature intensity algorithm_1,obj, wherein n is total thermocouple number;

Step 6, it is judged that whether the object function in step 5 is less than setting threshold epsilon₁, thermal conductivity λ (the T)=a of testing sample that if so, then will assume in step 3₁+a₂·T+a₃·T²Export as a result, otherwise adopt Particle Swarm Optimization correction sample thermal conductivity and three key band absorptances, return step 4;

The unit of λ (T) is W/ (m K);

Step 7, three key band absorptances that step 6 exports, as final sample thermal conductivity, are by warm Thermal Conductivity Varying With Temperature step 6 exportedWithInitial value as Particle Swarm Optimization;

Step 8, solves radiation transfer equation and Heat Conduction Differential Equations, it is thus achieved that calculate the radiant intensity q of testing sample upper surface different directions_est(θ_j), lower footnote est represents value of calculation;

Step 9, utilizes the different directions radiant intensity q (θ measuring gained in step 2_j) with step in corresponding value of calculation q_est(θ_j), in conjunction with formula:

$F_{2,obj}=\frac{1}{2}\underset{j=1}{\overset{m}{\Σ}}{\[q_{est}\left({\mathrm{\θ}}_j\right)/q\left({\mathrm{\θ}}_j\right)-1.0\]}^2---\left(2\right)$

Obtain the object function F in Particle Swarm Optimization_2,obj; Wherein, m is the direction number measuring the radiant intensity obtained;

Step 10, it is judged that whether the object function in step 8 is less than setting threshold epsilon₂, if so, then by the band absorption coefficient of the testing sample of acquisitionWithExport as a result, the method completing simultaneously to obtain high temperature translucent sample temperature Thermal Conductivity Varying With Temperature and band absorption coefficient, otherwise adopt three key band absorptances of Particle Swarm Optimization correction, and warm Thermal Conductivity Varying With Temperature step 6 exported is as sample thermal conductivity, returns step 8.

2. the method simultaneously measuring high temperature translucent medium thermal conductivity and absorptance based on intrinsic light and heat information according to claim 1, it is characterised in that step 4 and step 8 obtain the method in the temperature field in computational fields and be: utilize Heat Conduction Differential Equations

$\mathrm{\λ}\left(T\right)\frac{{\∂}^{2}T}{\∂z^2}-\frac{\∂q^r}{\∂z}=0---\left(3\right)$

$q_{w1}^r-\mathrm{\λ}\frac{\∂T}{\∂z}{|}_{z=0}=h_{w1}(T_{\mathrm{\∞}}-T_{w1})---\left(4\right)$

$-\mathrm{\λ}\frac{\∂T}{\∂z}{|}_{z=L}=h_{w2}(T_{\mathrm{\∞}}-T_{w2})---\left(5\right)$

Realize, wherein ρ and c_pRepresenting density and the specific heat capacity of testing sample respectively, T and h represents sample temperature and convection transfer rate respectively; q^rRepresenting radiant heat flux density, wherein footnote w1 and w2 represents coboundary and the lower boundary of testing sample respectively; T_∞Represent the temperature of surrounding fluid.

3. the method simultaneously measuring high temperature translucent medium thermal conductivity and absorptance based on intrinsic light and heat information according to claim 1 and 2, it is characterized in that, step 4 and step 8 obtain the method for the radiation field intensity in computational fields and are: utilize radiation transfer equation

$\frac{{\mathrm{dI}}_k(z,\mathrm{\θ})}{dz}=-{\mathrm{\κ}}_{\mathrm{\λ}k}{I}_k\left(z\right)+{\mathrm{\κ}}_{\mathrm{\λ}k}{I}_{b,\mathrm{\λ}k}\left(z\right)---\left(6\right)$

Realize, κ in formula_λkRepresenting the kth band absorption coefficient of testing sample, k represents different bands of a spectrum, and I represents radiant intensity, I_bRepresenting the radiant intensity of black matrix at identical temperature, z represents thickness of sample; I_b,λkRepresent the radiant intensity of the kth bands of a spectrum of black matrix.

4. the method simultaneously measuring high temperature translucent medium thermal conductivity and absorptance based on intrinsic light and heat information according to claim 3, it is characterised in that the method obtaining the heat flow density in Heat Conduction Differential Equations is: utilize equation

$q_{w1}^r=\mathrm{\ϵ}\[{\mathrm{\σT}}_{w1}^4-\underset{k=1}{\overset{3}{\mathrm{\Σ}}}{\∫}_{cos\mathrm{\θ}>0}2{\mathrm{\πI}}_k(L,\mathrm{\θ})\left|cos\mathrm{\θ}\right|sin\mathrm{\θ}d\mathrm{\θ}\]---\left(7\right)$

$\frac{\∂q^r}{\∂z}=\underset{k=1}{\overset{3}{\Σ}}4\mathrm{\π}\·{\mathrm{\κ}}_{\mathrm{\λ}k}\·\[I_{b,\mathrm{\λ}k}\left(z\right)-\frac{1}{2}{\∫}_0^{\mathrm{\π}}{I}_k(z,\mathrm{\θ})d\mathrm{\θ}\]---\left(8\right)$

Realizing, in formula, ε represents the emissivity of testing sample upper surface, and σ represents this special fence-Boltzmann constant, I_b,λkRepresent the radiant intensity of the kth bands of a spectrum of black matrix.