Next Article in Journal
Raman 2D Peak Line Shape in Epigraphene on SiC
Next Article in Special Issue
Low Thermal Conductivity Materials and Very Low Heat Power: A Demanding Challenge in the Detection of Flaws in Multi-Layer Wooden Cultural Heritage Objects Solved by Pulse-Compression Thermography Technique
Previous Article in Journal
Virtual Environments and Augmented Reality Applied to Heritage Education. An Evaluative Study
Previous Article in Special Issue
Feature Point Matching Method Based on Consistent Edge Structures for Infrared and Visible Images
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Thermal Resistance Field in Layered Materials by Analytical Asymptotic Method

by
Marie-Marthe Groz
1,
Mohamed Bensalem
1,
Alain Sommier
2,
Emmanuelle Abisset-Chavanne
2,
Stéphane Chevalier
2,
Arsenii Chulkov
3,
Jean-Luc Battaglia
2,
Jean-Christophe Batsale
2 and
Christophe Pradere
2,*
1
Epsilon-Alcen, Esplanade des arts et métiers, 33405 Talence CEDEX, France
2
I2M UMR CNRS-ENSAM-UB 5295, Esplanade des arts et métiers, 33405 Talence CEDEX, France
3
National Research Tomsk Polytechnic University, Tomsk 634050, Russia
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(7), 2351; https://doi.org/10.3390/app10072351
Submission received: 29 January 2020 / Revised: 10 March 2020 / Accepted: 10 March 2020 / Published: 30 March 2020
(This article belongs to the Special Issue Infrared Imaging and NDT)

Abstract

:
In this paper, the problem of the quantitative characterization of thermal resistance fields in a multilayer sample is addressed by using the classical front face flash method as the thermal excitation and infrared thermography (IRT) as the monitoring sensor. In this challenging problem, the complete inverse processing of a multilayer analytical model is difficult due to the lack of sensitivity of some parameters (layer thickness, depth of thermal resistance, etc.) and the expansive computational iterative processing. For these reasons, the proposed strategy is to use a simple multilayer problem where only one resistive layer is estimated. Moreover, to simplify the inverse processing often based on iterative methods, an asymptotic development method is proposed here. Regarding the thermal signal reconstruction (TSR) methods, the drawback of these methods is the inability to be quantitative. To overcome this problem, the method incorporates a calibration process originating from the complete analytical quadrupole solution to the thermal problem. This analytical knowledge allows self-calibration of the asymptotic method. From this calibration, the quantitative thermal resistance field of a sample can be retrieved with a reasonable accuracy lower than 5%.

1. Introduction

