| Home | E-Submission/Review | Sitemap | Editorial Office |  
Korean Journal of Metals and Materials > Volume 56(6); 2018 > Article
Chang, Kwon, and Lee: Phase-field Modeling of Precipitate Behavior in RPV Steel Using CALPHAD Database


We developed a multi-scale modeling framework using the thermodynamic UW1 database of the Fe-Mn-Ni-Si quaternary system. To obtain input data for the phase-field model, such as mobility and interfacial energy data, we used DICTRA software and the TC-PRISMA package with the commercial TCFE8 and MOBFE3 databases. Using the developed framework, we evaluated the stability of single precipitates and compared the results to those from the thermodynamic model. We also investigated the precipitate interactions of paired particles of identical and different phases.


Predicting multi-particle behavior is an extremely complicated task. The Ostwald ripening of single-phase precipitates has been extensively investigated using both analytical [1-8] and numerical methods [8-15]. The role of the coarsened phase volume fraction has been explored [3,5,6,9,10] and morphological deviations from the spherical (circular) coarsened phases at high volume fractions have already been observed in previous studies [6,9,10,15]. Ardell evaluated the role of the coarsened phase fraction using a theoretical approach and he found that the peak of the particle size distribution decreased as the particle fraction increased [3]. Brasilsford and Wynblatt [7] also found that the particle size distribution peak position decreased as the volume fraction increased, and they claimed that the effect was much less sensitive to volume fraction than proposed by Ardell [3]. Although the microstructural evolutions driven by curvature minimization have been studied extensively for many decades, current understanding of Ostwald ripening remains limited. Because we have no analytical solutions for the diffusion equations in most cases and ripening is extremely complex, multi-particle behaviors are usually predicted based on the mean-field assumption [1-8]. or numerical approaches [6,8-10,15]. Previous simulations of particle coarsening have predicted the behaviors of single phases using phenomenological models [9,10,15]. In this study, we performed a simulation of multiphase precipitates using a CALPHAD-type database to predict the behaviors of realistic precipitates. The major elements of low-Cu reactor pressure vessel (RPV) steel are Fe, Mn, Ni and Si [16-18]; therefore, we selected the Fe–Mn–Ni–Si (hereafter FMNS) quaternary system to investigate multi-phase coarsening, using the FMNS database. The phase-field method has been extensively applied to investigate the microstructural evolution of low-alloy Fe-based steel [19-21] assuming only one precipitate phase within the system. T3, T6, and T7 Mn–Ni–Si-rich (MNS) precipitates are generally observed in RPV steel [16]; therefore, we assumed the presence of three types of MNS precipitates in the system. We adopted the Kim–Kim–Suzuki (KKS) model [22] to perform multi-phase and multi-component phase-field modeling. The FMNS system, a prototype of RPV steel, is of great interest in nuclear engineering applications [16-18,23,24]. Because MNS precipitates are among the main sources of late-stage hardening for FMNS [16-18,23,24], we examined the stability of the MNS precipitates. The T3, T6, and T7 precipitate phases, major types of MNS precipitates, were considered in the UW1 database [16] in this study; we also simulated the evolution of pairs of T3–T3, T7–T7, and T3–T7 precipitates to quantify the interactions between particles.


In the UW1 database [16], one bcc phase representing the matrix and 12 MNS precipitate phases are considered. In our study, we selected the matrix phase and three MNS precipitate phases for simplicity. The thermodynamic parameters we used were taken from the supplementary material of ref. [16] as follows:
for bcc FMNS phase,
LFe,Mnbcc0=-2759 + 1.23T
LFe,Nibcc0=956.63  1.28726Τ
LFe,Nibcc1=1789.03  1.92912Τ
LFe,Sibcc0=153138.56 + 46.86T
LFe,Sibcc1=92352, LFe,Sibcc2=62240
LMn,Nibcc0=3508.43  23.7885T
LMn,Sibcc0=89620.7 + 2.9410T
TFe,Mnbcc0=123.0, TFe,Sibcc0=504.0
for the T3: Mn6⁄29Ni16⁄29Si7⁄29 phase,
for the T6: Mn1/3(Ni,Si)2/3 phase,
LMn,Ni,Sibcc0=159474.81, LMn,Ni,Sibcc1=172110.47
and for the T7: Mn1⁄2Ni1⁄3Si1⁄6 phase,
gMn,Ni,SiT70=32434.25  5T+1/20gMncbcc+1/30gNifcc+1/60gSidiamond


