The effect of the TiC particle size on the preferred oxidation temperature for self-healing of oxide ceramic matrix materials

The effect of particle size on the oxidation kinetics of TiC powders is studied. Different sizes of TiC powder ranging from nanometre to submillimetre sizes are investigated. The samples are heated at different heating rates from room temperature up to 1200 °C in dry synthetic air. The Kissinger method for analysis of non-isothermal oxidation is used to estimate the activation energy for oxidation of the powders and to identify the active temperature window for efficient self-healing. The master curve plotting method is used to identify the model which best describes the oxidation of TiC powders, and the Senum and Yang method is used to approximate the value for the Arrhenius constant. The oxidation of TiC proceeds via the formation of oxycarbides, anatase and then finally the most stable form: rutile. The activation energy is found to be a strong function of the particle size for particle sizes between 50 nm and 11 µm and becomes constant at larger particle sizes. The data demonstrate how the minimal healing temperature for oxide ceramics containing TiC as healing particles can be tailored between 400 and 1000 °C by selecting the right average TiC particle size.


Introduction
Embedded titanium carbide (TiC) particles are being considered as a potential healing agent to autonomously repair crack damage in oxide ceramics used in high-temperature applications [1]. For such (extrinsic) self-healing ceramics, microscopic cracks in the material will also intersect the sacrificial healing particle and allow atmospheric oxygen to reach it.
The local supply of oxygen via the open cracks to the intersected TiC particles will cause them to oxidize, leading to the formation of TiO 2 which both fills the crack and adheres strongly to the faces of the crack [2]. Hence, the so-called healing reaction leads to a partial or complete recovery of the mechanical properties of the material and an extension of the lifetime of the component. To be a successful strategy, it is important that the oxidation reaction takes place at the prevailing conditions (in particular the right temperature) and with the right kinetics. A general analysis of the required properties of the healing particles leading to autonomous self-healing at high temperature can be found in [2].
However, the cited analysis did not specify the minimal temperature window required to trigger the healing reaction. Of course, the required minimal temperature to trigger the healing reaction depends primarily on the chemical stability of the healing material in an oxygen-containing environment. This dependence has already been demonstrated (but apparently not been appreciated) in the combined first-generation studies on self-healing of Al 2 O 3 or Si 3 N 4 filled with SiC or Ni particles [3][4][5][6][7]. In these studies, the crack healing ability was studied as a function of both temperature (from 900 to 1400°C) and time (between 1 and 300 h) and widely different values were obtained for different particle sizes.
While the optimal healing temperature has been shown to vary significantly with the chemical composition of the healing particle and to a lesser degree of that of the matrix material, very few studies, with the exception of [8], performed a systematic study on the effect of particle size on the healing kinetics given a fixed composition for both the particle and the matrix. In their study, Nakao and Abe [8] measured the strength recovered by high-temperature oxidation of alumina containing 18 vol% SiC particles of size 270, 30 and 10 nm. They reported that reducing the size of SiC reduced the oxidation/healing temperature from 1200 to 950°C, due to the faster oxidation kinetics at lower particle sizes.
Recently, the ability of 2-lm-sized TiC particles to heal surface cracks in alumina has been studied by Yoshioka et al. [1]. Healing of two different composites containing 15 and 30 vol% of TiC was studied as a function of temperature ranging from 400 to 800°C in air and a fixed time of 1 h. Full strength was obtained for a composite containing 30 vol% of TiC upon annealing in air at 800°C for 1 h. Also, in situ healing of these alumina composites was demonstrated [9]. Surface cracks were healed when the material was exposed in a combustion chamber with a high-velocity exhaust gas mixture at approx. 1000°C and a low oxygen partial pressure.
The oxidation behaviour of bulk TiC has already been investigated as a function of temperature (from 600 to 1200°C) and oxygen partial pressure (13-100 kPa) [10,11]. The observable trend in these studies is that oxidation proceeds parabolically from 600 to 800°C and switches to linear kinetics at higher temperatures and higher oxygen partial pressures. Four steps were identified during the transformation of TiC powders to TiO 2 [12,13]. The first step involves the substitution of atomic oxygen for the carbon present at the interstitial vacancies of the TiC lattice leading to the formation of the oxycarbides/ titanium suboxide layer. In the next step, the titanium suboxides oxidise into amorphous titanium dioxide and further crystallize to anatase. In third stage, the crystallization of anatase continues and the resulting volume expansion leads to the cracking of the oxide layers, and this provides fast short-circuit diffusion paths oxidation. In the fourth and final step, anatase is formed continuously and rather rapidly transformed into rutile. The rate-determining step is the diffusion of oxygen through the already formed oxide layers.
With TiC having been identified as a potential healing agent for alumina [2], the current work focusses on the determination of the effect of TiC particle size on its oxidation kinetics in order to explore the possibility to lower the minimal healing temperature of alumina containing TiC to temperatures below 1000°C. To this aim, differential thermal and thermogravimetric analyses are performed for TiC powders of different particle sizes over the temperature range from room temperature to 1200°C. The kinetic triplet, i.e. the activation energy, the reaction model and the Arrhenius constant, is determined for all powders and such information, in principle, allows the prediction of the optimal healing temperature for any inert ceramic containing TiC particles as the discrete, extrinsic healing agent.

