Quantifying Solid Solution Strengthening in Nickel-Based Superalloys via High-Throughput Experiment and Machine Learning

Zihang Li,Zexin Wang,Zi Wang,Zijun Qin,Feng Liu,Liming Tan,⋆,Xiaochao Jin,Xueling Fan and Lan Huang

1State Key Laboratory of Powder Metallurgy,Central South University,Changsha,410083,China

2AECC Commercial AircraftEngine Co.,Ltd.,Shanghai,200241,China

3State Key Laboratory for Strength and Vibration of Mechanical Structures,School of Aerospace Engineering,Xi’an,710049,China

ABSTRACT Solid solution strengthening(SSS)is one of the main contributions to the desired tensile properties of nickel-based superalloys for turbine blades and disks.The value of SSS can be calculated by using Fleischer’s and Labusch’s theories,while the model parameters are incorporated without fitting to experimental data of complex alloys.In this work,four diffusion multiples consisting of multicomponent alloys and pure Ni are prepared and characterized.The composition and microhardness of single γ phase regions in samples are used to quantify the SSS.Then,Fleischer’s and Labusch’s theories are examined based on high-throughput experiments,respectively.The fitted solid solution coefficients are obtained based on Labusch’s theory and experimental data,indicating higher accuracy.Furthermore,six machine learning algorithms are established,providing a more accurate prediction compared with traditional physical models and fitted physical models.The results show that the coupling of highthroughput experiments and machine learning has great potential in the field of performance prediction and alloy design.

KEYWORDS Multicomponent diffusion multiples;solid solution strengthening;strengthening models;machine learning

1 Introduction

Superalloys,especially Ni-based superalloys,are crucial materials used in the high temperature components,including turbine blades and turbine disks of aircraft engines [1].For example,engine turbine blades are mainly made of Ni-based single crystal superalloys due to their superior mechanical properties at elevated temperatures,such as creep resistance and high temperature strength [1–3],which can be attributed to the strengthening effects of solid solution matrixγphases and coherent precipitationsγ′phases [4–13].Regarding the tensile properties,the strengthening mechanisms of nickel-based superalloys mainly include solid solution strengthening(SSS),precipitation strengthening(PS),and grain boundary strengthening(GBS)in polycrystalline alloys[4–6,9].The quantification of these strengthening effects is significant in alloy design[14–16].SSS plays an important role in creep resistance and tensile properties[11,16,17],describing the hindrance of solute atoms to the movement of mobile dislocations attributed to elastic misfit and atomic size misfit[17].

SSS was firstly proposed in binary or ternary alloy systems[18–21].Gypen et al.[22,23]further developed a suitable model for multicomponent alloys by integrating the strength increments of various alloying elementsiin the Ni-ibinary systems.This model was extended by Roth et al.[24]and Galindo-Nava et al.[7]to estimate SSS in nickel-based superalloys based on Fleischer’s[19]and Labusch’s[20]theories,respectively.However,the SSS coefficients of traditional models are related to simple alloy systems without fitting the experimental data of multicomponent alloys[7,24].

In addition,diffusion multiple technologies as a high throughput way has been widely used in the fields of phase diagram research[25–29],diffusion kinetics[30,31],and alloy design[32,33]et al.Compared with traditional alloy diffusion couples,diffusion multiple techniques can effectively integrate composition,microstructure,and even property data in a simple sample.Ur-Rehman [34] used diffusion couples to research the SSS in binary nickel-and cobalt-based superalloys,but few works on SSS have been studied for multicomponent systems.

