EN FR
EN FR
##### I4S - 2016
Application Domains
New Software and Platforms
Bilateral Contracts and Grants with Industry
Partnerships and Cooperations
Bibliography
Application Domains
New Software and Platforms
Bilateral Contracts and Grants with Industry
Partnerships and Cooperations
Bibliography

## Section: Research Program

### Thermal methods

#### Infrared thermography and heat transfer

This section introduce the infrared radiation and its link with the temperature, in the next part different measurement methods based on that principle are presented.

Infrared is an electromagnetic radiation having a wavelength between $0.2$ $\mu m$ and 1 $mm$, this range begin in uv spectrum and it ends on the microwaves domain, see Figure 1.

Figure 1. Electromagnetic spectrum - Credit MODEST, M.F. (1993). Radiative Heat Transfer. Academic Press.

For scientific purpose infrared can be divided in three ranges of wavelength in which the application varies, see Table 1.

 Band name wavelength Uses $\setminus$ definition Near infrared (PIR, IR-A, NIR) $0.7-3\mu$m Reflected solar heat flux Mid infrared (MIR, IR-B) $3-50\mu$m Thermal infrared Far infrared (LIR, IR-C, FIR) $50-1000\mu$m Astronomy

Our work is concentrated in the mid infrared spectral band. Keep in mind that Table 1 represents the ISO 20473 division scheme, in the literature boundaries between bands can move slightly.

The Plank's law, proposed by Max Planck en 1901, allow to compute the black body emission spectrum for various temperatures (and only temperatures), see Figure 2 left. The black body is a theoretical construction, it represents perfect energy emitter at a given temperature, cf Equation (20).

 ${M}_{\lambda ,T}^{o}=\frac{{C}_{1}{\lambda }^{-5}}{{exp}^{\frac{{C}_{2}}{\lambda T}}-1}$ (20)

With $\lambda$ the wavelength in m and $T$ as the temperature in Kelvin. The ${C}_{1}$ an ${C}_{2}$ constant, respectively in W.m${}^{2}$ and m.K are defined as follow:

 $\begin{array}{ccc}\hfill {C}_{1}& =& 2h{c}^{2}\pi \hfill \\ \hfill {C}_{2}& =& h\frac{c}{k}\hfill \end{array}$ (21)

with

• $c$ The electromagnetic wave speed (in vacuum $c$ is the light speed in m.s${}^{-1}$).

• $k=1.381{e}^{-23}$ J.K${}^{-1}$ The Boltzmann (Entropy definition from Ludwig Boltzmann 1873). It can be seen as a proportionality factor between the temperature and the energy of a system.

• $h\approx 6,62606957{e}^{-34}$ J.s The Plank constant. It is the link between the photons energy and their frequency.

Figure 2. Left: Plank's law at various temperatures - Right: Energy spectrum of the atmosphere

By generalizing the Plank's law with the Stefan Boltzmann law ( proposed first in 1879 and then in 1884 by Joseph Stefan and Ludwig Boltzmann) it is possible to address mathematically the energy spectrum of real body at each wavelength dependent of the temperature, the optical condition and the real body properties, which is the base of the infrared thermography.

For example, Figure 2 right presents the energy spectrum of the atmosphere at various levels, it can be seen that the various properties of the atmosphere affect the spectrum at various wavelengths. Other important point is that the infrared solar heat flux can be approximated by a black body at 5523,15 K.

##### Infrared Thermography

The infrared thermography is a way to measure the thermal radiation received from a medium. With that information about the electromagnetic flux it is possible to estimate the surface temperature of the body, see section 3.2.1.1. Various types of detector can assure the measure of the electromagnetic radiation.

Those different detectors can take various forms and/or manufacturing process. For our research purpose we use uncooled infrared camera using a matrix of microbolometers detectors. A microbolometer, as a lot of transducers, converts a radiation in electric current used to represent the physical quantity (here the heat flux).

This field of activity includes the use and the improvement of vision system, like in [3].

#### Heat transfer theory

Once the acquisition process is done, it is useful to model the heat conduction inside the cartesian domain $\Omega$. Note that in opaque solid medium the heat conduction is the only mode of heat transfer. Proposed by Jean Baptiste Biot in 1804 and experimentally demonstrated by Joseph Fourier in 1821, the Fourier Law describes the heat flux inside a solid, cf Equation (22).

 $\begin{array}{cc}\hfill \varphi =k\nabla T& X\in \Omega \end{array}$ (22)

Where $k$ is the thermal conductivity in W.m${}^{-1}$.K ${}^{o}$, $\nabla$ is the gradient operator and $\varphi$ is the heat flux density in Wm${}^{-2}$. This law illustrates the first principle of thermodynamic (law of conservation of energy) and implies the second principle (irreversibility of the phenomenon), from this law it can be seen that the heat flux always goes from hot area to cold area.

