Phase field modeling microstructural evolution of Fe-Cr-Al systems at thermal treatment

Abstract

A phase field model to study dynamics of microstructure transformations and the evolution of defect structure during heat treatment of Fe-Cr-Al systems is developed. Statistical and kinetic properties of evolving microstructure and defect structure in alloys with different content of alloying elements and at different temperatures were studied. Point defects rearrangement during precipitation is discussed in details. Universality of statistical distributions over precipitate size is revealed for considered class of alloys.

1 Introduction

It is known that Fe-Cr-Al alloys are successfully used in many industrial applications where high temperature oxidation resistance is needed, for example, in fuel energy plants (Gulbransen and Andrew, 1959; Wukusick and Collins, 1964). Nuclear grade iron-based Fe-Cr-Al alloys have been developed as candidate materials for accident tolerant fuel cladding since the 1960s due to their corrosion resistance, radiation response, fuel-cladding chemical interaction, good mechanical properties for extended periods which do not significantly penalize the overall nuclear reactor operation [see discussion in Ref. (Yamamoto et al., 2015)].

By studying an influence of alloying elements on properties in these alloys it was found that the higher Cr content in these alloys would be beneficial for corrosion resistance in the high temperature water reactors and also high-temperature steam oxidation resistance (Stott et al., 1995). At the same time higher Cr content would result in potential embrittlement of these alloys at operating temperatures in the interval 288°C–320°C due to the formation of Cr-rich α′-phase (Mathon et al., 2003; Kobayashi and Takasugi, 2010; Ejenstam et al., 2015). The higher Al content is beneficial for high-temperature steam oxidation resistance, but it may raise the ductile-brittle transition temperature (Qu et al., 2013). It was shown that more than six at% of Al addition would be effective in reducing the stability of α′ precipitates, a further Al addition results in increasing the solution limit of Cr in the matrix (Kobayashi and Takasugi, 2010; Li et al., 2013; Ejenstam et al., 2015).

Extensive studies of ternary Fe-Cr-Al alloys were started from 1960s for use in nuclear applications for a wide range of compositions experimentally and theoretically. Essential progress was in studying material properties in such systems was achieved by using different methods of multiscale modeling scheme. Among theoretically obtained results one can issue ab initio calculations of energetic parameters for atoms relocations (Kresse and Hafner, 1993; Kresse and Hafner, 1994; Kresse and Furthmüller, 1996), kinetic Monte Carlo simulations of the long-term microstructural stability of Fe-Cr-Al alloys (Ejenstam et al., 2015), molecular dynamics studying the dependence of elastic, bulk and shear moduli on alloy composition (Dai et al., 2022), classical phase field modeling precipitation kinetics [see for example, Refs. (Chen et al., 2019; Lee et al., 2020)]. A correlative theoretical and experimental research strategy has been developed in Ref. (Chang et al., 2019) to investigate the Fe-Cr-Al system, which includes an integration of ab initio calculations, CALPHAD, and experiments. The CALPHAD approach has been successfully applied to obtain phase diagrams at thermodynamic equilibrium by exploiting the Gibbs energy functions (Dinsdale, 1991). It serves as a powerful tool for component selection, composition design and optimization of materials for specific applications.

Recently successful progress in understanding microstructure transformation in Fe-Cr-Al systems separation kinetics of the nano-size α′-phase in Fe-Cr-Al alloys was achieved. At the same time some open questions related to study microstructure evolution of long aging time and predicting the mechanical properties of these alloys remain actual. Indeed, an influence of alloy composition and aging temperature onto volume fraction and the mean precipitate size affecting the hardness and embrittlement in Fe-Cr-Al alloys is still under discussion.

In this paper we describe properties of nano-size α′-phase precipitation in such ternary systems by considering a phase decomposition with local rearrangement of point defects. We exploit CALPHAD approach and take into account interaction of atomic subsystem with subsystem of defects. The derived model includes elastic anisotropy and can be used to study alloys with different content of doping. Main efforts are concentrated onto studying an influence of alloying elements onto kinetics of phase decomposition at different temperatures and statistical properties of growing α′ precipitates. In our approach we study dynamics of precipitate size, their number density and volume fraction of precipitates at thermal treatment of solid solution for alloys with different compositions. Our study provides an insight into the physics of phase decomposition which can be generalized next in studying radiation induced precipitation and precipitates dissolution with defects clustering.

The work is organized in the following manner. In Section 2 we describe formalism used in phase field modeling. Section 3 is devoted to stability analysis and discussion of numerical modeling results. Here the main dynamical laws for precipitate size and their number density are shown and analysis of phase decomposition kinetics at different content of alloying elements and annealing temperatures is provided. We conclude in Section 4.

2 Phase field model

