EN FR
EN FR


Section: Research Program

Thermal methods

Infrared vision and measurement Method

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 radiation

Infrared is an electromagnetic radiation having a wavelength between 0.7 μm and 1 mm, this range begin early after the visible 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.
IMG/SpectreElectro.png

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

Table 1. Wavelength bands in the infrared according to ISO 20473:2007
Band name wavelength Uses definition
Near infrared (PIR, IR-A, NIR) 0.7-3μm Reflected solar heat flux
Mid infrared (MIR, IR-B) 3-50μm Thermal infrared
Far infrared (LIR, IR-C, FIR) 50-1000μ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 transmitter at a given temperature, cf Equation (20 ).

M λ , T o = C 1 λ - 5 exp C 2 λ T - 1 (20)

With λ the wavelength in m and T as the temperature in Kelvin. The C1 an C2 constant, respectively in kg.m4.s-3 and m.K are defined as follow:

C 1 = 2 h c 2 π C 2 = h c k (21)

with

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

  • k=1.381e-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.

  • h6,62606957e-34 J.s The Plank constant. It's the link between the photons energy and their frequency.

Figure 2. Left: Plank's law at various temperatures - Right: Energy spectrum of the atmosphere
IMG/Planck.png IMG/Solar_spectrum_en.png

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's 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, cf Equation (22 ).

ϕ λ = σ s T o 4 Balck body emission - ϵ λ σ s T 4 Real body emission - 1 - ϵ λ σ s T o 4 reflexion Real body radiosity ϕ λ = ϵ λ σ s T o 4 - T 4 (22)

where

  • σs=5,67.10-8 W.m-2.K-4 The Stephan constant

  • ϵλ The emissivity ratio of the real body at a given wavelength (For the black body ϵλ=1 )

  • ϕλ The energy flux in Wm-2 at a given wavelength

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.

The above demonstrations are here to illustrate how energy and temperature can be linked, this can only be used as an general purpose introduction for the next sections.

Infrared Thermography

The infrared thermography is a way to measure the thermal radiation emitted or reflected by a medium. With that information about the electromagnetic flux it's possible to compute the surface temperature of the body, see section 3.2.1.1 . Various types of detector can assure the measure of the electromagnetic radiation. Table 2 illustrates for various spectral bands in infrared some of the detectors which can be used to measure the thermal radiation.

Table 2. Infrared sensors
Bands Wavelength Detectors
Near Infrared 0 . 7 - 1 μ m End of the human eye response, Aluminium gallium arsenide detectors (InGaAs)
Short Wave Infrared - SWIR\Band I 1 - 3 μ m Germanium detector ans InGaAs
Mid Wave Infrared - MWIR\Band II 3 - 5 μ m Mercury cadmium telluride detectors (HgCdTe), Lead selenide (PbSe) and InSb
Long wave Infrared - LWIR\Band III 7 - 14 μ m Microbolometers (VOx) and HgCdTe
Very long wave infrared 12 - 30 μ m silicon detectors
Figure 3. Left: AX5 low cost camera - Right: Example of measurement process [6]
IMG/flir_a65.jpg IMG/BilanRadi.png

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). Figure 3 shows an example of FLIR AX5 low cost mini infrared camera that we use for long term SHM purpose, and also a schematic representation of the measurement principle.

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

Heat transfer theory

Once the acquisition process is done, it's useful to model the heat conduction inside the cartesian domain Ω. 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 (23 ).

ϕ = k T X Ω (23)

Where k is the thermal conductivity in W.m-1.K o, is the gradient operator and ϕ 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 Ω, cf Equation (24 ). This equation has been proposed by Joseph Fourier in 1811.

ρ C T ( X , t ) t = · k T + P X Ω (24)

With .() the divergence operator, C the specific heat capacity in J.kg-1.oK-1, ρ the volumetric mass density in kg. m-3, X the space variable X={x,y,z} and P a possible internal heat production in W.m-3.

To solve the system (24 ), it's 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's possible, for example, to express the thermal radiation and the convection phenomenon which can occur at Ω the system boundaries, cf Equation (25 ).

ϕ = k T · n = h T f l u i d - T B o u n d a r i e Convection + ϵ σ s T e n v i r o n e m e n t 4 - T B o u n d a r y 4 Radiation + ϕ 0 X Ω (25)

Equation (25 ) is the so called Robin condition on the boundary Ω, where n is the normal, h the convective heat transfer coefficient in W.m-2.K-1 and ϕ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's 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 (26 ). For example if A represents the heat transfer, b can be the temperature evolution.

A P = b (26)

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 (27 ).

P = A - 1 b (27)

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

A P (28)

To do that it's 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's why we dont solve directly the system (28 ) but we minimise the quadratic coast function (29 ) which represents the Legendre-Gauss least square algorithm for linear problems.

m i n P A P - 2 = m i n P (29)

Where can be a product of matrix.

= A P - T A P - (30)

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 (31 ) and find where it's equal to zero..

( P ) = 2 - A P T P A P - = 2 J ( P ) T A P - (31)

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's the case when you want to estimate the external heat flux, ϕ0 in equation 25 . 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 (32 ).

P k + 1 = P k + [ ( J k ) T J k + μ k Ω k ] - 1 ( J k ) T - A ( P k ) (32)

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