A numerical algorithm for the solution of a phase-field model of polycrystalline materials

被引：41

作者：

Dorr, M. R. ^{[1
]}

Fattebert, J. -L. ^{[1
]}

Wickett, M. E. ^{[1
]}

Belak, J. F. ^{[1
]}

Turchi, P. E. A. ^{[1
]}

机构：

[1] Lawrence Livermore Natl Lab, Livermore, CA 94550 USA

来源：

JOURNAL OF COMPUTATIONAL PHYSICS | 2010年 / 229卷 / 03期

关键词：

Phase-field model; Polycrystalline microstructure; Method of lines; Newton-Krylov methods; GRAIN-GROWTH; SIMULATION; SOLIDIFICATION; DIFFUSION; PLUTONIUM; EVOLUTION;

D O I：

10.1016/j.jcp.2009.09.041

中图分类号：

TP39 [计算机的应用];

学科分类号：

081203 ; 0835 ;

摘要：

We describe an algorithm for the numerical solution of a phase-field model (PFM) of microstructure evolution in polycrystalline materials. The PFM system of equations includes a local order parameter, a cluaternion representation of local orientation and a species composition parameter. The algorithm is based on the implicit integration of a semidiscretization of the PFM system using a backward difference formula (BDF) temporal discretization combined with a Newton-Krylov algorithm to solve the nonlinear system at each time step. The BDF algorithm is combined with a coordinate-projection method to maintain quaternion unit length, which is related to an important solution invariant. A key element of the Newton-Krylov algorithm is the selection of a preconclitioner to accelerate the convergence of the Generalized Minimum Residual algorithm used to solve the Jacobian linear system in each Newton step. Results are presented for the application of the algorithm to 2D and 3D examples. (C) 2009 Elsevier Inc. All rights reserved.

引用

页码：626 / 641

页数：16

共 36 条

[1] MICROSCOPIC THEORY FOR ANTIPHASE BOUNDARY MOTION AND ITS APPLICATION TO ANTIPHASE DOMAIN COARSENING [J].