We describe ternary systems Fe-Cr-Al with molar concentration of Fe, Cr and Al (coinciding with atomic ones) denoted as x_μ = N_μ/N with μ = {Fe, Cr, Al}, where N_μ is the number of atoms of the sort μ, N is the total number of atoms, the mass conservation law is satisfied ∑_μx_μ = 1. Concentrations of point defects c_d with d = {i, v} denotes interstitial (i) and vacancies (v) is defined in the same manner. We consider introduced quantities as continuous fields evolving in space r and time t: x_μ = x_μ(r, t), c_d = c_d (r, t). We introduce point defects as additional substances allowing us to distinguish atomic and point defects subsystems. These equilibrium defects do not contribute to mass conservation due to their small concentrations.

2.1 Gibbs energy for ternary system with defects and elastic anisotropy

The total Gibbs functional for the class of considered systems has the form:where the term G_ch relates to the chemical Gibbs energy depending on the concentrations of the alloy components and point defects; G_∇ is a gradient energy part and the last term G_el corresponds to an elastic contribution; V_m is the molar volume; integration is taken over the whole volume V.

The free energy density G_ch can be decomposed into two terms, describing atomic and defect subsystems as follows: G_ch = G_Fe−Cr−Al + G_d. The molar Gibbs energy of the atomic subsystem G_Fe−Cr−Al is defined as: , where is the reference Gibbs energy, denotes the ideal Gibbs energy contribution due to random mixing of atoms, G^ex is the excess term defining a deviation from ideality. The reference Gibbs energy is written in the standard form . The corresponding expressions for the reference energies can be obtained with the help of phase diagrams by using CALPHAD method (Dinsdale, 1991). The ideal part G^id relates to entropic contribution G^id = RT∑_μx_μ ln x_μ, where R is the universal gas constant, T is the temperature. The excess term G^ex is defined by temperature dependent interaction coefficients L_μν = L_μν(T) in the following form: G^ex = ∑_μ≠νx_μx_νL_μ,ν. The Gibbs energy G_d for a point defect subsystem d = {i, v} can be written as follows: . Here is the defect formation energy. It is defined through the formation energies of defects in pure materials Fe, Cr and Al and mean concentrations of atoms in the alloy as , where relates to nominal concentration of μ. The corresponding entropic contribution is . Interactions between defects with host atoms and alloying elements ( and ) can be described by the component , where for the defect-atom interaction energies we use definitions (Was, 2016): . Here denotes the corresponding cohesive energy, Z relates to the coordination number. Point defects formation energies for pure materials and the corresponding cohesive energies can be taken from the literature or can be calculated by ab initio calculations or molecular dynamics simulations.

By combining all above expressions, the total Gibbs energy for the ternary alloy Fe − Cr − Al with point defects takes the form:

The gradient term G_∇ has the form: . The energy constants {κ_μ} are defined according to the Hilliard approach Ref. (Hilliard, 1970). Next we assume κ_μ = κ, where κ = L_Fe,Cra²/6, a = ∑_μa_μx_μ, a_μ is the lattice constant for pure materials. In order to avoid negative gradient energy and ensure the disorder of the system, for the energy constants κ_i and κ_v we use the formalism discussed in Refs. (Liang et al., 2018; Yan et al., 2020) and define them as κ_d = κ_v = κ_i.

The α − α′-interface can be considered as a semi-coherent due to small difference in lattice parameters of two main elements Fe and Cr, as usual (Ribis and Lozano-Perez, 2012). To make our description more general we consider the elastic contribution G_el in the Gibbs energy has the standard form , here C_ijkl denotes the elastic constants tensor, is the elastic strain tensor, where u = (u_x, u_y, u_z) is the displacement vector. By using the Vegard law one defines the lattice parameters misfit caused by introducing alloying element μ. Hence, the eigenstrain tensor acquires the general form , here δ_ij is the Kronecker delta-symbol, . By taking into account that α′-phase is hard phase, one can consider only Chromium contribution to elastic effects. In such a case eigenstrain and corresponding stress tensors acquire the form , with . The elastic energy can be decomposed into four parts: related to eigenstrain components (proportional to squared concentration of alloying elements); related to elastic components defined through elementary strains e_s; is responsible for coupling between deformation elastic deformations e₁ and alloying element concentration; corresponds to cubic symmetry.

As was shown previously (Cahn, 1962; Onuki, 2002; Khachaturyan, 2013) in BCC metals/alloys the minimum free energy is along the soft direction ⟨100⟩. The contribution to elastic energy coming from the lattice mismatch takes the form (Cahn, 1961): , , where Y_[100] is the averaged stiffness parameter for cubic crystals related to the free energy minimum along ⟨100⟩ direction.

By using standard relations between stiffness components and elastic moduli one can write a contribution coming from elastic continuum as follows (Onuki, 2002):where the stretching (Φ) and shear (Ψ) contributions are defined by

Here elementary strain components are: e₁ = e_xx + e_yy + e_zz, , e₅ = e_yz, e₆ = e_zx, . For small strains we get , and , as usual.

The component describing coupling between deformation and concentration has the form (Onuki, 2002): , where . Elastic misfit results in anisotropy with a parameter ξ_a = (C₁₁ − C₁₂ − 2C₄₄)/C₄₄ defining anisotropy contribution Onuki (2002).where ∇_i∇_jw = n_in_jδx_Cr (r, t), n_i denotes component of the interface normal n.

