Advanced Computational Analysis of Cobalt-Based Superalloys through Crystal Plasticity

This study introduces an advanced computational method aimed at accelerating continuum-scale processes using crystal plasticity approaches to predict mechanical responses in cobalt-based superalloys. The framework integrates two levels, namely, sub-grain and homogenized, at the meso-scale through crystal plasticity finite element (CPFE) platforms. The model is applicable across a temperature range from room temperature up to 900 °C, accommodating various dislocation mechanisms in the microstructure. The sub-grain level explicitly incorporates precipitates and employs a dislocation density-based constitutive model that is size-dependent. In contrast, the homogenized level utilizes an activation energy-based constitutive model, implicitly representing the γ′ phase for efficiency in computations. This level considers the effects of composition and morphology on mechanical properties, demonstrating the potential for cobalt-based superalloys to rival nickel-based superalloys. The study aims to investigate the impacts of elements including tungsten, tantalum, titanium, and chromium through the homogenized constitutive model. The model accounts for the locking mechanism to address the cross-slip of screw dislocations at lower temperatures as well as the glide and climb mechanism to simulate diffusions at higher temperatures. The model’s validity is established across diverse compositions and morphologies, as well as various temperatures, through comparison with experimental data. This advanced computational framework not only enables accurate predictions of mechanical responses in cobalt-based superalloys across a wide temperature range, but also provides valuable insights into the design and optimization of these materials for high-temperature applications.


