A Material-by-Design Approach to Develop Ceramic- and Metallic-Particle-Reinforced Ca-α-SiAlON Composites for Improved Thermal and Structural Properties

α-SiAlON is commonly used to machine superalloys owing to its desirable thermal and structural properties. α-SiAlON is among the crystalline forms of SiAlON and has more favorable properties than β-SiAlON. However, it becomes fragile during the machining of hard-to-cut materials due to its low fracture toughness and machinability. Recent research efforts focus on improving the thermal and structural properties of α-SiAlON using suitable dopants, nano-sized precursors, and the addition of metallic/ceramic reinforcement particles. The present study presents a material-by-design approach to designing and developing ceramic and metal-particle-reinforced Ca-α-SiAlON composites with properties tailored for the cutting tool applications. The mean-field homogenization theories and effective medium approximations implemented in an in-house code are used to effectively optimize the thermal and structural properties of the Ca-α-SiAlON composite by varying essential parameters such as inclusion material, volume fraction, porosity, particulate size, and thermal interface resistance. Individual properties of the matrix and reinforcements are considered in the computations of effective properties such as thermal conductivity, thermal expansion coefficient, modulus of elasticity, and fracture toughness. The main objective of the study is to enhance the thermal conductivity and fracture toughness of Ca-α-SiAlON, while lowering its thermal expansion coefficient. At the same time, the elastic modulus and hardness/strength must be maintained within an acceptable range. As a validation, Ni/Ca-α-SiAlON and SiC/Ca-α-SiAlON composites are synthesized from the nano-sized precursors, CaO dopant, and Ni/SiC microparticles via spark plasma sintering (SPS) process. The thermal conductivity, coefficient of thermal expansion, and elastic modulus of the composites are measured and compared with the computational predictions. The computational predictions are found to be comparable to that of the experimental measurements. Moreover, the studies show that WC, SiC, and Cr can be suitable reinforcement materials for enhancing the thermal and structural properties of Ca-α-SiAlON material for the cutting tool inserts.


Introduction
Owing to its excellent thermomechanical and tribological characteristics, SiAlON has become a suitable ceramic material for a wide range of industrial applications, such as bearings, heat insulators, sealing elements, engine parts, cutting tools, refractories, etc. [1][2][3]. The use of this material originated in the early 1970s when Oyama and Kamigaito [4] and Jack and Wilson [5] reported that Al 3+ could be infused into the lattice of Si 3 N 4 without altering the crystal structure, i.e., by replacing Si 4+ such that the N 3− are also replaced by O 2− , to form the solid solution, SiAlON. SiAlON is thus formed from the chemical reaction wear resistance in SiAlON. However, the strength of the nanocomposite was considerably reduced due to the high interfacial mismatch strain developed [31,32]. Similarly, the hardness and fracture toughness of SiAlON can be improved by adding cBN inclusions [33] as well as WC particles [34]. As the thermal conductivity of SiC is significantly higher than that of SiAlON, adding 10-25 wt% SiC improves the thermal diffusivity of α-SiAlON by~30-70%. At the same time, the thermal diffusivity of β-SiAlON improves by about 19% [35]. The thermal diffusivity increases with increasing SiC particles due to increasing the path for phonon transport. A recent study showed that the mechanical properties of SiAlON become further improved with increasing SiC particle size [36][37][38]. On the other hand, the addition of metallic inclusions (such as Ni) can improve the thermal conductivity and fracture toughness, even though the hardness and the thermal expansion coefficient become affected due to higher porosity and partial melting of metallic particles during sintering [39].
Although the SiAlON-based composite has been extensively studied in previous research, there are no comprehensive studies investigating the influence of a wide range of particle reinforcement types, volume fractions, and particle sizes to date. This impedes the progress of utilizing SiAlON-based composites as a replacement for pure SiAlON or other ceramic materials for various industrial applications. Therefore, there is a need for a critical evaluation of the influence of metallic and ceramic particles on the effective properties of SiAlON-based composites. For a more systematic study, computational design tools have become necessary to minimize the costs associated with the trial-and-error approach utilized experimentally.
Previous studies show that the CaO additive resulted in a single-phase α-SiAlON ceramic with a relatively more stable/dense microstructure and a suitable range of thermomechanical properties compared to when the other additives were used [15]. In addition, CaO is highly soluble, cheaper, and widely available from mineral resources such as fly ash; hence, Ca-α-SiAlON is chosen for the current study. The present study used a computational design approach to tailor Ca-α-SiAlON composite properties for cutting tool applications by reinforcing it with a wide range of metal and ceramic particles. Effective medium and mean-field homogenization schemes implemented in our in-house codes were used to optimize the effective properties of Ca-α-SiAlON composites to obtain superior properties that suit the intended application. The predicted properties included thermal conductivity, coefficient of thermal expansion, modulus of elasticity, and fracture toughness. The effect of the inclusion volume fraction, particle sizes, thermal interfacial resistance between the inclusion and matrix, and porosity are studied. The parameters are optimized for the thermal conductivity, fracture toughness, and thermal expansion coefficient of the Ca-α-SiAlON composite while maintaining the modulus of elasticity at an acceptable level.
The most suitable inclusion materials were selected based on the optimized predictions from the computational studies. To validate the predictions, the designed particle volume fractions and particle sizes were used to synthesize Ni/Ca-α-SiAlON and SiC/Caα-SiAlON composites from the nano-sized precursors, CaO dopant, and Ni/SiC microparticles via the spark plasma sintering (SPS) process. The microstructural morphologies of the synthesized samples were assessed via scanning electron and optical microscopy, while phase characterization was conducted via X-ray diffraction analysis. Porosity was measured via the Archimedes technique, while the hardness and elastic modulus were measured using the instrumented indentation. The thermal conductivity and coefficient of thermal expansion were measured using the thermal and thermo-mechanical analyzers, respectively. The measured properties of the composite materials were then compared with that of the computational predictions.

