摘要
A fully implicit numerical method,based upon a combination of adaptively refined hierarchical meshes and geometric multigrid,is presented for the simulation of binary alloy solidification in three space dimensions.The computational techniques are presented for a particular mathematical model,based upon the phase-field approach,however their applicability is of greater generality than for the specific phase-field model used here.In particular,an implicit second order time discretization is combined with the use of second order spatial differences to yield a large nonlinear system of algebraic equations as each time step.It is demonstrated that these equations may be solved reliably and efficiently through the use of a nonlinear multigrid scheme for locally refined grids.In effect this paper presents an extension of earlier research in two space dimensions(J.Comput.Phys.,225(2007),pp.1271–1287)to fully three-dimensional problems.This extension is validated against earlier two-dimensional results and against some of the limited results available in three dimensions,obtained using an explicit scheme.The efficiency of the implicit approach and the multigrid solver are then demonstrated and some sample computational results for the simulation of the growth of dendrite structures are presented.
基金
The authors would like to thank Jan Rosam and James Green for their contributions to the development of the implicit software used in this work.The PARAMESH software used in this work was developed at the NASA Goddard Space Flight Center and Drexel University under NASA’s HPCC and ESTO/CT projects and under grant NNG04GP79G from the NASA/AISR project.In addition we wish to thank Jonathan Dantzig for providing the explicit software used for comparison.The parallel simulations have been run on HECToR,the UK’s National Supercomputing Service,for which we gratefully acknowledge the compute time provided to us.