Introduction
For over half a century, nickel-based superalloys have stood as the preferred materials in the aerospace industry, particularly in the high-temperature zones of turbine engine components, such as blades and disks [1].The exceptional mechanical properties, including remarkable yield stress and creep resistance at elevated temperatures, can be attributed to the distinctive γ-γ ′ microstructure characteristics [2].Presently, nickel-based superalloys find application not only in aircraft engines, but also in power generation systems, where operating temperatures can reach up to approximately 1100 • C. The continuous evolution of nickel-based superalloys has been reached in compositions closely approaching their melting temperatures, prompting the need for the development of new high-temperature alloys [3].The recent identification of a stable ternary system featuring an L12 phase in the Co-Al-W ternary by Sato et al. [4], capable of coexisting with an FCC solid-solution cobaltbased system, has opened new avenues for creating higher-temperature resistant materials in gas turbine applications.The γ-matrix phase is predominantly an alloy of Co having a face-centered cubic (FCC) lattice structure, while the precipitates or γ ′ phase is a coherent, ordered inter-metallic of an L12 crystal structure.The revelation of γ-γ ′ microstructure compositions akin to those in nickel-based superalloys [5] has sparked considerable research efforts [3,5] towards establishing a novel class of high-temperature, high-strength materials.
Considerable investments have been made to boost the efficiency of turbine engines, resulting in substantial demands for materials employed in their hot sections.The prevailing trend leans towards materials with higher melting points and extended service lifetimes, necessitating enhanced mechanical properties, particularly in terms of creep response.This is especially crucial in turbine blades, where temperatures approach approximately 90% of the material's melting temperature [6].In response to these advancements, computational approaches have been devised to accommodate these alterations through minor compositional adjustments.These approaches aim to explore the impacts of various elements such as W, Ta, Ti, Cr, etc., each with different concentrations, where empirical testing is conducted to assess their effects on crucial mechanical properties, including yield stress, plastic deformations, and creep properties.This is predominantly attributed to the intricate dislocation mechanisms and stress states prevalent in the hot sections, specifically in turbine engine components like blades and disks.Addressing these complexities necessitates the utilization of precise, physics-based models coupled with detailed constitutive models.This combination is essential for the development of robust computational methods, facilitating efficient analysis and design within a shorter time frame.These models must strike a balance between appropriate fidelity and practicality in terms of computational time and resources, given the intricacies involved.The fundamental dislocation mechanisms are heavily dependent on the microstructure's morphology, where precipitates serve as hindrances to the continuous movement of dislocations.The dislocations either loop around or shear the precipitates, and the specific outcomes depend on factors such as temperature and stress level.
Creep tests on single crystals of cobalt-based superalloys have been conducted, with a focus on both cobalt (Co)-based and cobalt-nickel (CoNi)-based alloys [7,8].The creep resistance observed in these alloys is comparable to that of single crystals belonging to the first-generation nickel-based superalloys.These first-generation alloys are characterized by a substantial volume fraction of 70% for cubic precipitates.While there have been significant studies on nickel-based superalloys, both experimentally and computationally, at low and high temperatures, and dislocation mechanisms, experimental data on hightemperature deformation mechanisms in cobalt-based superalloys are limited, and the corresponding computational methods are even less investigated.
Comprehensive analyses of the deformation behavior of nickel-based superalloys under various loading and temperature conditions have been conducted, especially for single crystals [9,10] and polycrystals [11][12][13].Cobalt-based superalloys have been reported to exhibit similar deformation processes [14].When the temperature varies, the dislocation processes in both scenarios changes.Octahedral slip systems are claimed to be activated in both phases at lower temperatures, and mechanical property anomalies are seen as the temperature increases.
In the current study, a comprehensive computational method using a crystal plasticity finite element (CPFE) approach is developed to explore the potential mechanical properties of cobalt-based superalloys.In doing so, several elements of Ta, Ti, and Cr are added to the ternary system of Co-Al-W to investigate the effects on the 0.2% flow stress as well as creep deformations.The computational framework needs to balance fidelity and performance at the same time, making multi-scale approaches a suitable choice.Multi-scale modeling techniques are effective computational tools whenever a problem of interest cannot be modeled on a single scale.The two-phase cobalt-based superalloys are one example where the average size of precipitates is much smaller than the microstructure's dimensions.It is, as of now, computationally almost impossible to simulate the microstructure in the presence of all precipitates.Hence, multi-scale modeling techniques should be utilized in order to bridge the scales [15][16][17][18][19]. Homogenization theory, one of the techniques used in multi-scale modeling, is employed to develop constitutive models, including parametric forms of microstructural variables and their evolution [20,21].The parametric forms and their evolution can accurately represent the effects of microstructural morphology and deformation mechanisms [22].Meso-scale analyses of superalloys, incorporating precipitate distributions as well as grain structures, have been conducted using phenomenological viscoplastic constitutive laws in [23][24][25].Hardening parameters in many of the constitutive models have been expressed as assumed functions of the average precipitate size.
The proposed multi-scale scheme is in the process of being implemented into the Object-Oriented Finite Element platform, OOF, developed by the scientists at the National Institute of Standards and Technology (NIST) [26], a publicly accessible platform intended to assist materials scientists and engineers in undertaking computational investigations of structure-property relations in a large variety of systems, including mechanical systems, whose behavior is dominated by crystal plasticity.

Multi-Scale Framework
The multi-scale method offers a substantial efficiency advantage, especially in simulating polycrystals, as there is no longer a need to consider explicit representation of the precipitates in the microstructure.Figure 1 illustrates the general framework of the multi-scale approach, encompassing the sub-grain scale with an explicit representation of the morphology of the γ ′ phase, the homogenized single crystal scale with parametric constitutive equations, and the polycrystalline grain scale comprising various grains.The sub-grain scale involves detailed dislocation mechanisms dependent on critical microstructural and mechanical features in plasticity models.Microstructural features encompass the morphology and crystallography of materials, while mechanical features include load intensity, load rate, and temperature.In the realm of single crystals of superalloys, two distinct length scales can be delineated.The smaller scale, known as the sub-grain scale, involves an explicit representation of the γ-γ ′ phases, where the larger scale, or the homogenized grain scale, has implicit effects of morphology through parametric constitutive equations [27,28].Two constitutive models are proposed for each scale, tailored to address dislocation mechanisms occurring either below or above the peak temperature, typically around 850 • C, where the maximum flow stress is observed.These models correspond to distinct dislocation mechanisms operative within their respective temperature ranges.At the sub-grain scale, the constitutive models are physics-based or dislocation density-based.For temperatures below the peak temperature, the constitutive model incorporates non-Schmid components to account for cross-slip mechanics.Conversely, for temperatures exceeding the peak temperature, the model employs the glide-climb creep model [29].At the homogenized grain scale, the constitutive models take into account the implicit effects of both morphology and composition as constitutive parameters.The proposed constitutive model at this scale demonstrates the capability of encompassing temperature variations ranging from room temperature to 900 • C for various crystal orientations.Crystal plasticity models are employed in a hierarchical manner to integrate information at each scale, facilitating the implementation of the proposed constitutive models.This approach proves valuable for establishing microstructure-property relations and guiding microstructure design.This study focuses on developing a multi-scale model for cobalt-based superalloys to predict mechanical responses, utilizing a homogenized grain-scale crystal plasticity model suitable for design purposes.The method offers a significant efficiency advantage, wherein parametric hardening variables related to slip system resistances are expressed in terms of the γ-γ ′ morphology and composition.The overall model framework spans from the sub-grain γ-γ ′ level to the homogenized polycrystalline scale, which includes homogenized grains with an implicit presence of microstructural features.