Computational Material Design
The computational design of composite materials has become feasible with the recent advancement in computing and homogenization theories [40]. Composite materials contain certain microstructural constituents (such as inclusions, fibers, defects, pores, and Nanomaterials 2022, 12, 2176 4 of 24 interfaces) that must be selected and dispersed to attain the desired properties. Computer simulations can be used to assess the thermomechanical response of composite materials under various combinations and forms of microstructural constituents before their synthesis and development; hence, they serve as cost-effective design tools for tailoring the desired properties.

Mean-Field Homogenization Using Effective Medium Theories
Homogenization theories are used to estimate the effective properties of heterogeneous materials through computational modeling or effective medium approximations. Effective medium theories (EMTs) are used to predict effective properties by utilizing various analytical or theoretical relations that describe the macroscopic behavior of composite materials. With EMTs, acceptable relations (or approximations) are used to predict the microscopic behavior of a composite system based on the individual properties and volume fractions of individual microconstituents. Subsequently, the macroscale response of the composite material is obtained by averaging the microscale results obtained with the theoretical relationships.
In the present study, multi-inclusion effective medium theories were used to predict the thermo-mechanical properties of various composite materials ( Figure 1). The properties of the matrix and inclusion materials, particle size and shape, particle orientation, volume fraction, and interfacial thermal resistance were adopted in the calculations. The effective thermal conductivity of the composite material was predicted using the multi-inclusion effective medium approximation (EMA) model by Siddiqui and Arif [41] as well as theoretical models by Hasselman and Johnson [42]. The elastic modulus and coefficient of thermal expansion were predicted using the Mori-Tanaka mean-field homogenization scheme. In contrast, fracture toughness was predicted using the energy-based semi-empirical model previously developed by Li and Zhou [43,44].
The Hasselman and Johnson model [42] is among the earliest analytical models used to compute the effective thermal conductivity (κ e ) of a composite material by taking into account the thermal boundary resistance between the matrix and reinforcement phases as well as the other essential composite parameters such as particle shape, particle radius, and volume fractions (Equation (1)). The model has been proven to give accurate predictions for a composite with dilute fractions of the reinforcement material [45]. However, to apply the model to a hybrid composite with the second inclusion and/or any porosity, it is necessary to use the two-step calculation procedure where the porosity or second inclusion is used to arrive at the effective conductivity for the matrix phase [46]: where κ e f f mat = 2k m 1−ϕ p 2+ϕ p is the effective thermal conductivity of matrix material, k m is the matrix thermal conductivity, ϕ i is the volume fraction of the first inclusion, ϕ p is the porosity/second inclusion fraction, κ is the effective thermal conductivity of the first inclusion, k i inc is the thermal conductivity of the inclusion, R i TB is the thermal boundary resistance, and a i is the inclusion radius.
Similarly, the multi-inclusion effective medium approximation (EMA) model [41] assumes that: (i) the inclusion volume fraction is small and sufficiently dilute, (ii) heat transfer across the imperfect mating interfaces is affected by the interfacial thermal resistance, (iii) porosity is the second phase inclusion with a specific size, distribution, shape, and volume fraction. The advantages of the EMA model over the Hasselman and Johnson model are the considerations for (i) nanosize effect, (ii) multiple inclusions, (iii) different particle orientation and shape, and iii) non-uniform dispersion of inclusions. Nevertheless, both theories are developed based on the concept of Kapitza radius, which represents the mini-mum threshold particle radius beyond which the thermal boundary resistance becomes negligible in the calculations of the effective thermal conductivity. The multi-inclusion EMA model can be expressed as [40]: where: is the orientation factor for the inclusion i, p i = a i 3 /a i 1 is the particle aspect ratio, l i , a i k are the radii that define the particle shapes, a kap = R i TB κ mat is the Kapitza radius, and α i k = a kap /a i k is a dimensionless quantity that represents the influence of thermal boundary resistance on heat flow, R i TB is thermal interface resistance, κ mat is the matrix thermal conductivity, κ e are the overall thermal conductivities along the three orthogonal axes, κ i inc is the thermal conductivity of the inclusion i and ϕ i is the volume fraction of inclusions.
Eshelby's research in 1965 provides the fundamental framework for predicting effective stiffness for inhomogeneous elastic media. According to Eshelby, the microstrain distribution developed within inhomogeneous elastic media relates to the macroscale strains via the strain localization equation. Based on Eshelby's theory, several schemes are used to predict the effective properties of composite materials. Among other schemes, the Mori-Tanaka scheme optimally fits heterogeneous materials with inclusions (of volume fraction <30%), which may interact with one another. Using the Mori-Tanaka mean-field homogenization scheme [47], the effective elastic modulus and thermal expansion coefficient are estimated using the following equations [48]: A i = B a : A mat (10) where C refers to the elastic modulus, α is the thermal expansion coefficient, A is the strain localization tensor, ϕ is the volume fraction, S is the Eshelby tensor, I 2 is a second-order identity tensor and I 4 is the fourth-order identity tensor. The inclusion and matrix are designated with the subscripts i and mat, respectively.
where refers to the elastic modulus, is the thermal expansion coefficient, is the strain localization tensor, is the volume fraction, is the Eshelby tensor, is a second-order identity tensor and is the fourth-order identity tensor. The inclusion and matrix are designated with the subscripts and , respectively. In their previous works, Doghri and Tinel [49] showed that the direct application of the Mori-Tanaka scheme for composites with dual inclusions could lead to misleading results. Hence, a modified scheme where the composite material is subdivided into pseudo-grains that contain only one type of inclusion was used in the current work. As shown in Equation (13), the required properties are more accurately estimated by applying the Voigt method to carry out homogenization across the entire grains: The dominant fracture mechanism in composite materials includes particle cracking, matrix cracking, and interface debonding. The effective fracture toughness ( ) of the ceramic composite is predicted by the effective elastic moduli ( ̅ ) and critical energy release rates ( ) of the ceramic composite based on Equations (14) and (15). The critical In their previous works, Doghri and Tinel [49] showed that the direct application of the Mori-Tanaka scheme for composites with dual inclusions could lead to misleading results. Hence, a modified scheme where the composite material is subdivided into pseudograins that contain only one type of inclusion was used in the current work. As shown in Equation (13), the required properties are more accurately estimated by applying the Voigt method to carry out homogenization across the entire grains: The dominant fracture mechanism in composite materials includes particle cracking, matrix cracking, and interface debonding. The effective fracture toughness (K IC ) of the ceramic composite is predicted by the effective elastic moduli (C e ) and critical energy release rates (J ic ) of the ceramic composite based on Equations (14) and (15). The critical energy release rate is predicted using the analytical model by Akhtar et al. [50], which was developed based on the original model by Li and Zhou [36,37] for the prediction of fracture toughness in ceramic composites. The developed models assume fracture by quasi-static crack growth, and the ultra-low crack growth rate is predicted through some dynamic calculations. Hence, the evaluated fracture toughness depends on essential microstructural features such as inclusion material type, inclusion size and distribution, individual fracture Nanomaterials 2022, 12, 2176 7 of 24 toughness of matrix and inclusion materials, interfacial fracture energy, etc. Details about the fracture toughness model are given in the work of Akhtar et al. [50]: where v is the effective Poisson's ratio of composite system. Due to cracking by multiple modes, the critical energy release rate is expressed as the average energy release rate due to particle-cracking, matrix-cracking, and interface-debonding: is the crack length multiplication factor, D 1 and D 2 are functions of inclusion particle sizes, ϕ is the fracture energy, H is the proportional crack length, in denotes the inclusion material, mat denotes the matrix, and i denotes the inclusion particle. The proportional crack lengths can be estimated from the empirical relations: where D is the total crack length, p is the probability of crack deflection, Q is the strength ratio and P ij is probabilistic functions for matrix and particle cracking. The Ca-α-SiAlON properties are enhanced by adding a single inclusion from various materials, sizes, and distributions to the ceramic matrix to attain better cutting performance and tool life. The desired properties considered include high thermal conductivity, low coefficient of thermal expansion (CTE), moderate elastic modulus, and high fracture toughness. Table 1 highlights the properties of Ca-α-SiAlON and the various combinations of metallic and ceramic inclusion materials considered for the present study. Three metallic inclusion materials, including Ni, Co, and Cr, were selected due to their high fracture toughness and high electrical and thermal conductivities. However, compatibility with the ceramic material becomes an issue due to the high structural mismatch at mating interfaces resulting from the considerable variation in CTE between the metallic reinforcement and ceramic matrix materials. With the use of ceramic filler materials (such as the those provided in Table 1), the structural mismatches at the mating interfaces do not usually occur; hence, ceramic inclusion materials that possess relatively higher electrical and thermal conductivity are selected for the present study (i.e., WC, SiC, ZrB 2 , cBN, and TiN). The main issue with the ceramic class is their low fracture toughness, even though their dispersion within the ceramic matrix is expected to increase the crack growth resistance.

Experimental
Following the computational design of the composite materials, a few Ca-α-SiAlONmatrix composites were reinforced with the SiC and Ni particles that were synthesized for validation. The details describing the materials, synthesis, testing, and characterization are discussed in the following sections.
To validate the computational predictions, SiC (as a ceramic reinforcement) and Ni (as a metal reinforcement) were selected to synthesize Ca-α-SiAlON matrix composites with varying contents (0 wt%, 10 wt%, 20 wt%, 30 wt%) in line with the computational predictions. To prepare the powder mixtures for each sample, the starting precursors were carefully weighed. Table 2 shows the details of the prepared composites. Each sample's powder blends were placed in a methanol medium and then dispersed for around 20 min by ultrasonic probe sonicator (VCX 750, Sonics, Newtown, CT, USA). After homogenization, the samples were placed in a muffle furnace for around 24 h at a temperature of 80 • C for drying. After they had dried, the homogenized mixtures were crushed using a granite pestle and mortar for ease of consolidation. Subsequently, the homogenized mixtures were consolidated in a 20 mm graphite die for 30 min at 1500 • C, 100 • C/min rate, and a uniaxial pressure of 50 MPa by SPS (HP D5, FCT Systeme, Gewerbepark, Frankenblick, Germany). Experiments were conducted at a pressure of 5 × 10 −2 mbar. To remove the graphite embedded on the surface of the synthesized samples, the samples were ground using SiC abrasive paper with grit sizes ranging from 60 to 1000 grit sizes. Moreover, an alumina suspension (particle size of 9 µm) was used to obtain mirror-like surfaces.

