Next Article in Journal
Synthesis and Investigation of Physicochemical Properties and Biocompatibility of Phosphate–Vanadate Hydroxyapatite Co-Doped with Tb3+ and Sr2+ Ions
Next Article in Special Issue
Numerical Investigation of GaN HEMT Terahertz Detection Model Considering Multiple Scattering Mechanisms
Previous Article in Journal
Amphiphilic Bowl-Shaped Janus Particles Prepared via Thiol–Ene Click Reaction for Effective Oil–Water Separation
Previous Article in Special Issue
Classical Analog and Hybrid Metamaterials of Tunable Multiple-Band Electromagnetic Induced Transparency
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Elastic Strain Relaxation of Phase Boundary of α′ Nanoscale Phase Mediated via the Point Defects Loop under Normal Strain

1
School of Materials Science and Engineering, Nanjing University of Science and Technology, Nanjing 210094, China
2
MIIT Key Laboratory of Advanced Metallic and Intermetallic Materials Technology, Nanjing 210094, China
*
Author to whom correspondence should be addressed.
Nanomaterials 2023, 13(3), 456; https://doi.org/10.3390/nano13030456
Submission received: 21 November 2022 / Revised: 30 December 2022 / Accepted: 19 January 2023 / Published: 22 January 2023

Abstract

:
Irradiation-induced point defects and applied stress affect the concentration distribution and morphology evolution of the nanophase in Fe–Cr based alloys; the aggregation of point defects and the nanoscale precipitates can intensify the hardness and embrittlement of the alloy. The influence of normal strain on the coevolution of point defects and the Cr-enriched α′ nanophase are studied in Fe-35 at.% Cr alloy by utilizing the multi-phase-field simulation. The clustering of point defects and the splitting of nanoscale particles are clearly presented under normal strain. The defects loop formed at the α/α′ phase interface relaxes the coherent strain between the α/α′ phases, reducing the elongation of the Cr-enriched α′ phase under the normal strains. Furthermore, the point defects enhance the concentration clustering of the α′ phase, and this is more obvious under the compressive strain at high temperature. The larger normal strain can induce the splitting of an α′ nanoparticle with the nonequilibrium concentration in the early precipitation stage. The clustering and migration of point defects provide the diffusion channels of Cr atoms to accelerate the phase separation. The interaction of point defect with the solution atom clusters under normal strain provides an atomic scale view on the microstructure evolution under external stress.

Graphical Abstract

1. Introduction

An irradiation cascade reaction produces point defects in crystal materials. Point defects are either intrinsic, as with vacancies, interstitial atoms, and their small clusters, or extrinsic, as with impurities and dopants; they play a major role in materials’ properties and microstructure evolution. Some properties of these point defects, such as their formation and migration energies, are mainly determined by the region in the immediate vicinity of defects where the crystal structure is strongly disturbed [1]. The formation of point defects not only promotes the phase separation [2,3], but also produces the lattice misfit and the elastic energy [4] to affect the morphology of the precipitated phase.
Furthermore, point defects can induce a long-range perturbation of the host lattice, leading to an elastic interaction with other structural defects, impurities, or an applied elastic field. Following the work of Eshelby et al. [5], the elastic model of a point defect can be regarded as a spherical inclusion that is forced into a spherical hole of slightly different size in an infinite elastic medium. The variations in the “size” and “shape” of point defects with the strain field is described by assigning different elastic constants to the inclusion. In the study of Clouet et al. [1], the elastic model was used to describe the evolution of a point defect in an external elastic field, and the influence of a defect was discussed using the dipole tensors.
The strain produced by the thermal or external stress has a substantial effect on the phase separation and spatial morphology [6,7]. In the study of a monolayer alloy, the strain can induce a structural phase transformation [8]. Furthermore, applied tensile or compressive stresses can induce the orientation morphologies of precipitates [9,10], and directional coarsening is observed when annealed with tensile and compressive stresses [11,12]. Prikhodko and Ardell found that applied compressive stress retards the growth of the precipitates and widens the distribution of particle size [13,14]. However, the interaction between the point defects and the nanophase under the applied strain is still not clear; also, the morphology of the defects cluster and its effect on the evolution of nanophase are not studied.
Fe–Cr-based stainless steels are used as the critical material in nuclear power plants for their excellent resistance to oxidation, corrosion resistance, and mechanical properties; both the Fe-enriched α and the Cr-enriched α′ phase with bcc structure are formed via phase separation under high temperature, irradiation, or stress. Tong et al. studied the morphology and kinetics evolution of the Cr-enriched α′ nanophase under multiaxial strain by utilizing the three-dimensional phase-field simulation; they found that phase separation is accelerated under unidirectional normal strain, and is only slowed down under triaxial compressive strain [15]. When the precipitate’s size is at the nanoscale, the irradiation point defects, such as vacancy and interstitial atoms, will essentially affect the atoms’ migrations, and then the formation of a Cr-enriched α′ phase will be faster [16,17]. These Cr-enriched α′ nanoscale precipitates can prevent the motion of dislocations, which causes an increase in hardness [18]; thus, unfavorable thermal embrittlement will be produced [19]. However, the evolutions of nanoparticles and point defects in Fe–Cr alloys under normal stress and radiation are not clarified. Therefore, the effect of strain accompanied with irradiated defects on the precipitation is an interesting question.
In addition, temperature influences the phase separation under strain and irradiation. Zhu et al. found that the effect of applied strain is more obvious at high temperature [20], indicating that the initial stage phase separation is faster and that the elongation of the particle is obvious at high temperature under applied strain. The strain can induce the splitting of coherent precipitates in Ni-based superalloys close to the solvus temperature [21,22], and the particle splitting based on the competition of the coherency elastic energy and the interfacial energy was discussed by Lee [23]. This splitting of the γ′ particle is caused by the stress relaxation of the internal misfit strain between the γ’ and γ phases, and an applied load affects the splitting direction [24]. However, the splitting of the particle with the influence of point defects under the difference temperature is not well studied. Therefore, the effect of temperature on the evolution of defects and the splitting of nanoparticles needs further study in the strained system.
In view of previous studies regarding the irradiation defects and strains in phase separation, the evolution of point defects and nanophase under the normal strain will be studied in this work. The phase-field model is developed coupling with the production and recombination of vacancy and interstitial atoms, and the external strains with the coherent strain are incorporated. The effects of normal strain on the concentration evolution of point defects and Cr in the Cr-enriched nanophase, and the shape change of the nanophase around the defects loop are revealed. In addition, the splitting of nanoparticles with point defects is presented with the temperature change.

2. Model and Methods

2.1. Elastic Strain Including Vacancy and Interstitial Atoms