Materials and Methods
The mechanical properties involving plastic deformation have been studied for both γ ′ precipitates [30] and two-phase γ-γ ′ [14,31].Single crystals of ternary Co-Al-W-based alloys exhibit poorer mechanical properties compared to single crystals of nickel-based superalloys.However, the addition of elements such as Ta and Ti shows promise in enhancing high-temperature strength [14,31].The mechanical properties of the Ni 3 Al compound in nickel-based superalloys exhibit anomalous behaviors from room temperature up to the peak temperature around 750 • C.This leads to a substantial increase in flow stress.In contrast, this region in γ ′ in cobalt-based superalloys is confined to a narrow temperature range between 650 • C to 850 • C, as illustrated in Figure 2. Furthermore, this has minimal effects on the flow stress for ternary Co-Al-W.

Constitutive Models for the Sub-Grain Level
The constitutive models at the sub-grain level are formulated using large deformation crystal plasticity.This involves a multiplicative decomposition of F = F e F p where the total deformation gradient F is composed of an inelastic, incompressible part F p , associated with slip without rotation, and an elastic part F e , which accommodates rigid-body rotations and elastic stretching.For plastic velocity gradient L p , the plastic shear strain rate γα on the slip system α (including the slip direction m α 0 and slip plane normal n α 0 in the reference configuration) and the Schmid tensor s α 0 = m α 0 ⊗ n α 0 can be employed to calculate the evolution of plastic deformation as The stress-strain relation invokes the second Piola-Kirchoff stress S and its work conjugate, Green-Lagrange strain tensor E e , in the intermediate configuration where C is a fourth order anisotropic elasticity tensor, σ the Cauchy stress tensor, and I the identity tensor.
The constitutive model utilized at the sub-grain scale draws inspiration from [13], wherein the plastic shear strain rate on a slip system adheres to the Orowan equation, γα = ρ α m b α v α , which is contingent upon the mobile dislocation density, ρ α m , and the Burgers vector, b α .The average velocity of the dislocations, influenced by the parallel and perpendicular slip resistances, τ α pass and τ α cut , respectively, can be represented by this equation.
And with the hardening parameters as where G is the shear modulus.The hardening parameters for passing and cutting are determined by the densities of parallel, ρ P , and forest, ρ F , dislocations, which are dislocation densities projected along and orthogonal to the slip planes.Statistically stored dislocations (SSDs), which are scalar and rate-dependent, contribute to the hardening through dislocation generation and annihilation processes.Additionally, due to the presence of plastic deformation gradients, particularly at the γ-γ ′ interface where misfit is a concern, slip resistance is further influenced by geometrically necessary dislocations (GNDs) [28].GNDs are included to address this phenomenon, particularly noticeable at the γ-γ ′ boundaries.These superalloys exhibit strain gradient effects due to the dimensions of γ ′ , which are captured by GNDs through morphological parameters such as shape, average size, and volume fraction of γ ′ .A higher volume fraction for the same average size leads to more precipitates within a unit cell, resulting in a larger gradient of plastic deformation.Conversely, a larger average size for the same volume fraction increases channel width, thereby reducing the gradient of plastic deformation.As a result, the densities of parallel and forest dislocations can be expressed in terms of scalar SSDs and vectorial GNDs, influencing the mobile dislocation density ρ α m as described [27]. (5)