Characterization and Testing Details
The 18 mm synthesized samples had their thermal conductivities measured at room temperature by applying continuous heat using a one-sided interfacial heat reflector (TCi, C-Therm, Fredericton, NB, Canada). A diamond cutter was used to cut each of the four samples into three pieces, one of which was utilized to measure the coefficient of thermal expansion (using TMA/SDTA 1 L.F./1100, Mettler Toledo, Columbus, Ohio, Switzerland). The Vickers hardness of the SiAlON samples was determined using instrumented indentation (ZHU250, ZwickRoell, Horsham, Pennsylvania, Germany) at a load of 10 kg. The backscattered mode of a field emission scanning electron microscope (FESEM, FEI, Hillsboro, OR, USA) was used to examine the microstructure of the samples. It was equipped with an energy-dispersive X-ray spectrometer (EDX, Malvern Panalytical, Worcester, MA, USA), which aided in collecting compositional data. The Archimedes principle was used to determine the density of the sintered samples.
The Evans and Charles Equation (Equation (19)) was used to compute the fracture toughness (K IC ), where H is the Vickers hardness, a is the length from the center of the indentation to its edge, and c is the length from the end of the vertical crack line to the end of the horizontal crack line.

Computational Prediction of Effective Properties
Computational prediction of the effective properties is based on effective medium approximations for various combinations of inclusion materials, volume fractions, and porosity utilizing the various material properties obtained from the literature. This section summarizes the findings of the computational study conducted with the mean-field homogenization codes. Firstly, the effective properties were predicted for selected reinforcement materials and several combinations of inclusion volume fraction, size, and porosity to identify the threshold for real improvement in the required properties. The porosity and inclusion particle size were optimized, while the volume fraction range was determined. The inclusion particle size is solely optimized utilizing the thermal conductivity code. Simultaneously, the porosity is fixed at a particle level based on the results of effective thermal conductivity, elastic moduli, and CTE. The final stage of the computational study involved running the simulations for varying volume fractions, fixed porosity, and particle sizes to determine the effective properties resulting from the addition of each inclusion material shown in Table 1. The data generated at this stage was used to hone the effective composite properties and identify the most suitable inclusion materials to improve cutting tool life and performance. Figure 2 illustrates the effective thermal conductivities predicted for the Ca-α-SiAlON composites as a function of the volume fractions of several ceramic and metal reinforcements following Table 2. Overall, it can be observed that the thermal conductivity increases proportionally to the volume fractions, i.e., as the volume fraction of the inclusions increases, the thermal conductivity of the composite increases because of the filler material having higher thermal conductivity. The diamond-reinforced composite had the most significant rise in thermal conductivity, followed by the cluster of SiC, Co, Cr, Ni, WC, ZrB 2 , and finally, cBN and TiN. The rise in the thermal conductivity due to the addition of diamond inclusion is the highest because it has the highest thermal conductivity (~2000 W/m.K). TiN, predictably, displays the lowest rise in the composite's thermal conductivity, which can be attributed to it having the lowest thermal conductivity (~19 W/m.K) among the remaining inclusions. Hence, by considering the highest and the lowest thermal conductivity rises as exceptional cases or outliers, the region of the aforementioned cluster is deemed the most appropriate to aid in the selection of the most optimal effective thermal conductivity (~8 W/m.K) with respect to its volume fraction (~16%).
It is worth mentioning that the aforementioned thermal conductivity trends were obtained at a fixed interfacial resistance (ITR) value of 5.0 × 10 −8 K.m 2 /W, a porosity of 2%, and a particle size of 10 µm. The average value for the ITR between the inclusions and Ca-α-SiAlON matrix is 5.0 × 10 −8 K.m 2 /W as obtained from the literature [51]. The ITR value between the inclusion and SiAlON matrix depends on the wettability of the two phases, which can be controlled by proper functionalization for any desired value. Hence, it is possible to further the thermal conductivity by appropriate surface treatment of particles if deemed necessary. Although an average ITR value was selected for our model, the relationship between ITR and particle size was essential in governing the effective thermal conductivity of metal/ceramic-reinforced ceramic matrix composites. As the inclusion particle size decreased, the ITR increased because a reduction in particle size increased contacted the surface area. The Kaptza radius for the Ca-α-SiAlON composite is about 0.29 µm; thus, the thermal boundary resistance created by the 10 µm size inclusions does not notably influence the thermal conductivity values. However, the porosity and volume fractions directly affected the effective thermal conductivity. The thermal conductivity decreased with increasing porosity, which was associated with the high thermal resistance due to the reduction of gas-phase heat conduction. Figure 2b shows that the predictions made with the model by Hasselman and Johnson slightly overestimated the effective thermal conductivities of the composite materials. This can be attributed to the incorporation of porosity (as the second inclusion) via the two-step approach with the model by Hasselman and Johnson. On the contrary, the effective medium approximation (EMA) model implicitly considered multiple inclusions and porosity in the formulation; hence, it could lead to a better estimation of the effective thermal conductivities. Nevertheless, the predictions made with the two theories considered here were in very close range for the 10 µm particle size adopted here. Figure 3 was plotted in accordance with the determined volume fraction of 16% along with the average ITR value of 5.0 × 10 −8 K.m 2 /W to predict the change in the effective thermal conductivity with respect to inclusion size and to determine the threshold particle size value. The general trend observed is an increase from the plain Ca-α-SiAlON thermal conductivity value of 7.6 W/mK only after an inclusion size of approximately 1 µm. This result was taken as the minimum threshold particle size post from which the ascension of thermal conductivity can be observed. Additionally, the saturation points for all the reinforced composites coincided at roughly 130 µm. Although increasing the volume fractions and reducing the porosity results in higher thermal conductivity, other effective properties, such as the elastic modulus and CTE, need to be considered for benchmarking inclusion volume fractions. Figure 4a depicts the predicted effective elastic modulus for various ranges of volume fractions and the reinforcement materials stated in Table 2. It can be observed that the addition of reinforcements from the cluster of diamonds through SiC increases the elastic modulus, with the diamonds providing the sharpest rise. This can be attributed to them showcasing higher elastic moduli than the plain Ca-α-SiAlON matrix. In contrast, the grouping of Cr through Ni shows a decline in the elastic modulus for the opposite reason (comparatively lower elastic modulus than plain Ca-α-SiAlON). Moreover, cBN's elastic modulus is akin to that of plain SiAlON, so the rise in elastic modulus is only marginal. Figure 4b shows the change in the effective coefficient of thermal expansion (CTE) with respect to the volume fraction of the nine previously mentioned reinforcements. The composite must maintain a low CTE and high thermal conductivity to ensure that the thermal loads are low and heat dissipation is high, respectively, for improved cutting tool performance. In this regard, SiC, cBN, and Cr were deemed to be the most appropriate reinforcers as they do not cause the CTE of the composite to upsurge dramatically with an increase in inclusion volume fraction. Diamond would be the best candidate in this regard, but it is ignored due to its high cost. Ni through ZrB 2 , on the other hand, can be considered unsatisfactory for our target application due to an increase in CTE with the addition of more inclusion material. It is already established that the CTE and material elastic modulus are interrelated. Materials having strong interatomic interactions, such as ceramics, have low CTEs and high elastic moduli, making them stiffer. Hence, by adding metal inclusions such as Ni and Co in Ca-α-SiAlON, which has high CTE, but low elastic modulus, the CTE is expected to increase rise, as predicted by the model, and verified experimentally.
At high cutting speeds, increased fracture toughness can be beneficial to the cutting tool's durability and lifetime. The computational model is initially performed for various potential ceramic reinforcements to determine the most appropriate candidates likely to improve the effective K IC of the resulting Ca-α-SiAlON composite. Figure 5 shows the effective K IC as a function of various ceramic inclusions in the Ca-α-SiAlON matrix, assuming a particle size of 10 µm in each combination. The shape of the inclusions was assumed to be perfectly spherical, and the aspect ratio was kept constant. It is worth noting that the K IC of generic Ca-α-SiAlON [52] (~4.67 Mpa.m 1/2 , which corresponded to fracture energy of 60 J/m 2 ) was set as a benchmark to assess the selected reinforcement types. Metal reinforcements are excluded from the analysis as the existing computational model for fracture toughness is valid only for ceramic particles in ceramic matrices. It can be observed that the effective K IC increases as the volume fraction of inclusions increases in general. However, K IC drops after around 20% volume for inclusions having low fracture resistance, such as diamond and cBN. This peak value of K IC shifts towards a higher volume percentage value for inclusions with higher fracture resistance. Hence, the WC with the highest fracture energy (85 J/m 2 ) leads to the highest K IC with a continuously increasing trend, followed by TiN and ZrB 2 . The decreasing trend in K IC beyond a specific volume concentration (in the case of reinforcements with low fracture energies) is attributed to the dominant particulate cracking mode, unlike interface cracking at high volume percentages in the Ca-α-SiAlON matrix. As a result, this tendency shows that crack deflection, rather than particle cracking, can be fostered. Therefore K IC can be improved by including inclusions with high fracture energies for a well-bonded interface, as considered in the current model and reported by Kumai et al. [53]. As a result, particles with high fracture energy can be employed to increase the K IC of SiAlON composites.

