Abstract
The flash method is a commonly used technique for estimating the thermal diffusivity of solid materials. In this work, we develop a new efficient formula for the thermal diffusivity calculation to process the flash method data in the case of heat losses at the rear and front faces. First, the transient heat equation is analytically solved using the Green functions under time-dependent convection boundary conditions. The property of the descending part of the analytical solution, containing series which converges very quickly for long times, allows to express the thermal diffusivity as a function of the slope of the natural logarithm of the temperature rise and the Biot numbers related to heat losses at the rear and front faces. These Biot numbers are calculated from the ratio between the integrals of the complete temperature trends of the two faces. For this purpose, we developed a novel way to evaluate these integrals using the analytical solutions in Laplace domain. To verify the accuracy of these formulas, we simulated experimental data adding a Gaussian noise to the theoretical temperature results of the flash method data.
Similar content being viewed by others
Abbreviations
- a:
-
Thermal diffusivity (m2·s−1)
- Cp :
-
Specific heat (J·k−1·kg−1)
- h1 :
-
Heat transfer coefficient at front face (W·m−2·k−1)
- h2 :
-
Heat transfer coefficient at rear face (W·m−2·k−1)
- λ:
-
Thermal conductivity (W·m−1·k−1)
- e:
-
Length of the sample (m)
- Q0 :
-
Amount of heat absorbed through the front surface of the sample per unit area (J·m−2)
- Bi1 :
-
\(\frac{{eh_{1} }}{\lambda }\) the Biot number at front face
- Bi2 :
-
\(\frac{{eh_{2} }}{\lambda }\) the Biot number at rear face
- ρ:
-
Density of the sample (kg·m−3)
- Tma :
-
\(\frac{{Q_{0} }}{{e \cdot \rho c_{p} }}\) adiabatic limit temperature (°C)
- \(\tilde{T}\) :
-
Simulated front-surface temperature rise data (°C)
- T:
-
Simulated rear-surface temperature rise data (°C)
- t:
-
Time (s)
- td :
-
Pulse duration (s)
- ti :
-
ith sampled time (s)
- \(\tilde{t}_{N}\) :
-
Final sampled time at front face (s)
- t Z :
-
Final sampled time at rear face (s)
- x:
-
Spatial coordinate (m)
- r:
-
Normally distributed random number (°C)
- σ:
-
Standard deviation
- N:
-
Number of temperature samples recorded at front face (excluding initial temperature)
- Z:
-
Number of temperature samples recorded at rear face (excluding initial temperature)
- \(\varepsilon\) :
-
Relative error (%)
- T(0,t):
-
The theoretical rise temperature at front face (°C)
- T(e,t):
-
The theoretical rise temperature at rear face (°C)
- S:
-
\(\int\limits_{0}^{ + \infty } {T\left( {x,t} \right)} dt\) (°C·s)
References
W.J. Parker, R.J. Jenkins, C.P. Butler, G.L. Abbott, J. Appl. Phys. 32, 1679–1684 (1961)
L. Vozar, W. Hohenauer, High Temp. High Press. 35, 253–264 (2003)
J. Gembarovic, R.E. Taylor, Rev. Sci. Instrum. 65, 3535–3539 (1994)
E.J. Carr, Chem. Eng. Sci. 199, 546–551 (2019)
R.C. Heckman, J. Appl. Phys. 44, 1455–1460 (1973)
J.A. Cape, G.W. Lehman, J. Appl. Phys. 34, 1909–1913 (1963)
Y. Tao, L. Yang, Q. Zhong, Z. Xu, C. Luo, Rev. Sci. Instrum. 86, 124902 (2015)
T. Azumi, Y. Takahashi, Rev. Sci. Instrum. 52, 1411–1413 (1981)
R.E. Taylor, J.A. Cape, Appl. Phys. Lett. 5, 212–213 (1964)
M. Potenza, A. Cataldo, G. Bovesecchi, S. Corasaniti, P. Coppa, S. Bellucci, AIP Adv. (2017). https://doi.org/10.1063/1.4995513
R.D. Cowan, J. Appl. Phys. 34, 926–927 (1963)
K. Beedham, I.P. Dalrymple, The measurement of thermal diffusivity by the flash method. An investigation into errors arising from the boundary conditions. Revue Internationale des Hautes Temp. et des Refractaires 7, 278–283 (1970)
S. Bellucci, G. Bovesecchi, A. Cataldo, P. Coppa, S. Corasaniti, M. Potenza, Materials 12, 696 (2019). https://doi.org/10.3390/ma12050696
B. Czél, K.A. Woodbury, J. Woolley, G. Gróf, Int. J. Thermophys. 34, 1993–2008 (2013)
L. Chen, A.M. Limarga, D.R. Clarke, Comput. Mater. Sci. 50, 77–82 (2010)
H.M. James, J. Appl. Phys. 51, 4666–4672 (1980)
Y.H. Zhao, Z.K. Wu, S.L. Bai, Int. J. Heat Mass Transf. 101, 470–475 (2016)
E.J. Carr, C.J. Wood, Int. J. Heat Mass Transf. 144, 118609 (2019)
L. Duquesne, C. Lorrette, C. Pradere, G.L. Vignoles, J.C. Batsale, Quant. InfraRed Thermogr. J. 13, 182–194 (2016)
R.E. Taylor, L.M. Clark, High Temp. High Press. 6, 65 (1974)
Heckman, R. C. Error analysis of the flash thermal diffusivity technique, Thermal conductivity, vol 14 (Springer, Boston, 1976), pp. 491–498
A. Degiovanni, Correction de longueur d’impulsion pour la mesure de la diffusivité thermique par méthode flash. Int. J. Heat Mass Transf. 30, 2199–2200 (1987)
L.M. Clark III, R.E. Taylor, J. Appl. Phys. 46, 714–719 (1975)
A. Degiovanni, M. Laurent, Revue de physique Appliquée 21, 229–237 (1986)
D.L. Balageas, Revue de Physique Appliquée 17, 227–237 (1982)
L. Vozár, J. Gembarovič, V. Majerník, Int. J. Heat Mass Transf. 34, 1316–1318 (1991)
L. Vozár, J. Gembarovič, V. Majerník, High Temp. High Press. 23, 397–402 (1991)
G.F. Bocchini, G. Bovesecchi, P. Coppa, S. Corasaniti, R. Montanari, A. Varone, Int. J. Thermophys. 37, 1–14 (2016). https://doi.org/10.1007/s10765-016-2050-4
J.A. McKay, J.T. Schriempf, J. Appl. Phys. 47, 1668–1671 (1976)
T. Yamane, S. Katayama, M. Todoki, Int. J. Thermophys. 18, 269–290 (1997)
Y. Takahashi, K. Yamamoto, T. Ohsato, T. Terai, Proc. 9th Japan Symp. Thermophysical Properties, 1988, pp. 175–178
J. Gembarovič, L. Vozár, V. Majerník, Int. J. Heat Mass Transf. 33, 1563–1565 (1990)
Y. Chihab, S. Raefat, M. Garoum, N. Laaroussi, R. Bouferra, 6th international renewable and sustainable energy conference, IRSEC, 2018. https://doi.org/10.1109/IRSEC.2018.8702890
Y. Chihab, S. Raefat, M. Garoum, N. Laaroussi, R. Bouferra, in IOP conference series: materials science and engineering, Vol 446, No 1 (IOP Publishing, 2018), p. 012008
Y. Chihab, M. Garoum, N. Laaroussi, Int. J. Innov. Technol. Explor. Eng. (2019). https://doi.org/10.35940/ijitee.k1801.0981119
N. Laaroussi, G. Lauriat, M. Garoum, A. Cherki, Y. Jannot, Constr. Build. Mater. 70, 351–361 (2014)
M. Boussaid, Y. Boulaadjoul, M. Hedibel, T. Loulou, Heat Transf. Eng. 35, 1201–1208 (2014)
R.P. Brent, Chapter 4: An algorithm with guaranteed convergence for finding a zero of a function, Algorithms for minimization without derivatives (Prentice-Hall, Englewood Cliffs, 1973). ISBN 0-13-022335-2
R.L. Burden, J.D. Faires, 2.1 The bisection algorithm, Numerical analysis, 3rd edn. (PWS Publishers, Boston, 1985). ISBN 0-87150-857-5
M.B. Allen, E.L. Isaacson, Numerical analysis for applied science (Wiley, Hoboken, 1998), pp. 188–195. ISBN 978-0-471-55266-6
Weisstein EW. “Inverse quadratic interpolation.” From MathWorld—a wolfram web resource. http://mathworld.wolfram.com/InverseQuadraticInterpolation.html
L.I. Turchak, P.V. Plotnikov, Fundamentals of numerical methods (Nauka, Moscow, 1987)
K.E. Atkinson, An introduction to numerical analysis, 2nd edn. (Wiley, Hoboken, 1989). ISBN 0-471-50023-2
A. Degiovanni, Diffusivité et méthode flash. Revue générale de thermique 185, 420–442 (1977)
J.A. Nelder, R. Mead, Comput. J. 7, 308–331 (1965)
Y. Jannot, A. Degiovanni. Mesure des propriétés thermiques des matériaux. ISTE Editions, 1st edn, April 2018
A. Tadeu et al., Energy Build. 68, 611–619 (2014)
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendix A. Analytical Expression of the Temperature Rise
Appendix A. Analytical Expression of the Temperature Rise
The system (A1) can be reduced to a problem with homogeneous boundary conditions.
We suggest that
where K(x,t) is a function satisfying the boundary conditions (Eq. A3)
We assume that K(x, t) is written as follows:
The coefficients A1 and A2 are calculated using the boundary conditions (Eq. A3):for x = 0,
for x = e,
K(x,t) can be written as follows:
Then, the function F(x,t) must be the solution of the problem:
And
With
The first step in resolving the system (Eqs. A8 and A9) is to consider the PDE (Eq. A8) homogenous:
Using the method of separating the variables, we obtain
The eigenvalues fn are the solutions of the following transcendent equation:
For this problem, the normalized eigenfunctions are
The orthogonality properties of the normalized eigenfunctions An (x):
We achieve
We express the unknown solution F(x, t) as a generalized eigenfunctions series of Fourier with time-dependent coefficients:
After replacement Eq. A17 in Eq. A8 we obtain
We have
\(\Rightarrow\)
Replacement of Eq. A19 in Eq. A18:
Then we also develop p(x, t) as a generalized eigenfunctions series of Fourier with a time-dependent coefficient.
We multiply Eq. A21 by An(x):
We integrate Eq. A22 from 0 to e:
Assuming the uniform convergence:
The orthogonality relation of the normalized eigenfunctions leads to:
Replacement of Eq. A21 in Eq. A18 gives
This leads to
The time-dependent ODE solution mentioned above (Eq. A28) is presented as follows:
In the equation (Eq. A29), τ is the auxiliary time and G1(t, μ) is the first-order Green function defined as follows:
We use the initial condition to obtain the initial coefficient Bn(0):
Assuming uniform convergence of the series and the orthogonality property of the eigenfunctions An(x), we obtain
The equation (Eq. A29) becomes
Finally, the solution of F(x, t) can be written as follows:
The temperature T(x,t) can be written as follows:
The temperature rise on the front face(x = 0) is written as
The rise temperature on the rear face (x = e) is written in the following form:
Rights and permissions
About this article
Cite this article
Chihab, Y., Garoum, M. & Laaroussi, N. A New Efficient Formula for the Thermal Diffusivity Estimation from the Flash Method Taking into Account Heat Losses in Rear and Front Faces. Int J Thermophys 41, 118 (2020). https://doi.org/10.1007/s10765-020-02704-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10765-020-02704-w