Analysis of reactions
In solid-gas reactions, the solid phase is usually a granular medium through which the gaseous reactant circulates and reacts at the solid interface [14]. Due to the reactions occurring at a microscopic level, the solid material undergoes chemical and structural change as the reaction proceeds. Several reaction models are used in describing gas-solid reactions [15]. The unreacted shrinking core model is commonly used. In this model, the reaction product forms a solid layer that allows diffusion of reactant gas towards the interface between the product layer and unreacted core. A continuous reaction leads to the advancing of a reaction front, or equivalently, a shrinking core. The entire reaction kinetics involves three steps: transport and dissolution of the gas reactant at the solid surface, diffusion of the gas reactant through the product layer towards the unreacted core and finally chemical reaction of the gas with the solid reactant at the product-core interface [16].
Thermal analysis is employed to study the heterogeneous reactions. In (simultaneous) differential thermal analysis (DTA), the heat flow and, in thermogravimetric analysis (TGA), the mass change are measured while the material is heated at a constant the heating rate [17]. From the measured data, the kinetic triplet, which describes the conversion of a material (a) as a function of temperature (T) and time (t), is derived. The kinetic triplet parameters are the activation energy E A , the Arrhenius or frequency constant A and the reaction model f(a) [18]. For nonisothermal experiments, the rate of transformation of any material is given by: where the fraction converted at any time is given by: in which m 0 , m t and m ? denote the mass at time t 0, t and ? (i.e. after full conversion), respectively. dT/dt is the heating rate denoted as b, and da/dt is the isothermal conversion rate and expressed as: Combining Eqs. (1) and (3) leads to the nonisothermal conversion rate:

RT fðaÞ ð 4Þ
Integration results in the non-isothermal transformation rate law given by: This temperature integral has no analytical solution [17,19], but is usually approximated by a series expansion or converted in an approximate form that can be integrated. A frequently used approximation is that by Senum and Yang [20]. Consider x = E A / RT, then: With pðxÞ ¼ R 1 x e Àx x 2 dx, Eq. (6) is approximated with the fourth rational of Senum and Yang [20], i.e.
Successively, the activation energy, the reaction model and finally the Arrhenius constant will be determined from non-isothermal experimental data. In the above analysis, the kinetic parameters are considered to be independent of temperature and fraction converted.

Activation energy
Regardless of the reaction order, the activation energy of the conversion reaction E A can be obtained directly from non-isothermal TGA or DTA data using the so-called Kissinger method [21,22]. If the temperature rises during the reaction, the reaction rate will rise to a maximum value, but return to zero as the reactant is exhausted. The maximum rate is associated with either a maximum mass change or heat evolution, which occurs at the peak temperature. Then, the governing equation for the determination of the activation energy reads: where T p is the peak temperature, corresponding to the temperature at which the reaction rate is at maximum. Then, the highest heat flux or mass change rate is recorded. R is the gas constant, and C is a constant. Thus, the activation energy (E A ) is obtained by plotting ln (b/T p 2 ) versus 1/T P for a series of measurements at different heating rates b.

