Ada 344715
Ada 344715
U:.j> ^ST^K^NT A
Papers presented at the 85th Meeting of the AGARD Structures and Materials Panel, held in
Aalborg, Denmark 15-16 October 1997.
19980520 102
North Atlantic Treaty Organization
Organisation du Tratte de l'Atlantique Nord
According to its Charter, the mission of AGARD is to bring together the leading personalities of the NATO nations in the
fields of science and technology relating to aerospace for the following purposes:
— Recommending effective ways for the member nations to use their research and development capabilities for the
common benefit of the NATO community;
— Providing scientific and technical advice and assistance to the Military Committee in the field of aerospace research
and development (with particular regard to its military application);
— Continuously stimulating advances in the aerospace sciences relevant to strengthening the common defence posture;
— Improving the co-operation among member nations in aerospace research and development;
— Providing assistance to member nations for the purpose of increasing their scientific and technical potential;
— Rendering scientific and technical assistance, as requested, to other NATO bodies and to member nations in
connection with research and development problems in the aerospace field.
The highest authority within AGARD is the National Delegates Board consisting of officially appointed senior
representatives from each member nation. The mission of AGARD is carried out through the Panels which are composed of
experts appointed by the National Delegates, the Consultant and Exchange Programme and the Aerospace Applications
Studies Programme. The results of AGARD work are reported to the member nations and the NATO Authorities through the
AGARD series of publications of which this is one.
Participation in AGARD activities is by invitation only and is normally limited to citizens of the NATO nations.
* AGARD merged with the Defence Research Group of NATO (DRG) on 1 January 1998 to form the Research and Technology Organization (RTO)
of NATO. However, both AGARD and DRG will continue to issue publications under their own names in respect of work performed in 1997.
®
Printed on recycled paper
ISBN 92-836-1073-3
Executive Summary
Thermal barrier coatings are an emerging technology which will allow either increasing the inlet
turbine temperatures or on the other hand decreasing the working temperature of the metal of the blades
and consequently increasing their life-time.
The Workshop allowed a survey of the state of the art, a description of the existing technologies or of
the technologies under development, a review of the present knowledge of damage mechanisms,
including microstructural, mechanical and thermal aspects, and an account of the advantages and
drawbacks of the various families as perceived by engine manufacturers and users. A final discussion
was held to identify the needs for further R&D.
Les revetements anti-mur de chaleur
(AGARD-R-823)
Synthese
Les barrieres thermiques sont une technologie emergente qui permettra soit d'augmenter les
temperatures d'entree des turbines, soit de reduire la temperature ä laquelle est soumis le metal des
aubes et par consequent d'en accroitre la duree de vie.
L'atelier a eu pour objet de presenter l'etat de l'art, de donner une description des technologies
existantes ou en cours de developpement, de donner une revue des connaissances actuelles sur les
mecanismes d'endommagement includant les aspects microstructuraux, mecaniques et thermiques, de
rendre compte des avantages et des inconvenients des diverses families tels que les percoivent les
fabricants et les utilisateurs de moteurs.
Une discussion finale a permis d'identifier les besoins futurs en R & D.
Contents
Page
Synthese IV
Preface Vll
Preface Vlll
Reference
SESSION I: INTRODUCTION
Paper 2 Withdrawn
Paper 19 Withdrawn
Thermal barrier coatings are a technology used to protect the metal of combustion chambers, turbines blades or
veins. Since these parts are air-cooled low-thickness elements, it is possible to reduce the actual temperature of
the metal introducing a large temperature gradient though the ceramic coating. This technology has been used for
several decades in the combustion chambers, and more recently on modern highly-cooled blades.
Since it can allow an increase of the temperature difference between the metal part and the gas flow, up to 100°C,
it can be considered as one of the bottle-neck technologies giving access to higher temperatures in aeronautical
engines. Alternatively, it can allow also, for the same turbine inlet temperature, a reduction of the working
temperature of the metal of the blades, leading to Longer lifetimes of the parts, or a reduction of the cooling air
flow, with an increase in engine efficiency.
Up to now zirconia-base ceramics in particular yttried-zirconia have been used; they are compatible with metal
temperatures up to 900°C for combustion chambers and 1050°C for blades, under fairly low thicknesses.
Two deposition processes are presently in use: plasma spray, which gives relatively thick, low cost, low
conductivity coatings, well adapted to combustion chambers, and for blades Physical Vapour Deposition
techniques which give columnar microstructure which lead to more acceptable damage and to a surface of limited
rugosity, which does not perturb the aerodynamic qualities of the blades but with unfortunately lead to lower
deposition rates, and also higher conductivities.
Research is in progress in several directions: improved or new deposition techniques, new composition of the
ceramics.
1. The presentation of families and technologies under development with reference to new ideas and limits.
2. The description of the advantages and drawbacks of the various families of thermal coatings as perceived by
engine manufacturers or users.
4. Modeling the behaviour and the lifetime of the coatings with special reference to the underlying bond coat.
P. COSTA
Chairman of the Subcommittee
Preface
Les barrieres thermiques sont une technologie qui est utilisee pour proteger le metal des chambres de combustion,
des aubes de turbine mobiles et fixes. Comme ces pieces sont des elements de faible epaisseur fortement
refroidies par un courant d'air arriere, il est possible de reduire la temperature effective du m6tal du fait d'une
importante chute de temperature au travers du revetement ceramique. Cette technologie est utilisee depuis
plusieurs decennies pour les chambres de combustion, et depuis plus recemment pour des aubes mobiles
fortement refroidies.
Dans la mesure oü eile permet d'accroitre la difference entre la piece metallique et le flux gazeux de pres de
100°C, eile peut etre considered comme une technologie goulet pour qui veut accroitre les temperatures dans les
moteurs aeronautiques. En parallele, eile permet aussi, pour une meme temperature d'entr^e de turbine, de require
la temperature du metal des pieces travaillantes, et d'en accroitre ainsi la dur6e de vie, ou encore de require les
flux de refroidissement, avec un accroissement du rendement des moteurs.
Jusqu'ici, ce sont des ceramiques ä base de zircone qui ont ete utilisees, en particulier la zircone-yttriee. Elles sont
compatibles avec des temperatures atteignant 900°C pour le metal des chambres de combustion, de 1050°C pour
les aubes mobiles, et ceci pour des epaisseurs de revetement relativement faibles.
Deux precedes de depot sont actuellement utilises, la pulverisation plasma qui permet d'obtenir des revetements
relativement epais, de coüt modere et de faible conductibilite thermique, bien adaptes au probleme des chambres
de combustion, et pour les aubes des depots en phase vapeur, qui donnent des structures colonnaires qui sont
faiblement endommageables, et un 6tat de surface de bonne qualite ne perturbant pas l'aerodynamique de l'aube,
avec toutefois des vitesses de depot, done des coüts plus eleves, et egalement une conductibilite plus eleväe. Les
recherches se poursuivent dans plusieurs directions: methodes de depots ameliorees ou nouvelles, nouvelles
compositions de la c6ramique.
2. La presentation des families et des technologies en cours de developpement, avec une reference particuliere
aux idees nouvelles et ä leurs limites.
3. La description des avantages et des inconv6nients des diverses families de barrieres thermiques tels que les
percoivent fabricants de moteurs et utilisateurs.
4. La comprehension des mecanismes d'endommagement en correlation avec les proprietes des barrieres et leur
microstructure.
5. La modelisation du comportement et de la duree de vie des revetements, avec une reference particuliere pour
le revetement interm6diaire.
P. COSTA
President du Sous-Comite
Structures and Materials Panel
SUB-COMMITTEE MEMBERS
PANEL EXECUTIVE
ABSTRACT. In this paper, the importance of the the well-known t' phase [5] and its toughening
challenge associated with the control of the thermal mechanisms [6]. Important materials studies have also
conductivity of thermal barrier coatings for turbine been devoted to finding an optimal metallic bond
engines hot stages is being reviewed (yttria stabilised coating to support the ceramic topcoat [7,8]. The
zirconia mostly). It is firstly illustrated by the development of TBC has also been closely linked to
description of a practical aeronautic coated and that of the corresponding deposition processes: plasma
uncoated turbine blade design exercise. The various spraying and electron beam physical vapour deposition
contributions to TBC thermal conductivity are then (EBPVD). A number of studies [4, 9, 10] have thus
reviewed. Their low conductivity finds its source not been devoted to optimise the deposition process
only in the nature of the ceramic layer (highly parameters to obtain favourable structures in order to
disordered material), but also in the morphology of achieve longer life cycles. The general idea there, was
the insulating layer, closely linked to its fabrication to try to relate coatings microstructure to their
process parameters. The influence of various factors thermomechanical resistance.
(such as yttria content, temperature, porosity content
and distribution, etc.) on the thermal conductivity is This necessity to address in the very first place the
examined, largely based on a literature review. In this feasibility and reliability problems for such coatings
field of investigation, the modelling tool should allow has somewhat hindered, until recently, the extent of
to predict, to a certain extent, which morphology research activities directly dealing with their thermal
would lead to the lowest values of thermal conductivity (K). The implicit notion there, was that
conductivity. Eventually attempts are made to identify by using stabilised zirconia as TBC would lead to a
research domains where further understanding is value of K for the coating not significantly different
needed, and to formulate several suggestions from one coating to the other, lower in any case than
concerning possible ways to lower the thermal that of the dense material having the same
conductivity. composition. The scope of this paper is to emphasise
the importance of the challenge associated with the
1. INTRODUCTION control of coating thermal conductivity. It is first
Yttria stabilised zirconia based thermal barrier illustrated through a practical aeronautic turbine blade
coatings are viewed today as the new generation of design. The various contributions to thermal barrier
protective coatings that will allow to take the next coating thermal conductivity are then reviewed. The
quantum step forward to increase turbine engine low conductivity of thermal barriers comes on the one
efficiency, without having to pay the price of hand from the nature of the yttria stabilised zirconia
replacing nickel-base superalloys by some other more ceramic layer, a highly disordered material indeed, and
refractory material [1]. They are also considered by on the other, from the morphology of this layer, closely
diesel engine manufacturers to increase the efficiency linked to coating processes: pores, cracks, etc. The
of their engines. These coatings have been studied influence of various factors (such as yttria content,
and developed for more than 20 years in the western temperature, porosity distribution, etc.) on the thermal
countries and even longer in former USSR. conductivity is examined, largely based upon a
literature review. In this field of investigation, the
It can be said that one of the first concern of materials modelling tool should allow to predict, to a certain
scientists in the TBC field was to optimise the extent, which type of morphology would lead to the
coatings from a materials science standpoint. In the lowest values of conductivity. Eventually, attempts are
early days, important studies have been carried out to made to identify research domains where further
determine for instance which was the most adequate understanding is needed and to formulate several
stabilising agent of zirconia for a given application, suggestions concerning possible ways to lower the
and in what amount [1-3]. As far as very high thermal conductivity.
temperatures applications are concerned (mostly
driven by the aero industry), there was a rapid and 2. A REAL CASE DESIGN STUDY
general move towards 6-8 wt. % yttria partially In this case we consider a high pressure turbine vane
stabilised zirconia (Y-PSZ) because these coatings airfoil of an aero-engine. The component is meshed for
showed the highest degree of resistance to spallation 2D finite element calculation ; for one airfoil, about
in thermal fatigue solicitations [3] and an excellent 350 elements and 1200 nodes are being generated. The
thermal stability. Later, explanations of the thermal barrier coating is represented as a single
phenomenon appeared through the identification of element in the foil thickness.
^1300 N>7BC
A mission profile is defined for the component, as the
variation with time of relevant engine parameters such Ü1200 I.JVVMII.IV
Examination of the leading edge metal temperature 100 1000 10000 100000
profile (figure 1) shows that the application of a TBC
(1.9 W/m.K) reduces the maximum temperature Component life in cycles ■ 1.1 W/m.K
(1090°C) of "the metal by 50°C. Decreasing the thermal H1.9W/m.K
conductivity down to 1.1 W/m.K permits to gain a D No TBC
further 32°C. Considering the whole leading edge
profile, the maximum gain even reaches 113°C. Figure 3 : Component life duration (creep andfatigue)
If we consider now that the engine is run so that the
components keep the same life duration, but the gain is
concentrated entirely on reducing the cooling rate and
by way of consequence decreasing the fuel
1-3
consumption, table 1 shows that up to 0.7% specific cooling, as a result of the limited elastoplastic
consumption can be saved in the cases envisaged. deformation capability of the ceramic. The other is
mainly oriented parallel to the surface and corresponds
Thermal conductivity Specific to poor contacts between different splats. The latter,
(W/m.K) Cooling rate consumption being perpendicular to the heat flux are particularly
1.9 -15% - 0.4% effective in reducing the thermal conductivity of the
1.1 - 25% - 0.7% coating. In addition, a dispersed, spheroidal porosity
contributes to lowering the conductivity.
Table 1 : Reduction in cooling rate and specific
consumption owing to the application of a TBC of In the case of EBPVD coatings, the reduction in
thermal conductivity K on a vane airfoil (for a given thermal conductivity comes mainly from the
lifetime). intracolumnar fine porosity, the intercolumnar porosity
being much less effective as it is arranged
These few results clearly demonstrate the importance perpendicular to the surface. Schematically, the
of lowering the thermal conductivity in terms of columnar coating can be considered, in an electrical
economy and performance. Moreover, a decrease in analogy, more as a set of parallel resistances than a
the thermal conductivity permits to lower the series set which gives a closer image of the
thickness of the coating. Given the fact that a TBC (microcracked) sprayed coating.
represents an added stress on a rotating blade of 10%
(for a 125 urn thick ceramic layer), lowering its
thickness (provided the density is kept constant) will System Thermal conductivity Comment
result in a significant' gain on the stress applied on the Reference
rotating parts, blades and disks. 300K 1300 K
Zr02-7%Y2O3 1.7 2 [4]
EBPVD
Of course, a rigorous computation would have to take Zr02-7%Y2O3 0.9 [12]
into account the radiation contribution to the APS
conductivity within the coating. This is an effect Zr02-7-9%Y2O3 0.9 to 1.1 [13]
which depends on the actual environment seen by the APS Various
airfoils within the turbine stage [11] and difficult to porosity levels
assess at present as the optical properties of zirconia- Zr02-20%Y2O3 0.4 0.55 As sprayed,
based plasma-sprayed and EBPVD coatings are APS 7% porosity
poorly known. ri4i
Zr02-20%Y2O3 1.5 1.7 lOOh, 1480°C,
APS 7% porosity
3. THERMAL CONDUCTIVITY OF ZIRCONIA- ri4i
Zr02-6%Y2O3 1.1 [15]
BASED COATINGS APS
Although the thermal conductivity of thermal barrier Zr02-8%Y2O3 1. [15]
coatings has not been a major research and APS
development topic compared to work on Zr02-12%Y2O3 0.6 [15]
thermomechanical behaviour for instance, a number APS
of characterisations have been published, in particular Zr02-7%Y2O3 1.3 1.3 [16]
to provide designers with data for lifetime prediction EBPVD
models. Looking through some of the data available Zr02-7%Y2O3 0.8 0.8 [16]
and reported in table 2, several points are worth APS
noticing: Zr02-7%Y2O3 0.8 0.8 [16]
APS
- Thermal conductivity values for zirconia based Zr02-8%Y2O3 1.5 1.3 [17]
EBPVD as
coatings are relatively low (within the range 0.5 to 2 coated
W/m.K). This is significantly lower than values Zr02-8%Y2O3 1.9 1.5 [17]
determined on bulk zirconia [18, 19], typically in the EBPVD
2 to 4 W/m.K range, depending on the composition Zr02-7%Y2O, 2.9 . 2.7 [16]
and crystallographic structure. These values lie in the bulk
lower range of thermal conductivity for oxides.
Table 2 : Typical values of the thermal conductivity of
- Thermal conductivities of EBPVD coatings are zirconia-based coatings (EBPVD : electron beam
significantly higher than values measured for plasma physical vapour deposition, APS: air plasma spray)
sprayed coatings having the same composition. Composition percentages in wt%.
This can be explained by the structure of the coatings.
Plasma sprayed coatings exhibit a microcrack
network, with basically two crack populations : one
oriented perpendicular to the coating surface and
originating from cracking of individual splats during
1-4
(2680°C), and the tetragonal to monoclinic actually easier starting from an expression proposed by
transformation being accompanied by a large volume Slack [25] for non-metallic crystals and for
variation. The only available data seem to have been temperatures higher than the Debye temperature 6D
derived by Hasselman et al. [18] from experimental and derived form the same principles as the one
determination on various bulk zirconia-based system, followed by [22]:
knowing the relative proportion of the different
crystallographic forms and assuming that the thermal K = B.<M>.8.eD3/T.N2/3.y2
conductivity of a mixture can be expressed simply as a
function of the individual thermal conductivity. As where B is a constant, <M> is the average atomic mass
reported in table 3, the monoclinic form would in the crystal, 83 is the average volume occupied by an
present a conductivity in the range 4 to 5 W/m.K. atom and y is the Griineisen's constant (y = oc/Cv.%
with a the thermal expansion coefficient, % the
Another approach, based on purely theoretical
compressibility and Cv the specific heat). Bearing in
considerations, has been proposed by Klemens [22]
for calculating an «intrinsic thermal conductivity » of mind that 6D is inversely proportional to <M>1/2, it is
cubic zirconia, supposed to contain no stabiliser (and clear that low conductivity crystals correspond to :
therefore no associated vacancies). In a most general
way, the thermal conductivity K can be expressed as: - high atomic mass and weak interatomic bonding
(low 9D, which is a dominant term),
K=(l/3)JC'(a>).v.l((D)dffl - complex crystal structure (large N),
- high anharmonicity (large y).
where C'(ra) dm is the contribution of the specific heat
in the frequency range oo, oo+dro, v is the transverse According to these considerations, other ceramic
materials could then be envisaged as a basis for
phonon velocity and l(m) is the phonon mean free thermal barrier coatings, keeping in mind that other
path.
criterias have to be fulfilled (thermal expansion
coefficient, thermal stability, etc.). Among the
In the high temperature limit and in the case of no parameters which can have a large influence on K,
defect present, the mean free path of phonons is <M> is certainly a privileged one, and to a lesser
limited by three phonon interactions, owing to the
extent N. The Griineisen's constant (y) does not vary
cubic anharmonicities of the lattice forces, and this
expression leads to the following expression of the much from one material to another. An interesting
intrinsic thermal conductivity [23] : candidate, for example, is Ce02 (the atomic mass of
Ce02 is 40% higher than that of zirconia) a material
which has been tested recently by Schulz et al. [26]
K = (3/2)y2 . (n VVODD)!'-1 . N"2'3 and Maloney et al.[27]. However, if EBPVD Ce02-
9wt.%Y203 coatings present a low thermal
where y is the Grüneisen constant, |x is the shear conductivity (around half that of standard EBPVD
modulus, v is the transverse sound velocity, coD the Zr02-8wt.%Y203 coating [27]), they seem to exhibit a
Debye frequency, T the temperature and N the poor erosion resistance, a shortcoming which might be
number of atoms per primitive unit cell. Recently, alleviated though, if a multilayer architecture is
Klemens [22] applied such an approach to zirconia adopted, with a more erosion resistant external layer
and derived for the intrinsic thermal conductivity of added on top of the coating.
zirconia (i.e. that of an hypothetical «pure cubic
zirconia »): As far as compositions for the basis are concerned,
other ways are explored. Observing that only the
K(W/m.K) = 1700/T(K) thermal conductivity at high temperature is relevant,
and that for most ceramics, the thermal conductivity is
According to this expression, the intrinsic thermal a decreasing function of temperature (in principle T"1
conductivity of zirconia: law above), Padture and Klemens [28] are considering
- would amount to 5.7 W/m.K at room temperature, a ceramic materials with complex lattice cells. An
value surprisingly close to the values derived by [18] additional idea in this case is to try to develop an
for pure monoclinic zirconia, oxygen-impermeable dense ceramic layer to avoid or to
- decreases with a classical T1 variation, classical for minimise the oxidation of the bondcoat with formation
non-metallic materials [24], down to 1.3 W/m.K at of a thermally grown alumina at the ceramic/bondcoat
1300K. This last value appears excessively low, interface (locus of the final failure in the most
considering that dense PSZ polycrystals have a advanced systems). It is far from certain however, that
thermal conductivity higher than 2 W/m.K [18] and such a ceramic material can be deposited as a reliable
this represents a strict lower limit of pure zirconia. and dense coating and stay there in service,
(thermomechanical cycles). Moreover, the garnet
Even if absolute values cannot be derived from such compositions these authors propose still have a
simplified theoretical approaches, it is possible to relatively high thermal conductivity, even at 1000°C
infer trends concerning the effect of different (typically above 2.5 W/m.K).
parameters on the thermal conductivity. This is
1-6
In conclusion, possibilities exist to change the average atomic mass, and relative difference in ionic
chemistry of the basis material to decrease thermal size between solute and solvent. Diverse expressions
conductivity. A balance has, of course, to be carefully have been proposed to describe their influence on the
evaluated between the expected gains and the thermal conductivity [see for example 31, 32] and a
consequences on other properties (density, thermal critical study would be necessary to clarify the case of
stability, thermomechanical resistance, erosion zirconia-based materials.
resistance, etc.) to evaluate the final interest of these
potential new candidates. In the case of Zr02-Y203, the replacement of
zirconium ions (M = 91.22, r = 0.080 nm) by yttrium
4.2. Defects : ions (M = 88.90, r = 0.093 nm) should have a minor
An efficient way to decrease the thermal conductivity effect, their atomic mass in particular being relatively
of an insulator is to introduce structural defects which close. However, a significant effect may be expected
constitute obstacles to the propagation of phonons. from other substituents, in particular heavy rare earth
Several types of defect may be envisaged : vacancies, elements (at the expense though of the density, an
substitution ions, dislocations, interfaces (grain important drawback for mobile components). Very
boundaries), etc. few elements lighter than zirconium would satisfy
these criteria, with the exception though of (expensive)
The role of vacancies is particularly important in this scandium. As far as divalent ions are concerned, the
respect for zirconia-based materials. In fact, as pure choice is much more limited. Calcium or magnesium
zirconia cannot be of practical use, due to its cannot be envisaged in high temperature applications
undesirable transformation at around 1100°C, it is due to structural stability problems.
stabilised (partially or completely) by adding
heteroelements (most often Y, Mg, Ca) which Grain boundaries and interfaces can also scatter
stabilise the cubic (or the tetragonal phase). Thus the phonons and reduce the thermal conductivity. These
addition of 20wt% Y203 fully stabilises the cubic processes are most effective at low temperature for
form. Addition of 6 to 8wt. % Y203 permits to obtain ceramics in general [33], and only slightly at high
the so-called metastable t' phase (nontransformable temperature. It is to be noted that in TBCs, the grain
up to around 1200°C). The incorporation of size is much larger (typically 0.1 to 1 urn) than the
heteroelements is accompanied by the introduction of mean free path calculated according to :
vacancies to maintain the electrical neutrality of the
ionic lattice. Thus, introducing two yttrium ions is K~(l/3)C.v.l~0.5nm
accompanied by the introduction of one vacancy
according to the following reaction, using the and it is doubtful that, in these conditions,
nomenclature by Kroger and Vink : conventional grain boundaries may have a significant
effect.
Y203 ->2Yir'+Vo°°+3 0ox
As a conclusion, it can be said that the most efficient
where Y^' is an yttrium ion on a zirconium site ways to decrease the thermal conductivity are to
(single negative charge), V0°° is an oxygen vacancy introduce point defects :
twice positively charged and 00x is an oxygen ion on - vacancies associated with stabilising dopants. In
an oxygen site (no charge). zirconia-based coatings, up to 12 mol.% ROi 5 for a
trivalent dopant R, corresponding to the stabilised
Effectively, incorporating heteroelements such as Ca, cubic form, seem possible. However a fully stabilised
Mg or Y in zirconia lowers its thermal conductivity plasma sprayed material may not present a satisfactory
and a comparison of diffusivity measurements on thermomechanical behaviour.
zirconias with different stabilisers indicates that - for a maximum effect, the dopant should have a very
similar defect populations result in comparable different atomic mass than zirconium, provided that
thermal diffusivity [19]. It is to be noted though that elastic factor does not counterbalance the mass one.
this effect is observed only for stabiliser contents less
than a critical value [29]. For high vacancy contents, Finally, one could wonder what the minimal thermal
they cannot be considered as independent defects and conductivity would be with completely disordered
ordering may appear (as demonstrated by [30] in the materials, in other words what the minimal thermal
case of cubic zirconia stabilised with yttria and conductivity for a refractory material could be. Two
magnesia). It is therefore likely that phonons interact attempts have been proposed to calculate this minimal
less strongly with a population of ordered vacancies thermal conductivity.
than with a population of randomly distributed
defects. Slack [25] derives the theoretical minimum
conductivity assuming that the minimum mean free
Another source of efficient point defects acting as path for a phonon cannot be less than one phonon
phonon scatter centres is substitution ions. This acts wavelength. For pure zirconia, the value calculated is
essentially through two factors : AM/<M>, the ratio of 2.1 W/m.K. In a more recent approach, based on an
the mass difference between solute and solvent to the adaptation of the Einstein model for highly disordered
1-7
materials, Cahill et al. [34], calculate a minimum example) while keeping the computations tractable on
thermal conductivity of 1.3 for dense zirconia-yttria microcomputers.
systems. These values show that a gain may be
obtained by increasing the atomic disorder, but this This type of approach should offer guidelines for
gain will remain limited if the base material (in designers and coaters, provided the relationship
occurrence zirconia) remains the same. between process parameters and coating morphology is
known.
4.3. Optimising coating architecture :
a. Pores, crack network. b. Thin films.
The fact that the thermal conductivity of EBPVD, and It has been reported by several authors [38-42] that
in particular plasma sprayed coatings is much lower ceramic thin films may exhibit a lower thermal
than for bulk materials of the same composition is conductivity than the corresponding bulk materials.
explained by the presence of pores and microcracks Some examples are reported in table 4.
originating from the coating deposition process, as
described above. It is important to remark that up to Thermal Thermal
now, only few studies have been devoted to optimise Material conductivity conductivity of Ref.
the porosity distribution and crack network with of thin films bulk material
respect to thermal conductivity. In fact, most of the (W/m.K) (W/m.K)
PECVD Si02 0.05 (100 nm) [38]
effort in this field has been devoted instead to the 0.4 (1.4 urn)
influence of these features on the thermomechanical Glassy Si02 0.12 1.3 (amorphous) [391
resistance of the systems, with mainly empirical Ti02 1.5-5 7 [40]
approaches. (depending on
deposition
method)
As a first step, it is necessary to understand and Ti02 0.6 (evaporated) 5.5 - 7.6 (rutile) [41]
describe the relationship between thermal 0.5 (sputtered)
conductivity and morphological features (pores, Zr02(evaporated, 0.8 2-5 [42]
annealed In,
cracks,..). McPherson [35] was probably the first to 730°C)
propose a model for the microstructure of plasma A1N 16 70-220 [42]
sprayed coatings involving regions of good and poor (polycrystalline)
contact between lamellae to explain the low thermal
conductivity of coatings with respect to bulk Table 4 : Examples of thermal conductivity of thin
materials. The regions of poor contact act as thermal films compared with bulk values.
resistances. Considering an electrical analogy, and
assuming a very schematic geometry for the coating, Although this effect is clear, the values reported have
the ratio of thermal conductivity of the coating (Kc) to to be taken with care as measurements of thermal
that of solid material (K0) can be expressed as : conductivity of submicronic films are very delicate and
the error margins reported may be important.
Kc/Ko = 2.f.5/7r.a Different explanations have been proposed for this
effect. Most likely, this effect is due to a combination
where f is the fraction of «true contact», 8 is the of grain-boundary scattering (the grains are
lamellae thickness and a is the radius of the particularly small) and of atomic disorder during
individual contact areas. With a value of 0.2 for f and growth of thin films. Non equilibrium structures may
assimilating a with 5/2, this relation gives Kc/K0 be obtained by PVD or PECVD, in particular for
approximately 1/4, a right order of magnitude. It is to relatively low substrate temperatures (the diffusion of
be noted however that the predictive capacity of such species on the surface of the growing coating, and bulk
an approach is extremely limited due to the simplicity diffusion within the coating being too slow to reach an
of the geometry assumed and the experimental equilibrium state). In this respect, it is to be noted that
difficulty to determine the parameters involved. This the diffusion of cations within zirconia is quite slow.
approach has been employed by Moreau et al. [36] to With the diffusion data (D = D0 . exp-Q/RT with D0 =
describe the thermal diffusivity of plasma-sprayed 0.031 cm2/s and Q = 391 kJ/mol) published by [43],
tungsten based coatings, the geometrical quantities one can estimate that it takes about 8 years for a cation
being determined from image analysis. to travel 1 urn at 1100°C by solid state diffusion .
As described in Dorvaux et al. [37], the most To fully exploit this thin film effect, one could imagine
promising approach for developing a tool to describe thermal barrier coatings constituted of multiple
the thermal conductivity in relation with nanolayers of different materials, such as zirconia and
microstructure seems to be via a finite difference alumina as in [44], with the added advantage of
calculation carried out on digitised images of real incorporating boundary thermal resistances. The
coatings. This type of approach presents the materials of each layer and the process parameters
advantage of taking into account the complexity of the should have to be selected so that prolonged exposure
geometry involved (interconnected porosity for at service temperature would not result in structural
transformations towards equilibrium structures
1-8
(crystallisation of amorphous material for example), ■ designing new coating architectures (nanolayers).
activated by solid state diffusion and accompanied by
an undesirable increase in thermal conductivity. In It is felt that, up to now, several points have received
this respect, it is interesting to note that, as reported little attention and would be worth exploring in more
by Lee et al. [40], annealing a Hf02 thin film at details :
900°C enhances its thermal conductivity by only 20%. ■ characterising the optical properties of the real
It is to be noted also that the fabrication of such coatings, in order to evaluate the possible influence
coatings would certainly require serious modifications of radiation to the heat transfer through the
of techniques employed at present in production, or ceramic coating, and develop ways to decrease this
adapting new techniques to this application contribution.
(sputtering with multiple targets for instance).
■ assessing the thermal stability of the new systems;
with TBCs having low thermal conductivity, the
4.4. Concluding remarks
surface temperatures of the ceramic coating are
Several possibilities exist, in principle to lower the
likely to rise, compared to the temperatures reached
thermal conductivity of TBCs : incorporation of point
with present systems, and provoke a
defects (vacancies, substitution cations), optimising
microstructural evolution of the outmost region of
the microcrack network and, to a lesser extent,
this ceramic layer, through sintering, or phase
porosity distribution, and introducing structural
transformations, which in turn will increase the
disorder via the deposition process. It must be kept in
thermal conductivity.
mind though that, as exposed in the introduction,
TBCs must fulfil a set of requirements, and thermal
Finally, it is important to insist on the fact that
conductivity is only one of them. In particular, it
designing TBCs for engine airfoils requires a
would be highly questionable to develop TBCs that
multidisciplinary approach taking into account various
could not be deposited with processes such as plasma
aspects : materials (ceramic behaviour, oxidation of
spraying or EBPVD (or alternative processes
bondcoat, thermal conductivity, etc.), processes
currently under development). And the design of low
(ceramic morphology), life prediction modelling, blade
thermal conductivity coatings must take into account
design, quality control, commercial (patents),..
the requirements associated with the processes. For
Thermal conductivity is only one aspect of the design
instance, depositing materials having complex
and changing the ceramic composition to lower the
compositions, in particular with elements having
intrinsic conductivity, for example, is likely to have
widely differing vapour pressures, may prove to be too
important repercussions on the other properties. The
difficult to control in production with EBPVD
design of the final coating may have to fully exploit
processes ; similarly, a stack of numerous very thin
the flexibility and capabilities offered by present
layers presenting a highly disordered structure may be
coating processes, by developing multilayer coatings,
possible to fabricate by a PVD (or a CVD) process,
each layer (corresponding to a temperature range and a
but not by plasma spraying. Additional work is
particular environment in service) having a specific
needed however to study the thermal (meta)stability of function.
these defective structures.
Acknowledgements : Particularly fruitful discussions
5. CONCLUSION
with Odile Lavigne, during the preparation of this
Lowering the thermal conductivity of the ceramic
paper, are gratefully acknowledged.
coating of a TBC applied on airfoils in the hot stages
of a turbine represents an important challenge for
References
designing improved performance systems. This paper
has modestly indicated and listed a few paths to tackle
1. Miller R.A., «Thermal barrier coatings - History
this topic, some of which have already been explored :
and directions», in «Thermal Barrier Coating
■ changing the chemistry of the ceramic material in Workshop, NASA CP3312, 1995, pp 17-34.
order to incorporate point defects (vacancies, 2. Stecura S., « Optimisation of NiCrAl/Zr02-Y203
substitution cations) by playing both on the base Thermal barrier system, NASA TM 86905, 1985.
element (zirconia at the .moment) and on the 3. Bose S., DeMasi-Marcin J., «Thermal barrier
stabiliser, the main difficulty being to obtain coating experience in the gas turbine engine at
sufficiently stable phase structures and staying Pratt & Whitney», in NASA CP 3312, 1995, pp
within composition ranges such that the disorder 63-78.
is maximum as far as phonon scattering is 4. Meier S.M., Nissley D.M., Sheffler K.D., Cruse
concerned. T.A., « TBC life prediction model development -
■ optimising the coating morphological features phase II» - PWA Final report; NASA CR 189111,
(porosity, crack network) so that the heat transfer 1991.
perpendicular to the external surface is minimised. 5. Miller R.A., Garlick R.G., Smialek J.L., « Phase
This requires the development of modelling tools stability in plasma-sprayed zirconia-yttria», Am.
such as in the approach described in [36] and a Ceram. Soc. Bull., 62, 1983, pp 1355-1358.
good knowledge of the relationship between 6. Lelait L. (1991) : Etude microstructurale fine de
processes parameters and coating structures. revetements ceramiques de type barriere thermique
1-9
; incidence sur la resistance thermomecanique de stabilised (yttria) zirconia single crystals», J. Am.
ces revetements. These de l'Universite, de Paris Ceram. Soc. 71(4), 1988, 255-260.
Sud (Orsay). 20. Pawlowski L., Lombard D., Mahlia A., Martin C,
7. Alperine S., Lelait L., «Microstructural Fauchais P., « Thermal diffusivity of arc plasma
investigations of plasma-sprayed yttria partially sprayed zirconia coatings», High Temp. High
stabilized zirconia TBC, in relation with Press. 16, 1984, pp 347-359.
thermomechanical resistance and high 21. Taylor R.E., « An overview of the measurements of
temperature oxidation mechanisms», in Proc. thermophysical properties and some results of
37th ASME international Gas Turbine & molten superalloys and semiconductors», in
Aeroengine Congress (Köln) 92-GT-317 (1992). NASACP10121, 1993,ppl29-142.
8. Brindley W.J., «Properties of plasma sprayed 22. Klemens P.G., «Thermal conductivity of
bond coats», in TBC Workshop, March 27-29, zirconia», in «Thermal conductivity 23»,
1995, NASA CP3312, pp 189-202. Technomics, Lancaster PA (1996), pp 209-220.
9. Pint B.A., Wright LG., Lee Y.J., Zhang Y., 23. KlemensP.G., «Theory of the thermal
Pruzner K., Alexander K.B., « Substrate and bond conductivity of solids», in «Thermal
coat compositions : factors affecting alumina scale conductivity», vol.1, Acad. Press, 1969, 1-68.
adhesion », to be published in Materials Science 24. Kittel C, «Introduction to solid state physics »,
and Engineering. 1986, Wiley.
10. Schulz U., Fritscher K., Ratzer-Scheibe H.J., 25. Slack G.A., «The thermal conductivity of
Kaysser W.A., Peters M. , in Proc. 4th Internat. nonmetallic crystals», Solid State Physics 34,
Symp. on High temp. Corrosion & Protection (Les 1979, pp 1-71.
Embiez, France, 20-24 May 1996), to be 26. Schulz U., Fritscher, Peters M., « EBPVD Y203-
published. and Ce02/Y203-stabilized zirconia thermal barrier
11. Siegel R., Spuckler CM., «Analysis of thermal coatings - crystal habit and phase composition»,
radiation effects on temperatures in turbine engine Surf. Coat. Technol. 82, 1996, pp 259-269.
thermal barrier coatings», in Proc. TBC 27. Maloney M.J., Achter H.S., Barkalow B.K.,
Workshop 97, Cincinnati, OH, 19-21 May 1997, « Development of low thermal conductivity thermal
pp 17-26. barrier coatings». Communication presented in
12. DeMasi J.T., Ortiz M., Sheffler K.D., « Thermal Thermal Barrier Coating Workshop (Cincinnati,
barrier coating life prediction model development. may 1997).
Phase 1 », NASA CR 182230, 1989. 28. Padture N.P., Klemens P.G., «Low thermal
13. Miller RA, Leissler G.W., Jobe J.M. (1993) : conductivity in garnets», J. Am. Ceram. Soc.
Characterization and durability testing of plasma- 80(4), 1997, pp 1018-1020.
sprayed zirconia-yttria and hafnia-yttria thermal 29. Fournier D. et al., to be published (1997).
barrier coatings. Part I. NASA TP 3295. 30. Z.R. Dai, Z.L. Wang, Y.R. Chen, H.Z. Wu, W.X.
14. H.E. Eaton, J.R. Linsey, R.B. Dinwiddie (1994) : Lu (1996) : Local ordering of oxygen vacancies in
The effect of thermal aging on the thermal cubic zirconia stabilized with yttria and magnesia.
conductivity of plasma sprayed fully stabilized I. Electron diffuse scattering study. Phil. Mag.
zirconia. Thermal Conductivity 22, 289-300. 14b. 73(2) 415-430.29b. Z.R. Dai, Z.L. Wang, W.X. Lu
Dinwiddie R.B., Beecher S.C., Porter W.D., (1996) : Local ordering of oxygen vacancies in
Nagaraj B.A., « The effect of thermal aging on the cubic zirconia stabilized with yttria and magnesia.
thermal conductivity of plasma-sprayed and II. Determination of local ordering parameters of
EBPVD thermal barrier coatings», ASME 96- oxygen vacancies. Phil. Mag. 73(6) 1685-1698.
GT-982, pp 1-7. 31. Klemens P.G., «Theory of heat conduction in
15. Morrell P. , Taylor R., « Thermal diffusivity of nonstoichiometric oxides and carbides», High
thermal barrier coatings of Zr02 stabilized with Temp. High Press., 17, 1995, pp 41-45.
Y203 », High Temp. High Press. 17, 1985, pp 79- 32. Fukushima S., Ohmichi T., Handa M., « The effect
88. of rare earths on thermal conductivity of uranium,
16. Staniek G., Marci G, in Proc. Werkstoff- plutonium and their mixed oxide fuels», J. Less-
Kolloquium DLR, 10 dec. 1996, eds M. Peters et Common Metals, 121, 1986, pp 631-636.
al., pp 50-53. 33. Klemens P.G, « Phonon scattering and thermal
17. Portal R, Etude de la conductivite thermique de resistance due to grain boundaries», Internat.
couches minces de Zr02-Y203 deposees par Journal of Thermophysics 15(6) , 1994, ppl345-
EBPVD », Rapport de stage SNECMA,1997. 1351.
18. Hasselman D.P.H., Johnson L.F., Bentsen L.D., 34. Cahill D.G., Watson S.K., RO. Pohl R.O.,
Syed R, Lee H.L., Swain M.V., «Thermal «Lower limit to the thermal conductivity of
diffusivity and conductivity of dense disordered crystals », Phys. Rev. B 46(10), 1992,
polycrystalline Zr02 ceramics : a survey», Am. pp 6131-6140.
Ceram. Soc. bull. 66(5), 1987, 799-806. 35. McPherson R, «A model for the thermal
19Youngblood G.E., Rice R.W., Ingel R.P., conductivity of plasma-sprayed ceramic coatings»,
«Thermal diffusivity of partially and fully Thin Solid Films, 112, 1984, 89-95 .
1-10
The main criteria for power generation use of - a columnar structure under certain conditions
TBCs are the time and temperature effects on bond of evaporation,
coat and substrate, coating densification and changes - a fast tapering at the opening of cooling holes
in thermal or mechanical properties of the coatings. while plasma sprayed coatings have a greater
As underlined in reference [3] for advanced tendency to built up coating at that location,
land-based gas turbines, the goals that are to be - possibility to achieve much more smoother
achieved together with issues and challenges are surfaces,
summarized in fig. 1. - high adhesion/cohesion values (~ 70 MPa),
For a long period (more than 20 years) TBCs - low residual stress at room temperature
were essentially plasma-sprayed in air (APS) [4-7] or (~ 70 MPa),
controlled atmosphere [8]. Since a few years Physical - reduction of temperature of 38-66°C in stage of
Vapor Deposition assisted by Electron Beam (EB- High Pressure Turbine Blade,
PVD) TBCs, have been available on the market [9,10] - high number of cycles in engine testing : 180
especially for aircraft engines. The advantages for blades and 20000 for engine nozzle.
claimed relatively to APS TBCs are the following :
Reduced spallation
at the interface
V
Understanding damage
mechanisms
Oxidation,Fatigue, Creep, Phase
& Property changes, Hot
Improved corrosion, Erosion, Diffusion, etc
oxidation/corrosion Reliable life prediction
Protection
f
• PtAI, MCrAlY, Desulfurization
• VPS, HVOF, CVD etc.
However, it does not mean necessarily that APS spraying to achieve a much better reproducibility of
TBCs have no future. For example they have a the plasma spray process.
significant advantage in being able to accommodate
the large size components found in modern large 2. PLASMA SPRAYED COATING FORMATION
power generation machines. Each component
application of a TBC requires careful study and All the results presented in this paper concerned
review to determine whether APS or EB-PVD has to fused and crushed powders with 8 wt% Y203. They
be used [3]. Moreover, a better understanding of the are called in the following YSZ.
APS coating formation could result in coatings with
similar properties as those obtained by EB-PVD. The 2.1. Characteristic times and cooling rates
aim of this paper is first to review what is our Table 2 summarizes characteristic time values
knowledge in plasma sprayed coating formation. We for particles in the size range 22-70 urn sprayed with
will first were discuss the splat formation, their d.c. torches and impacting at velocities between 80
layering and the resulting coating adhesion/cohesion and 260 m/s on smooth surfaces (Ra < 0.1 urn) [11-
as well as residual stresses. We will second underline 14].
which parameters have to be controlled during
3-3
Flattening time Starting of nucleation Solidification time Time between two Time between two
if (us) tn(Hs) ts (us) impact events at the successive passes
same location t; (|is) U (Us)
1-3 0.3-2.5 <10 500-5000 106-107
Table 2 : Characteristic times in d.c. plasma spraying.
tf = 2.dp.Reo.2/(3.vp) (2)
Anon -n
3500 -
* 3000 - \Splat surface
§ 2500
g 2000 -
Q) Vsplat bottom
a 1500 -
\ "■^SSSte*™-,,,
o 1000 -
Substrate surface 500 -f
500 - Substrate surface
0
n
50 100
() 5 10
Time (ps) Time (|JS)
2
a) R,,, = 10 s K.m /W b) R^ = 5.HM K.m2/W
Fig. 2 : Splat surface and bottom as well as interface temperature evolutions for two values of R^,
(1-um thick Zirconia splat on a smooth stainless steel substrate at 600 K).
3-4
hot substrate, the mean Shape Factor SF is close to The latter exhibit on their whole surface, a
one, both in the spray cone central part and its homogeneous microstructure with a columnar growth
periphery even if the mean splat diameter D is perpendicular to the surface [20] (see fig. 5.a) except
slightly smaller in the spray cone periphery, as it in the splat rim (see fig. 5.b) where the columnar
could be expected. SF is defined as 4 n S/P2 where structure is parallel to the substrate. Indeed, due to
P is the splat perimeter and S its surface area (SF = 1 splat curling [24] the rim is not in contact with the hot
for a disk). For cold substrates D is by far smaller substrate and cooling occurs through the already
due to material splashing and SF ~ 0.67 against 0.94. solidified columnar structure of the central part of the
The dispersion is also higher. However, if the splat resulting in a columnar structure parallel to the
substrate is oxidized (due to a too long preheating) substrate. As the cooling is much slower, the size of
splats are again extensively fingered. The oxidation the columns is bigger. On cold substrates (Tph < Tc)
stage depends strongly on the preheating time with the only a very few AFM measurements could be done,
plasma jet and the substrate material [20,22]. Splat the splats being pulled off very easily. In this case, the
cooling rates, measured by fast pyrometry [12-14], are central part of the splat consists of a few islands
3-6 times higher with the hot substrate. For the disk where the contact is good (columnar structure),
shape splats obtained at 573 K, the splat cooling however with bigger columns than on hot substrates,
theoretical curve [23] of a 1-D model has been surrounded by rims where the crystal growth is
adjusted to the experimental curve by varying the parallel to the substrate [20]. The "fingers" exhibit the
thermal contact resistance. Values of R,h ~ 10-7 - same morphology as that of the rims. As a general
10"8 m2.K/W were obtained corresponding to an rule, the size of the columns depends very strongly on
almost perfect contact as confirmed by the crack the cooling rate as illustrated in Table 4 for two
network due to quenching stress relaxation observed different substrate materials and spraying equipments.
in fig. 4.b. On cold substrates, the contact resistance With the RF plasma torch, the particle impact velo-
Rth calculated assuming that the splats have a disk- cities are lower than 80 m/s resulting in splat thick-
shape, were found to be of the order of 10"6 m2K/W. nesses es in the range 2-2.5 urn against 0.9-1.2 um
This poor contact resistance is confirmed by the pul- with d.c. plasma torches. These thicker splats cool
ling off of splats by the tip of a perthometer. This does down slower (in a ratio 6-10) than the thinner splats.
not occur when splats are collected on hot substrates.
3-6
Fig, 5 : AFM view of an YSZ plasma sprayed with a d.c. torch (vp = 210 m/s, Tp = 3400 K)
(a) splat central part (note the cracks) (b) splat rim.
Plasma torch d.c. torch with a 7 mm i.d. nozzle r.f. torch with a 50 mm i.d.
substrate material stainless steel 304 YSZ stainless steel 304L YSZ
splat thickness es (|xm) 0.95 2.2
CR cooling rate (K/us) 645 115 105 16
ts solidification time (us) 0.74 1.23 4.54 7.6
column width (nm) 125 250 320 >400
Table 4 : YSZ column width with respect to cooling conditions [19].
The decreasing in A/C values when the preheating 2.3.2. Thick passes
time increases is explained by the time-growth and In order to limit oxidation problems 300 urn
time-composition change of the oxide layer at the thick coatings sprayed on In 100 substrates were
304L surface, as shown by recent studies [22] : the preheated to different temperatures and the pass
plasma plume contains atomic oxygen and oxidation thickness was varied from 25 to 300 um/pass. With 25
is somewhat different from that obtained by urn passes for preheating temperature Tph over 873 K,
preheating substrates in a furnace. Similar results are aq is more than doubled : 20 MPa against 8 MPa. For
obtained on FT25, XC38 and In 100 substrates for Tp ~ 973 K, the fractured cross section of coatings
different temperatures and preheating times. The exhibits a columnar structure growing through the
coating Young's modulii obtained by comparing the whole coating (see fig. 8.a). The Vickers hardness
deflections of XC38 beams measured during spraying with a 5N load is increased to 1200 ± 80 MPa against
[25,28] with values calculated using a 1-D code based 800 ± 70 MPa with the 3-um thick pass (with Tp =
on the equations proposed by Tsui and Clyne [30], are 773 K). Unfortunately when spraying on In 100
in the range 20-30 GPa. However, the precision is substrates preheated at 973 K, the coating detached
poor (± 30%) with such low values. upon cooling probably due to oxidation problems.
3-8
a) b)
Fig. 8 : Fractured cross section of an YSZ coating :
a) at x = 0.4 mm from substrate and sprayed on an already deposited YSZ coating (300 urn on In 100)
preheated at 973 K ep = 25 um, b) sprayed on In 100 preheated at 973 K, x = 0.005 mm.
resistance. For the lower limit of the particle results in an increase of the jet instability and a
distribution the injection problems are critical. Most decrease in jet length coupled with a decrease in
of small particles, especially those below 10 urn, coating quality. When a tungsten insert is used this
hardly penetrate" within the plasma jet, travel in its erosion problem may become even more critical when
periphery and are sucked down farther downstream tungsten particles are ejected and imbedded within
resulting in defects in the coating when they stick to coatings. Therefore, it is important to monitor the
it, which is especially the case when the coating anode erosion by the statistical analysis of the voltage
surface is kept over 773 K during spraying. signal [43] or its power spectrum [44], and determine
when this wear will induce poor quality coating, long
3.1.2. Powder injection in the plasma-jet before its failure.
Internal injection allows a better penetration of This type of monitoring makes it possible to
the particles within the plasma jet compared to check rapidly a defective gun assembly resulting for
external injection [37]. However with small particles, example in a poor electrode centering ; a non
requesting a higher carrier gas flow rate, the plasma symmetrical arc attachment influencing deeply the
flow perturbation is more important with internal restrike mode frequency and the voltage signal
injection [38]. The position of the injection port and evolution.
its distance z from the torch axis is very critical. A
variation of z of a few tenths of mm modifies deeply 3.2.2. Real time control
the particles mean trajectory. Thus, any change in the The outputs that have to be set and maintained
injector port location and its wear have to be checked are the particles mean temperature, velocity, trajectory
very carefully and regularly. It is also the case of the and the substrate and then coating surface temperature
internal wear of the injector, a rough internal surface before (preheating), during and after (cooling down)
modifying deeply the powder distribution at the spraying. Among these parameters the most
injector exit compared to that obtained with a new important, at least to our opinion, are the two last
injector with a smooth internal surface. More ones.
precisely (see section 3.2.2) the trajectory distribution Substrate and coating temperature control
of particles within the plasma jet has to be followed depends strongly on the size, shape and thickness of
continuously because it controls acceleration and the parts to be coated, the relative movements
heating of the particles. between the torch and the substrate and powder feed
rate, linked to the pass thickness, and cooling means.
3.2. Plasma spray process If the surface temperature can be followed by IR
3.2.1. Wear of the electrodes pyrometry, the close-loop controller has to be
The arc behavior at the anode and cathode is designed according to substrate and coating heating
quite different. Thoriated tungsten cathode tip wear is and cooling, if possible by using simplified 1-D or
rather fast (~ 1 h). It is due to the fast diffusion of 2-D computer codes easily adaptable to different
thoria which increases drastically the cathode tip geometries and pass thicknesses.
temperature (by about 1000 K) and induces a neck The mean trajectory of the particles within the
erosion below the tip which is blown up after a while plasma jet is very sensitive to the powder injection
[39,40]. The transition of a sharp cathode tip to a conditions which have to be adapted to the spraying
rounded one, which extremity becomes slowly and conditions [37]. With a 1-D photodiode array, it is
progressively wider, induces, within the first working possible to follow continuously the hot particles
hour, a reduction of the flow velocity of about 20% distribution and, when coupling it to a laser sheet, to
and a reduction of the arc voltage of 4 to 6% [41]. The have simultaneously the trajectory distribution of the
widening of the cathode tip reduces progressively the cold particles. The hot particles distribution is very
voltage by 5-8% after 20 h, and the jet velocity by sensitive to any variation of the injection parameters
about 10%. [46]. For example with Zr02 particles (- 45 + 22 urn)
The arc root fluctuations (restrike mode) [41,42] a variation of the powder feed rate of 50 g can be seen
which condition the arc root life time and thus the and the signal intensity is multiplied by 1.9 when the
anode erosion, evolve with operation time and powder flow rate is doubled (from 1 kg/h to 2 kg/h).
number of torch ignitions. This erosion is distributed Only 1.9 and not 2 due to the starting of the load
in an area which length is about 1.5 x d (d being the effect. Any change of the optimum carrier gas flow
internal diameter of the anode). The depth of the rate m°cg modifying the trajectory is immediately
eroded area as well as the life time of the arc root detected. For example a change of 10% relatively to
increase with the operation time. Anode erosion
3-10
its optimum value induces a 10% change of the signal control the induced residual stresses, temperature of
amplitude and position. substrate and coating have to be monitored
Of course, it is also possible to control on line continuously by adjusting the pass thickness and the
the surface temperature of the particles prior to their cooling devices to the coated part dimensions with the
impact [47] or together the temperatures, velocities help of simplified 1-D or 2-D computer codes.
and trajectories of the particles [48]. These robust At last, to improve coatings reproducibility,
easy-to-use optical sensors, now commercialized, strongly linked to the molten state and velocity of
allow us, after positioning the torch in front of the impacting particles, a rather simple set-up has been
sensor head for 1 minute, to check in industrial developed giving in real time the trajectory
conditions the optimization of the spray parameters distribution of hot particles within the plasma jet and
and their time evolution. For this last point after a allowing to check any change in powder flow rate,
certain spraying time the torch has to be positioned carrier gas injection flow rate, injector wear...
again in front of the sensor head. Such devices make Compared to commercial set-ups following particle
them possible to determine when the load effect velocity, surface temperature and flux density it is by
becomes important, measure the time necessary far much simpler and seems to be as efficient. The
before reaching the stability of particle spraying erosion of the torch nozzle has also to be followed by
conditions, increase the reproducibility of the analyzing the voltage fluctuations which vary
spraying process. drastically when conditions are reached where coating
quality is affected and allows also the detection of
Conclusion defective gun assembly.
rn
CD
3
Figure 4 Time evolution of YSZ particle temperature
00 2600 - for 3 guns used in production : (T) before and after
(A) deposition on different parts.
t
L_
CO
Q.
E
CD 2400 - The number of points differs for each gun as the guns were
i
used for spraying different powders and only the
i measurements on YSZ particles are shown. It is worth
??nn i
i I i noting that there is a difference of typically 20°C between
temperature measured before and after deposition on each
50 100 150 200 250
part. Even if this difference is small, it is significant,
Time (sec) indicating that the conditions of the particles have slightly
3000 changed during spraying.
The guns 1,2,3 have been changed after 46, 58 and 48 hours
9 2800 of spraying, respectively. For gun 3, the reason of change is
not imposed by consideration of stability of the process but
CD
is due to a problem of waterproofness in the torch. With this
-I—«
öS
CO 2600 - gun, the process was considered stable according to the
L_
CD PWC criteria and also according to the measurements on the
Q. particles (Fig. 4). The situation is different for guns 1 and 2
E as the temperature and also velocity of the particles (not
CD 2400
shown) vary significantly during the use of these guns in
production. It appears that these variations are often related
2200 to variations of voltage between electrodes which is one of
the parameters recorded by the operators. In fact, when the
4200 4250 4300 4350 4400 4450 voltage variations are too large, the usual procedure is to
stop spraying as it is the case for gun 1 after 46 hours of
Time (sec) operation. Consequently, the two triangles at 46 hours
correspond to a beginning of spraying but no deposition was
Figure 3 Example of particle temperature evolution made as the voltage variation was too large. Even if gun 1
during coating production with YSZ particles, (a) before was used for coating deposition during 46 hours, Fig. 4
deposition, (b) after deposition.
4-4
shows that the spraying conditions of the particles have metallic (NiCoCrW) powders, 27 urn average diameter.
constantly and significantly changed after 35 hours of Figure 5 shows the evolution of the particle jet intensity
operation. It is worth noting that the process was controlled before and after deposition. In this case, the process didn't
by keeping the input power constant. So, as the electrodes reach a stable operation point even after 3 minutes of
become worn, the voltage tends to decrease and the spraying when the torch normally moves to spray the part to
controller reacts by increasing the current in the gun. In the be coated. Such a long time to obtain stable spray
case of gun 1, it clearly appears that such a correction was conditions is likely related to the feeding conditions of the
too large leading to a significant increase of the temperature fine powder into the plasma.
of the sprayed particles. A corresponding increase in the
velocity of the particles was also noted.
The values for the average particle parameters for each gun
are given in Table 1. The values for each parameter are
average values obtained during the first hours of use of each
gun, that means during the period of time they give
relatively stable conditions for the particles. From Table 1, n 1 r
it appears that the different particle parameters are very
similar for guns 1 and 3 except the difference of 50°C in 600 650 700 750 800 850 900
temperature. The differences with gun 2 are larger as the
temperature difference reaches 82°C. From Table 1, it is Time (sec)
clear that in spite of the control established on the process,
there are still significant differences in the particle Figure 5 Time evolution of the particle jet intensity
characteristics. These differences can be explained by a collected by the CCD camera.
variation in the way the particles are injected in the plasma
jet. This is especially important if the plasma characteristics
change as we can expect if the plasma gun is changed.
4.2 Effect of the wear of the electrodes
Consequently, a control of the particle state appears more
In order to better understand the effect of the wear of the
sensitive to variations in the process conditions than the
electrodes during spraying, a series of measurements of the
current approach based on the control of plasma power.
particle characteristics was carried out in the laboratory
[21]. A torch was mounted with a new set of electrodes
It was observed that, in the case of YSZ particles, the
and the state of the particles was monitored as a function of
spraying conditions are relatively stable during the stability
the spraying time. At regular intervals, YSZ coatings about
periods (Fig. 3) even if the temperature can slightly vary
500 urn thick were sprayed.
during deposition. This is not true in the case of finer
4.2.1 Two-Dimensional Profiles
Four 2-D profiles were measured after 1, 8, 22 and 37
Table 1 Particle parameters measured during the first hours of spraying. Mean temperature, velocity and
hours of operation for 3 guns used in production. diameter of the sprayed particles were monitored during the
scans. These profiles were collected by moving the sensor
GUN1 GUN 2 GUN 3 head with a X-Y unit, scanning in this way the
measurement volume of the DPV2000 across the particle
Temperature (°C) 2687 2769 2734
jet at 65 mm from the gun. The number of detected events
Velocity (m/s) 156 161 152 which is related to the particle flux at the measurement
Diameter (urn) 55 67 60 volume location was also counted for a fixed period of time
Flux (part/min.) 1065 1290 1123 (25 ms). This number can't be taken as the true particle flux
Jet position (mm) -0.9 0.0 -0.5 but is qualitatively related to it in the sense that both follow
Jet width (mm) 4.7 4.6 4.7 the same trends.
Jet intensity (a.u.) 1.0 1.3 1.1
4-5
-20 -10 0 10 20
4.2.2 Coatings and Deposition Efficiency
Coatings were sprayed after 3, 10, 22 and 37 hours of X(mm)
spraying in order to evaluate the effects of the sprayed
particle state on the microstructure of the deposited 10 -
material through the long-term experiment. Image analysis 2700 S 2900 27l 0
>
showed that the total porosity was the same for all
coatings, at about 16%. E ^
After three hours of spraying, the deposition efficiency
was 55% while at the end of the experiment (37 hours of >• -10
spraying) it fell to 41% reflecting important changes in the
state of the sprayed particles. Figure 8 shows the -20
microstructure of the coatings deposited after 3, 10 and 37 I After 22 hours of spraying
i i i i i i i i i i i i i i i ii i i
hours of spraying. As seen in the figure, a clear evolution
of the coating microstructure took place with increasing -20 -10 0 10 20
spraying time. At the beginning, after 3 hours of spraying,
X(mm)
the deposited coating contained a relatively large amount
of cracking. The cracking pattern, though still present, is
not so apparent in the coating sprayed after 10 hours. The 10 ,2700-
mean particle temperature was nearly the same for these
two coatings whereas the particle velocity differed by
nearly 30 m/s. Coatings sprayed after 22 (not shown in E °
Figure 8) and 37 hours were very similar in their E,
microstructure. However, in these coatings cracking
>
patterns like those observed after 3 and 10 hours of -10 "I 27U0
spraying did not form. These microstructure differences
can be attributed to the different particle state before their -20 :
impact on the substrate. Indeed, in a previous After 37 hours of spraying
i i i i l i i i i l i i i i l i i i i
communication [22] it has been demonstrated with a
similar ceramic powder that a mean particle temperature -20 -10 0 10 20
variation of 200°C, such as the one measured in the present X (mm)
experiment between 3 and 37 hours of spraying, can result
in a more than 10% thermal diffusivity variation as well as
significant microstructure changes. It should be reminded Figure 6 Temperature distribution (°C) of the
that between 3 and 37 hours of spraying the gun power sprayed particles after different spraying times. The
gun axis corresponds to X = 0 mm and Y = 0 mm.
4-6
>" -10
-20 -
After 1 hour of spraying
i i i i i i i i i | i i i i | i i i i
-20 -10 0 10 20
X (mm)
10
0
E
E after 10 hours of spraying
> -10
-20
After 8 hours of spraying
i i i I i i i i I i i i i I i i i i
-20 -10 0 10 20
X(mm)
10 :
:
E
E
0
'ffi^s >*
> -10
-20 -10 0 10 20
X(mm)
10
E °
E.
> -10
-20
THIpI I I | I I I ! |I II I
varied by only 1 kW. It would have been difficult, if not with higher hydrogen proportion and gas flow rate is
impossible, to predict any microstructure changes from this shown in Fig. 9. It can be seen from this figure that the
power variation. Present results show that coatings cracking pattern is rather similar to that of the coating
prepared after different spraying times with the same input sprayed after 3 hours, as seen in Fig. 8.
spray parameters can result in materials having significant
differences in their microstructures and properties. So after a long spraying time (about 40 hours) by changing
some input parameters it was possible to recover the
4.2.3 Attempts to Recover Initial Particle State sprayed particle parameters obtained after three hours of
After the forty-hour experiment, the electrodes were spraying and to recover the same type of coating structure
disassembled from the plasma gun. They were mounted observed before significant electrode wear took place.
again several months later when attempts were made to try
to recover the initial particle parameters. In a first attempt, 4.3 Controlling the particle state and coating structure
the gun current was set to 900 A. However, no significant In order to implement a feedback loop to control the state of
changes in the mean particle temperature and velocity the in-flight particles, it is necessary to know how the input
were observed neither when gun current was increased spray conditions influence the state of the particles. It is also
from 800 A to 900 A and to 1000 A, nor when the powder necessary to know how the temperature and velocity of the
carrier gas flow rate was changed. Nevertheless, when impinging particles affect the coating structure. For TBCs
hydrogen proportion was changed to 5% (as compared to coatings, one key characteristic is the crack network within
3.9% vol. during the forty-hour experiment) the mean the coating. Indeed, it seems advantageous to maximize the
particle temperature increased by about 150°C (from number of cracks perpendicular to the substrate while
3000°C to 3150°C) whereas the mean particle velocity minimizing the number of parallel cracks [24].
increased by more than 15% (from 148 m/s to 175 m/s).
Moreover when the total gas flow rate was increased from An example of how the arc current and total gas flow rate
52.0 1/min to 60.6 1/min, keeping the hydrogen proportion influence the state of YSZ particles is shown in Fig. 10 [22].
to 5% vol., the mean particle temperature increased from The powder used in this study was a fused and crushed YSZ
3150°C to 3200°C while the mean particle velocity powder (-45 +11 urn). After setting the plasma current and
increased from 175 m/s to 190 m/s.). These velocity and total gas flow rate, the powder carrier gas flow rate was
temperature are similar to those measured after two hours adjusted to position the spray jet axis at a constant angle of
of spraying. The measured gun power was 35 kW, which is about 3 degrees below the torch axis. As shown in the
about 5 kW higher than its value after two hours of process control plot of Fig. 10, an increase in plasma current
spraying. increases both the temperature and velocity of the particles
whereas an increase in the gas flow rate increases their
Using these spraying conditions, 2-D profiles of the velocity but decreases their temperature. Subsequently,
particle jet were measured. They are rather similar to those conditions for spraying particles at an exact temperature and
measured after eight hours of spraying (see Figs. 6 and 7). velocity may be obtained by carefully regulating both the
The numbers of detected events are similar to those arc current and gas flow rates. In this way it is possible to
observed after 8 hours of spraying (Fig. 7). A coating as regulate the spray process in order to spray coatings in a
well as a deposition efficiency test were performed using reproducible manner.
the new input parameters. The deposition efficiency was
measured to be 53% which is essentially the same as what After having established the data shown in Fig. 10, it is
had been measured after three hours of spraying (55%) and possible to optimize the coating structure by varying
much higher than the 41% value measured after 37 hours independently the temperature and velocity of the sprayed
of spraying. The microstructure of the coating prepared
3050
3000
current / /
2950 increase /? *^^ /
/
O
v
S- 2900 /gas flow /
/ increase/
f 2850 - / /
| 2800
h-
2750
2700
-Y
2650 i i i i
3050 3000
#2 II
• 2995 °C Vparticle=250 m/s
3000 - #7 • : «#1 2700 " \ Tparticle=2850 °C
2950
250 m/s 2400
ü
2900
2840 °C .
2100 V •/
I 2850 - * #3
CD
V■
#4 *
CD K 1800 - \•
|2800 */ /
200 m/s
1500 -
2750 •/ ■ /
■^ #5 1200
2700 • ■ ■
#6 900 -Tparticle=3014°C
2650 i i i i i i i i i
140 160 180 200 220 240 260 280 300 320 20 40 60 80
Velocity (m/s) Angle (degrees)
Figure 11 Average temperature and velocity of the Figure 12 Angular crack distribution in coatings 4 and
sprayed particles for coatings sprayed in the study of 5 sprayed with particles at 250 m/s and (•)3014°C
Prystay et al. and (■) 2850 °C, respectively.
particles instead of varying the input spray parameters. For the sprayed coatings increase, the process control becomes a
example, in the work of Prystay et al. [22], the coatings key issue in order to produce, day after day, high quality
were sprayed in the 7 conditions shown in Fig. 11. The coatings in a reproducible manner. Monitoring the
velocity and temperature effects were decoupled by temperature, velocity and trajectory of in-flight particles
selecting operating parameters to generate data points that before their impact on the substrate appears as an efficient
fall along a staircase pattern. For example, the influence of way to characterize the process and to perform feedback on
particle velocity at about 2995 °C is seen by comparing the input spray parameters. Moreover, it is also an efficient
coatings 1 and 7. (The average temperature between the tool to optimize the spraying conditions by varying
particle temperatures of samples 1 and 7 is used. Other parameters that have the most crucial effects on the coating
values shown on the graph were determined in the same properties.
manner.) The effect of velocity at particle temperatures of
2840 °C is seen by comparing samples 3 and 4. The results 5. ACKNOWLEDGEMENTS
of changing particle temperature for fixed particle velocities The authors are grateful to many people that contributed in
are obtained by comparing sample 4 with 5 (200 m/s) and this work. In particular, we would like to thank M.
sample 3 with 2 (250 m/s). Lamontagne and S. Belanger from NRC Canada for their
help in spraying coatings and developing the optical
In this work the angular crack distribution was characterized diagnosis system. Moreover, we are grateful to F. Richard
by image analysis using a dedicated algorithm developed in for taking charge of the project at Pratt & Whitney Canada.
our laboratory. This algorithm treats each pore or crack as
an individual feature of equal weighting and no results 6. REFERENCES
regarding their width are included. The orientation of each
feature is determined and the total length of the features (in 1. Smith, R.W., "Plasma Spray Deposition: A Need for
pixels) as a function of their orientation is calculated. Direct Process Control", in "Review of progress in
quantitative nondestructive evaluation," Vol 8B, D. O.
An example of results is given in Fig. 12 where the crack Thomson and D. E. Chimenti (Eds), Plenum Press,
distribution of coatings 4 and 5 are compared. Both coatings New York, 1989, pp 1169-1183.
were sprayed at particle velocity of 250 m/s but with 2. Fincke, J.R., "Diagnostics and Sensor Development fir
different particle temperatures. In this case, the increase in Thermal Spray Technologies", in "Thermal Spray:
temperature has the effect of sharply decreasing the number International Advances in Coatings Technology," Proc.
of horizontal cracks (0 degree) and increasing the vertical Inter. Thermal Spray Conference, Orlando, ASM
cracking (at 90 degrees). In fact, coating sprayed with International, Materials Park, OH (1992) pp 1-9.
particles at 250 m/s and 3014 °C exhibits only a weak 3. McPherson, R., "The Relationship Between the
preferred angular crack orientation. In addition, although Mechanism of Formation, Microstructure and
this coating is not the densest in this study, it is the coating Properties of Plasma-Sprayed Coatings", Thin Solid
with the highest thermal diffusivity [22]. This suggests that Films, 83, 1981, pp 297-310.
the significant reduction in horizontal cracking decreased 4. Mishin, J., Vardelle, M., Lesinski, J. and Fauchais, P.,
the thermal insulating properties of the coating. "Two-colour Pyrometer for the Statistical
Measurement of the Surface Temperature of Particles
5. CONCLUSION under Thermal Plasma Conditions", J. Phys. E:Sci.
Plasma-sprayed coatings are used in a wide range of new Instrum., 20, 1987, pp 620-625.
and demanding applications. As the quality requirements of
4-9
5. Vardelle, M., Vardelle, A., Fauchais, P., Boulos, M.I., On-line Diagnostics of Particles under Thermal
"Plasma-Particle Momentum and Heat Transfer: Spraying Conditions", in "Advances in Powder
Modelling and Measurements", AIChE Journal, 29, Metallurgy & Particulate Materials", Vol.6, Advanced
1983, pp. 236-243. Processing Techniques, C. Lall and A. J. Neupaver
6. Sakuta, T., Ohtsuchi, T., Sakai, K., and Takashima, T., (Eds.), Proc. Inter. Conf. on Powder Metallurgy and
"Measurement of In-Flight Particle Parameters in Particle Materials, Toronto, Canada, Metal Powder
Thermal Plasma Spaying Process by the Light Pulse Industries Federation, 1994, pp 199-210.
Analysis Method", in "Proc. Jpn. Symp. Plasma 19. Tecnar Automation Ltee, 3502, First Street, St-Hubert,
Chem.", Vol. 4, 1991, pp. 175-180. Quebec, Canada, J3Y 5K2.
7. Fincke, J.R., Swank, W.D., Jeffery, C.L. and Mancuso, 20. Gougeon, P., Moreau, C. and Richard, F., "On-line
CA., "Simultaneous Measurement of Particle Size, Control of Plasma Sprayed Particles in the Aerospace
Velocity and Temperature", Meas. Sei. Technol. 4, Industry", in "Advances in Thermal Spray Science and
1993, pp 559-565. Technology", C. C. Bemdt and S. Sampath (Eds.),
8. Guselnikov, S.M., Zavarzin, A.G., Lyagushkin, V.P., Proc. National Thermal Spray Conference, Houston,
Mikhalchenko, A.A. and Solonenko, O.P., "Laser Texas, ASM International, 1995, pp 149-155.
Optical Diagnostics Complex For Investigation of 21. Leblanc, L., Gougeon, P. and Moreau, C,
High-Temperature Heterogeneous Jets", in "Plasma "Investigation of the Long-Term Stability of Plasma
Jets", Solonenko and Fedorchenko (Eds.), VSP, 1990, Spraying by Monitoring Characteristics of the Sprayed
pp 163-170. Particles", in "Proc. United Thermal Spray Conf,
9. Coulombe, S., Boulos, M.I. and Sakuta, T., Indianapolis, Indiana, ASM Inter., Sept. 1997, to be
"Simultaneous Particle Surface Temperature and published.
Velocity Measurements under Plasma Conditions", 22. Prystay, M., Gougeon, P. and Moreau, C, "Correlation
Meas. Sei, Technol., 6, 1995, pp. 383-390. Between Particle Temperature and Velocity and the
10. Kudora, S., Fujimori, H., Fukushima, T. and Kitahara, Structure of Plasma Sprayed Zirconia Coatings", in
S., "Measurement of Temperature and Velocity of "Thermal Spray: Practical Solutions for Engineering
Thermally Sprayed Particles using Thermal Radiation", Problems", C.C. Berndt (Ed.), Proc. National Thermal
Trans. Jpn. Welding Soc, 22, 2, 1991, pp. 10-17. Spray Conf, Cincinnati, Ohio, ASM Inter., 1996, pp
11. Frind, G., Goody, C.P. and Prescott, L.E., 517-523.
"Measurement of Particle Velocity in Two Low- 23. Lethtinen, T., Knuuttila, J., Vattulainen, J., Mäntylä, T.
Pressure Plasma Jets", in "6th Inter. Symposium on and Hernberg, R., "Correlations Between In-flight
Plasma Chemistry (ISPC-8)", Montreal, 1983, pp 120- Particle Concentrations and Coating Properties in
125. Atmospheric Plasma Spraying on Alumina", in
12. Mayr, W. and Henne, R., "Investigation of a VPS "Thermal Spray: Practical Solutions for Engineering
Burner with Laval Nozzle using an Automated Laser Problems", C.C. Berndt (Ed.), Proc. National Thermal
Doppler Measuring System", in "First Plasma-Technik Spray Conf, Cincinnati, Ohio, ASM Inter., 1996, pp
Symposium", Vol. 1, 1988, pp 87-97. 525-530.
13. Söelter, H.-J., Müller, U. and Lugscheider, E., "High- 24. Bengtsson, P., Jonannesson, T. and Wigren, J., "Crack
Speed Temperature Measurement for On-line Process Structures in Plasma Sprayed Thermal Barrier
Control and Quality Assurance during Plasma Coatings as a Function of Deposition Temperature", in
Spraying" Powder Metal. International, 24, 1992, pp. "Thermal Spraying - Current Status and Future
169-174. Trends", Proc. Inter. Thermal Spray Conf, Kobe,
14. Hencken, K.R., Tichenor, D.A. and Wang, C.F., Japan, A. Ohmori (Ed.), 1995, pp 347-351.
"Optical Double-Slit Particle Measuring System", U.S.
Patent 4,441,816,1984.
15. Swank, W.D., Fincke, J.R. and Haggard, D.C., "A
Particle Temperature Sensor for Monitoring and
Control of the Thermal Spray Process", in "Advances
in Thermal Spray Science and Technology", Proc. 8th
National Thermal Spray Conf, ASM International,
1995, pp. 111-116.
16. Li, K.-L, Vardelle, M. and Fauchais, P., "Comparison
Between Single and Double Flow Injectors in the
Plasma Spraying Process", in "Advances in Thermal
Spray Science and Technology", Proc. 8,h National
Thermal Spray Conf., ASM International, 1995, pp.
45-50. ;
17. Moreau, C, Gougeon, P., Lamontagne, M., Lacasse,
V. Vaudreuil, G. and Cielo, P., "On-line Control of the
Plasma Spraying Process by Monitoring the
Temperature, Velocity, and Trajectory of In-flight
Particles", in "Thermal Spray Industrial Applications",
C. C. Berndt and S. Sampath (Eds.), Proc. National
Thermal Spray Conference, Boston, ASM Intern.,
Boston, 1994, pp 431-437.
18. Gougeon, P., Moreau, C, Lacasse, V., Lamontagne,
M., Powell I., and Bewsher, A., "A New Sensor for
5-1
T. Haubold
H. Gans
BMW Rolls-Royce GmbH
Material Technology
Hohemarkstr. 60-70
D-61440 Oberursel, Germany
D. Schwingel,
R. Taylor
University of Manchester
Materials Science Institute(UMIST)
Grosvenor Street
Manchester, M1 7HS, UK
SUMMARY
1. INTRODUCTION
Basic features for future advanced gas turbines
will be low emission rates and increased A basic feature for the development of advanced
efficiency. To fulfil these requirements a reduction gas turbines will be low emission rates of Nitrogen
of cooling air in the combustor chamber Oxides (NOx), Carbon monoxides (CO) and
demanding more efficient cooling technology Hydrocarbons (HC) to suit expected lower limits
and/or an increasing inner wall temperature is for pollution in the future. One way to fulfil these
needed. The development of thicker thermal requirements is a reduction of cooling air in the
barrier coatings for combustor tile applications is combustor chamber demanding more efficient
one concept followed to enable a reduction in cooling technology and/or increasing inner wall
cooling air. temperatures.
Thermal barrier coatings with a thickness up to Beside the development of more efficient cooling
2mm and improved thermal cycling life have been technologies there are several materials
developed as well as mechanical and technology concepts that are followed at BMW
thermophysical data determined. Rolls-Royce to reduce cooling air requirements in
the combustor:
Combustor segment rig testing demonstrated a i) a development of improved superalloys for
potential of 25% cooling air reduction by the use of combustor application with higher temperature
thick thermal barriers for the tile design capability,
investigated. ii) investigation of the potential of ceramic matrix
composite materials for combustor tiles and
iii) development of thicker thermal barrier coatings
(TBCs) for combustor tile applications.
Typical current applications of TBCs are thin has been performed at the BMW Rolls-Royce
coatings on guide vanes, blades and combustors APS spray facility (Sulzer Metco A3000S, 9MB
of gas turbines often introduced during overhaul in gun) using different spray parameters. During
order to extend lifetime of parts. Current spraying temperature has been measured by
technology is the use of plasma sprayed bond means of a pyrometer at the coating surface and
coat (BC) such as CoNiCrAlY mainly for adhesion by means of thermocouples at the backside of the
and oxidation protection and partially stabilised substrate. Substrates have been 2mm thick flat
ZrC>2 (about 8wt%Y2C>3) with a thickness below plates of Hastalloy X.
0.5mm. The benefit is a component temperature
reduction up to 150°C resulting in reduced cooling Microstructural investigations are performed at
airflow requirements and/or component durability BRR to determine typical coating features such as
improvement (1, 2, 3). Thickness and lifetime of thickness, porosity and crack density.
the barriers are limited by mechanical and thermal
stresses in the coatings due to the mismatch of Thermophysical and mechanical Properties
thermal expansion coefficients of the ceramic and
the metals and/or oxidation of the bond coat. Four point bend testing has been performed at
UMIST at room temperature and elevated
In a current Brite Euram project (BE7287)1 plasma temperature to determine the nominal strength,
sprayed thermal barrier coatings up to 2mm in the failure strain and Young's modulus. For testing
free standing TBCs have been used by dissolving
thickness have been developed to meet the
away the substrate. This has been possible due to
requirements of improved thermal insulation (4).
the thickness of the coating.
Different microstructures of coatings has been
addressed from different partners. Volvo Aero For the measurement at elevated temperatures, a
Corporation has investigated segmented TBCs high temperature rig has been designed in
(4), BMW Rolls-Royce 'conventional' highly porous accordance with European standards. All parts of
TBCs and ANSALDO Richerce has developed a the rig within the hot zone of the furnace are made
process route to spray highly porous TBCs using out of Alsint 99,7 which can theoretically withstand
a mixed ceramic-polymer powder (5). a temperature of 1700 BC. The capacity of the
furnace enables temperatures up to 1500 BC.
This paper will report on the development of the Figure 1 shows the design of the apparatus. The
highly porous coatings, the investigation of sample is supported by two rollers 20 mm apart.
thermophysical and mechanical properties and In order to minimise the torsional loading, one of
thermal cycling experiments as well as combustor these rollers can also rotate about an axis parallel
sector rig tests performed to determine the to the length of the sample. The specimen is
potential for cooling air reduction of thicker TBCs loaded by a second set of rollers which are
applied on combustor tiles. The sector rig test has isocentric to the first ones. They are 10 mm apart
been performed in the course of a national and both of them allow a rotation about their own
research programme. axis as well as about the axis of the sample.
' "V"";
the conditions are raised to take off rating with an
inlet air pressure of 20 bar and inlet temperature '. ■■■>!
of 885K. The facility limits the inlet pressure to 20 '"• ■.■:.-'•-•.■:£ ■„" ,■■■ u **'•', V- ;
. ,#."■*.-> V:-
bar. ■ ■,''•■ , " * •tf:rc-*."i '•■,■," ' ' ■ *';'".'■'.' i» '.'.% ■''•■ ~'T
•., ■ \ ' "Jx1*f .»■*''•- • .'' ■ ■•' •"' ■.-.'■'.-■ ; ,
w;> . •*■•»',■'.■ m-J ,-.v, >'■* ..*L •*.■•-..•*•»'■
Porosity values obtained for the different specimen Fig. 4: micrograph of TOP 7,15%porosity, x 150
are in the range of 8.5 to 20% (see fig. 4). For a
number of specimen segmentation cracking has
been found , which has been related to porosity Thermophvsical and mechanical properties
values in the range of 8.5 -12.5%.
Since the structure of the coatings is expected to
vary through thickness two different loading types
have been considered; where the top is either in
the state of compression or of tension.
TOP 8 surface temperature
k
Top 12 I8.8 j 30.4 iO.37
150
Top 13 I 7.8 125.1 I 0.35
100
Table 1: Mechanical properties of thick TBCs
Fig. 5: Failure stress versus Youngs modulus with p and cp being the density and the specific
heat capacity respectively (6). Values ranged from
The expansion coefficient of the samples tested 0.8 W/mK to 1.5 W/mK for the TBCs investigated
seems to have a rather constant value of 10.7*10" here.
6
K'1, which is typically for plasma sprayed 7%Y203
stabilised Zr02 (8, 11). Obviously, the porosity has
rather little influence. Over the whole temperature Thermal cycling
range, the thermal expansion coefficient is almost
constant and differences during heating and Thermalshock test at BRR (perrformed at
cooling, a typical indicator of phase changes, University Karlsruhe, IKM) used a cycling time of
cannot be determined. 60 sec for heating up (TTBC = 1250°C, Tsubstrate =
800°C) and cooling down to room temperature
a) and has been performed for several specimen up
0.006 to 2000 cycles yet in comparison to a cycling life
below 10 cycles for the first tests. Some specimen
tested to 2000 cycles did not fail at all. Fig 7 was
an overview about number of cycles obtained.
....T.!?.?.r.!53)..!?y!*!?.5.
ACKNOWLEDGEMENT
REFERENCES
Hot Gas systems by alloying with ceria and yttria [9]. This paper
Temperature discusses the theoretical models of heat transfer in solids
with regard to the reduction of the thermal conductivity,
Temperature Profile
without TBC and reviews these against the measured values on TBC
coatings engineered for low conductivity.
k.r =—an
3
T L [2]
'—'S.
\
.
1 0.1
«c
S
H
0.01
0.1 1 10 100
5 wt% Y203 20 wt% Y203
Wavelength (micrometers)
Figure 3 Emittance and transmittance of zirconia as a function Figure 4 Theoretical conductivity of EB-PVD zirconia ceramics
of wavelength. as a function of yttria content.
Whilst the thermal conductivity of air plasma sprayed Figure 7 Differences in the column density for the nucleate
zone (top) and the steady state region (bottom)
6-5
and distribution of columnar boundaries is responsible process used to produce the thermal barrier coating
for the reduction in thermal conductivity. More strongly affects the coatings resultant thermal
boundaries mean additional centres capable of scattering conductivity. Lamellar interfaces, parallel to the
thermal waves. Thus the thermal conductivity of an EB- coating/substrate interface as produced in plasma
PVD TBC can be modelled successfully using a two- sprayed ceramic, act to strongly scatter thermal waves
layer model, shown schematically in Figure 8 [19]. The reducing the coatings thermal conductivity. This
inner "nucleate" zone is dominated by defect/grain behaviour is observed to a lesser extent in the inner
boundary scattering and this results in a lower growth zone of EB-PVD deposited coatings due to the
conductivity of around 1.0 W/m.K. With increasing inclined nature of many of these inner boundaries.
thickness, the structure is characterised by a dominant Defect scattering must also play a role in this early
crystallographic texture and the thermal conductivity deposition zone of EB-PVD ceramics. However, due to
increases as the outer part of the coating becomes more the competitive nature of coating growth for vapour
perfect. In this outer zone the thermal conductivity deposited coatings a dominant growth texture is
approaches that of bulk zirconia at circa 2.2 W/m.K. developed. This vertically aligned microstructure of the
The overall coating thermal conductivity results from the outer layer of an EB-PVD coating results in thermal
combined effect of these two layers and can be predicted conductivities close to that for the bulk ceramic.
using a rule of mixtures. For coating thicknesses (x) less
than Xj (lOOum for the schematic illustrated in Figure 8) These observations would suggest that to further reduce
then the thermal conductivity of EB-PVD TBC's one should
introduce layers and/or additional phonon scattering
k = k, [5a] centres into the coating.
while for coating thickness (x) greater than x, 4. REDUCTION OF THERMAL CONDUCTIVITY IN
EB PVD COATINGS BY LAYERING
x
k_ \k\ (x~xi) j.
+
[5b] The former section discussing heat transfer into and
through the ceramic has highlighted a number of
Figure 9 provides a comparison between predicted and methods capable of reducing the thermal conductivity of
measured thermal conductivities for EB-PVD zirconia the zirconia ceramic. To reduce radiative transport
thermal barrier coatings, illustrating the good agreement 'colouring' and 'layering' are possibilities, while
to experiment of this two layer thermal conductivity introducing 'atom dimension defects', 'local changes in
model. ceramic density' and 'local variation in ceramic elastic
modulus' may be effective in reducing phonon
It should be expected that the thermal conductivity of conduction.
chemical vapour deposited (CVD) ceramic should show
a similar thickness dependence to EB-PVD zirconia as As indicated earlier the plasma sprayed TBC structure is
both coating morphologies are grown from the vapour such that the contribution to total thermal conductivity
phase. A study on the thermal conductivity of CVD by radiation is minimised by the lamellar splat
deposited synthetic diamond confirms this hypothesis boundaries at a typical spacing of 1 um. If this
[19,20]. attenuation of the radiative transport mechanism could
be applied to EB-PVD it could make a significant
From the foregoing, it is evident that the manufacturing contribution to the reduction of the thermal conductivity,
6-6
especially at high temperatures. Thus layering of the resistance. The concept is the subject of a European
ceramic appears a most promising route to lower the patent [22].
ceramics thermal conductivity, especially if it can be
coupled with charging the local density or elastic This paper focuses on the use of ion bombardment to
modulus of the material. A suitable periodicity of the create layers in the EB-PVD structure with varying
layers should be 0.2-2.0um (between X and Xf4 for an density.
incident radiation of 0.5-5.0 urn), which is able to scatter
photons and thus reduce radiative transport. This 4.1 Layering by Plasma Assisted Deposition
concept of layering within each column of an EB-PVD
coating is illustrated schematically in Figure 10. This method uses a glow discharge plasma to vary the
density of the ceramic during deposition. Figure 11
To reduce phonon conduction layer periods would have illustrates a typical micrograph of the morphological
to be of near atomic dimensions (0.3-2.0nm) and this is changes that can be introduced.
not considered feasible, although it should be possible to
locally change either the elastic modulus or density of The layers were produced by switching the D.C. bias
the ceramic by introducing atomic level defects into the applied to the substrate between high and low levels
lattice. during deposition. This has the effect of periodically
changing the degree of ion bombardment and thus
The PVD process offers a number of ways to tailor altering the density of the layers produced. The
coating microstructures and thereby introduce multi- micrograph illustrated in Figure 11, was produced by
layers, since both the column angle and morphology switching between a D.C. bias with a peak voltage of
may be influenced by variations in the flux density and 2500V and ground. The degree of bombardment that is
the direction this flux subtends to the substrate surface attained depends on the applied bias and local current
[21]. • density, current densities typically greater than
O.lmA/cm2 are required to achieve a significant density
Three alternative EB-PVD processing routes could be change [22].
used to generate layering within the columns.
The measured thermal conductivity for this
• Route 1 varying the speed of rotation of the parts microstructure is significantly lower than that for a
during coating coating produced without ion bombardments (Figure
12). Reductions of the order of 37-45% compared to
• Route 2 rocking the substrate between two fixed tilt state-of-the-art EB-PVD TBC's have been measured for
angles. these layered structures that are approaching the values
for thermally sprayed TBC systems.
• Route 3 varying the degree of ion bombardment to
change the density of the ceramic.
1.2
1.00
I3 1
T3
g 0.8
o 0.67
| 0.6 ■ 0.55 ^B 0 53
o
jE 0.4
a
I 0.2
0
EB-PVD Coloured PAPVD PAPVD Plasma
TBC EB-PVD Layered Layered Sprayed
TBC TBCver.1 TBCver.2 TBC
Clearly, the combination of layering at micron followed by a region with bias and high ion current
dimensions and the introduction of density change from density to produce maximum density changes between
layer to layer work in combination to significantly layers. This was followed by a region of low current
reduce the thermal conduction of the coating. As density giving low density differences between layers
discussed earlier the layering periodicity has been and finally a region of high current density giving a
selected to significantly reduce photon transport, while denser structure. (Figure 14).
local changes in layer density act to scatter phonons and
thus reduce thermal conduction by lattice vibrations. 4.2 Addition of dopants to reduce thermal conductivity
Using this process route it has been possible to introduce Tamarin and co-workers [17] have examined the
layers into individual columns, without change the addition of divalent transition metal oxides to reduce the
overall vertical columnar microstructure required for thermal conductivity of zirconia-yttria TBC's. As eluded
good strain, tolerance and good erosion resistance. For to earlier this has two effects, firstly to reduce the
example, it can be seen from Figure 13 that layering phonon transport in the material and secondly a
produces no significant deterioration to the erosion reduction in the radiative transport mechanism. The use
performance of an EB-PVD thermal barrier coating. of a divalent transition metal oxide introduces vacancies
as well as strain centres into the lattice both of which
Also it is possible by programming the sequencing of will reduce the phonon mean free path. In addition this
bias application to tailor the position of the layers in the doping changes the colour of the TBC material to dark
coating. This potentially has benefits in allowing the grey, reducing radiation transport in the visible range
properties of the coating to be varied from the surface to and by inference in the near infra-red.
the interface. For example more erosion resistant
structures could be generated at the outer surface of the Following this concept of dopant additions, further
coating with the application of bias in this region, whilst reduction in phonon thermal conductivity could be
the inner zone of the coating is layered to lower the achieved by the selection of low levels of dopant which
thermal conductivity. To demonstrate this a coating was maximise lattice strain and lattice anharmonicity. This is
produced where the early deposition has no bias applied, best achieved by the addition of transition metal ions
with high mass and large ionic radii.
7.00
and metallic bond coat of the TBC. The coat is deposited by VPS and the top coat
rough surface of the TBC, inherent to this by APS. NASA's Marshall Space Flight
approach and particularly of the weld Center [5] has reported a five fold increase
regions, leads to an undesirable change in in thermal cycling over the conventional
flow pattern of the hot gas. Moreover, TBC, when both the bond and top coats are
because the current fabricating process applied via VPS. Furthermore,
consists of mechanically forming sections PyroGenesis has experimented with
of the component, there is a limitation on advanced TBCs which can be much thicker
the choice of applicable superalloys (only than conventional zirconia and offer
superalloys with relatively high elongation significantly improved thermal protection.
can be used). Finally, the cost of drilling
the thousands of small holes required for PyroGenesis has committed itself to
cooling is exorbitant and the use of cooling addressing the problems associated to
air significantly adds to the size of the current combustion system components and
compressor. their fabrication method, by combining the
advantages of spray forming and advanced
There have been numerous studies on spray TBCs with those of the VPS technology.
forming of metals and ceramics via VPS
[1, 2, 3, 4]. However, to the authors'
knowledge, there has not been any PROCESS DESCRIPTION
published work on spray forming complex
shape components with a ceramic inner A suitable mold material was selected and
layer, followed by multiple metallic layers. surface machined to the inner geometry of
A closed component with an inner ceramic the desired component. The mold was then
layer has potential applications for thermal cleaned and placed inside the VPS
barrier or wear protection components. In chamber, ready for service. The following
the absence of film cooling, the use of a unit operations were then used to fabricate
conventional thermal barrier coating does multilayered combustion system
not provide adequate insulation to protect components by VPS spray forming:
the component. The temperature drop (i) mold preconditioning; (ii) preheating
across conventional TBCs is approximately the mold to a suitable temperature and
120°C. In order to protect the combustion depositing the TBC layer; (iii) preheating
liners and transitions ducts of an advanced the mold to a higher temperature and
turbine operating without film cooling, a depositing the bond coat and the structural
temperature drop across the barrier of superalloy; and (iv) removing the mold
approximately 350°C is required. Thus, from the near net-shape component.
both a new fabrication technique and a Figure 1 illustrates the process.
dramatically improved TBC is needed for
the proposed application. The design of the mold is of particular
importance in the VPS spray forming
PyroGenesis' experience in VPS-applied process. The mold is used both to define
zirconia based TBC systems (both top and the shape of the component and to control
bond coats) for gas turbine engine the energy balance of the process, a critical
applications has revealed superior thermal issue in controlling the stresses within the
cycling performance at a lower cost over spray formed structure.
the commonly used TBCs, where the bond
7-3
Figure 3 is a micrograph of IN-738LC after • • •- -, ■*, M'Wz • - it, it," " iti', ?*"•> J , • I ' :■ * ' *» ---" -"' '■* «F ■
50 microns
50 microns
300 microns
300 microns
Thermal cycling was performed on the The results from the thermal conductivity
PSZ/Ca2Si04 TBC by heating Hat coupons study are presented in Figure 8. Figure 8
to 1,100°C with an oxyacetylene torch and shows that the AT between the hot and cold
cooling rapidly with compressed air. faces of a uncoated stainless steel coupon is
Following 300 cycles no delamination was approximately 200°C when the hot phase
observed. reaches 1,000°C. By comparison, a
conventional PSZ TBC structure offers a
A number of PSZ/Ca2Si04 TBCs were AT of 275°C (Sample Z-10). Under the
produced on flat stainless steel substrates same operating condition an advanced 1.8
and tested for their thermal barrier mm thick PSZ/Ca2Si04 TBC offers a AT of
properties. The test was performed by approximately 420°C (Sample ZC-9).
heating the TBC face with two high Thus, the thermal barrier properties of the
intensity infrared lamps and measuring the PSZ/Ca2Si04 developed in this study can be
temperature of both the hot TBC surface almost 3x better than a conventional TBC,
and the cooler stainless steel surface. A which should be sufficient to meet the
list of the samples evaluated in this study is requirements of the proposed application.
shown in Table 3.
300 ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' '
Sample Z-10 represents a conventional PSZ 650 700 750 800 850 900 950 1,000
TBC sprayed by VPS. Sample ZC-9
Temperature onlhe hot surface, (°C)
represents a TBC which theoretically
should provide a AT in excess of 350 °C
and thus meet the performance
requirements for advanced combustion Figure 8. Thermal properties of
liners and transition ducts without film PSZ/Ca2Si04 advanced TBCs.
cooling.
7-8
SUMMARY
REFERENCES
Y.Jaslier, S.AIperine
SNECMA - Materials and Processes Department
Site de Villaroche - Bätiment 41
77550 Moissy-Cramayel
France
The early application of TBC's was the protection of In the Eastern countries, the EB-PVD technology was bred in
combustor parts. This application involved applying the prestigious institutes such as the Paton Welding Institute in
ceramic coat on top of an MCrAlY coating used as a bond Ukraine and the VIAM in Russia, with industrial recipients
coat, both deposited by air plasma spraying. One advantage of such as Saturn and Nickolai Kuznetzov (NK). The first
the air plasma spaying technique is to provide the MCrAlY Structure Zone Model for physical vapour deposited coatings
coating with a rough surface onto which the ceramic coating was first published in 1963 by Movchan and Demchishin5
can mechanically bond. Another advantage of air plasma and has been universally referred to since. The design
spraying lies in the heavily microcracked ceramic approach of Russian jet engine manufacturers was such that it
microstructure it produces which is beneficial in terms of was not felt there was a need for coating HP vane aerofoils
strain tolerance and thermal shock resistance. with EB-PVD TBC's. As a result Russian experience in
production was bound to coating blades, not vanes6. This
Since the potential benefits of thermal barrier coatings were detail bears some implications on the development of tool
additive to the gains in high temperature strength capability design and masking technology (or lack of) in the former
of superalloys and in internal and film cooling technologies,
level of a few 10"4 Torr. There is no intentional bleed of a although poorly understood, is known to be moisture related.
reactive gas into the chamber. A typical coating cycle consists Other investigators have reported the effect of moisture on
first in mounting the carrousel holding the parts on the the spallation resistance of thermally grown oxides on bare
driving shaft inside the loading chamber. After the loading superalloys15. In these tests, thermal cycling was interrupted
chamber has been evacuated, the locking vane is opened and every 20 cycles with the test pieces left at room temperature
the carrousel with the cold parts is introduced into the for a time up to 4 hours, after what the specimens were
deposition chamber by translation of the shaft. The parts are inspected and the cycling resumed. A dwell time of 4 hours at
pre-heated inside the deposition chamber by the heating room temperature was chosen so as to let enough chance for
electron beam gun while the evaporating guns are kept shut the moisture to have an effect, if any. The failure criterion
down. Evaporation is started after the work pieces have was a ceramic spalled area of at least 20% of the specimen
reached the right temperature regime. Once the deposition is coated surface.
completed the carrousel is driven back into the loadlock
where it is left to cool down before it is unloaded. RESULTS
It should be noted at this point that Western coaters are Characterization of as-deposited TBC's
generally designed to be dedicated to one type of coatings i.e. The overall TBC systems are presented in Figure 3. It may be
metallics or ceramics whereas, in the Eastern art, the coaters seen that each ceramic coating system exhibits a specific
are designed and used to deposit both types of coatings. ceramic morphology. For a same nominal YSZ composition,
Whether the latter approach is better or worse from the the RT33 process produces a coarser, more columnar shape
production standpoint than the dedicated type coater approach morphology compared to the finer, feather-like, KDPi
is not the debate of the present discussion. The point is that ceramic morphology. Moreover, the KDPI feathered columns
the flexibility of the Eastern type coaters permits to evaporate coarsen towards the outer surface. This effect is not as
metallics concurrently with ceramics or to deposit both in dramatic in the RT33 morphology.
sequence in the same coating cycle and this is essential to
their art14.
wi
Ni Cr Fe Mo Co W
bal 22 18 9 1 0.2
Characterization
The morphology of the ceramic layers was evaluated using a
field emission gun scanning electron microscope (FEG-
SEM). Energy Dispersive X-Ray analysis coupled with SEM
Acc.V Spot Magn Det WD
was used to perform chemical analysis in both the as- 115.0 kV 4.0 360x BSE 10.0 VTP21(SDP?-Ätathni«V7fl1
deposited and as-tested state. Ceramic phase composition was Figure 3: SEM back-scattered electron image of the as-
determined using CuKa X-Ray diffraction. deposited RT31ART33 (a) and SDP2/KDP1 (b) TBC
systems
TBC performance was evaluated using a thermal cycling test.
Thermal cycling was conducted on 6 and 10mm diameter bars High magnification examination of the ceramic layers at 50
in a furnace with a cycle consisting of 55 minutes at 1100°C micron height from the bond coat gives some insight into the
with a 5 minute heat up and 15 minute cool down to room column sub-structure (Figure 4). Back-scattered electron
temperature using forced air convection cooling. The imaging reveals internal striation consisting of alternate dark
probability for the ceramic to spall off after many cycles and white layers. This effect is a lot less pronounced in the
increases with dwell time at room temperature. This effect,
8-4
case of the KDP1 ceramic where only minor striation may be Overall bond coat composition was determined using semi-
distinguished in the core centre line of the feathered columns. quantitative EDS analysis over an 80x80microns window
centred within the bond coat depth. The analysis results are
shown in the table below:
Co Ni Cr Al Y Fe
RT31 wt% bal 38 18 8.9 <1 <1
at% 34 18 17
SDP2 wt% - bal 14 12 <1 <1
at% 14 22
1BBBB
[counts]
9BBB
-
KDP1
,
(28B)
5806
4aee
3BBB
(488)
RT33 (311)
. \
e
£V 38
JB_ 4B 5B
SB 68
bö 78
YB 1i'iei
4B I B'B
BB
ilPfPiiW \\l ■o
"55
o
a
W)
rh
/
41
0 - In
SDP2/KDP1
, Ilk
RT31/RT33
Acknowledgements: The work concerning the study of nickel-base superalloy », Materials science and Engineering
Russian TBC's was supported by the SPAE.The authors A203 (1995) 388-398
gratefully acknowledge the contributions of E.Celerier and
M.Scher for the thermal cycling work as well as 16.D.V.Rigney, R.Viguie, D.J.Wortman, D.W.Skelly, «PVD
J.M.Duchemin and J.M.Thubert for the SEM work at thermal barrier coating applications and process development
Snecma. The X-ray diffraction studies were performed by for aircraft engines», TBC Workshop 1995, NASA
C.Diot at Onera. Fruitful contacts and discussions with Pr Conference publication 3312
Y.Tamarin at the VIAM Institute in Moscow were much
appreciated. Last but not least, the authors would like to 17.Yann Jaslier, «Development of EB-PVD TBC's: the role
thank Remy Mevrel at Onera for returning his useful of deposition temperature and plasma assistance», PhD
comments on the first draft of this paper. thesis, Cranfield University, 1995
2. S.M.Meier, D.K.Gupta, K.D.Sheffler, Ceramic TBC's for 20. Maricocchi, Rigney, US Patent 5,350,599
commercial gas turbine engines, JOM, March 1992 Sep. 1994
3. K.D.Sheffler, D.K.Gupta, Current status and future trends 21.LLelait, S.Alperine, C.Diot «Microstructural
in turbine applications of TBC's, Gas Turbine Aeroengine investigations of EB-PVD TBC's », Journal de Physique IV,
Congress and Rxp., Amsterdam, Sep. 30 1988 ASME 88-GT- Colloque C9, supplement au journal de physique HI, vol 3,
286 dec 1993
4. Comments from W.Goward, D.Boone, W.Hainan, 22.M.G.Scott, «Phase relationships in the zirconia-yttria
TBC Workshop 1997, May 19-21, Cincinatti, Oh system », J.Mat.Sci, 10 (1975), pp.1527-1535
W. A. Kaysser
M. Peters
K. Fritscher
U. Schulz
DLR
German Aerospace Research Establishment
Institute of Materials Research
D-51170 Köln
Germany
techniques or by advanced high temperature materi- lives by avoiding hot spot formation and subsequent
als [4] and in particular through the introduction of failure by thermal fatigue. These materials were
electron-beam physical vapour deposition (EB-PVD) broadly accepted for 20 years until their replacement
thermal barrier coatings (TBCs) [5]. by PYSZ (partially yttria-stabilised zirconia) coatings
was initiated. The APS (air-plasma sprayed) PYSZ
TBCs offer outstanding mechanical, chemical and
thermal properties. Two generations of high tem-
perature (>1000°C) TBCs were based on these
unique ceramics (generation I is APS MCrAlY / APS
PYSZ and generation II is LPPS (low-pressure
plasma sprayed) MCrAlY / APS PYSZ) which allowed
in the early eighties to introduce TBCs on highly
thermally loaded parts like vane platforms and vane
airfoils. The need for higher operating temperatures
in turbines, however, inspired material scientists for
half a decade to search for TBCs which could also be
applied on blade platforms and airfoils which are
200
1950 1960 1970 1980 1990 2000 2010 2020 subjected to high thermal and additional mechanical
loads.
Service Date of Gas Turbine Aero Engines
Figure 1: Increase of turbine inlet temperature of Here the electron-beam physical vapour deposition
aircraft gas turbine engines [2] (EB-PVD) technology offered the opportunity to gen-
erate TBCs with vastly superior strain tolerance. This
is due to their specific coating structure which grows
2. THERMAL BARRIER COATINGS from the vapour phase in a columnar form with indi-
vidual ceramic columns being weakly bonded to their
TBCs consist of thin ceramic layers of low thermal neighbour columns as e.g. shown in Figures 8a or
conductivity - typically partially stabilised zirconia - 13.
which are applied on aerofoil surfaces that just have
a metallic corrosion resistant coating. The coating Coating all these high performance parts with reliable
imparts good adhesion of the ceramic to the sub- TBCs would allow to alter the initial component de-
strate. Application of the TBCs enables increasing sign and to operate the TBC bearing components
engine performance/thrust by either increasing the with the coatings as an integral part of them. The
gas temperature or reducing the cooling airflow. exploitation of strain-tolerant TBCs initiated a renais-
Alternatively the lifetime of the turbine blades can be sance of industrial EB-PVD technology in 1987 which
extended by decreasing metal temperatures as sche- until then had been utilised for deposition of metallic
matically outlined in Figure 2 [6, 7]. MCrAIY-type coatings for 20 years. The novel EB-
PVD PYSZ ceramic layers with advanced LPPS
MCrAlY bond coat, designated now as generation III
TBC, offered a roughly 3 times improvement in blade
life or a surface temperature increase on the hot
sections of airfoils of approximately 150K.
4. PROCESSING OF TBCS
Since the early sixties PS (plasma sprayed) calcia- During EB-PVD processing a high energy electron
and magnesia-stabilised zirconia TBCs have been beam melts and evaporates a ceramic source ingot in
used extensively as ceramic top coat on combustion a vacuum chamber. Ingots are bottom fed into the
chamber walls and on burner cans to prolong their crucibles during evaporation to ensure continuous
9-3
TBC growth. To achieve defined stoichiometry of the the coating systems. Failure in TBC systems almost
zirconia a controlled amount of oxygen is bled into always occurs by TBC spallation due to stresses in
the deposition chamber. Preheated substrates are the ceramic or in the bond coat. But failure of EB-
positioned in the vapour cloud above where the va- PVD processed TBCs happens, unlike those for PS
pour is deposited on substrates at deposition rates of based TBCs, at the ceramic - bond coat interface.
100 to 250 nm/s [9]. Typical columnar microstruc- There is a thermally grown oxide (TGO) of a few urn
tures and aerodynamically smooth surfaces are ob- in thickness which plays a key role in the adhesion of
tained without the need for final polishing or condi- EB-PVD TBCs since it is considered to be the weak
tioning of cooling holes. Due to the columnar micro- link in the system. Its strength and adhesion to the
structure the life time of the TBCs is prolonged and bond coat governs the true location of failure whether
the damage tolerance improved. Typical characteris- first cracks will propagate within or along the TGO -
tics and major advantages of thermally sprayed and bond coat interface. Here spallation is initiated by
on-evaporated TBCs are schematically outlined in accommodation of stresses surmounting a critical
Figure 3 [10]. Figure 4 shows a selection of EB-PVD value.
TBCs on aero engine and stationary gas turbine
blades produced at DLR on semi-commercial single-
source 60kW Leybold and dual-source 150kW von
Ardenne EB-PVD coaters.
PS EB-PVD
TBC»-
BC »-
[HIM11 Figure 4: EB-PVD TBCs on turbine blades for use in
stationary gas turbines, civil aero engines and heli-
copter engines (decreasing in size)
Substrate
merits the modified high creep resistance phase 7. INFLUENCE OF STABILISERS ON TBC
structure offers further benefits which refer to a low PERFORMANCE
diffusivity and a reduced coefficient of thermal ex-
pansion (CTE). The latter property is of uppermost Partially yttria-stabilised zirconia (PYSZ) is the cur-
importance for the amount of residual stresses ac- rent state-of-the-art material for TBCs. Unfortunately,
cumulating in the TGO layer. They can be reduced if the material shows insufficient phase stability and
the differences between the CTEs of the respective accelerated sintering at temperatures above 1200°C.
ceramic and metallic partners of the whole system Therefore, further increases of gas inlet temperature
are minimised. require alternative stabilisers with improved phase
Finally the involvement of superalloy substrates in stability. Another argument for new stabilisers in
the interplay of the physical properties of the respec- zirconia or for completely new ceramics is the insuffi-
tive material components and their influence on the cient resistance of current TBCs against chemical
compositional behaviour of TBC systems has to be attack by pollutants in the combustion gas. Hot cor-
addressed. Ni-base superalloys have gained a high rosive decay of TBCs by Na2S04 and vanadates are
degree of technical maturity in the early nineties. reported which involves leaching out of stabilisers
They have passed a lot of development stages, e. g. from parent zirconia. The degradation occurs by
hardening by precipitated y phases, extensive solid transformation of high temperature phases to mono-
solution strengthening, directional solidification and clinic zirconia on cooling. This failure mode may be
finally single crystals to trust in more creep-resistant found mainly in heavy duty engines and off-shore
crystal orientations and continuous grain structures service [15].
[11]. The most creep-resistant crystal orientation
(001), however, provides the least Young's modulus One alternative stabiliser for zirconia is Ce02- The
of only 120 MPa. Obviously it is a drawback for the benefits of ceria-stabilized (CeSZ) TBCs are good
adhesion of TBCs as thin-walled airfoils can be bent corrosion resistance [16-17] and excellent phase
more rigorously under same loads in service than the stability at high temperature [18]. Furthermore, the
conventionally cast counterparts. Three generations thermal conductivity is found to be lower than for
of single crystal superalloy development, however, PYSZ and some benefits for lifetime and thermocyclic
resulted in a reduction of the CTE which compares resistance are reported, too. La203 is another candi-
with the recent low expansion bond coats.
date for replacing Y2O3 in zirconia.
above. Other effects that affect the microstructure Burner rig tests as well as furnace tests revealed that
include variations in phase composition, ion radii, and composition and phase structure of the variously
surface energy aspects during condensation [22]. stabilised zirconias are closely related to the cyclic
lifetime of the coating system. The only phase found
Phase analyses by X-ray diffraction (XRD) revealed in PYSZ was metastable non-transformable t'. The
that the compositions of the three binary TBCs, equilibrium phase diagram [23] predicts a two-phase
PYSZ, FYSZ, and LaSZ, were close to the ingot mixture consisting of t containing 4 % Y2O3 and c
compositions. For PYSZ exclusive tetragonal (f) containing 16 % Y2O3 at deposition temperature.
phase was identified while FYSZ only contained The c phase, however, can only be maintained at RT
cubic (c) phase and LaSZ exhibited a mixture of under the preconditions that a very low critical grain
mainly tetragonal, minor cubic, and substantial size of c is established and/or the c particles are
amounts of monoclinic (m) phase. In the case of the subjected to high compressional stresses by the
ternary composition Zr02-Ce02-Y2C>3, however, surrounding t matrix. Both mechanisms would cause
analysis showed no constant but fluctuating compo- the half widths of the XRD peaks to be broadened
sitions across the TBC thickness. Due to these varia- which, however, is not observed. Therefore, forma-
tions the surface content of ceria varied among test tion of the t' phase is considered to be caused via
pieces of different deposition runs between 13 and rapid quenching on PVD processing as can similarly
38 wt%. For CeSZ a mixture of cubic, tetragonal, and be assumed for solidifying TBCs during PS deposi-
occasionally monoclinic phases was found [22]. tion.
300- PYSZ
FYSZ
OJ250-
E, /V LaSZ
8 200- CeSZ^-
0
£ 150-
cn
Q)
g 100-
50-
01 k—O rO—-rO-—rO- 1— i
10 20 30 40 50 60 70
Time [h]
temperatures up to 1150°C, they transfer after 100 the fluctuation of composition across the thickness of
hrs annealing at 1400°C to a mixture of 48% tetrago- the TBC. This may cause alternating layers of "weak"
nal + 48% cubic + 4% monoclinic (mol%) phases phases including monoclinic and "strong" phases like
[24]. The transformation into three phases suggests a t' [24]. Apparently, the coating will first break in
sequential mechanism consisting of yttrium cation weaker regions. Lattice misfits between the various
diffusion out of t' and destabilisation of t' into c and t layers may further contribute to reduced adhesion
with subsequent phase transformation to m during between the layers. Higher erosion rates of CeSZ
cooling. At temperatures below the stability point t' is found for PS TBCs [28, 31] may also account for this
still the most effective phase for durability in thermal failure mode. The partial loss of this TBC may act as
barrier coatings as can be seen in Figure 5. The t' a strain/stress relief mechanism that allows the re-
microstructure is characterised by outstanding bend- mainder part of the coating, which is thinner due to
ing strength, high crack propagation energy, high partial spallation, to adhere longer. CeSZ apparently
fracture toughness and highly tolerant thermoshock offers some potential in thermal barrier coatings es-
behaviour. Lattice distortion due to tetragonality, a pecially if a columnar microstructure can be utilised.
domain structure and the well-described tweed micro- Multiple source evaporation is however required to
structure inside the t' grains [25] are supposed to be bring about reproducible CeSZ TBCs of stable com-
responsible for the excellent performance of PYSZ. position.
often vary in diameter from root to top or from one need of appropriate column features in particular
column to the other. Columns at the root section are TBC thickness zones. Several aspects may contrib-
much thinner than at the top and are enlarging dis- ute to the observed effects including surface tem-
continuously from root to top conically. Increasing perature fluctuations, change of diffusion kinetics and
both temperature and rotational speed improves the enlarged shadowing due to rotation. For instance,
regularity and parallelity of the microstructure and measurements of actual surface temperatures on flat
enlarges the column diameter. After competitive samples during deposition gave temperature differ-
selection of favoured crystal orientations in the thin ences of 60 to 80°C for a low rotational speed com-
root area during the first stages of coating growth, pared to 15 to 25°C for a three times higher speed.
columns start immediately to grow in width to their
terminal diameter. Although the coating density
measured is higher for high temperature/high rota-
tional speed TBCs than in low temperature deposited
TBCs, the columns appear to be less densely
packed.
30 min-'
12 min'
\AAA/
The good performance of TBCs on highly loaded In hot corrosion testing of EB-PVD TBCs no evidence
parts in severest environments, achieved in less than was found for hot corrosion of the bond coat and for a
a decade of development, is a good indication for the chemical reaction between the TBC and corrosion
versatility of the columnar structure concept of TBCs. accelerators like molten salts and sulphur com-
The unfavourable crystallographic relations between pounds [37, 38]. However, TBC failure can be initi-
the columns and the direction of highest straining in ated by mechanical attack of solidified salt com-
service, hypothetically addressed before, stand for a pounds. Spallation of small areas of TBC in a funnel-
positive view into the future of TBCs: they still bear a shaped manner with some thinner coating still adher-
potential for further improvements. ent on the bondcoat indicate this failure mechanism.
For columns that posses larger voids between the
interior column faces, salts may penetrate deeper
11. STRUCTURE PROPERTY RELATIONSHIPS into the coating than in denser columnar microstruc-
tures. Consequently, the weight loss is higher for a
Some representative microstructures of TBCs pro- coarser microstructure deposited at higher rotational
duced by single source evaporation were chosen for speed compared to a microstructure deposited at
cyclic furnace testing (Figure 12). All samples had lower temperature and rotational speed [37].
the same EB-PVD NiCoCrAlY bond coat and surface
600
treatments. Samples with a microstructure consisting
of tapered and discontinuous columns due to multiple ^- SDallation
500-
branching on growth exhibit the longest cyclic lifetime •■**&■; i»
(12rpm/1020°C). TBCs with more regular columns of
400- -«.-.
larger diameters appear less tolerant against thermo-
cyclic loading. Burner rig tests confirm these results. (/)
| 300
Obviously, microstructure of TBCs and thermocyclic >.
lifetime are closely interlinked. Figure 12 evidences Ü
200 - ws*$$
the superiority of TBCs with a non-regular micro-
structure consisting of columns with variable diame-
ters over TBC with large uniform column diameters
100 SÜ
and a regular microstructure.
0
12/1020 30/920 30/1060
Revolutions [RPM] /Temperature [°C]
Figure 12: Cyclic lifetime in furnace test between
1100°Cand130°C
In Figure 10 two EB-PVD PYSZ coatings are com- This question in mind the habit morphology and ori-
pared with vapour incline angles of 25-40° (Figure entation of a TBC will be detailed first. EB-PVD TBCs
10a) and < 10° (Figure 10b) deposited under con- form a preferred {100} texture not only in normal but
stant rotating conditions. It is obvious that the colum- also in in-plain direction which will establish along the
nar structure closely follows the incline angle of the main axis of rotation and perpendicular to it (Figure
vapour cloud to the substrate surface. This shows 11). This axis usually coincides with the main axis of
that the angle of the columns is directly related to the a turbine blade which also is most heavily subjected
vapour incline angle, which is, however, different for to mechanical stressing in service. In essence, the in-
the crystallographic orientation. Texture measure- plain {100} texture orientation of TBCs is the most
ments have revealed that the deposition conditions heavily strained one. Any other crystal orientation
clearly influence the type of texture of EB-PVD TBCs should be better. A look over a TBC in lateral direc-
[36]. For rotated substrates a four-fold symmetry with tion (Figure 7b) clearly shows the longitudinal and
{100} planes nearly parallel to the surface evolves often diamond-shaped cross section of crystals which
independent of vapour incline angle, as indicated in spread with their elongated main diagonals in {100}
Figure 10. The technical importance of textures in direction. So less stress-depressant boundaries will
TBCs is related from the fact that the associated accumulate in this direction compared to all others.
anisotropy in elastic moduli as well as thermal ex-
pansion within the t' phase may very well have a
profound influence on the lifetime of the coating sys- (113) (100)
tems.
14. REFERENCES
[16] B.A. Nagaraj, A.F. Maricocchi, D.J. Wortman, Coatings", Surface and Coatings Technology,
J.S. Patton, and R.L.Clarke, "Hot Corrosion Re- 50(1992)141-149
sistance of Thermal Barrier Coatings", ASME [32] R.E. Demaray, W.K. Hainan, and S. Shen,
paper 92-GT-44 (1992) "Development of Electron Beam Physical Vapor
[17] B.A. Nagaraj and D.J. Wortman, "Burner Rig Deposition of Ceramic Coatings", Proc. 2nd
Evaluation of Ceramic Coatings With Vanadium- Conf. on Advanced Materials for Alternative
Contaminated Fuels", ASMEJ. Eng. Gas Tur- Fuel-Capable Heat Engines, (Palo Alto, CA:
bine Power112 (1990) 536-542 Electric Power Research Institute 1981), pp.
[18] P. Vincenzini, G. Appiano, F. Brossa, and S. 151-184
Meriani, "Stability of Thermal Barrier Coatings", [33] D.V. Rigney, R. Viguie, D.J Wortman, and D.W.
Proc. 3rd Int. Symp. "Ceramic Materials and Skelly, "PVD Thermal Barrier Coating Applica-
Components for Engines", ed. V.J. Tennery tions and Process Development for Aircraft En-
(1989), pp. 201-210 gines", NASA Conference 3312 (1995) pp. 135-
[19] U. Schulz, K. Fritscher, M. Peters, "Thermocyclic 149
Behavior of Variously Stabilized EB-PVD Thermal [34] K. Fritscher and U. Schulz, "Burner-Rig Perfor-
Barrier Coatings", Journal of Engineering for Gas mance of Density-Graded EB-PVD Processed
Turbines and Power, 119 (1997) in print Thermal Barrier Coatings" in Ceramic Coatings,
[20] J.A. Thornton, "Influence of Substrate Tempe- ed. K.Kokini, ASME MD-Vol.44 (1993), New
rature and Deposition Rate on Structure of York, pp. 1 -8
Thick Sputtered Cu Coatings", J. Vac. Sei. [35] U. Schulz, K. Fritscher, C. Leyens, M. Peters,
Technol., 12 (1975) pp. 830-835 W.A. Kaysser, Thermocyclic Behavior of Differ-
[21] B.A. Movchan and A.V. Demchishin, "Study of ently Stabilized and Structured EB-PVD Ther-
the Structure and Properties of Thick Vacuum mal Barrier Coatings, Materialwissenschaft und
Condensates of Nickel, Titanium, Tungsten, Werkstofftechnik 28 (1997) 370-376
Aluminium Oxide and Zirconium Dioxide", Fiz. [36] U. Schulz, H. Oettel, W. Bunk, "Texture of EB-
Metal. Metalloved., 28 (1969), pp. 83-90 PVD thermal barrier coatings under variable
[22] U. Schulz, K. Fritscher, and M. Peters, "EB-PVD deposition conditions", Zeitschrift für Metallkun-
Y2O3 and Ce02/Y2C>3 Stabilized Zirconia de 87 (1996) 488-492
Thermal Barrier Coatings - Crystal Habit and [37] U. Schulz, K. Fritscher, H.-J. Rätzer-Scheibe,
Phase Composition", Surface and Coatings M. Peters, and W.A. Kaysser, "Thermocyclic
Technology, 82 (1996) 259-269 Behaviour of Microstructurally Modified EB-
[23] H.G. Scott, "Phase Relationships in the Zirco- PVD Thermal Barrier Coatings", Materials Sci-
nia- Yttria System. J. Mat. Sei., 10 (1975) 1527- ence Forum, in press
1535 [38] B.A. Nagaraj, A.F. Maricocchi, D.J. Wortman,
[24] U. Schulz, K. Fritscher, and M. Peters, J.S. Patton, R.L Clarke, "Hot corrosion resis-
"Thermocyclic Behavior of Variously Stabilized tance of thermal barrier coatings", ASME paper
EB-PVD Thermal Barrier Coatings", ASME 92-GT-44(1992)
Turbo Expo, 10-13.6.1996 Birmingham, ASME- [39] P. Lenk, J. Senf, B.-D. Wenzel:" Von Ardenne
paper 96-GT-488 guns and their suitability for PVD processes for
[25] L. Lelait, S. Alperine, C. Diot, and R. Mevrel, thermal barrier coating", Proc. conf. on electron
"Thermal Barrier Coatings: Microstructural Inve- beam melting and refining-state of the art, 1996;
stigation after Annealing", Mater. Sei. Eng., ed. R. Bakish, Bakish Materials Corporation,
A121 (1989)475-482 Engelwood, USA, 182-190
[26] R.A. Miller and C.C. Berndt, "Performance of
Thermal Barrier Coatings in High Heat Flux En-
vironments", Thin Solid Films, 119 (1984) 195-
202
[27] S. Stecura, "Effects of Compositional Changes
on the Performance of a Thermal Barrier Coat-
ing System", NASA Technical Memorandum
78976(1979)
[28] N.P. Anderson and K.D. Sheffler, "Development
of Strain Tolerant Thermal Barrier Coating Sys-
tems", NASA Contract NAS3-22548, Report No.
NASA-CR-168251 (1983)
[29] O. Unal, T.E. Mitchell, and A.H. Heuer,
"Microstructures of Y203-Stabilized ZrÜ2 Elec-
tron Beam-Physical Vapor Deposition Coatings
on Ni-Base Superalloys", J. Am. Ceram. Soc,
77(1994)984-992
[30] N.S. Jacobson, "Thermodynamic Properties of
Some Metal Oxide-Zirconia Systems", NASA
TM 102351 (1989), pp. 1-63
[31] R. Taylor, J.R. Brandon, and P. Morrell,
"Microstructure, Composition and Property Re-
lationships of Plasma-Sprayed Thermal Barrier
11-1
S. Chevillard
S. Drawin
M.H. Vidal-Setif
2.1 Zr and Y source materials ZrCl4 powder is commercially available at high purity
levels (more than 99.5%) and presents high vapor
Since neither mineral nor metal-organic zirconium and pressures (315 to 4000 Pa) in the range 200-250°C. This
yttrium compounds are gaseous at room temperature, allows the design of a vaporization furnace working at
condensed source materials have to be used and relatively low temperature. In practice, a porous silica
vaporized in a specific manner (bubbler, furnace, etc). crucible is half filled with ZrCl4 and placed in the
The precursors for Zr and Y were chosen from literature temperature-controlled furnace (70 cm3 inner volume).
data and specific investigations were conducted on some Argon carrier gas is then forced through the powder and
11-3
the ZrCl4-enriched gas flows to the mixing chamber. This deposition experiment and is continuously kept under
percolation technique ensures a good solid to gas mass vacuum. No excessive corrosion is noticed after three
transfer. years operation.
The metal-organic compound Y(thd)3 is available The whole gas generation part, from the furnaces
commercially at slightly lower purity levels (more than downwards to the injection nozzle, is heated at 250°C to
98%). It melts at 174°C and presents a high enough avoid condensation of the gaseous precursors.
vapor pressure at temperatures greater than ca. 160°C. A
2.2 Gas injection
silica boat with a small amount of Y(thd)3 is placed in a
second temperature-controlled furnace. Argon circulates The injection of the precursor gases into the deposition
over the powder surface, which is pasty at the selected chamber is of primary importance as coatings with
temperatures. The Y(thd)3 + Ar gas mixture is then homogeneous thickness are required. This is especially
introduced into the mixing chamber. the case if real complex-shaped turbine blades are to be
coated. But the present laboratory-scale deposition
Oxygen and argon are added to the gas flow in the
reactor has been designed for investigations on small
mixing chamber. Argon is used as an inert carrier gas but
samples, so that this requirement is less severe.
also to facilitate the ignition of the plasma, to improve its
Nevertheless, both to get insight in the fluid dynamics
stability and to adjust the partial and total pressures in the
around the sample and to prepare an up-scaling of the
deposition chamber.
reactor, flow simulations in the deposition chamber have
Only a few grams of Zr and Y precursors are vaporized at been performed [24].
each deposition experiment, as given by crucible mass
The small cylindrical nozzle used by Bertrand [22] has
measurements. This corresponds, in the standard
been replaced by a conical convergent-divergent nozzle
conditions (see Table 1), to approximate mass loss rates
as shown on Figure 2. The nozzle has been designed
of 1.4 g/h and 0.25 g/h and gaseous precursor mole
using one-dimensional isentropic flow relations [26] so
fractions in the mixing chamber of about 2 10"3 and
that a shock wave is located in the divergent part. The
1.5 10"4, respectively. After deposition, the crucibles with
flow at the nozzle exit remains thus subsonic in a wide
the remaining powder are removed from the furnaces,
range of experimental conditions. The flow properties at
after having cooled down to room temperature, and
stored until the next utilisation. Vaporization conditions that location were then used as input boundary conditions
for both furnaces were investigated, and charts of for finite element calculations [27] in the deposition
chamber, neglecting the presence of Zr and Y precursors
precursor mass loss vs. furnace temperature and argon
flow rate were drawn. (less than 1 mol%), as illustrated on Figure 3. For some
calculations, the heating of the gas by the microwave
Although both precursors are moisture sensitive, they can plasma was taken into account, but the numerous
be stored for several months in appropriate dessicators, chemical reactions induced by the plasma were up to now
without significant variation of their properties. To neglected. These simulations will be used to optimize the
prevent corrosion of metallic parts by hydrochloric acid deposition process, with respect to both deposition
(produced by the reaction of moisture with chlorides), the efficiency and coating homogeneity on complex-shaped
whole apparatus is flushed with argon after each substrates.
Table 1 - Experimental conditions for the deposition of
Zr02-Y203 TBC coatings.
Standard
Parameter Range
conditions
02 flow rate (seem) 25 - 500 200
Ar total flow rate (seem) 50 - 1400 800
for Zr furnace (seem) 20 - 700 100
for Y furnace (seem) 10-140 50
Zr furnace temperature (°C) 25 - 250 210
Y furnace temperature (°C) 25 - 250 170
Total pressure in deposition
chamber (Pa) 13-1330 133
Microwave input power (W) 0 - 1200 750
Substrate temperature (°C) 300 - 750 400
Nozzle exit - substrate
distance (mm) 25-90 50
Deposition duration (h) 0.25 - 3 1
Figure 2 - Schematic view of the gas injection nozzle.
11-4
XJ
B
§■
Table 2 - Measured specific mass peXp and pore volume diffusion becomes more important, by ca. one order of
fraction P, for various substrate temperatures. magnitude between 400°C and 750°C if a leV activation
energy is assumed, so that less porous structures are
Substrate temperature Pexp P allowed to grow. At high temperatures, bulk diffusion
(°Q (kg/m3) can also contribute to the elimination of porosity.
400 2760 54% The porosity is not homogeneous throughout the coating,
600 3500 42% the first micrometers appearing to be more dense.
Annealing at 1100°C in air for 24 hours produces the
750 4420 26%
coalescence of the intergranular porosity which becomes
observable by SEM. This sintering takes place
üü^i throughout the coating, as seen on Figure 6 (bulk of
coating and column base).
The deposition of highly porous coatings could have
been interesting with respect to lowered thermal
conductivities, if the obtained morphology had been
thermally stable up to temperatures corresponding to gas
turbine applications. This is clearly not the case here.
3.2 Deposition rate
The deposition rate is about 55 um/h (48 mg/h) in the
standard conditions. Figure 7 shows that the deposited
zirconia amount (here at Ts = 400°C) is proportional to
the vaporized precursor mass and thus to the partial
pressure of zirconium compounds in the gas phase.
Figure 8 shows that the substrate temperature has very
little effect on the deposition rate. This indicates that the
process is not controlled by surface kinetics. Further
experiments are necessary to determine the limiting step
of the deposition process (mass transport, gas phase
reactions, ...) for various deposition conditions. Note that
the analysis of the deposition mechanism is more
complicated than for conventional thermal CVD, as the
influence of the plasma must be taken into account.
The deposition rate depends both on the gaseous
precursors flow rates (at least for ZrCl4) and the overall
deposition efficiency. The former is related to the
vaporization capability of the ZrCl4 and Y(thd)3 furnaces.
The ZrCL, furnace has been tested up to mass flow rates
corresponding to deposition rates of 150 um/h (about 200
mg/h Zr02 - Y203 deposited on the substrate, at
Figure 6 - SEM cross-sectional view of a coating Ts = 750°C); up-scaling of both furnaces should not
deposited in the standard conditions (Ts = 400°C) and present any problems. The latter is related first to the
annealed for 24 hours at 1100°C in air : bulk of coating physical and chemical processes in the gas phase and at
(top) and substrate-coating interface (bottom). the gas/coating interface and secondly to the gas injection
conditions and to the geometrical configuration of the
activated particle surface diffusion [29]. In the present deposition chamber, i.e. nozzle geometry, substrate size
case, the time interval between ad-particles arrival is and shape, nozzle-substrate distance, etc. Hence, the
about 20 ms, being defined as c/2R with c the lattice deposition efficiencies, defined as the number of
parameter (c«5 10"'° m) and R the deposition rate. At deposited Zr moles to the number of vaporized Zr moles
low deposition temperatures, the distance travelled by the ratio, that can be calculated from Figure 8, i.e. 8-10%,
adsorbed particles by surface diffusion during that time apply only to the current non-optimized deposition
interval is too small to allow the construction of a configuration. For instance, the deposition rate for a
« perfect» crystal so that many voids (vacancy clusters nozzle exit - substrate distance of 75 mm is the half of
or pores) are incorporated into the coating. It is seen on that at 50 mm : this is due to the fact that the flow
the SEM coating cross-sections (Figure 4) that undergoes a conical expansion in the nozzle and in the
intercolumnar porosity cannot account for the whole deposition chamber. An efficiency of at least 20% would
measured 54% pore volume fraction, showing that be achievable after process optimization.
intracolumnar porosity exists which cannot be resolved Most coatings were deposited at rates lower than about
by SEM. At higher deposition temperatures, surface 80 um/h, because beyond that value thickness
11-6
70 20 f)
75 26 O 76
Figure 9-XRD diagrams of as-deposited PECVD Zr02 - Figure 10 - {400} region ofXRD diagrams for a ZrÖ2 -
XY2O3, with x = 7wt% (top), Uwf/o (middle) and 7wt% Y2O3 coating. Experimental diagrams and
16wt% (bottom). AI2O3 peaks are due to spoiled areas calculated contributions of the « t'» and « c » phases :
(edges) on an alumina substrate. as-deposited (top), annealed at 1200°C in air for 24
hours (middle) and 120 hours (bottom).
coatings, which appears after about 250 hours at 1200°C. substrates, with high rates (up to 150 um/h) at
This different behaviour of the PECVD sample may be temperatures ranging from 400 to 750CC. ZrCl4 and
explained by the high porosity level of the coating Y(thd)3 were used as source materials; vaporization,
deposited at low temperature (7^ = 400°C): the high handling and storage procedures have been established
amount of internal surfaces thus generated may favour, for both precursors.
comparatively to volume diffusion, the surface diffusion
of Y3+ cations, which are responsible for the evolution of Morphology, impurity content and crystallographic
structure of the coatings were investigated. At all
the «t' » metastable system towards thermodynamic
equilibrium. deposition temperatures, the coatings exhibited columnar
morphology and {100} crystallographic texture. At low
4. CONCLUSION substrate temperature (400°C), the pore volume fraction
was high {ca. 50%) and the coatings incorporated some
Plasma-assisted chemical vapor deposition has been carbon and chlorine impurities, but at higher temperature
successfully used to deposit thick (20 - 150 urn) zirconia- (750°C), the pore volume fraction reduced to about 25%
(5-20wt%) yttria coatings on metallic and alumina and little carbon and no chlorine were detected. It thus
11-8
appears that medium to high temperatures are temperatures », Surf. Coat. Technol. 68/69, 81-85 (1994).
unavoidable to synthesize suitable coatings for thermal [9] Kreutz E.W., Lemmer O., Wesner D.A., Alunovic M.,
barrier applications. Erkens G., Leyendecker T., Voss A., « Electron and laser
Metastable tetragonal «t'» phase was the major radiation as sources of zirconia film deposition », Surf.
constituent of coatings with 7 and llwt% yttria, while Coat. Technol. 74/75, 1005-1011 (1995).
16wt%-yttria coatings exhibited only the cubic phase, in [10] Movchan B.A., « EBPVD technology in the gas
accordance to the phase diagram. Annealing for 24 and turbine industry : present and future », JOM, 40-45 (nov.
120 hours of Zr02 - 7wt% Y203 coating deposited at 1996).
400°C showed that cubic phase with an yttria content of
about llwt% is allowed to grow, at the expense of the [11] Grünling H.W., Mannsmann W., « Plasma sprayed
« t' » phase whose yttria content diminishes to 5wt% and thermal barrier coayings for industrial gas turbines :
4wt%. This relatively fast evolution towards morphology, processing and properties », J. Phys. IV,
thermodynamic equilibrium may be related to the high Colloque C7 suppl. J. Phys. Ill, 3, 903-912 (1993).
porosity level. [12] MeVrel R., « Barrieres thermiques pour aubages de
moteurs aeYonautiques. Etat de 1'art et perspectives », La
Work is in progress to further investigate the applicability Recherche Aerospatiale 5-6, 381-392 (1996).
of PECVD to coat complex-shaped substrates and to
study the relationships between deposition parameter and [13] Meier S.M., Gupta D.K., « The evolution of thermal
coating properties. barrier coatings in gas turbine engine applications»,
ASME J. Eng. Gas Turbines Power 116, 250-257 (1994).
ACKNOWLEDGEMENTS [14] Osyka A.S, Rybnikov A.I., Leontiev S.A., Nikitin
N.V., Malashenko I.S., « Experience with metal/ceramic
This study has benefited from a partial financial support
coating in stationary gas turbines », Surf. Coat. Technol.
from the Delegation G6n6rale de l'Armement (DSP and
76/77,86-94(1995).
SPA6). Thorough XRD investigations performed by C.
Diot and T. Ochin, as well as fruitful discussions, have [15] Parker D.W., «Improving thermal efficiency with
been are greatly appreciated. ceramic thermal barrier coatings», Proceed, of the
Internat. Symp. on Developments and Applications of
REFERENCES Ceramics and New Metal Alloys, Quebec City, Canada,
aug. 29 - sept. 2, 1993, 603-613.
[1] Stecura S., «Optimization of NiCrAl/Zr02-Y203 [16] Minh N.Q., «Ceramic fuel cells», J. Am. Ceram.
thermal barrier systems », NASA TM 86905 (1985). Soc. 76(3), 563-588 (1993).
[2] Maleto M.I., Solovjeva L.I., Turesvskaya E.P., [17]«Thermal Barrier Coating Workshop '95»
Vorolitov K.A., Yanovskaya M.I., «Alkoxy-derived (Cleveland, Ohio, USA, mar. 27-29, 1995), NASA-CP-
Y203-stabilized Zr02 thin films », Thin Solid Films 249, 3312(1995).
1-5 (1994). [18] «Thermal Barrier Coating Workshop 1997»
[3] Yamane H., Hirai T., «Yttria stabilized zirconia (Cincinnati, USA, may 19-11, 1997)
transparent films prepared by chemical vapour [19] DeMasi-Marcin J.T., Gupka D.K., «Protective
deposition », J. Crystal Growth 94, 880-884 (1989). coatings in the gas turbine engine », Surf. Coat. Technol.
[4] Kim J.S., Marzouk H.A., Reucroft P.J., « Deposition 68/69, 1-9(1994).
and structural characterization of Zr02 and yttria- [20] Rossnagel S.M., Cuomo J.J., Westwood W.D., Eds.,
stabilized Zr02 films by chemical vapour deposition », « Handbook of plasma processing technology », Noyes
Thin Solid Films 254, 33-38 (1995). Publications, ISBN-0-8155-1220-1 (1990).
[5] Cao C.B., Wang J.T., Yu W.J., Peng D.K., Meng [21] Seiberras G., «Realisation et caracterisation de
G.Y., « Research on YSZ thin films prepared by plasma- revetements de zircone obtenus par döpöt chimique en
CVD process », Thin Solid Films 249, 163-167 (1994). phase vapeur assist^ par un plasma miro-onde », These
[6] Brinkman H.W., Meijerink J., de Vries K.J., de Doctorat, University Paris 11-Orsay, mar. 24th, 1994.
Burggraaf A.J., «Kinetics and morphology of [22] Bertrand G., M6vrel R., « Zirconia coatings realized
electrochemical vapour deposited thin zirconia/yttria by microwave plasma-enhanced chemical vapor
layers on porous substrates », J. Eur. Ceram. Soc. 16, deposition », Thin Solid Films 292, 241-246 (1997).
587-600(1996). [23] Gavillet J., Belmonte T., Hertz D., Michel H, « Low
[7] Carter W.B., Godfrey S., «Combustion vapor temperature zirconia thin film synthesis by a chemical
deposited partially stabilized zirconia coatings», in vapour deposition process involving ZrCl4 and 02-H2-Ar
« Elevated temperature coatings : Science and technology microwave post-discharges. Comparison with a
I », Dahotre N.B., Hampikian J.N., Stiglich J.J. Eds., The conventional CVD hydrolysis process», Thin Solid
Minerals, Metals & Materials Society, 103-111 (1995). Films 301, 35-44 (1997).
[8] Andritschky M., Rebouta L., Teixeira V., « Corrosion [24] Chevillard S., «Etude d'un procedö de döpöt
and adherence of stabilized Zr02 coatings at high chimique en phase vapeur assists par un plasma micro-
11-9
L. T. Kabacoff
Materials Division
Office of Naval Research
800 N. Quincy Street
Arlington, VA 22217-5660
stabilized and exhibit no transformation dramatically reduced as the grain size goes
toughening. The toughening mechanism below 100 nm, the traditional definition of a
reappears if the percentage of yttria is reduced nanostructured material. However, the drop in
to below 2 percent (Ref 2). Another important thermal conductivity at operating temperatures
feature of nanostructured ceramics is reduced within gas turbine engines becomes significant
thermal conductivity (Ref 3). This effect only when the grain size goes below about 20
results from the grain boundaries being closer nm. Several unpublished studies have
-together than the typical distance between point demonstrated that this is indeed the case. Not
defects in a conventional ceramic. Thus, the shown in the calculations is the photon
grain boundary contribution to phonon contribution to thermal conductivity. This is of
scattering becomes anomalously high. The major importance since nanostructured zirconia
calculated phonon contribution to thermal is quite transparent to infrared. Therefore, a
conductivity is illustrated for zirconia (figure 2) nanostructured TBC would have to include
and YAG (figure 3). In both cases, the thermal micron size inclusions, either pores or a
conductivity at room temperature is suitable second phase.
>^
•a
u
3
T3
C -
o 2
u
Intrinsic
S-l
D
X!
5 nm
0
0 200 400 600 800 1000 1200
Temperature, T (°C)
Figure 2. Thermal conductivity of 7YSZ versus grain size and temperaure (Ref. 3)
12-4
10 j— i i i i
YAG
,__,
- 8
U
E
^ 6 A \
>->
*> \>lQpnmv\0
N
3 4
^ v^ o
N
^ \>
o
S v ^
^ »x "N-
^
^^w o o
U V* ^
"3
i 2 ■* *■
— —
- —.
5 nm
i i i i i
u
0 200 400 600 800 1000 1200
Temperature, T (°C)
Figure 3. Thermal con ductivity of YAG versus grain size and temperaure (Ref. 3)
12-5
of the sprayed powder during processing. The TBC. In spite of the lack of a large body of
grain size was limited by the quench rate and definitive data, two facts emerge.
the phase separation of zirconia and alumina Nanostructured ceramics deform
during solidification. The critical size for superplastically at moderate temperature (Ref
stabilization of rutile was larger in the 6). This temperature is typically less than half
composite than that found in individual pure the melting temperature and can be even lower
particles. Typical grain size in these studies for very small grain size. It has also been
was 30 to 50 nm. observed that intrinsic stress can be relieved at
low temperature, perhaps even room
The mechanical properties of nanostructured temperature. Room temperature relaxation has
ceramics depend critically on composition and definitely been observed in nanostructured
microstructure. Almost all of the studies have cermets (Ref 7). It was also observed in
been performed on "fully dense" material (in nanostructured titania by Gleiter (Ref 8).
reality, unintentionally porous). Very little Gleiter attributed this to very high room
work has been done on intentionally porous temperature diffusion (more than six orders of
materials, as one would prefer to have in a magnitude higher than that found in microscale
400-
ZrOr 3mol%Y203 +
+ commercial submicron
(30 nm starting grain size;
200 nm agglomerate size) /+
300- ■ G47% nano
6 D G58% nano
(nano: 15 nm starting grain size;
• 15-30 nm agglomerate size)
200- A
Q
■I A
Ö 100- 97.8%
-56.7%
. E 7^82.2%
5 °C/min heating rate
0 2 hr. hold at tempr
, . , —i ■ r -
600 800 1000 1200 1400
materials). It is clear that these coatings ceramic powder does not exist. Several
should not be used in situations where creep companies have the capability to provide
resistance is required. Another property production capacity, but will need to see an
affected by the ultrafine microstructure is adequate market before doing so. Research
fracture toughness. In the case of thermally quantities are easily obtained from numerous
sprayed nanostrucured zirconia, the biggest sources.
effect may depend on the way the coating is
formed. If the feedstock consists of Development of techniques for fabricating
agglomerated small particles, such as colloidal coatings by thermal spray are at a very early
material, and no melting (or only partial stage. It was noted above that improvement of
melting) occurs, then the coating will lack the thermal resistivity will require an extremely
splat boundaries present in conventional small grain size, of the order of 10 - 20 nm,
coatings. These splat boundaries are a major while improvement of mechanical properties
factor in the failure of conventional thermal can be achieved by reducing grain size only to
sprayed coatings. It is possible to form about 30-70 nm. The difficulty in retaining an
coatings without melting because of the very ultrafine microstructure in the coating during
high sinterability of the small particles. processing an subsequent exposure to elevated
Moreover, use of such small particles is temperature increases dramatically with
necessary if one wants to achieve a grain size reduced grain size. One reason for this is that
small enough to effect thermal conductivity. If the driving force for grain coarsening is highly
melting does occur, the coating will still be grain size dependent. Another reason is that,
nanostructured provided the sprayed material is while it is possible to produce a 30 - 70 nm
a composite of immiscible phases. In this case, grain structure through melting and rapid
the thermal conductivity will not be greatly solidification, a 10 - 20 nm grain structure can
reduced, but the mechanical properties will still only be achieved by avoiding melting, relying
be enhanced. Several companies are pursuing instead on rapid sintering of very fine particles.
both strategies, but the information is This may, if fact, be quite feasible, but it has
considered proprietary at this time. yet to be demonstrated. Clearly, improving
mechanical properties will be much easier than
Discussion and Conclusions reducing thermal conductivity. In both cases,
control of coating microstructure will require
The first requirement for successful fabrication very tight control of the thermal history of the
of nanostructured TBC's be to obtain suitable sprayed particles. This, it turn, will require
material for thermal spray. The synthesis of excellent control of both the size and
suitable sprayable material has not been morphology of sprayed particles, and spray
difficult. Nanoscale ceramic powder can be conditions. It may, in fact, require real time
synthesized by gas condensation, by sol gel and process control at a scale not currently
colloidal processes, by various flame and available. Thus, it may be necessary to extend
plasma processes, and mechanical attrition. the state of the art in thermal spray. This much
Several techniques have also been developed more likely to be true for improvement in
for reprocessing raw powder into a suitable thermal conductivity than for improvement of
form for thermal spray. These reprocessing mechanical properties. Finally, a great deal of
techniques are proprietary and cannot be research still needs to be done in order to
disclosed at this time. At present, commercial determine exactly what composition and
manufacturing of large quantities of nanoscale microstructure to attempt to achieve. This will
12-7
References
Zirconia-based thermal barrier coatings exhibit already low As a general rule, the thermal conductivity X of a porous solid
thermal conductivities (typically between 1 and 2 W/m.K). derived with these models can be expressed as :
These low values of thermal conductivity arise both from the
intrinsic value of dense partially (and fully) stabilised zirconia X/Xo = f (s, nf)
materials and from the heterogeneous microstructure of the
coatings deposited by plasma spraying or EBPVD (electron where Xo is the thermal conductivity of the pore-free solid, s is
beam physical vapour deposition). the volume fraction of the pores and n( are adjustable
parameters, to be determined empirically, related to the shape
The introduction of a stabiliser, required to avoid the detrimental and/or orientation of pores.
effect of monoclinic to tetragonal phase transformation at typical
service temperature, is accompanied by the incorporation of a To evaluate the validity of these models, several authors have
substantial amount of vacancies. Thus, in the case of yttria examined the analytical results given by some of these models in
stabilised zirconia, the substitution of two tetravalent zirconium light of experimental values on different systems [19, 20, 21, for
ions by trivalent yttrium ions is accompanied by the example]. Examining experimental data on alumina, graphite,
incorporation of one oxygen vacancy. For PSZ containing 8wt% uranium dioxide, sandstone, silica brick and limestone, and with
porosity up to 75%, Rhee [19] found a satisfactory correlation and measured values. This is attributed to the complex pore
with the equation proposed by Aizanov and Domashnev [22]. structure, which cannot be properly taken into account by the
For El-Fekey et al. [20] who evaluated the Maxwell, Loeb and limited areas on the cross sectional micrographs serving as input
Murabayashi models on thoria compacts with porosities up to information into the calculation.
46%, the Murabayashi model [14] yields the best result among
the three. More recently, Jackson et al. [21] found that In a study on the thermal conductivity of U02 pellets, Bakker
McLachlan's equation [18] comes closest to modelling the [30] used a finite element method to compute the conductivity of
experimentally measured thermal conductivity of A1N-Ln203 a matrix containing a dispersed phase (pores). Starting from a
systems. photograph of a cross section, the microstructure of the material
is described as a triangular mesh used as input data. Arbitrary
Remarking that the experiments considered may not be suited to temperatures are imposed on the upper and lower boundaries of
determine the quality of theoretical calculations, due to various the corresponding area (other boundaries are taken adiabatic)
experimental inaccuracies and simplifications, Bakker et al [23] and the FEM program computes the conductivity from the
have computed the influence of pores of different shapes thermal flux profile. As explained later, it remains doubtful
(cylinders, spheres, ellipsoids) with the FEM technique. The whether this type approach can be extended to large areas,
comparison with the equations describing the effect of keeping memory storage reasonable.
inclusions on the overall thermal conductivity shows that the
analytical equation of Schulz [15] gives a good description of In the work presented here, the objective was to develop and
such materials. evaluate an alternative approach, based on a finite-difference
calculation for computing the thermal conductivity of thermal
A full review of all the models is beyond the scope of this barrier coating, using as input data digitised images of the real
article. However, the puzzling diversity of conclusions derived material. With this method, the morphology of the ceramic
by these authors point out to the fact that all these models with coating, as complex as it can be, is properly taken into account.
simple algebraic expressions assume that the porosity can be
experimentally assessed with a very limited number of 2. CHARACTERISATION OF TBCs POROUS
parameters (the pore volume fraction e and adjustable STRUCTURE
parameters), which singularly restricts their predictive capacity. In order to provide data to the thermal conductivity calculation
software, reliable procedures have been developed to quantify
To overcome these shortcomings, several authors have the real TBCs porous structure. The success of the modelling
developed more sophisticated theoretical approaches as well as strongly depends on the accuracy and reliability of the results of
finite element calculations to estimate the effective thermal these procedures. The complete microscopic characterisation of
conductivity of multiphase materials. the TBC morphology involves numerous steps as sample
preparation, image acquisition, image analysis and validation
Tzou [24] estimated the thermal conductivity starting from a with macroscopic porosity determination.
general field theory (supposing uniformly distributed internal
cavities), deriving tractable analytical expressions for simple Sample preparation
pore geometries, i.e. insulated spherical cavities and penny- A correct and reproducible metallographic preparation of plasma
shaped cracks. sprayed Zr02 coatings is not straightforward but is essential for
obtaining a representation of the real porous structure. Porosity
In a more sophisticated approach, Furmanski [25] developed an determination by quantitative image analysis and numerical
effective macroscopic description for heat conduction in calculation of the thermal conductivity will give reliable results
heterogeneous materials using an averaging technique and provided that no error is introduced at the sample preparation
Green's function method and applied it to describe a model stage.
composite with randomly oriented parallelepiped inclusions
distributed in an epoxy matrix. Materials
All investigations are performed on a 1 mm thick free standing
The expressions obtained with these approaches become readily TBC plasma sprayed by SNECMA, initially on aHastelloyX
complex and apparently they have been used up to now only in substrate without any bondcoat. The powder is a commercial 8
the case of uniform distributions involving rather simple weight% yttria partially stabilised zirconia spray dried powder
geometrical shapes for the pores (ellipsoid, cubes). For materials (HCST Amperit 827.423).
having complex pore microstructures, such as those met in
plasma sprayed coatings, and in order to be closer to the real Sample preparation
structure (interconnected porous structure), however, numerical The following procedure has been set up. The free standing
schemes appear to be the most promising approaches. coating is vacuum impregnated (Struers/Epovac, P-100 mbar)
with a superfiuid epoxy dye (Struers/Epofix 301) before
With the objective of extending the rather simple analytical sectioning. Sample slices are then cut off with a precision saw
approaches (and poorly predictive) proposed by McPherson [26] (Struers/Accutom-2) using a thin diamond wheel and the lowest
and by Moreau et al. [27] on the thermal conductivity of plasma rotation and sample forward speeds to avoid any damage. The
sprayed coatings (and later improved by Bjorneklett et al.[28]), slices are then manually polished according to the following
Hollis [29] developed a numerical scheme in which the actual steps : 1200 then 4000 SiC papers for grinding, polishing using
pore structures of vacuum (VPS) and air (APS) plasma sprayed a 3 urn diamond spray first on a hard cloth (Struers/Pan W, blue
tungsten coatings are used as the basis of finite-element models lubricant) then on a smoother cloth (Struers/DP-Mol, pink
to calculate the effect of pores on the thermal conductivity of lubricant), final polishing with a !4 urn diamond paste.
plasma sprayed coatings. If, for VPS coatings, whose pore Polishing quality and absence of pull-outs are checked with light
distribution and shapes are relatively simple, this approach gives microscopy after each step. For scanning electron microscopy
a reasonable agreement between calculated and experimental observations a -20 nm thick conducting carbon layer is vacuum
values (respectively 70% and 60% of the bulk value), for APS sputtered within a Balzers MED010 unit.
coatings, instead, a large discrepancy exists between calculated
13-3
Image acquisition
Cross sections of TBC prepared as previously described have JÄÄiBiteMSJ!ÄiÄiJÄ*«l
been observed in a digital scanning electron microscope
(Zeiss/DSM960). Backscattered electron mode is chosen
fill "J 1 T 1 ['
•-jsoa»
■
obtain a reproducible grey level histogram (Fig. 1). Two
I I 1 iJ «ax»
magnification levels (180x180 um2; 45x45 um2) and two image I'll
[ill 1 i,„
1 k^ ■ lOOOO
resolutions (512x512 or 1024x1024 pixels) have been used for .. 70 ?1 72 73 U ?6 ?S 77 ?8 7SJÖ *i S3 S3 &* 85 8g S? S3 83 SO
characterising the TBCs morphology. Their influence will be orgy twes
Image analysis
Image analysis procedures have been developed with two Figure 2: Thresholding.
objectives:
• provide binary images of TBC cross sections which will be
used directly as input data for the thermal conductivity Porosity separation
calculation software, This procedure is based on opening and reconstruction image
• quantify the porous structure in such a way as to be able to analysis operations. It creates three new images from the initial
determine the contribution of each morphological feature to thresholded image : one containing the dispersed globular pores,
the thermal conductivity of a thermal barrier coating. one with the crack network and one with the globular pores and
the cracks directly connected to them.
This will be the starting point for modelling thermal
conductivity as a function of the morphology. This fully automated procedure will be widely used to determine
the contribution of each morphological feature to the
Thresholding conductivity (see § 4). A still more precise description can be
The transformation of SEM grey level images into binary images achieved using the following procedure.
is tht key step on which depends the reliability of the subsequent
results. The threshold level is manually fixed on images with Porosity reconstruction
optimised and reproducible grey level histogram (Fig. 1). From an initial parent binary image it is possible to derive a
series of images (about 300 or more according to the
This last one contains two peaks: a very sharp and intensive morphology complexity) with evolving porosity. This is done by
peak at the 0 level (black level) corresponding to the majority of applying a sequence of erosion image analysis operations to the
the globular pores and a wider one centred at middle grey levels parent image which leads to a final image containing only
which includes cracks and zirconia matrix. «ultimate erodeds ». These pixels exactly correspond to the
globular pore centroi'ds. Starting from this new image, a step by
step reconstruction of the porous features creates the derived
images. The operation is ended when the initial total porosity is
«MCi—T i
j completely reconstructed (Fig. 5).
I i 1 1
1 1 1 An application of this procedure in relation with the thermal
j f
conductivity computation will be shown in the last section.
rl
1 zJ1
i
I
40000
Image assembling
30C00 1 This algorithm is developed to provide large size images of the
[ 1 samples which are more representative for the conductivity
20000
11 I calculation. This is done by an automatic search of maximal
tOCOO
a
20 40 60 80
LL _XJ_J
100 120 140 160 1W 203 220 24Q 260 ■
overlap between two binary images from contiguous and slightly
overlapping fields.
Gtey lev& |
|0*tilavf|j«Jrpx^|
Reconstruction to globs
Figure 3: Globular pore distribution (% of particles
versus area")
• •
Analysis of the crack network
This procedure is based on the «skeletonization» (image
analysis terminology) of the crack network, followed by a
decomposition in juxtaposed segments of minimal fixed length.
A distribution of these segments with regards to their orientation
is then available (Fig. 4).
- *
Reconstruction to globs and linked cracks
i e 5 S Si 5*?'
I i ! 3 3 3 ! ! 3 i
3. NUMERICAL COMPUTATION OF THE THERMAL properties of each node (e.g. the local thermal conductivity) and
CONDUCTIVITY are stored as coefficients of the partial derivative equation to be
It has been shown that the morphology of a zirconia plasma solved. The schematic representation of our problem is shown
sprayed coating is relatively complex in so far as it involves on figure 6:
dispersion of pores with various sizes and shapes and connected
crack network. As we are concerned with the influence of
morphology on thermal conductivity and we want all the
morphological information obtained with SEM to be taken into
account, the modelling method must face this structure
complexity without oversimplification.
Specific implementation of the finite difference method Four first order distinct schemes can be constructed on this basis
The grid construction may be reduced to a very simple process if (the four corner implementations). It can be shown that they lead
we use a regular uniform mesh grid where nodes are simply to a space transformation in the neighbourhood of
defined by the pixels (square pixels) of the input image. With discontinuities as they operate a small translation along bisectrix
such a grid no geometric internal boundary information is of ij axes. For random geometries they are not strictly
associated with meshes. These information are related to the equivalent as they produce small variations (about 10"4 or less)
13-6
of the resulting computed conductivities, but as we will see They can be implemented without extra storage so our criterion
below (Fig. 9) this is not the most critical source of uncertainties is satisfied. Provided that we are able to find some well suited
in the present problem. (Fig. 7) convergence acceleration parameter (over-relaxed
Gauss-Seidel method), they can outperform matrix methods
For the main boundaries (e.g. the boundaries where Dirichlet or when applied on scalar or superscalar computer [33], but they
Neumann conditions are applied) we use the standard method, are notably less efficient when deep vectorisation is allowed.
that is the image method based on the same scheme for
Neumann condition (the Dirichlet condition does not need any So we have implemented an over-relaxed Gauss-Seidel solver
computation as T is imposed). which allows the solution of problems with sizes up to about
6000x6000 on our 512Mo RAM workstation. With the
When the temperature is known for each node, the average or commonly used SEM magnification this theoretically allows the
equivalent thermal conductivity of the whole domain is computation of TBC conductivity for thick samples (1.2x1.2
computed from : mm2), or extended areas of relatively thin coatings (for example
0.20x36mm2).
\A—-dx
>L = Gauss-Seidel over-relaxation factor
/ Ti-T0
The solver
The proposed scheme leads to a set of Nx-2*Ny linear equations
to be solved. As we are concerned with large values of Nx and
Ny (computation with Nx~4200 and Ny~5100 has been
performed), the solver has to be as efficient as possible.
With the definition of A,B,C,D we have seen before, it can For huge problems, the main limitation is the duration of
easily be shown that, for a simple diphasic problem, we have to computation. The following table shows computation duration in
store only 23 + 22 + 23 different configurations. So the storage for seconds for a standard 1024x1024 problem with 10"4 precision
A,B,C,D can be drastically reduced provided that we can store criterion (e.g. the convergence criterion on heat flux in relative
some case index for each node. form), for three different superscalar computers :
As this case number is in the range of tens, an array of bytes is HP9000/780 DEC Alpha 500 Pentium Pro
convenient. So the minimal storage required is : 260 405 1164
• one Nx*Ny double precision array for T,
• one Nx*Ny byte array for case markers, The computation was achieved with 3000 iterations using a
• 4*20 double precision scalars for coefficients, convergence acceleration parameter of 1.997. The maximal
• one Nx*Ny byte array for input data. residual error on heat flux (relative value) was about 9.4.10"5.
This is about 10 Moctets for a 1024x1024 problem and up to One can consider this case as the simplest one because the input
210 Moctets for a 4200x5100 problem. So if the solver does not data image was of medium complexity. For higher complexity
need intermediate storage, we expect that such problems can be input the maximal computation time may be twice the indicated
solved on common 256Mo RAM computers. The solution time. For the 4200x5100 problem we have solved only on the
method may be chosen with respect to this criterion. The two FTP workstation, the computation time was about 10 hours in the
common approaches for this kind of solver are iterative matrix most difficult case. Notice that the complexity (which can be
methods and iterative by points methods. evaluated as some function of the number of porous objects)
increases roughly as the area of the domain. Despite that, the
The first kind is well represented by the conjugate gradient relation between the size of the domain and the computation
method [32] which, as other matrix methods, does not satisfy time which is a power function of the maximal dimension of the
our criterion. It needs, at least, two or three times the minimal domain exhibits an exponent close to 3, which is the theoretical
storage defined above. So the maximal dimensions of the value for the over-relaxed Gauss-Seidel method.
problems we could solve with such methods would be reduced
by a factor of V2 or V3. Alternatively, we would need some The initial value for T may also strongly influence the speed of
supercomputer.Therefore, matrix methods are very efficient on convergence. Our standard initialisation is the temperature
computers with high vectorisation capabilities (as Cray for distribution for the homogeneous medium. A possibly more
example), but seems to be less attractive on scalar or superscalar efficient initial distribution may be constructed from the
computers (such as common workstations) [33]. Richardson's extrapolation concept [34]. Let Tm be the solution
The second kind is commonly known as the Gauss or Gauss- of the problem for the domain reduced by a factor 2 in each
Seidel method [32]. dimension (e.g. we compute only one node out of two). This
13-7
K - Ki + N
For the porous one, we can use the thermal conductivity of air
Such representations performed on images with separated which is well known. The major problem is for the zirconia
porosity give a qualitative useful description of the effect of each matrix thermal conductivity. At present only bibliographical
morphological feature on the heat transfer through a TBC. values are available. They fall between 2 and 3 W.nf'.K"1 [5], so
the related uncertainties on computed thermal conductivities are
4. APPLICATIONS large. Some in-situ evaluations of the matrix conductivity are in
Before applying these tools to various TBCs, three major points progress by means of microscopic thermal diffusivity
have to be considered. measurements on real TBCs. They may lead to more realistic
values of matrix properties in a near future. In the meantime, we
The resolution problem have used an arbitrary value of 2 W.itf'.K"1 for all
The first one concerns the optimal choice for input image computations.
resolution. As figure 9 shows, the computed conductivity
strongly depends on the input image dimensions.
The 2D to 3D effects
Images of the same physical field with various resolutions (e.g. The third point is related to restriction on actual geometry. All
256x256, 512x512, 1024x1024, 2048x2048 pixels) lead to four the data and results we are able to obtain with these methods are
inputs which have been used for conductivity computation. An
13-8
two-dimensional. Three-dimensional effects can heavily change correspond to the thermal conductivity computed from the
the conductivity if inclusions have geometrical properties which global 4234x5148 image resulting from the assembly of the 35
cannot be evaluated with two-dimensional sections. As Bakker fields according to §2.
[29] mentionned, the three-dimensional problem is unreachable The first important result is that the average values of the 35
with current computers. Therefore the construction of the fields computed conductivities are very close to those obtained
relevant three-dimensional input data sets would be for the global field as it can be seen on figure 10.
impracticable with the resolution we use. So the present work is
limited to the two-dimensional approach. Fractographies Thermal conductivity computations
MT70A serf« (cracks) <>,, 2.0 W.m"'.k ')
performed on plasma sprayed TBCs suggest that the cracks,
which are responsible of the main conductivity reduction, have
depth/width ratios far exceeding one and a two-dimensional
section is an acceptable representation for such geometries. But
this is not the case for globular porosity.
Muls> [.
Owing to these various uncertainties the computation results
have to be considered indicative rather than exact.
Results
These tools were first applied to a free standing TBC plasma
sprayed by SNECMA. Extensive analysis and computations
were done on such coating in order to :
Figure 12: Crack contribution.
• test and validate the different procedures,
• study the influence of different classes of porosity on
thermal conductivity. Therefore the decomposition in small (180x180 urn2) domains
with somewhat arbitrary boundary conditions leads to a correct
estimation of the equivalent conductivity of a larger element of
Thermal conductivity computations
MP70A strie (Ai, 2.0 W.m"'.Kr') TBC.
1.9 •
1.7
»obi
1.3
Figure 10: Plasma sprayed TBC computed thermal
Lilted
conductivity versus porosity. 1.1 Creeks
Cracks \
A set of 35 adjacent (1024x1024 pixels) fields was obtained
0.000 0.023 0.050 0.075 0.100 0.125 0.1:
with the SEM. The thermal conductivity has been computed on £
each field and is represented versus porosity on figure 10
Figure 13: Influence of porosity class.
Thermal conductivity computations
MP70A serfe (gobukr) (*.„ ZO W.nf'.K."')
5. CONCLUSION
A software has been developed for the computation of the 12. Rüssel H.W., J. Am. Ceram. Soc. 18, 1939, pi.
thermal conductivity of porous coatings (plasma sprayed TBCs)
from binary images of real material cross sections. This 13. Bruggeman D.A.G., «Berechnung Verschiedener
approach based on a finite difference method takes directly into Physikalischer Konstanten von Heterogenen Substanzen », Ann.
account the actual complex morphology of the ceramic which is Phys. 24(7), 1935, pp. 636-679.
mostly original with regard to existing models. This numerical
model is able to determine the contribution of each 14. Murabayashi M. et al, « Effect of porosity on the thermal
morphological feature to the thermal conductivity and therefore conductivity of Th02», J. Nucl. Sei. Tech. 6,1969, p. 47
compare different microstructures or coating architectures. In a
near future it will be associated to a morphology generator for 15.Schulz B., in High Temperature-High Pressure 13, 1981, pp
building a predictive tool. This one will be used in particular by 649-653.
engine manufacturers as guideline for modifying coating
spraying conditions in order to obtain specific morphologies 16. Koh J.C., Fortini A., «Thermal conductivity and electrical
leading to optimised coating thermal properties. resistivity of porous material », NASA report NAS3-12012, CR
It is important to note that this approach, developed on TBCs is 1200854(1971)
most general and can be applied to a variety of multiphase media
(refractories, composites, etc.). 17. Cunningham M.E., Peddicord K.L., «Heat conduction in
spheres packed in an infinite regular cubical array », Internat. J.
Acknowledgements : Part of the work reported has been carried Heat and Mass Transfer 24, 1991, ppl081-1088.
out within a Brite Euram project [4] and the support of the
European Commission is gratefully acknowledged. Special 18. McLachlan D.S., Blaszkiewicz M., Newnham R.E.,
thanks are due to SNECMA who provided the plasma sprayed « Electrical resistance of composites », J. Am. Ceram. Soc.
coatings. 73(8), 1990, pp 2187-21203.
7. Kingery W.D., Bowen H.K., Uhlmann D.R., in 25. Furmanski P., « Effective macroscopic description for heat
« Introduction to ceramics ». J. Wiley (NY, 1976). conduction in heterogeneous materials», Int. J. Heat Mass
Transfer 35(11), 1992, pp 3047-3058.
8. Madarasz F.L., Klemens P.G., « Reduction of lattice thermal
conductivity by point defects at intermediate temperatures », 26. McPherson R., « A model for the thermal conductivity of
Internat. J. Thermophysics. 8(2), 1987, pp. 257-262. plasma-sprayed ceramic coatings», Thin Solid Films 112,
1984, pp 89-95.
9. Maloney M.J., Achter H.S.,Barkalow B.K., « Development of
low thermal conductivity thermal barrier coatings». 27. C. Moreau, S. Boire-Lavigne, R.G. Saint-Jacques, «The
Communication presented in Thermal Barrier Coating relationship between the microstructure and thermal diffusivity
Workshop (Cincinnati, may 1997). of plasma-sprayed tungsten coatings», in Proc. of the 7th
National Thermal Spray Conference 20-24 June 1994, Boston
10. Maxwell J.C., «Treatise on Electricity and Magnetism». (MA). Edited by C.C. Berndt and S. Sampath (ASM, 1994), pp
Oxford University Press. 621-626.
11. Eucken A., « Thermal conductivity of ceramic refractory 28. Bjorneklett A., Haukeland L., Wigren J., Kristiansen H,
materials : calculation from thermal conductivity of «Effective medium theory and the thermal conductivity of
constituents », VDI-Forschungsh. 353 (3-4), 1932, pp. 1-16. plasma-sprayed ceramic coatings », J. Mat. Science 29, 1994, pp
4043-4050.
13-10
30. Bakker K., « Using the finite element method to compute the
influence of complex porosity and inclusion structures on the
thermal and electrical conductivity », Int. J. heat Mass Transfer
40(15), 1997, pp 3503-3511.
reliability of the conductivity data is assessed in sprayed condition and on the detached coatings
terms of the sensitivity of the laser flash in the heat-treated condition. Measurements
technique to the uncertainties in the coating and were made from room temperature to a
the substrate parameters. temperature of 1000°C. For the high
temperature measurements, the samples were
heated in a vacuum chamber in 100°C steps; a
2. EXPERIMENTAL PROCEDURE thermal diffusivity measurement was made at
each step. Specific heat measurements of A1203
Plasma-sprayed coatings were obtained by and 8YSZ were also made using a standard
spraying 8YSZ and A1203 onto 3 mm thick Rene Perkin-Elmer Model DSC-2 Differential
95 superalloy substrates measuring 62.5 mm X Scanning Calorimeter with sapphire as the
12.5 mm. Prior to deposition, the substrate reference material. Powders, scrapped from the
surfaces were grit blasted to improve coating as-sprayed coatings were used for this purpose.
adherence. Powders of 8YSZ (Metco 204NS; The standard and the sample were subjected to
average particle size: 10 urn) and A1203 the same heat flux as a blank and the differential
(Plasmalloy Al-1010; average particle size: 5 powers required to heat the sample and standard
Jim) were used. A Plasma Technik Spray system at the same rate were determined. From the
with a single spray nozzle and dual powder feeder masses of the sapphire standard and the sample,
was employed. Both monolithic and multilayer the differential power, and the known specific
coatings involving alternating layers of A1203 heat of sapphire, the specific heat of the sample
and 8YSZ were prepared. The coatings was computed. The thermal conductivities of
exhibited residual porosity. Porosity levels were the coatings were determined using the
determined using measurements of coating mass relationship:
and volume as well as by image analysis and k = ocCpp ...(1)
point counting techniques on optical in which k is the thermal conductivity, a is the
micrographs. To determine if residual porosity thermal diffusivity, Cp is the specific heat and p
could be closed by sintering, heat treatment of is the density of the coating. X-ray diffraction
the coatings was performed at 1300°C for 50 hrs analyses of the coatings were performed using
in a furnace under flowing argon. For this Siemens D5000 unit, with CuKa radiation, to
purpose, 10 mm2 size samples, cut using a
diamond wafering blade, were employed. The identify the phases.
coatings detached from the substrate as units
after the heat treatment. The surfaces of the 3. RESULTS:
detached coatings that corresponded to the PLASMA SPRAYED COATINGS
coating-substrate interfaces were
metallographically polished to remove the metal (a) Monolithic Coatings:
oxide layers formed during the heat treatment. Microstructures of the as-plasma sprayed
EB-PVD coatings were made by electron beam coatings are shown in Fig. l(a&b); those after
evaporation of high purity A1203 and 8YSZ heat treatment are illustrated in Fig. 2(a&b).
sputtering targets on to CMSX-4 single crystal The coatings exhibited a porous structure,
superalloy substrates. The substrate temperature typical of air-plasma-sprayed TBCs. A
was about 300-500°C during the deposition comparison of Figures 1 and 2 reveals little
process. Both monolithic and multilayer change in porosity in A1203 and YSZ layers after
coatings involving alternating layers of A1203 heat treatment. Therefore, porosity levels in
and 8YSZ were prepared by this technique. the as-sprayed condition were assumed for the
heat treated coatings.
Thermal diffusivity measurements were made
using the laser flash technique [10]. The laser Table I. Data on Monolithic Coatin es
flash technique involves heating one side of the Property A1,0, 8YSZ
sample with a laser pulse of short duration and Coating thickness (mm) 0.34 0.37
measuring the temperature rise on the other side Porosity (%) 19 12
with an infrared detector. The thermal Density of monolithic (g/cc) 3.9 5.74
diffusivity is determined from the time required Expt. coating density (g/cc) 3.17 5.06
to reach one-half of the peak temperature and a Calc. coating density (g/cc) 3.28 5.13
transient heat conduction analysis of a two-layer
body. From this analysis, the thermal diffusivity
of a single layer coating on a substrate can be The porosity, density and the thickness of each
independently determined. Measurements were coating are presented in Table I. As mentioned
made on the coatings with substrate in the as- above, porosity levels were determined by a
"mass/volume" (direct) measurement, or by
14-3
,Jft
J
y. «•*
200 fxm
p». *'*'$'.'±->'*--. V*.:-'- '.:'..'■* •».;
Fig. 1(a) Microstructure of PS A1203 in the as- ftV!
sprayed condition
1.4
1.2 <•».
Diff.-Rene95
tt
E
O)
1 e»
—D-Sp. Ht.-Rene 95 5
^ 200 fim -o--Sp.Ht.-AIO
2 3 0,8 4*<
-*--Sp.Ht.-8YSZ X
Q.
Fig. 1(b) Microstructure of PS 8YSZ coating in ,j, &•—t,—*—*
0.6 (0
the as-sprayed condition
0.4
0 200 400 600 800 1000 1200
Temperature (°C)
2.5
u 2 200 |nm
3
■o
C 1.5
o Fig. 6(a) Microstructure of PS AZ4 coating in
u
75 1 the as-sprayed condition
€ 0.5
0)
■c
0 _l_
200 400 600 800 1000 1200
Temperature (°C)
— Dense Al O
*S^^Ä^» 2
Dense 8% YSZ
3
AI203(19%P) + 8%YSZ(12%F|)
(Calc. From Fully Dense Data)
o AZ1 (As Sprayed)
100 I I | I I'" a AZ2 (As Sprayed)
A AZ4 (As Sprayed)
o AZ8 (As Sprayed)
6 Al 20+8%
3
YSZ
Fig. 7(a) Microstructure of PS AZ4 coating in (Calc. From Expt. Coating Data)
o
3
•o
c
5 1
75
E
0.1
200 400 600 800 1000 1200
Temperature (°C)
total number of splat interfaces in each coating interface thermal resistance [5, 11]. This effect
can be expected to be nearly the same. It due to internal interfaces can be appreciated
appears that due to this similarity, the thermal from the data shown in Fig. 9. The predictions
conductivity levels of multilayer coatings for a bilayer, calculated from the porosity-
differed little and thus can be predicted with corrected data of fully dense materials as well as
reasonable accuracy from the experimental data from the experimental data for heat-treated
of monolithic coatings. monolithic coatings using Equation (3), agree
well with the experimental data for multilayer
While the measured conductivity data for all the coatings. This suggests that heat treatment may
as-sprayed multilayer coatings showed a have eliminated the effects due to the imperfect
significant decrease with temperature, the data splat interfaces in the coatings and that the
after heat treatment (Figure 9) were largely primary microstructural factor that influences
independent of temperature. Several thermal conductivity in this condition is the
microstructural factors should be considered in porosity within the layers. However, since the
understanding thermal conductivity changes coatings detached after heat treatment, it was
after heat treatment. Porosity and thermal not possible to determine the relative
resistance at interfaces can significantly contribution of the coating-substrate interface
influence the thermal conductivity in solids versus that of the splat interfaces.
[4, 5, 10, 11, 15]. Since porosity levels changed
only a little, the effect on thermal conductivity 5. RESULTS: EB-PVD COATINGS
due to this change can be considered negligible.
On the other hand, interfaces between the splats Details of EB-PVD coatings, including total
in thermal sprayed coatings have been suggested coating thickness, individual layer thickness and
to reduce the thermal conductivity due to the density are presented in Table III.
Table III. Description of Coatings, Layer Thickness and Density Values for EB-PVD
Coatings with CMSX Single Crystal Substrates.
wm m
^itfr^iiski frl Vi
L'j!ry. »':-j6i
BSr-i-'lB'.'Ä
H&M
Al203 (Monolithic)
8YSZ (Monolithic) 30000
— -N203 (0.54) + 8YSZ (0.46) (series Calc I
AI2Oa (0.59) + 8YSZ (0.41) (series Calc
Al203 (0.35) + 8YSZ (0.65) (series Calc
8YSZ
1 layer each AI203/8YSZ
4 layer each AI203/8YSZ
-350 alt. layers:
*
Al2 0 3 /8YSZ
-1000 alt. layers: Al203 /8YSZ
V — *■---»-
—»
100
_!__!__!__!__,_
400 600 800 1000 1200
Temperature (°C) Fig. 16. X-ray diffraction pattern for the 8YSZ
coating
Fig. 14. A comparison of the experimental
thermal conductivities of EB-PVD coatings with First, it is of interest to see how the relative
the data for solid materials and the calculated data proportion of phases in YSZ, influenced by Y203
for multilayer coatings level, would influence the conductivity of YSZ.
In Fig. 17, the conductivity data of various dense
YSZ specimens reported in literature [16,20,21],
6. DISCUSSION: EB-PVD COATINGS as well as that of solid 8YSZ evaluated in the
present study, are compared. The data differed
First of all, the thermal conductivity data of significantly among them, owing to the variation
8YSZ coating is higher than the bulk 8YSZ data in the amount of Y203 stabilizer. For example,
by a factor of 1.5. In order to understand this increasing Y203 from 5.3wt.% to 9wt.%, a
discrepancy, some microstructural factors are to significant decrease in the conductivity level can
be considered. The solid 8YSZ and the coating be seen, over the entire temperature range.
deposited using it, were found to be not under the However, at llwt.% Y203, the conductivity
same microstructural condition. Whereas the appears to increase, relative to that of 8-9wt.%
solid consisted of a mixture of monoclinic and Y203 stabilized Zr02 specimens. This trend is
tetragonal phases (Fig. 15), the coating consisted consistent with the variation in thermal
almost entirely of tetragonal phase (Fig. 16). The diffusivity as a function of Y203 content, as
type of phase present in the YSZ microstructure observed by Youngblood et. al. [16].
is known [20] to affect the thermal conductivity
and therefore, this factor should be considered in 5
—o— 9YSZ-Hasselman et. al. [20]
the comparison of the monolithic material and E 4.5 --A--Solid 8YSZ-Present study A
the coating. —o—11YSZ-Mg. Elekt.[21]
/
-«-5.3YSZ-Mg. Elekt. [21]
4
-+- -5.3YSZ-Hasselman et. al. [20]
2000 —a--7YSZ-Youngblood et. al.
3.5 (Single Crystal) [16] /'
u
3 3
■o
C
O
Z> O 2.5
« 2
E
1.5
1
200 400 600 800 1000 1200
M. I
& PS
Temperature (°C)
nLyvAv-ju^wji~
60 80 100
Angle
Fig. 17. Thermal conductivity data for YSZ with
various amount of stabilizer content, compiled
Fig. 15. X-ray diffraction pattern for solid 8YSZ from literature
material
14-10
[22] and (iii) a change in defect structure with -50 -40 -30 -20 -10 0 10 20 30 40 50
composition [23]. In Fig. 18, the thermal % A in Coating / Substrate Parameter
conductivity comparisons employed in Fig. 14
have been revised using the solid 5.3wt% YSZ Fig. 19. Sensitivity of thermal conductivity of
conductivity data, since, this material entirely PS 8YSZ coating to the changes in
consisted of tetragonal phase. As can be seen, coating/substrate parameters
the agreement between the conductivities of solid
8YSZ and 8YSZ coating as well as between the For the PS 8YSZ coating, in the order of
experimental data and the calculated data for decreasing sensitivity, the coating conductivity is
multilayer coatings has improved. It is to be
noted that the variations between the sensitive to: substrate thickness (ts), half-time
conductivity levels of multilayer coatings appear (S0.5), substrate diffusivity (ocs), coating thickness
to be due to differences in the relative layer (tc), coating density (pc) or specific heat (Cpc)
thicknesses of 8YSZ and A1203 layers.
and substrate density (ps) or specific heat (Cps).
7. ANALYSIS OF SENSITIVITY OF For example, for PS 8YSZ coating (Fig. 19), a -
THERMAL CONDUCTIVITY 10/+10 change in the parameter caused changes
MEASUREMENT in the coating conductivity as: -24/+46 for
substrate thickness, +35/-20 for half-time, +22/-
13 for substrate diffusivity, -21/+21 for coating
14-11
thickness, -11/+11 for coating density or specific substrate diffusivity, coating thickness, coating
heat and +1/-1 for substrate density or specific density or specific heat and substrate density or
heat. specific heat. For example, a -10/+10 change in
the parameter caused changes in the coating
For the present PS coatings, the accuracy of the conductivity as: -44/+°« for substrate thickness,
conductivity data principally depends on the +105/-35 for half-time, +86/-29 for substrate
parameters that cause the most variability in the diffusivity, -22/+27 for coating thickness, -
present coatings. Since substrate thickness could 15/+16 for coating density or specific heat and
be measured to the accuracy of 0.1 mm, the +5/-3 for substrate density or specific heat (Fig.
uncertainty in substrate thickness (0.03%) is too 20). The order of parameters in sensitivity
small to cause any change in coating ranking is similar to that of the PS coatings.
conductivity. The half-times are recorded using
digital oscilloscopes, so the uncertainty in their It is to be noted that for the EB-PVD coatings
measurement can be considered to be better than the variability encountered in coating parameters
1%. The substrate diffusivity was independently are about the same as that of the PS coatings.
measured using the free-standing sample and However, the EB-PVD coating thickness is a
hence the measurements can be considered factor of 4 smaller compared to the PS coatings.
accurate within ±1%. The density of the The substrate thickness was 1 mm. The sample
substrate were determined within ±1% using the used for thermal conductivity measurement was
Archimedes principle. The specific heat values prepared by mechanically grinding the substrate
of the coating and the substrate can also be from its original thickness of 3 mm to 1 mm.
determined to a level of ±1% or better. This After grinding, the average variation in substrate
leaves the coating thickness and the density as thickness was found to be about +/- 0.04 mm.
the most sensitive parameters involved in the This would suggest that the uncertainty in
present study. The coatings thicknesses were of substrate thickness is +A4%. Therefore, the
the order of 400 |J.m, with the coating surface increased sensitivity of conductivity data in the
roughness of the order of 40-50 urn, indicating case of EB-PVD coating, relative to PS coatings,
that the variability in coating thickness is about seems to arise from the reduced thicknesses of
10%. Additionally, the measurements of density the coating and the substrate.
of plasma sprayed coatings could be in error to
some degree, due to the penetration of liquid 8. CONCLUSIONS
medium through the open pores. Therefore, the
present thermal conductivity data for PS coatings Plasma Sprayed Coatings:
can be expected to be subject to these (1) The reductions in the thermal conductivities
uncertainties. of single layer as well as multilayer plasma-
sprayed A1203 and 8YSZ coatings are brought
100 about by porosity and thermal resistance at the
interfaces in the coating. The thermal resistance
appears to arise from the interfaces between the
splats as well as the coating/substrate interface.
The interlayer interfaces in the multilayer
appears to play only a minor role in influencing
the coating thermal conductivity.
i s pas» RPÄHRiwii. m\
0491S6 20KV 979321 29KV
Figure 1 Representative microstructures of TBC's applied by (a) plasma spraying in air (b) vacuum plasma
spraying and (c) EBPVD (S: substrate, B: bond coat, C: ceramic overlay layer, N: electroless Ni film used to
protect ceramic for metallography).
spraying (VPS) and electron beam physical vapour known to fail above this oxide layer. The mechanisms of
deposition (EBPVD) are shown in Figure 1. APS the reactions of the coatings with the environment in the
coatings (Fig. la) usually exhibit a layered structure, with engine and the ultimate coating failure are discussed in
a swirl pattern and a large amount of voids and pores. detail, elsewhere (8-10).
The pores help to increase the insulating capacity of the
oxide. VPS coatings (Fig. lb) are also layered but are Evaluating Thermal Barrier Coatings
denser than APS coatings and the interface between the The effectiveness of a TBC coating is assessed from its
bond coat and the ceramic layer is generally more ability to stay in place without spalling. The best way to
uniform than in air sprayed products. EBPVD coatings assess coating durability is from actual engine tests.
(Fig. lc) have a dense columnar structure, which can These, however, are expensive and time consuming.
accommodate the strains arising from differential thermal Alternatively, the performance of a coating can be
expansion between the metallic and ceramic layers. evaluated in a burner rig, where engine operating
conditions can be closely simulated. Engine
Failure Modes of Thermal Barrier Coatings manufacturers have used burner rig tests in the past to
The life of the overlay ceramic is limited by its evaluate and compare the performance of coatings for
susceptibility to thermal fatigue cracking, while the life of screening purposes, prior to actual engine tests. In
the bond coat is limited by its oxidation and corrosion Canada, military operational groups have also relied on
behaviour. Hot corrosion of the ceramic overlay may also burner rig tests to assess the potential of emerging
be life limiting in marine environments. The combined coatings for use in their engines, under various
effects of stress, time at high temperatures and oxidizing component retrofit programs. However, there is no
environment are complex, causing changes in the standard procedure established for burner rig tests.
microstructure of the coatings as well as the substrates. Manufacturers and users adopt different cycling times
These changes eventually will lead to failure of the and temperatures (heating times from a few minutes up to
coatings. Transverse and longitudinal cracks form at the an hour, at temperatures from 1000°C to 1500°C,
surface and inside the ceramic layer, due to thermally followed by cooling for a few minutes) tailored to suit
induced stresses. Oxidizing gases penetrate through the their end use.
transverse cracks and oxidize the bond coat. The layer of
A1203, which forms at the interface between the bond Objectives of the present study
coat and the ceramic (and possibly between the bond coat At NRC, several studies have been performed using
and substrate), thickens during service. In this layer, burner rig testing to compare TBC's produced by
oxidation products such as Cr203, NiO, Hf02 and spinels different techniques and different processing conditions,
(NiO-Cr203-Al203) are formed by loss of ß-NiAl from and to optimize some of the coatings for use in Canadian
the bond coat. In the case of EBPVD coatings, cracks Forces engines. The objective of this paper is to discuss
form in this layer of oxides, causing the top ceramic the testing procedures and methodologies employed in
coating to spall, whereas plasma sprayed coatings are these studies with a view to evaluating the suitability and
15-3
practicality of burner rig testing as a laboratory ranking of combustor chemistry are readily reproduced and
method for TBC's. The influence of key test parameters controlled. This is particularly important when assessing
on modes and rates of damage accumulation is discussed. the effects of gas/coating reactions on the durability of
Ways to ensure that the modes known to be life limiting coatings (14,15). Fuel (JP-4), along with compressed air,
in engines are reproduced during rig testing are also is introduced through a pressure-atomizing nozzle into
addressed. the primary section of the combustor, within which
combustion occurs in a swirl stabilized recirculation
2. EXPERIMENTAL MATERIALS AND zone. The expanding combustion gases are cooled by
METHODS dilution within a secondary section, and are exhausted
through a nozzle in front of which the test pieces are
2.1 Materials and Test Coupons placed for exposure to high velocity hot gases. By
The TBC's investigated in this work included independently adjusting the primary and secondary air
experimental coatings produced by EBPVD, APS and flows, along with the fuel flow, it is possible to vary the
VPS, under three development programs, as well as gas velocity and gas temperature independently of each
TBC's produced by established techniques by other. With a standard 2-inch axisymmetric circular
commercial vendors. In the first series of tests, the exhaust nozzle, gas velocity can be varied from Mach 0.2
performance of experimental EBPVD TBC's having to Mach 0.8 and gas temperature from 540°C (1000°F) to
compositionally and functionally graded interfaces 1650°C (3000°F). The latter is controlled by a
between the metallic bond coat and the ceramic layer, programmable microprocessor and measured by a
was evaluated. These coatings were applied to nickel base thermocouple located near the throat of the exhaust
superalloy pins, 0.25" diameter and 4" long, made from nozzle. Nozzles of different sizes can be used to vary the
IN625 and a proprietary single crystal (7). In the second cross sectional geometry of the hot gas jet to
series of tests, zirconia and calcia based TBC's were accommodate different test piece geometries. In the
plasma sprayed onto IN625 pins, 0.5" diameter and 5.5" present studies, the standard 2-inch diameter nozzle was
long, using different conditions of particle velocities and used when testing the Vi" diameter pins and a 3-inch
temperatures, and studied to assess the effects of oversized nozzle was used for the Vi" diameter pins. In
composition and processing conditions on product both cases, the pins were positioned 1.5" away from the
performance (11,12). In the third series of tests, the exhaust nozzle.
performances of a non-graded EBPVD TBC, TBC's with
VPS bond coats and APS ceramic overlays and TBC's The rig in the standard 2-inch nozzle configuration is
With VPS bond coats and VPS overlays, were also shown in Figure 2, with a set of uncoated Vi inch
studied. These latter coatings were applied to the same diameter pins positioned in front of the exhaust nozzle.
pins as those used for the first series of tests. The test pins are mounted on a rotating carrousel, used to
translate the pins back and forth between the hot gases
and a jet of room temperature air from a 2-inch diameter
pipe located beside the hot nozzle. In this way, the pins
can be alternately heated and cooled, simulating thermal
cycling in engines. Rotation of the pins, typically at 100
to 200 rpm, is used to promote gas mixing and ensure
uniform exposure of the pins to the hot and the cold
gases. An electrical slip ring assembly mounted on the
drive shaft of the rotating carrousel can be used for
instrumentation of 8 of the 12 test pins.
surface temperature distribution, as a function of gas while a close-up view of the spring loaded mechanism is
temperature and gas velocity and the changes in metal shown in Figure 6, with the test pins removed from the
temperature at the core of the pins, as a function of carrousel, but with the thermocouples still in place. The
number of thermal cycles, were established in this spring-loaded device ensured that the thermocouple beads
manner. For surface temperature measurements, holes were kept in contact with the bottom ends of the holes,
drilled radially and along the axes of the pins were used irrespective of the degree of thermal expansion and
as feedthroughs for the thermocouple leads, to avoid contraction experienced by the pins during thermal
disturbing the heat transfer boundary layer around the cycling. The position of the bottom of the hole was
pins. Several pins with one thermocouple each, and with chosen to correspond to the position within the hot gas jet
thermocouple beads regularly spaced along the axes of known to experience the highest temperature, in order to
the pins, were used for these measurements. A schematic record the maximum metal core temperatures.
cross section of the instrumented Vi inch diameter pins,
where all the bead positions have been combined in a
single drawing for simplicity, is shown in Figure 4. The
surface temperature measurements were performed on
bare pins.
r jr 7 j £-i h _T
Internal temperatures were only measured in the case of Burner rig tests involved up to 2150 cycles of 4 minutes
the Vi inch diameter pins. The intent was to track the of heating in the high velocity hot gas jet followed by 2
degradation in thermal protection efficiency of their minutes of forced air cooling. In some cases, the gas
TBC's as a function of number of thermal cycles, or test temperature was increased in a stepwise manner, from a
time, through measurements of drifts in heating and minimum of 1300°C to a maximum of 1500°C, to
cooling rates during thermal cycling. To measure the increase the severity of the tests and ensure that failure of
internal metal temperatures, thermocouples (one per pin) the coatings occurred before the end of the tests. The
were mounted with their beads pushed against the bottom surface conditions of the pins and their weight to the
of a blind hole drilled along the axes of the pins, using a nearest 0.5 mg were recorded before testing and after
spring-loaded mechanism. A schematic of the every 50 to 150 cycles. The weight changes provided a
instrumented Vi inch diameter pin is shown in Figure 5, measure of the spalling resistance of the various TBC's to
15-5
the combined effects of thermal loads and hot gas However, the coated pins reach higher steady state core
erosion. Changes in the heating and cooling rates of the temperatures than do the bare pins, as the rates of heat
pins during thermal cycling were also extracted on-line loss from the free ends of the coated pins are reduced
from internal temperature measurements. Finally, the relative to bare pins. Also, the internal temperatures of
microstructures of the TBC's were compared before and both bare and coated pins are higher than their surface
after testing, and the extent of thermally induced temperatures. This difference increases with gas
microstructural changes, including phase transformation temperature and is higher in coated pins, as shown in
as revealed by x-ray diffraction, was assessed. Table 1.
SFT3-150 O 1200
1400
' '"\ r
;
1400 "i" iooo *.$) JWlv' ■■■:^I t
1100 ü S^ i
1 ! ' I
I
O • 1500 ^V;
2^ 1000 3
S> N |600
<1400 • 1200 C Flame, PI- 800 P2* 50.0 • 1100CFI8TO. P1-47,5 P2-20.0
3 0)
a 900 Q. • 1200 C Flame, PI- 47.5 P2- 20 0 1300 C Flame. P1 - 47.5 P2- 20.0
(a) (b)
Figure 10 Profiles of steady state pin (bare) surface temperatures for (a) Vz inch diameter pins, and (b) %
inch diameter pins, measured at different gas temperatures and gas velocities.
exposed to the hot gases. The cracked surface of an spalling initiation varied with the deposition method, and
EBPVD coating on a Vi inch pin, after some thermal even more significantly, with processing conditions for
cycling had been applied, is shown in Figure 11. each of the deposition methods. Some coatings from all
four forms of TBC's investigated (EBPVD, APS,
VPS/APS and VPS) were found to survive temperatures
up to 1500°C for several hundred cycles, without
spalling. A series of APS TBC's, after rig testing are
shown in Figure 13. The 1.2 mm thick calcia based
coating, (pin (a) in Figure 13), shows no evidence of
spalling, whereas the zirconia based TBC's (all about 300
um thick, pins (b) to (f),) exhibit varying degrees of
delamination.
surface of the bond coat, as shown in Figure 14. The trends are due to the rate of bond coat oxidation
vertical cracks may relieve the stresses at the interface, increasing with temperature.
thereby reducing the driving force for delamination. One
of the pins, (b) in Figure 13, shows evidence of bond The data shown in Figure 16 also reveal that processing
coat oxidation, where the ceramic topcoat has been conditions during deposition can be adjusted to achieve
entirely eroded away by the hot gases. These observations improved performance with VPS coatings. The optimized
are generally consistent with expected modes of failure VPS coating showed minimum weight changes for more
for APS coatings (8-10). than 1000 cycles at gas temperatures up to 1400°C,
indicative of good resistance to ablation and spalling.
3.3 Damage Monitoring and Quantification However, it spalled off rapidly after the temperature was
raised to 1500°C, as did the APS coatings.
Weight Change Measurements
Changes in the weight of some coated V% inch pins,
measured over 1150 thermal cycles, are presented in
Figure 15. The weight losses observed for all the TBC's
are due to gas-induced ablation and spalling of the
topcoats, during the test. The thick APS calcia-based
coating appeared to have incurred the least damage in this
test, while the performance of APS zirconia-based
coatings varied greatly with processing conditions. The
relationships between processing conditions,
microstructures and properties of APS coatings are being
investigated at NRC and the results will be published
elsewhere (11,12).
Numbw of CyclM
E
O)
\ p~Vs \
Calcla bind
i
Figure 16 Changes in the mass of a series of
0)
01
c ■::,
V
\
TBC
1
-^
EBPVD, APS and VPS coatings applied to % inch
n ZlrconEa b* ml
1
pins as a function of number of rig thermal cycles.
TBC
\ The thermal cycling schedule is shown in the figure
\ "■-•'•
D)
Heating Rates Measurements were used in the present experiments tended to oxidize
The maximum heating rates experienced by coated !4 and fail by thermal fatigue. Type S thermocouples should
pins, as measured by thermocouples located inside the perform better in this application and ought to be tested.
pin's axial blind holes, are plotted in Figure 17 as a Also, the measurements should be made on pins that are
function of the number of thermal cycles. Each data point internally cooled to better simulate the thermal gradients
in Figure 17 corresponds to the maximum slope of the and levels of thermally induced stresses experienced by
rising portions of the temperature/time curves obtained TBC's applied to internally cooled components.
for each of the thermal cycles, examples of which are
shown in Figures 7 and 8. There are several features in CONCLUSIONS
Figure 17 that warrant comment. Firstly, the This work has shown that the modes and rates of
discontinuities in the data arise from interruptions in the degradation of TBC's during burner rig testing of coated
test sequence required to weigh the pins, following blocks solid round pins are influenced by many test parameters,
of 5, 50 or 150 thermal cycles. The initial sharp drop in besides being influenced by differences in deposition
heating rates at the onset of each block is due to the pins method and processing conditions.
being at room temperature at the onset of the first cycle in
a block, and at greater than room temperature for Pins, V* inch in diameter, heat up and cool down much
subsequent cycles, because the pins do not cool down faster than Vi inch pins. They also reach a higher
fully to room temperature at the end of a cycle (see temperature and spend more time at maximum
Figure 7). temperature over the duration of a thermal cycle. Thus,
coatings applied to Vt inch pins will be subjected to
TBC Efficency Breakdown higher transient and steady state stresses, and their bond
coat will oxidize more in the course of a test, as
compared to Vi inch pins.
i ♦ ►
/
12 V& Over the steady state portion of a thermal cycle, the
"
sir K
2 1 internal temperature of a solid pin is higher than its
2 11 *l V'
.«. ¥ t
surface temperature. This difference is larger in a TBC
w
coated pin than in a bare pin. High internal temperatures
► in coated pins cause high shear stresses near the interface
Zirco lia between the ceramic and metallic layers of the TBC,
owing to differential thermal expansion between the two
0 200 400 600 800 1000 1200
media, thereby promoting delamination near this
interface. Because of their highly strain tolerant columnar
Cycle Number
microstructure, EBPVD coatings generally fared better
than thermally sprayed coatings under such testing
Figure 17 Changes in heating rates of a V2 pin conditions. Using internally cooled pins, to better
coated with zirconia as function of number of simulate service conditions, would reduce the magnitude
thermal cycles applied to the pin, showing change in of these stresses and provide a more realistic laboratory
heat transfer characteristics of coated system. ranking of TBC coatings.
As the temperature at the start of a cycle increases, the The performance of each type of TBC examined, is a
apparent heating rate inside the pin decreases. It takes strong function of processing conditions. Analysis of
several cycles for heating rates to become stable because damage modes arising from thermal cycling in a rig
of the complex heat transfer transients affecting the pins confirms that plasma sprayed coatings spall by cracking
and their supporting carrousel. Secondly, the general of the ceramic overlay, above the interface with the bond
upward trend in the data, observed within each block and coat, and in the case of EBPVD coatings, by cracking
from block to block across the plot, may be viewed as an within the oxide layer that builds up at the surface of the
indication of the breakdown in thermal protection bond coat. The failure modes differ because of
efficiency of the TBC system. As the coating ablates, or differences in the degree of strain tolerance of the two
spalls off the test pin, the internal heating rates increase, types of coatings.
as expected, because heat is transferred more quickly into
the pins. Thirdly, there appears to be a threshold number ACKNOWLEDGEMENTS
of applied thermal cycles beyond which the increase in This work was performed at IAR and IMI under projects
heating rate accelerates sharply. This threshold is of the involving collaboration with a number of companies
order of 500 cycles, and corresponds to the point beyond (Cametoid Advanced Technologies Inc., Pyrogenesis
which significant weight loss was recorded for this Inc., Orenda Aerospace Corporation, SURFTEC
coating, as shown in Figure 15. Thus, monitoring the Industrial Consortium), with financial support from the
internal temperatures of coated pins could be used as an Chief Research and Development of the Department of
alternative method for the quantification of damage National Defense. Special thanks are due to Mr. Ray
accumulation in TBC's. However, further work is needed Dainty of IAR, for his support in SEM analysis.
to develop this technique. The type K thermocouples that
15-10
5. S. Durham, S. Manning Maier, D.K. Gupta and K.D. 14. A.K. Gupta, T. Terada, P.C. Patnaik and J-P.
Schefler, "Ceramic Thermal Barrier Coatings", in Immarigeon, "Evaluation of High Temperature
Advances in High Temperature Structural Materials Protective Coatings for Gas Turbine Engines Under
and Protective Coatings. Eds. A.K. Koul et al, Pub. Simulated Service Conditions", Proc. AGARD-
NRC, Canada, Ottawa, 1994, pp.226-236 NATO Conf. on High Temperature Surface
Interactions, Ottawa, May 1989, AGARD-CP-461,
6. Anonymous, "Vacuum plasma spray coating protects 1989, pp. 14-1/14-31.
turbine blades", a note in Advanced Materials and
Processes, March 1997, p. 12. 15. A.K. Gupta, J-P. Immarigeon and P. C. Patnaik, "A
Review of Factors Controlling the Gas Turbine Hot
7. F. Jamarani, M. Korotkin, R.V. Lang, M.F. Ouelette, Section Environment and their Influence on Hot Salt
K.L. Yan, R.W. Bertram and V.R. Parameswaran, Corrosion Test Methods", High Temp. Technology,
"Compositionally Graded TBC for High Vol. 7, No. 4, November 1989, pp. 173-186.
Temperature Aero Gas Turbine Components",
Surface and Coatings Technology, 54/55, 1992, pp. 16. S. Gendron, N.J. Marchand, C. Korn, J-P.
58-63. Immarigeon and J.J. Kacpryzynski, "An
Experimental Investigation of Convective Heat
8. M. Gell, "Application Opportunities for Transfer at the Leading Edge of a Turbine Airfoil",
Nanostructured Materials and Coatings", Materials Int. Gas Turbine and Aeroengine Congress and
Science and Engineering, Vol. A204, 1995, pp. 246- Exposition, Cologne, Germany, June 1992,
251. Publication ASME 92-GT-248, 1992
16-1
1 mm 10 um
Fig. 2: Cross section of a thick alumina coating Fig. 3: Surface of sputtered alumina layer on top
obtained by EB-PVD at low pressure ofNiCoCrAlYbondcoat
16-4
of deposition temperature on phase formation of alu- In Fig. 5 chemical composition is plotted as a function
mina is shown. of sample position above the ingots. Experiments
clearly show that dual source evaporation of two ceram-
XRD peaks were broadened for alumina films deposited ics with widely differing vapor pressures and melting
at low temperatures. This can be attributed to the pres- points is possible using the jumping-beam technique. A
ence of areas within the coating that were not fully mixing of both ceramic vapor clouds takes place and
crystallized. Coatings deposited at the highest tempera- both species condense on the substrates thereby forming
ture showed the presence of only the stable oc-Al203 mixed layers. While the portion of A1203 decreases more
phase. Additionally it was found that chamber pressure rapidly with decreasing distance to the ZrOz-ingot, the
and deposition rate had only little influence on phase portion of Zr02 rises slower and shows a flatter gradi-
formation of alumina which is contrary to results found ent. The zirconia content is predominant in all coatings.
for the morphology. The reason for that is the difference of vapor pressures
Phase transformation was also studied after annealing in and evaporation behavior between Zr02 and A1203 that
the temperature range between 1000°C and 1200CC for leads to a different geometry of the vapor clouds. Fur-
1 to 9 hrs. Annealing at 1000°C had little effect on any thermore it can be assumed that the re-evaporation rate
of the coatings; only a sharpening of diffraction peaks of alumina is higher than that of zirconia because of its
was observed. For all coatings, transformation to ot- lower melting and boiling point.
A1203 started upon annealing at 1100°C. Short term (1
h) annealing at 1100CC resulted in partial transforma-
tion to a-Al203. The higher the deposition temperature
was, the lower was the tendency for phase transforma-
tion. This can be explained by the fact that the driving
force for transformation (i.e. difference between coating
and annealing temperature) is lower. After annealing at
1200°C for 9hrs all coatings consisted mainly of oc-
A1203 with only small amounts of other phases /21, 22/.
Zirconia
A120, wt%
39
24
%uttr "Wir^-r Wirtr *HinrT"r KM c
H 1 1 1 1 1 1 1
253035 40 455055606570 75 80859095 100
2ttetain°
structural features are in good agreement with observa- calculation of the stress state effects like creep, plastic
tions on mixed layers described in chapter 3.3. deformation, oxidation, and sintering have to be taken
into account. It must contain also stresses after fabrica-
tion, in the transient stages, in the temperature field
during service and after cooling down to room tempera-
ture. All these aspects described above have to taken
into account in order to draw the right conclusions from
PYSZ TBC calculations.
In the following one special aspects of the design of
graded TBCs will be explained in detail.
R(ZJQZ)
deviations because of small temperature gradients.
Compared to a non-graded TBC of 200 urn Zr02 (curve
a) the heat flux of a linear graded coating with the same
thickness (curve b) is more than twice as high. In this
case the graded layer shows higher interfacial tensile
stress, which should decrease performance. Based on
(1) and (2) the overall thickness of the graded layer was
increased to same heat flux conditions (curve d), result-
ing in the lowest compressive interface stresses in the
ceramic but in a dramatic higher thickness too. Alterna-
tively, a 150 um thick zirconia layer was applied on top
of a 200 urn thick graded layer (curve c), resulting in
0.1 0.5 1 5 10 50 100 the same heat flux as the conventional system. Stresses
Profile Parameter p in the bondcoat are much lower than for all other sys-
tems, but interfacial stresses are higher compared to a
Fig. 12: Effective thermal resistivity for different
and d. From these examples one can see that the grading
rules of mixtures and profile parameters
profile significantly influences the stress situation. For
optimization, knowledge of critical stresses that are
this requirement, thickness has to be increased by a
factor of 4.6 in case of a linear graded layer compared to responsible for failure is essential. In case of coating
pure PYSZ. Increase of coating thickness is limited to a thickness of same heat flux the FGM concept offers a
certain level for aero engines. Therefore application of potential for improvement of coating performance.
an additional zirconia layer on top of the graded layer is
preferred to get same heat fluxes.
5. Conclusions
For a loading situation similar to a turbine blade stresses The potential of chemically graded thermal barrier
of different coating systems have been calculated (Fig. coatings to improve lifetime and reliability of turbine
13). In a biaxial stress state the radial stresses of an blades was evaluated. It has been demonstrated, that
infinite plate have been calculated assuming linear elas- graded ceramic coatings can be produced by dual
tic behavior /32. The model involves a 1mm thick su- source EB-PVD and jumping beam technology. Various
peralloy substrate that possesses a NiCoCrAlY bondcoat mixtures of alumina-PYSZ as well as pure alumina and
of 150 urn and the TBC on top. All material properties graded coatings have been produced and analyzed.
were assumed to follow the linear rule of mixture and to 1. Morphology and phases of PVD alumina films de-
be not temperature dependend which causes only little pend on deposition parameters, mainly on deposition
rate, chamber pressure and atmosphere, film thickness,
and substrate temperature. Alpha alumina can be ob-
Bond TBC tained either by high substrate temperatures or by post
300 -Coat .-, Tsurface =1200PC
coat annealing.
\ Tbutton =90(rC
2. Morphology of alumina-zirconia mixtures is more
- 'X \ compact than that of pure YPSZ TBC but it depends to
\ \
a less extend on the amount of alumina. Zirconia is very
u 100 - L \\b
fine grained in an amorphous alumina matrix.
\ " x
\ 3. By defining an effective thermal resistivity coating
£ • \ \
\\
\
\
thickness can be adapted to constant heat flux conditions
.2 400 \ \ x With a graded transition zone and with a modified coa-
->d
-oet ting thickness a variation of the stress field can be
**■----*:
\ \c achieved.
P4 -200 a ^
I.I.I
0 0L2 0.4 0.6 a8 1 1.2 Optimization of coating conditions and the design of the
Distance to Substrate [mm] graded zone with respect to chemical composition pro-
file, phases, microstructure, morphology, and optimized
curve Zr02 FGM heat flux lifetime remains to be carried out in future research.
um urn MW/m2
a 200 0 2
b 0 200 4.2 ACKNOWLEDGEMENTS
c 150 200 2 The authors gratefully acknowledge technical support
d 0 920 2 by J. Brien, C, Kröder, H. Mangers and H. Schurmann..
The work was partially funded by the Deutsche For-
Fig. 13: Effect of corrected coating thickness on schungsgemeinschaft.
stress distribution
16-9
1 Kaysser et al. AGARD this issue 13 Ulion, N.E.; Ruckle, D.L.: Columnar grain ceramic
thermal barrier coatings on polished substrates. European
2 M. Peters, K. Fritscher, G. Staniek, WA. Kaysser, U. Patent EP0044329B1
Schulz: Design and properties of thermal barrier coat-
ings for advanced turbine engines. Materialwissenschaft 14 Alaya, M.; Grathwohl, G.; Musil, J.: A comparison of
und Werkstofftechnik 28(1997) 357-362 thermal cycling behaviour of graded and duplex Zr02-
thermal barrier coatings. Proc. FGM'94, eds. B.Ilschner et
3 J.T. DeMasi-Marcin and D.K. Gupta, "Protective al. Lausanne (1994) 405-411
Coatings in the Gas Turbine Engine", Surface and
Coatings Technology 68/69 (1994) 1 -9 15 Lee, W.Y.; Stinton, D.P.; Berndt, C.C.; Erdogan, F.;
Lee, Y-D.; Mutasim, Z.: Concept of functionally graded
4 S.M. Meier, D.M. Nissley, K.D. Sheffler, T.A. Cruse: materials for advanced thermal barrier coating applica-
"Thermal barrier coating life prediction model develop- tion. J. Am. Ceram. Soc. 97 (1996) 12, 3003-3012
ment", ASME J. Eng. Gas Turbine Power 114(1992) 258-
263 16 Vincenzini, P.: Zirconia thermal barrier coatings for
engine applications. Industrial Ceramics 10(1990)3,113-
5 U. Schulz, K. Fritscher, C. Leyens, M. Peters, WA. 126
Kaysser: Thermocyclic Behavior of Differently Stabi-
lized and Structured EB-PVD Thermal Barrier Coat- 17 Schmitt-Thomas, Kh.G; Dietl, U.: Surface and coat-
ings. Material Wissenschaft und Werkstofftechnik ings technology, 68/69 (1994) 113-115
28(1997) 370-376
18 Sun, J.H., Chang, E.; Wu, B.C.; Tsai, C.H.: The prop-
6 Schulz, U.; Fritscher, K.; Rätzer-Scheibe, H.-J.; Peters, erties and performance of (Zr02-8wt%Y203)/(chemically
M.; Kaysser, WA.: Thermocyclic Behaviour of Micro- vapour-deposited Al203)/(Ni-22wt%Cr-10wt% Al-
structurally Modified EB-PVD Thermal Barrier Coatings. lwt%Y) thermal barrier coatings. Surface and Coatings
4th International Symposium on High Temperature Corro- Technology 58(1993) 93-99
sion and Protection of Materials. Les Embiez 20-24.5.1996
19 Strangman, T.E; Solfest, P.A.: Ceramic thermal bar-
7 Sohn, Y.H.; Biederman, R.R.; Sisson Jr, R.D.: Micro- rier coating with alumina interlay er. US Patent 5 015
structural development in physical vapour-deposited 502(14.5.1991)
partially stabilized zirconia thermal barrier coatings. Thin
Solid Films 250(1994)1-7 20 T. Krell, U. Schulz, M. Peters, WA. Kaysser: Influ-
ence of various process parameters on morphology and
8 Maricocchi, A.; Bartz, A.; Wortman, DJ.: PVD TBC phase content of EB-PVD thermal barrier coatings. Poc.
experience on GE Aircraft engines. Proc. TBC workshop EUROMAT 97, eds. L.A.J.L. Sarton and H.B. Zeedijk,
Cleveland 1995, ed. W.J. Brindley, NASA Conf. Publi- Netherlands Soc. for Materials Science, Vol. 3(1997)
cation 3312 (1995) 79-89 3/29-3/32
9 Schmücker, M.; Fritscher, K; Schulz, U.: Haftmecha- 21 U. Leushake, U. Schulz , T. Krell, M. Peters, WA.
nismen in ausgewählten EB-PVD Wärmedämmschicht- Kaysser: A1203 - Zr02 Graded Thermal Barrier Coatings
systemen. Fortschrittsberichte der DKG, Werkstoffe, by EB-PVD - Concept, Microstructure and Phase Sta-
Verfahren, Anwendung 10(1995)4, 379-384 bility. Proc. 4th Int. Symp. on FGM'96.
10 Fritscher, K; Schmücker, M.; Leyens, C; Schulz, U: 22 U. Leushake, T. Krell, U. Schulz, M. Peters, WA.
TEM Investigation on the Adhesion of EB-PVD Thermal Kaysser: Microstructure and Phase Stability of EB-PVD
Barrier Coatings. 4th International Symposium on High Alumina and Alumina/Zirconia for Thermal Barrier
Temperature Corrosion and Protection of Materials. Les Coating Application. ICMCTF '97, Surface and Coat-
Embiez 20-24.5.1996 ings Technology (in press)
11 Goward, G.W.; Grey, DA.; Krutenat, R.C.: US Pat- 23 Fritscher, K; Schulz, U: Burner-rig performance of
ent No. 4 248 940 density-graded EB-PVD processed thermal barrier coat-
ings, in "Ceramic Coatings" ed. K.Kokini, New York,
12 Jamarani, F.; Korotkin, M.; Lang, R.V.; Ouelette, ASME MD-Vol.44 (1993) 1-8.
M.F.; Yan, K.L.; Bertram, R.W.; Parameswaram, V.R.:
Compositionally graded thermal barrier coatings for high
16-10
Introduction 2. Methodology
It is universally acknowledged that bond coats for This study used NIKE 2-D stress models [9] to describe
plasma sprayed TBCs have a strong effect on the a burner rig test of a 2.54 cm (1.0 in.) diameter solid
thermal fatigue life of the ceramic layer, but the Waspaloy rod coated with 0.013 cm (0.005 in.) of
specific mechanisms by which the bond coat influences NiCrAlY bond coat, which was in turn coated with a
TBC life have not been established [1, 2, 3, 4]. As an "top coat" of 0.025 cm (0.010 in.) Zr02-8wt.% Y203
example, bond coat oxidation has been clearly linked to (Figure 1(a)). It was assumed that the coating was stress
spallation of the ceramic top coat but the mechanism by free after deposition and cooling to room temperature a
which oxidation causes failure of the top coat has not reasonable assumption in view of the low in-plane
been clearly identified [1, 4, 5]. Similarly, bond coat compressive stresses (approximately 3 MPa) measured
coefficient of thermal expansion (CTE), roughness and on TBCs fabricated with appropriate temperature
creep properties have been proposed to influence TBC control [10]. The model geometry was a thin
failure, but the magnitude of the effects and the axisymmetric slice of the substrate coating system, as
physical mechanisms by which these effects act, and was used previously (Figure 1(a)) [8], with a condition
interact, have not been established [1, 2, 6]. of generalized plane strain. Unlike the Petrus model
[8], the current model used refined elements for the
Part of the difficulty in determining which mechanisms bond coat and top coat at the interface region to allow
and interactions are involved is due to the difficulty in for growth and analysis of an oxide layer. The
geometry of the bond coat/ceramic interface was no oxidation and a high equal to the oxidation rate
modeled as a sine wave to simulate circumferential found empirically for the alumina former NiCoCrAlY
surface roughness conditions typical for a plasma [11]. The oxide scale was assumed to be of uniform
sprayed TBC (Figure 1(b).) [6, 7, 8].
STEADY STATE
HEATING HOLD COOLING
BOND 1200
INTERFACE
900
m
DIMENSIONS E
comments on the effects of roughness and CTE, the The effects of roughness and bond coat CTE will first
remainder of the paper will concentrate on the effects be discussed to demonstrate general effects that are
of oxidation and creep of the constituents. important to the detailed analysis of the effects of
oxidation and creep that are presented later.
The output of the model included stress as a function of
time, location and direction. It was found that the radial Roughness
stresses in the ceramic (perpendicular to the interface)
provided an adequate indication of the stresses most Interfacial roughness has been the subject of study for
likely to promote delamination cracking. Therefore, TBCs [6, 8] and is a well known consideration for
radial stresses will be presented as a function of time stresses in thin oxide films [14]. The case that is usually
and location to depict the active stress generating treated is for cooling of a two layer system that is
mechanisms in a TBC system and to try to define the assumed stress free at some high temperature, generally
likely failure scenario. It was also found (see Figures 8 the maximum temperature of a thermal cycle. Petrus et
and 9, discussed later) that the highest tensile radial al. [8] started with the assumption of a stress free
stresses were located at the bond coat peak, off-peak condition at room temperature and then modeled one
and valley regions. Therefore, the large volume of data entire thermal cycle for a plasma sprayed TBC. They
generated during FEM modeling will be narrowed to demonstrated that the assumption of a stress free
these locations in order to facilitate comparison of the
various cases presented. RADIAL STRESS CONTOURS (MPa)
"*"11 ^^v-i'o
3. Results and Discussion
condition at high temperature is valid for a TBC if the bond coat for these conditions. This large effect
bond coat and / or top coat relax sufficiently during a emphasizes the importance of roughness to the stresses
high temperature hold. When this holds true, Petrus et in the near interface region.
al. achieve qualitatively similar results as previous
studies: a rough interface causes an increase in While it is tempting to suggest that a high CTE bond
delamination stresses (tensile stresses normal to the coat favors cracking at bond coat peaks and should
interface) at the peaks of the metallic asperities while therefore shorten TBC life, a high bond coat CTE also
increasing compression normal to the surface in the favors crack arrest over the valleys at this small number
valleys between asperities, as compared to a smooth of cycles. As will be shown later, oxidation at larger
surface. This finding was reexamined for the current numbers of cycles may have a dramatic effect on this
study using four thermal cycles, instead of one, and situation. The combined effect of CTE and higher
including oxidation. It was found that, even when oxidation at large numbers of cycles has not yet been
oxidation is considered, a smooth surface had stresses investigated, but it is known that CTE interacts with
normal to the interface ranging from 5 MPa at the some of the other parameters studied [15].
interface to 0.5 MPa at the surface. A rough surface,
however, exhibited tensile delamination stresses as high The remainder of the discussion will focus on the
as 67 MPa at bond coat peaks, and normal compressive generation of stresses due to the operation of oxidation
stresses as high as 17 MPa in the valleys (Figure 3(a)). and creep. The effects of the various mechanisms are
Tensile stresses were also developed in the valleys due best understood when presented separately, starting
to oxidation, as will be discussed later. Several points with the simple, all elastic case.
of interest arise from this finding: 1) roughness tends to
cause a substantial increase in stresses in the near Elastic Case
interface region as compared to a smooth interface, and
2) the stress is periodic and should tend to limit In the case where all the constituents are elastic and
cracking to the bond coat peak regions. This result is there is no oxidation, the radial stress profiles for the
qualitatively the same as for previous models that three locations are the same for all four cycles (Figure 4
considered only cooling and only one thermal cycle. A (a)). At the bond coat peak, the ceramic is driven into
third point is noted that the properties of the bond coat compression on heating and returns to zero stress on
have nearly no influence on TBC stresses for the case cooling. The stress in the off-peak location also cycles
of a smooth interface. However, bond coat properties between a compressive stress and zero stress. However,
become extremely important when the interface is the stress in the valley cycles between a zero cool-down
rough [7, 8]. Note that roughness profiles, and hence stress and a tensile stress of approximately 70 MPa
the effect of roughness on stress, are extremely during steady state heating. The change in sign of stress
complex and at an early stage of investigation. from the peak to the valley is a result of the change in
sign of the radius of curvature [14].
Bond Coat CTE
Oxidation is introduced incrementally in the model,
For the case of a smooth interface, it has been with the first increment of oxidation starting with the
demonstrated that the properties of the bond coat are start of the hold time in the third cycle. Oxidation
relatively unimportant. Rather, the substrate dominates changes the stresses significantly from the case of no-
stress generation in the coating. The strength of an oxidation (compare cycles 1 and 2 to cycles 3 and 4 in
effect in the presence of a rough bond coat is Figure 4(b)). Oxidation pushes the peak radial stress
exemplified by the effect of bond coat CTE. Figure 3 toward slightly greater compression for all portions of
shows the radial ceramic layer stresses for low CTE the cycle and reduces the amplitude of the stress per
(Figure 3(a)) and high CTE (Figure 3(b)) bond coats cycle. Concurrently, the off peak cool-down stress
and the standard conditions of oxidation and creep of becomes tensile while the off-peak steady state stress
the top coat and bond coat. The tensile stresses at the remains constant. This results in an increase in the
bond coat peak increase from 67 to 240 MPa with the cyclic stress amplitude for the off-peak location, where
increase in CTE from 14 x 10"6 °C_1 to 20 x 10"6 °C1. an increase in the stress amplitude may have an effect
Also, the valley stresses for the high CTE case are on the fatigue life of this material [16]. The most
compressive, up to -221 MPa, while the valley stresses dramatic effect is on the stress in the valley. The valley
for the low CTE case tend to be mildly tensile. The stress increases rapidly with growth of the oxide. The
factor of approximately three change in stress level at linear, nearly vertical increase in stress reflects the
the peak scales well to the approximately factor of three linear growth rate imposed on the oxide. This rapid
change in CTE mismatch between the ceramic and the increase in stress is exaggerated due to the fact that all
17-5
200
CYCLE 1
HEATINQ
have a significant effect on stress when no other
_ 150
HOLD mechanisms are active. Interestingly, the four cycle
/ COOLING
stress history in Figure 4 b would result in a system that
S 100 -+- could not support the commonly observed peak-to-peak
ÜJ 50
V cracking. This is due to the fact that the radial stresses
w ffB
DC
^ are zero or compressive over the bond coat peaks for
0
-50
-OFF-PEAK all 4 cycles. A real system, however, includes plasticity
-PEAK and creep (or, identically, stress relaxation) behavior
-100
that can modify these stresses substantially.
-150
500 1000 1500 2000 2500
TIME (seconds) Creep Effects
(a)
Before examining the interactions of top coat and bond
coat creep with oxidation, it is informative to revisit the
effect of these mechanisms on stresses in the absence of
oxidation [8] and to examine how these stresses evolve
with thermal cycling. Figures 5 and 6 show the peak,
(b)
Figure 4. Ceramic layer radial stresses at the peak, off-
peak and valley locations for the case of (a)
no oxidation and (b) oxidation starting at cycle
3.
of the oxide growth for the first three cycles is made to 1,000 1,500 2,000 2,500
occur during a short time in the third cycle in the TIME (seconds)
model. As is shown later, this exaggerated increase in Figure 5. Variation in ceramic radial stress with
stress is unimportant when creep is allowed to occur. position and time for the case of a creep in
The valley cool-down stress increased from zero to 120 the ceramic and an elastic bond coat.
MPa and the steady state stress increased from roughly
50 to 170 MPa, causing a substantial oxidation induced off-peak and valley stresses for systems with creeping
increase in the mean tensile stress. Figure 4(b) provides top coat-elastic bond coat and creeping bond coat-
an indication of the stress levels that are expected if elastic top coat, respectively. While the stress always
stresses due to oxide growth are not relieved. The returns to zero in the all-elastic case (without
change in stress amplitude for all locations is a result of oxidation), creep of the top coat (Figure 5) causes a
the change in the system properties with the addition of deviation from a zero cool-down stress at all locations.
a very thin oxide layer with a drastically different At the peak, the cool down stress progressively
modulus and CTE than the coating materials. increases from the zero stress of the all-elastic case to
roughly 70 MPa during the course of four thermal
This simple case illustrates a trend found throughout cycles. This gradually increasing cool-down stress is a
the results: oxidation tends to cause radial cool down result of partial relaxation of the top coat compressive
stresses at the bond coat peak to become more negative, stresses at the peak during each steady state hold
and off-peak and valley stresses to become more (Figure 5). As with any relaxation process, the
positive. Oxidation also causes a decrease in the peak relaxation rate decreases with time (apparent as a
and valley cyclic stress amplitudes and a concurrent decreasing slope during the hold time) and the cool-
increase in the off-peak stress amplitude. down stress tends toward a steady state value.
However, relaxation is not complete by the end of the
The overall view in Figure 4(b) is that oxidation can fourth cycle. The stresses at the off-peak and valley
17-6
\ r \U1r.. \U-1r.
200 coo .ING
1—
OFF-PEAK 250 E NO-OXIDATION
1 150
■ s~- U.. * s
200
150
■ OXIDATION
iJI-rii
0) +" • PEAK
U ■J
«2 loo EC £• 100
i U «
= (O
i
afft
J 3 w 50
50 --- 1 r 1 < oc
III </)
0
u
V ^VLA
^V ILA
r\
-50
0 /\_A -100
-50 V
,1
ALL
ELASTIC
TOP
COAT
BOND
COAT
ALL
CREEP
CREEP CREEP
-100 ' ■ i i
500 1,000 1,500 2,000 2,500 (b)
TIME (seconds)
g£
Bond coat creep has a similar effect to that of top coat 2*
o
creep, but the rate of the change is much higher simply <aDC
because the bond coat creeps more readily (Figure 6)
[2]. The magnitude of the change is also greater within TOP BOND ALL
the first four cycles. This is due to the fact that COAT
CREEP
COAT
CREEP
CREEP
be expected of either mechanism separately. This result Figure 3(a)), shows that the area of the peak, off-peak
illustrates the difficulty in looking at only one and especially the valley tensile stress regions are small
mechanism to characterize the stresses in a system that and not connected. Also, a compressive region exists
has more than one mechanism operating. The various between the peaks (note Region 3 in Figure 8) that
creep cases do, however, remain consistent in the sign would prevent planar cracking from peak top to peak
of the stress at each location. top, as is usually observed in failed coatings. Similarly,
the high tensile stress regions are not in the same plane,
Comparing oxidized to non-oxidized stresses shows so that cracks moving from one region to another
there is interaction among mechanisms at the valley would have to follow a path well off of parallel to the
location but that oxidation and creep are roughly interface. Since cracks typically appear planar and
additive for the other locations. In the valley, oxidation parallel to the interface, this configuration does not
reverses the stress from moderately negative in all three agree with observations.
creep cases to highly positive. An interesting point is
that the magnitude of the stress change is for the RADIAL STRESS CONTOURS (MPa)
creeping bond coat and combined creep case is nearly
double that for the ceramic creep case.
67 -10
"X
~n'
REGION 2 v.
REGION 3 \>
67s. y//
1.1 !
REGION 2
h
10 um
122 -~-~-Srr 39
-\— L!
H 10 jim
Figure 8. Contour plot of the radial stress in the h
ceramic after cool-down following 4 thermal (b)
cycles for the case of combined oxidation Figure 9. Contour plots of the radial ceramic stress at
and creep. cool-down following 51 cycles for the case of
combined creep and (a) without oxidation
the interface. Examination of a contour plot for the and (b) with oxidation.
general creep case with oxidation, Figure 8 (same as
17-8
This difficulty was resolved by running an additional form a uniform tensile region a short distance from the
combined creep case, with and without oxidation, to a bond coat. Cracking in this large, general tensile region
greater number of thermal cycles. Figures 9 compares then becomes possible. Thus, it appears that early
contour plots of the radial stresses after 51 thermal cracking should occur at bond coat peak regions, aided
cycles for non-oxidizing and oxidizing cases. For 51 by ceramic and bond coat creep, and continued
cycles without oxidation, a compressive stress region cracking is driven by the changes in stress over the
extends from the valley well into the ceramic layer valleys driven by oxidation.
(Region 2 in Figure 9 (a)). These compressive stresses
would limit cracking to the region above the bond coat Summary and Conclusions
peaks. With oxidation, it is evident in both the 4 and
51-cycle cases (Figures 8 and 9(b), respectively) that The failure mechanisms of thermal barrier coatings
the tensile stress at the peak has become compressive. have been examined through a FEM model of residual
In the 51 cycle case the tensile off-peak region has stress generation due to oxidation, top coat creep and
grown in area and become more tensile. Additionally, bond coat creep. The results indicate that top coat and
the compressive radial stress regions over the valley bond coat creep act to generate tensile stresses at bond
after 4 cycles with oxidation (Region 3 in Figure 8), coat peak and off-peak locations, while generating
and after 51 cycles without oxidation (Region 2 in compressive stresses in the valley regions. The primary
Figure 9(a)), have completely disappeared after 51 effect of creep on the cool down stress appears to be
cycles with oxidation (Figure 9b, Region 2). Thus, the limited to short times at temperature. Oxidation acts
stresses at 51 cycles with oxidation is much more opposite to creep, causing increasing tensile stresses in
favorable to crack growth over the valleys than is the the valleys while pushing the peak regions into
four cycle case. However, the stresses at 51 cycles are compression. Continuing oxidation results in an
unfavorable for crack growth at the peak of the bond increase in the size of the tensile region associated with
coat. The histograms in Figure 10 provide a direct the valley, growing to the point that the tensile stresses
comparison of the stresses at each position for each in the valley extend beyond the peak regions.
case. Continued oxidation also causes a slight increase in the
maximum tensile stresses at the off-peak and valley
locations.
References
Carlos A. Arana
Aero Propulsion & Power Directorate
Air Force Wright Laboratory
Wright Patterson AFB OH USA
liner system, TBCs act as thermal insulators reduc- • Combustor Inlet Temperature: As the inlet
ing the amount of heat transmitted to the combustor temperature to the combustor increases, gas tem-
wall from the products of combustion. The heat perature also increases and, therefore, the liner
flux gets dissipated through radiation back into the temperature increases proportionally requiring
hot gas path, conduction across the coating and more air to maintain a given liner temperature.
substrate, and convection through back side cooling TBCs can ameliorate this issue by enabling a re-
air. duction in the amount of radiant heat flux from the
flame to the liner (figure 2).
2.0 THE PROBLEM
Fuel Injector
O)
3
O
V
\ x^ E
_o
IE
o
3
Figure 1. Internal view of a typical annular m
0
combustor system 0.01 0.03 0.05 0.07
Fuel-Air Ratio
Two major factors will significantly impact the Figure 3. Increased combustor exit
ability to cool the combustor liner: temperatures reduce available
dilution and cooling air
18-3
From the practical combustor design standpoint, • 99% efficiency at all operating conditions
the objective is not to use improved cooling meth-
ods to lower wall temperatures below the levels Other combustor system design considerations:
achieved with present methods, but rather to
achieve the same wall temperatures despite reduc- • 2000 engine flight hours and 4000 Low Cycle
tions in the amount of coolant flow as shown in Fatigue TACS - hot parts
figure 4. • 4000 engine flight hours and 8000 Low Cycle
Fatigue TACS - cold parts
Impingement/Film
CMC • 10 hour of operation at any steady-state flight
Uncooled CMC o
Hybrid Liner condition
+900
mat 5 • 1 hour of operation at simultaneous T3/T4
■f
Q.
a Quasi-Transpiration
Ü
ODS <* condition
9>
Single Cristal Effusion Cooled
+600 ^- Impingement/Film CMC go« c
• 350 hour of test cell operation
0)
a. '^ou
B 3
■o
a>
Other requirements:
Quasi-Transpiration a
•a
a +300 Double Wall 10' ■a
m £
a
a Film-Cooled 3
• Buckling margin of safety as dictated by com-
Effusion Cooled o-
Double Wall
Single Wall 0)
DC
bustor geometry (shape and wall thickness) shall
exceed 50 percent under maximum pressure and
TIME temperature
Figure 4. Combustor liner material/cooling
progression for the past 15 years. • 5% creep limit without loss of mechanical
alignment, as dictated by structural geometry to
4.0 GENERAL REQUIREMENTS AND DE- mitigate the interaction of thermal and mechanical
SIGN CONSIDERATIONS stresses
The primary technical objective of combustor re- • Component high cycle fatigue life greater than
search and development is to provide the technol- 107 cycles, as determined by both material and
ogy that will enable the designers attain combustor structure
performance, operability, weight and cost objec-
tives for a particular mission and cycle definition. • -3a material properties (average for crack
propagation)
As systems become more efficient and probably
more complex, the next generation of combustor • 500 0-maximum-0 cycles at maximum power
designs will need to balance conflicting demands
such as increased temperature while improving du- 5.0 COMBUSTOR TBC EVOLUTION
rability, maintainability, and affordability. In order
to achieve these demands an integrated prod- Since the birth of thermal barrier coatings in the
uct/process must be implemented. This process early 1960's, combustors and augmentors have
should enable the production of a robust design benefited from a series of incremental growths in
capable of operating efficiently in a wide variety of this technology. This has resulted in an increased
environments dictated by the intended air vehicle. ability to withstand higher temperatures with a de-
To satisfy the requirements of future cycle parame- crease in the severity of oxidation damage, and an
ters, technology improvements are required. Ad- improved cracking resistance to severe thermal
vanced planning allow us to establish combustor cycling for longer periods of time (Harris,
performance objectives that will support overall 1991) [3]. Early combustor coatings consisted of an
engine system improvements. Some combustor outer ceramic insulative layer of air plasma depos-
specific performance objectives are: ited, 22 weight percent MgO stabilized Zr02, and
an inner metallic bond coat of air plasma deposited
• 50% weight reduction Ni-Cr or Ni-Al (Gupta, et al.,1994)[4].
• 100% increase in volumetric heat release rate
• +900° F maximum combustor exit temperature As these temperatures in the combustor have been
• +400° F combustor inlet temperature increased, the durability of this early combustor
• 300% increase in fuel/air turndown ratio coating has become inadequate. The maximum use
18-4
temperature of the plasma sprayed magnesia- duced through the liner walls downstream of the
stabilized zirconia coating is on the order of dilution jets may not be effective in the combustion
1255°K since magnesia-stabilized zirconia crystal- process, therefore, contributing to potential com-
lographically destabilizes above 1227°K. bustor inefficiency.
£ 4
D)
Electron Beam Physical
Vapor Deposition £
O 3
J
Generation 3 ^ Nano-layers z
for Reduced Thermal O) 2|—
^ Conductivity
Flame Spray-Air Plasma Deposited Figure 6. Efficiency Benefits for Idle Engine
for Improved Component Durability
Operating Condition
Figure 5. TBC Progression for Hot Section
Applications 6.2 Exhaust Smoke Emissions
way between the maximum gas temperature and the was performed in order to estimate the thermal per-
temperature of the cooling flow (T3). Also shown formance of TBCs and compare the capability of
in this figure is the relationship of the allowable TBC improvements against cooling effectiveness
design metal temperature levels to the incipient improvements and T3 reductions. The heat transfer
metal temperatures. As shown, the allowable metal model outlined in this section closely follows the
temperature levels are near the incipient melt point development approach as detailed by Lefevbre, et
of nickel based superalloys. Therefore, any break- al., 1960 [11].
down in the cooling layer would be detrimental.
Experimental Data
40 r —Uncoated Panel
6. 360 —Coated Panel
E Calculated Data
v. 32 CD
<D *-> 320 o Uncoated Panel
n
■ Coated Panel
i_
E co
3
■i—■ ?80
Z 24 ID
JS c
o
E Test Condition Fuel 1 240
0) 16
LU / A Cruise Jet A
Q.
< O Takeoff Jet A E 200
CO / CD
•4—>
_L j 1 1 1 CD
+■ 1100 1180 1260 1340 1420
>
co
120
Stabilized Zirconia
Ceramic Coating
Ceramic Filled Bond
Coating Surface
w E Xl
Hot Combustion °° • Metal
Substrate
Gas Radiation Wie
lel
and Convection
1
Ceramic, Tmetal
Tg Surface
TBC
Figure 10. Schematic of radiant heat flow through the TBC at isothermal conditions.
7.1 Internal Radiant Heat Flux 7.2 External Radiant Heat Flux
The radiation heat flux was calculated for the case Radiant heat flow from the combustor liner to the
of a TBC-coated liner using the results obtained by casing is given by:
Liebert, 1978 [12] where emittance of both the ce- £
ramic translucent layer and the material substrate Wec[TWC~T3 )
<1 RAD = ° A {5}
were taken into consideration.
ec+ew(l-ec){ yAc)
<1"RAD = <y(0.5)egT*-5(Tg2-5 -TWH25) {1}
7.3 Internal Convective Heat Flux
The radiation heat flux was calculated for the case Convective heat flow to the combustor liner from
of an uncoated liner using reference [11]: the gases is given by:
Eg=\- exp(-290P LV(/ / a)lb Tfh5) {3} Convective heat flow from the combustor liner to
the annulus air is given by:
* 800,000
* r%in
where T) is the thermal effectiveness, which is a i 600,000 1 * •
R-in
function of the internal heat transfer coefficient I 400,000
within the liner. This heat flux is absorbed by the C-out
200,000 R-ouk *
coolant and therefore reduces the heat load that 'K
0 ''— F —
must be removed by q'coNv and q"RAD on the out-
0 0.2 0.4 0.6 0.8 1
side liner. % TBC thickness (3.81x10-4)m
The following figures summarize the results of the Figure 12. The Effect of TBC Thickness on
investigation: heat fluxes
> .N \ 1S1400
1500 <u
2
3
as 1400
~—-j
!i=»*_.
-i
vT- rneiai nut siue £1350
<u
T- metal cold sic le_
03 x\ =1300
xX
Q.
E to
1300
x> Sv ^s,
51250
5 1200
\^ 1200
700 750 800 850 900 950 1000
1100 1 1 ■ Coolant Air Temperature ~ K
0 0.02 0.1 0.2 0.5 0.75 1
% TBC thickness (3.81 x 10-4)m
1600
tem. TBC reliability must be equal or greater than
the substrate since the integrity of the coating be-
T-TBC
* 1400 comes prime reliable as loss of the coating can re-
§ sult in rapid deterioration of the liner. Also, the
I 1200 ability to apply TBCs to components with complex
(D
geometries should be improved. As combustor in-
11000 — i -meiai noi siae
L~-~- ._-••■•
let air temperature approaches the maximum al-
I 800 —
C5— — lowable wall temperature, there simply will not be
T-metal cold side
enough cooling potential left to cool the wall ade-
600 quately. This leaves the possibility of no metal wall
700 750 800 850 900 950 1000 structure, no matter how effective, will maintain
Coolant Air Temperature - K
reasonable metal temperatures at those conditions.
Figure 14. The Effect of Coolant Air An analytical investigation was conducted to de-
Temperature Variation on TBC termine, through qualitative trends rather than
Liner Wall Temperature quantitative values, the effect of a ceramic coating
on the average metal temperatures of a full annular,
• Figure 15 shows the impact of liner thermal semi-transpiration cooled combustor liner. This
effectiveness variation on wall temperatures. As investigation was conducted at a pressure of 3040
would be expected, the hot side metal temperature KPa, inlet air temperatures of 1144°K and an
is reduced as the ability of the liner to extract heat overall fuel air ratio of 0.67. The insulating effects
becomes greater for the same amount of coolant as well as increased reflectivity of the ceramic
flow rate coating were responsible for the reduction in heat
transfer through the liner walls.
1700
^ 1600 T-tbc
To solve the heat flux equations a non linear algo-
rithm was used. For that purpose a Newton-
2 1500
Raphson non-linear solver was developed to solve
JS
Q-
1400 the equations in an iterative manner. It should be
1 1300 T-me tal he t side
mentioned that the energy balance equations repre-
h, „„.. ^*™~i h-«^z; sent a gross approximation to the heat transfer
=5 1200 process taking place between the liner structure and
T ■meta cold side- —=-—,
1100 1i
i 1
its surroundings. A more rigorous formulation for
0.5 0.55 0.6 0.65 0.7 the heat transfer process may be developed if spe-
Thermal Effectiveness cific details are required.
2. Turner D., "Enhanced Performance Extended 13. Odgers, J., Kretschmer, D., "A Simple Method
Life Combustor", WL-TR-2136, November for the Prediction of Wall Temperatures in a
1996 Gas Turbine", ASME Paper No.78-GT-90,
April 1978.
3. Harris, D., "Composite Thermal Barrier
Coatings", Fourth National Spray
Conference, CF, May 1991.
1350
1300
1250 Thermal Barrier Coatings
1200
o 1150
£3 1100
Directionally Solidified
2at 1050
Alloys
ingle Crystal
a 1000 Alloys
Ea> 950
900 Conventionally Cast
Alloys
850
800 Wrought Alloys
750
1940 1950 1960 1970 1980 1990 2000 2010
Year
Figure 1 Increase in operational temperature of turbine components made possible by alloy, manufacturing technology and thermal
barrier development.
Hot Gas
Temperature
Temperature Profile
without TBC
1400 1600
Temperature (*C)
Temperature Profile
with TBC
Coolant Figure 4 Relative phase compositions of aged 8% PYSZ TBC's.
Temperature Note the initial coatings were all T face-centred tetragonal
phase.
* R
■2 ^
E Cooling i
Isothermal
ageing
5 Heating
I2 --
\-
n
200 400 600 800 1000 1200 1400
Temperature (°C) Figure 5 An example of a combustion can with plasma sprayed
TBC on the inner can walls.
Figure 3 Thermal diffusivity of CaSZ TBC showing thermal
instability due to precipitation of CaO out of solution.
performance are discussed in the following sections.
20-3
200 400 600 800 The relatively high modulus associated with the plasma
Temperature (°C)
sprayed TBC can lead to high thermal strains in the
coating, especially around tight radii of curvature which
Figure 11 Young's moduli of solid sintered 8% PYSZ, plasma
sprayed and EB-PVD 8% PYSZ TBC systems measured by
are found around the leading edges of high-pressure
dynamic techniques. turbine aerofoils. As an example Figure 12 compares the
microstructures of a plasma and an EB-PVD TBC which
have undergone thermal cyclic testing. In the case of the
For TBC applications in large turbine components plasma sprayed coating, there has been significant
(particularly in the power generation industry) the cracking within the ceramic leading to gross failure of the
relative roughness of the plasma sprayed TBC compared TBC system, whereas the EB-PVD system has suffered
to the blade chord dimensions becomes more favourable little mechanical damage within the ceramic layer.
and the effects of increase friction coefficient become
less important. Similarly, the local gas velocities in the
combustion zones are lower whcih can lead to low
coefficient of friction losses and use of plasma sprayed
TBC systems in such circumstances are not so
detrimental to performance.
Strain Compliance
0.16
Si02 ERODANT AT ROOM
TEMPERATURE '
„0.12
£ \ EB-PBD TBC
l
20 40 60 80
Impingement Angle (")
Plasma Sprayed EB-PVD Substrate Figure 14 Room temperature erosion rate of plasma and EB-PVD
TBC systems as a function of impingement angle.
Figure 13 Erosion resistance of plasma sprayed and EB-PVD
TBC's, compared to that of a metallic substrate.
Coating (Elastic)
EB-PVD MCrALY
v bondcoat
BONDCOAT DEVELOPMENT
Table of contents
one coating manufacturer to the other. So, great care 1.2.4 Plasma assisted CVD
should be taken when comparisons are made! A. Lasalmonie (SNECMA) - This new deposition
There is place for both processes (PS and EBPVD) in technique seems to present good potentials ; its
aeroengines and other parameters than durability application to coat real parts has not still to be
have to be taken into account. demonstrated and the costs have to be evaluated.
XXX (??) - Have other studies on reproducibility of PS V.R. Parameswaran (NRC) - The use of chlorides (ZrCl4)
coatings been performed (c/paper N°4)? may lead to environmental issues.
A. Lasalmonie (SNECMA) - Equipments are now XXX (DLR) - The ability to coat compex-shaped parts
availavable and installed, for a good control of the remains questionable.
reproducibility.
S. Drawin (ONERA) - The problems with chlorides can
C. Moreau (NRCC) - A good control for a better be solved even for industrial plants ; other zirconium
reproducibility can be achieved by controlling the precursors are currently studied. A larger deposition
particle velocity and temperature and the substrate reactor will be built next year, that will be able to
temperature. The latter depends of the type of thermal coat larger parts, as well coupons as blade and vanes.
contact between splats.
J.R. Nicholls (CRANFIELD Univ.) - The use of a plasma
XXX (??) - How much variability can be tolerated? will make it possible to deposit layered systems,
similarly to what has been achieved in our laboratory.
C. Moreau (NRCC) - For variations of 100°C for the
substrate temperature, different coating micro-
structures were observed. 1.3 Lifetime of TBCs
P. Morrell (ROLLS-ROYCE) - From the point of view of 1.3.1 Role of bond coat, TGO and oxidation
the manufacturer, the shape of the component to be
coated is of prior importance. The now available A. Lasalmonie (SNECMA) - Differences in the coating
processes are a compromise of what the engineers properties are observed when different bond coats are
want and what they can get. used : how can they be explained? The role of the bond
coat surface preparation (roughness) has to be assessed.
1.2.2 EBPVD Various failure modes are encountered in or at the TGO
(Thermal Grown Oxide): the failure mechanisms have to
A. Lasalmonie (SNECMA) - Different coating structures be understood more precisely.
and properties can be obtained by EBPVD. If we refer to
the presentation by G. Marijnissen (INTERTURBINE, M.R. Winstone (DERA) - Are the failure mechanisms the
paper N°10), this is a process that has to be further same for TBCs having experienced 1000 one-minute
controlled : the influence of a lot of parameters has to be cycles and 1000 one-hour cycles? TBC
studied. manufacturers have large databases on the possible
failure modes...
V.R. Parameswaran (NRC) - Grading should have a good
influence on the properties of the EBPVD coatings. W. Brindley (NASA) - Yes, TBCs cycled at maximum
temperature for six or sixty minutes present different
failure modes. It is thus important to design the right
1.2.3 Layered structures
durability tests : maximum temperature, coolant
A. Lasalmonie (SNECMA) - Interesting results have temperature, time at maximum temperature,...
been presented (cf. paper N°6), but no thermal stability
and lifetime values are available for these TBCs. For highly oxidized TBCs, failure occurs at locations
where the experienced temperatures were the highest.
Furthermore, are these structures industrially feasable?
J.R. Nicholls (CRANFIELD Univ.) - These are still too About the surface preparation, the presented paper
new systems : they have only been tested at 1100°C (N°17) adressed PS coatings with a simplified (sine-
wave) roughness. Can the conclusions be generalized
for 50 hours.
to all PS coatings?
A. Lasalmonie (SNECMA) - The layer thickness appears
to be a fundamental parameter. 1.3.2 Testing methodology
J.R. Nicholls (CRANFIELD Univ.) - Yes, the layers A. Lasalmonie (SNECMA) - As shown by paper N°15,
should have a thickness comprised between A/4 and X burner rig tests may be very difficult to interpret: there is
(k: wavelength of the incident radiation) to a need to establish a testing methodology for the
efficiently reduce the radiative contribution to the durability tests.
thermal conductivity. Note that the layers deposited
by K.S. Ravichandran (UTAH Univ., paper N°14) 1.3.3 Modelling
were either too thick (15-49 urn) or too thin (0.06 - A. Lasalmonie (SNECMA) - Thermomechanical
0.15 um) to induce a noticeable reduction of the modelling of the TBC system should allow a better
thermal conductivity. understanding of the failure mechanisms.
D-4
1.4.1 Erosion
A. Lasalmonie (SNECMA) - Can the erosion behaviour
of PS coatings be improved?
C. Moreau (NRCC) - Yes, by using for instance powders
with smaller particle size ... but several problems will
then occur (feeding of the powder into the plasma
torch, moisture sensitivity, cost,...).
W. Brindley (NASA) - Another way to solve this
problem would be to find an erosion-resistant
material to spray on the leading edge (and having a
lower thermal conductivity than the superalloy
substrate).
C. Moreau (NRCC) - Here too, feedback from the users
is needed.
P. Morrell (ROLLS-ROYCE) - Currently, TBCs behave
well on combustor cans and vane platforms.
W. Brindley (NASA) - The problem is also to get the
right informations from the TBC coated parts : this
needs periodic inspections, inspections means, and
inspection procedures.
P. Morrell (ROLLS-ROYCE) - And how will the
metallic substrates be inspected if there is a TBC on
it.
8. Author(s)/Editor(s) 9. Date
Multiple April 1998
10. Author's/Editor's Address 11. Pages
Multiple 188
12. Distribution Statement There are no restrictions on the distribution of this document.
Information about the availability of this and other RTO
unclassified publications is given on the back cover.
13. Keywords/Descriptors
Protective coatings Plasma spraying
Heat resistant materials Coating processes
Barrier coatings Vapor deposition
Combustion chambers Research projects
Turbine blades Electron beams
Ceramic coatings Thermal conductivity
Temperature gradients Bonded coatings
High temperature Reviews
Aircraft engines
14. Abstract
Thermal barrier coatings are an emerging technology which will enable either an increase in
turbine inlet temperatures, or a reduction in the duty temperature of the blades and as a result,
an extension of their life-cycle.
This Workshop presented the state of the art, gave a description of the existing or emerging
technologies, summarised our present knowledge of damage mechanisms, including
microstructural, mechanical and thermal aspects, and reported on the benefits and drawbacks of
the various families as perceived by engine manufacturers and users.
A final discussion attempted to identify future R&D requirements.
NATO -0- OTAN
7 RUE ANCELLE • 92200 NEUILLY-SUR-SEINE DIFFUSION DES PUBLICATIONS
FRANCE AGARD NON CLASSIFIEES
Telecopie 0(1)55.61.22.99 • Telex 610 176
L'AGARD detient un stock limite de certaines de ses publications recentes. Celles-ci pourront eventuellement gtre obtenus sous forme de
copie papier. Pour de plus amples renseignements concernant l'achat de ces ouvrages, adressez-vous ä l'AGARD par lettre ou par t61ecopie
ä l'adresse indiquee ci-dessus. Veuillez ne pas UUphoner.
Des exemplaires supplementaires peuvent parfois etre obtenus aupres des centres de diffusion nationaux indiques ci-dessous. Si vous
souhaitez recevoir toutes les publications de l'AGARD, ou simplement celles qui concernent certains Panels, vous pouvez demander d'etre
inclus sur la liste d'envoi de l'un de ces centres.
Les publications de l'AGARD sont en vente aupres des agences de vente indiquees ci-dessous, sous forme de photocopie ou de microfiche.
Certains originaux peuvent egalement etre obtenus aupres de CASI.
m
Imprime par le Groupe Communication Canada Inc.
(membre de la Corporation St-Joseph)
45, boul. Sacre-Cceur, Hull (Quebec), Canada K1A 0S7
jMBs^m
NATO -®- OTAN
7 RUE ANCELLE • 92200 NEUILLY-SUR-SEINE DISTRIBUTION OF UNCLASSIFIED
ISBN 92-836-1073-3