In this work,the SSS effects in multicomponent nickel-based superalloys based on highthroughput experiment techniques were studied,which yielded vast experimental data.However,how to analyze the data becomes a big challenge for material scientists.Recently,data-driven machine learning (ML) algorithms played a key role in the field of material sciences,including phase identifying and mechanical performance prediction [35–38].Qin et al.[39] coupled the highthroughput experiment and ML to quickly identify the precipitation phase of superalloys with superior accuracy.Wang et al.[40]accelerated the compositional design of nickel-based superalloys with target mechanical properties by integrating high-throughput experiments and artificial neural networks(ANN).Liu et al.[41]designed a high-strength and tough NiCoCrAl eutectic high entropy alloy by combining machine learning and CALPHAD method.Wen et al.[42]formulated a materials design strategy that combines ML surrogate models with experimental design algorithms to search for the high-entropy alloys (HEAs) with excellent hardness,and the high hardness HEAs were obtained by only seven iterative experiments.Liu et al.[43]developed a divide-and-conquer self-adaptive(DSCA)learning method incorporating five material descriptors to accelerate the prediction of the creep rupture life of nickel-based superalloys,which shows superior predictive performance.Hence,coupling high-throughput experiments,calculations and ML exhibit great potential in performance prediction.

Here,four multicomponent diffusion multiples were prepared to rapidly obtain the composition-SSS hardness data of nickel-based superalloys.Thereafter,the existing SSS models were evaluated by experimental data,and related model parameters were fitted.Furthermore,several ML algorithms were used to predict the SSS values and their accuracy was compared with physical models.

2 Experiment and Method

2.1 Diffusion Multiple Experiment

The base alloys were prepared by vacuum induction melting (VIM) technology,and the experimental compositions are listed in Table 1.The multicomponent diffusion multiples were assembled from base alloy ingots and a pure Ni can,as shown in Fig.1.Seven diffusion couples (DCs) were obtained in one diffusion multiple and denoted as DC1 to DC7,as highlighted by the black lines in Fig.1.Four assembled diffusion multiples were hot isostatic pressed(HIPed)at 1180°C/150 MPa for 8 h to achieve good interfacial contact among the base alloys.Thereafter,the diffusion multiples were sealed in a vacuum quartz jar to avoid oxidation of the sample surface during heat treatment process.Subsequently,the samples underwent solution heat treatment at 1180°C for 1000 h to form a continuous composition gradient.Finally,four kinds of aging treatment 700°C/2000 h,800°C/2000 h,900°C/1000 h,and 1000°C/1000 h were used on four same diffusion multiples,respectively.Considering the different in aging heat treatment,the prefix of WX is used to distinguish the DCs with same combination but different aging heat treatment,where X is the aging temperature.For example,W800 DC1 refers to the DCs aged at 800°C,and composed of pure Ni and W1 alloy.

Table 1: Chemical compositions of base alloys for preparing multicomponent diffusion multiples(wt.%)

Figure 1:Schematic diagram of multicomponent diffusion multiple

The compositions of DCs were detected using Micro-XRF (X-ray fluorescence) device Bruker M4 Tornado due to high characterization efficiency compared with Energy Dispersive Spectrometer(EDS)and Electron Probe Micro Analyzer(EPMA).Then,the microstructure was characterized by a ZEISS SUPRA 55 scanning electron microscope(SEM)equipped with an automated photographic platform.The specimens for SEM observation were prepared by polishing and etching in the solution of 33 ml H2O + 33 ml CH3CH2OH + 33 ml HNO3+ 1 ml HF at room temperature for 10–20 s.Ultimately,the hardness of different zones was measured by the automatic microhardness tester Qness Q10a,and indentation pitch was set at 100 μm to avoid the influence of adjacent indentations.

2.2 SSS Models

The yield stressσyof alloy can be calculated by

whereσbis the basic strength (in MPa),including the contributions of pure matrix and extra boundary strengthening effect,σSSis the strength increment caused by SSS.However,the yield stress of micro-zones cannot be directly measured.Thus,following empirical relationship between the Vickers hardness and yield stress was used to quantify the SSS,which can be expressed as[6]

whereκequal to 2.46 MPa/Hv,His the measured Vickers hardness(in Hv).

Since the hardness of the micro-zones at distance equaling to 0 μm has almost no the SSS effect,the basic hardnessHbof each DC can be defined as

whereH0is the hardness of the micro-zone at distance equaling to 0 μm,σSS0is the SSS value of microzone at distance equaling to 0 μm.

To evaluate the effect of the SSS,Fleisher’s theory[17,19]employed the average value of interaction force between solute atoms and mobile dislocation.Hence,the increment of critical stress caused by SSS is given by

whereCis a constant,Fis the force of the solute atoms to hinder the movement of the dislocation,bis the Burgers vector,Tis the dislocations line tension,cis the concentration of solute atoms.However,Labusch’s theory [17,20] considered the effect of average field distribution of obstacles on mobile dislocation,and the increment of critical stress is defined as

wherewis also a constant.Eqs.(4)and(5)can be simplified to

whereβiis the SSS coefficient of elementiin the Ni-ibinary systems,xiis the atomic fraction (in at.) of elementiin the alloy,qis equal to 1/2 and 2/3 for Fleisher’s theory and Labusch’s theory,respectively.To further modeling SSS in multicomponent alloys,Gypen et al.[22] proposed an equivalent concentration method,and the increment of critical stress derived from SSS is expressed as

SSS coefficientsFβi(in MPa/at.1/2) andLβi(in MPa/at.2/3) for Fleisher’s theory and Labusch’s theory,respectively,are listed in Table 2.Therefore,the increment of yield stress caused by SSS in the nickel-based superalloys can be given by

or

The fitting of SSS coefficients in Section 3.2 is implemented with leastsq function in the scipy.optimize package within Python 3.7.7 environment.

2.3 ML Models

There are many effective models in the machine learning community,and different models exhibit respective advantages and disadvantages.Therefore,several commonly used ML algorithms were employed to quantify the SSS of nickel-based superalloys,including Linear Regression(LR),Support Vector Regression (SVR),K-Nearest Neighbor (KNN),Random Forest (RF),Gradient Boosting Decision Tree(GBRT),and ANN.

LR is a simple linear model,which is widely used in the prediction of objective properties owing to it can obtain specific expressions to clearly illustrate the relationship between input features and target performance.Especially,LR model can effectively reduce the risk of overfitting by introducing a regularization term based on simple multiple linear regression.

SVR is an important nonlinear regression model in a machine learning community,which could predict the objective properties by minimizing the cost functions.Especially,different kernel functions can be used in the SVR model to map the feature variables from low dimensional space to high dimensional space,including radius basis function (RBF) kernel,linear kernel,polynomial kernel,and sigmoid kernel.Compared to LR,SVR can effectively solve nonlinear regression problems and has a stronger generalization ability.In this paper,the RBF kernel is adopted during the experiments.

K-Nearest Neighbor is a non-parametric supervised learning method first developed by Evelyn Fix and Joseph Hodges.It solves the regression problems by finding the k nearest neighbors of a sample and assigning the average value of an attribute of these neighborhoods to the sample.

RF is one of the ensemble regression models,which consists of many trees.Firstly,the bootstrap resampling technology is used to extract some samples from the original sample space to construct the decision trees.Then,each of decision trees is conducted to predict the objective properties.Finally,the predicted values were calculated by using the average value of all decision tree predictions.In addition,the RF can provide the interpretability of ML model by evaluating the importance of each input feature’s contribution to the output feature variable on each decision tree.

GBRT is another ensemble model with a regression tree as the base learner.Different from the RF,each tree of GBRT learns from the residuals of all previous trees and the negative gradient value of the loss function was used to approximate the residual.Finally,the GBRT model accumulates the predictions of all regression trees as the final prediction result.

ANN is a traditional model in the machine learning community and has been widely used to quickly construct the relationship of composition,microstructure and properties in the field of material science.Generally,the ANN consists of three major components,including input layers,hidden layers,and output layers.Input layer is commonly the feature variable;the hidden layer is used to extract and process the input features into another dimension space by the neurons;the output layer is the objective performance.