The strain ε i j a caused by the external stress influences the morphology of the nanophase, which is introduced as part of the elastic strain. Therefore, the elastic strain can be calculated by [25]
ε i j el = ε i j a + ε i j ε i j 0
where ε i j is the heterogeneous strain, and ε i j 0 is the eigenstrain caused by the compositional heterogeneity and point defects and can be calculated by [4]
ε i j 0 = δ i j ( ε 0 Δ c + ε V 0 Δ c V + ε I 0 Δ c I )
where δ i j is the Kronecker-delta function. ε 0 = 1 a ( d a d c Cr ) , ε V 0 = 1 a ( d a d c V ) and ε I 0 = 1 a ( d a d c I ) are the expansion coefficients of the lattice parameter due to the introduction of concentration inhomogeneity, vacancy, and interstitial atoms [4], Δ c = c Cr c 0 , where c Cr is the concentration of Cr, c 0 is the initial average concentration of Cr, Δ c m = c m c m 0 ( m = V , I ) , c m is the concentration of vacancy and interstitial atoms, and c m 0 is the initial average concentration of vacancy and interstitial atoms and can be calculated via [26]
c m 0 = exp ( E m f R T )
where E m f is the formation energy of defect. The vacancy and interstitial atoms are assumed to occur concurrently; therefore, the concentrations of solution atoms and point defects satisfy c Cr + c V + c I + c Fe = 1.0 ; c Fe is the concentration of Fe. The variation of the stress-free lattice parameter is a = ( 1 c Cr c V c I ) a Fe + c Cr a Cr , where a Cr and a Fe are the lattice parameters of Cr and Fe atoms. The elastic strain can be calculated using Hook’s law [27]
σ i j el = C i j k l ε i j el = ( C i j k l 0 + Δ C i j k l Δ c ) ( ε i j a + ε i j ε i j 0 )
where C i j k l 0 is the average elastic modulus tensors and can be calculated by
C i j k l 0 = λ C i j k l P + ( 1 λ ) C i j k l M
where λ is the volume fraction of precipitates, and C i j k l P and C i j k l M are the elastic modulus tensors of precipitates and matrix, respectively. Δ C i j k l is the difference of the elastic modulus tensors and is given by
Δ C i j k l = C i j k l P C i j k l M
In the Fe–Cr alloy system, the elastic constants of the precipitates are C 11 P = 3.5 × 10 11 Pa , C 12 P = 0.678 × 10 11 Pa and C 44 P = 1.008 × 10 11 Pa , and those of the matrix are C 11 M = 2.33 × 10 11 Pa , C 12 M = 1.354 × 10 11 Pa and C 44 M = 1.178 × 10 11 Pa [28]. The elastic strain is related to the displacement
ε k l = 1 2 ( u k x l + u l x k )
Since the mechanical equilibrium with respect to elastic displacements is established much faster than through any diffusion processes, for any given distribution of composition, the system is always at mechanical equilibrium [29].
σ i j el x j = 0
Then, the elastic strain energy density can be calculated by
E e = 1 2 C i j k l ( ε i j a + ε i j ε i j 0 ) ( ε k l a + ε k l ε k l 0 )

2.2. Free Energy

The total free energy [30] includes the chemical free energy, gradient energy, and elastic energy
F ( c Cr , c V , c I ) = V { 1 V m [ G ( c Cr , c V , c I , T ) + k Cr 2 | c Cr | 2 + G mag ] + k V 2 | c V | 2 + k I 2 | c I | 2 + E e ( c Cr , c V , c I , T ) } d V
where G ( c Cr , c V , c I , T ) is the molar Gibbs energy, chemical free density is E c = G / V m , and V m is the molar volume. k Cr is the gradient energy coefficients for the α/α′ phase, where k Cr = 1 6 r 0 2 L FeCr , and L FeCr is the interaction parameter between Fe and Cr, where L FeCr = 20500 9.68 T J mol−1 [31]. r 0 is the interatomic distance and changes with concentration. k V and k I are the gradient energy coefficients for vacancy and interstitial atoms, where k V = k I = 6.91 × 10 9 J   m 1 [4]. Then, the molar Gibbs energy can be written as
G ( c Cr , c V , c I , T ) = E V f c V + E I f c I + c Fe G Fe 0 + c Cr G Cr 0 + L FeCr c Cr c Fe + R T [ c Cr ln c Cr + c V ln c V + c I ln c I + ( 1 c Cr c V c I ) ln ( 1 c Cr c V c I ) ]
where G Fe 0 and G Cr 0 are the molar Gibbs free energies of the element of Fe and Cr, respectively. E V f and E I f are the respective formation energies of vacancy and interstitial atoms, which are taken as 2.24 eV and 3.03 eV [32]. Gmag represents the magnetic ordering contribution to the Gibbs free energy [33]
G mag = R T ln ( β + 1 ) f ( τ )
where β is the atomic magnetic moment in the Bohr magneton, dependent on the concentration of Cr, β = 2.22 c Fe 0.008 c Cr 0.85 c Fe c Cr , and f ( τ ) satisfies the polynomial function as follows [33]
f ( τ ) = { 0.90530 τ 1 + 1.0 0.153 τ 3 6.8 × 10 3 τ 9 1.53 × 10 3 τ 15 , τ < 1 0.06417 τ 5 2.037 × 10 3 τ 15 40278 × 10 4 τ 25 , τ 1

2.3. Phase-Field Model

The effect of irradiation on phase separation is accomplished via the defects action. So, the field variations in the phase-field model are described by the concentration field of Cr c Cr ( r , t ) , vacancy c V ( r , t ) , and interstitial atoms c I ( r , t ) through the Cahn–Hilliard diffusion equation [34,35]
c Cr ( r , t ) t = V m ( M δ F δ c Cr ) + ξ Cr
c V ( r , t ) t = V m ( M V δ F δ c V ) + ξ V + P V R VI S V
c I ( r , t ) t = V m ( M I δ F δ c I ) + ξ I + P I R VI S I
where r and t are the spatial coordinate and time, respectively, and M is the chemical mobility of the Fe–Cr alloy and is given by Darken’s equation [36]
M = ( c Cr M Fe + c Fe M Cr ) c Cr c Fe
where M i ( i = Fe , Cr ) is the atomic mobility of Fe or Cr, and can be expressed as M i = D i R T , where D i ( i = Fe , Cr ) is the diffusion coefficient of Fe or Cr, R is the gas constant, and T is the absolute temperature. M V and M I are the chemical mobilities of vacancy and interstitial atoms, respectively, which can be described by M j = c j D j R T ( j = V , I ) , where D j ( j = V , I ) is the diffusion coefficient of vacancy or interstitial atoms. ξ i ( i = Cr , V , I ) are the thermal fluctuations to introduce the random distribution of concentration and point defects to trigger the nucleation, and their magnitudes are set as 10 8 . P V and P I are the production rates of vacancy and interstitial atoms, respectively, which are taken as 1 × 10 8 dpa / s [37]. R VI is the recombination rate of the vacancy and interstitial atoms [4,38], and is given by
R VI = 4 π r VI ( D V + D I ) Ω c V c I
where r VI is the recombination radius and is equivalent to the atomic radius, r VI = 3 a 4 . Ω is the atomic volume. S V and S I are the sink rates of the vacancy and interstitial atoms, respectively.

2.4. Numerical Calculation

By substituting Equation (10) into Equations (14)–(16) respectively, the evolution equations of the solution atoms and point defects are given by
c Cr ( r , t ) t = [ M ( δ G δ c Cr k Cr 2 c Cr + V m δ E e δ c Cr ) ] + ξ Cr
c V ( r , t ) t = [ M V ( δ G δ c V V m k V 2 c V + V m δ E e δ c V ) ] + ξ V + P V R VI S V
c I ( r , t ) t = [ M I ( δ G δ c I V m k I 2 c I + V m δ E e δ c I ) ] + ξ I + P I R VI S I
For solving Equations (19)–(21) numerically, the following normalization parameters are introduced, r = r / l , t = t D / l 2 , M i = R T 0 M i / D , G = G / R T 0 , k Cr = k Cr / R T 0 l 2 , k V = V m k V / R T 0 l 2 , k I = V m k I / R T 0 l 2 , E c = V m E c / R T 0 , E e = V m E e / R T 0 , = r = ( r / l ) , ξ i = l 2 D ξ i , P i = l 2 D P i , R i = l 2 D R i , S i = l 2 D S i , where l is the grid length of the simulation cell and is set as the average lattice parameter of Cr and Fe, l = 0.29 × 10 9   m . D = 10 26   m 2   s 1 is the diffusion coefficient and is used for dimensionless, and T 0 = 900 K is the critical temperature of spinodal decomposition of the Fe–Cr alloy [39]. The dimensionless grid size is Δ x = Δ y = 1.0 , and the 2D simulation cell is 64 Δ x × 64 Δ y with the grid length Δ x = Δ y = l . The time step is Δ t * = 0.001 . The semi-implicit Fourier spectrum method is used to solve Equations (19)–(21) with a periodic boundary condition. Barker et al. [40] concluded that the qualitative morphology is very similar between the 2D and 3D simulations when they discussed the effect of concentration on the spinodal decomposition of the Fe–Cr alloys. Zhu et al. also showed that the morphology and the dynamic evolution law of the α′ phase are similar for the 2D and 3D simulations under the uniaxial strain [41], although the elasticity is inherently a 3D problem. In addition, in the microstructure evolution of Ni-base superalloys during [110] creep loading, the 3D simulation focuses on the microstructure formation, and more quantitative analyses are performed in 2D simulation [42]. Therefore, the 2D simulation can provide reasonable results and save more time in simulation.

3. Results and Discussion

The simulation starts from a particle with an initial concentration of 35 at.% Cr in a circle region with a radius of 5 l and 25 at.% Cr in the other regions of the simulation cell, which is to ensure that there is only one particle nucleation in the system, and that the vacancy and interstitial atoms are randomly distributed in the simulation cell. The input concentrations of vacancy and interstitial atoms are c V i = c I i = 0 and c V i = c I i = 1 × 10 7 for comparing the effects with and without defects.

3.1. α′ Phase Evolution with a Point Defects Loop under Normal Strain

In this section, the effects of normal strain on the concentration evolution of Cr in the Cr-enriched α′ phase with the vacancy and interstitial atoms at 750K are studied. The tensile and compressive strains show different results on the elongation of the α′ nanophase, and are explained by the distribution of the elastic energy density in the system with defects loop and nanophase.
Figure 1 shows the morphology of the Cr-enriched α′ nanophase and vacancy loop at 750 K under normal trains ε x x a = 0 , ε x x a = 0.025 , and ε x x a = 0.025 . Under the effect of compressive strain ε x x a = 0.025 along the x*-axis, the α′ phase particle elongates along the x*-axis, as shown in Figure 1b. When the tensile strain ε x x a = 0.025 is along the x*-axis, the α′ phase particle is elongated along the y*-axis, as shown in Figure 1c. This orientation of the α′ phase depends on the sign of the eigenstrain and the applied strain [43]. The elongation direction under the tensile or compressive strain depends on whether the phase is hard or soft, and the lattice mismatch of the precipitates is either positive or negative [44]. In the Fe–Cr alloy, the α′ phase is the soft phase and the lattice mismatch is the positive swelling coefficient [20], and so the orientation of the α′ phase is perpendicular to the direction of tensile strain and parallel to the compression strain direction [15,45].
When the point defects loop forms at the α/α′ phase boundaries, as shown in Figure 1g–i, the coherent relationship between the α and α′ phases is released, and the egenstrain is reduced, so that the elongation of the Cr-enriched particle along the tensile direction is weakened under normal strain, as shown in Figure 1e,f. In the elastic model of point defect [1], the relation of defect and elastic dipole is established via Kanzaki forces. Kanzaki forces are applied to the atoms in the neighborhood of the point defect. Taking vacancy as an example, when the atom is added in the vacancy reversely, a force is produced between the host atom and the neighbor atoms, which is contrary to the Kanzaki forces. So, the point defect offsets the force between the neighbor atoms of the normal crystal lattice, and the strain between the neighbor atoms of Fe and Cr is released. As a result, the elongation of the particle around the defects loop is restricted.
To identify the particle shape change under the influences of normal strain and defects loop, the distributions of local Cr concentration, cCr, are shown in Figure 2 along x* = 32 and y* = 32, where, the Lx* and Ly* delegate the concentration width for cCr = 0.5 in the x*-axis and y*-axis directions, respectively. Then, the aspect ratio k of the α′ phase is calculated using the ratio of the particle′s length to width, i.e., k = Lx*/Ly* and k = Ly*/Lx* under the compressive and tensile strains, respectively. When t* = 50, the concentration of Cr reaches the equilibrium value 0.84 and 0.83 with and without strain, as shown in Figure 2.
As the initial concentration of point defect increases from 0.0 to 1.0 × 10–7, when ε x x a = 0 , the aspect ratio k = 1 of the α′ phase, as shown in Figure 2a,d. When ε x x a = 0.025 , the aspect ratio decreases from 4.54 to 2.22 as the raising of the point defect concentration, as shown in Figure 2b,e. The value of k decreases from 1.71 to 1.37 with the effect of the tensile strain ε x x a = 0.025 , as shown in Figure 2c,f. Therefore, when there is a defects loop at the interface of the α/α′ phase, the deformation of the Cr-enriched nanoparticle is weakened under the compressive or tensile strain.
According to the results of rafting in the superalloys, the elastic energy is dominant in the formation of the rafting morphology [11]. Therefore, the distribution of elastic energy can clearly show the influence of strain on particle deformation with the point defects. Figure 3 shows the dimensionless elastic energy density E e distribution along the direction of the particle elongation through the center of the α′ phase under different strains at 750 K. Points A and B show the peak values of the dimensionless elastic energy density E eA and E eB at the interface regions of the α/α′ phase, respectively.
When ε x x a = 0 , the E e with defects is higher than that without defects, and the peak value of the dimensionless elastic energy density increases from E eA = 1.4 × 10 3 to E eB = 2.7 × 10 3 at the interface region, as shown in Figure 3a. The result indicates that, when there is no normal strain, the E e at the interface region is elevated when there is a defects loop, even if the coherent eigenstrain between the α/α′ phase is reduced. However, the vacancy or interstitial atom produces the negative eigenstrain to increase the elastic strain. The large E e results in a higher concentration of Cr under ε x x a = 0 , as shown in Figure 4a,b.
Under the compressive strain ε x x a = 0.025 , when the concentration of point defects increases to 1.0 × 10–7, the peak value of the dimensionless elastic energy density of the interface regions also increases from E eA = 3.7 × 10 2 to E eB = 5.2 × 10 2 , as shown in Figure 3b. Compared with Figure 3a and Figure 3c, the E e in the α′ phase is higher than that without normal strain and lower than that with tensile strain. In addition, the E e in the matrix is higher than that of the α′ phase under the compressive strain ε x x a = 0.025 . For soft precipitates, the energy required for a sheet or a lathy shape deformation is lower than that of spherical precipitates, and the matrix will prevent them from contracting along a sheet or a lathy shape [9]. In the Fe–Cr alloy, the α′ phase is softer than the matrix, so the lamellar and lathy shapes are formed obviously under uniaxial compression strain, and the aspect ratio k of the α′ phase is greater than that of the tensile strain (see Figure 2).
Under tensile strain ε x x a = 0.025 , the peak values E eA and E eB at the interface of the α/α′ phase are similar with and without the defects loop, as shown in Figure 3c. The E e inside the particle is decreased as the concentration of point defects increases to c V i = c I i = 1 × 10 7 , which is contrary to that of ε x x a = 0 and ε x x a = 0.025 . Under tensile strain, the low E e inside the α′ phase weakens the deformation of the particles.

3.2. Concentration Evolution of the α′ Phase with the Effects of Point Defects, Strains and Temperature

The concentration evolution of Cr in the Cr-enriched α′ phase with the vacancy and interstitial atoms are studied under normal strain at 700 K and 750 K; the influences of temperature, strain, and vacancy on equilibrium concentration are analyzed. The concentration distributions of Cr in the α′ phase and vacancy at 750 K are shown in Figure 4. The concentration changes of Cr at 750 K under the effects of compressive and tensile strains are listed in Table 1. When there are no point defects, the Cr concentration in the α′ phase increases from c Cr 0 = 0.48 to c Cr c = 0.49 under the effect of compressive strain ε x x a = 0.025 , and to c Cr t = 0.58 under tensile strain ε x x a = 0.025 , as shown in Figure 4a. The concentration changes of Cr under the effect of compressive and tensile strain are Δ c Cr = c Cr c c Cr 0 = 0.01 and Δ c Cr = c Cr t c Cr 0 = 0.10 , respectively. The relative concentration differences, Δ c Cr / c Cr 0 , are 0.021 and 0.208, respectively. The results indicate that the effect of compressive strain on initial concentration clustering is smaller than that of tensile strain when there are no point defects.
When the initial point defect concentration is c V i = c I i = 1 × 10 7 , the concentration of Cr in the α′ phase increases from c Cr 0 = 0.55 to c Cr c = 0.64 and c Cr t = 0.73 under the effects of compressive strain ε x x a = 0.025 and tensile strain ε x x a = 0.025 , as shown in Figure 4b. The concentration changes of Cr in the α′ phase are Δ c Cr = 0.09 and Δ c Cr = 0.18 , and the relative concentration differences Δ c Cr / c Cr 0 are 0.164 and 0.327, respectively. Compared with no defects, the Δ c Cr increases by 0.08 under both ε x x a = 0.025 and ε x x a = 0.025 , and the relative concentration difference Δ c Cr / c Cr 0 increases by 0.143 and 0.119 under the effects of compression and tensile strains. The results show that the point defects promote the Cr atoms cluster, and the effect of compressive strain on accelerating the phase separation is enhanced with the point defects. However, the effect of strain is slight on the concentration change of point defects, as shown in Figure 4d.
A previous study of Fe–Cr alloy shows that the effect of point defects on the phase separation is reduced at low temperature [17]. The concentration distributions of Cr and vacancy at 700 K under different strains are displayed in Figure 5. The concentration changes of Cr with and without defect at 700 K under the effects of compressive and tensile strains are listed in Table 2, where, the peak values of cCr are defined as c Cr i 0 (i = c, t) for c V i = c I i = 0 and c Cr i 1 for c V i = c I i = 1 × 10 7 under the effect of compressive (i = c) ε x x a = 0.025 and tensile strains (i = t) ε x x a = 0.025 , respectively. When the point defects concentration changes from 0.0 to 1.0 × 10–7, the Cr concentration increases from c Cr c 0 = 0.77 to c Cr c 1 = 0.83 , and from c Cr t 0 = 0.76 to c Cr t 1 = 0.84 under the effects of compressive and tensile strains, as shown in Figure 5a,b. The concentration differences Δ c Cr i = c Cr i 1 c Cr i 0 , of Cr in the α′ phase for compressive and tensile strains are Δ c Cr c = 0.06 and Δ c Cr t = 0.08 , respectively.
With the increase in initial point defects concentration, the concentration difference of Cr in the α′ phase shows a very slight change. Therefore, point defects have little influence on the Cr concentration at low temperature (700 K), whether there is a normal strain or not. At 700 K, the effect of tensile strain on the concentration evolution is similar to that of compressive strain, while the compressive strain shows an obvious influence at 750 K (see Figure 4a,b). Under the effect of normal strain, the concentration change of vacancy is very small, as shown in Figure 5d. Similar to 750 K, the effect of the normal strain is very small for the concentration of point defects.
In order to further investigate the effect of point defects and strains at different temperatures, the equilibrium concentration of Cr, c Cr eq , in the α′ phase with temperature change is shown in Figure 6. It can be seen that the c Cr eq decreases with the increased temperature with or without point defects. Under the influence of external strains, c Cr eq is raised, as shown in Figure 6a,b. Zhu et al. studied the free energy change under the effects of external strain and found that the increase in strain makes the free energy decrease and leads to a higher Cr concentration c Cr eq [20]. When the point defects concentration changes from 0.0 to 1.0 × 10–7, the c Cr eq in the α′ phase shows an obvious decline under the effect of the compressive strain ε x x a = 0.025 , as shown in Figure 6b. In addition, the Cr concentration differences among the three strain states become small. Therefore, the defects can reduce the difference of c Cr eq , especially under normal strain.
In the vacancy diffusion mechanism, the atom migrates via the position exchange with vacancy [46]. To study the evolution of vacancy and its influence on the phase separation under strain, the morphology evolution of the Cr-enriched α′ phase and concentration distributions of Cr and vacancy under ε x x a = 0.025 are shown in Figure 7. In the stage of nucleation and growth, as time changes from t* = 0.001 to t* = 2, the concentration of Cr increases from 0.35 to 0.44 with point defects, and to 0.41 without point defects, as shown in Figure 7d. At the same time, the Cr atoms cluster and form a small sized α′ phase, as shown in Figure 7a–d. The peak of vacancy concentration increases from 1.8 × 10–7 to 1.7 × 10–4 and moves from the initial boundary position to the center, as shown in Figure 7g–j. The rate of the concentration change of Cr in the α′ phase and vacancy is calculated by Δ c k = ( c k t = 2 c k t = 0.001 ) / c k t = 0.001   ( k = Cr ,   V ) , Δ c Cr = 0.257 in the α′ phase, and Δ c V = 944.9 for vacancy. The change rate of vacancy concentration is faster than that of the Cr-enriched α′ phase. The early stage evolution of vacancy provides more position to facilitate the mobility of Cr atoms, which makes the nucleation and growth of Cr-enriched nanoscale particles faster. With the increase in time, the vacancy is moving from the center to the interface of the α/α′ phases (see Figure 7k,l), and the exchange process of the vacancy and the Cr atoms promotes the growth of the Cr-enriched α′ phase [17]. In general, the cluster and mobility of point defects can induce a high diffusion path to promote the initial phase separation [47]. With the growth of the particle, the elongation of a Cr-enriched nanoscale particle becomes obvious under the compressive strain, and is weakened with the formation of the defects loop, as shown in Figure 7e,f. Similarly, the evolution of the Cr-enriched α′ phase under tensile strain is the same as that under compression strain.

3.3. Splitting of the α′ Phase with Defects Loop under Normal Strain

The splitting of nanoparticles under the influence of strain and temperature, and the relationship between the initial particle size and particle splitting are studied. The main factors of splitting are revealed by studying the elastic energy and interfacial energy. Figure 8 shows the concentration distribution in 2D and 3D for the Cr-enriched α′ phase at 700 K, 720 K and 750 K with c V i = c I i = 1 × 10 7 in an initial circle region under ε x x a = 0.025 . At 700 K, the Cr concentration increases with time and then splits into two peaks, as shown in Figure 8a–d. It can be seen from the experiment that the particle splitting requires specific conditions [21], and that the phenomenon of particle splitting was also found during the growth and coarsening of γ′ precipitates in Ni-based superalloys [48,49], which is caused by the mismatch of the interface energy and the elastic energy. The particle splits into smaller ones to relax the elastic energy when its size reaches a critical value [48]. Therefore, the splitting is not a generic phenomenon. In the simulation of single particle splitting, the initial component of the particle is non-equilibrium, and with the growth of the particle, the particle splits [50]. When the concentration of the particle reaches a higher value or equilibrium, the particle itself does not continue to split, but the coarsening will continue [50]. Similarly, in our study, under a non-equilibrium initial concentration, the particle splits at low temperature and the splitting stops when the concentration reaches 0.54, as shown in Figure 8c. Therefore, the splitting is only a phenomenon under certain simulation conditions, and in our study it occurs at the initial phase separation.
At 720 K, the particle shows a concentration clustering and a tendency to split, as shown in Figure 8f,g, while the particle is merged again into one with an equilibrium concentration 0.862 in Figure 8h. Then the particle is elongated under the compressive strain. At 750 K, the concentration of Cr clusters and evolves into a particle with small size, then the particle grows up and elongates rapidly, as shown in Figure 8k–o. Because the atomic diffusion is faster at a high temperature, the lower interface energy and elastic energy in the center of particle will further inhibit the splitting of particle, as shown in Figure 8k,l, so that the Cr-enriched nanoparticle grows up and coarsens individually without splitting at 750 K.
In addition, when there is no point defect, the particle still splits at 700 K under a compressive strain ε x x a = 0.025 , as shown in Figure 9a–e. When the strain decreases to ε x x a = 0.019 , the split particles will be recombined, as shown in Figure 9h–j, which is similar to the state of 720 K (see Figure 8f–j) with a compressive strain ε x x a = 0.025 . Therefore, the particles splitting can be eliminated by increasing the temperature and decreasing the strain. The results show that the point defects are not a determining factor for the splitting of particle, and that the particle splitting is affected by the interface energy and the strain energy [50]. When the ratio of the interfacial energy to the elastic energy is increased, the splitting becomes difficult [24]. Therefore, the relationship between the interfacial energy and the elastic energy needs to be discussed. With the increase in the applied strain, the α′ phase splits into two or more particles during its growth. Therefore, the particle’s splitting occurs under high strain.
To reveal the α′ phase splitting with a defects loop under normal strain, the dimensionless gradient energy density E g , the dimensionless chemical free energy density E c , the dimensionless elastic energy density E e , and the Cr concentration cCr distributions along the y*-axis through the center of the α′ phase are shown in Figure 10. On the whole, the distribution of E c (see Figure 10b,f,j) is similar to that of cCr (see Figure 10a,e,i), while the distribution of E e inside of the α′ phase is contrary to that of cCr.
At 700 K, the Cr concentration shows two peaks at t* = 41, which means that the splitting of the α′ phase particle occurs. At the same time, the E g has a valley in the center position, as signed by the arrow in Figure 10d; while it is not matched with the increased E e in the center position, as shown by the arrows in Figure 10c, this energy that is different inside of the α′ phase causes the particle to split, as shown in Figure 8a–c. In the Ni-Al alloy with negative mismatch, the split of the γ′ precipitates results in the rapid increase in interface energy and a reduction in elastic strain energy, which is caused by the inverse precipitation of the γ phase in the center of the γ′ precipitates [48,50]. As the temperature increases, the change of E g and E e is gradually consistent in the center of particle at t* = 41; as shown in Figure 10g,h,k,l, the particles coarsen independently without splitting.
In addition to the temperature, the initial particle size also affects the splitting of the particle in the single-particle system [50]. Figure 11 shows the concentration distribution in 2D and 3D for the Cr-enriched α′ phase at 700 K and 750 K in two-particle systems with ε x x a = 0.025 and c V i = c I i = 1 × 10 7 , where the critical radius of the initial particle is R= 3.5 l at 700 K and 8 l at 750 K, respectively, as shown in Figure 11a,c. At 700 K, when the radius of the initial particle is less than R = 3.5 l, the particles coarsened independently without splitting under ε x x a = 0.025 , as shown in Figure 11b. When the temperature increases to 750 K, the particles split at a larger particle radius of 8 l, as shown in Figure 11c,d. Therefore, the critical radius of particle splitting increases with the increase of the temperature. The split morphology of the two particles at 750 K is similar to that of single particles at 700 K in Figure 8e. Therefore, the particle size and their distance affect the splitting behavior in the multi-particle system [50]. Generally, when the size of the initial non-equilibrium particles is larger, the particles are easy to split at the early stage of phase separation under a higher normal strain (Supplementary Materials).

4. Conclusions

The phase-field model coupled with defects and strain is developed to study the evolutions of point defects and α′ nanophase in the Fe-35at.% Cr alloys. The effects of the point defect and normal strain on the concentration evolution and the shape change of the nanoparticles are presented. It is shown that a coherent strain between the α/α′ phases is relaxed by forming a defects loop at the interface regions of the α/α′ phase, thus reducing the eigenstrain to weaken the deformation of the nanophase. The point defects enhance the effect of normal strain in accelerating the phase separation, which is obvious under the compressive strain at high temperature. Under the normal strain, the small initial radius of the Cr-enriched nanophase splits in the early-stage precipitation at 700 K, while the splitting disappears with the increased temperature and the decreased strains. In addition, the large initial size of the particles is more likely to split.
The equilibrium concentration of Cr in the α′ phase increases with the increased normal strain and the concentration of the point defects. However, under compression strain, the increased concentration of point defects can decrease the equilibrium concentration of Cr in the α′ phase. With the reconciliation of the point defects, the concentration differences become small for the strain-free, tensile, and compressive states. This study offers a unique perspective on the formation of the irradiation point defect loop in the phase separation alloy, and its influence on the evolution of nanoparticles under normal strain, and the relationship between point defects and normal strain are clarified.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nano13030456/s1. The ppt is a gif of the evolution of the Cr-enriched α′ phase and point defect under compressive strain.

Author Contributions

Formal Analysis and writing-original draft preparation, Z.Y.; Validation, S.S.; Visualization, P.S.; Data Curation, K.L.; Investigation, Q.Q.; Conceptualization, software, methodology, funding acquisition, writing-review, and editing, Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the [National Natural Science Foundation of China] grant number [52275342] and the [Fundamental Research Funds for the Central Universities] grant number [30921013107].

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Clouet, E.; Varvenne, C.; Jourdan, T. Elastic modeling of point-defects and their interaction. Comput. Mater. Sci. 2018, 147, 49–63. [Google Scholar] [CrossRef] [Green Version]
  2. Soisson, F.; Jourdan, T. Radiation-accelerated precipitation in Fe–Cr alloys. Acta Mater. 2016, 103, 870–881. [Google Scholar] [CrossRef]
  3. Grönhagen, K.; Ågren, J.; Odéna, M. Phase-field modelling of spinodal decomposition in TiAlN including the effect of metal vacancies. Scr. Mater. 2015, 95, 42–45. [Google Scholar] [CrossRef]
  4. Liang, L.Y.; Mei, Z.G.; Kim, Y.S.; Anitescu, M.; Yacout, A.M. Three-dimensional phase-field simulations of intragranular gas bubble evolution in irradiated U–Mo fuel. Comput. Mater. Sci. 2018, 145, 86–95. [Google Scholar] [CrossRef]
  5. Eshelby, J.D. The continuum theory of lattice defects. Solid State Phys. 1956, 3, 79–144. [Google Scholar] [CrossRef]
  6. Sugathan, S.; Bhattacharya, S. A phase-field study of elastic stress effects on phase separation in ternary alloys. Comput. Mater. Sci. 2020, 172, 109284. [Google Scholar] [CrossRef] [Green Version]
  7. Tian, J.Y.; Xu, P.; Liu, Q.B. Effects of stress-induced solid phase transformations on residual stress in laser cladding a Fe–Mn–Si–Cr–Ni alloy coating. Mater. Des. 2020, 193, 108824. [Google Scholar] [CrossRef]
  8. Apte, A.; Kochat, V.; Rajak, P.; Krishnamoorthy, A.; Manimunda, P.; Hachtel, J.A.; Idrobo, J.C.; Amanulla, S.A.S.; Vashishta, P.; Nakano, A.; et al. Structural phase transformation in strained monolayer MoWSe2 alloy. ACS Nano 2018, 12, 3468–3476. [Google Scholar] [CrossRef]
  9. Fratzl, P.; Penrose, O.; Lebowitz, J.L. Modeling of phase separation in alloys with coherent elastic misfit. J. Stat. Phys. 1999, 95, 1429–1503. [Google Scholar] [CrossRef] [Green Version]
  10. Tien, J.K.; Copley, S.M. The effect of orientation and sense of applied uniaxial stress on the morphology of coherent gamma prime precipitates in stress annealed nickel-base superalloy crystals. Metall. Trans. 1971, 2, 543–553. [Google Scholar] [CrossRef]
  11. Wang, D.; Li, Y.S.; Shi, S.J.; Tong, X.W.; Yan, Z.W. Phase-field simulation of γ’ precipitates rafting and creep property of Co-base superalloys. Mater. Des. 2020, 196, 109077. [Google Scholar] [CrossRef]
  12. Sakaguchia, M.; Okazaki, M. Distinctive role of plastic and creep strain in directional coarsening of a Ni-base single crystal superalloy. Mater. Sci. Eng. A 2018, 710, 121–128. [Google Scholar] [CrossRef]
  13. Prikhodko, S.V.; Ardell, A.J. Coarsening of γ′ in Ni–Al alloys aged under uniaxial compression: I. Early-stage kinetics. Acta Mater. 2003, 51, 5001–5012. [Google Scholar] [CrossRef]
  14. Ardell, A.J.; Prikhodko, S.V. Coarsening of γ′ in Ni–Al alloys aged under uniaxial compression: II. Diffusion under stress and retardation of coarsening kinetics. Acta Mater. 2003, 51, 5013–5019. [Google Scholar] [CrossRef]
  15. Tong, X.W.; Li, Y.S.; Yan, Z.W.; Wang, D. Phase-field simulation of effects of normal strain on the morphology and kinetics evolution of nanoscale phase. J. Mater. Res. Technol. 2020, 9, 2063–2071. [Google Scholar] [CrossRef]
  16. Maqbool, S.; Li, Y.S.; Muhammad, S.; Yan, Z.W. Phase-feld simulation of Cu enriched nanoparticles with variation of defects migration energy under neutron irradiation. Model. Simul. Mater. Sci. Eng. 2021, 29, 085011. [Google Scholar] [CrossRef]
  17. Yan, Z.W.; Shi, S.J.; Li, Y.S.; Chen, J.; Maqbool, S. Vacancy and interstitial atom evolution with the separation of the nanoscale phase in Fe–Cr alloys: Phase-field simulations. Phys Chem. Chem. Phys. 2020, 22, 3611–3619. [Google Scholar] [CrossRef] [PubMed]
  18. Chen, L.; Huang, X.F.; He, K.; Wang, H. Precipitation behavior and age hardening effect of the precipitates in a Fe-13Cr-4Al-2Mo-1.2Nb alloy. Mater. Charact. 2021, 173, 110918. [Google Scholar] [CrossRef]
  19. Kim, H.; Subramanian, G.O.; Kim, C.; Jang, H.; Jang, C. 400 °C aging embrittlement of FeCrAl alloys: Microstructure and fracture behavior. Mater. Sci. Eng. A 2019, 743, 159–167. [Google Scholar] [CrossRef]
  20. Zhu, L.H.; Li, Y.S.; Liu, C.W.; Chen, S.; Shi, S.J.; Jin, S.S. Effect of applied strain on phase separation of Fe–28 at. % Cr alloy: 3D phase-field simulation. Model. Simul. Mater. Sci. Eng. 2018, 26, 035015. [Google Scholar] [CrossRef]
  21. Doi, M.; Miyazaki, T.; Wakatsuki, T. The effect of elastic interaction energy on the morphology of γ′ precipitates in nickel-based alloys. Mater. Sci. Eng. 1984, 67, 247–253. [Google Scholar] [CrossRef]
  22. Yang, M.; Wei, H.; Zhang, J.; Zhao, Y.; Jin, T.; Liu, L.; Sun, X.F. Phase-field study on effects of antiphase domain and elastic energy on evolution of γ′ precipitates in nickel-based superalloys. Comput. Mater. Sci. 2017, 129, 211–219. [Google Scholar] [CrossRef]
  23. Lee, J.K. A study on coherency strain and precipitate morphologyvia a discrete atom method. Metall. Mater. Trans. A 1996, 27, 1449–1459. [Google Scholar] [CrossRef]
  24. Li, D.Y.; Chen, L.Q. Shape evolution and splitting of coherent particles under applied stresses. Acta Mater. 1999, 47, 247–257. [Google Scholar] [CrossRef]
  25. Khachaturyan, A.G. Theory of Structural Transformations in Solids; John Wiley and Sons: New York, NY, USA, 1983. [Google Scholar]
  26. Smirnova, D.E.; Kuksin, A.Y.; Starikov, S.V. Investigation of point defects diffusion in bcc uranium and U–Mo alloys. J. Nucl. Mater. 2015, 45, 304–311. [Google Scholar] [CrossRef]
  27. Hu, S.Y.; Chen, L.Q. A phase-field model for evolving microstructures with strong elastic inhomogeneity. Acta Mater. 2001, 49, 1879–1890. [Google Scholar] [CrossRef]
  28. Li, Y.S.; Li, S.X.; Zhang, T.Y. Effect of dislocations on spinodal decomposition in Fe–Cr alloys. J Nucl Mater. 2009, 395, 120–130. [Google Scholar] [CrossRef]
  29. Wang, J.J.; Yang, T.N.; Zorn, J.A.; Wang, E.; Irwin, J.L.; Lindemann, S.; Rzchowski, M.S.; Hu, J.M.; Eom, C.B.; Chen, L.Q. Strain anisotropy and magnetic domain structures in multiferroic heterostructures: High-throughput finite-element and phase-field studies. Acta. Mater. 2019, 176, 73–83. [Google Scholar] [CrossRef]
  30. Cahn, J.W. Phase separation by spinodal decomposition in isotropic system. J. Chem. Phys. 1965, 42, 93–99. [Google Scholar] [CrossRef]
  31. Andersson, J.O.; Sundman, B. Thermodynamic properties of the Cr–Fe system. Calphad 1987, 11, 83–92. [Google Scholar] [CrossRef]
  32. Wang, Y.Y.; Ding, J.H.; Huang, S.S.; Zhao, J.J.; Wang, Y.Z. Mesoscale modeling of irradiation damage evolution in bcc iron and vanadium: A comparative study. Fusion Eng. Des. 2018, 137, 303–311. [Google Scholar] [CrossRef]
  33. Dinsdale, T. SGTE data for pure elements. Calphad 1991, 15, 317–425. [Google Scholar] [CrossRef]
  34. Millett, P.C.; Azab, A.E.; Rokkam, S.; Tonks, M.; Wolf, D. Phase-field simulation of irradiated metals Part I: Void kinetics. Comput. Mater. Sci. 2011, 50, 949–959. [Google Scholar] [CrossRef]
  35. Jiang, Y.B.; Liu, W.B.; Li, W.J.; Sun, Z.Y.; Xin, Y.; Chen, P.H.; Yun, D. Phase-field simulation of the interaction between intergranular voids and grain boundaries during radiation in UO2. Comput. Mater. Sci. 2021, 188, 110176. [Google Scholar] [CrossRef]
  36. Mohanty, R.R.; Guyer, J.E.; Sohn, Y. Diffusion under temperature gradient: A phase-field model study. J. Appl. Phys. 2009, 106, 034912. [Google Scholar] [CrossRef]
  37. Okita, T.; Sekimura, N.; Garner, F.A. Effects of dpa rate on swelling in neutron-irradiated Fe–Cr and Fe–Cr–Mo alloys. J. Nucl. Mater. 2011, 417, 944–948. [Google Scholar] [CrossRef]
  38. Gao, Y.P.; Zhang, Y.F.; Schwen, D.; Jiang, C.; Sun, C.; Gan, J. Formation and self-organization of void superlattices under irradiation: A phase field study. Materialia 2018, 1, 78–88. [Google Scholar] [CrossRef]
  39. Xiong, W.; Selleby, M.; Chen, Q.; Odqvist, J.; Du, Y. Phase equilibria and thermodynamic properties in the Fe–Cr system. Crit. Rev. Solid State Mater. Sci. 2010, 35, 125–152. [Google Scholar] [CrossRef]
  40. Barkar, T.; Höglund, L.; Odqvist, J.; Ågren, J. Effect of concentration dependent gradient energy coefficient on spinodal decomposition in the Fe–Cr system. Comput. Mater. Sci. 2018, 143, 446–453. [Google Scholar] [CrossRef]
  41. Zhu, L.H.; Li, Y.S.; Shi, S.J.; Yan, Z.W.; Chen, J.; Maqbool, S. Morphology and Kinetics Evolution of Nanoscale Phase in Fe–Cr Alloys under External Strain. Nanomaterials 2019, 9, 294. [Google Scholar] [CrossRef]
  42. Cottua, M.; Appolaire, B.; Finel, A.; Bouar, Y.L. Microstructure evolution under [110] creep in Ni-base superalloys. Acta Mater. 2021, 212, 116851. [Google Scholar] [CrossRef]
  43. Mughrabi, H. The importance of sign and magnitude of γ/γ′ lattice misfit in superalloys with special reference to the new γ′-hardened cobalt-base superalloys. Acta Mater. 2014, 81, 21–29. [Google Scholar] [CrossRef]
  44. Wang, J.C.; Osawa, M.; Yokokawa, T.; Harada, H.; Enomotoc, M. Modeling the microstructural evolution of Ni-base superalloys by phase field method combined with CALPHAD and CVM. Comput. Mater. Sci. 2007, 39, 871–879. [Google Scholar] [CrossRef]
  45. Li, Y.S.; Zhu, H.; Zhang, L.; Cheng, X.L. Phase decomposition and morphology characteristic in thermal aging Fe–Cr alloys under applied strain: A phase-field simulation. J. Nucl. Mater. 2012, 429, 13–18. [Google Scholar] [CrossRef]
  46. Anento, N.; Serra, A.; Osetsky, Y. Effect of nickel on point defects diffusion in Fe–Ni alloys. Acta Mater. 2017, 132, 367–373. [Google Scholar] [CrossRef]
  47. Rachbauer, R.; Massl, S.; Stergar, E.; Holec, D.; Kiener, D.; Keckes, J.; Patscheider, J.; Stiefel, M.; Leitner, H.; Mayrhofer, P.H. Decomposition pathways in age hardening of Ti–Al–N films. J. Appl. Phys. 2011, 110, 023515. [Google Scholar] [CrossRef] [Green Version]
  48. Banerjee, D.; Banerjee, R.; Wang, Y. Formation of split patterns of γ′ precipitates in Ni-Al via particle aggregation. Scr. Mater. 1999, 41, 1023–1030. [Google Scholar] [CrossRef]
  49. Yoo, Y.S.; Yoon, D.Y.; Henry, M.F. The effect of elastic misfit strain on the morphological evolution of γ′-precipitates in a model Ni-base superalloy. Met. Mater. 1995, 1, 47–61. [Google Scholar] [CrossRef]
  50. Liu, L.; Chen, Z.; Wang, Y.X. Elastic strain energy induced split during precipitation in alloys. J. Alloy Compd. 2016, 661, 349–356. [Google Scholar] [CrossRef]
Figure 1. Morphology of the Cr-enriched α′ nanophase (af) and vacancy (gi) at 750 K for t* = 50. (ac) c V i = c I i = 0 , (di) c V i = c I i = 1 × 10 7 . (a,d,g) ε x x a = 0 , (b,e,h) ε x x a = 0.025 , (c,f,i) ε x x a = 0.025 .
Figure 1. Morphology of the Cr-enriched α′ nanophase (af) and vacancy (gi) at 750 K for t* = 50. (ac) c V i = c I i = 0 , (di) c V i = c I i = 1 × 10 7 . (a,d,g) ε x x a = 0 , (b,e,h) ε x x a = 0.025 , (c,f,i) ε x x a = 0.025 .
Nanomaterials 13 00456 g001
Figure 2. The concentration distributions of Cr in the α′ phase along x* = 32 and y* = 32 with and without normal strain at 750 K for t* = 50. (ac) c V i = c I i = 0 , and (di) c V i = c I i = 1 × 10 7 .
Figure 2. The concentration distributions of Cr in the α′ phase along x* = 32 and y* = 32 with and without normal strain at 750 K for t* = 50. (ac) c V i = c I i = 0 , and (di) c V i = c I i = 1 × 10 7 .
Nanomaterials 13 00456 g002
Figure 3. The dimensionless elastic energy density distribution through the center of α′ phase at t* = 50 with and without point defects. (a) ε x x a = 0 , (b) ε x x a = 0.025 , and (c) ε x x a = 0.025 .
Figure 3. The dimensionless elastic energy density distribution through the center of α′ phase at t* = 50 with and without point defects. (a) ε x x a = 0 , (b) ε x x a = 0.025 , and (c) ε x x a = 0.025 .
Nanomaterials 13 00456 g003aNanomaterials 13 00456 g003b
Figure 4. The concentration distributions of Cr and vacancy across the Cr-enriched α′ phase and vacancy at 750 K under different strains at t* = 3. (a) c V i = c I i = 0 , (bd) c V i = c I i = 1 × 10 7 , and (d) local magnification of the concentration in Figure 4c.
Figure 4. The concentration distributions of Cr and vacancy across the Cr-enriched α′ phase and vacancy at 750 K under different strains at t* = 3. (a) c V i = c I i = 0 , (bd) c V i = c I i = 1 × 10 7 , and (d) local magnification of the concentration in Figure 4c.
Nanomaterials 13 00456 g004
Figure 5. The concentration distributions of Cr and vacancy across the Cr-enriched α′ phase, and vacancy at 700 K under different strains at t* = 71. (a) c V i = c I i = 0 , (bd) c V i = c I i = 1 × 10 7 , and (d) local magnification of the concentration in Figure 5c.
Figure 5. The concentration distributions of Cr and vacancy across the Cr-enriched α′ phase, and vacancy at 700 K under different strains at t* = 71. (a) c V i = c I i = 0 , (bd) c V i = c I i = 1 × 10 7 , and (d) local magnification of the concentration in Figure 5c.
Nanomaterials 13 00456 g005
Figure 6. The equilibrium concentration of Cr in the α′ phase with temperature change at t* = 791 under different strains. (a) c V i = c I i = 0 , (b) c V i = c I i = 1 × 10 7 .
Figure 6. The equilibrium concentration of Cr in the α′ phase with temperature change at t* = 791 under different strains. (a) c V i = c I i = 0 , (b) c V i = c I i = 1 × 10 7 .
Nanomaterials 13 00456 g006
Figure 7. The morphology evolution of the Cr-enriched α′ phase and concentration distributions of Cr in the α′ phase with and without point defects (af), the morphology evolution and concentration distributions of vacancy (gl) with c V i = c I i = 1 × 10 7 , T = 750 K, ε x x a = 0.025 .
Figure 7. The morphology evolution of the Cr-enriched α′ phase and concentration distributions of Cr in the α′ phase with and without point defects (af), the morphology evolution and concentration distributions of vacancy (gl) with c V i = c I i = 1 × 10 7 , T = 750 K, ε x x a = 0.025 .
Nanomaterials 13 00456 g007
Figure 8. The concentration surface evolution of the Cr-enriched α′ phase from the initial circle regions with 35 at.% Cr and c V i = c I i = 1 × 10 7 under ε x x a = 0.025 , (ae) 700 K, (fj) 720 K and (ko) 750 K.
Figure 8. The concentration surface evolution of the Cr-enriched α′ phase from the initial circle regions with 35 at.% Cr and c V i = c I i = 1 × 10 7 under ε x x a = 0.025 , (ae) 700 K, (fj) 720 K and (ko) 750 K.
Nanomaterials 13 00456 g008aNanomaterials 13 00456 g008b
Figure 9. The concentration surface evolution of the Cr-enriched α′ phase from the initial circle regions with 35 at.% Cr at 700 K. (ae) ε x x a = 0.025 ,   c V i = c I i = 0 , and (fj) ε x x a = 0.019 ,   c V i = c I i = 1 × 10 7 .
Figure 9. The concentration surface evolution of the Cr-enriched α′ phase from the initial circle regions with 35 at.% Cr at 700 K. (ae) ε x x a = 0.025 ,   c V i = c I i = 0 , and (fj) ε x x a = 0.019 ,   c V i = c I i = 1 × 10 7 .
Nanomaterials 13 00456 g009
Figure 10. The dimensionless energy density and concentration distributions of Cr across of the Cr-enriched α′ phase with c V i = c I i = 1 × 10 7 at 700 K (ad), 720 K (eh) and 750 K (il) under ε x x a = 0.025 . (a,e,i) concentration of Cr, (b,f,j) E c , (c,g,k) E e , (d,h,l) E g .
Figure 10. The dimensionless energy density and concentration distributions of Cr across of the Cr-enriched α′ phase with c V i = c I i = 1 × 10 7 at 700 K (ad), 720 K (eh) and 750 K (il) under ε x x a = 0.025 . (a,e,i) concentration of Cr, (b,f,j) E c , (c,g,k) E e , (d,h,l) E g .
Nanomaterials 13 00456 g010
Figure 11. The concentration surface evolution of the Cr-enriched α′ phase from the initial circle regions with 35 at.% Cr under ε x x a = 0.025 and c V i = c I i = 1 × 10 × 7 . (a,b) 700 K, R = 3.5 l, (c,d) 750 K, R = 8 l.
Figure 11. The concentration surface evolution of the Cr-enriched α′ phase from the initial circle regions with 35 at.% Cr under ε x x a = 0.025 and c V i = c I i = 1 × 10 × 7 . (a,b) 700 K, R = 3.5 l, (c,d) 750 K, R = 8 l.
Nanomaterials 13 00456 g011
Table 1. The concentration changes of Cr at 750 K under the effects of compressive and tensile strains.
Table 1. The concentration changes of Cr at 750 K under the effects of compressive and tensile strains.
c V i   =   c I i   =   0 c V i   =   c I i   =   1 × 10 7
ε x x a 0−0.0250.0250−0.0250.025
c Cr c Cr 0 = 0.48 c Cr c = 0.49 c Cr t = 0.58 c Cr 0 = 0.55 c Cr c = 0.64 c Cr t = 0.73
Δ c Cr = c Cr i c Cr 0   , ( i = c , t ) -0.010.10-0.090.18
Δ c Cr / c Cr 0 -0.0210.208-0.1640.327
Table 2. The concentration changes of Cr with and without defect at 700 K under the effects of compressive and tensile strains.
Table 2. The concentration changes of Cr with and without defect at 700 K under the effects of compressive and tensile strains.
c V i   =   c I i   =   0 c V i   =   c I i   =   1 × 10 7
ε x x a −0.0250.025−0.0250.025
c Cr c Cr c 0 = 0.77 c Cr t 0 = 0.76 c Cr c 1 = 0.83 c Cr t 1 = 0.84
Δ c Cr i = c Cr i 1 c Cr i 0   , ( i = c , t ) Δ c Cr c = 0.06 Δ c Cr t = 0.08
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yan, Z.; Shi, S.; Sang, P.; Li, K.; Qin, Q.; Li, Y. Elastic Strain Relaxation of Phase Boundary of α′ Nanoscale Phase Mediated via the Point Defects Loop under Normal Strain. Nanomaterials 2023, 13, 456. https://doi.org/10.3390/nano13030456

AMA Style

Yan Z, Shi S, Sang P, Li K, Qin Q, Li Y. Elastic Strain Relaxation of Phase Boundary of α′ Nanoscale Phase Mediated via the Point Defects Loop under Normal Strain. Nanomaterials. 2023; 13(3):456. https://doi.org/10.3390/nano13030456

Chicago/Turabian Style

Yan, Zhengwei, Shujing Shi, Peng Sang, Kaiyue Li, Qingqing Qin, and Yongsheng Li. 2023. "Elastic Strain Relaxation of Phase Boundary of α′ Nanoscale Phase Mediated via the Point Defects Loop under Normal Strain" Nanomaterials 13, no. 3: 456. https://doi.org/10.3390/nano13030456

APA Style

Yan, Z., Shi, S., Sang, P., Li, K., Qin, Q., & Li, Y. (2023). Elastic Strain Relaxation of Phase Boundary of α′ Nanoscale Phase Mediated via the Point Defects Loop under Normal Strain. Nanomaterials, 13(3), 456. https://doi.org/10.3390/nano13030456

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop