Values of lattice, diffusion, thermodynamic and elastic constants.
Abstract
This chapter is focused on the application of the phase-field method to the analysis of phase decomposition during the isothermal aging of alloys. The phase-field method is based on a numerical solution of either the nonlinear Cahn-Hilliard equation or the Cahn-Allen equation. These partial differential equations can be solved using the finite difference method among other numerical methods. The phase-field method has been applied to analyze different types of phase transformations in alloys, such as phase decomposition, precipitation, recrystallization, grain growth, solidification of pure metals and alloys, martensitic transformation, ordering reactions, and so on. One of the main advantages of phase-field method is that this method permits to follow the microstructure evolution in two or three dimensions as the time of phase transformations progresses. Thus, the morphology, size, and size distribution could be determined to follow their corresponding growth kinetics. Additionally, the evolution of chemical composition can also be followed during the phase transformations. Furthermore, both Allen-Cahn and Cahn-Hilliard equations can be solved simultaneously to analyze the presence of ordered phases or magnetic domains in alloys.
Keywords
- phase field method
- phase decomposition
- alloys
- aging
1. Introduction
In the field of materials science, it is important to analyze different moving free boundary problems in order to understand its effect on the phase transformations that may occur in materials. For instance, the solidification and grain growth, and martensitic transformation are diffusion-controlled and diffusion-less phase transformations, respectively, with this characteristic. One way to overcome this need is the use of the diffuse interface model [1]. Furthermore, the application of the phase-field method to this sort of problem permits the use of the order parameter and phase-field variable which takes into account the composition gradient energy present in a diffuse interface model. For instance [2], the order parameter
The phase-field method is based on the equations proposed by Cahn-Hilliard [1], Allen-Cahn [2], or Ginzburg-Landau [3]:
The composition gradient energy can be defined, for instance, for the field variable, concentration
One of the main advantages of phase-field method is that this method permits to follow the microstructure evolution in two or three dimensions as the time of phase transformations progresses. Thus, the morphology, size, and size distribution could be determined to follow their corresponding growth kinetics. Additionally, the evolution of chemical composition can also be followed during the phase transformations.
To solve either of the partial differential equations, Eqs (1) or (2), several numerical methods have been used such as finite difference method, difference volume method, and finite element method [5]. The use of explicit finite difference method is simple and good alternative to solve this type of differential equations. For instance, the finite difference method can be used to solve a simple partial differential equation such as the simplified one-dimension equation of the second Fick’s law:
This chapter is mainly focused on the application to the phase decomposition by the spinodal decomposition mechanism during the isothermal aging of hypothetical binary alloys using the nonlinear Cahn-Hilliard equation [7]. The effect of main parameters such as the atomic mobility of alloy and elastic-strain energy on the microstructure evolution and growth kinetics is analyzed. To conclude, the application of phase-field method is extended to the analysis of spinodal decomposition during isothermal aging of real binary and ternary alloy systems, such as Cu-Ni, Fe-Cr, and Cu-Ni-Fe. A comparison of simulated results with experimental ones is also included.
2. Phase decomposition in alloys
The formation of phases in alloys usually takes place by nucleation mechanism, growth mechanism, or spinodal decomposition mechanism, which is followed by the coarsening of phases in alloy systems. These three processes can be analyzed using the phase-field method, and their results can also be compared with the fundamental theories of phase transformations such as the Cahn-Hilliard spinodal decomposition theory [4] and the Lifshitz-Slyozov-Wagner (LSW) diffusion-controlled coarsening theory [8]. The phase decomposition that takes place by the spinodal decomposition mechanism is distinguished from the phase separation that occurs by nucleation and growth mechanism by the formation of an initial composition modulation, which shows an increase in the modulation amplitude with time. In contrast, the phase formation by nucleation and growth predicts that the formed phase has a composition very close to the equilibrium one from the start to the finish of phase transformation [9]. Besides, the spinodal decomposition is usually associated with the presence of a miscibility gap in the equilibrium phase diagram, as shown in Figure 1. The miscibility gap is the equilibrium line and there is only one α phase for compositions and temperatures above this line, whereas a mixture of two phases, α1 and α2, is present inside the miscibility gap. This figure also shows the existence of the chemical spinodal located within the miscibility gap. A supersaturated αsss phase is expected to decompose spinodally into a mixture of A-rich α1 and B-rich α2 phases for an alloy composition after heating at a temperature higher than that of the miscibility gap and then quenching and heating or aging at a temperature lower than the chemical spinodal. The miscibility gap in Figure 1 is usually related to the plot of free energy versus composition shown in Figure 2. This figure shows the free energy curve shape changes as the temperature decreases. This type of curve is known as the spinodal curve, and it indicates that any alloy composition is in unstable state and thus it is expected to decompose into a mixture of A-rich α1 and B-rich α2 phases. The minimum and saddle points at each temperature of the spinodal curve correspond to the equilibrium and chemical spinodal shown in Figure 1.
The Cahn-Hilliard theory of spinodal decomposition [4] was developed by the modification of the second Fick’s law equation in order to allow only the growth of the modulation amplitude in composition with a wavelength larger than a critical value. Furthermore, the nonlinear Cahn-Hilliard equation, used in the simulations of this work, has its origin in this theory.
3. Simulation of hypothetic binary alloys
In the numerical simulation of the phase decomposition for a hypothetic A-B binary alloy, the nonlinear Cahn-Hilliard equation [7] can be rewritten as follows:
The local chemical free energy
This equation represents a spinodal curve similar to those shown in Figure 2. In the first case of simulation, the mobility
In the second example, the mobility
The numerical solution of the former partial differential equation conducted to the following concentration profiles (Figure 5). The same characteristics described in the previous example are also observed for this case. However, the amplitude of the composition modulation for this case increases faster with time than that for the former case. The microstructure evolution for this case is different from the one shown for the previous case (Figure 6). That is, the decomposed phases form a lamellar structure instead of the irregular and interconnected morphology of the previous case. This may be attributed to the variable mobility of the decomposed phases [4].
To complete the numerical simulation of hypothetic binary alloys, the next case includes the presence of an isotropic elastic-strain energy
Figure 7 illustrates the concentration profiles for this case and the same characteristics, observed in the others, are also present. That is, the modulation amplitude increases with time. Besides, a mixture of A-rich α1 and B-rich α2 phases is formed after aging. The microstructure evolution is similar to that of the first case for short times. That is, the morphology of the decomposed phases is irregular and interconnected, which is known as percolated structure [4] (Figure 8). Nevertheless, one of the decomposed phases takes a cuboid shape for the longest times, which is attributed to the isotropic elastic-stain energy [9].
4. Simulation in real binary alloys
In the next part of this chapter, the numerical simulation of the phased decompositions of Cu-Ni and Fe-Cr alloys after aging at different temperatures for different times is shown. Simulated results are compared to the experimental ones. To begin, the Cu-Ni alloys are widely used in different industrial applications. The equilibrium phase diagram is shown in Figure 9. This diagram has a miscibility gap located at temperatures lower than 350°C [11]. Thus, a supersaturated solid solution, formed by heating above 350°C and then quenching, is expected to decompose spinodally into a mixture of Cu-rich and Ni-rich phases after aging at a temperature lower than 350°C. Nevertheless, the growth kinetics of spinodal decomposition is very slow due to the low atomic diffusivity at these temperatures [7]. Thus, the application of the phase-field method to analyze the spinodal decomposition seems to be a good alternative because of the slow kinetics.
The nonlinear Cahn-Hilliard equation, Eq. (7), was solved to analyze the phase decomposition in these alloys. The local energy
The interdiffusion coefficient
The gradient energy coefficient
The thermodynamic, diffusion, crystal lattice, and elastic parameters for the microstructure simulation were obtained from the literature [12–15] and these are shown in Table 1. The effect of coherency elastic-strain energy was considered to be present during the phase decomposition of Cu-Ni alloys in spite of the similar lattice parameters of copper and nickel [13]. This elastic-strain energy was introduced into Eq. (7), according to the simple definition proposed by Hilliard [4]:
Parameter | Ni-Cu alloys | |
---|---|---|
Crystal lattice parameter |
0.360 | |
0.0016 | ||
(8366.0 + 2.802T) + (−4359.6 + 1.812T)( |
||
Diffusion coefficient |
Cu 1.5–2.3 exp (−230,000–260,000 J mol−1)/RT | |
Ni 17–35 exp (−270,000–300,000 J mol−1)/RT | ||
Cu/Ni [15] |
The elastic constants,
Considering the elastic strain energy,
The microstructural simulation was carried out using the finite difference method with 101 × 101 points square grid with a mesh size of 0.25 nm and a time-step size of 10 s. The simulations were performed for the Ni-30at.% Cu alloy at temperatures between 250 and 322°C for different times.
Figure 10 shows the numerically calculated plots of Cu concentration versus distance for the Ni-30at.%Cu alloy solution treated (0 h) and aged at 300°C for different times. There is an increase in the modulation amplitude with aging time. The increase in amplitude at this temperature confirms that the phase decomposition occurs spinodally in this alloy. The long simulated aging times also confirm that the growth kinetics of phase decomposition is very slow in this alloy system.
The simulated microstructures of the Ni-30at.%Cu alloy aged at 300°C for 0, 222, 278, 361, and 444 h are shown in Figure 11 (a–e), respectively. The black and white regions correspond to the Cu-rich and Ni-rich phases, respectively. It can be seen that the morphology of the decomposed phases is irregular and interconnected. The volume fraction of the Ni-rich phase increased with aging time.
The experimental Ne-gas field ion microscopy (FIM) images of the Ni-30at.%Cu alloy aged at 300°C for 500 h is shown in Figure 12. There is a good agreement between the calculated and experimental morphologies of the decomposed phases.
In the next paragraphs, a second example about the numerical simulation of real alloys is presented. The Fe-Cr alloys are a very important alloy system since this is used as the basis for different industrial alloys, such as the family of stainless steels. The Fe-Cr equilibrium phase diagram [11] also shows a miscibility gap found at temperatures lower than 500°C (Figure 13). Thus, the phase decomposition of the supersaturated solid solution into a mixture of Cr-rich and Fe-rich phases is also expected as a result of aging at temperatures lower than 500 °C.
The phase decomposition simulation was based on a numerical solution of the nonlinear Cahn-Hilliard equation, Eq. (7). The formulation for this case is very similar to that described previously for the numerical simulation of the phase decomposition in Ni-Cu alloys.
The crystal lattice, thermodynamic, diffusion, and elastic constants for the microstructural simulation in Fe-Cr alloys were taken from references [13–16] and these parameters are shown in Table 2. The simulation of phase decomposition was pursued using the explicit finite difference method with 101 × 101 and 201 × 201 points square grids with a mesh size of 0.1 and 0.25 nm and a time-step size up to 10 s. The numerical simulation was performed for the Fe-40at.%Cr alloy aged at 470°C for times from 0 to 1000 h. It is important to mention that the initial composition modulation corresponding to the solution-treated sample was calculated using a random number generator.
Parameter | Values |
---|---|
Crystal lattice parameter |
0.2866 |
0.00614 | |
(18600.0 + 0.1T) | |
Diffusion coefficient |
|
Fe Cr |
The plots of Cr concentration versus distance, concentration profiles, for the Fe-40at.%Cr alloy aged at 470°C for different times are shown in Figure 14. These concentration profiles indicate clearly that the supersaturated solid solution decomposed spinodally into a mix of Cr-rich and Fe-rich phases since the modulation amplitude increases as the aging time increases.
Figure 15 shows the simulated microstructural evolution of the phase decomposition in the Fe-40-at.%Cr alloy aged at 470°C for times from 10 to 750 h. The white and gray zones represent the Fe-rich and Cr-rich phases, respectively. It can be observed and irregular and interconnected morphology of the decomposed phases in the alloy aged for times up to 10 h. This morphological characteristic is known as percolated structure, and it has been commonly observed to occur during the early stages of aging in the spinodally decomposed alloys. The HR-TEM micrographs of this alloy aged at 470°C for 250 h shows clearly the presence of spheres corresponding to the Cr-rich phases imbedded in the ferrite phase matrix (Figure 16). The decomposed phases present a coherent interface. This shape of decomposed Cr-rich phase is in good agreement with the simulated microstructure (Figure 15(e) and (f)) [17].
5. Simulation in real ternary alloys
The numerical simulation of phase decomposition of ternary alloys can be also conducted using the phase-field method. The equilibrium Cu-Ni-Fe phase diagram [18] is shown in Figure 17 at different temperatures. The presence of a miscibility gap is evident and thus the phase decomposition can also be simulated by the phase-field method.
The Cahn-Hilliard nonlinear equation for a multicomponent system with a constant mobility can be used for the present simulation, Eq. (7).
The local free energy
Parameter | Cu-Ni-Fe alloys | |
---|---|---|
(9534.49+2.839T) + (−424.255−0.629T)( (48206.0 − 8.446T) + (−5918.0 + 5.017T)( (−18298.8 + 5.149T) + (14313.6 − 7.659T)( −35982.0 − 12.0T |
||
−8.65T−22.64TlnT−3.13×10−3 T2−7023.9 93.23T−12.54TlnT+1.23×10−3 T2−532.3 39.0T−26.61TlnT+1.23×10−3 T2−4154.5 |
||
Diffusion coefficient |
||
Crystal lattice parameter |
0.360 | |
0.0016 | ||
Cu/Ni [15] |
c44= 7.54×1010 |
c11= 24.65 × 1010 c12= 14.73 × 1010 c44 = 12.47 × 1010 |
The atomic mobility
The atomic mobility
The gradient energy coefficient
Thus, the gradient energy coefficient
The effect of coherency elastic-strain energy was introduced into Eq. (7), according to the simple definition proposed by Hilliard [4]. It was shown in Eq. (17).
The elastic constant
The elastic constants,
Considering the elastic strain energy,
Equations (28) and (29) were solved numerically using the explicit finite difference method with 101 × 101 points square grid with a mesh size of 0.36 nm and a time-step size up to 10 s. The simulations were performed for the Cu-46at.%Ni-4at.%Fe alloy at temperature of 400°C for times from 0 to 200 h. This composition was selected for comparison of the morphology and kinetics of the phase decomposition. It is important to mention that the initial composition modulation corresponding to the solution treated sample was calculated using a random-number generator as proposed in reference [19].
Figure 18 shows the calculated concentration profiles of the Cu-46at.%Ni-4at.%Fe alloy aged at 400°C for times from 0 to 200 h. An increase in the amplitude of the composition modulation with aging time can be noticed in both cases. This fact has been associated with the phase decomposition via the spinodal decomposition mechanism [4]. This behavior is also in good agreement with the experimental evidence reported in the literature [18] for the aging of Cu-Ni-Fe alloys. The calculated Cu and Fe concentration profiles are shown in Figure 19 for the alloys aged at 400°C for 200 h. These plots indicate clearly that the supersaturated solid solution decomposes into two phases: a Cu-Ni-Fe-rich phase with a poor content of Fe and a Ni-Cu-Fe-rich phase with a lower content of Cu and a higher content of Fe. The decomposed Cu-Ni-Fe and Ni-Cu-Fe phases are in agreement with the miscibility gap of the equilibrium Cu-Ni-Fe phase diagram [17].
The simulated microstructural evolution for the Cu-46at.%Ni-4at.%Fe alloy aged at 400°C for different times is shown in Figure 20. The white and black zones correspond to the Ni-Cu-Fe-rich and Cu-Ni-Fe-rich phases, respectively. The morphology is irregular and interconnected, and it has no preferential alignment in a specific crystallographic direction (Figure 20(a–c)). This type of microstructure has been named isotropic [4]. The volume fraction of phases is similar because of the small difference in chemical composition. For higher temperatures, the isotropic morphology at the early stages of aging, as the aging progresses, changes to a cuboid or plate shape, crystallographically aligned along the <100> directions because of the low coherency-strain energy associated with these directions [9].
Figure 21 shows the FIM microstructural evolution of the Cu-46at.%Ni-4at.%Fe alloy aged at 400°C for 50 h. FIM image of the solution treated and quenched sample shows the characteristic concentric ring of the (001) plane in both alloys. In the case of the aged samples, brightly imaged zones correspond to the Ni-Cu-Fe-rich phase and dark zones to the Cu-Ni-Fe-rich phase (matrix). The morphology of the decomposed phases was also cuboids or plates aligned in the <100> directions as the aging progressed.
6. Summary
The
Acknowledgments
The authors wish to acknowledge financial support from Tecnológico Nacional de México/Instituto Tecnológico de Pachuca, and SIP-IPN and Conacyt.
References
- 1.
Cahn J.W., Hilliard J.E., Free energy of a nonuniform system I. Interfacial free energy. J. Chem. Phys. 1958; 28, 258–267. - 2.
Allen S.M., Cahn J.W., A microscopic theory for antiphase boundary motion and its application to antiphase domain coarsening. Acta Metal. 1979; 27, 1085–1095. - 3.
Ginzburg V.L., Landau L.D. On the theory of superconductivity. Zh. Eksp. Teor. Fiz. 1950; 20, 1064–1082. - 4.
Hilliard J.E., Phase transformations, ASM, 1970. - 5.
Raabe D. Computational materials science, Wiley-VCH, 1998. - 6.
Rappaz M., Bellet M., Deville M. Numerical modelling in materials science and engineering, Springer, 1998. - 7.
Avila-Davila E.O., Lopez-Hirata V.M., Saucedo-Muñoz M.L., Gonzalez-Velazquez J.L. Microstructural simulation of phase decomposition in Cu-Ni alloys. J. Alloys Comp. 2008; 460, 206–212. - 8.
Lifshitz I.M., Slyozov V.V. The kinetics of precipitation from supersaturated solid solution. J. Phys. Chem. Solids. 1961; 19, 35–50. - 9.
Kostorz G. Phase transformations in materials, Wiley-VCH, 2001. - 10.
Avila-Davila E.O., Lezama-Alvarez S., Saucedo-Muñoz M.L., Lopez-Hirata V.M. Gonzalez-Velazquez J.L., Simulación numérica de la descomposición espinodal en sistemas de aleación hipotéticos A-B y A-B-C. Revista de Metalurgia. 2012; 48, 223–236. - 11.
Baker H. Alloy phase diagrams handbook, Vol. 3, ASM International, 1992. - 12.
Moser Z., Zakulski W., Spencer P., Hack K. Thermodynamic investigations of solid Cu-Ni and Fe-Ni alloys and calculation of the solid state miscibility gap in the Cu-Fe-Ni system. CALPHAD. 1985; 9, 257–269. - 13.
Pearson W.B. A handbook of lattice spacing and structures of metals and alloys, ASM International, 1986. - 14.
Mehrer H. Diffusion in solid metals and alloys, Springer-Verlag, 1990. - 15.
Dieter G.E. Mechanical metallurgy, Mc Graw Hill, 2001. - 16.
Soriano-Vargas O., Avila-Davila E.O., Lopez-Hirata V.M., Dorantes-Rosales H.J., Gonzalez-Velazquez J.L. Spinodal decomposition in an Fe-32 at% Cr alloy during isothermal aging. Mater. Trans. JIM. 2009; 50, 1753–1757. - 17.
Lopez-Hirata V.M., Sano N., Sakurai T., Hirano K. A study of phase decomposition in Cu-Ni-Fe alloys. Acta Metall. Mater. 1993; 41, 265–271. - 18.
Avila-Davila E.O., Melo-Maximo D.V., Lopez-Hirata V.M., Soriano-Vargas O., Saucedo-Muñoz M.L., Gonzalez-Velazquez J.L. Microstructural simulation in spinodally-decomposed Cu-70 at.%Ni and Cu-46 at.%Ni-4 at.%Fe alloys. Mater. Characterization. 2009; 60, 560–567. - 19.
Honjo M., Saito Y. Numerical simulation of phase separation in Fe-Cr binary and Fe-Cr-Mo ternary alloys with use of the Cahn-Hilliard equation. ISIJ Int. 2000; 40, 914–919.