Constitutive Models for the Grain Level
There are two primary dislocation mechanisms observed from room temperature up to the peak temperature, which corresponds to the maximum flow stress at around 850 • C, and temperatures higher than that.Before reaching the peak temperature, the activation of octahedral slip systems is reported in both phases.As the temperature increases, anomalies in flow stress are observed, as depicted in Figure 2.These anomalies are attributed to dislocation behaviors in the intermetallic γ ′ phase, where screw dislocations tend to become locked in a Kear-Wilsdorf (KW) configuration, as illustrated in Figure 3, due to cross-slip [9].The rate of lock formation rapidly increases as the temperature rises, reaching a maximum around the peak temperature.At low to mid-temperatures, particularly at larger strains, dislocations in the channel tend to shear γ ′ precipitates, and this rate increases with rising temperature.Similar to nickel-based superalloys, superdislocations have been observed instead of regular dislocations in cobalt-based superalloys [14].These superdislocations exhibit a behavior where they change planes and cross-slip along the anti-phase boundary (APB).However, at elevated temperatures, the dislocations lack sufficient energy to cut through the precipitates.Instead, they glide just inside the channel of the γ phase [14,33,34].This glide in the channel typically continues by climbing along the precipitates, as illustrated in Figure 4, giving rise to a glide-climb dislocation mechanism [27,29].Accordingly, two constitutive models are employed to address the distinct dislocation mechanisms.The first model accommodates the cross-slip mechanism, primarily occurring at low to mid-temperatures.The second model addresses the glide-climb mechanism, predominant at higher temperatures.

Constitutive Model for Cross-Slip Mechanism
The grain level, or homogenized single crystal grain-scale, for cobalt-based superalloys is obtained from a parametrized representation of the morphology of the precipitates, following the method explained in [35].This model is a function of the microstructure's morphology, especially the precipitates, containing their shape, size, and volume fraction of the γ ′ phase.The homogenized constitutive model is an Orowan-based model utilizing hardening parameters of thermal shear resistance, τ α * , and athermal or passing shear resistance, τ α pass , which reflect the effects of forest and parallel dislocation densities, as employed in the sub-grain scale, stated as Dislocations become active when the effective shear stress, represented by the difference between the resolved shear stress, τ α , and the passing stress, τ α pass , are expressed as τ α e f f = |τ α − τ α pass |, which is greater than zero.This implies that the resolved shear stress on a slip system must exceed the passing stress for dislocation motions to occur.In general, τ α * represents the cutting stress, τ α * = τ α cut , as a morphology-dependent parameter in the homogenized model [13], though additional terms are needed to be modified for the crossslip mechanism and also to account for the composition effects.The reference shear strain rate, γα * , is a function of plastic strain and morphological parameters [35].Both cutting and passing shear resistances undergo evolution due to self and latent hardening, and this evolution is defined as The total slip resistance rate can be written as ṡα = ( τα pass ) 2 + ( τα * ) 2 .The hardening matrix, h αβ , can be calculated using the interaction matrix, which accounts for interactions between slip systems by The resistance parameter, h β , is a function of the saturation shear resistance, s β sat , and total slip system resistance.h 0 and r are material constants, and q αβ = q + (1 − q)δ αβ is the interaction coefficient matrix, which includes q as a latent-hardening factor taken to be 1.4.

Modification of τ α
* for KW Lock Formation Due to Cross-Slip Mechanism As mentioned in the previous section, an additional term is required to integrate into the thermal shear resistance, τ α * , to accommodate cross-slip resistance resulting in τ α * = τ α cut + τ α cross .Locked cross-slipped dislocations, as depicted in Figure 3, increase the hardening by acting as obstacles to further dislocation motions.The enthalpy term for the cross-slip mechanism, developed in [13] for nickel-based superalloys, is used and modified for cobalt-based superalloys.In this modification, the non-Schmid components of τ α pe , τ α se , and τ α cb are responsible for dislocation dissociation during the cross-slip mechanism.However, there are no dislocation core effects for cube slip systems.Therefore, the cross-slip resistances are considered only for the octahedral slip systems as a function of non-Schmid factors and temperatures.For cube slip systems, it will be solely a function of temperature.
on cube slip systems (9) Following [13], the cross-slip shear resistance can be stated as where The temperature parameter ξ = ξ 0 exp( A θ−θ c ) signifies the strength of the pinning obstacles and diminishes as temperature rises.The rate of decrease for this parameter accelerates notably as the upper boundary temperature θ c is approached, at which climb mechanism ensues.Initially, this parameter has a value exceeding unity for temperatures below the peak temperature, but it decreases to less than unity for temperatures beyond that.The yield stresses corresponding to three temperatures of 300 K, 900 K, and 1200 K are utilized to calibrate the three parameters ξ 0 , A, and θ c , which results in 5.5, 50, and 880, respectively.G is the shear modulus, and ρ 0 is the initial cross-slip dislocation density set to be 5.0 * 10 15 .The enthalpy term has material constants k 1 and k 2 calibrated as 0.2, and 0.7, Γ 111 , and Γ 010 (J/m 2 ) are the APB energies per unit area on the octahedral and cube planes, respectively.t α xx is the normalized non-Schmid component of the resolved shear stress as t α xx = τ α xx Γ 111 /b , B = Gb 2 2πΓ 111 , and c H = Gb 3 4π .The cross-slip shear stress increases as the dislocation densities of the cross-slip mechanism rise.Conversely, the obstacle strength diminishes with rising temperature.Therefore, there is a competition between the increasing strength with the formation of KW locks and the decreasing obstacle strength with increasing temperature.

Homogenized Constitutive Model from the Sub-Grain Representative Volume Element (RVE) Model
The constitutive parameters described in Equation (6) for the homogenized single crystal grain-scale model are influenced by macro-micro energy equivalence.In this framework, the micromechanical analysis is carried out using the sub-grain RVE model.Functional parameters within the constitutive model are expressed in relation to morphological variables.
The sub-grain microstructural RVE comprises a two-phase γ-γ ′ system, characterized by three morphology variables to account for the volume fraction, shape, and size of the γ ′ phase.These variables are designated as follows: (1) The constitutive parameters at the single crystal grain scale are denoted as τ α cut0 (V p , n, w), τ α cross (V p , n, w), and s α sat (V p , n, w), as expressed in Equation ( 6) by comparing the stress-strain response obtained from the homogenized constitutive model with the volumeaveraged response from simulations of the sub-grain RVE model.This calibration process involves generating various RVE microstructures with different morphologies, including distinct volume fractions (V p ), shape factor (n), and channel width (w) [35].Subsequently, dislocation density-based models at the sub-grain scale are applied to each of these microstructures to simulate volumetric stress-strain curves up to 10% true strain.To adhere to the Hill-Mandel principle, the same volumetric stress-strain response is replicated at the grain scale using the homogenized constitutive model by adjusting four constitutive parameters.This comprehensive set of simulations yields the functional forms of the single crystal constitutive parameters through the least square minimization method.The final form of these constitutive parameters can be found in [28].Ultimately, the homogenized single crystal grain-scale model with parametric constitutive relations is validated against findings from the sub-grain scale and experimental data, as documented in [13].This validation demonstrates a significant acceleration in the processing speed by approximately four orders of magnitude, all while maintaining a high level of accuracy.

Modification of τ α
* for the Effect of Composition Experimental data indicate that the addition of certain elements, such as Ta and Ti [7,8,14], has positive effects on mechanical properties, including yield stress and creep resistance.Yet, the inclusion of materials like Cr can have negative effects on these mechanical outcomes.In other words, the addition of these materials alters the slip system resistances.Generally, they can be incorporated into the total slip system resistance as where x m is the atomic concentration of element m and τ α m is the temperature-dependent shear resistance to account for compositions.Now, the total slip system resistance τ α * can be stated as