We utilize the phase-field model to simulate the microstructural evolution of the FMNS system by solving the Cahn–Hilliard [25] and Allen-Cahn (Ginzburg-Landau) equations [26] to simulate the microstructural evolution. We will denote the composition (i = 1, 2, 3, and 4 for Fe, Mn, Ni, and Si, respectively) in phase θ using ciθ(r,t) at position r and time t. θ indicates the T3, T6, and T7 phases. We introduce four non-conserved order parameters (ϕi) to indicate the regions of the four precipitated phases. The composition ci(r,t) is given as follows [22]:
where [27],
where np represents the number of θ precipitate.
The local free energy density G(ciθ,t) of the system is expressed as follows:
where the units of all free energy densities (expressed by g in this study) in this study are joules per mole: the unit of total free energy of the system G is J. We assumed that density ρ is 7.9 g cm3 and the atomic molar mass M is 56.
with c1α=1.0-c2α-c3α-c4α. We chose ωp = 200.0 and γ = 5.0 in Eq. 5. From ref. [28], we obtain the free energy of each element of the α phase. The free energies of the T3, T6, and T6 precipitates is given as follows:
where y3II and y4II denote the site fraction of Ni and Si at the second sub-lattice of the T6 phase, respectively.
We solve the Ginzburg-Landau equation in Eq. 9 and the Cahn-Hilliard equation in Eq. 10.
ϕiθr,tt=-Lθδgδϕiθr,t θ= T3, T6, and T7
We implemented the forward Euler scheme to solve Eq. 9 and 10. To determine Di in Eq. 10 we adopted the relation in refs. [20,21].
We use the relation DiθT=Diθ0exp-QiθRT to determine the diffusivity; the parameters used to determine this value are listed in Table 1.
In principle, one has to know the diffusivity data of each element of each precipitated phase to perform phase-field modeling. However, the diffusivity data of T3, T6, and T7 phases are not available. Since a precipitate larger than 17 nm has the FCC crystal structure [30], we adopted the diffusion data of each element in the Austenite (FCC) phase for all precipitated phases.
We could not find diffusivity data for Si in the γ(fcc) phase; therefore, we calculate it from a kinetic simulation. We performed a DICTRA simulation using Thermo-Calc 2017a software with the TCFE8 and MOBFE3 databases. At T = 550K, in the Fe–Si binary system (γ(fcc) phase) the Si mole fraction is always 0.03 at x = 0; the value is 0.005 initially, except for the point of x = 0 in Fig. 1. From the diffusion length depending on the simulat ion time, we calculated the diffusion coefficient of Si in the γ(fcc) phase.
We restated Eq. 9 as Eq. 12.
ϕiθr,tt=-LθGϕiθr,t-ωθ2ϕiθ(θ = T3, T6, and T7)
In our simulations, we use the energy normalized by RT, where R = 8.3144598 J/mol·K and the temperature T = 550 K. The diffusivity values are normalized by D2α(T). For the length of the grid spacing, Δx = Δy was taken as 0.125 nm [21]. The time in our study is normalized by 0.01 s.
To determine ωθ in Eq. 12, the interfacial energies between the matrix and precipitate phases are needed. So far, we have not found any literature that has evaluated the interfacial energies relevant for our study; therefore, we measure the interfacial energy using the extended Becker’s model function [31]. To perform the interfacial energy estimation, we used TC-PRISMA software (included in Thermo-Calc 2016a) with the TCAL3 and MOBAL3 databases. However, evaluating the interfacial energy for the whole Si fraction in the T6 phase using TC-PRISMA is not possible. We find that the Si fraction in the T6 phase is generally smaller than 0.3; therefore, we assume that the interfacial energy is given by Eq. 13. 0.2496 J/m2 [22]. For the T3 and T7 stoichiometric compounds, we also evaluate the interfacial energies using TC-PRISMA as 0.3311 and 0.1791 J/m2. The interfacial energy is given by Eq. 13 when the gradient coefficient is ε2.
Since we used the normalized variables in this study, we evaluated the normalized interfacial energy.
The unit of ε is J⁄m2⋅mol. We chose ε*,T6 = 1.0 as a non-dimensional parameter which corresponds to the interfacial energy 0.26 J/m2 which is quite comparable to the value 0.2496 J/m2 which we obtained from the TC-Prisma calculation. We also determined ε*,T3 = 1.33 and ε*,T7 = 0.72 using a consistent way.


We performed the simulations with a simulation cell size of 8 × 8 nm2. The initial precipitate radius was 1.5 nm. The discretized time step was selected as given in Table 2.
We chose Δt after convergence and accuracy testing. The simulations obtained results in 1 to 12 h, depending on the Δt of the system. For the matrix, the initial compositions in the α phase were c2α = 0.008, c3α = 0.008, and cα = 0.00808. For the T6 precipitate, c3T6=0.459234067 and c4T6=0.2074326.
For the T3 and T7 cases, the compositions of the precipitates were fixed and the initial compositions of the α phase were assumed to be equal to those of the T6 case. For the case with a single precipitate, we placed a circular precipitate with an initial radius r = 1.5 nm at the center of the simulation cell. For multiple precipitate phases, we placed two precipitates with initial radii r = 1.5 nm at two symmetrical positions with the distances L = 3.25, 3.50, and 3.75 nm relative to the center of the system.


We performed the simulation to examine the stability of a single precipitate as shown in Fig. 2. In Fig. 3(a), the T6 precipitate undergoes minor accommodation at a very early stage. After 0.0002 s, the precipitate radii remains the same as it was during the modeling in Fig. 3(b). This means that the T6 precipitate is thermodynamically stable. Because the T3 and T7 phases are unstable, these precipitates shrink with time; the T7 precipitate radius decreases more rapidly than that of the T3 precipitate. We performed CALPHAD modeling using Thermocalc software and the UW1 database [16] with c2α = 0.008, c3α = 0.008, c4α = 0.00808, and T = 550 K. The equilibrium phases were bcc(A2) (98.940 mol%) and T6 (1.363 mol%). Therefore, we concluded that the phase-field simulation with a single precipitate was consistent with the thermodynamic model. As shown in Fig. 4, the diffusion field of the solutes overlapped between the particles, therefore, the solute concentration is relatively higher than in other matrix regions. As a result, the dissolution rate decreases along the direction of the other particle. As a result, the particle shape becomes asymmetric in Fig. 4.
Therefore, the T3 particles remain longer for two interacting T3 particles compared to a single T3 particle. In addition, when the distance between particles L is 3.75 nm, the particle size evolution curve shown in Fig. 5 approaches convergence with that curve of the single T3 particle. Consistent results are observed for the T7 precipitate in Fig. 6.
We performed a set of simulations to investigate the interactions between T3 and T7 precipitates, with T3 and T7 precipitates located as described in Fig. 7 and 8.
In Fig. 9 and 10, we find that the T3 and T7 precipitates survive longer under T3–T7 interactions. Even though the precipitates are of different phases, particle–particle interactions increase the lifetime of the precipitates, i.e., enhance their stability.


We developed a multi-scale modeling framework to simulate the microstructural evolution of MNS precipitates in the quaternary FMNS system representing RPV steel. The UW1 quaternary database was implemented and the DICTRA and TC-PRISMA packages were used to evaluate the mobility of solutes and the interfacial energies between the matrix phase and precipitate phases. The results of the single-precipitate simulations were consistent with the thermodynamic model. We also evaluated the interactions between precipitates and found that the precipitates survived longer when two particles were present, whether of the same or different phases, compared to the particle lifetime in the single-precipitate case.


This work was supported by a National Research Foundation of Korea(NRF) grant funded by the government of Korea (Ministry of Science and ICT) (NRF-2017M2A8A4015157).

Fig. 1.
Profile of Si mole-fraction at four different times.
Fig. 2.
ϕT3 profile when t = 6 s. The figure is visualized by mapping a summation of the ϕT3 value to a gray scale (maximum 1.0, minimum 0.0).
Fig. 3.
(a) Evolution of radii of T6 precipitate at very early stage (<0.0002 s), (b) precipitate radii of T3, T6, and T7 phases with respect to time for the case of single-precipitate simulation.
Fig. 4.
ϕ1T3 and ϕ2T3 profiles when t = 6 s. The distance between the particle centers L = 3.25 nm when t = 6 s. The figure was visualized by mapping a summation of the ϕ1T3 + ϕ2T3 value to a gray scale (maximum 1.0, minimum 0.0).
Fig. 5.
Precipitate radius of T3 with respect to time for a single T3 precipitate versus two T3 precipitates separated by different distances L.
Fig. 6.
Precipitate radius of T7 precipitate with respect to time for a single T7 precipitate and two T7 precipitates with varying distances between particles L.
Fig. 7.
Profile of Mn concentration when t = 1.5 s. The distance between the particle centers L is 3.25 nm. The figure was visualized by mapping the c2 value to a gray scale (maximum 0.5, minimum 0.05).
Fig. 8.
Profile of Ni concentration when t = 1.5 s. The distance between the particle centers L is 3.25 nm. The figure was visualized by mapping the c3 value to a gray scale (maximum 0.6, minimum 0.0).
Fig. 9.
Precipitate radius of T3 for a single T3 precipitate and paired T3–T7 precipitates with varying distances between particles L.
Fig. 10.
Precipitate radius of T7 for a single T7 precipitate and paired T3–T7 precipitates with varying distances between particles L.
Table 1.
The values of the maximal diffusion coefficient and activation energy of each element [29].
Elements Phase 0Diθ (m2/s) Qiθ (J/mol)
Mn α(bcc) 1.5 × 10–4 2.34 × 105
Mn γ(fcc) 1.6 × 10–5 2.62 × 105
Ni α(bcc) 4.2 × 10–3 2.68 × 105
Ni γ(fcc) 7.7 × 10–5 2.81 × 105
Si α(bcc) 1.7 × 10–4 2.29 × 105
Table 2.
Normalized Δt for six precipitate configurations.
Precipitate Δt
Single T3 1.0 × 10–4
Single T6 1.0 × 10–4
Single T7 1.0 × 10–4
T3-T3 1.0 × 10–4
T7-T7 1.0 × 10–4
T3-T7 1.0 × 10–5


1. I. M. Lifshitz and V. V. lyozov, J. Phys. Chem. Solids. 19, 35 (1961).
2. C G Garay-Reyes, S. E. Hernández-Martínez, J. L. Hernández-Rivera, J. L. Hernández-Rivera, E. J. Gutiérrez-Castañeda, H. J. Dorantes-Rosales, J. Aguilar-Santillan, and R. Martínez-Sánchez, Met. Mater. Int. 23, 298 (2017).
crossref pdf
3. A. Ardell, Acta Metall. 20, 61 (1972).
4. P. W. Voorhees, J. Stat. Phys. 38, 231 (1985).
5. C. Davies, P. Nash, and R. N. Stevens, Acta Metall. 28, 179 (1980).
6. P. Voorhees and M. Glicksman, Metall. Mater. Trans. A. 15, 1081 (1984).
crossref pdf
7. A. Brailsford and P. Wynblatt, Acta Metall. 27, 489 (1979).
8. J. H. Yao, K. Elder, H. Guo, and M. Grant, Phys. Rev. B. 47, 14110 (1993).
9. D. Fan, S. Chen, L. Q. Chen, and P. W. Voorhees, Acta Mater. 50, 1895 (2002).
10. S. G. Kim, Acta Mater. 55, 6513 (2007).
11. J. A. Warren and B. T. Murray, Modelling Simul. Mater. Sci. Eng. 4, 215 (1996).
12. T. Kraska, J. Phys. Chem. B. 112, 12408 (2008).
13. G. Pévot, Phys. Rev. B. 84, 045434 (2011).
14. P. W. Voorhees and M. Glicksman, Acta Metall. 32, 2013 (1984).
15. P. W. Voorhees, G. McFadden, R. Boisvert, and D. Meiron, Acta Metall. 36, 207 (1988).
16. W. Xiong, H. Ke, R. Krishnamurthy, P. Wells, L. Barnard, G. R. Odette, and D. Morgan, MRS Commun. 4, 101 (2014).

17. D. Sprouster, J. Sinsheimer, E. Dooryhee, S. Ghose, P. Wells, T. Stan, N. Almirall, G. Odette, and L. Ecker, Scripta Mater. 113, 18 (2016).
18. S. Hong, C. Lee, M.-C. Kim, and B.-S. Lee, Korean J. Met. Mater. 55, 752 (2017).
19. T. Koyama and H. Onodera, Mater. Trans. 46, 1187 (2005).
20. T. Koyama, K. Hashimoto, and H. Onodera, Mater. Trans. 47, 2765 (2006).
21. S. Biner, W. Rao, and Y. Zhang, J. Nucl. Mater. 468, 9 (2016).
22. S. G. Kim, W. T. Kim, and T. Suzuki, Phys. Rev. E. 60, 7186 (1999).
23. G. R. Odette and G. Lucas, JOM. 53, 18 (2001).
24. G. R. Odette and G. Lucas, Radiat Eff. Defects Solids. 144, 189 (1998).
25. J. W. Cahn and J. E. Hilliard, J. Chem. Phys. 28, 258 (1958).
26. S. M. Allen and J. W. Cahn, Acta Metall. 27, 1085 (1979).
27. J. Zhu, T. Wang, A. Ardell, S. Zhou, Z. Liu, and L. Chen, Acta Mater. 52, 2837 (2004).
28. A. Dinsdale, Calphad. 15, 317 (1991).
29. Metals data book, Japan Institute of Metals, Maruzen, Tokyo, Japan (2004).

30. P. J. Othen, M. L. Jenkins, and G. D. W. Smith, Phil Mag A. 70, 1 (1994).
31. R. Becker, Anna. Phys. 424, 128 (1938).

PDF Links  PDF Links
PubReader  PubReader
ePub Link  ePub Link
Full text via DOI  Full text via DOI
Download Citation  Download Citation
Related article
Editorial Office
The Korean Institute of Metals and Materials
6th Fl., Seocho-daero 56-gil 38, Seocho-gu, Seoul 06633, Korea
TEL: +82-70-4266-1646   FAX: +82-2-557-1080   E-mail: metal@kim.or.kr
About |  Browse Articles |  Current Issue |  For Authors and Reviewers
Copyright © The Korean Institute of Metals and Materials. All rights reserved.                 Developed in M2Community