An energy balance with respect to the first principle drives to the expression of the heat conduction in all point of the domain $\Omega$, cf Equation (23). This equation has been proposed by Joseph Fourier in 1811.

 $\begin{array}{cc}\hfill \rho C\frac{\partial T\left(X,t\right)}{\partial t}=\nabla ·\left(k\nabla T\right)+P& X\in \Omega \end{array}$ (23)

With $\nabla .\left(\right)$ the divergence operator, $C$ the specific heat capacity in J.kg${}^{-1}$.${}^{o}$K${}^{-1}$, $\rho$ the volumetric mass density in kg. m${}^{-3}$, $X$ the space variable $X=\left\{x,y,z\right\}$ and $P$ a possible internal heat production in W.m${}^{-3}$.

To solve the system (23), it is necessary to express the boundaries conditions of the system. With the developments presented in section 3.2.1.1 and the Fourier's law it is possible, for example, to express the thermal radiation and the convection phenomenon which can occur at $\partial \Omega$ the system boundaries, cf Equation (24).

 $\begin{array}{cc}\hfill \varphi =k\nabla T·n=\underset{\text{Convection}}{\underbrace{h\left({T}_{fluid}-{T}_{Boundarie}\right)}}+\underset{\text{Radiation}}{\underbrace{ϵ{\sigma }_{s}\left({T}_{environement}^{4}-{T}_{Boundary}^{4}\right)}}+{\varphi }_{0}& X\in \partial \Omega \end{array}$ (24)

Equation (24) is the so called Robin condition on the boundary $\partial \Omega$, where $n$ is the normal, $h$ the convective heat transfer coefficient in W.m${}^{-2}$.K${}^{-1}$ and ${\varphi }_{0}$ an external energy contribution W.m${}^{-2}$, in cases where the external energy contribution is artificial and controlled we call it active thermography (spotlight etc...) in the contrary it is called passive thermography (direct solar heat flux).

The systems presented in the different sections above (3.2.1 to 3.2.2) are useful to build physical models in order to represents the measured quantity. To estimate key parameters, as the conductivity, one way to do is the model inversion, the next section will introduce that principle.

#### Inverse model for parameters estimation

Lets take any model $A$ which can for example represent the conductive heat transfer in a medium, the model is solved for a parameter vector $P$ and it results another vector $b$, cf Equation (25). For example if $A$ represents the heat transfer, $b$ can be the temperature evolution.

 $AP=b$ (25)

With $A$ a matrix of size $n×m$, $P$ a vector of size $m$ and $b$ of size $n$, preferentially $n>>P$. This model is called direct model, the inverse model consist to find a vector $P$ which satisfy the results $b$ of the direct model. For that we need to inverse the matrix $A$, cf Equation (26).

 $P={A}^{-1}b$ (26)

Here we want find the solution $AP$ which is closest to the acquired measures $M$, Equation (27).

 $AP\approx ℳ$ (27)

To do that it is important to respect the well posed condition established by Jacques Hadamard in 1902

• A solution exists.

• The solution is unique.

• The solution's behavior changes continuously with the initial conditions.

Unfortunately those condition are rarely respected in our field of study. That is why we dont solve directly the system (27) but we minimise the quadratic coast function (28) which represents the Legendre-Gauss least square algorithm for linear problems.

 $mi{n}_{P}\left({\parallel AP-ℳ\parallel }^{2}\right)=mi{n}_{P}\left(ℱ\right)$ (28)

Where $ℱ$ can be a product of matrix.

 $ℱ={\left[AP-ℳ\right]}^{T}\left[AP-ℳ\right]$ (29)

In some case the problem is still ill-posed and need to be regularized for example using the Tikhonov regularization. An elegant way to minimize the cost function $ℱ$ is compute the gradient, Equation (30) and find where it is equal to zero.

 $\nabla ℱ\left(P\right)=2\left[-\frac{\partial A{P}^{T}}{\partial P}\right]\left[AP-ℳ\right]=2J{\left(P\right)}^{T}\left[AP-ℳ\right]$ (30)

Where $J$ is the sensitivity matrix of the model $A$ to its parameter vector $P$.

Until now the inverse method proposed is valid only when the model $A$ is linearly dependent of its parameter $P$, for the heat equation it is the case when you want to estimate the external heat flux, ${\varphi }_{0}$ in equation 24. For all the other parameters, like the conductivity $k$ the model is non-linearly dependant of its parameter $P$. For such case the use of iterative algorithm is needed, for example the Levenberg-Marquardt algorithm, cf Equation (31).

 ${P}^{k+1}={P}^{k}+{\left[{\left({J}^{k}\right)}^{T}{J}^{k}+{\mu }^{k}{\Omega }^{k}\right]}^{-1}{\left({J}^{k}\right)}^{T}\left[ℳ-A\left({P}^{k}\right)\right]$ (31)

Equation (31) is solved iteratively at each loop $k$. Some of our results with such linear or non linear method can be seen in [4] or [2], more specifically [1] is a custom implementation of the Levenberg-Marquardt algorithm based on the adjoint method (developed by Jacques Louis Lions in 1968) coupled to the conjugate gradient algorithm to estimate wide properties field in a medium.