In the domain of material and microelectronic device characterizations, many groups and studies have emerged with the goal to estimate both the thermal properties and resistance of these thin multilayers. One of the most common procedures is based on the 3 ω technique [1], which has been extensively used since it offers absolute measurements of the heat flux and the temperature and is also well suited for characterizations at low temperatures. At high temperatures, contactless photothermal methods, such as the thermoreflectance [2,3] and infrared radiometry [4,5,6] techniques, have been implemented and allow measuring the relative change in temperature and heat flux. Nevertheless, calibration is a complex task within those experimental configurations. In this community, very few studies are based on field characterizations. The notion of fields introduces the requirement to propose new methods for which the computation time is reasonable, passing from only one time vector, to find the inverse of a matrix. On the other hand, within the research groups involved in quantitative infrared thermography, many authors have developed a technique for qualitative in-depth defect detection by using flash thermography and thermal signal reconstruction (TSR) [7,8] or the lock-in technique [9]. Few authors, e.g., Tranta et al. [10] and Roche et al. [11], have linked the TSR method with quantitative estimation of the defect depth. The same kind of studies have been performed by Ibarra et al. [12] with a technique based on pulsed phase transformation (PPT) and thermographic signal reconstruction similar to [13]. As depicted in [14] one alternative to these qualitative methods is the use of parametrized 1D simulations with the quadrupole method as presented in the book of Maillet et al. [15]. Some works were performed by coupling infrared thermography (IRT) and the quadrupole approach, such as the one of Bendada et al. [16,17,18]. In fact, with these analytical methods, it becomes easy to calculate a multilayer sample with several thermal properties as well as thermal resistance. Then, by using nonlinear regression methods as proposed by Muller et al. [14], the parameters of the model can be adjusted by least square minimization between the model output and the experimental data. This approach is difficult when the initial values are poorly known and could be time consuming when addressing a large amount of data, such as temperature fields, with IRT. In the publication of Feuillet et al. [19], artificial delaminations were manufactured by introducing PTFE (Teflon) inlays in carbon-fiber-reinforced plastics, and they attempted to reduce the inversion by using singular value decomposition (SVD) [20,21,22,23]. Nevertheless, in this challenging problem, the complete inverse processing of a multilayer analytical model is difficult due to the lack of sensitivity of some parameters (layer thickness, depth of thermal resistance, etc.) and the expansive computational iterative processing. For these reasons, the proposed strategy is to use a simple multilayer problem where only one resistive layer is estimated. Moreover, to simplify the inverse processing often based on iterative methods, an asymptotic development method is proposed. The main novelty comes from the calibrating method between the analytical calculation of a theoretical base and the estimated parameters of the asymptotic method. Then with the knowledge of this calibration curve, the estimated parameters based on the measured one can be turn to quantitative ones.
Finally, in this work, a well-known flash setup is presented. A method is developed to model the heat transfer within the multilayer sample by using the quadrupole formulation. Then, in the inverse processing section, the autoregressive asymptotic development method is presented as well as the calibration procedure. Finally, the results and discussion section presents the characterization of a standard sample as a validation of the method.

2. Experimental Setup and Materials

The setup illustrated in Figure 1a and presented in Figure 1b is the classical front face flash system associated with an IR camera. The flash lamp comes from Uniblitz and has an energy E of 1600 J. Each lamp is synchronized with the IR camera by using an analogical TTL link to perform a pretrigger mode before the flash. The IR camera is an MCT long-wave ( λ = 9–11 μ m) FLIR SC7600 with a matrix sensor of 320 × 256 pixels and a pitch of 30 × 30 μ m 2 . With the used lenses and the sample distance, the resulting spatial resolution is approximately 300 × 300 μ m 2 per pixel.
For this study, a homemade reference sample was machined to calibrate the depth and thickness of the thermal resistance, represented by an air square hole, as illustrated in Figure 2. The sample is composed of polycarbonate layer of thickness e P C and the air thickness of the thermal resistance is noted e R t h .The geometrical and thermal properties are reported in Figure 2 and extracted from reference [24] for aluminum and [25] for the polycarbonate. In this simple case, only a bilayer sample was designed to demonstrate the capacity of the proposed method. However, the method can be extended to multilayer stack. It is important to note that the minimum tolerance guaranteed by the supplier is approximately +50 μ m of the desired thickness (see Figure 2).

3. Methods

3.1. Direct Problem

As illustrated in Figure 3, the direct problem of the heat transfer within a multilayer sample can be solved by using the quadrupole method [15]. It is important to note that according to the multilayer used for characterization, the complete stack can be overlaid onto a substrate with different thermal properties to also control the boundary condition imposed at the rear face of the sample. Usually, a foam bulk is used to mimic an adiabatic condition, whereas a metallic plate is used to mimic the imposed temperature.
Using the quadrupole formulation [15], the thermal response of any assembly (see Figure 3) can be generalized according to the Equation (1). The quadrupole method enables taking into account many layers as well as thermal resistances. A MATLAB homemade application (https://fr.mathworks.com/matlabcentral/fileexchange/74199-1dt-multilayer-thermalquadrupole-solver-and-builder) was developed to generate any assembly and to realize a sensitivity study for any parameters of the multilayer. The model is then expressed as:
θ e ( p ) Ψ e ( p ) h θ e ( p ) = i = 1 N M 1 R 1 R i 1 M i R i R N 1 M N θ s ( p ) h θ s ( p )
In the Equation (1), M i is the matrix relative to layer i and R i is the matrix associated with the thermal resistance between layers i and i + 1. They are both expressed in Equation (2). Moreover, θ e ( K ) and θ s ( K ) represent the Laplace transforms of the front and rear face temperatures, respectively. Ψ e ( p ) ( W . m 2 ) is the Laplace transform of the heat flux at the front face. Usually, the heat flux waveforms can be Dirac, heaviside or a pulse with a duration linked with the flash lamp discharge time. Finally, h ( W . m 2 . K 1 ) represents the convective heat losses.
M i = cosh p a i e i sinh p a i e i p a i λ i p a i λ i sinh p a i e i cosh p a i e i R i = 1 R t h i = e i λ i 0 1
In Equation (2), p ( s 1 ) is the Laplace variable, e i ( m ) is the thickness of layer i, a i (m 2 .s 1 ) is the thermal diffusivity of layer i, λ i (W.m 1 .K 1 ) is the thermal conductivity and R t h i ( m 2 . K . W 1 ) is the interfacial thermal resistance. In the case of the sample represented in Figure 2, the front face temperature can be solved for a wide range of layer 1 (i.e., the polycarbonate layer) thicknesses e P C = [ 100 , 1000 ] μ m (9 steps of 100 μ m) and corresponding thermal resistance thicknesses e R t h = [ 0 , 1000 ] μ m (10 steps of 100 μ m) with a thermal conductivity of λ a i r = 0.025 W.m 1 .K 1 to give the final generated thermal resistance R t h = [0,0.036] m 2 .K.W 1 with 10 steps of 0.004 m 2 .K.W 1 . Since the heat flux amplitude is unknown, the data are normalized according to any time step included in the semi-infinite behavior of the temperature relaxation inside the first cross layer according to the following formulation: T n ( t ) = T ( t t r e f ) T ( t r e f ) . To calculate the inverse Laplace transform of the temperature, the Stehfest algorithm [26] is used. The dimensionless temperatures as a function of time and for different polycarbonate layer thicknesses and thermal resistances are plotted in Figure 4. The thermal properties of the polycarbonate layer and the aluminum layer are given in Figure 2.

3.2. Autoregressive Asymptotic Method

The asymptotic method consists of the development of the amplitude and time delay between a reference temperature signal and a measured or generated one. The main assumption is that this pseudothermal contrast is sensitive only to the intensity variations A (due to the presence of the thermal resistance) and a delay τ ( s ) proportional to both the depth of the thermal resistance and its thickness. If both parameters are estimated only the amplitude will be used to quantitatively estimated the thermal resistance.
T m ( t ) = A T r e f ( t ) + τ d T r e f ( t ) d t
From the asymptotic Equation (3), and to avoid the influence of noisy data, an autoregressive inverse algorithm [27] is applied for a long time, as described by the Equation (4). The inverse processing is based on the Moore-Penrose inversion method [28,29]. The advantage of using this autoregressive method is the ability to estimate the parameters as a function of time to create a strong link with the maximum sensitivity corresponding to the early "separation" between the reference signal and the signal to process. With the above remarks under consideration, the inverse of the system can be written as follows
0 t 1 T m ( t ) d t 0 t i T m ( t ) d t 0 t N T m ( t ) d t = 0 t 1 T r e f ( t ) d t T r e f ( 0 , t 1 ) 0 t i T r e f ( t ) d t T r e f ( 0 , t i ) 0 t N T r e f ( t ) d t T r e f ( 0 , t N ) A 1 A i A N τ 1 τ i τ N Y = X β
where A ( ) is the amplitude of the asymptotic development, τ ( s ) is the time delay, t ( s ) is the time, and the indices m and r e f correspond to the measured temperature and the reference one. From Equation (4), the pseudoinverse X is obtained by using the singular value decomposition (SVD) methods [20,21,22,23]. The principle of SVD is to decompose a matrix I of size m by n into a product of three matrices, where U is an orthogonal matrix of size m by m, V is an orthogonal matrix of size n by n, and S is a “diagonal” matrix of size m by n. The diagonal elements of S are called singular values. They are positive and ranked from largest to smallest. If only the first p singular values are different from zero (p < n), the SVD can then be further reduced. Thus, the pseudoinversion of Equation (4) gives:
β = k = 1 p U k T Y s k
From the reference sample and the temperature response calculated and plotted in Figure 4, the autoregressive method is applied to all the data, with the R t h = 0 plot as the reference, corresponding to an assembly without any thermal resistance. The estimated parameters A ( ) is reported in Figure 5 as a function of the generated thermal resistance R t h . The time required to calculate all the time steps of the autoregressive method is approximately 6 s with a 2.3 GHz quad-core i5 MacBook Pro with 16 GB of RAM. The data can be fitted by using polynomial regression for the thermal resistance (Figure 5b). These results can be qualified by using a calibrating curve to link the nonquantitative estimation given by the asymptotic development method to the absolute value of he thermal resistance.

4. Results and Discussion

In this part, the results on the sample described in Section 2 are measured under the following parameters: (i) flash duration of 1 ms, (ii) camera frequency acquisition of 200 Hz, and (iii) 4000 images corresponding to a final time of 20 s and a pixel area of 300 × 300 μ m 2 for the spatial resolution. The obtained results are plotted in Figure 6:
From these measurements, 5 regions of interest (ROIs) appeared. For each ROI, the corresponding temperature evolution of one pixel can be extracted: 4 of them in the center of the thermal resistance square hole and the fifth in the center of the image. All these temperatures plot are normalized (baseline subtraction and reference division T n * ( t ) = T * ( t t r e f ) T * ( t r e f ) with T n * ( t ) = T n * ( t t 0 ) T 0 * ( t 0 ) , where t 0 ( s ) is the time of the flash and t r e f ( s ) any time during the -1/2 slope) and superposed onto the "healthy" analytical model (calculated for the bilayer of Figure 2) assuming perfect contact between all the layers. This normalization also avoid the non-uniformity of the global temperature coming from the non-homogeneity of the flash lamp. The results represented in Figure 7 show a -1/2 slope for a short time and the influence of the thermal resistance as well as the different depth effects as a function of time. In the Figure 7, it can be clearly observed that for the deepest ROI the temperature evolution quit the -1/2 slope further. For higher thermal resistance, the behavior is closer to adiabatic case with variable inflexion tendencies. Moreover, it is very interesting to note that the ROI extracted in the center of the image seems to be affected by a small thermal resistance (no exact superposition with the supposed healthy model can be achieved).
From this, we can estimate in Figure 8 the corresponding parameters. First, in Figure 8a, the obtained amplitude is plotted. Here, one can see that the estimated values are in good agreement with the ones estimated using the analytical solution represented in Figure 5a. Then, from the calibrated curves calculated and represented in Figure 5b, the absolute and quantitative fields of the amplitude value can be represented, as shown in Figure 8b. Here, the time required to calculate the parameters with the inverse process (4) is approximately 7 s with a 2.3 GHz quad-core i5 MacBook Pro with 16 GB of RAM. This computational time includes the inversion of a 51,012 × 4000 matrix representative of the number of pixels in space (218 × 234) and the number of time steps of the thermograms.
In this validation case, the geometry and thermal properties of the studied sample are well known (see Figure 2). Then, from the results in Figure 8, the agreement with the value is good, and the respective theoretical and estimated values are reported in Table 1.
From the data Table 1 one can see that the error is always lower than 5%, with the maximum uncertainty for ROI number 3 corresponding to a hole of 500 μ m. One can also note that a threshold has been performed before the identification in order to treat only the pertinent points belonging to the patches. Indeed, a certain number of points in the ROIs borders are not pertinent because of the strong gradient present here due to the 3D effect of the diffusion not taken into account in the identification process (see Figure 8). For example, for ROI3, 1500 points over 1900 have been used for the identification.

5. Conclusions

In this paper, a methodology for quantitative thermal resistance imaging estimation is proposed. This method is based on asymptotic development and quantitative calibration by using an analytical quadrupole formulation of the problem. The main advantages of this method are as follows: (i) the possibility to address any type of multilayer (high numbers of layers, etc.), (ii) a fast inverse process that can be adapted for online monitoring (inverse processing in less than 7 s), and (iii) the ability to overcome the knowledge of the initial excitation heat flux as well as the absolute measured temperature. The main drawbacks of the proposed methods are (i) the 1D character of the model and (ii) the sensitivity to only one thermal resistive layer. Nevertheless, the complete method applied on a heterogeneous bi-layer sample reveals a very good result with a complete estimation of the thermal resistance fields in fast time (complete calculation time lower than 7 s) and accuracy with maximum thermal resistance errors around 5%.

Author Contributions

Conceptualization, M.-M.G., M.B., A.S., E.A.-C., S.C., A.C., J.-L.B., J.-C.B. and C.P.; methodology, M.-M.G., M.B., A.S., E.A.-C., S.C., A.C., J.-L.B., J.-C.B. and C.P.; software, M.-M.G., M.B., A.S., E.A.-C., S.C., A.C., J.-L.B., J.-C.B. and C.P.; validation, M.-M.G., M.B., A.S., E.A.-C., S.C., A.C., J.-L.B., J.-C.B. and C.P.; formal analysis, M.-M.G., M.B., A.S., E.A.-C., S.C., A.C., J.-L.B., J.-C.B. and C.P.; investigation, M.-M.G., M.B., A.S., E.A.-C., S.C., A.C., J.-L.B., J.-C.B. and C.P.; resources, M.-M.G., M.B., A.S., E.A.-C., S.C., A.C., J.-L.B., J.-C.B. and C.P.; data curation, M.-M.G., M.B., A.S., E.A.-C., S.C., A.C., J.-L.B., J.-C.B. and C.P.; writing–original draft preparation, M.-M.G., M.B., A.S., E.A.-C., S.C., A.C., J.-L.B., J.-C.B. and C.P.; writing–review and editing, M.-M.G., M.B., A.S., E.A.-C., S.C., A.C., J.-L.B., J.-C.B. and C.P.; visualization, M.-M.G., M.B., A.S., E.A.-C., S.C., A.C., J.-L.B., J.-C.B. and C.P.; supervision, M.-M.G., M.B., A.S., E.A.-C., S.C., A.C., J.-L.B., J.-C.B. and C.P.; project administration, M.-M.G., M.B., A.S., E.A.-C., S.C., A.C., J.-L.B., J.-C.B. and C.P.; funding acquisition, M.-M.G., M.B., A.S., E.A.-C., S.C., A.C., J.-L.B., J.-C.B. and C.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cahill, D.G. Thermal conductivity measurement from 30 to 750 K: The 3ω method. Rev. Sci. Instrum. 1990, 61, 802–808. [Google Scholar] [CrossRef]
  2. Dilhaire, S.; Pernot, G.; Calbris, G.; Rampnoux, J.M.; Grauby, S. Heterodyne picosecond thermoreflectance applied to nanoscale thermal metrology. J. Appl. Phys. 2011, 110, 114314. [Google Scholar] [CrossRef]
  3. Baba, T.; Ishikawa, K.; Yagi, T.; Taketoshi, N. Measurements of thermophysical property of thin films by light pulse heating thermoreflectance methods. arXiv 2007, arXiv:0709.1845. [Google Scholar]
  4. Horny, N.; Chirtoc, M.; Fleming, A.; Hamaoui, G.; Ban, H. Kapitza thermal resistance studied by high-frequency photothermal radiometry. Appl. Phys. Lett. 2016, 109, 033103. [Google Scholar] [CrossRef]
  5. Depriester, M.; Hus, P.; Delenclos, S.; Sahraoui, A.H. New methodology for thermal parameter measurements in solids using photothermal radiometry. Rev. Sci. Instrum. 2005, 76, 074902. [Google Scholar] [CrossRef]
  6. Kusiak, A.; Battaglia, J.L.; Noe, P.; Sousa, V.; Fillot, F. Thermal conductivity of carbon doped GeTe thin films in amorphous and crystalline state measured by modulated photo thermal radiometry. J. Phys. Conf. Ser. 2016, 745, 032104. [Google Scholar] [CrossRef]
  7. Shepard, S.M. Flash thermography of aerospace composites. 2007. Available online: https://www.ndt.net/article/panndt2007/papers/132.pdf (accessed on 29 January 2020).
  8. Balageas, D.; Chapuis, B.; Deban, G.; Passilly, F. Improvement of the detection of defects by pulse thermography thanks to the TSR approach in the case of a smart composite repair patch. Quant. InfraRed Thermogr. J. 2010, 7, 167–187. [Google Scholar] [CrossRef]
  9. Karpen, W.; Wu, D.; Steegmuller, R.; Busse, G. Depth profiling of orientation in laminates with local lockin thermography. Proc. QIRT 1994, 94, 23–26. [Google Scholar]
  10. Oswald-Tranta, B.; Maier, A.; Schledjewski, R. Defect depth determination in a CFRP structure using TSR technique. In Proceedings of the 12th International Conference on Quantitative Infrared Thermography (QIRT 2014), Bordeaux, France, 7–11 July 2014. [Google Scholar]
  11. Roche, J.M.; Balageas, D. Imagerie quantitative TSR-RVB: De la détection de défauts au D-scan thermique. Available online: https://www.sft.asso.fr/Local/sft/dir/user-3775/documents/actes/Congres_2015/Communications/41753-fichier4.pdf (accessed on 29 January 2020).
  12. Ibarra-Castanedo, C.; Maldague, X. Pulsed phase thermography reviewed. Quant. Infrared Thermogr. J. 2004, 1, 47–70. [Google Scholar] [CrossRef]
  13. Shepard, S.M.; Lhota, J.R.; Rubadeux, B.A.; Wang, D.; Ahmed, T. Reconstruction and enhancement of active thermographic image sequences. Opt. Eng. 2003, 42, 1337–1343. [Google Scholar] [CrossRef]
  14. Müller, J.P.; Krankenhagen, R. Optimizing thermographic testing of thick GFRP plates by assessing the real energy absorbed within the material. Compos. Struct. 2019, 215, 60–68. [Google Scholar] [CrossRef]
  15. Maillet, D. Thermal Quadrupoles: Solving the Heat Equation Through Integral Transforms; John Wiley & Sons Inc.: Hoboken, NJ, USA, 2000. [Google Scholar]
  16. Bendada, A.; Erchiqui, F.; Lamontagne, M. Pulsed thermography in the evaluation of an aircraft composite using 3D thermal quadrupoles and mathematical perturbations. Inverse Probl. 2005, 21, 857. [Google Scholar] [CrossRef]
  17. Winfree, W.P.; Zalameda, J.N.; Howell, P.A.; Cramer, K.E. Simulation of thermographic responses of delaminations in composites with quadrupole method. Thermosense Ther. Infrared Appl. XXXVIII Int. Soc. Opt. Photonics 2016, 9861, 98610N. [Google Scholar]
  18. Bernegger, R.; Altenburg, S.J.; Roellig, M.; Maierhofer, C. Applicability of a 1D analytical model for pulse thermography of laterally heterogeneous semitransparent materials. Int. J. Thermophys. 2018, 39, 39. [Google Scholar] [CrossRef]
  19. Feuillet, V.; Ibos, L.; Fois, M.; Dumoulin, J.; Candau, Y. Defect detection and characterization in composite materials using square pulse thermography coupled with singular value decomposition analysis and thermal quadrupole modeling. NDT E Int. 2012, 51, 58–67. [Google Scholar] [CrossRef]
  20. Stewart, G.W. On the early history of the singular value decomposition. SIAM Rev. 1993, 35, 551–566. [Google Scholar] [CrossRef] [Green Version]
  21. Jordan, C. Sur la réduction des formes bilinéaires. Comptes Rendus De l’Académie Des Sciences 1874, 78, 614–617. [Google Scholar]
  22. Beltrami, E. Sulle funzioni bilineari. Giornale di Matematiche ad Uso degli Studenti Delle Universita 1873, 11, 98–106. [Google Scholar]
  23. Sylvester, J.J. Sur la reduction biorthogonale d’une forme lineo-linéaire a sa forme canonique. Comptes Rendus 1889, 108, 651–653. [Google Scholar]
  24. Dean, J.A. Lange’s Handbook of Chemistry; Mcgraw-Hill, Inc.: New York, NY, USA; London, UK, 1999. [Google Scholar]
  25. Touloukian, Y.S.; Powell, R.W.; Ho, C.Y.; Klemens, P.G. Thermophysical Properties of Matter-the tprc Data Series. Volume 1. Thermal Conductivity-Metallic Elements and Alloys; Thermophysical and Electronic Properties Information Analysis Center: Lafayette, IN, USA, 1970. [Google Scholar]
  26. Stehfest, H. Algorithm 368: Numerical inversion of Laplace transforms [D5]. Commun. ACM 1970, 13, 47–49. [Google Scholar] [CrossRef]
  27. Zheng, W.X. A least-squares based method for autoregressive signals in the presence of noise. IEEE Trans. Circuits Syst. II Analog Digit. Signal Process. 1999, 46, 81–85. [Google Scholar] [CrossRef]
  28. Moore, E.H. On the reciprocal of the general algebraic matrix. Bull. Am. Math. Soc. 1920, 26, 394–395. [Google Scholar]
  29. Penrose, R. A generalized inverse for matrices. In Mathematical Proceedings of the Cambridge Philosophical Society; Cambridge University Press: Cambridge, UK, 1955; Volume 51, pp. 406–413. [Google Scholar]
Figure 1. (a) Schema of the experimental setup and (b) photograph of the bench.
Figure 1. (a) Schema of the experimental setup and (b) photograph of the bench.
Applsci 10 02351 g001
Figure 2. (a) Photograph of the polycarbonate layer and (b) schema of the reference sample used in this study.
Figure 2. (a) Photograph of the polycarbonate layer and (b) schema of the reference sample used in this study.
Applsci 10 02351 g002aApplsci 10 02351 g002b
Figure 3. Schema of multilayer with thermal resistance and a boundary condition.
Figure 3. Schema of multilayer with thermal resistance and a boundary condition.
Applsci 10 02351 g003
Figure 4. Loglog dimensionless temperature response of the front face to a pulse of 1 ms duration and under adiabatic boundary conditions as a function of the thermal resistance depth and thickness of the reference sample composed of polycarbonate and aluminum.
Figure 4. Loglog dimensionless temperature response of the front face to a pulse of 1 ms duration and under adiabatic boundary conditions as a function of the thermal resistance depth and thickness of the reference sample composed of polycarbonate and aluminum.
Applsci 10 02351 g004
Figure 5. Results of the autoregressive inverse method based on the simulated analytical temperature response from Figure 4: (a) estimated amplitude A as a function of time, the solid lines corresponds to the different R t h of the Figure 4 and (b), calibration curve obtained from the extraction of the maximum amplitude A as a function of the thermal resistance and amplitude at any time step of 19.7 s.
Figure 5. Results of the autoregressive inverse method based on the simulated analytical temperature response from Figure 4: (a) estimated amplitude A as a function of time, the solid lines corresponds to the different R t h of the Figure 4 and (b), calibration curve obtained from the extraction of the maximum amplitude A as a function of the thermal resistance and amplitude at any time step of 19.7 s.
Applsci 10 02351 g005
Figure 6. Temperature of the sample at (a) t = 0.87 s, (b) t = 1.12 s, (c) t = 1.67 s and (d) t = 5.12 s.
Figure 6. Temperature of the sample at (a) t = 0.87 s, (b) t = 1.12 s, (c) t = 1.67 s and (d) t = 5.12 s.
Applsci 10 02351 g006
Figure 7. Loglog dimensionless temperature response of the front face to a pulse of 1 ms duration and adiabatic boundary conditions as a function of the thermal resistance for a sample composed of one layer of polycarbonate, one layer of thermal resistance and one aluminum layer.
Figure 7. Loglog dimensionless temperature response of the front face to a pulse of 1 ms duration and adiabatic boundary conditions as a function of the thermal resistance for a sample composed of one layer of polycarbonate, one layer of thermal resistance and one aluminum layer.
Applsci 10 02351 g007
Figure 8. Fields of: (a) estimated amplitude and (b), calculated thermal resistance based on the calibration curve in Figure 5a.
Figure 8. Fields of: (a) estimated amplitude and (b), calculated thermal resistance based on the calibration curve in Figure 5a.
Applsci 10 02351 g008
Table 1. Comparison between the theoretical and estimated values.
Table 1. Comparison between the theoretical and estimated values.
e Rth λ air Rth theo Rth est Error ***
( μ m)(W.m 1 .K 1 )(m 2 .K.W 1 )(m 2 .K.W 1 )(%)
* ROI 1[200–250] **0.025[0.008–0.01]0.00873.4
* ROI 2[300–350] **0.025[0.012–0.014]0.01263.1
* ROI 3[500–550] **0.025[0.02–0.022]0.02215.2
* ROI 4[800–850] **0.025[0.032–0.034]0.03280.85
*, The ROI order can be found in the schema in Figure 2 or the results in Figure 8; **, Systematic maximum machining uncertainty of +50 μ m given by the supplier; see Section 2; ***, Error 100 × 1 R t h e s t R t h t h e o ¯ calculated with the average value R t h t h e o ¯ of the R t h t h e o gap.

Share and Cite

MDPI and ACS Style

Groz, M.-M.; Bensalem, M.; Sommier, A.; Abisset-Chavanne, E.; Chevalier, S.; Chulkov, A.; Battaglia, J.-L.; Batsale, J.-C.; Pradere, C. Estimation of Thermal Resistance Field in Layered Materials by Analytical Asymptotic Method. Appl. Sci. 2020, 10, 2351. https://doi.org/10.3390/app10072351

AMA Style

Groz M-M, Bensalem M, Sommier A, Abisset-Chavanne E, Chevalier S, Chulkov A, Battaglia J-L, Batsale J-C, Pradere C. Estimation of Thermal Resistance Field in Layered Materials by Analytical Asymptotic Method. Applied Sciences. 2020; 10(7):2351. https://doi.org/10.3390/app10072351

Chicago/Turabian Style

Groz, Marie-Marthe, Mohamed Bensalem, Alain Sommier, Emmanuelle Abisset-Chavanne, Stéphane Chevalier, Arsenii Chulkov, Jean-Luc Battaglia, Jean-Christophe Batsale, and Christophe Pradere. 2020. "Estimation of Thermal Resistance Field in Layered Materials by Analytical Asymptotic Method" Applied Sciences 10, no. 7: 2351. https://doi.org/10.3390/app10072351

APA Style

Groz, M. -M., Bensalem, M., Sommier, A., Abisset-Chavanne, E., Chevalier, S., Chulkov, A., Battaglia, J. -L., Batsale, J. -C., & Pradere, C. (2020). Estimation of Thermal Resistance Field in Layered Materials by Analytical Asymptotic Method. Applied Sciences, 10(7), 2351. https://doi.org/10.3390/app10072351

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop