Modelling the coefficient of thermal expansion in Ni-based superalloys and bond coatings

The coefficient of thermal expansion (CTE) of nickel-based superalloys and bond coat layers was modelled by considering contributions from their constituent phases. The equilibrium phase composition of the examined materials was determined using thermodynamic equilibrium software with an appropriate database for Ni-based alloys, whereas the CTE and elastic properties of the principal phases were modelled using published data. The CTEs of individual phases were combined using a number of approaches to determine the CTE of the phase aggregate. As part of this work, the expansion coefficients of the superalloy IN-738LC and bond coat Amdry-995 were measured as a function of temperature and compared with the model predictions. The predicted values were also validated with the published data for the single-crystal superalloy CMSX-4 and a number of other Ni-based alloy compositions at 1000 K. A very good agreement between experiment and model output was found, especially up to 800 ∘\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^\circ $$\end{document}C. The modelling approaches discussed in this paper have the potential to be an extremely useful tool for the industry and for the designers of new coating systems.


Introduction Background
Gas turbine components used in aero-and land-based engines are designed to operate at high temperatures and pressures in corrosive environments. These days, such conditions are satisfied by coated nickel-based superalloys [1,2]. Over the years, Ni-based superalloy compositions have evolved such that they are complex alloys for the most demanding of applications containing more than ten elements, and often cast in single-crystal form. Superalloys derive their strength and creep resistance from a precipitation-hardened microstructure comprising an ordered Ni 3 Al-based precipitate phase (c 0 ) in a solution-hardened fcc Ni (c) matrix (e.g. [3]). The mechanical properties of a modern superalloy are optimised for high-temperature creep and fatigue properties over corrosion protection, and hence, there is an increasing emphasis on coating systems for environmental and thermal protection. A typical coating system for a turbine blade for the high pressure section of a gas turbine consists of a ceramic thermal barrier top coat (TBC) and an intermediate metallic bond coat (BC) between the TBC and substrate. The ceramic TBC is normally made from yttria-stabilised zirconia (YSZ), and the overlay BC is of the MCrAlY type but, more generally, its composition can be represented by the MCrAlX type where M mainly refers to Ni and / or Co, and sometimes heavy elements such as Ta. A reactive element is represented by X, such as Y, Hf and/or Si.
The role of the BC is to provide a sacrificial oxidation barrier for the substrate, whilst also acting as a binding interface for the TBC (e.g. [4]). During service, the BC develops a surface layer of thermally grown oxide (TGO), which consists primarily of a-Al 2 O 3 . The TGO protects against further oxidation of the BC by the combustion gases which are believed to permeate the porous TBC freely. The BC contains a relatively high amount of Al and Cr to maintain a TGO layer which has to withstand many cycles of breakage (spallation) and regrowth. The high Al concentration in the BC gives rise to the NiAl (b) phase, which is commonly present with the c, and sometimes with c 0 , phases. Some BC compositions which have a large amount of Cr also contain the bcc a-Cr phase in significant amounts.
In the service environment, the Al and Cr concentrations of the BC decrease as a result of interdiffusion with the substrate and also due to consumption by the scale formation process. This results in the initial BC, which generally has a high fraction of b, transforming into a microstructure consisting of c 0 and c phases. Furthermore, as Al is removed from the BC surface by scale formation, an outer layer depleted of b is created, the extent of which is often an indicator of the remaining life of the BC. At the same time, the heavy elements present in the substrate migrate into the BC, giving rise to an interdiffusion zone (IDZ) near the BC/substrate interface. The IDZ is generally rich in c 0 and is often populated by the topologically close-packed (TCP) phases, which are considered harmful for their potential for weakening the matrix by removing solution hardening elements.
The complexity of the layered coating system can give rise to in and out of plane thermal stresses. Due to the variations in their chemistry and atomic bonding, the TBC, TGO, BC and substrate have considerably different thermal expansion characteristics. These CTE differentials can cause substantial mechanical stresses between the layers during production and service, which with additional stresses created due to the growth of the TGO, can cause de-lamination and/or cracking. Furthermore, the expansion properties of the layers gradually change over time with the evolving microstructure. Fortunately, the thermal stresses created can relax at high temperature by processes such as creep [5]. On the other hand, the structure is more vulnerable to mechanical failure during cooling due to the lack of flow processes at low temperatures and also because of the increase in brittle phases such as b in the BC. The TGO is particularly susceptible to failure, and stresses caused by thermal expansion difference are acknowledged to be the general cause for TGO spallation. Thermal stresses also set upper limits for the thickness of layers of TBC and TGO.
Therefore, it is necessary to minimise thermally induced stresses by reducing the thermal expansion mismatch between layer interfaces when designing coating systems. The BC compositions in particular can be tailored to reduce such stresses. Although the influence of phase constitution on the CTE of substrates and bond coat has been acknowledged (e.g. [6]), most traditional models are mainly empirical in nature and do not explicitly take the phase composition into account. However, there are now thermodynamic modelling tools available which can, along with a properly assessed thermodynamic parameter database, predict the equilibrium phase structure of bond coat compositions and substrates. In the modelling work described here, thermodynamic computation tools have been used to generate a more accurate prediction of the CTE of a material by accounting for individual phase volumes and their CTE characteristics at different temperatures.

Coefficient of thermal expansion (CTE)
Most materials expand when heated and shrink when cooled. The associated change in length with temperature is measured using the coefficient of linear thermal expansion a which is defined by where l is the expanded length of a specimen at temperature T and l 0 is the initial length at temperature T 0 , which is generally taken to be room ambient 20 C (293 K). The CTE is also defined as the fractional increase in length per unit temperature increment as follows: where l is the length of the test piece at temperature T and a is commonly identified as expansivity or physical expansion coefficient. The expansivity a is determined by the slope of tangent in the fractional length change versus temperature plot, whereas the mean coefficient a is measured by the slope of the chord between two points on the plot. The work here is primarily concerned with a, which is the most technologically useful form and the focus of many studies, although some of the works referenced have reported values for a. The term CTE is used in place of a throughout the article to avoid possible confusion with the a-Cr phase.

Thermal expansion coefficient of phases
Nickel-based superalloys have a phase structure consisting mainly of the c and c 0 phases. The b phase occurs in high Al-containing MCrAlY bond coatings along with the c and sometimes, c 0 phases. The a-Cr phase can also be found in high Cr-containing MCrAlY coatings. The other phases found in the systems are TCP phases, which normally include r, l and P phases.
The reported CTE values of pure c(Ni), c 0 (Ni 3 Al), b(NiAl) and a-Cr, which are the main phases found in a typical superalloy-MCrAlY coating system, are shown as a function of temperature in Fig. 1. The thermal expansion coefficients vary according to CTE Ni [ CTE Ni 3 Al [ CTE NiAl [ CTE aÀCr over the temperature range 500-1100 C as described in more detail below.

CTE of Ni
The graph of CTE of pure Ni shown in Fig. 1 has been plotted using the equation where the constants a 0 , a 1 , a 2 , a 3 and n are given as À3:810, 1.633, À1:746 Â 10 À1 , 6:510 Â 10 À3 and 1 / 3, respectively, which have been derived by regression in [7] using published data from a number of sources.
There have been a number of studies about the effect of common alloying elements on the CTE of Ni. Giamei et al. [8] have reported the effect of Cr, W, Re and Mo on CTE of c-Ni at 1000 C using Ni-10Cr, Ni-10W, Ni-8Re, Ni-20Mo and Ni-25Mo ( at.%) alloys. They found that Cr has virtually no effect on expansion, whereas W, Mo and Re decreased the CTE substantially in that order on a per atom basis, and the higher Mo containing Ni-25Mo alloy has a lower CTE than the Ni-20Mo composition.
In a study conducted at lower temperatures up to 538 C, Hull et al. [9] observed that in a 20Cr-80Ni (wt%) c nonmagnetic alloy, Mo had the most dramatic effect in reducing the a, and the effect of reduction decreases in the order of Ti, W, Si, Al, (Nb ? Ta). The elements which were found to increase CTE were Mn and Fe, in that order. However, unlike the work of Giamei et al. [8] at 1000 C, Cr was found to increase the CTE up to 23.6 wt% and then decrease.
Morrow et al. [10] measured the effects of Mo in c up to 1050 C. In ternary alloys containing approximately 33 wt% Co, 35 wt% Cr, and 32 wt% Ni, they found that addition of Mo up to 6.5 wt% clearly decreased the CTE significantly, in agreement with earlier studies.
Pavlovic et al. [11] examined the CTE of fcc Ni binary alloys Ni-Cr, Ni-Mo and Ni-Re, for alloying concentrations up to 25 at.% at temperatures up to 900 K (627 C) and reported that for all three systems, CTE decreased with increasing concentration of the second element, the addition of Cr having the least effect on CTE whereas Re the most. The findings have been compared with molecular dynamics (MD) simulation predictions carried out by Mei et al. [12], and the best agreement was found to be in the case of Mo and Re the least. A rule of mixtures equation derived by the authors shows a linear decrease of CTE with concentration for all three of these alloys, in qualitative agreement with the experimental results.
Yamamoto et al. [13] found that in Fe-free Ni alloys, Cr increased the CTE, whereas W, Mo, Al and Ti decreased it from room temperature to 700 C. The effect of Cr and W between 20 and 800 C was examined using a series of Ni-W, Ni-Cr, and Ni-Cr-W alloys by Dosovitskiy et al. [14]. The work revealed that W decreased the CTE considerably. The CTE could be suppressed with Cr at low temperature, but at high temperatures, Cr increased it.
Using a first-principles study, Kim et al. [15] have shown that CTE is reduced by the alloying elements Al, Cr, Hf, Pt, Y and Zr up to 1600 K (1327 C) and that higher Pt concentrations further decreased CTE, whereas higher Al was ineffective in altering the CTE.
These studies reveal that the CTE of Ni is decreased by most common alloying elements, with the exception of Fe and Mn. The refractory elements W, Re, Mo and Ta are the most influential in reducing the CTE, but in the case of Cr, one of the most regular additions in superalloy/coating chemistries, the behaviour seems to be more complex with respect to concentration and temperature dependency.

CTE of Ni 3 Al
The CTE curve of Ni 3 Al up to 1300 K shown in Fig. 1 has been derived from data in [16] which are based on two sources [17,18]. The ordered Ni 3 Al phase accommodates off-stoichiometric variations within a narrow concentration band, and the influence of Al on the CTE has been investigated by several authors. Williams et al. [19] have found no significant variation in well-annealed Ni 3 Al containing either 24 or 25 at.% Al. However, the measurements by Ramesh et al. [20], up to 1000 C, have shown the CTE of  Fig. 1 The variation of CTE with temperature of pure Ni [16], Ni 3 Al [16], NiAl [23] and a-Cr [24]. Ni-24 at.% Al to be slightly lower than that of stoichiometric Ni-25 at.% Al. However, Rao et al. [21] report that the CTE shows an increase with lower Al concentrations up to 1000 C, and attributes the increase to a larger number of vacancies at off-stoichiometric compositions. There have been some studies on the effect of a third alloying element. Giamei et al. [8] report that CTE is reduced by the elements Mo, Ta, Nb, W and Hf up to 1000 C. Arbuzov and Zelenkov [17] have studied the effect of Ta, Zr, V, Mn and Ti and concluded that the variation in CTE depends on the CTE of the alloying elements themselves. For example, B and Zr, which have lower values of a, decrease the CTE of Ni 3 Al, whereas Hf and Ti do not alter it significantly. Porter and Maziasz [22] have measured the CTE of the commercial Ni 3 Al alloy IC-221M which has the composition 15.9Al-8Cr-0.8Mo-1Zr-0.04B at.% and find the expansion of IC-221M to be very similar to published data for the commercial nickel-based superalloy Inconel 713C [16], possibly due to their similarity in composition.

CTE of NiAl
The graph of CTE of NiAl shown in Fig. 1 has been plotted using the equation where T is temperature in K and is valid between 300 and 3000 K. The equation was derived by Noebe et al. [23] considering a number of published data for polycrystalline materials. It is noted that the CTE of NiAl is a strong function of temperature but not of composition. To the authors' best knowledge, there are no reported investigations on the variations of CTE with off-stoichiometric compositions in NiAl, nor on third or higher element additions.

CTE of a-Cr
The a-Cr phase, mostly found in MCrAlY coatings containing a high amount of Cr, and sometimes also in superalloys, is based on the BCC a phase of the element Cr. The graph shown in Fig. 1 is derived from data in [16] for pure Cr and is in good agreement with reported figure for a-Cr by Fritscher et al. [24] who observe that the CTE remains fairly constant at 10:0 Â 10 À6 C À1 from room temperature up to 1000 C.

Modelling of CTE
Appropriate modelling techniques can be divided broadly into three categories as (a) rule of mixture of elements, (b) empirical models and (c) phase-structure-based approaches.

Rule of mixture of elements
In this approach, the CTE of an alloy is composed of individual contributions of constituent elements, and in the simplest form, a simple rule of mixtures formulation can be given as where a and a i are the CTEs of the alloy and constituent element i, respectively, and V i is the atomic volume fraction of the element i.
The ability of heavy refractory elements such as Re, W, Mo, Ta and W to lower the CTE of nickel-based alloys is well documented as discussed earlier (e.g. [8,11,14]), and many workers have also identified these elements to have much lower CTEs themselves compared to Ni. Therefore, at least at first sight, a rule of mixtures model seems very plausible. The graphs in Fig. 2a, b illustrate this point and show how CTE varies between 500 and 1150 C for a number of pure elements often found in the nickel-based superalloy and bond coat compositions.
Among the refractory elements, W and Mo have the lowest thermal expansion characteristics with a CTE smaller than 6 Â 10 À6 C À1 . Figure 2b concentrates on the lower half of the graph, with elements having melting points above 2200 C and CTE between 6.0 and 8.5 C À1 .
The group of elements Cr, Ti and V, which are adjacent in the periodic table, show close expansion coefficients which vary between 9.0 and 11.0 C À1 . The drop observed of CTE in the case of Ti at 882 C is due to the phase transition from hcp (a) to bcc (b). The CTE of Co is similar to and slightly higher than that of Ni. Aluminium has the highest CTE. Since thermal expansion is related to atomic bonding, which also affects the melting point (MP), it is perhaps not surprising that W with the highest melting point 3420 C has the lowest CTE of the group and Al with the lowest MP (660 C) has the highest CTE.
Several authors have commented that although the presence of low expansion elements in the alloy composition reduces the CTE of an alloy, a simple rule of mixtures approach cannot be adopted universally to predict the expansion behaviour of an alloy. Pavlovic et al. [11] in their study on the effect of Cr, Re and Mo on the thermal expansion of Ni argue that the simple rule of mixtures seems to work best when the atomic volumes of the atoms in the alloy are not too different. They found that the experimental decrease in CTE with the at.% of Mo and Re is three to four times larger than the prediction of the simple rule of mixtures equation, whereas for Cr, the experimental and predicted values are quite comparable. Furthermore, the authors argue that the predicted decrease in CTE is the result of the lower values of CTE for Cr, Re and Mo relative to that of Ni, offset by the atomic volume contributions of Cr, Re and Mo relative to Ni which has the effect of increasing CTE. The simple rule of mixtures approach in Eq. (5), however, becomes more complicated when predicting the CTE of a multi-phase alloy, as a i and V i are likely to be different for each phase.

Empirical models
These models have been developed by fitting polynomials to a set of experimentally measured CTE values using regression analysis. The effects of Cr, Fe, Mo, W, (Nb ? Ta), Mn, Si, Ti and Al on CTE of Ni-20Cr wt% alloys were modelled by Hull et al. [9] using a linear form up to 1000 F (538 C) where b i are the regression coefficients and C i are concentrations of alloying elements i. The contributions of second order terms were either too small (or not considered) except in the case of Cr, which has shown a maximum at around 23.6 wt% through the C 2 Cr interaction. The above model has been validated against a number of commercial alloys.
For Ni-based alloys containing Cr, W, Mo, Al and Ti, Yamamoto et al. [13] derived the following formulation for the CTE of Fe-free Ni alloys from room to 700 C: where element concentrations are given in wt%. The results of work by Dosovitskiy et al. [14] on elements Cr and W to 800 C using a series of Ni-W, Ni-Cr, and Ni-Cr-W alloys are shown by the following derivation: where concentrations are given in at.%. Sung and Poirier [7] have derived a temperature-and concentration-dependent empirical CTE model for Ni-based superalloys comprising c and c 0 phase mixtures by fitting experimental data determined up to 1300 K (1027 C) to the following form of linear relationship using regression where C i are the concentrations of element i in wt% and d j ,   Fig. 2 The variation of CTE with temperature of a number of common elements a and b, found in nickel-based superalloys and MCrAlY bond coat compositions with b concentrating on elements with a low expansion coefficients [16], reported by [16], except Ni data which were taken from [7]. where T is temperature in K, b i , i ¼ 1. . .3 are fitting coefficients and n is 1 / 3 as for Ni in Eq. (3). The b i up to 1300 K are given as b 0 ¼ À5:9431, b 1 ¼ 2:5805, b 2 ¼ À0:30383 and b 3 ¼ 1:2181 Â 10 À2 . However, since extrapolation beyond 1300 K with these b i gives rise to too much deviation from the CTE of c-Ni given in Eq. (3), the following set of b i has been derived for T ! 1300 K by the authors: b 0 ¼ 6:9411 Â 10 À2 , b 1 ¼ 0:50692, b 2 ¼ À6:5013 Â 10 À2 and b 3 ¼ 3:0825 Â 10 À3 . The derivation is based on the observation made by Morrow et al. [10] that Ni-based superalloys containing c and c 0 have a smaller CTE compared to a single phase c. When c and c 0 alloys are heated and transform to c entirely, the CTE increases by an increment that depends on the fraction of c 0 . Morrow et al. [10] quantified the increase of CTE for alloys containing Al and Mo as 0.6Â10 À6 C À1 .
Sung and Poirier [7] have also determined the temperature dependence of d i (i [ 0) by fitting experimental data at 400, 700, 1000 and 1100 K to the following equation, and the fitted coefficients p i and q i are shown in Table 1.
Amongst the empirical models, the scheme of Sung and Poirier [7] encompasses the widest range of components, compositions and temperatures applicable for Ni-based superalloys. However, their model is derived for superalloy compositions, consisting mainly of c/c 0 mixtures and hence care has to be taken when the scheme is extended to bond coat compositions which normally contain the phase b and sometimes other phases such as a-Cr in high fractions.

Mixtures of phases
The CTE of an alloy can also be modelled by taking contributions of the thermal expansions of individual phases into account. A simple rule of mixtures approach can then be given by where a is the CTE of the multi-phase body, a i are the CTEs of individual phases i, and V i are the volume fraction of phases i. The method can be used to recognise the micro-stresses, termed tessellated stresses, which are known to develop between differentially expanding phases. The most recognised formula which accounts for the elastic stresses is given by the Turner's equation [25,26] where the additional K i term represents the bulk modulus of phases i. The formula assumes only uniform hydrostatic stresses to exist within the phases. If the Poisson's ratios of the phases are identical, Young's moduli E i can replace bulk moduli K i in Eq. (13) which gives a ¼ In general, Turner's equation yields lower CTE values compared to the simple rule of mixtures method because the phases with low expansion coefficients are likely to have high elastic moduli [26]. Hermosilla et al. [27] calculate the CTE of a phase aggregate using a modified version of an expression derived by Wakashima et al. [28] given by a ¼ where a i are CTEs of individual phases and term w given as where m is the Poisson's ratio of the aggregate. For an isotropic phase aggregate, the elastic properties are defined by the bulk K and shear moduli G which can be calculated using bulk K i and the shear G i moduli of constituent phases i following the expressions originally derived by Budiansky [29], which are in turn based on Eshelby's inclusion technique [30,31].
where / and u are defined by following expressions: Using the following standard relationship available for isotropic properties of phases equations (15)- (21) can be used to calculate the CTE of the phase aggregate using a numerical root finding scheme for solving a non-linear systems of equations (e.g. [32]). Although the phase mixture approaches are powerful, they have not been used widely in the past due to the difficulties in accurately determining phase volume fractions and elastic properties at high temperatures. Fortunately, due to the advancements in thermodynamic equilibrium calculation software such as MTDATA [33] and the wide availability of assessed thermodynamic parameter databases such as Ni DATA for Ni-based alloys [34], the task of determining phase fractions as a function of temperature and composition has now become possible for Ni-based superalloys and coatings (e.g. [35,36]).
However, for the use of methods by Turner [25] and Hermosilla et al. [27], calculation of elastic properties of phases as a function of composition and temperature remains a problem for Ni-based alloys, and especially coatings, as a result of general lack of high-temperature data and suitable models.

Modelling strategy
In this work, the capability of the model introduced by Hermosilla et al. [27] (''Mixtures of phases'' section) has been further developed to include temperature-dependent elasticity, density and phase composition data to calculate the CTE of a number of Ni-based superalloys and bond coats. The model predictions have then been compared with experimental data, both new data determined during the course of this work and those available in literature. The model predictions have also been compared with two other modelling schemes, a simple rule of mixtures method and Turner's scheme [25] (''Mixtures of phases'' section), both of which take into account the phase constitution of the materials. The simple rule of mixtures method adopted here is illustrated in Fig. 3. The phase constitution of the material was calculated using thermodynamic equilibrium calculation software (''Determination of phase constitution'' section) which returned the equilibrium weight fractions of phases for a given combination of chemical composition, temperature and pressure. The weight fraction x i of each phase i predicted to be present at the temperature was used to calculate the volume fraction of each phase present using the densities of each phase q i . The volume-weighted CTE values of individual phases a i were combined using a simple rule of mixtures method given in Eq. (22) to obtain a value for the CTE of the material. a ¼ Determination of phase constitution The equilibrium phase structure as a function of material composition, temperature and pressure was calculated using the Version 5.02 of the Application Program Interface (API) of MTDATA software (Linux Version) [33] together with the Ni-DATA thermodynamic parameter database Version 3.1 [34].

Modelling of the bulk modulus of phases
For the implementation of the methods of Turner [25] and Hermosilla et al. [27], the elastic properties of the constituent phases of the material are needed. Systematically derived elastic properties of multi-component alloy phases are unavailable, and hence, data for pure materials are used in this work. The values for bulk modulus (K) of c-Ni up to 762 K have been derived in [37] using experimental data determined by Alers et al. [38], and these have been adopted here by fitting a third-degree polynomial of the following form: where K i are given in GPa and T in K. The fitted polynomial coefficients a i are given in Table 2 which also shows the coefficients obtained in a similar manner for c 0 up to 1100 K derived from work of Tanaka and Koiwa [39] and for b up to 875 K by Rusović and Warlimont [40]. For higher temperatures, isentropic bulk moduli calculated using a first-principles approach by Wang et al. [37] Concentration,Temperature were used, after fitting the third-degree polynomial given in Eq. (23). The fitted polynomial coefficients are given in Table 2.
For the c and c 0 phases, data calculated without the quasistatic approximation were used, whereas for the b phase, those obtained with the quasistatic approximation were used, as these sets showed best correspondence with lower temperature experimental data in [38][39][40]. Due to the unavailability of moduli data for the a-Cr and r phases, the bulk modulus of 161 GPa reported for pure Cr at room temperature [41] was used for both of these phases as an approximation, due to the high amount of Cr likely to be present in these phases.

Modelling of Poisson's ratios of the phases
The Poisson's ratios of c, c 0 and b phases were calculated using modelled stiffness coefficients c 11 and c 12 in [37], which showed reasonable agreement with low-temperature experimental data, using the following relationship which holds for isotropic materials: The stiffness parameters as well as the Poisson's ratios m i of phases i were modelled using a third-degree polynomial of the form given in Eq. (25), where temperature T is given in K. The coefficients of best fit are given in Table 3.
The Poisson's ratios of a-Cr and r phase have yet to be determined, and hence, the reported value of 0.21 for pure Cr at room temperature [41] was used for both of these phases here as an approximation.

Modelling of the density of phases
The density q i of each phase i of the material was used to calculate the volume fraction v i of the phase using the standard expression where w i is the mass fraction of phase i returned from the thermodynamic calculations. The density q of the c and c 0 phases was determined using the standard expression for fcc crystals where M j and X j are the atomic weight and fractional concentration of element j, N 0 is Avogadro's number and a is the lattice parameter of the alloy. The lattice parameter was modelled using the Vegard's Law [42] a where a 0 are 3.5219 and 3.5691 Å for c and c 0 , respectively, and oa=oX j for various elements in c and c 0 were obtained from work by Kablov et al. [43] and are shown in Table 4.
The density q of a phase at temperature T was determined using the general relationship where q 0 is the density of the phase at room temperature T 0 and a its CTE. In the absence of a model for calculating density as a function of composition in the c phase, the reported q of pure Ni was used in Eq. (29). The density of b-NiAl was given as a function of Ni concentration by Noebe et al. [23] using data from a number of sources. The density decreases linearly with decreasing Ni content, although a change in slope of the line occurs at the stoichiometric composition which has a density of 5.9 g cm À3 . Due to the discontinuity at the stoichiometric concentration, 50 at.% Al, two schemes are used to determine the density q NiAl above and below this composition which are modelled as follows: for 40\C Ni \50 at.% and for 50:0\C Ni \55:0 at.% Implementation of the models The models were coded using the C language with the MTDATA software library being linked via its Application Program Interface (API) [44]. In order to solve the nonlinear system given in Eqs. (15)- (21) for the scheme of Hermosilla et al. [27], a C routine based on that given in [32] for the multidimensional downhill-simplex method of Nelder and Mead [45] was used.

Experimental method
In order to provide additional validation, the thermal expansion in an IN-738LC superalloy and Amdry-995 bond coat material was measured in a Linseis mechanical dilatometer equipped with an alumina test piece holder and push-rod. The nominal chemical compositions of the samples are given in Table 5.
The instrument was calibrated in accordance with ASTM E228/CEN EN821-1. The transducer sensitivity was determined using a calibrated drum micrometer, and the instrument baseline shift was determined using a piece of push-rod alumina as a test piece. The apparatus expansion correction was determined using a platinum test piece and reference data for platinum taken from [46]. For each test piece, at least three thermal cycles were employed at a heating/cooling rate of 2 min À1 using a flushed argon atmosphere to restrict oxidation. The first thermal cycle was used for settling and annealing of the test piece, and the thermal expansion data were taken from second and third cycles as the average of heating and cooling. The accuracy of the determined expansion coefficients is considered to be about AE0.1Â10 À6 C À1 for temperature ranges from room to above 400 C.
The measurements were carried out in IN-738LC samples cut in the shape of a bar with approximate dimensions 28 Â 8 Â 3 mm 3 to 1200 C. A delaminated layer of sample was used in the case of Amdry-995 and reliable measurements could be collected only up to 800 C due to the fragility of the sample.

Results and discussion
The graphs in Figs. 4a-c show a comparison between the measured and calculated CTEs of superalloys IN-738LC and CMSX-4 and also for the bond coat Amdry-995. The experimental measurements of CTEs in the IN-738LC and Amdry-995 samples were made during the current study, and those in CMSX-4 were extracted from the study in [47]. The nominal compositions of the materials are listed in Table 5, which were used as input for the modelling (except Y which was not supported by the thermodynamic database). The calculated CTE profiles shown in the Fig. 4 have been produced using three schemes: (1) that proposed by Hermosilla et al. [27] which has been extended in this work with the incorporation of temperature-dependent elastic, density and phase composition data, (2) Turner's method [25] described in ''Mixtures of phases'' section using temperature-dependent elastic, density and phase composition data and (3) the simple rule of mixtures (ROM) method described in ''Modelling strategy'' section with temperature-dependent density and phase composition data.
As the graphs reveal, the model predictions are in very good agreement with experiment measurements. At higher temperatures, above 950 C in IN-738LC and above 850 C in CMSX-4, the model starts to underpredict the CTE slightly. In Amdry-995, again an excellent match could be found within the temperature range 400-800 C but above 650 C, the calculations over-predicts the measured CTE by a small margin.
The discrepancies seen between experimental and modelled profiles can be related to the phase structure  variation between samples. Figure 5 shows the predicted phase structure variation with temperature for each sample calculated in this work. For the superalloy substrates IN-738LC and CMSX-4, the phase structure at lower temperatures consists of c, c 0 and r phases. As temperature is increased, first the r and then c 0 is dissolved from the structure, and the temperature at which the predicted CTEs start to deviate from measured values coincides with the disappearance of the r phase. Therefore, the modest disagreement observed between the experimental and modelled curves at higher temperatures in Fig. 4a, b could be due to a delay in dissolution of the r phase in particular or due to an underestimation of the CTE values of c and c 0 phases at higher temperatures. The model however is still capable of estimating the CTE with a high degree of accuracy in the temperature range where these materials are commonly used in practice.
In the case of the Amdry-995 bond coat sample, the phase structure is predicted to contain also the b phase in addition to c, c 0 and r phases, as Fig. 5c reveals. Above 500 C, the amount of r and c 0 gradually decreases while the c and b content does the opposite, although beyond 700 C, the b concentration declines. Despite the complexity in the phase transitions, the model predicts the CTE of this sample reasonably accurately as evident in Fig. 4c, with only a slight deviation occurring at around 620 C which overlaps with the substantial phase structure change observed in Fig. 5c.
The CTE calculation scheme of [27] requires the determination of elastic properties of the specimen material using the properties of its constituent phases, as outlined in Eqs.    [27]. Turner and ROM refer to the models by Turner [25] and the simple rule of mixtures method described in ''Modelling strategy'' section and Fig. 3, respectively. the moduli of those two phases, as evident indeed in the case of CMSX-4. However, the presence of a fairly high fraction of a, r or other TCP phases at lower temperatures, such as in IN-738LC, can cause the modulus to be lower than that of c at low temperatures. It is noted here that due to the lack of data, the temperature variation of the r and a is assumed to be constant in the current model. In Amdry-995 and other such coating compositions, the high proportion of the b phase present, which can be stable up to high temperatures, can increase the modulus of the coating over that of a superalloy composition, as illustrated in Fig. 6. Figure 7 shows the variation of calculated shear modulus of c, c 0 and b phases using the model given in ''Modelling of the bulk modulus of phases and Modelling of Poisson's ratios of the phases'' sections and also that of Amdry-995 computed using the moduli of phases via Eqs. (17)- (21). As shown, the shear modulus G i of phases i shows the relationship G r [ G c [ G c 0 [ G b . At temperatures lower than 650 C, the shear modulus of Amdry-995  is higher than that of c as the phase structure has a high proportion of r phase which is shown to have a high shear modulus. As temperature is increased, the shear modulus of Amdry-995 goes lower than that of c because of the dominance of the c and b phases at these temperatures which have lower moduli compared to the r phase. The overall agreement shown in Fig. 4 between experimental measurement and predictions by a number of models is very good. All three models account for the phase constitution of the materials but the schemes provided by Hermosilla et al. [27] and Turner [25] additionally consider the elastic properties of constituent phases.
At above 850 C in the case of superalloys and IN-738LC in particular, the model predictions deviate from measured values. The method of simple rule of mixtures gives much better predictions, up to 1100 C, than the models which incorporates elastic properties of phases, which tend to underpredict the CTEs at higher end of the temperature scale. It could be possible that at high temperatures, the modelled elastic properties are less accurate and further work is needed to develop that part of the model.
It should also be mentioned that during experimental measurement at a typical heating/cooling rate, the equilibrium amounts of phases are never present, possibly more during heating and less during cooling from above the solvus. There is often hysteresis in thermal expansion curves in the apparent solvus temperature range. These artefacts may well have affected the agreement between model predictions and experimental data observed.
Interestingly, both methods that incorporate elastic properties of phases give very similar predictions for the examined alloy compositions over the temperature range considered. While the approach by Hermosilla et al. [27] provides a rigorous method of taking account of interactions between phases, in practice, almost exactly the same results have been achieved for the specific materials considered here using simpler methods which do not need so much data. The main purpose of the method by Hermosilla et al. [27] is therefore to provide a rigorous validation of what are effectively semi-empirical phase mixture models which might or not apply to other multiphase alloys.
The graph in Fig. 8 shows how the predicted CTE values of a number of Ni-based superalloys and coatings compare with the measured CTE figures available in literature at 1000 K (726.85 C). There is in general a good agreement between experiment and modelling with many points overlying on the dashed line which represents exact one to one correspondence. The calculated root mean square value of the difference between the model and measurement figures of all the points shown in Fig. 8 is 0.6103 Â 10 À6 C À1 , indicating that the model can predict the CTE with a high degree of accuracy in this case.
It should be mentioned, however, that the CTE of a material can also be affected by factors such as grain size, homogeneity and crystalline texture which the current model does not take into account. Furthermore, the thermodynamic modelling procedure used in this work has limitations, for example, some elements and phases may not be supported. Although the predictions are very good as discussed, some of these factors may have contributed to the discrepancies observed in Figs. 4 and 8.

Conclusions
Most of the models currently available for predicting CTEs of nickel-based superalloys are empirical in nature and therefore are restricted by the temperature or the  Amdry-995 Fig. 8 Comparison between experimentally measured CTE values and those predicted using the model described in [27] that has been extended in this work. The values shown are for 1000 K (726.85 C) for a number of superalloys. The experimental data are mainly from [7] (open markers). composition range in which the experiments were carried out. Models based on foundations of materials physics are very limited. The recent advances in modelling thermodynamic properties based on CALPHAD techniques and the availability of reliable assessed databases for Ni-based alloys make it possible to predict phase constitution of such alloys accurately. As material properties such as CTE are very likely to be dependent on the phase structure and the CTEs of individual phases in an alloy, it is possible to exploit advances made in such thermodynamic modelling to develop better CTE prediction models. In this work, a model originally proposed by Hermosilla et al. [27] has been further developed to incorporate temperature-dependent elasticity, density and phase composition data, and the CTE predictions by the extended model have been compared against measured CTEs of superalloys IN-738LC, CMSX-4 and bond coat Amdry-995. The predictions agree very well at lower temperatures, but the model slightly under-predicts the CTE in the superalloys considered at higher temperatures. For the Amdry-995, IN-738LC and CMSX-4 compositions examined, the method of [27] gave close results to those predicted by the model proposed by Turner [25]. As both methods are based on the same elastic properties and phase prediction models, comparable results can be expected. The work also found that simple volume-weighted rule of mixtures method, which does not take elastic properties into consideration, can also predict experimental measurements of CTEs appreciably well, and even give better agreement at high temperatures. It has also been identified that further work is still needed to improve the elastic parameter models.