In order to improve the interpretability and generalization ability of ML models,physical features based on empirical knowledge are usually introduced as input features.As is known to all,the SSS of nickel-based superalloys originates from the elastic interactions between the local stress fields of solute atoms and mobile dislocations,including the elastic misfit and atomic size misfit[7].Hence,the atomic misfit and modulus misfit were added as the physical descriptors to the SSS dataset to enhance the performance of ML models.The atomic misfit and modulus misfit can be expressed as

wherenis the number of constituent elements,xi,ri,andGiare the atomic fraction,atomic radius(in nm),and shear modulus(in GPa)of theith element,respectively.andare the mean atomic radius and mean shear modulus,respectively.The atomic radius and shear modulus used in this work refer to reference[7].

In the following,to avoid the differences in orders of magnitude,each feature of the input data is scaled and normalized by the Min-Max normalization method,which is defined as

wherexandx′are the original value and normalized value of each input feature,respectively.

Furthermore,to evaluate the performance of these ML models,three different indicators,including root mean squared error(RMSE),mean absolute error(MAE)and R-square(R2)are proposed.Givenmobserved valuesyj,j=1,2,···,mand predicted values,j=1,2,···,m,these indicators can be calculated by the following equations:

Here a lowerRMSEorMAEand a higherR2mean a better ML model.

In this work,169 data points after selection (mentioned in Section 3.2) from W800,W900,and W1000 MCDMs were used as the train set,49 data points from W700 multicomponent diffusion multiples were used as test set.Moreover,to prevent the overfitting of ML algorithms,the grid search and ten-fold cross validation (CV) were used to tune the optimum hyperparameters of ML models.Finally,the ML models were constructed and compared with the physical models.All the ML algorithms mentioned above are implemented in the Scikit-learn package within Python 3.7.7 environment.

3 Results and Discussion

3.1 Composition-Microstructure-Microhardness

The detected compositions of DCs in MCDM W700,W800,W900,and W1000 were shown in Fig.2.To obtain smooth composition-distance profiles,the superimposed Boltzmann function was used

wherenis the number of the Boltzmann functions used,x−andx+are the terminal compositions,pjandqjare the fitted parameters.As shown in Fig.2,the fitted curves are in good agreement with experimental data.Thus,the fitted compositions were used in the following sections.

The existence of theγ′phases,the geometric close-packed (GCP) phases,and the topological close-packed(TCP)phases will interfere with the quantification of SSS.Hence,the hardness of singleγphase regions was selected to avoid the influence of other phases.Fig.3 shows the singleγphase region of W800 DC1 from 0 to 1360 μm,and the rest is aγ+γ′dual phase region.Besides,the boundary betweenγandγ+γ′region of W700 DC1-7 is not defined due to the unexpected corrosion.Therefore,the regions of W700 DC1-DC7 from 0 to 1200 μm are assumed to be a singleγphase region.

Figure 2: (Continued)

Figure 2: Comparison between the experimental (denoted in symbols) and fitted (denoted in lines)composition-distance profiles of DCs:(a)W700 DC1-DC7;(b)W800 DC1-DC7;(c)W900 DC1-DC7;(d)W1000 DC1-DC7

Figure 3:Secondary electron(SE)images of W800 DC1:(a)distance from 1600 to 3000 μm(b)γ +γ′dual phase region(c)distance equal to 1360 μm(d)single γ phase(e)distance from 0 to 1400 μm

The singleγphase regions of each diffusion couple were detected via the SEM images,as shown in Fig.4.Thereby,the microhardness values of singleγphase regions were obtained.The difference in microhardness can be attributed to the SSS caused by the difference in composition.Hence,the composition and microhardness can be used to quantify the SSS in the following sections.

Figure 4: Distribution of microhardness in single γ phase regions of DCs: (a) W700 DC1-DC7 (b)W800 DC1-DC7(c)W900 DC1-DC7(d)W1000 DC1-DC7

3.2 Quantification Based on SSS Models