Constitutive Model for Climb and Glide Mechanism
At elevated temperatures and under lower applied stress, dislocations are unable to cut through the precipitates.Instead, they climb along the interface of γ-γ ′ towards the matrix, where glide takes place in the channel before reaching the next precipitates [29], as illustrated in Figure 4b.
The rate of plastic shear strain for a given slip system can be calculated in terms of the morphology of γ-γ ′ and the composition of γ-γ ′ using the Orowan-based constitutive model, as outlined in [27] by V p , d, and w represent the γ-γ ′ morphology volume fraction, average size, and average distance between the precipitates, respectively, while Γ α e is the released frequency containing the effects of composition and ρ α m is the mobile dislocation density.Utilizing Gibbs free energy [36], the released frequency can be rewritten as where the effective diffusion coefficient due to the climb mechanism is Hence, the crystal plasticity constitutive model for the climb and glide mechanism can be stated with The effective diffusivity in Equation ( 16) is considered by the Arrhenius representa- ), where D α e f f 0 and Q α e f f are the effective pre-exponential coefficient and activation energy, respectively.
The effective activation energy is assessed in terms of the activation energy [29] of the constituents in the composition.This assumes that the impact of solute elements on the solvent alters the effective activation energy in pure cobalt, considering an average weighted percentage of the various solutes by where Q α Co is the activation energy for self-diffusion of cobalt and Q α m,Co the activation energy for the solute with the atomic concentration of m.The effective pre-exponential coefficient can be calculated in terms of different elements in the composition and the weighted percentage; , where D α 0m,Co is the pre-exponential factor of solute concentration in cobalt.

Validation of the Computational Model with Experiments
The computational model at the grain scale for cobalt-based superalloys incorporates two dislocation mechanisms over a temperature range from 25 • C to 1000 • C. At lower temperatures up to 800 • C, the yield stress experiences almost a linear decline from room temperature to around 600 • C, followed by an anomaly between 600 • C to 800 • C.This anomaly is primarily attributed to the cross-slip mechanism following KW locking formations.At temperatures higher than 800 • C, a sharp decline is observed, indicating that dislocations do not have enough energy to cut through the precipitates; instead, they glide inside the channel and climb along the precipitates.To validate both constitutive models, yield stresses from room temperature up to 1000 • C are computed and compared to the experimental data.Hence, simulations were conducted for three ternaries with the nominal chemical compositions (mole fraction × 100) of Co-9.2Al-9W, referred to as 9 W; Co-9.4Al-10.7 W, referred to as 11 W; and Co-8.8Al-9.8W-2Ta,referred to as 2 Ta.The microstructures after heat treatments are shown in Figure 5.According to [14], the microstructures were homogenized at temperatures between their γ ′ solvus (around 1300 K) and solidus (around 1700 K) temperatures, and then aged at temperatures about 60 K to 130 K below the γ ′ solvus temperatures.Subsequently, the simulation results were compared with the experimental data reported in [14].The volume fraction in 9 W is 60%, which is increased by adding tungsten content in 11 W to 80% and in 2 Ta to 70%.The variation of 0.2% flow stress vs. temperature is illustrated in Figure 6.All tests depicted in Figure 6 entail tensile tests conducted under a constant strain rate of 10 −4 1/s.The volume fraction of precipitates varies: 60% for Co-9Al-9W, 80% for Co-9Al-11W, and 70% for Co-9Al-9W-2Ta.According to the data, the changes in flow stress due to the addition of tungsten content from 9 W to 11 W are negligible.However, adding 2% Ta results in an increase in flow stress of almost 20% for temperatures below 800 • C. The deviation between experimental and computational data at very high temperatures may be related to the rafting mechanism, which is not considered in this study.