In our calculations we will follow approach developed by Onuki [see Ref. (Onuki, 2002) and citations therein] and results of molecular dynamics (Dai et al., 2022) and assume that all elastic constants depend on concentration of alloying elements: . A condition of mechanical equilibrium implies .

By combining all above expressions, the elastic contribution G_el to the total Gibbs energy for the ternary alloy Fe − Cr − Al becomes the form:

In our approach we use concentration dependent stiffness constants K = K₀ + K₁x_Cr, μ₂ = μ₂₀ + μ₂₁x_Cr, μ₃ = μ₃₀ + μ₃₁x_Cr, , where bulk, shear moduli and anisotropy components are collected in Table 1. For stiffness parameter Y_[100], next, we use elastic constants defined through the nominal compositions, i.e., . By using the Hook law for elastic stress tensor one finds for diagonal components σ_ii: σ_xx = σ_yy = σ_zz = (K₀ + K₁x)e₁, for off–diagonal ones we get σ_xy = (μ₃₀ + μ₃₁x_Cr) sin (2πe₄)/2π, σ_yz = (μ₃₀ + μ₃₁x_Cr) sin (2πe₅)/2π, σ_xy = (μ₃₀ + μ₃₁x_Cr) sin (2πe₆)/2π.

2.2 Evolution equations

By taking into account the mass conservation law and defining the concentration of the Iron as the host element as: x_Fe = 1 − x_Cr − x_Al dynamics of the complete system will be governed by the following equations of the Cahn-Hilliard type (Cahn, 1962; Hilliard, 1970):

Mobility coefficients M_Cr,Cr, M_Al,Al and M_Cr,Al are as follows (Huang et al., 1995; Wu et al., 2001; Wu et al., 2004; Lee et al., 2020):where M_Fe, M_Al and M_Cr are the corresponding mobilities for pure elements, which are defined in the standard manner: M_μ = D_μ/RT, D_μ is the corresponding atomic diffusivity. The mobility for the ensembles of vacancies and interstitials we choose in the standard manner: L_d = D_dc_d/RT, defined through the corresponding point defect diffusivity D_d. Here we use a simplified model for point defects mobilities.

In the Cahn-Hilliard-Cook equation Eq. 7 with a variable mobility Eq. 8, it is not convenient to accurately discretize the gradient operator and divergence operator by using the finite-difference method. However, the Fourier spectral method developed in Refs. (Chen and Shen, 1998; Zhu et al., 1999; Kim, 2007; Canuto et al., 2012; Biner, 2017) allows us to solve numerically Eq. 7 without much more difficulty than the case with a constant mobility. In the framework of the Fourier spectral method the system of coupled non-linear differential equations takes the form:

Here, x_Cr,Al (k, t), c_v,i (k, t), and represent the Fourier transforms of x_Cr,Al(r, t), c_v,i(r, t), and in r-space, respectively. k is the vector in the Fourier space with the magnitude ; { }_k is the Fourier transform and [ ]_r means the inverse Fourier transform to real space.

In further modeling, we use thermodynamic and elastic parameters characterizing ternary systems collected in Table 1.

3 Results and discussions

3.1 Stability diagram

In order to define ranges of concentrations of alloying additions Cr and Al and temperature, when phase decomposition is possible we exploit the linear stability analysis of a homogeneous state to inhomogeneous perturbations. According to obtained kinetic equations one can find phase diagram defining critical values of concentration of alloying elements at fixed temperatures shown in Figure 1. It follows that phase separation related to instability of the solid solution to inhomogeneous perturbations is possible at fixed range of Chromium and Aluminum concentrations. With the temperature growth the instability domain bounded by spinodal curves shrinks. Outside of the corresponding spinodal the system is stable to homogeneous perturbations. In this domains the solid solution is realized. It is clear from this diagram that in the case of high values of the annealing temperature (see, for example, the curve at T = 850 K), when precipitation is impossible in the binary Fe-Cr systems, the addition of Aluminum with a concentration induces such processes. At high temperatures (see curve at T = 900 K) with the temperature increase phase decomposition related to formation of α′-precipitates enriched by Cr is possible only at elevated concentration of Al and with content of Cr around 50 at%.

FIGURE 1

3.2 Numerical analysis

To study process of phase decomposition at annealing numerically we choose values of nominal concentration of two alloying elements and fix the annealing temperature according to stability diagram. We take solid solution as initial microstructure of the corresponding alloy with homogeneously distributed concentration fields of all components including point defects. As initial configurations we use: , ; , , where relate to equilibrium point defect concentration. All simulations were done on square lattice L × L of the linear size L = NΔx with N (N = 128) nodes in each direction and effective lattice constant for simulations Δx = 1 by exploiting Fourier spectral method (Chen and Shen, 1998; Zhu et al., 1999; Canuto et al., 2012; Biner, 2017). We take Δt′ = 10⁻³ as dimensionless time step of integration. Our simulations were done in dimensionless time t′ = tℓ²/D_Al, for spatial coordinate we use r′ = r/ℓ, where ℓ is of the lattice parameter a. Boundary conditions were periodic. We discuss dynamics of α′-phase formation and growth and analyze effect of alloying elements and annealing temperature onto main statistical properties of precipitates.

According to provided stability analysis in order to study α′-phase precipitation next we consider alloy Fe − 30%Cr − 5%Al as a model system with annealing temperature T = 710 K. The main indicator of phase separation is the growth of the dispersion of the concentration field, determined in the standard way: . This quantity plays the role of an effective order parameter during phase decomposition. Indeed, if grows in time, it means that enriched and depleted domains of μ-component are formed and can be distinguished by concentration values, therefore phase decomposition occurs (García-Ojalvo and Sancho, 2012).

Typical scenario of phase decomposition is shown in Figure 2. Snapshots of the concentration field of Cr and Al are shown in Figure 2A at time instants related to markers indicated at protocols of effective order parameters in Figure 2B. Form Figure 2A it follows that after a certain incubation period (t < t_c) domains enriched and depleted in the concentration of impurities emerge from the initial homogeneous distribution (solid solution). Thus, after 22 h of annealing, fine Chromium precipitates can be observed. At the same time, Aluminum is homogeneously distributed in the matrix (outside the precipitates). The number of precipitates and their size increase with time. During this nucleation stage, all formed precipitates are characterized by almost the same size. With further annealing one passes to the coarsening stage, where precipitates continue to grow and interact with each other: small precipitates dissolve, while large ones take material from the matrix and continue to grow following Ostwald ripening scenario. Detailed analysis of α′-phase precipitation can be seen from protocols of the order parameters shown in Figure 2B. One finds that after a certain incubation period dispersion of both x_Cr and x_Al grow rapidly, meaning a passage to the nucleation stage. At the coarsening stage the growth rate of dispersions and decreases, they continue to increase slowly and the system goes to a quasi-stationary regime. Here the total amount of precipitates decreases with time, and the dispersion in their sizes increases.

FIGURE 2

Results of statistical analysis of precipitation are shown in Figure 3. Protocols for the Chromium concentration in the α′-phase and the volume fraction of precipitates are shown in Figure 3A. Both values grow rapidly after the incubation period t_c at the nucleation stage whilst their slow growth is observed at coarsening stage. At the physical time of hours almost 70% of Chromium is concentrated in the α′-phase, whilst the volume fraction of precipitates is about 20%. Obtained data relate well to experimental observations and results of numerical modeling Yang et al. (2020), Chen et al. (2019), Lee et al. (2020). Figure 3B illustrates dynamics of the mean precipitate radius and their number density recalculated from obtained data. It can be seen that the number of precipitates and their mean radius rapidly increase with time at the nucleation stage. At the coarsening stage the mean precipitate size (radius) ⟨R_p⟩ increases according to the power-law ⟨R_p⟩∝ tⁿ with the exponent n close to the classical Lifshits-Slyozov-Wagner (LSW) value 1/3; the total amount of precipitates at this stage decreases exponentially (see power-law and exponential decay fitting shown by dash and dash-dot lines, with fitting parameter τ). Data obtained for the mean precipitate size (around 1–1.2 nm at annealing time 300 h) and precipitate number density (1.5 × 10¹⁹cm⁻³) for the studied alloy relate well to experimental and numerical data for class of system related to Fe-Cr-Al alloys (Messoloras et al., 1984; Chen et al., 2019; Lee et al., 2020; Yang et al., 2020). Scaling dynamics of mean precipitate size in Fe-Cr alloys was observed experimentally [see, for example, Refs. (Brenner et al., 1982; Messoloras et al., 1984; Bley, 1992; Hyde et al., 1995; Soriano-Vargas et al., 2010)], where the scaling exponent n depends on the nominal composition and annealing temperature.

FIGURE 3

Size distribution α′-phase precipitates f (R_p/⟨R_p⟩) is shown in Figure 3C at different time instants by markers. It can be seen that data obtained from the numerical experiment at different times relate to each other meaning an universality of the size distribution of precipitate sizes at the coarsening stage. The solid curve shows the classical universal LSW distribution (Lifshitz and Slyozov, 1961; Wagner, 1961). The obtained numerical data are in good agreement with the LSW distribution. Hence, LSW theory can be used to predicate particle size distribution in Fe-Cr-Al systems at thermal treatment.

Let us consider a behavior of defect subsystem during precipitation. From Figure 4A it follows that the dispersion of the vacancy concentration decreases with time, while the dispersion of the concentration of interstitials slowly increase after the decreasing. Spatial rearrangement of point defects is well seen from snapshots in Figure 4A (see insets). Here vacancies are quasi-homogeneously distributed in a bulk, while interstitials are mostly localized in small precipitates (R_p < ⟨R_p⟩) and near the phase interface in large precipitates (R_p > ⟨R_p⟩). Detailed information about defect distribution around precipitates can be found in profiles of alloying elements and defects shown in Figure 4B. Here one finds the following: the interstitial concentration is higher in precipitates compared to its value in the matrix phase; they are mostly localized at regions of precipitates with elevated curvature (see the snapshot of c_i(r) in Figure 4A); vacancy concentration is locally increased around interfaces and has lower values in precipitates comparing to their values in matrix that relates to thermodynamics; α′-phase precipitates contain extremely small amount of Al atoms comparing to matrix phase.