As mentioned in Section 2.2,the basic strength is treated as a whole contributed by the strength of pure matrix and the effect of grain boundary strengthening.The contribution of grain bo undary strengthening to the total yield stress can be given by the Hall–Petch relationshipσGB=kHP/[7]wherekHPis Hall–Petch constant andDis mean grain size.It is assumed that the effects of grain boundary strengthening in the same DC are equal so that the basic hardness of position X=0 defined by Eq.(3) can be used for all positions in the same DC.Basic hardness of all DCs were listed in Tables 3 and 4.It can be found that nearly all DCs’basic hardness values in Table 3 are lower than these in Table 4 except for W900 DC3 and DC5,and W1000DC 3-5.This means the solid solution strength based on Fleisher’s theory is higher than that based on Labusch’s theory.

Table 3: Basic hardness value(Hv)with σSS0 based on Fleisher’s theory

Table 4: Basic hardness value(Hv)with σSS0 based on Labusch’s theory

The microhardness of MCDM W800,W900,and W1000 singleγphase regions DCs was used to estimate the traditional SSS models.In addition,14 sets of data with a relative standard deviation (RSD) of hardness larger than 0.06 were filtered out,as shown in Figs.5a–5c show the comparisons between experimental and fitted SSS values based on Fleisher’s theory and Labusch’s theory,respectively.It can be observed that the predicted SSS values agree well with experimental values.TheR2was used to compare the accuracy of prediction.As shown in Figs.5b and 5c,Labusch’s theory can better describe the SSS of nickel-based superalloys with higher accuracy.

Figure 5: (a) The relative standard deviation (RSD) of the hardness dataset.Comparisons between experimental data and hardness increment based on(b)Fleisher’s theory and(c)Labusch’s theory

To further improve the performance of the model,the data from DCs in MCDM W800/W900/W1000 is used to refit the SSS coefficients based on Labusch’s theory.It should be pointed out that the SSS coefficients of Co,Mo,and Cr were set to be the same as Table 4 due to excessive deviation,and the remaining fitting coefficients were listed in Table 5.It can be seen that W shows a stronger SSS effect,while Al,Ti,Ta,and Nb show a weaker SSS effect than the SSS coefficients in Table 4.Fig.6a shows the comparison between experimental SSS values and predicted values using the fitted SSS coefficients.Moreover,the model with fitted coefficients achieved a higher accuracy withR2of 0.847.

The data from DCs in MCDM W700 were used to evaluate the generalization ability of the fitted SSS coefficients.The comparisons between experimental and predicted values based on Labusch’s theory using original coefficients and the fitted coefficients were performed,as shown in Figs.6b and 6c.As can be seen,the model with the fitted coefficients exhibits higher accuracy compared with original coefficients.

Table 5: Fitted SSS coefficients

3.3 QUantification Based on ML Models

All the 10 features are considered to be associated with SSS.To identify the importance of the feature variables,RF is conducted on the SSS train set with the two calculated physical descriptors.The feature importance for each feature is shown in Fig.7.It is found that the descriptors,δrandδG,play the most importance role in SSS.As for the alloy composition,the atomic fraction of alloying elements Co and Ti are important factors to the SSS.This result may due to the fact that the element Ti can increase the atomic misfit and modulus misfit.

Figure 6:Comparison between experimental and predicted SSS values:(a)fitted SSS coefficients with experimental data from MCDM W800,W900,and W1000;(b)original SSS in Table 4;(c)fitted SSS coefficients with experimental data from MCDM W700

Figure 7:The ranking of feature importance of composition features by RF

To further confirm the effect of descriptors on SSS,six ML models(LR,SVR,KNN,RF,GBRT and ANN)were used to predict the SSS train set with and without descriptors,respectively.The average RMSE and MAE of 10 iterations of six ML models on the train set with and without descriptors are shown in Table 6.It can be seen that the prediction accuracy of the six ML models on the train set with two descriptors is only slightly better than that without descriptors.

Table 6:The comparison of prediction performance(with standard deviation)of six ML models on the SSS dataset with and without two microstructural descriptors

Fig.8 shows the comparison between the experimental SSS values and predicted values by six ML models.It can be seen that the predicted SSS value is in good agreement with the experimental values,indicating the ML models are adequate for predicting the SSS values of nickel-based superalloys.Compared to the fitted physical model withR2of 0.753,the ML models exhibit equivalent or higher accuracy.

