A combined finite element and finite difference scheme for computer simulation of microstructure evolution and its application to pore-boundary separation during sintering
S. Kucherenko et al., A combined finite element and finite difference scheme for computer simulation of microstructure evolution and its application to pore-boundary separation during sintering, COMP MAT SC, 18(1), 2000, pp. 76-92
This paper presents two new developments to the numerical technique previou
sly developed by Pan and Cocks [cf. J. Pan, A.C.F. Cocks, A numerical techn
ique for the analysis of coupled surface and grain-boundary diffusion, Acta
Metall. 43 (1995) 1395-1406; J. Pan, A.C.F. Cocks, S. Kucherenko, Finite e
lement analysis of coupled grain-boundary and surface diffusion with grain-
boundary migration, Proc. R. Sec. London A 453 (1998) 2161-2184] for the co
mputer simulation of microstructure evolution of materials. Coupled grain-b
oundary diffusion, surface diffusion and grain-boundary migration are consi
dered as the underlying mechanisms for the evolution. The first development
is that a set of "link elements" are developed to link the finite differen
ce scheme of Pan and Cocks (1995) with the finite element scheme of Pan et
al. (1998) for the surface diffusion and grain-boundary migration parts of
the problem respectively. Unlike the method used in Pan and Cocks (1995), t
hese link elements are designed to link the two discretisation schemes away
from the interface junction so that the dihedral angle can be maintained i
n the variational sense at the junction while the finite difference scheme
can still be used for most of the interface network. Such a combined scheme
is more efficient than the full finite element scheme because most of the
degrees of freedom for surface diffusion and grain-boundary migration no lo
nger contribute to the global linear simultaneous equations. The second dev
elopment is that an implicit time integration method is implemented. In gen
eral the implicit time integration method allows much larger timesteps to b
e used than that allowed by an explicit method. The two new developments to
gether significantly improved the efficiency of the numerical scheme. Sever
al test cases are provided to verify the numerical scheme. As an example of
application, the effect of pore shrinkage on pore-boundary separation is i
nvestigated using the numerical scheme. It is shown that the existing separ
ation criteria significantly over-predict separation. (C) 2000 Elsevier Sci
ence B.V. All rights reserved.