Composition Effects on Flow Stress
Experimental data indicate that certain elements, such as tungsten, tantalum, and titanium, have positive effects on mechanical properties, including flow stress and creep resistances, whereas elements like chromium exhibit negative effects on these properties.The impact of these elements depends on the testing temperature, with the positive or negative effects diminishing as the temperature increases.These effects can be expressed in Equation (11) as a modification to the slip system resistances.In this equation, τ α m will have a general form, particularly under constant strain rate by τ α m0 is the initial shear resistance, which is positive for Ta and Ti and negative for Cr.T re f is the reference temperature, which is calibrated to be 1100 Kelvin for all three elements.The equation is valid for the temperatures under T re f .The initial shear resistance along with the reference temperature is provided in Table 1.To investigate the effects of elements such as tungsten, tantalum, and titanium on flow stress, the microstructures provided in [37] are considered.The base microstructure has a nominal chemical composition (mole fraction × 100) of Co-10Ni-(9 − x)Al-(9 − x)W-2xTi.The flow stress at different temperatures is similar to the ternary case when x is 0. The microstructures for x = 1 to 4 are shown in Figure 7 as 2 Ti, 4 Ti, 6 Ti, and 8 Ti.According to [37], the volume fraction of the precipitates increases from 56% to 68% when the content of titanium in mole fraction × 100 rises from 0 to 8%.The variation of flow stresses with respect to temperature for these five microstructures is depicted in Figure 8.It is evident that adding titanium, even just 2%, increases the yield stress by more than 20% at the peak temperature.All microstructures with titanium content exhibit an increase in flow stress at temperatures below 900 • C. The flow stress for 8 Ti shows competitive values with nickel-based superalloys.

Validation of the Creep Model with Experiments
At elevated temperatures, gliding inside the channel and climbing along precipitates will be the dominant mechanism.The plastic shear strain or creep properties can be measured using Equation (16).The first study will be on the creep behaviors of ternary cobalt-based superalloys and the effect of tungsten content.As learned from Figure 6, tungsten will have small effects on flow stress, which is the onset of plastic strain.To understand the effect of W on the creep properties, two ternary microstructures are taken from [34], as shown in Figure 9, with the composition of Co-9Al-9W and Co-9Al-11W having 41% and 73% of the γ ′ volume fraction, respectively.
Creep tests on each microstructure are conducted, with the first test under 275 MPa at 900 • C and the second one under 460 MPa at 850 • C. The experimental data are compared with the computational results in Figure 10.As noted, adding more tungsten increases the creep resistance in cobalt-based superalloys, which can be explained by two reasons.First, tungsten is a slow diffuser, and second, adding more W increases the volume fraction of γ ′ , which is the main morphology factor strengthening creep properties in cobaltbased superalloys.To investigate the effect of volume fraction on creep responses, three simulations are performed for tests with the applied stress of 275 MPa at 900 • C based on the constitutive model stated in Equation ( 16).The volume fraction of γ ′ in ternary cobalt-based superalloys can be changed by altering the content of tungsten.The change in volume fraction of γ ′ from 40% to 80% will significantly decrease the plastic strain, as shown in Figure 11.The plastic strain in the microstructure with 40% of the γ ′ volume fraction will be more than 10% in less than 200 h, while the plastic strain in the microstructure with 80% of the γ ′ volume fraction will show less than 2% plastic strain.

