An adaptive, multilevel scheme for the implicit solution of three‐dimensional phase‐field equations

Phase-field models, consisting of a set of highly nonlinear coupled parabolic partial differential equations, are widely used for the simulation of a range of solidification phenomena. This article focuses on the numerical solution of one such model, representing anisotropic solidification in three space dimensions. The main contribution of the work is to propose a solution strategy that combines hierarchical mesh adaptivity with implicit time integration and the use of a nonlinear multigrid solver at each step. This strategy is implemented in a general software framework that permits parallel computation in a natural manner. Results are presented that provide both qualitative and quantitative justifications for these choices.

[1]  N. Goldenfeld,et al.  Adaptive Mesh Refinement Computation of Solidification Microstructures Using Dynamic Data Structures , 1998, cond-mat/9808216.

[2]  K. Olson,et al.  PARAMESH: A Parallel, Adaptive Grid Tool , 2006 .

[3]  Peter K. Jimack,et al.  A fully implicit, fully adaptive time and space discretisation method for phase-field simulation of binary alloy solidification , 2007, J. Comput. Phys..

[4]  Andrew M. Mullis,et al.  Quantification of mesh induced anisotropy effects in the phase-field method , 2006 .

[5]  Hidehiro Onodera,et al.  Parallel Computer Simulation of Three-Dimensional Grain Growth Using the Multi-Phase-Field Model , 2008 .

[6]  R. Cochrane,et al.  A phase field model for spontaneous grain refinement in deeply undercooled metallic melts , 2001 .

[7]  R. Cochrane,et al.  Microstructural evolution and growth velocity-undercooling relationships in the systems Cu, Cu-O and Cu-Sn at high undercooling , 2000 .

[8]  Ruo Li,et al.  A multi-mesh adaptive finite element approximation to phase field models , 2009 .

[9]  Wheeler,et al.  Phase-field models for anisotropic interfaces. , 1993, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[10]  Michael Griebel,et al.  Parallel multigrid in an adaptive PDE solver based on hashing and space-filling curves , 1999, Parallel Comput..

[11]  Martin Berzins,et al.  A comparison of some dynamic load-balancing algorithms for a parallel adaptive flow solver , 2000, Parallel Comput..

[12]  William L. George,et al.  A Parallel 3D Dendritic Growth Simulator Using the Phase-Field Method , 2002 .

[13]  A. Mullis.  Effect of the ratio of solid to liquid conductivity on the stability parameter of dendrites within a phase-field model of solidification. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[14]  N. Goldenfeld,et al.  Phase field model for three-dimensional dendritic growth with fluid flow. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[15]  R. Cochrane,et al.  Deformation of dendrites by fluid flow during rapid solidification , 2001 .

[16]  A. Karma,et al.  Quantitative phase-field modeling of dendritic growth in two and three dimensions , 1996 .

[17]  A. Karma,et al.  Phase-field method for computationally efficient modeling of solidification with arbitrary interface kinetics. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[18]  Brener.  Needle-crystal solution in three-dimensional dendritic growth. , 1993, Physical review letters.

[19]  Roy D. Williams,et al.  Performance of dynamic load balancing algorithms for unstructured mesh calculations , 1991, Concurr. Pract. Exp..

[20]  J. Warren,et al.  Prediction of dendritic growth and microsegregation patterns in a binary alloy using the phase-field method , 1995 .

[21]  P. Jimack,et al.  Quantitative phase-field modeling of solidification at high Lewis number. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[22]  Peter K. Jimack,et al.  An adaptive, fully implicit multigrid phase-field model for the quantitative simulation of non-isothermal binary alloy solidification , 2008 .

[23]  Wouter-Jan Rappel,et al.  Phase-field simulation of three-dimensional dendrites: is microscopic solvability theory correct? , 1997 .

[24]  Stefan Lang,et al.  Parallel adaptive multigrid methods in plane linear elasticity problems , 1997, Numer. Linear Algebra Appl..

[25]  D. Brandt,et al.  Multi-level adaptive solutions to boundary-value problems math comptr , 1977 .

[26]  Peter K. Jimack,et al.  An adaptive multigrid tool for elliptic and parabolic systems , 2005 .

[27]  Barbieri,et al.  Predictions of dendritic growth rates in the linearized solvability theory. , 1989, Physical review. A, General physics.