Reaction model
The reaction model mathematically describes the conversion process of a material as a result of an occurring chemical or physical reaction. Different models have been derived for solid-state reactions [18], which comprises nucleation, geometrical contraction, diffusion and reaction-order models. To determine the appropriate model to analyse the experimental data, the so-called master curve plotting method is used. This method is based on the concept of generalized time h [23], which is defined by Eq. (9) for isothermal or non-isothermal reactions, respectively, as: h denotes the reaction time taken to attain a particular fraction converted at infinite temperature. Now, Eqs. (3) and (5) can be rewritten in terms of generalized time, respectively, as: As a reference point for the master curve plotting analysis, generally the level a = 0.5 is chosen and using Eq. (10) it can be derived that: From the experimentally determined conversion rate, with respect to the generalized time, as a function of fraction converted, f(a)/f(0.5) versus a can be obtained, since By comparing the experimental result with theoretical reaction models [23], the reaction model that best fits the measured data can be determined. From this, the key reaction mechanisms can be determined and this knowledge can be used to optimise the material performance.

Arrhenius constant
The last of the kinetic triplet, i.e. the pre-exponential factor (A), can be evaluated from Eq. (11), once the reaction model is determined; cf. ''Reaction model'' section. The generalized time h for non-isothermal reactions [Eq. (9)] resembles Eq. (6) and hence can be written as: while p(x) is approximated using Eq. (7).

Experimental procedure
The five TiC powders used in this study are all commercially available, and their main features are summarized in Table 1. The first and second powders, having an average particle size of 140 and 48 lm, respectively, were obtained by sieving a TiC powder with a broad size distribution (Goodfellow, UK; product number TI546030). The third TiC powder with an average particle size of 5 lm was supplied by Alfa Aesar, Germany (product number TiC 40178). The fourth TiC powder with an average particle size of 0.8 lm was supplied by American Elements, USA (product number TI-C-03M-NP.800NS). Finally, the fifth TiC powder with an average particle size of 50 nm was supplied by Chempur, Germany (Product number CFF-NP059).
The purity of the powders studied was determined with X-ray fluorescence spectroscopy (XRF) using an Axios Max WD-XRF instrument (PANanalytical, The Netherlands), and the data were evaluated with SuperQ5.0i/Omnian software. The size distribution of the TiC particles was determined by laser diffraction using the Microtrac 3500 (Microtrac Inc., USA) equipped with a tri-laser diode detection system. The particle morphology of each TiC powder was recorded with scanning electron microscopy (SEM) using a JSM 6500F (JEOL, Japan). The instrument used was equipped with an energy-dispersive spectrometer (EDS; Ultra Dry 30 mm 2 detector, Thermo Fisher Scientific, USA) for X-ray microanalysis (XMA).
The oxidation kinetics of the TiC powders was studied non-isothermally using simultaneous thermogravimetric (TG) and differential thermal analysis (DTA) using a high-performance symmetrical double furnace balance (TAG 16/18, Setaram, France). 20 ± 1 mg of the powder was put into a 100-lL Al 2 O 3 crucible and placed onto a TG-DTA rod with an S-type thermocouple. A similar rod but without crucibles was placed as a counterweight on the other side of the balance to aid automatic correction of buoyancy effects. Then, the dual furnaces were heated up from room temperature to 1200°C with various heating rates (between 1 and 10°C min -1 ). In each measurement, both furnaces were purged with dry synthetic air (N 2 ? 20 vol% O 2 and H 2 O \ 10 ppm) at a total flow rate of 50 sccm. The evolved gas at the sample location was collected and analysed with an Omnistar GSD 301 system having a quadrupole mass spectrometer (QMS 200 M3 from Pfeiffer Omnistar, Germany) in order to observe the evolution of CO 2 . The thermal analysis system used was equipped with Calisto acquisition and processing software (v1.43 DB v1.44, Setaram, France) which records both the furnace and the sample temperature. For accurate results, the temperature measured by the thermocouples at the sample position is used for the analysis as this represents the true sample temperature. Prior to admitting the gases to the furnaces, each gas, i.e. Ar, N 2 and O 2 (all with 5 N vol% purity), was filtered to remove any residual hydrocarbons, moisture and oxygen (for N 2 and Ar). The Accosorb (\ 10 ppb hydrocarbons), Hydrosorb (\ 10 ppb H 2 O) and Oxysorb (\ 5 ppb O 2 ) filters (Messer Griesheim, Germany) were used. The flow of each gas was regulated and monitored with LabView (version 12, National Instruments, USA) using calibrated mass flow controllers (Bronkhorst, The Netherlands).
The oxidation products and lattice parameter were identified and measured by X-ray diffraction using a Bruker D8 Advance diffractometer (Bruker, Germany) equipped with a graphite monochromator. Diffractograms were recorded with Co Ka radiation in the 2h range of 10°-130°and a step size of 0.034°. The data were processed with the Diffrac EVA 4.1 Bruker software.