Characterization and Morphology of Synthesized Samples
The interaction of reinforcements with the SiAlON matrix directly impacts the final properties such as thermal conductivity, expansion coefficient, stiffness, and fracture toughness. The accuracy of the computational predictions was validated through an experimental phase. This included determining the composite samples' characteristics and characterization. To evaluate the computational design accuracy, a variety of SiC and Ni-reinforced Ca-α-SiAlON composites were developed in line with predictions using spark-plasma sintering. A comprehensive microscopic study of the sintered sample was performed to describe distinct composite combinations.
FESEM images of the polished surfaces of the sintered 10% SiC, 20% SiC, and 30% SiC Ca-α-SiAlON composites are shown in Figure 6a-c, respectively. The homogenous dispersion of the white SiC particles in the black Ca-α-SiAlON matrix demonstrates that the ultrasonic probe sonication was sufficient to achieve appropriate SiC particle homogeneity. FESEM micrographs of fractured surfaces of monolithic Ca-α-SiAlON and 30%SiC/Ca-α-SiAlON composite samples are shown in Figure 6d,e, respectively. As can be observed, a typical microstructure of single-phase α-SiAlON with equiaxed grains are evident in the case of a monolithic SiAlON sample. An intergranular form of fracture in monolithic SiAlON is evident without any major pull-out of particles. On the contrary, sharp step-like structures inside the SiC/SiAlON composite indicate cleavage or delamination in the grains. Grain pull-out and cleavage stages have been identified as essential variables in improving the fracture toughness of ceramic composites with increasing SiC content [54].  Figure 8 shows the XRD pattern of the 30 wt% Ni/Ca-α-SiAlON composite at the various sintering temperatures. The mixture of the precursors was heated up to 1500 • C during the high-temperature XRD under an argon atmosphere. The crystalline phases were continually noted during the heating process. Si 3 N 4 is the dominant constituent of the powder mixture, and a major peak of Ni was visible at 100 • C and a diffraction angle of 44.5 • . The diffractogram shows that the presence of Ni reinforcements particles did not react to form any new compounds such as the aluminides, silicides, etc. However, a notable change in the XRD pattern is observed at about 1200 • C and 1400 • C with the formation of the liquid phase Ca-α-SiAlON and solution-diffusion-precipitation phenomenon, respectively. The AlN and the intermediate phase(s) disappear at temperatures above 1400 • C, while the formation of the α-sialon increases during the heating process until reaching completion at about 1500 • C.