FIGURE 4

By studying phase decomposition of alloys with different content of alloying elements and different temperatures one finds that an increase in one of the component can control time scale for precipitation kinetics. By studying α′-phase precipitation it was found that the incubation time t_c for precipitation decreases with content of Cr and Al as is shown in Figures 5A, B. Here we denote critical value for concentration of alloying elements as (from phase diagram in Figure 1). It means that an increase in alloying element content accelerates driven forces for phase decomposition decreasing the incubation time t_c needed for emergence of precipitates. From theory of phase decomposition it is known that the quantity t_c depends on the annealing temperature (Khachaturyan, 2013). The corresponding dependence is shown in Figure 5C for the Fe-30%Cr-5%Al system (as an example), where one gets an infinite time needed to initiate phase decomposition in the vicinity of the critical temperature related to spinodal depicted in the phase diagram in Figure 1.

FIGURE 5

Effect of alloying elements and temperature influence onto dynamics of mean precipitate size, volume fraction of precipitates and precipitate size distribution is shown in Figure 6. In Figures 6A, B we plot the corresponding protocols in physical time scale according to the time related to phase decomposition starting from t_c. From obtained results it follows that an increase Al content results in slight increase in precipitate size and volume fraction (cf. curves 1 and 2), whereas Cr addition increases the corresponding values of above quantities crucially (cf. curves 1 and 3). These results are related well to previous studies [see for example, Ref. (Lee et al., 2020)]. With the temperature decrease from 710 K down to 690 K (cf. curves 1 and 4) one gets low values of the mean precipitate size and low values for volume fraction of α’ precipitates that is constituent with results discussed in Ref. (Chen et al., 2019). Distributions over precipitate size in alloys with different content of Al and Cr and two different temperatures are shown in Figure 6C. One finds that obtained data are well predicted by classical LSW distribution. Therefore, it follows that additional introduction of Cr or Al does not change universality of precipitate size distribution at coarsening stage at different temperatures. By comparing dynamics of the precipitate size and precipitate size distributions in alloys with different content of doping one concludes that both dynamics of the system and statistical properties of growing precipitates remain universal. The difference in measured quantities (mean precipitate size, precipitate number density, precipitate volume fraction) is related mostly to scaling factors determined by nominal composition of the studied alloys.

FIGURE 6

4 Discussions

We need to stress out that the used model for microstructure transformations in Fe-Cr-Al alloy is based on equilibrium CALPHAD free energy, however Aluminum and Chromium concentrations will be out of thermal equilibrium at cooling. Therefore, to describe more precisely free energy in metastable zone one should further incorporate corrections to the free energy related to non-equilibrium states. At the same time this generalization will not affect the final state related to equilibrium as was shown in previous studies (Chen et al., 2019; Lee et al., 2020).

In our study concentration-dependent model for stiffness constant was used for solid systems with cubic symmetry. Obtained results relate well to most of experimental observations. However, for more precise description of the system dynamics the measurement of the interface thickness obtained experimentally allows one to fix the corresponding values and thus the space scale.

For the considered class of concentrated alloys there is no formal theory for the mobilities exploited in the present study. Therefore here the model assumptions for concentration-dependent mobilities were used on the basis of previous works (Huang et al., 1995; Wu et al., 2001; Wu et al., 2004; Lee et al., 2020). This approach allows one to study transient dynamics of phase decomposition process in details, and can play a significant effect on the coarsening kinetics (Lacasta et al., 1992; Lacasta et al., 1993; Bray and Emmott, 1995; Zhu et al., 1999; Kharchenko and Dvornichenko, 2008; Kharchenko and Kharchenko, 2014; Kharchenko et al., 2022). A constant mobility model can be used in order to modeling the late time microstructure.

Nevertheless used in the present study assumptions and simplifications do not affect significantly the obtained results related well to most of experimental and theoretical predictions (Messoloras et al., 1984; Chen et al., 2019; Lee et al., 2020; Yang et al., 2020).

5 Conclusion

We have developed a phase field model to study dynamics of microstructure transformations and the evolution of defect structure during heat treatment of Fe-Cr-Al systems. Statistical and kinetic properties of evolving microstructure and defect structure in alloys with different content of alloying elements and different temperatures were studied.

In the framework of the stability analysis we obtain a phase diagram illustrating the range of concentrations of Chromium/Aluminum and annealing temperature when phase separation in Fe-Cr-Al systems accompanied by α′-phase precipitation can be realized. We exploited numerical simulations to discuss dynamics of precipitates formation, their growth, evolution of their number densities, precipitate mean size, and point defects rearrangement at precipitation. It is shown that mean precipitate size grows according to the power-law, whereas precipitate number density manifests exponential decay at coarsening stage. Precipitate size distribution function is of universal character.

It is shown that addition of Aluminum up to 8% initiates phase decomposition decreasing the incubation time for development of concentration waves. It slightly increases the mean precipitate size and volume fraction of precipitates. An increase in the nominal concentration of Chromium from 24% up to 40% increases content of Cr inside precipitates and essentially increases precipitate size with their volume fraction and accelerates phase decomposition processes. Dynamics of growing precipitates and the corresponding distribution functions over precipitate size manifest universality with scaling factors depending on the nominal compositions of alloys.

It was found that the annealing at elevated temperatures results in a decrease in the precipitates number density and in an increase their size. It was found that at high and low temperatures, precipitation processes are slowed down compared to the case of moderate temperatures.

By studying kinetics of equilibrium point defects it was shown that during annealing of solid solution vacancies are mostly homogeneously distributed in a bulk, while interstitial atoms are mostly localized in small α′ precipitates and inside large precipitates near the interfaces with large curvature.

We expect that obtained results can be used to predict materials properties change at different content of alloying elements and will be useful as basic data for microstructure optimization of Fe-Cr-Al fuel cladding. Results can be used to describe radiation induced effects by incorporating ballistic mixing and dynamics of non-equilibrium defects produced by irradiation.

References

• 1

  BinerS. B. (2017). Programming phase-field modeling. Switzerland: Springer International Publishing.

  • Google Scholar
  • View reference in article

• 2

  BleyF. (1992). Neutron small-angle scattering study of unmixing in fe-cr alloys. Acta metallurgica materialia40, 1505–1517. 10.1016/0956-7151(92)90094-u

  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  BrayA.EmmottC. (1995). Lifshitz-slyozov scaling for late-stage coarsening with an order-parameter-dependent mobility. Phys. Rev. B52, R685–R688. 10.1103/physrevb.52.r685

  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  BrennerS.MillerM.SoffaW. (1982). Spinodal decomposition of iron-32 at.% chromium at 470 c. Scr. Metall.16, 831–836. 10.1016/0036-9748(82)90239-3

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  CahnJ. W. (1961). On spinodal decomposition. Acta metall.9, 795–801. 10.1016/0001-6160(61)90182-1

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  CahnJ. W. (1962). On spinodal decomposition in cubic crystals. Acta metall.10, 179–183. 10.1016/0001-6160(62)90114-1

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  CanutoC.HussainiM. Y.QuarteroniA.ThomasA.Jr (2012). Spectral methods in fluid dynamics. Berlin, Heidelberg: Springer Science & Business Media.

  • Google Scholar
  • View reference in article

• 8

  ChangK.MengF.GeF.ZhaoG.DuS.HuangF. (2019). Theory-guided bottom-up design of the fecral alloys as accident tolerant fuel cladding materials. J. Nucl. Mater.516, 63–72. 10.1016/j.jnucmat.2019.01.002

  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  ChenL. Q.ShenJ. (1998). Applications of semi-implicit Fourier-spectral method to phase field equations. Comput. Phys. Commun.108, 147–158. 10.1016/s0010-4655(97)00115-x

  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  ChenS.LiY.ShiS.JinS. (2019). Quantitative phase-field simulation of composition partition and separation kinetics of nanoscale phase in fe-cr-al alloy. J. Nanomater.2019, 1–11. 10.1155/2019/6862390

  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  DaiH.YuM.DongY.SetyawanW.GaoN.WangX. (2022). Effect of cr and al on elastic constants of fecral alloys investigated by molecular dynamics method. Metals12, 558. 10.3390/met12040558

  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  DieterG. E.BaconD. (1976). Mechanical metallurgy. New York: McGraw-Hill.

  • Google Scholar
  • View reference in article

• 13

  DinsdaleA. T. (1991). Sgte data for pure elements. Calphad15, 317–425. 10.1016/0364-5916(91)90030-n

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  EjenstamJ.ThuvanderM.OlssonP.RaveF.SzakalosP. (2015). Microstructural stability of fe–cr–al alloys at 450–550 c. J. Nucl. Mater.457, 291–297. 10.1016/j.jnucmat.2014.11.101

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  García-OjalvoJ.SanchoJ. (2012). Noise in spatially extended systems. New York, NY: Springer Science & Business Media.

  • Google Scholar
  • View reference in article

• 16

  GulbransenE. A.AndrewK. F. (1959). Oxidation studies on the iron-chromium-aluminum heater alloys. J. Electrochem. Soc.106, 294. 10.1149/1.2427333

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  HilliardJ. E. (1970). Spinodal Decomposition, in Phase Transformations. American Society of Metals, Metals Park, 497.

  • Google Scholar
  • View reference in article

• 18

  HuangC.de La CruzM. O.SwiftB. (1995). Phase separation of ternary mixtures: Symmetric polymer blends. Macromolecules28, 7996–8005. 10.1021/ma00128a005

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  HydeJ.MillerM.HetheringtonM.CerezoA.SmithG.ElliottC. (1995). Spinodal decomposition in fe-cr alloys: Experimental study at the atomic level and comparison with computer models-ii. development of domain size and composition amplitude. Acta metallurgica materialia43, 3403–3413. 10.1016/0956-7151(95)00041-s

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  KhachaturyanA. G. (2013). Theory of structural transformations in solids. New York, NY: Dover Publications, Inc.,.

  • Google Scholar
  • View reference in article

• 21

  KharchenkoD.DvornichenkoA. (2008). Phase separation in binary systems with internal multiplicative noise. Phys. A Stat. Mech. its Appl.387, 5342–5354. 10.1016/j.physa.2008.05.041

  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  KharchenkoV. O.KharchenkoD. O. (2014). Abnormal grain growth in nonequilibrium systems: Effects of point defect patterning. Phys. Rev. E89, 042133. 10.1103/physreve.89.042133

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  KharchenkoV.XinT.WuL.KharchenkoD.KupriienkoV.ShudaI. (2022). Phase stability and precipitation modeling in neutron irradiated zr–2% nb alloy. Model. Simul. Mater. Sci. Eng.30, 075006. 10.1088/1361-651x/ac8fad

  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  KimJ. (2007). A numerical method for the Cahn–Hilliard equation with a variable mobility. Commun. Nonlinear Sci. Numer. Simul.12, 1560–1571. 10.1016/j.cnsns.2006.02.010

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  KimS.BuyersW. (1978). Vacancy formation energy in iron by positron annihilation. J. Phys. F Metal Phys.8, L103–L108. 10.1088/0305-4608/8/5/001

  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  KittelC.McEuenP. (2018). Introduction to solid state physicsGlobal ed., 9th ed.New Jersey: Wiley, 692.

  • Google Scholar
  • View reference in article

• 27

  KobayashiS.TakasugiT. (2010). Mapping of 475 c embrittlement in ferritic fe–cr–al alloys. Scr. Mater.63, 1104–1107. 10.1016/j.scriptamat.2010.08.015

  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  KresseG.FurthmüllerJ. (1996). Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B54, 11169–11186. 10.1103/physrevb.54.11169

  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  KresseG.HafnerJ. (1993). Ab initio molecular dynamics for liquid metals. Phys. Rev. B47, 558–561. 10.1103/physrevb.47.558

  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  KresseG.HafnerJ. (1994). Norm-conserving and ultrasoft pseudopotentials for first-row and transition elements. J. Phys. Condens. Matter6, 8245–8257. 10.1088/0953-8984/6/40/015

  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  LacastaA.Hernández-MachadoA.SanchoJ. M.ToralR. (1992). Domain growth in binary mixtures at low temperatures. Phys. Rev. B45, 5276–5281. 10.1103/physrevb.45.5276

  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  LacastaA.SanchoJ. M.HernaA.ToralR. (1993). Effects of domain morphology in phase-separation dynamics at low temperature. Phys. Rev. B48, 6854–6857. 10.1103/physrevb.48.6854

  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  LeeJ.ParkK.ChangK. (2020). Effect of al concentration on the microstructural evolution of fe-cr-al systems: A phase-field approach. Metals11, 4. 10.3390/met11010004

  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  LiW.LuS.HuQ. M.MaoH.JohanssonB.VitosL. (2013). The effect of al on the 475 c embrittlement of fe–cr alloys. Comput. Mater. Sci.74, 101–106. 10.1016/j.commatsci.2013.03.021

  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  LiangL.MeiZ. G.KimY. S.AnitescuM.YacoutA. M. (2018). Three-dimensional phase-field simulations of intragranular gas bubble evolution in irradiated u-mo fuel. Comput. Mater. Sci.145, 86–95. 10.1016/j.commatsci.2017.12.061

  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  LifshitzI. M.SlyozovV. V. (1961). The kinetics of precipitation from supersaturated solid solutions. J. Phys. Chem. solids19, 35–50. 10.1016/0022-3697(61)90054-3

  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  MathonM.De CarlanY.GeoffroyG.AvertyX.AlamoA.De NovionC. (2003). A sans investigation of the irradiation-enhanced α–α′ phases separation in 7–12 cr martensitic steels. J. Nucl. Mater.312, 236–248. 10.1016/s0022-3115(02)01630-6

  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  MessolorasS.PikeB.StewartR.WindsorC. (1984). Precipitation in iron–chromium–aluminium alloys. Metal Sci.18, 311–321. 10.1179/030634584790419999

  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  NagasakiS. (2004). Metals data book. Tokyo: Japan Institute of Metals.

  • Google Scholar
  • View reference in article

• 40

  OgorodnikovV.RakitskiiA.RogovoiY. I. (1988). Calculation of the vacancy formation energy of metals. Sov. Powder Metall. Mater. Ceram.(Engl. Transl.)27, 55–60. 10.1007/bf00799739

  • CrossRef
  • Google Scholar
  • View reference in article