TiC powder oxidation
Oxidation of the 5 lm powder was studied first by heating from room temperature to 1200°C at 5°C min -1 in air, and the heat flow and the conversion curves are presented in Fig. 1. Transformation starts just before 400°C and ends around 900°C via two exothermic peaks at about 470 and 710°C. Therefore, to fully describe the oxidation of TiC, two separate experiments are performed on the same powder under the same conditions till the end of each peak (i.e. till 530 and 900°C) and the oxidation products analysed.
XRD confirms that part of the TiC powder is converted into anatase and rutile in the ratio of about 3:2 in the first peak and the remaining TiC, oxycarbide and anatase are converted into rutile in the second peak. The morphology of the 5 lm TiC powder at different stages of non-isothermal oxidation with a heating rate of 5°C min -1 in dry synthetic air is shown in Fig. 2. The starting TiC particles are irregularly shaped with sharp edges; see Fig. 2a. After oxidation up to 530°C, some anatase was formed at the surfaces; see Fig. 2b. Finally, after full transformation into rutile the particles became more rounded at the edges as shown in Fig. 2c.
As the experiment in Fig. 1 proceeded, the gases evolved were collected and analysed. CO 2 was the only gas detected, and the recorded pressure is Table 1 Size and composition of TiC powders with D 50 as the average particle size and D 10 and D 90 as the particle size at 10 and 90% cut-off, respectively, in the size distribution  presented in Fig. 3. It is noticeable that the removal of C via CO 2 release and the precipitation of TiO 2 occur simultaneously as the heat flow peaks in Fig. 1 and the CO 2 signal in Fig. 3 closely overlap at the same temperatures.
A zoom in at about 320°C of Fig. 3 shows a small peak (inset of Fig. 3) which is not immediately identified in the heat flow signal in Fig. 1. Hence, the mass change signal in Fig. 1 is further differentiated in Fig. 4, and then a kink coinciding with the peak at 320°C in Fig. 3 is identified.
Another experiment using a fresh sample was run under the same conditions until 360°C to determine the reaction occurring around 320°C. No oxide was identified by XRD after controlled cooling of the reaction product obtained although the sample weight had increased and CO 2 had been released. The combination of observations shows that the lowtemperature peak is due to the substitution of atomic oxygen at carbon sites of the TiC lattice [12], leading to the formation of TiC y O 1-y which is evidenced by the reduction in the lattice parameter of TiC [24]. In this specific case, the lattice parameter of the starting material reduced from 4.3275 ± 0.0001 to 4.3260 ± 0.00005 Å after heating up to 360°C. This oxycarbide phase exists throughout the oxidation   process [12,25] and thus may act as a precursor for the formation of anatase and rutile.
The oxidation of the 5 lm powder can be summarized as follows. A first oxycarbide formation stage takes place over the interval from about 300-360°C. A second oxidation stage which leads to the formation of both rutile and anatase sets in after the first reaction and peaks around 470°C and ends just before 530°C; see Fig. 1. (This stage is hereafter referred to as the anatase formation stage.) Although a single peak in the heat flow signal characterizes this step, the mass change signal clearly shows two consecutive stages. The first of these two stages starts from 430 to 470°C, and this proceeds relatively fast showing a very steep conversion versus temperature curve. The second stage which proceeds at a relatively slower rate continues from 470°C and ends at about 520°C. A last oxidation stage starts from 530°C and ends at around 960°C, where the remaining TiC, oxycarbide and anatase are converted into rutile. A similar reaction scheme was observed and reported earlier by [12] when studying the oxidation behaviour of TiC powder at various oxygen partial pressures. However, in that study the second stage was split up into the formation of anatase and rutile; in addition, it was mentioned that these two stages may overlap and cannot be clearly distinguished from each other.
To continue, the other powders were all heated at 5°C min -1 in air in separate experiments. The results are presented in Fig. 5. As observed for the 5 lm powder, multiple peaks are identified for the different particle sizes. The stages of oxidation were resolved by running separate experiments up to each identified peak in the heat flow or the differentiated mass change signal, and then the products were analysed by XRD and/or XMA.
For the coarsest powders having an average particle size of 140 and 48 lm, only the peaks pertaining to the formation of anatase and rutile can be identified in the heat flow and mass change signals (see Fig. 5a, b, respectively). The stage of oxycarbide formation appears to be overshadowed. For both powders, the peak occurring between 400 and 600°C corresponds to the formation of anatase. The formation of rutile sets in at 800°C and ends at 1195°C for the 140 lm powder, and for the 48 lm powder it starts from 700°C and ends around 1050°C. It is noteworthy that the transformation into rutile of both coarse powders occurs via multiple steps while the mass gain signal for this range reflects a single peak event (i.e. a non-changing slope).
For the 0.8 lm TiC powder, all three stages observed for the 5 lm powder can be identified in the heat flow signal. The formation of the oxycarbide is shown in Fig. 5c peaking at about 310°C, as well as in the differentiated mass gain signals shown in Fig. 6. Anatase is formed in the sharp peak around 420°C followed immediately by conversion into rutile which also occurs via two peaks.
For the nano-powder the mechanism somewhat reverses to that of the coarser powders as no trace of the oxycarbide formation can be detected. A peak corresponding with the formation of anatase can be identified at 400°C, and two peaks associated with the formation of rutile can be identified at about 430 and 550°C, respectively, in Fig. 5d. The various stages of oxidation, as identified for the different particle sizes, are summarized as a function of temperature in Fig. 7. Depending on the particle size, some stages may overlap or is overshadowed and hence not always clearly resolved in the mass gain and heat flow signals as a function of temperature.