Effects of Composition on Creep Properties in Cobalt-Based Superalloys
Composition has significant effects on the creep properties of cobalt-based superalloys.In this section, we inspect the effects of Ti, Ta, and Cr and combinations of these elements in the composition on creep properties.The first composition is Co-9.5Al-5W-16Ni-6Cr-xTa, reported by [38].The microstructures with 1.8 Ta and 2.8 Ta are taken from [38] and shown in Figure 12.As reported in [38], increasing Ta from 1.8% to 2.8% has slight effects on the volume fraction of γ ′ ; however, it results in coarser precipitates, where the average size increases from 241 µm to 301 µm.The variation of the plastic deformation over time for these two microstructures under 275 MPa at 900 • C is shown in Figure 13.The results indicate that increasing Ta by 1% improves creep resistance notably, with the plastic strain for this composition being around 2% after 200 h, while for 1.8 Ta, it will be more than 10%.The simulation results and the experimental data have the same trend; however, there are some differences between them, which are possibly due to the high content of nickel in the compositions, which are not investigated in the current research.The second example from [7] provides better insights into how these three elements in the composition will change the creep behaviors.The base microstructure is ternary cobalt-based superalloys Co-9Al-11W with a volume fraction of 56%, as shown in Figure 14.As shown in Figure 15, the plastic deformation for the ternary will be around 8% after 200 h.By adding 2% Ta to the ternary, the volume fraction will increase to 70%, as illustrated in Figure 14b, where the plastic strain will be around 1% after 200 h and less than 8% after 400 h for this microstructure, which exhibits the positive effects of adding tantalum.Now, if we add 4.5% Cr to the previous microstructure, the volume fraction will decrease to 43%, as demonstrated in Figure 14c; however, it diminishes the creep resistances dramatically, with the plastic deformation being more than 10% after 200 h.In fact, the creep strength for this microstructure is less than the ternary one.A total of 6% Ti is added to the ternary for the last tested microstructure, as shown in Figure 14d, to increase the volume fraction to 70%.The remarkable increase in creep resistance can be seen in Figure 15, where the plastic strain is less than 3% after 400 h.Based on the current experimental data obtained from the high content titanium samples, it seems that the simulation is providing lower predictions.With access to more experimental data, it would be possible to better evaluate the accuracy of the computational approach's convergence.

Conclusions
A comprehensive computational framework is proposed to predict the mechanical properties of cobalt-based superalloys.The framework employs a multi-scale approach to capture the effects of both morphology and composition.Applicable from room temperature to 900 • C, the model accounts for two dominant dislocation mechanisms: at lower temperatures, precipitates are cut under sufficient loading, while at higher temperatures, they glide inside the channel and climb around precipitates.The homogenized constitutive models for simulating both mechanisms consider the effects of morphology and composition.The mechanical properties, such as flow stress and creep strength, in ternary cobalt-based superalloys can be notably improved by adding certain elements to the composition, making them competitive with nickel-based superalloys.Simulations indicate that increasing tungsten content significantly enhances the volume fraction of γ ′ and creep resistances, with minimal effects on flow stress.Additionally, adding tantalum or titanium, even in small amounts, has substantial increasing effects on both flow stress and creep resistance.The weakening effects of chromium on both flow stress and creep strength are also noted.Additionally, the model holds promise for simulating material behavior under extreme conditions, offering insights into its potential applications across various industrial settings.Exploring these avenues further could provide valuable insights into the practical utility of our model.
As part of future work, the constitutive models will be integrated into the Object-Oriented Finite Element (OOF) Framework developed by scientists at the National Institute of Standards and Technology (NIST).This publicly accessible platform aims to assist materials scientists and engineers in conducting computational investigations of structure-property relations, particularly in systems dominated by crystal plasticity in mechanical behavior.Funding: There was no funding for this project however we acknowledge Material Measurement Laboratory of National Institute of Standards and Technology.

Figure 1 .
Figure 1.Multi-scale framework from single crystal to polycrystalline microstructures of superalloys: (a) Turbine blades with single crystalline microstructure and cubic precipitates.(b) Discretized finite element mesh of a single crystal description with a two-phase material.(c) Homogenized single crystal grain-scale finite element model.(d) Turbine disk with polycrystalline mictostructure.

1 ;
for volume fraction of the γ ′ phase; (2) n, for the shape factor of the γ ′ phase, which is obtained by considering γ ′ as a super-ellipsoid with the same dimension in principle directions a = b = c in the equation (3) w, the distance between precipitates.

Figure 6 .
Figure 6.Comparison of yield stress vs. temperature between simulations and experimental data for cobalt-based superalloys 9 W, 11 W, and 2 Ta [14].

Figure 11 .
Figure 11.Creep resistances for different volume fractions of the γ ′ phase in ternary cobalt-based superalloys.
All microstructures are subjected to creep under 310 MPa at 900 • C.

Figure 15 .
Figure 15.Creep resistances for different compositions of cobalt-based superalloys.

Table 1 .
Nominal chemical composition (mass fraction) for four generations of single-crystal nickelbased superalloys.