• 41

  OnukiA. (2002). Phase transition dynamics. Cambridge: Cambridge University Press.

  • Google Scholar
  • View reference in article

• 42

  QuH.LangY.YaoC.ChenH.YangC. (2013). The effect of heat treatment on recrystallized microstructure, precipitation and ductility of hot-rolled fe–cr–al–rem ferritic stainless steel sheets. Mater. Sci. Eng. A562, 9–16. 10.1016/j.msea.2012.11.008

  • CrossRef
  • Google Scholar
  • View reference in article

• 43

  RibisJ.Lozano-PerezS. (2012). Orientation relationships and interface structure of α′-cr nanoclusters embedded in α-fe matrix after α − −α′ demixing in neutron irradiated oxide dispersion strengthened material. Mater. Lett.74, 143–146. 10.1016/j.matlet.2012.01.115

  • CrossRef
  • Google Scholar
  • View reference in article

• 44

  Soriano-VargasO.Avila-DavilaE. O.Lopez-HirataV. M.Cayetano-CastroN.Gonzalez-VelazquezJ. L. (2010). Effect of spinodal decomposition on the mechanical behavior of fe–cr alloys. Mater. Sci. Eng. A527, 2910–2914. 10.1016/j.msea.2010.01.020

  • CrossRef
  • Google Scholar
  • View reference in article

• 45

  StottF.WoodG.StringerJ. (1995). The influence of alloying elements on the development and maintenance of protective scales. Oxid. metals44, 113–145. 10.1007/bf01046725

  • CrossRef
  • Google Scholar
  • View reference in article

• 46

  TerentyevD.Hafez HaghighatS.SchäublinR. (2010). Strengthening due to cr-rich precipitates in fe–cr alloys: Effect of temperature and precipitate composition. J. Appl. Phys.107, 061806. 10.1063/1.3340522

  • CrossRef
  • Google Scholar
  • View reference in article

• 47

  TerentyevD.OlssonP.KlaverT.MalerbaL. (2008). On the migration and trapping of single self-interstitial atoms in dilute and concentrated fe–cr alloys: Atomistic study and comparison with resistivity recovery experiments. Comput. Mater. Sci.43, 1183–1192. 10.1016/j.commatsci.2008.03.013

  • CrossRef
  • Google Scholar
  • View reference in article

• 48

  WagnerC. (1961). Theorie der Alterung von Niederschlägen durch Umlösen (Ostwald‐Reifung). Z. Elektrochem65, 581–591. 10.1002/bbpc.19610650704

  • CrossRef
  • Google Scholar
  • View reference in article

• 49

  WasG. S. (2016). Fundamentals of radiation materials science: Metals and alloys. Berlin, Heidelberg, New York: Springer.

  • Google Scholar
  • View reference in article

• 50

  WuK.MorralJ.WangY. (2001). A phase field study of microstructural changes due to the kirkendall effect in two-phase diffusion couples. Acta mater.49, 3401–3408. 10.1016/s1359-6454(01)00257-9

  • CrossRef
  • Google Scholar
  • View reference in article

• 51

  WuK.MorralJ.WangY. (2004). Movement of kirkendall markers, second phase particles and the type 0 boundary in two-phase diffusion couple simulations. Acta mater.52, 1917–1925. 10.1016/j.actamat.2003.12.031

  • CrossRef
  • Google Scholar
  • View reference in article

• 52

  WukusickC. S.CollinsJ. F. (1964). An iron-chromium-aluminum alloy containing yttrium. Mat. Res. std.4 (1964), 637.

  • Google Scholar
  • View reference in article

• 53

  YamamotoY.PintB.TerraniK.FieldK.YangY.SneadL. (2015). Development and property evaluation of nuclear grade wrought fecral fuel cladding for light water reactors. J. Nucl. Mater.467, 703–716. 10.1016/j.jnucmat.2015.10.019

  • CrossRef
  • Google Scholar
  • View reference in article

• 54

  YanZ.ShiS.LiY.ChenJ.MaqboolS. (2020). Vacancy and interstitial atom evolution with the separation of the nanoscale phase in fe–cr alloys: Phase-field simulations. Phys. Chem. Chem. Phys.22, 3611–3619. 10.1039/c9cp06247e

  • CrossRef
  • Google Scholar
  • View reference in article

• 55

  YangZ.WangZ.XiaC.OuyangM.PengJ.ZhangH.et al (2020). Aluminum suppression of α precipitate in model fe–cr–al alloys during long-term aging at 475 c. Mater. Sci. Eng. A772, 138714. 10.1016/j.msea.2019.138714

  • CrossRef
  • Google Scholar
  • View reference in article

• 56

  ZhuJ.ChenL. Q.ShenJ.TikareV. (1999). Coarsening kinetics from a variable-mobility Cahn-Hilliard equation: Application of a semi-implicit Fourier spectral method. Phys. Rev. E60, 3564–3572. 10.1103/physreve.60.3564

  • CrossRef
  • Google Scholar
  • View reference in article