Mechanics & Industry
Volume 21, Number 5, 2020
Scientific challenges and industrial applications in mechanical engineering
|Number of page(s)||12|
|Published online||11 August 2020|
Fast simulation of grain growth based on Orientated Tessellation Updating Method
Navier, CNRS, École des Ponts ParisTech,
6 & 8 Ave Blaise Pascal,
Marne La Vallee,
2 LMS, École Polytechnique, CNRS, Université Paris-Saclay, 91128 Palaiseau, France
* e-mail: email@example.com
Accepted: 28 April 2020
This work is part of a more general idea consisting in developing a macroscopic model of grain growth whose state variables contain for each material point the statistical descriptors of the microstructure (e.g., disorientation, grain size and shape distributions). The strategy is to determine macroscopic free energy and dissipation potentials on the basis of a large number of computations at the scale of the polycrystal. The aim is to determine enriched macroscopic evolution laws. For sake of simplicity, this contribution only deals with grain growth of a single phased metal without diffusion or segregation of alloying elements. In order to test this upscaling strategy it is necessary to establish a simulation tool at the scale of the polycrystal. It should be sufficiently simple and fast to enable a large number of simulations of various microstructures, even if it leads to neglect some phenomena occurring at this scale. Usual grain growth models relying on mobile finite element modeling, level set functions, phase field or molecular dynamics are too computationally costly to be used within the proposed framework. Therefore, this paper focuses on the development of a “toy” model. Tessellation techniques are usually used to approximate polycrystalline microstructures. Therefore, one can approximate the real evolution of the microstructure as a succession of tessellation approximations. It then becomes quite natural to attempt to establish the evolution law of the microstructure directly on the parameters defining the tessellation. The obtained model is very light in terms of computational cost and enables to compute a large number of evolutions within the framework of the proposed statistical upscaling method.
Key words: Grain growth / Voronoi-Laguerre tessellation / grain boundary energy / dissipation / grain mobility
© D. Weisz-Patrault et al., published by EDP Sciences 2020
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This work is part of a more general idea consisting in developing a macroscopic model of grain growth whose state variables contain for each material point the statistical descriptors of the microstructure (e.g., disorientation, grain size and shape distributions). Statistical distributions are too much detailed to be processed at each material point, thus reduced information is considered instead: means and variances at each material point. This ambition arises as very few information about microstructure is usually processed at the macroscopic scale. The general framework is the standard generalized media characterized by a free energy per unit mass and a dissipated power potential. These two potentials arise in the balance energy equation combining the first and second laws of thermodynamics. The balance energy equation should be verified for all possible virtual evolution. The two potentials depend on macroscopic state variables and their time derivatives. The determination of the macroscopic free energy and dissipation potentials enables to establish the evolution laws of the state variables characterizing the microstructure evolution at macroscale.
The strategy is to determine macroscopic free energy and dissipation potentials not axiomatically (with parametric functions) and calibration with experiments at macroscale, but on the basis of a large number of computations at the scale of the polycrystal. The aim is to quantify for each microstructure evolution the total grain boundary energy and the total dissipated power as a function of the macroscopic state variables that characterize the statistical distribution of the microstructure in order to determine enriched macroscopic evolution laws. For sake of simplicity, this contribution only deals with grain growth of a single phased metal without diffusion or segregation of alloying elements. In order to test this upscaling strategy it is necessary to establish a simulation tool at the scale of the polycrystal. It should be sufficiently simple and fast to enable a large number of simulations of various microstructures, even if it leads to neglect some phenomena occurring at this scale. Usual grain growth models relying on cellular automaton [1–4], mobile finite element modeling , level set functions [6,7], phase field [8–11] or molecular dynamics [12–14] have been intensively studied and very interesting results have been obtained. However, the computational cost of such approaches is usually incompatible with an intensive use as suggested within the proposed framework. Vertex methods [15–18] consist in establishing the evolution law directly at the triple junctions and are sufficiently simple in two dimensions to reach short computation time. However, the extension in three dimensions is very difficult.
Evolution of polycrystalline microstructures (such as grain growth) involves coupled mechanisms at different scales: (i) the atomic scale (crystal lattice), (ii) the microscale including mechanisms at grain boundaries (iii) the mesoscale revealing the grain structure (i.e., polycrystalline Representative Volume Element) and (iv) the macroscopic scale modeled as a continuous medium. Thus, this paper focuses on the development of a “toy” model formulated at mesoscale and enabling to simulate grain growth within short computation time.
Voronoi-Laguerre tessellation techniques are usually used to approximate polycrystalline microstructures at mesoscale. Very efficient algorithms have been developed with the possibility of controlling statistical distributions of grain size and shape. Crystal lattice orientation can also be specified for each grain. The tessellation equipped with such an orientation field is called an Orientated Tessellation (OT). One can approximate the real evolution of the microstructure as a succession of OT approximations. It then becomes quite natural to attempt to establish the evolution law of the microstructure directly on the parameters defining the OT. This modeling strategy is called in this contribution the Orientated Tessellation Updating Method (OTUM).
As grain growth during annealing is essentially viscous, the time derivatives of the OT parameters are obtained as a function of the thermodynamic forces defined as the partial derivatives of the total grain boundary energy with respect to the OT parameters. This mesoscale evolution law is obtained through the energy balance equation by specifying two mechanisms at microscale: (i) the grain boundary energy (disorientation) and (ii) the dissipated power by crystal visco-plasticity through any grain boundary virtual motion. For sake of simplicity, crystal lattice rotation through crystal visco-plasticity is not considered in this contribution. In addition, plane polycrystals are considered with three slip directions in each grain. Thus, the crystal lattice is plane hexagonal. This configuration corresponds to the plane <1, 1, 1> of a face-centered cubic (fcc) crystal, as shown in Figure 1. Disorientation between two neighboring grains (characterized by five parameters in 3D) is characterized only by two parameters in 2D: the disorientation angle Δθ and the orientation of the grain boundary φ. Thus, the grain boundary energy considered in this paper is computed from fcc crystals sharing the same orientation <1, 1, 1> (asymmetrictilt boundaries).
The dissipative mechanism at microscale is calibrated to correspond to the classic notion of grain boundary mobility. The overall mesoscale grain growth model is validated by comparison with the classical curvature driven shrinkage of a spherical grain. The proposed model is very light in terms of computational cost and enables to compute a large number of microstructure evolutions, which is useful for the general statistical upscaling method.
FCC and slip directions in <1, 1, 1> plane.
A tessellation is defined by n seeds whose Cartesian coordinates are denoted by (xj, yj) and n weights denoted by wj (where 1 ≤ j ≤ n). Coordinates are dimensionless as well as weights. The tessellation is completely determined by the parameter vectors , , . The dimensionless tessellation is scaled by a physical length L0 that represents the length of the tessellation edge.
It is clear from the definition (1) that weights are defined up to a constant. Thus, the following constraint is added to obtain a univocal definition: (2)
Thus, weights w lie in an affine hyperplane of dimension n − 1 and denoted by , whose support is the hyperplane denoted by . In addition, it should be noted that a cell Cj may be empty as shown in Figure 2. This property will be intensively used as some grains should disappear during grain growth.
The crystallographic orientation is denoted by . Since the crystal lattice is plane hexagonal . The parameter set defining the OT (see Fig. 3) reads: (3)
A total energy is associated to the OT, where T is the temperature. In this contribution, the total grain boundary energy per unit depth is considered. The Orientated Tessellation Updating Method consists in computing a time evolution of the OT characterized by by deriving an evolution law of the form: (4)
where M(T, α) is a second order tensor to be defined by introducing a dissipative mechanism at microscale and is the driving force.
Seed and weight in tessellation.
Seed and weight in tessellation.
As mentioned in the introduction, the mesoscopic evolution law (4) is determined by introducing two local mechanisms at microscale: (i) the grain boundary energy and (ii) the dissipated power by crystal visco-plasticity through any grain boundary virtual motion. Grain boundaries are indexed by pairs (i, j) where i and j denote two neighboring grains. The set of pairs of neighboring grains defining grain boundaries is denoted by IGB : (6)
The condition j ≥ i is meant to count each grain boundary only once (otherwise if (i, j) ∈ IGB then (j, i) would also be an element of IGB).
In addition φij is the angle of the grain boundary as shown in Figure 4.
Several analytical formulations have been proposed for E. The simplest and most widely used function is the Read & Shockley formula : (9)
However, the range of validity of the Read & Shockley formula is limited to small disorientation angles and do not account for the energy cusps at certain disorientation angles. To overcome this difficulty, Wolf  introduced a piecewise interpolation function for the grain boundary energy per unit area of the form: (10)
This interpolation function is usually fitted on atomistic simulations. Other interpolations have been proposed within a three dimensional framework . Since the main goal of this contribution is to demonstrate the capability of the OTUM and not to present a study on a specific metal, a simplified grain boundary function is directly defined in the form of (10). Since the crystal lattice is plane hexagonal, a cusp is introduced at π∕3. In addition it is assumed that the dependance on φ is negligible. The grain boundary energy per unit area is defined as follows: (11)
where G(T) is the shear coefficient. In this paper, data for pure iron  are used to calibrate G(T): (12)
where bG ≈ 88134 MPa and aG ≈−24 MPa.K−1. The temperature dependence G(T)∕G(0) is a simple choice similar to . In addition: (13)
where a1 = a2 = 1 and J.m−2, which roughly corresponds to pure iron. The resulting simplified grain boundary energy per unit area is presented in Figure 5.
In addition to the grain boundary energy per unit area, the dissipated power per unit area through any virtual grain boundary motion should be specified. Consider a virtual normal velocity of the grain boundary (i, j) ∈ IGB. The dissipated power per unit area reads: (14)
Since grain growth is viscous, the analytical computation proposed by [24,25] is simply adapted for crystal visco-plasticity by considering that the critical shear stress τc linearly depends on : (17)
In other words, the driving force (proportional to γ(Δθ)) tends to zero when the disorientation tends to zero. However, the dissipative cost (proportional to X(Δθ)) that resists to the driving force, tends to zero much faster because of the logarithmic term in (19). Thus, for arbitrarily small disorientation, the grain boundary is unstable as the driving force is arbitrarily large in comparison to the resistive mechanism (the ratio is in ).
Grain boundary energy per unit area at T = 0 K.
The mesoscale evolution law is derived by considering the total grain boundary energy per unit depth and the total dissipated power per unit depth through any virtual grain boundary normal velocity field, which are defined as follows: (20)
where lij the joint length and . The dependance on α = (x, y, w, θ) comes from the grain boundary lengths lij that depend on (x, y, w) and that depend on θ. The number of grain boundary is . There exists a bijection between IGB and that definesthe numbering of grain boundaries. Thus, for each grain boundary denoted by (i, j) ∈ IGB there exists a unique . Therefore the total dissipated power per unit area reads: (21)
It is clear in (23) that for each grain boundary an arbitrary choice is made for the positive direction of the normal velocity , which has no consequence as the square of the virtual velocity arises in the dissipated power. Hence: (25)
The constraint (2) leads to: (27)
where is a dimensionless second order tensor of size n × n.
The “real” weight variation is determined for isothermal (i.e., Ṫ = 0) and homogenous temperature fields (i.e., ∇T = 0) by using the general energy balance equation (combining the first and second thermodynamic laws), which is established for “real” evolution (and not virtual): (29)
Thus, the energy balance (30) arises as a constraint on the “real” weight variation. The maximum dissipation principle , enables to determine the “real” weight variation as the maximum of the dissipation under the constraint (30): (31)
The evolution law (34) should be inverted. However R(α) is not invertible because of the constrain (27). The kernel of R(α) is . Thus, there are n − 1 strictly positive singular values of R(α) denoted by (λ1, ⋯, λn−1) and there exist two orthogonal tensors U and V such as: (35)
The singular value decomposition is computed in Scilab  and gives (λ1, ⋯, λn−1), U and V. The Moore-Penrose pseudo-inverse matrix is introduced as follows: (36)
It should be noted that a scale effect is demonstrated in the evolution law (37) as the physical dimension L0 arises.
The driving force involved in the evolution law (37) should be evaluated. An analytic solution is derived in this section, which contributes to reach short computation time. From (20) and (11) one obtains: (39)
Thus, the third order tensor β representing angles at the triple junctions as shown in Figure 6 are defined as follows: (43)
where the following symmetry rule holds βijk = βkji.
Consider a virtual weight variation (i.e., verifying the constrain (27)). One obtains: (44)
where δ is the Kronecker symbol.
where RI and RII are the radii of curvature along the principal directions (in three dimensions). For plane polycrystals RI = R and RII = ∞. A simple linear evolution law accounting for viscosity is usually used: (48)
where is the inward normal speed of the considered grain boundary and the classical grain boundary mobility for the curvature driven model (where the symbol refers to the curvature driven model). In the following, a circular grain shrinkage is compared to an hexagonal configuration as shown in Figure 7 in order to determine the mobility m(T) introduced in the dissipated power (17) as a function of . Thus for the curvature driven model the inward normal speed reads: (49)
The proposed approach is applied to the hexagonal situation. The partial derivative (46) reads: (50)
The driving force (40) reads: (52)
The second order tensor K introduced in (26) reads for this configuration: (53)
The second order tensor χ introduced in (22) reads: (54)
where I is the identity of size 6 × 6. The second order tensor R introduced in (28) reads: (55)
Thus, the second order tensor M introduced in (38) reads: (56)
Hence the evolution law (37) for the hexagonal situation: (57)
Thus, by using (23) one obtains the inward normal speed according to the proposed model: (58)
For the sake of simplicity it is assumed that depends on disorientation in 1∕X(Δθ) even though more realistic mobilities could be obtained by considering m(T, Δθ) instead of m(T) in (17). (60)
To estimate m(T) one can use the pure iron grain boundary mobility determined by : (61)
where the molar volume of fcc-Fe is Vm = 7.09 × 10−6 m3, the boundary thickness is δ = 1 nm, the burgers vector is b = 2.48 × 10−10 m, the gas constant is R = 8.3144621 J.mol−1.K−1. In addition, DFe(T) is the diffusivity of Fe atoms along the grain boundary (for fcc) and follows an Arrhenius law: (62)
where the pre-factor is m2.s−1 the activation energy is Q = 148 kJ.mol−1.
Validation by comparison to curvature driven spherical grain shrinkage.
The temperature kinetics is imposed in the following, thus T(t) is a known function. Voronoi-Laguerre tessellations (x, y, w) are generated with the free software Neper  and orientations θ are set in Scilab . Weights are initially set to a common value and are updated at each time step according to the evolution law (37) and the driving force (46). The numerical discretization consists in a simple explicit scheme with adaptive time step in order to capture accurately growth rate variations. The time step is determined as follows: (63) (64)
where δw > 0 is a constant value. As the simple model proposed in this paper is purely viscous, even though the temperature is maintained constant, grain growth never stops strictly speaking. The evolution becomes only slower and slower. Thus, a stopping criterion needs to be defined. The algorithm is simply presented in Figure 8.
Different Orientated Tessellations with different distributions of sphericity and grain size are generated. The orientation field correspond to a disorientation distribution very similar for each OT as shown in Figure 9. The temperature is set to T = 800 °C and grain growth is simulated by Orientated Tessellation Updating Method during approximately 4 h. The code is not optimized as Scilab has been used to compute the evolution law and Neper to generate tessellations. (Computation time is roughly 15 min for 1000 time steps and is mainly due to writing and reading text files between Scilab and Neper. Thus, a unifed code not relying on text files would significantly reduce computation time.)
Grain growth at 800 °C is presented for three OT (denoted by OT1, OT2, OT3) in Figures 10, 11 and 12 (orientations are in degrees). Simulations show grain coarsening and are consistent with experimental observations  on an other material (tantalum). Grain growth is heterogeneous depending on the local disorientation and grain size arrangement leading to heterogeneous grain size distribution during the transient state. Of course, grain size heterogeneity eventually reduce as bigger grain tends to grow at the expense of smaller grains. In addition, the disorientation distribution tends to favor small disorientations as shown in Figure 9. Sphericity tends to increase as shown in Figure 11 but this may be due to the fact that weights w are the only state variables. Future developments will include seeds x and y as state variables,which will enable to consider additional degrees of freedom to approximate more accurately grain shape evolution. In addition, orientations θ can also be added as state variables with an associated dissipative mechanism related to crystal plasticity in the grain bulk. Thus, reorientation is likely to occur for small grain sizes as suggested by molecular dynamics  because reorientation involves dissipated energy per unit volume although grain boundary migration involves dissipated energy per unit area.
Grain growth is a thermally activated process as captured by the Arrhenius law (62). This is illustrated by simulating the same OT as those presented in Figure 10 at 700 °C. The resulting evolution is much slower than at 800 °C as shown in Figure 13.
Disorientation distributions of three OT at 800 °C.
Grain growth of the OT1 at 800 °C during approximately 4 h.
Grain growth of the OT2 at 800 °C during approximately 4 h.
Grain growth of the OT3 at 800 °C during approximately 4 h.
Grain growth of the OT1 at 700 °C during approximately 4 h.
In this paper the Orientated Tessellation Updating Method has been introduced to simulate grain growth at mesoscale (i.e., at the scale of the polycrystalline structure). The method relies on Voronoi-Laguerre tessellation techniques and energetic contributions at microscale (i.e., at the scale of the grain boundary), namely the grain boundary energy due to the crystal disorientation and the dissipated power by crystal plasticity for any virtual motion of the grain boundary. The method is very light in terms of computation cost. Even though 2D structures have been studied in this contribution, the extension in 3D is straightforward and future works will focus on this aspect. In addition, the only state variables are the weights of the tessellation. Thus, further developments including seeds and orientations as state variables are needed to approximate more accurately structural evolution of polycrystals. Furthermore, tessellation techniques include options to generate subgrain structures that could also be introduced in the framework presented in this paper.
- D. Kandel, E. Domany, Rigorous derivation of domain growth kinetics without conservation laws, J. Stat. Phys. 58, 685–706 (1990) [Google Scholar]
- E.A. Holm, G.N. Hassold, M.A. Miodownik, On misorientation distribution evolution during anisotropic grain growth, Acta Mater. 49, 2981–2991 (2001) [Google Scholar]
- E.A. Holm, M.A. Miodownik, A.D. Rollett, On abnormal subgrain growth and the origin of recrystallization nuclei, Acta Mater. 51, 2701–2716 (2003) [Google Scholar]
- L. Zhang, A.D. Rollett, T. Bartel, D. Wu, M.T. Lusk, A calibrated monte carlo approach to quantify the impacts of misorientation and different driving forces on texture development, Acta Mater. 60, 1201–1210 (2012) [Google Scholar]
- J. Gruber, D.C. George, A.P. Kuprat, G.S. Rohrer, A.D. Rollett, Effect of anisotropic grain boundary properties on grain boundary plane distributions during grain growth. Scr. Mater. 53, 351–355 (2005) [Google Scholar]
- H. Hallberg, Influence of anisotropic grain boundary properties on the evolution of grain boundary character distribution during grain growth-a 2d level set study, Model. Simul. Mater. Sci. Eng. 22, 085005 (2014) [Google Scholar]
- B. Scholtes, R. Boulais-Sinou, A. Settefrati, D.P. Muñoz, I. Poitrault, A. Montouchet, N. Bozzolo, M. Bernacki, 3d level set modeling of static recrystallization considering stored energy fields. Comput. Mater. Sci. 122, 57–71 (2016) [Google Scholar]
- N. Ma, A. Kazaryan, S.A. Dregia, Y. Wang, Computer simulation of texture evolution during grain growth: effect of boundary properties and initial microstructure, Acta Mater. 52, 3869–3879 (2004) [Google Scholar]
- C.E. Krill Iii, L.-Q. Chen, Computer simulation of 3-d grain growth using a phase-field model, Acta Mater. 50, 3059–3075 (2002) [Google Scholar]
- L. Vanherpe, N. Moelans, B. Blanpain, S. Vandewalle, Bounding box framework for efficient phase field simulation of graingrowth in anisotropic systems, Comput. Mater. Sci. 50, 2221–2231 (2011) [Google Scholar]
- K. Chang, N. Moelans, Effect of grain boundary energy anisotropy on highly textured grain structures studied by phase-field simulations, Acta Mater. 64, 443–454 (2014) [Google Scholar]
- M. Upmanyu, D.J. Srolovitz, L.S. Shvindlerman, G. Gottstein, Misorientation dependence of intrinsic grain boundary mobility: simulation and experiment, Acta Mater. 47, 3901–3914 (1999) [Google Scholar]
- M. Upmanyu, D.J. Srolovitz, L.S. Shvindlerman, G. Gottstein, Molecular dynamics simulation of triple junction migration, Acta Mater. 50, 1405–1420 (2002) [Google Scholar]
- M. Upmanyu, D.J. Srolovitz, A.E. Lobkovsky, J.A. Warren, W.C. Carter, Simultaneous grain boundary migration and grain rotation, Acta Mater. 54, 1707–1719 (2006) [Google Scholar]
- F.J. Humphreys, Modelling mechanisms and microstructures of recrystallisation, Mater. Sci. Technol. 8, 135–144 (1992) [Google Scholar]
- F. Wakai, N. Enomoto, H. Ogawa, Three-dimensional microstructural evolution in ideal grain growth-general statistics, Acta Mater. 48, 1297–1311 (2000) [Google Scholar]
- M. Syha, D. Weygand, A generalized vertex dynamics model for grain growth in three dimensions, Model. Simul. Mater. Sci. Eng. 18, 015010 (2009) [Google Scholar]
- A. Vondrous, M. Reichardt, B. Nestler. Growth rate distributions for regular two-dimensional grains with read–shockley grain boundary energy, Model. Simul. Mater. Sci. Eng. 22, 025014 (2014) [Google Scholar]
- W.T. Read, W. Shockley, Dislocation models of crystal grain boundaries, Phys. Rev. 78, 275 (1950) [Google Scholar]
- D. Wolf, A read-shockley model for high-angle grain boundaries, Scr. Metal. 23, 1713–1718 (1989) [Google Scholar]
- V.V. Bulatov, B.W. Reed, M. Kumar, Grain boundary energy function for fcc metals, Acta Mater. 65, 161–175 (2014) [Google Scholar]
- A. Kagawa, T. Okamoto, H. Matsumoto, Young’s modulus and thermal expansion of pure iron-cementite alloy castings, Acta Metal. 35, 797–803 (1987) [Google Scholar]
- T. Cheng, D. Fang, Y. Yang, The temperature dependence of grain boundary free energy of solids, J. Appl. Phys. 123, 085902 (2018) [Google Scholar]
- M. Brocato, A. Ehrlacher, P. Tamagny, Détermination de la dissipation caractéristique dans la propagation d’un front de recristallisation, C.R. Acad. Sci. Ser. IIB Mech. Phys. Astron. 327, 179–184 (1999) [Google Scholar]
- M. Brocato, A. Ehrlacher, P. Tamagny, Stability of discontinuities in polycrystals, Waves and Stability in Continuous Media, World Scientific, Singapore, 1999, p. 57 [Google Scholar]
- K. Hackl, F.D. Fischer, On the relation between the principle of maximum dissipation and inelastic evolution given by dissipation potentials, Proc. R. Soc. A 464, 117–132 (2007) [Google Scholar]
- Scilab. Scilab: Free and open source software for numerical computation, Scilab Enterprises, Orsay, France, 2012 [Google Scholar]
- C.W. Sinclair, C.R. Hutchinson, Y. Brechet, The effect of nb on the recrystallization and grain growth of ultra-high-purity α-fe: a combinatorial approach, Metal. Mater. Trans. A 38, 821–830 (2007) [Google Scholar]
- R. Quey, P.R. Dawson, F. Barbe, Large-scale 3d random polycrystals for the finite element method: Generation, meshing and remeshing, Comput. Methods Appl. Mech. Eng. 200, 1729–1745 (2011) [Google Scholar]
- C. Kerisit, R.E. Logé, S. Jacomet, V. Llorca, N. Bozzolo, Ebsd coupled to sem in situ annealing for assessing recrystallization and grain growth mechanisms in pure tantalum, J. Microsc. 250, 189–199 (2013) [PubMed] [Google Scholar]
Cite this article as: D. Weisz-Patrault, S. Sakout, A. Ehrlacher, Fast simulation of grain growth based on Orientated Tessellation Updating Method, Mechanics & Industry 21, 513 (2020)
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.