Validation of Predictions
Experimental measurements of the composite properties such as thermal conductivity, CTE, hardness, and fracture toughness, were compared to the numerical predictions. The base properties of sintered Ca-doped α-SiAlONs were first determined experimentally and then used as input in the computational modeling for recalibration and realistic predictions.
The measured and predicted thermal conductivity, CTE, and hardness of monolithic Ca-α-SiAlON and 10% Ni-, 20% Ni-, and 30% Ni-reinforced Ca-α-SiAlON produced are shown in Table 3. The porosity of each sintered sample is also presented. The difference between the predicted and experimental values of thermal conductivity and CTE is less than 8%, which is deemed ab acceptable agreement. The discrepancy between experimental and numerical results can be attributed to several assumptions. The spherical shape of inclusions, uniform particle size distribution, and perfect dispersion of Ni inclusions inside the SiAlON matrix are among these assumptions. More importantly, the data on interface thermal resistance was collected from the published literature. With an increase in Ni content, the porosity of the composites was found to increase, which was also confirmed by the decrease in the relative density of the resulting composites. This decrease in density and increase in porosity is associated with the softening and/or melting of Ni particles during the composite's synthesis process. As a result, the mismatch of the CTE between Ca-α-SiAlON matrix and Ni results in the formation of voids at the interface. Furthermore, deposits of Ni were found on the exterior of the synthesized samples, demonstrating that some of the Ni particles melted and exited the graphite die during the synthesis. As the nickel content increased, the amount of nickel that melted increased,

Validation of Predictions
Experimental measurements of the composite properties such as thermal conductivity, CTE, hardness, and fracture toughness, were compared to the numerical predictions. The base properties of sintered Ca-doped α-SiAlONs were first determined experimentally and then used as input in the computational modeling for recalibration and realistic predictions.
The measured and predicted thermal conductivity, CTE, and hardness of monolithic Ca-α-SiAlON and 10% Ni-, 20% Ni-, and 30% Ni-reinforced Ca-α-SiAlON produced are shown in Table 3. The porosity of each sintered sample is also presented. The difference between the predicted and experimental values of thermal conductivity and CTE is less than 8%, which is deemed ab acceptable agreement. The discrepancy between experimental and numerical results can be attributed to several assumptions. The spherical shape of inclusions, uniform particle size distribution, and perfect dispersion of Ni inclusions inside the SiAlON matrix are among these assumptions. More importantly, the data on interface thermal resistance was collected from the published literature. With an increase in Ni content, the porosity of the composites was found to increase, which was also confirmed by the decrease in the relative density of the resulting composites. This decrease in density and increase in porosity is associated with the softening and/or melting of Ni particles during the composite's synthesis process. As a result, the mismatch of the CTE between Ca-α-SiAlON matrix and Ni results in the formation of voids at the interface. Furthermore, deposits of Ni were found on the exterior of the synthesized samples, demonstrating that some of the Ni particles melted and exited the graphite die during the synthesis. As the nickel content increased, the amount of nickel that melted increased, resulting in lower relative densities. This was confirmed when the actual amount of Ni in the final composite was compared to the initial value used in calculating theoretical density. Table 3. A comparative analysis of computational prediction and experimentally measured thermomechanical properties of Ni-reinforced Ca/α-SiAlON composites. ** The hardness of composite tends to increase with an increase in the elastic modulus as predicted. * Experimentally measured porosity from actual samples.