Effect of heating rate on oxidation
As the heating rate is increased, the amount of TiC transformed into anatase is reduced. For example, for the 140 lm powder the amount of TiC transformed decreases from 4 to 1% as heating rate is increased from 1 to 10°C min -1 . Similarly, there was a decrease of about 5-1.5% and 40-15% for the 48 and 5 lm powders, respectively, during the same stage, see Fig. 8.
For the same amount of powder, the surface area exposed to oxygen increases when the TiC particle size is reduced. Then, the conversion of TiC into TiO 2 is enhanced. This is evident from the heat flow signal recorded with the same heating rate, which increases for smaller particles. For example, the formation of rutile when oxidizing 5 lm TiC powder at 10°C min -1 corresponds with about 45 lV of heat evolved, while under the same conditions 140 lm TiC powder results in about 28 lV of heat evolved; cf. Fig. 8a, e. Since the oxidation of TiC powder is faster when the particle size is smaller, the final temperature for full conversion of TiC powder into TiO 2 (rutile) during non-isothermal oxidation (at the same heating rate) decreases with particle size; see Fig. 8.
A final observable trend is that the stages of oxidation are well separated from each other (in both the heat flow and the fraction converted signals) at larger particle sizes. However, as the particle size is reduced, these stages become indistinguishable in the c Figure 8 Heat flow and mass change in terms of converted fraction TiC (in %) as a function of temperature recorded at different heating rates in dry synthetic air of TiC powders with average particle size of: a, b 140 lm, c, d 48 lm, e, f 5 lm, g, h 0.8 lm and i, j 0.05 lm.

Figure 7
Chart showing the various stages of oxidation in dry synthetic air of TiC particles with different average particle sizes. fraction converted plots especially for the 50 nm powder.