To understand the influence of alloying elements on SSS,the feature importance of alloying elements,and Pearson correlation coefficients between alloying elements and SSS were calculated to quantify the influence of alloying elements on SSS,as shown in Fig.9.As can be seen from the Pearson correlation coefficients,the elements Al,Ti,Cr,Ta and Co have a strong relationship with SSS values.However,the Pearson correlation coefficients only reflect the linear relationship between alloying elements and SSS values.Therefore,the feature importance was adopted to further detect the nonlinear relationship between alloying elements and SSS values.A high feature importance value denotes that the feature plays a key role in the model,and the feature is significant in promoting the performance of the model.It can be found that the Ti,Cr,and Ta are important factors to the SSS values.

Therefore,the quick mapping between compositions and SSS values can be completed by the ML models.This demonstrates the role of ML methods in strengthening model,and the potential in the research of other mechanisms such as particle shearing(mainly in yield strength model),glidingclimbing(mainly in creep model)and so on.Moreover,the combination of various mechanism models can be the test-free method to obtain mechanical properties,which can reduce the cost in sample preparation and characterization.Hence,applying ML models in strengthening model is significant in accelerating the design of alloy for development of turbine blades and disks in aircraft engines.

Figure 8:Predictive performance of the six ML models on the SSS dataset:(a)RF(b)SVR(c)KNN(d)LR(e)GBRT(f)ANN without descriptors

3.4 Effect of Alloying Elements Co and W on the SSS of Nickel-Based Superalloys Based on ANN Model

Understanding the influence of alloying elements on SSS is significance for alloy design novel alloys with high SSS effect.Due to the Co and W are strong solid solution strengthening elements in multicomponent nickel-based superalloys.Therefore,we focus on the influence of Co and W elements on the SSS values.Based on the composition of W1 base alloys,the contents of alloying elements Co and W were adjusted to investigate the influence of alloying elements on the SSS effect.Fig.10 shows the variation of SSS values with different Co and W contents.It can be clearly observed that the SSS value increase with the Co content,and W element can also enhance the SSS effect when the Co content increase to a certain value.Hence,we could improve the content of Co to design novel solid solution strengthening superalloys with superior strength.

Figure 9:The Pearson correlation coefficients and feature importance of RF model without descriptors

Figure 10:Interaction effects of the elements Co and W on the SSS of nickel-based superalloys

In conclusion,compared with the physical model,the ML models show great potential for material property prediction.Moreover,the ML methods can be easily extended to other situations,such as the prediction of hardness,yield strength,creep life,etc.Furthermore,the ML models can be continuously enhanced by adding the new calculated or experimental data.

4 Conclusion

The composition-microhardness data of singleγphase regions in 28 DCs were high-throughput obtained,and the traditional SSS models of nickel-based superalloys were compared.The fitted SSS coefficients were fitted to quantify the SSS.Meanwhile,several ML algorithms were used to further quantify the SSS.Based on the experimental and modeling results,the following conclusions can be drawn:

1.The Labusch’s theory can better describe SSS of nickel-based superalloys compared with the Fleisher’s theory.

2.The fitted SSS coefficients were obtained by fitting experimental data in the framework of Labusch’s theory,which indicates that the accuracy is higher than original coefficients.

3.ML models exhibit a higher accuracy compared with the physical models,and Ti element plays a crucial role in predicting SSS of nickel-based superalloys.

Funding Statement:This work is supported by National Science and Technology Major Project(J2019-IV-0003-0070),the Natural Science Foundation of China(91860105,52074366),China Postdoctoral Science Foundation (2019M662799),Natural Science Foundation of Hunan Province of China(2021JJ40757),the Science and Technology Innovation Program of Hunan Province(2021RC3131),Changsha Municipal Natural Science Foundation(kq2014126),and Project Supported by State Key Laboratory of Powder Metallurgy,Central South University,Changsha,China.

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.