Extreme-Scale Phase Field Simulations of Coarsening Dynamics on the Sunway TaihuLight Supercomputer

Many important properties of materials such as strength, ductility, hardness and conductivity are determined by the microstructures of the material. During the formation of these microstructures, grain coarsening plays an important role. The Cahn-Hilliard equation has been applied extensively to simulate the coarsening kinetics of a two-phase microstructure. It is well accepted that the limited capabilities in conducting large scale, long time simulations constitute bottlenecks in predicting microstructure evolution based on the phase field approach. We present here a scalable time integration algorithm with large stepsizes and its efficient implementation on the Sunway TaihuLight supercomputer. The highly nonlinear and severely stiff Cahn-Hilliard equations with degenerate mobility for microstructure evolution are solved at extreme scale, demonstrating that the latest advent of high performance computing platform and the new advances in algorithm design are now offering us the possibility to simulate the coarsening dynamics accurately at unprecedented spatial and time scales.

[1]  Charles M. Elliott,et al.  The Cahn–Hilliard equation with a concentration dependent mobility: motion by minus the Laplacian of the mean curvature , 1996, European Journal of Applied Mathematics.

[2]  T. Takaki,et al.  GPU-accelerated phase-field simulation of dendritic solidification in a binary alloy , 2011 .

[3]  Marlis Hochbruck,et al.  Exponential Integrators for Large Systems of Differential Equations , 1998, SIAM J. Sci. Comput..

[4]  Chao Yang,et al.  A parallel domain decomposition-based implicit method for the Cahn-Hilliard-Cook phase-field equation in 3D , 2015, J. Comput. Phys..

[5]  J. E. Hilliard,et al.  Free Energy of a Nonuniform System. I. Interfacial Free Energy , 1958 .

[6]  Robert L. Pego,et al.  Front migration in the nonlinear Cahn-Hilliard equation , 1989, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[7]  Y. Saad.  Analysis of some Krylov subspace approximations to the matrix exponential operator , 1992 .

[8]  Desai,et al.  Numerical study of the late stages of spinodal decomposition. , 1988, Physical review. B, Condensed matter.

[9]  Jack J. Dongarra,et al.  A Portable Programming Interface for Performance Evaluation on Modern Processors , 2000, Int. J. High Perform. Comput. Appl..

[10]  Desai,et al.  Numerical study of late-stage coarsening for off-critical quenches in the Cahn-Hilliard equation of phase separation. , 1989, Physical review. B, Condensed matter.

[11]  Roger B. Sidje,et al.  Expokit: a software package for computing matrix exponentials , 1998, TOMS.

[12]  Marco Vianello,et al.  A massively parallel exponential integrator for advection-diffusion models , 2009, J. Comput. Appl. Math..

[13]  Ulrich Rüde,et al.  Massively parallel phase-field simulations for ternary eutectic directional solidification , 2015, SC15: International Conference for High Performance Computing, Networking, Storage and Analysis.

[14]  Jian Zhang,et al.  Fast and accurate algorithms for simulating coarsening dynamics of Cahn–Hilliard equations , 2015 .

[15]  C. Lubich,et al.  On Krylov Subspace Approximations to the Matrix Exponential Operator , 1997 .

[16]  Long-Qing Chen.  Phase-Field Models for Microstructure Evolution , 2002 .

[17]  Peter W Voorhees,et al.  Phase-field simulation of 2-D Ostwald ripening in the high volume fraction regime , 2002 .

[18]  T. Hughes,et al.  Isogeometric analysis of the Cahn–Hilliard phase-field model , 2008 .

[19]  Weidong Zhao,et al.  Fast High-Order Compact Exponential Time Differencing Runge–Kutta Methods for Second-Order Semilinear Parabolic Equations , 2016, J. Sci. Comput..

[20]  M. Hochbruck,et al.  Exponential integrators , 2010, Acta Numerica.

[21]  T. Takaki,et al.  Unexpected selection of growing dendrites by very-large-scale phase-field simulation , 2013 .

[22]  Milo R. Dorr,et al.  A numerical algorithm for the solution of a phase-field model of polycrystalline materials , 2008, J. Comput. Phys..

[23]  Qing Nie,et al.  Efficient semi-implicit schemes for stiff systems , 2006, J. Comput. Phys..

[24]  J. Cahn,et al.  A microscopic theory for antiphase boundary motion and its application to antiphase domain coasening , 1979 .

[25]  F. Otto,et al.  Upper Bounds on Coarsening Rates , 2002 .

[26]  S. Zwaag,et al.  Three-dimensional phase field modelling of the austenite-to-ferrite transformation , 2006 .

[27]  Qiang Du,et al.  Coarsening Mechanism for Systems Governed by the Cahn-Hilliard Equation with Degenerate Diffusion Mobility , 2014, Multiscale Model. Simul..

[28]  Emmott,et al.  Lifshitz-Slyozov scaling for late-stage coarsening with an order-parameter-dependent mobility. , 1995, Physical review. B, Condensed matter.

[29]  Harald Garcke,et al.  Transient coarsening behaviour in the Cahn–Hilliard model , 2003 .

[30]  S. Cox,et al.  Exponential Time Differencing for Stiff Systems , 2002 .

[31]  Jian Zhang,et al.  Fast Explicit Integration Factor Methods for Semilinear Parabolic Equations , 2015, J. Sci. Comput..

[32]  I. Moret,et al.  An interpolatory approximation of the matrix exponential based on Faber polynomials , 2001 .

[33]  T. Takaki,et al.  Two-dimensional phase-field simulations of dendrite competitive growth during the directional solidification of a binary alloy bicrystal , 2014 .

[34]  S. Krogstad.  Generalized integrating factor methods for stiff PDEs , 2005 .

[35]  T. Küpper,et al.  Simulation of particle growth and Ostwald ripening via the Cahn-Hilliard equation , 1994 .

[36]  Hui Lv,et al.  Cooperative Computing Techniques for a Deeply Fused and Heterogeneous Many-Core Processor Architecture , 2015, Journal of Computer Science and Technology.

[37]  Constantine Bekas,et al.  An extreme-scale implicit solver for complex PDEs: highly heterogeneous flow in earth's mantle , 2015, SC15: International Conference for High Performance Computing, Networking, Storage and Analysis.

[38]  Qiang Du,et al.  Motion of Interfaces Governed by the Cahn-Hilliard Equation with Highly Disparate Diffusion Mobility , 2012, SIAM J. Appl. Math..

[39]  Qiang Du,et al.  Analysis and Applications of the Exponential Time Differencing Schemes and Their Contour Integration Modifications , 2005 .

[40]  Olga Wodo,et al.  Computationally efficient solution to the Cahn-Hilliard equation: Adaptive implicit time schemes, mesh sensitivity analysis and the 3D isoperimetric problem , 2011, J. Comput. Phys..

[41]  Satoshi Matsuoka,et al.  Peta-scale phase-field simulation for dendritic solidification on the TSUBAME 2.0 supercomputer , 2011, 2011 International Conference for High Performance Computing, Networking, Storage and Analysis (SC).

[42]  Mahmood Mamivand,et al.  A review on phase field modeling of martensitic phase transformation , 2013 .

[43]  Jie Shen,et al.  Coarsening kinetics from a variable-mobility Cahn-Hilliard equation: application of a semi-implicit Fourier spectral method. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[44]  Cleve B. Moler,et al.  Nineteen Dubious Ways to Compute the Exponential of a Matrix, Twenty-Five Years Later , 1978, SIAM Rev..