Evaluation of the activation energy
The activation energy for oxidation of the different TiC powders is determined using Eq. (8). As more than one peak occurs during the transformation reaction, the peak temperature (T p ) is defined as the highest point in the differentiated mass gain signal. The results for the various TiC powders are presented in Fig. 9. With decreasing particle size, the lines for the maximum conversion of TiC into TiO 2 shifts to lower temperatures, and this defines the temperature window where a particle is activated. The slope of these lines corresponds to the activation energy of the reaction, and the figure clearly shows it decreases with decreasing particle size.
The values of the activation energy E A as a function of the TiC particle size are presented in Table 2. It is observed that E A significantly decreases with decrease in particle size for submicron and nanosized particles and is more or less constant for particle sizes [ 5 lm. The activation energies determined here are comparable to those determined in other studies (e.g. [11,25]), i.e. 192 and 245 kJ mol -1 . The (non-significant) difference may arise from the size of the TiC particles and the limited temperature window of the isothermal experiments used in the present analysis.
Effects of particle size on the activation energy have been reported for other materials as well [26,27]. It can be conceived that when the surface-tovolume ratio increases the activation energy decreases, since the atoms at the surface have fewer bonds and neighbouring atoms and thus need less energy to transfer into oxide [28].
Hence, the reduction of the bulk activation energy, E bulk A , can be related to the surface energy r according to [29]: in which V m is the molar volume and d the particle size. For TiC the surface energy equals 2.73 J m -2 [30] and the molar volume 12.18 cm 3 mol -1 [31]. From Fig. 10, the observed reduction in the activation energy for submicron or nano-sized particles is due to surface defects (like: dislocations, grain boundaries, steps and kinks, roughness), of which their number increases with surface area (i.e. with d -2 ) [32]. These surface defects promote the nucleation and growth of the oxide.   Hence, the activation energy as a function of particle size is presented in Fig. 10 and the relationship can be best described with Eq. (16).
where E A is the activation energy in kJ mol -1 and d is the particle size in nm.

Reaction model for TiC oxidation
Out of the many reaction models [18], the nonisothermal plots of the rate of conversion with temperature pre-empt that the oxidation of TiC may be either an Avrami-Erofeyev, diffusion, or an order reaction model. Hence, the experimental data are compared with the Avrami-Erofeyev (A3), 3D diffusion Jander (D3), Ginstling-Brounshtein (D4) and the third-order (F3) reaction models using the master curve plotting method [23]; see ''Reaction model'' section. To ascertain the validity of the identified model, master plots are made for all the powders at different heating rates, see Fig. 11. It is seen that for different heating rates the experimental data closely resemble the F3 (third-order reaction model) more than any other model. The integral and the differential form of the third-order reaction model is, respectively, defined as: The rate-determining step in the oxidation of TiC is the diffusion of oxygen from the gaseous phase through the formed oxide layers and the opposite diffusion of titanium and carbon from the material into the scale.

The Arrhenius pre-exponential factor
The last of the kinetic triplet, i.e. the pre-exponential factor (A), is evaluated from the integral form of the kinetic rate equation given in Eq. (11) combined with Eq. (17). The generalized time Eq. (9) is determined at every temperature from 25 to 1200°C using Eq. (14), the applied heating rate and the estimated E A . A linear fitting of Eq. (11) is made for the range of a pertaining to the third oxidation stage (i.e. transformation into rutile) to find the constant A. The calculated pre-exponential factors, as shown in Table 3, are relatively constant with the same order of magnitude. However, it significantly decreases for the nanometre-sized particles [33,34]. As the particle size decreases, the partial molar surface enthalpy and the partial molar surface entropy increase, leading to the decrease in the apparent activation energy and the pre-exponential factor [29,35]. There is no simple explanation for the high value of the Arrhenius parameter for the 5 lm particle, other than undetected impurities or structural defects.
Finally, the conversion a of TiC into TiO 2 by oxidation in air can be expressed as a function of temperature and time for each particle size by: Table 3 Evaluated preexponential factors (A) for TiC powders oxidized in dry synthetic air
a ¼ 1 À ½ð2Ah þ 1Þ ðÀ1=2Þ ð 19Þ with h calculated using Eq. (14). It can be seen in Fig. 12 that at higher temperatures the calculated results fit well with the experimental results with a slight deviation at lower temperatures. This deviation is associated with initial stages of oxidation which corresponds to the formation of the oxycarbides and/ or anatase.

Conclusions
Oxidation of TiC was found to proceed via the formation of oxycarbides ? anatase/rutile ? rutile and closely follows a third-order reaction path. The temperature window for activation of the oxidative decomposition of the TiC healing particles, i.e. for triggering the self-healing reaction, drops from 1000 to 400°C as particle size decreases from 140 lm to 50 nm. As the change in optimal healing temperature does not depend on the chemical composition of the matrix material as such, the present study shows in a generic manner that extrinsic self-healing in oxide ceramic matrices can be tailored to a specific temperature range by tuning the size of the healing particle. Further, the activation energy required for oxidation decreases significantly for particles ranging between 10 lm to 50 nm. Above 10 lm particle size, the activation energy for the thermal oxidation of TiC becomes constant.