Monolithic
Ca-α-SiAlON  Table 3 depicts the thermal conductivity variation as a function of Ni volume fractions within the Ca-α-SiAlON matrix which are all measured at room temperature. The thermal conductivity of the composites with 10% Ni, 20%, and 30% Ni content is somewhat higher than that of monolithic SiAlON (5.65 W/m K). However, as demonstrated, increasing the Ni level had no discernible influence on the thermal conductivity values. This is obviously due to the draining out of the Ni content and not achieving the desired Ni content, as more liquid Ni ran out of the graphite die at 1500 • C. Furthermore, the heat transfer across different phases, i.e., Ca-α-SiAlON matrix and Ni, and their interface all affect a composite's thermal conductivity, directly related to the particle size, voids, porosity, and interface resistance between the two phases [55,56]. A similar effect of porosity is observed on the effective CTE of Ni/Ca-α-SiAlON composites when the Ni content increases from 10% to 30%. The measured CTE of 10% Ni/Ca-α-SiAlON composite (3.18 × 10 −6 /K) was found to be larger than that of monolithic Ca-α-SiAlON (2.63 × 10 −6 /K), which is due to a higher intrinsic CTE of pure Ni (13.1 × 10 −6 /K). Nevertheless, the effective CTE decreased as the Ni content increased. Because porosity tends to accommodate the material's internal thermal expansion, CTE decreases with associated porosity. Considering measured thermal conductivity (5.65 W/m K) and CTE (2.63 × 10 −6 /K) of α-SiAlON together with observed porosity and Ni particle size of around 2µm, computational predictions are run and compared with the experimental values. A very close agreement can be found, as depicted in Table 3.
The Vickers indentation test was done on a monolithic Ca-α-SiAlON and sintered SiAlON composites with 10%, 20%, and 30 wt% Ni, with hardness values recorded decreasing from 21.1 GPa (for monolithic α-SiAlON) to the lowest value of 16.5 GPa in the case of 30%Ni/SiAlON composite, where the porosity is the highest. As expected, an increased porosity as the Ni loading increases tends to reduce the composites' hardness. Although the elastic modulus was not experimentally measured, a complementary trend is predicted numerically. The hardness of composite tends to increase with an increase in the elastic modulus, as predicted (Table 3). Due to the mismatch in the ability to deform a ductile Ni particle and a brittle SiAlON matrix, cracks are deflected when they reach a SiAlON/Ni interface, as can be observed in Figure 9. The presence of Ni in pores/voids in the matrix were expected to promote grain bonding/wettability of Ni with Ca-α-SiAlON matrix at the interface, resulting in increased fracture toughness values due to the decrease in hardness as observed in this study. Figure 10a compares the measured K IC values of sintered Ca-α-SiAlON composites (with 10%, 20%, and 30% volume of SiC) with that of the model predictions. It is worth noting that a SiC particle size of 3 µm was used in this model, which is in line with the sintered SiC/Ca-α-SiAlON samples. A very close agreement can be found mainly at low SiC content. To determine the impact of inclusion volume fraction and size on fracture toughness of Ca-α-SiAlON composites, SiC is used as a reinforcement in the model with various particle sizes and volume fractions, as depicted in Figure 10b. As discussed earlier in Figure 5, a similar trend of growing and subsequently lowering KIC can be observed for a specific particle size and varying volume fractions, which shows that fracture by particle-cracking can dominate the interface-debonding with increasing volume fractions.
In terms of the particle size effect on K IC , it is worth noting that decreasing particle size from 50 µm to 20 µm, 5 µm, then to 5 µm enhances K IC for any given SiC content. Furthermore, as the particle size decreases from 50 µm to 1 µm, the maxima of the curves corresponding to the highest K IC increasingly shift to the right. Peak shifting implies the embedment of smaller particles in the Ca-α-SiAlON matrix due to higher particle cracking resistance for any volume content, resulting in a more considerable K IC value due to increased crack deflection. In other words, the lower K IC values of composites embedded in the Ca-α-SiAlON matrix at larger particle sizes are due to an increased tendency for particle-cracking, as also reported by Evans [57]. Due to the large surface area of smaller particles contained in the Ca-α-SiAlON matrix, crack deflections through the interface between the matrix and particle tend to promote, resulting in increased K IC FESEM image shown in Figure 11 depicts the crack length and its propagation within 10% SiC/Ca-α-SiAlON composite complementing the aforementioned reasons for enhanced fracture toughness. the matrix were expected to promote grain bonding/wettability of Ni with Ca-α-SiAlON matrix at the interface, resulting in increased fracture toughness values due to the decrease in hardness as observed in this study. Figure 10a compares the measured KIC values of sintered Ca-α-SiAlON composites (with 10%, 20%, and 30% volume of SiC) with that of the model predictions. It is worth noting that a SiC particle size of 3 µm was used in this model, which is in line with the sintered SiC/Ca-α-SiAlON samples. A very close agreement can be found mainly at low SiC content. To determine the impact of inclusion volume fraction and size on fracture toughness of Ca-α-SiAlON composites, SiC is used as a reinforcement in the model with various particle sizes and volume fractions, as depicted in Figure 10b. As discussed earlier in Figure 5, a similar trend of growing and subsequently lowering KIC can be observed for a specific particle size and varying volume fractions, which shows that fracture by particle-cracking can dominate the interface-debonding with increasing volume fractions.
In terms of the particle size effect on KIC, it is worth noting that decreasing particle size from 50 µm to 20 µm, 5 µm, then to 5 µm enhances KIC for any given SiC content. Furthermore, as the particle size decreases from 50 µm to 1 µm, the maxima of the curves corresponding to the highest KIC increasingly shift to the right. Peak shifting implies the embedment of smaller particles in the Ca-α-SiAlON matrix due to higher particle cracking resistance for any volume content, resulting in a more considerable KIC value due to increased crack deflection. In other words, the lower KIC values of composites embedded in the Ca-α-SiAlON matrix at larger particle sizes are due to an increased tendency for particle-cracking, as also reported by Evans [57]. Due to the large surface area of smaller particles contained in the Ca-α-SiAlON matrix, crack deflections through the interface between the matrix and particle tend to promote, resulting in increased KIC FESEM image shown in Figure 11 depicts the crack length and its propagation within 10% SiC/Ca-α-SiAlON composite complementing the aforementioned reasons for enhanced fracture toughness.

