Abstract
Infrared thermography (IRT) and heat source reconstruction (HSR) were used in the study to measure two thermophysical properties of superelastic nickel-titanium (NiTi) shape memory alloy (SMA) wires, namely thermal conductivity and specific heat capacity. Since the values are potentially phase dependent, the identification was carried out at ambient temperature with and without mechanical loading, i.e., in the stress-induced martensite state and in the austenite state, respectively. A uniaxial testing machine was used to apply constant deformation while allowing temperature measurements by IRT during Joule heating and natural return the thermal equilibrium. The data processing by HSR, including preliminary filtering operations, was preliminarily evaluated from temperatures provided by a model (with added noise). It gave an error of ± 0.5% for the thermal conductivity and ± 1% for the specific heat capacity. The experimental analysis showed that the thermal conductivity of stress-induced martensite is 26% higher than that of austenite at the same temperature (here the ambient temperature). Regarding the specific heat capacity, that of stress-induced martensite is 4.7% lower than that of austenite. Explanations for these differences were proposed from solid-state physics theory. The measured values were also compared to data collected in the literature (obtained at zero stress at different temperatures).
Similar content being viewed by others
Data Availability
Data will be made available on request.
Abbreviations
- IRT:
-
Infrared Thermography
- HSR:
-
Heat Source Reconstruction
- SMA:
-
Shape Memory Alloy
- NiTi:
-
Nickel-Titanium
- DIC:
-
Digital Image Correlation
- TS1:
-
Heating state during the heating stage
- SS1:
-
Steady state during the heating stage
- TS2:
-
Heating state during the cooling stage
- SS2:
-
Steady state during the cooling stage
- \({A}_{\text{s}}\) :
-
Austenite-start temperature [\(^\circ \text{C}\)]
- \({A}_{\text{f}}\) :
-
Austenite-finish temperature [\(^\circ \text{C}\)]
- \({M}_{\text{s}}\) :
-
Martensite-start temperature [\(^\circ \text{C}\)]
- \({M}_{\text{f}}\) :
-
Martensite-finish temperature [\(^\circ \text{C}\)]
- \(Bi\) :
-
Biot number [-]
- \(h\) :
-
Thermal convection coefficient [\({\text{W}}\cdot {\text{m}}^{-2}\cdot {{\text{K}}}^{-1}\)]
- \(r\) :
-
Wire radius [\(\text{m}\)]
- \(\lambda\) :
-
Thermal conductivity [\({\text{W}}\cdot {\text{m}}^{-1}\cdot {{\text{K}}}^{-1}\)]
- \(C\) :
-
Specific heat capacity [\({\text{J}}\cdot {\text{kg}}^{-1}\cdot {{\text{K}}}^{-1}\)]
- \(T\left(z,t\right)\) :
-
1D spatio-temporal temperature distribution [\(^\circ \text{C}\)]
- \(\epsilon\) :
-
Thermal emissivity [-]
- \(\rho\) :
-
Density of the material [\(\text{kg}\cdot {\text{m}}^{-3}\)]
- \(\tau\) :
-
Exponential decay of the transient regimes [\(\text{s}\)]
- \({T}_{\text{ref}}\) :
-
Temperature of ambient air [\(^\circ \text{C}\)]
- \({\sigma }_{\text{SB}}\) :
-
Stefan-Boltzmann constant (\(5.67\times {10}^{-8}\text{ W}\cdot {\text{m}}^{-2}\cdot {{\text{K}}}^{-4}\))
- \({S}_{\text{el}}\) :
-
Heat source produced by the material [\({\text{W}}\cdot {\text{m}}^{-3}\)]
- \({R}_{\text{aust}}\) :
-
Radius of the wire in austenite phase [\(\text{m}\)]
- \({R}_{\text{mart}}\) :
-
Radius of the wire in stress-induced martensite phase [\(\text{m}\)]
- \({\rho }_{\text{el}}\) :
-
Electrical resistivity [\(\Omega \cdot \text{m}\)]
- \(U\) :
-
Voltage [\(\text{V}\)]
- \({S}_{0}\) :
-
Initial cross section of the wire [\({\text{m}}^{2}\)]
- \(I\) :
-
Electrical current [\(\text{A}\)]
- \({L}_{0}\) :
-
Initial gage length [\(\text{m}\)]
- \(\varepsilon\) :
-
Longitudinal strain [-]
- \(\mu\) :
-
Mean value
- \(\sigma\) :
-
Standard deviation
- \({\phi }_{\text{hs}}\) :
-
Storage term in the heat diffusion equation [\({\text{W}}\cdot {\text{m}}^{-3}\)]
- \({\phi }_{\text{cd}}\) :
-
Conduction term in the heat diffusion equation [\({\text{W}}\cdot {\text{m}}^{-3}\)]
- \({\phi }_{\text{cv}}\) :
-
Convection term in the heat diffusion equation [\({\text{W}}\cdot {\text{m}}^{-3}\)]
- \({\phi }_{\text{rd}}\) :
-
Radiation term in the heat diffusion equation [\({\text{W}}\cdot {\text{m}}^{-3}\)]
- \({S}_{\text{th}}\) :
-
Reconstructed heat source [\({\text{W}}\cdot {\text{m}}^{-3}\)]
- \({T}_{\text{calc}}\) :
-
Temperature calculated from a model [\({\text{K}}\)]
- \({T}_{\text{exp}}\) :
-
Experimental temperature [\({\text{K}}\)]
- \({\lambda }_{\text{L}}\) :
-
Lattice contribution to the thermal conductivity [\({{\text{W}}}\cdot {\text{m}}^{-1}\cdot {{\text{K}}}^{-1}\)]
- \({\lambda }_{\text{e}}\) :
-
Electronic contribution to the thermal conductivity [\({\text{W}}\cdot {\text{m}}^{-1}\cdot {{\text{K}}}^{-1}\)]
- \(L\) :
-
Sommerfeld value of the Lorentz number (\(2.45\times {10}^{-8}\text{ W}\cdot\Omega \cdot {{\text{K}}}^{-2}\))
- \({k}_{\text{b}}\) :
-
Boltzmann’s constant (\(1.380649\times {10}^{-23}{\text{ m}}^{2}\cdot \text{kg}\cdot {\text{s}}^{-2}\cdot {{\text{K}}}^{-1}\))
- \(\hslash\) :
-
Planck’s constant (1\(.054571817\times {10}^{-34}\text{ J}\cdot \text{s}\))
- \(c\) :
-
Phonon velocity [\(\text{m}\cdot {\text{s}}^{-1}\)]
- \(\Theta\) :
-
Debye temperature [\({\text{K}}\)]
- \(l\) :
-
Photon mean free path [\(\text{m}\)]
- \({C}_{\text{v}}\) :
-
Specific heat capacity at constant volume or deformation [\({\text{J}}\cdot {\text{kg}}^{-1}\cdot {{\text{K}}}^{-1}\)]
- \({C}_{\text{p}}\) :
-
Specific heat capacity at constant pressure or stress [\({\text{J}}\cdot {\text{kg}}^{-1}\cdot {{\text{K}}}^{-1}\)]
- \(R\) :
-
Gas constant \(\left(8.314\text{ J}\cdot {{\text{K}}}^{-1}\cdot {\text{mol}}^{-1}\right)\)
- \(m\) :
-
Number of atoms in crystal cell
- \({\theta }_{\text{ref}}\) :
-
Temperature change of ambient air with respect to the test beginning \([^\circ \text{C}]\)
References
J. Shaw, S. Kyriakides, J. Mech. Phys. Solids 43, 1243 (1995)
E. Alarcon, L. Heller, Mater. Des. 199, 109406 (2021)
K. Otsuka, X. Ren, Prog. Mater. Sci. 50, 511 (2005)
O. Molnárová, O. Tyc, L. Heller, H. Seiner, P. Šittner, Acta Mater. 218, 1 (2021)
P. Šittner, P. Sedlák, H. Seiner, P. Sedmák, J. Pilch, R. Delville, L. Heller, L. Kadeřávek, Prog. Mater. Sci. 98, 249 (2018)
L. Heller, H. Seiner, P. Šittner, P. Sedlák, O. Tyc, L. Kadeřávek, Int. J. Plast. 111, 53 (2018)
S.-M. Lee, D.G. Cahill, T.H. Allen, Phys. Rev. B 52, 253 (1995)
P. Cielo, L.A. Utracki, M. Lamontagne, Can. J. Phys. 64, 1172 (1986)
A.I. Kovalev, A.V. Logunov, N.V. Petrushin, L.S. Egorova, J. Eng. Phys. 32, 183 (1977)
J. F. Goff, J. Appl. Phys. 39, 2208 (1968)
M. Chirtoc, J. Gibkes, R. Wernhardt, J. Pelzl, A. Wieck, Rev. Sci. Instrum. 79, 093703 (2008)
J. Gibkes, B.K. Bein, J. Pelzl, A. Wieck, J. Frenzel, G. Eggeler, I. Delgadillo-Holtfort, M. Chirtoc, Acta Phys. Pol. A 114, 10 (2008)
C. Zanotti, P. Giuliani, G. Riva, A. Tuissi, A. Chrysanthou, J. Alloys Compd. 473, 231 (2009)
C. Zanotti, P. Giuliani, P. Bassani, Z. Zhang, A. Chrysanthou, Intermetallics 18, 14 (2010)
C. Zanotti, P. Giuliani, A. Chrysanthou, Intermetallics 24, 106 (2012)
V. Delobelle, D. Favier, H. Louche, J. Mater. Eng. Perform. 22, 1688 (2013)
V. Delobelle, D. Favier, H. Louche, N. Connesson, Exp. Mech. 55, 711 (2015)
M. Rohde, A. Schüssler, Sens. Actuators A Phys. 61, 463 (1997)
M.G. Faulkner, J. J. Amalraj, A. Bhattacharyya, Smart Mater. Struct. 9, 632 (2000)
Y. Terada, K. Ohkubo, K. Nakagawa, T. Mohri, T. Suzuki, Intermetallics 3, 347 (1995)
A. Jain, K.E. Goodson, J. Heat Transfer 130, 1 (2008)
B.D. Ingale, W.C. Wei, P.C. Chang, Y.K. Kuo, S.K. Wu, J. Appl. Phys. 110, 113721 (2011)
M. Kaya, A. Buğutekin, N. Orhan, Int. J. Thermophys. 32, 665 (2011)
M. Zhelnin, A. Iziumova, A. Vshivkov, O. Plekhov, Quant. Infrared Thermogr. J. 16, 74 (2019)
A. Jury, X. Balandraud, L. Heller, P. Šittner, M. Karlik, Exp. Mech. 61, 349 (2021)
L.R. Ingersoll, J. Opt. Soc. Am. 9, 495 (1924)
T. Ozawa, Jpn. J. Appl. Phys. 6, 1455 (1967)
J. H. Blackwell, 25 (1954).
E. García, I. Amaya, R. Correa, Appl. Therm. Eng. 129, 587 (2018)
W.J. Parker, R.J. Jenkins, C.P. Butler, G.L. Abbott, J. Appl. Phys. 32, 1679 (1961)
Á. Cifuentes, A. Mendioroz, A. Salazar, Int. J. Therm. Sci. 121, 305 (2017)
N. Tao, X.L. Li, J. G. Sun, Rev. Sci. Instrum. 88, 064903 (2017)
C. Boué, S. Holé, Infrared Phys. Technol. 55, 376 (2012)
A. Lara-Bernal, E. Marín, A. Calderón, O. Delgado-Vasallo, Int. J. Thermophys. 33, 1966 (2012)
M. Pilarczyk, B. Węglowski, Appl. Therm. Eng. 150, 1224 (2019)
D.G. Cahill, Rev. Sci. Instrum. 61, 802 (1990)
P. Cielo, L.A. Utracki, M. Lamontagne, Can. J. Phys. 64, 1172 (1985)
F. Cernuschi, P.G. Bison, A. Figari, S. Marinetti, E. Grinzato, Int. J. Thermophys. 25, 439 (2004)
K. Martínez, E. Marín, C. Glorieux, A. Lara-Bernal, A. Calderón, G. Peña Rodríguez, R. Ivanov, Int. J. Therm. Sci. 98, 202 (2015)
M.J. Assael, K.D. Antoniadis, I.N. Metaxa, S.K. Mylona, J. -A.M. Assael, J. Wu, M. Hu, Int. J. Thermophys. 36, 3083 (2015)
A.B. Donaldson, R.E. Taylor, J. Appl. Phys. 46, 4584 (1975)
R.A. Perkins, H.M. Roder, C.A.N. Decastro, J. Res. Natl. Inst. Stand. Technol. 96, 247 (1991)
S. Alvarado, E. Marín, A.G. Juárez, A. Calderón, R. Ivanov, Eur. J. Phys. 33, 897 (2012)
C. Boué, D. Fournier, Quant. Infrared Thermogr. J. 6, 175 (2009)
C. Ravey, C. Pradere, N. Regnier, J. -C. Batsale, Quant. Infrared Thermogr. J. 9, 79 (2012)
F. C. Malheiros, A. A. A. Figueiredo, L. H. da S. Ignacio, and H. C. Fernandes, Appl. Therm. Eng. 157, 113696 (2019).
J. Ospina-Borras, H. Benitez-Restrepo, J. Florez-Ospina, Appl. Sci. 7, 1253 (2017)
M. Zhelnin, A. Iziumova, and O. Plekhov, in AIP Conf. Proc. (2016), p. 040095.
C. Rodiet, B. Remy, A. Degiovanni, Int. J. Therm. Sci. 79, 67 (2014)
S.E. Gustafsson, Rev. Sci. Instrum. 62, 797 (1991)
I. Philippi, J. C. Batsale, D. Maillet, A. Degiovanni, Rev. Sci. Instrum. 66, 182 (1995)
V. Vavilov, D. Burleigh, V. Shiryaev, Quant. Infrared Thermogr. J. 4, 187 (2007)
H.A. Berman, E.D. West, A.G. Rozner, J. Appl. Phys. 38, 4473 (1967)
Z. Jiang, C.T. Imrie, J. M. Hutchinson, Thermochim. Acta 387, 75 (2002)
D. Delpueyo, X. Balandraud, M. Grédiac, Infrared Phys. Technol. 60, 312 (2013)
Jack P. Holman, Heat Transfer (2009).
P. Šittner, P. Sedlák, M. Landa, V. Novák, P. Lukáš, Mater. Sci. Eng. A 438–440, 579 (2006)
T.M. Tritt, Thermal Conductivity, 1st edn. (Springer, US, New York, 2004)
J. Vancea, G. Reiss, H. Hoffmann, Phys. Rev. B 35, 6435 (1987)
D. Bradley, J. Acoust. Soc. Am. 37, 700 (1965)
J. F. Goff, J. Appl. Phys. 35, 2919 (1964)
M. K. Stanford, Nasa Tech. Inf. Serv. 12 (2012).
C. Zanotti, P. Giuliani, A. Tuissi, S. Arnaboldi, and R. Casati, in ESOMAT 2009 - 8th Eur. Symp. Martensitic Transform. (EDP Sciences, Les Ulis, France, 2009), pp. 1–7.
I. Delgadillo-Holtfort, M. Chirtoc, J. Gibkes, B.K. Bein, J. Pelzl, Int. J. Thermophys. 35, 2308 (2014)
C. Kittel, R.W. Hill, Phys. Today 61, 1 (1980)
M. Fukuhara, M. Yagi, A. Matsuo, Phys. Rev. B 65, 224210 (2002)
R.J. Wasilewski, S.R. Butler, J. E. Hanlon, Met. Sci. J. 1, 104 (1967)
J. F. Smith, R. Lück, Q. Jiang, B. Predel, J. Phase Equilibria 14, 494 (1993)
M.J. Hadamard, Princet. Univ. Bull. XIII, 1 (1902)
Funding
The research leading to these results has received funding from the Czech Science Foundation (CSF) through the 22–15763 project (L. Heller).
Author information
Authors and Affiliations
Contributions
AJ performed the experiments and simulations. All authors participated to the analysis of the results. All authors contributed to the writing and reviewing of the manuscript.
Corresponding author
Ethics declarations
Competing interest
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendix
Appendix
The aim of this appendix is to describe the procedure for identifying the thermophysical properties involved in the heat equation, based on the reconstruction of heat sources. To test the robustness of the procedure, thermal data are first created from a model (Section A.1), to serve as input data for processing (Section A.2).
1.1 Creation of Synthetic Thermal Data
Synthetic thermal data were created using known properties of three non-SMA materials, namely glass, titanium and aluminum (see Table 3).These materials were chosen for their wide difference of thermophysical parameter values, although the Joule effect cannot obviously be applied to the glass material.
The simulations rely on the 1D expression of the heat diffusion equation and the experimental trends at the specimen boundaries (at \(z=0\) and \(z=L)\) observed in Figs. 3 and 4. Also, the trend for the ambient temperature \({T}_{\text{ref}}\) is inspired from the experimental observations. Thus,
where \({\theta }_{\text{ref}}\) is the temperature change of the ambient air with respect to the beginning of the test \([^\circ C]\), \({\Theta }_{1}\) and \({\Theta }_{2}\) are the temperature changes at the specimen ends at the end of the heating stage \([^\circ C]\) and \(\tau\) is a characteristic time \([s]\). This latter quantity was considered here to be defined by \(\tau =\rho Cr/2h\), where, as above, \(\rho\) is the material density \([{\text{kg}}\cdot {m}^{-3}]\), \(C\) is the specific heat capacity of the material \([{\text{J}}\cdot {\text{kg}} ^{-1}\cdot {\text{K}}^{-1}]\), \(r\) is the specimen radius \([m]\) and \(h\) is the heat convection coefficient [\({\text{W}}\cdot{\text{m}}^{-2}.{\text{K}}^{-1}\)].
The initial condition \(T\left(z,0\right)\) was obtained from the numerical solution of Eq. 13 below, where the heat source in the right-hand side is equal to zero, solved using the Matlab function ‘bvp5c’:
Then, the whole temperature evolution \(T\left(z,\text{t}\right)\) over the test duration was calculated using Eq. 2. Next, Gaussian noise with a standard deviation of 20 \({\text{mK}}\) is added to mimic experimental thermal data.
Finally, filtering operations are performed to smooth the data, as it is necessary from experimental data provided by an IR camera:
-
Step 1: time filtering is performed by sliding average operation over 10 values. This reduces the noise by \(\sqrt{10}\), leading to a standard deviation of \(20/\sqrt{10}=6.3\) mK while penalizing the temporal resolution from 0.1 \({\text{s}}\) (experimental acquisition frequency of \(10 {\text{Hz}}\)) to 1 \({\text{s}}\);
-
Step 2: spatial filtering is then performed by fitting the data by a polynomial function of order 6 at each time step. This order was defined after several tests, not reported here. It provides the best results for the identification of the thermophysical parameters in the next section.
Figure 8 shows an example of synthetic temperature evolution obtained by this procedure, for titanium. The filtered data are used in the following for the processing by HSR.
1.2 Identification Procedure of \({{C}}\), \({{\lambda}}\) and \({{h}}\) from the Thermal Data
The identification of parameters from experimental data is an inverse problem that is ill-posed in the sense of Hadamard [69] due to unavoidable noise in the inputs. It is here proposed to exploit the transient state (Eq. 1) and the steady state (Eq. 2) problems to reduce the computational burden and the number of “design” variables. A steady state 1D profile was created considering the time average from 125 \({\text{s}}\) to 300 \({\text{s}}\) corresponding to SS1 in Fig. 8a (while the heat storage term is equal to zero). This additional average operation further decreases the noise in the steady state data.
1.2.1 Identification of the Specific Heat Capacity \(C\)
The specific heat capacity \(C\) is estimated from Eq. 1 at the beginning of the heating stage: see Eq. 14. At the temporal resolution of the smoothed data recording (time step of 1 \({\text{s}}\)), the evolution can be considered as adiabatic, justifying the presence of the heat storage term only in the left-hand side of the equation. The time derivatives are estimated by finite differences.
Figure 9 shows the 1D profile of specific heat capacity obtained along the wire. It can be seen that the values are nearly constant in the central part of the wire. The higher values at both ends can be explained by a lack of adiabaticity due to contact with the steel jaws (boundary conditions have been defined on the basis of experiments, see section A.1). It can be also explained by the fitting by a polynomial function of degree 6, as described in Section A.1; see also Section A.2.3 for additional information. The values given in Tables 1, 2, and 3 for the specific heat capacity \(C\) are the average over the range \([10 mm;30 mm]\); see the red line in the graph.
1.2.2 Identification of \(\lambda\) and \(h\)
To facilitate the calculation of the two remaining design variables, \(\lambda\) and \(h\), a system of equations describing the steady state and the transient state was used. We thus obtained three Eqs. (15) for two design variables (\(\lambda\) and \(h\)).
The first and second equations correspond the heating and cooling stages of the test, the term \({S}_{\text{el}}\) being null in the latter case (natural cooling). The third equation corresponds to the steady regime of the heating stage, the heat storage term being null in this case. The optimization problem was carried out using the GlobalSearch solver of the Matlab Global Optimization toolbox to apply global optimization. In addition, the Matlab optimization algorithm fmincon with large bounds was used to obtain first estimates for \(\lambda\) and \(h\). Table 3 gives the results of the identification for the three materials considered. Following comments can be made from these results:
-
the relative error on the convection coefficient \(h\) is small in the three cases of material: less than 2%;
-
the relative error in the identification of the thermal conductivity \(\lambda\) decreases when \(\lambda\) increases. While it is about 0.05% for aluminum, that for titanium is about 0.32% and that for glass is 3.84%;
-
the relative error in the identification of the specific heat capacity \(C\) is low for glass and titanium, with relative errors of + 3.85% and + 0.65%, respectively. The largest errors was found for aluminum: + 13.10%.
As a conclusion, since the thermophysical parameters of NiTi are a priori close to those of titanium (according to the literature), the following errors can be reasonably declared for our experiments: ± 0.5% for the thermal conductivity \(\lambda\); and ± 1% for the specific heat capacity \(C\).
1.2.3 Additional Information
Figure 10a presents the evolution of the spatially averaged contributions of the left-hand side of Eq. 1 as a function of time, for titanium: heat storage term \({\phi }_{\text{hs}}\), conduction term \({\phi }_{\text{cd}}\), convection term \({\phi }_{\text{cv}}\), radiation term \({\phi }_{\text{rd}}\). The reconstructed heat source \({S}_{\text{th}}={\phi }_{\text{hs}}+{\phi }_{\text{cd}}+{\phi }_{\text{cv}}+{\phi }_{\text{rd}}\) and input heat source \({S}_{\text{el}}\) are also displayed, confirming the good reconstruction. The impact of the value of \(C\) on the heat storage term \({\phi }_{\text{hs}}\) along the 50 first seconds of the heating is shown in Fig. 10b. It can be clearly seen that the slope for \(C=528 {\text{J}}\cdot{\text{K}}^{-1}\cdot {\text{kg}^{-1}}\) (the input data) tends to the heat source \({S}_{\text{el}}\) when \({\phi }_{\text{hs}}\) are extended to \(t=0 s\) (time considered for the identification of \(C\); see Eq. 14).
Considering only spatially averaged quantities for the contributions of the left-hand side of Eq. 1 makes the analysis partial. For this reason, we used the steady state equation (Eq. 2) during the heating stage to visualize the contributions along the spatial dimension: see Fig. 10c. It can be seen that, while the temperature profile is convex, the convection and radiation terms are convex whereas the conduction term \({\phi }_{\text{cd}}\) is concave. In other words, they are opposite and complementary. The use of a polynomial expression of order 6 for the fitting of the noisy temperature profiles (see Section A.1) leads to a small discrepancy between \({S}_{\text{th}}\) and \({S}_{\text{el}}\) near the ends.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Jury, A., Balandraud, X. & Heller, L. Thermal Conductivity and Specific Heat Capacity of Austenite and Stress-Induced Martensite in Superelastic NiTi at Ambient Temperature. Int J Thermophys 44, 162 (2023). https://doi.org/10.1007/s10765-023-03279-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10765-023-03279-y