Conclusions
The present study uses a material-by-design approach to develop ceramic and metal particles-reinforced Ca-α-SiAlON composites with thermal and structural properties tailored for the cutting tools application. The mean-field homogenization and effective medium theories are used to predict the effective thermal conductivity computationally, coefficient of thermal expansion, elastic modulus, and fracture toughness of Ca-α-SiAlON composites by considering various reinforcement materials, volume fractions, and particle size. The simulations were run to identify the threshold inclusion size and volume fractions that could be used to achieve the desired level of improvements in the composite properties, which might lead to enhancement in cutting performance. The computational predictions are compared with the experimental measurements conducted on the Ni/Caα-SiAlON and SiC/Ca-α-SiAlON composite samples. The comparison between the model predictions and the experimental data shows good agreement with respect to the thermal and structural properties of the composite.
The overall thermal conductivity of Ca-α-SiAlON ceramic improved with the addition of the selected reinforcement materials (Diamond, SiC, Co, Cr, Ni, cBN, and TiN), depending on the inclusion volume fractions and particle sizes. A particle size of 10 µm could be sufficient to enhance the thermal conductivity of Ca-α-SiAlON to the threshold level. The most significant improvement was attained for the 15% volume fractions, with the most significant improvement noted for the diamond particles (i.e., about 43%). At the same time, SiC, Co, Cr, and Ni resulted in higher values of the effective thermal conductivities within a similar range (i.e., about 33-40%). The coefficient of thermal expansion (CTE) of Ca-α-SiAlON with the volume fractions increases for the Ni, Co, and TiN reinforcements. However, the build-up of high thermal stresses during cutting operations might restrict the use of the materials mentioned earlier as suitable reinforcers for the Caα-SiAlON composite. In this regard, SiC, cBN, and Cr become the most appropriate reinforcers for the ceramic material. The effective fracture toughness of the Ca-α-SiAlON composite also increases with increasing reinforcement volume fractions. However, for the low fracture-resistant materials, such as diamond and cBN, the improvement in fracture toughness became poor, beyond about 20% of the volume fractions. Consequently, WC, TiN, SiC, and ZrB2 particles yielded the most remarkable improvement in fracture toughness with a continuously increasing trend of fracture resistance. SiC and WC can result in

Conclusions
The present study uses a material-by-design approach to develop ceramic and metal particles-reinforced Ca-α-SiAlON composites with thermal and structural properties tailored for the cutting tools application. The mean-field homogenization and effective medium theories are used to predict the effective thermal conductivity computationally, coefficient of thermal expansion, elastic modulus, and fracture toughness of Ca-α-SiAlON composites by considering various reinforcement materials, volume fractions, and particle size. The simulations were run to identify the threshold inclusion size and volume fractions that could be used to achieve the desired level of improvements in the composite properties, which might lead to enhancement in cutting performance. The computational predictions are compared with the experimental measurements conducted on the Ni/Ca-α-SiAlON and SiC/Ca-α-SiAlON composite samples. The comparison between the model predictions and the experimental data shows good agreement with respect to the thermal and structural properties of the composite.
The overall thermal conductivity of Ca-α-SiAlON ceramic improved with the addition of the selected reinforcement materials (Diamond, SiC, Co, Cr, Ni, cBN, and TiN), depending on the inclusion volume fractions and particle sizes. A particle size of 10 µm could be sufficient to enhance the thermal conductivity of Ca-α-SiAlON to the threshold level. The most significant improvement was attained for the 15% volume fractions, with the most significant improvement noted for the diamond particles (i.e., about 43%). At the same time, SiC, Co, Cr, and Ni resulted in higher values of the effective thermal conductivities within a similar range (i.e., about 33-40%). The coefficient of thermal expansion (CTE) of Ca-α-SiAlON with the volume fractions increases for the Ni, Co, and TiN reinforcements. However, the build-up of high thermal stresses during cutting operations might restrict the use of the materials mentioned earlier as suitable reinforcers for the Ca-α-SiAlON composite. In this regard, SiC, cBN, and Cr become the most appropriate reinforcers for the ceramic material. The effective fracture toughness of the Ca-α-SiAlON composite also increases with increasing reinforcement volume fractions. However, for the low fractureresistant materials, such as diamond and cBN, the improvement in fracture toughness became poor, beyond about 20% of the volume fractions. Consequently, WC, TiN, SiC, and ZrB2 particles yielded the most remarkable improvement in fracture toughness with a continuously increasing trend of fracture resistance. SiC and WC can result in up to about 30% improvement in the overall fracture toughness. Therefore, the suitable reinforcement materials for α-SiAlON-based cutting tool inserts can be WC, SiC, and Cr due to the resulting improvement in the thermal and structural properties of the ceramic composites.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.