research-article

Extreme-scale phase field simulations of coarsening dynamics on the sunway taihulight supercomputer

Authors:

Zhao LiuAuthors Info & Claims

SC '16: Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis

Article No.: 4, Pages 1 - 12

Published: 13 November 2016 Publication History

Abstract

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 step-sizes 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.

References

[1]

L.-Q. Chen, "Phase-field models for microstructure evolution," Annu. Rev. Mater. Res., vol. 32, pp. 113--140, 2002.

[2]

J. Cahn and J. Hillard, "Free energy of a nonuniform system i. interfacial free energy," J. Chem. Phys., vol. 28, pp. 258--267, 1958.

[3]

S. Allen and J. W. Cahn, "A microscopic theory for antiphase boundary motion and its application to antiphase domain coarsening," Acta Metall., vol. 27, pp. 1084--1095, 1979.

[4]

J. Cahn, C. Elliott, and A. Novick-Cohen, "The Cahn-Hilliard equation with a concentration-dependent mobility: motion by minus the Laplacian of the mean curvature," European J. Appl. Math., vol. 7, pp. 287--301, 1996.

[5]

S. Dai and Q. Du, "Motion of interfaces governed by the Cahn-Hilliard equation with highly disparate diffusion mobility," SIAM J. Appl. Math., vol. 72, pp. 1818--1841, 2012.

[6]

S. Dai and Q. Du, "Coarsening mechanism for systems governed by the Cahn-Hilliard equation with degenerate diffusion mobility," Multi. Model. Simul., vol. 12, pp. 1870--1889, 2014.

[7]

R. Pego., "Front migration in the nonlinear Cahn-Hilliard equation," Proc. Royal Soc. London A, vol. 442, pp. 261--278, 1989.

[8]

A. Bray and C. Emmott, "Lifshitz-slyozov scaling for late-stage coarsening with an order parameter-dependent mobility," Phys. Rev. B, vol. 52, pp. R685--R688, 1995.

[9]

D. Fan, S. Chen, L.-Q. Chen, and P. Voorhees, "Phase-field simulation of 2-d Ostwald ripening in the high volume fraction regime," Acta Materialia, vol. 50, pp. 1895--1907, 2002.

[10]

H. Garcke, B. Niethammer, M. Rumpf, and U. Weikard, "Transient coarsening behaviour in the Cahn-Hilliard model," Acta Materialia, vol. 51, pp. 2823--2830, 2003.

[11]

J. Zhu, L.-C. Chen, J. Shen and V. Tikare, "Coarsening kinetics from a variable-mobility Cahn-Hilliard equation: Application of a semi-implicit fourier spectral method," Phy. Rev. E, vol. 60, pp. 3564--3572, 1999.

[12]

T. Kupper and N. Masbaum, "Simulation of particle growth and ostwald ripening via the Cahn-Hilliard equation," Acta.Metall. Mater., vol. 42, pp. 1847--1858, 1994.

[13]

T. Rogers and R. Desai, "Numerical study of late-stage coarsening for off-critical quenches in the Cahn-Hilliard equation of phase separation," Phys. Rev. B, vol. 39, pp. 11 956--11 964, 1989.

[14]

T. Rogers, K. Elder, and R. Desai, "Numerical study of the late stages of spinodal decomposition," Phys. Rev. B, vol. 37, pp. 9638--9649, 1988.

[15]

R. Kohn and F. Otto, "Upper bounds on coarsening rates," Comm. Math. Phys., vol. 229, pp. 375--395, 2002.

[16]

H. Gómez, V. M. Calo, Y. Bazilevs, and T. J. Hughes, "Isogeometric analysis of the cahn-hilliard phase-field model," Computer methods in applied mechanics and engineering, vol. 197, no. 49, pp. 4333--4352, 2008.

[17]

X. Zheng, C. Yang, X.-C. Cai, and D. Keyes, "A parallel domain decomposition-based implicit method for the Cahn-Hilliard-Cook phase-field equation in 3d," Journal of Computational Physics, vol. 285, pp. 55--70, 2015.

Digital Library

[18]

M. Mamivand, M. A. Zaeem, and H. El Kadiri, "A review on phase field modeling of martensitic phase transformation," Computational Materials Science, vol. 77, pp. 304--311, 2013.

[19]

M. Militzer, M. Mecozzi, J. Sietsma, and S. Van der Zwaag, "Three-dimensional phase field modelling of the austenite-to-ferrite transformation," Acta materialia, vol. 54, no. 15, pp. 3961--3972, 2006.

[20]

L. Ju, J. Zhang, and Q. Du, "Fast and accurate algorithms for simulating coarsening dynamics of Cahn-Hilliard equations," Computational Materials Science, vol. 108, pp. 272--282, 2015.

[21]

L. Ju, J. Zhang, L. Zhu, and Q. Du, "Fast explicit integration factor methods for semilinear parabolic equations," J. Sci. Comput., vol. 62, pp. 431--455, 2015.

Digital Library

[22]

T. Shimokawabe, T. Aoki, T. Takaki, A. Yamanaka, A. Nukada, T. Endo, N. Maruyama, and S. M. Clark, "Peta-scale phase-field simulation for dendritic solidification on the tsubame 2.0 supercomputer," in Proceedings of 2011 International Conference for High Performance Computing, Networking, Storage and Analysis (SC '11). Article No. 3, 2011.

Digital Library

[23]

A. Yamanaka, T. Aoki, S. Ogawa, and T. Takaki, "Gpu-accelerated phase-field simulation of dendritic solidification in a binary alloy," European J. Appl. Math., vol. 381(1), pp. 40--45, 2011.

[24]

T. Takaki, M. Ohno, T. Shimokawabe, and T. Aoki, "Two-dimensional phase-field simulations of dendrite competitive growth during the directional solidification of a binary alloy bicrystal," Acta Materialia, vol. 81(0), pp. 272--283, 2014.

[25]

T. Takaki, T. Shimokawabe, M. Ohno, A. Yamanaka, and T. Aoki, "Unexpected selection of growing dendrites by very-large-scale phase-field simulation," Journal of Crystal Growth, vol. 382, pp. 21--25, 2013.

[26]

M. Bauer, J. Hötzer, M. Jainta, p. Steinmetz, M. Berghoff, F. Schornbaum, C. Godenschwager, H. Köstler, B. Nestler, and U. Rüde, "Massively parallel phase-field simulations for ternary eutectic directional solidification," Proceedings of 2015 International Conference for High Performance Computing, Networking, Storage and Analysis (SC '15), 2015.

Digital Library

[27]

M. Dorr, J.-L. Fattebert, M. Wickett, J. Belak, and P. Turchi, "A numerical algorithm for the solution of a phase-field model of polycrystalline materials," Journal of Computational Physics, vol. 229, pp. 626--641, 2010.

Digital Library

[28]

O. Wodoa and B. Ganapathysubramanian, "Computationally effcient solution to the cahn-hilliard equation: adaptive implicit time schemes, mesh sensitivity analysis and the 3d isoperimetric problem," Journal of Computational Physics, vol. 230, pp. 6037--6060, 2011.

Digital Library

[29]

J. Rudi, A. C. I. Malossiy, T. Isaac, G. Stadlerk, M. Gurnisz, P. W. J. Staary, Y. Ineicheny, C. Bekasy, A. Curioniy, and O. Ghattas, "An extreme-scale implicit solver for complex pdes: Highly heterogeneous flow in earth's mantle," Proceedings of 2015 International Conference for High Performance Computing, Networking, Storage and Analysis (SC '15), 2015.

Digital Library

[30]

M. Hochbruck, C. Lubich, and H. Selhofer, "Exponential integrators for large systems of differential equations," SIAM J. Sci. Comput., vol. 19, pp. 1552--1574, 1998.

Digital Library

[31]

S. Cox and P. Matthews, "Exponential time differencing for stiff systems," J. Comput. Phys., vol. 176, pp. 430--455, 2002.

Digital Library

[32]

Q. Du and W. Zhu, "Stability analysis and applications of the exponential time differencing schemes and their contour integration modifications," J. Comput. Math., vol. 22, pp. 200--209, 2004.

[33]

S. Krogstad, "Generalized integrating factor methods for stiff pdes," J. Comput. Phys., vol. 203, pp. 72--88, 2005.

Digital Library

[34]

Q. Du and W. Zhu, "Analysis and applications of the exponential time differencing schemes," BIT Numer. Math., vol. 45, pp. 307--328, 2005.

Digital Library

[35]

Q. Nie, Y.-T. Zhang, and R. Zhao, "Efficient semi-implicit schemes for stiff systems," J. Comput. Phys., vol. 214, pp. 521--537, 2006.

Digital Library

[36]

M. Hochbruck and A. Ostermann, "Exponential integrators," Acta Numerica, vol. 19, pp. 209--286, 2010.

[37]

C. Moler and C. V. Loan, "Nineteen dubiousways to compute the exponential of a matrix, twenty-five years later," SIAM REVIEW, vol. 45, no. 1, pp. 3--49, 2003.

Digital Library

[38]

L. Zhu, L. Ju, and W. Zhao, "Fast high-order compact exponential time differencing Runge-Kutta methods for second-order semilinear parabolic equations," Journal of Scientific Computing, vol. 67, no. 3, pp. 1043--1065, 2016.

Digital Library

[39]

Y. Saad, "Analysis of some krylov subspace approximations to the matrix exponential operator," SIAM Journal of Numerical Analysis, vol. 29, pp. 209--228, 1992.

Digital Library

[40]

M. Hochbruck and C. Lubich, "On krylov subspace approximations to the matrix exponential," SIAM J. Numer. Anal., vol. 34, pp. 1911--1925, 1997.

Digital Library

[41]

R. Sidje, "Expokit: A software package for computing matrix exponentials," ACM Trans. Math. Software, vol. 24, pp. 130--156, 1998.

Digital Library

[42]

I. Moret and P. Novati, "An interpolatory approximation of the matrix exponential based on faber polynomials," J. Comput. Appl. Math., vol. 131, pp. 361--380, 2001.

Digital Library

[43]

A. Martinez, L. Bergamaschi, M. Caliari, and M. Vianello, "A massively parallel exponential integrator for advection-diffusion models," Journal of Computational and Applied Mathematics, vol. 231, pp. 82--91, 2009.

Digital Library

[44]

F. Zheng, H. Li, and H. L. et al., "Cooperative computing techniques for a deeply fused and heterogeneous many-core processor architecture," JOURNAL OF COMPUTER SCIENCE AND TECHNOLOGY, vol. 30(1), pp. 145--162, 2015.

[45]

S. Browne, J. Dongarra, N. Garner, G. Ho, and P. Mucci, "A portable programming interface for performance evaluation on modern processors," Int. J. High Perform. Comput. Appl, vol. 14, no. 3, pp. 189--204, 2000.

Digital Library

Cited By

Li XYe HZhang J(2020)Large-scale Simulations of Peridynamics on Sunway Taihulight SupercomputerProceedings of the 49th International Conference on Parallel Processing10.1145/3404397.3404421(1-11)Online publication date: 17-Aug-2020
https://dl.acm.org/doi/10.1145/3404397.3404421
Jiang LYang CMa W(2020)Enabling Highly Efficient Batched Matrix Multiplications on SW26010 Many-core ProcessorACM Transactions on Architecture and Code Optimization10.1145/337817617:1(1-23)Online publication date: 4-Mar-2020
https://dl.acm.org/doi/10.1145/3378176
Sun QZhang CWu CZhang JLi L(2018)Bandwidth Reduced Parallel SpMV on the SW26010 Many-Core PlatformProceedings of the 47th International Conference on Parallel Processing10.1145/3225058.3225074(1-10)Online publication date: 13-Aug-2018
https://dl.acm.org/doi/10.1145/3225058.3225074
Show More Cited By

Recommendations

Redesigning CAM-SE for peta-scale climate modeling performance and ultra-high resolution on Sunway TaihuLight
SC '17: Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis

The Community Atmosphere Model (CAM) is ported, redesigned, and scaled to the full system of the Sunway TaihuLight, and provides peta-scale climate modeling performance. We refactored and optimized the complete code using OpenACC directives at the first ...
Optimizing Convolutional Neural Networks on the Sunway TaihuLight Supercomputer

The Sunway TaihuLight supercomputer is powered by SW26010, a new 260-core processor designed with on-chip fusion of heterogeneous cores. In this article, we present our work on optimizing the training process of convolutional neural networks (CNNs) on ...
Performance Optimization of the HPCG Benchmark on the Sunway TaihuLight Supercomputer

In this article, we present some key techniques for optimizing HPCG on Sunway TaihuLight and demonstrate how to achieve high performance in memory-bound applications by exploiting specific characteristics of the hardware architecture. In particular, we ...

Comments

Information & Contributors

Information

Published In

cover image ACM Conferences

SC '16: Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis

November 2016

1034 pages

ISBN:9781467388153

Conference Chair:
John West
University of Texas at Austin

Sponsors

SIGARCH: ACM Special Interest Group on Computer Architecture
IEEE-CS: Computer Society

In-Cooperation

SIGHPC: ACM Special Interest Group on High Performance Computing, Special Interest Group on High Performance Computing

Publisher

IEEE Press

Publication History

Published: 13 November 2016

Check for updates

Qualifiers

Research-article

Conference

SC16

Sponsor:

SIGARCH
IEEE-CS

SC16: The International Conference for High Performance Computing, Networking, Storage and Analysis

November 13 - 18, 2016

Utah, Salt Lake City

Acceptance Rates

SC '16 Paper Acceptance Rate 81 of 442 submissions, 18%;

Overall Acceptance Rate 1,516 of 6,373 submissions, 24%

Contributors

Other Metrics

View Article Metrics

Bibliometrics & Citations

Bibliometrics

Article Metrics

6
Total Citations
View Citations
398
Total Downloads

Downloads (Last 12 months)15
Downloads (Last 6 weeks)1

Reflects downloads up to 30 Aug 2024

Other Metrics

View Author Metrics

Citations

Cited By

Li XYe HZhang J(2020)Large-scale Simulations of Peridynamics on Sunway Taihulight SupercomputerProceedings of the 49th International Conference on Parallel Processing10.1145/3404397.3404421(1-11)Online publication date: 17-Aug-2020
https://dl.acm.org/doi/10.1145/3404397.3404421
Jiang LYang CMa W(2020)Enabling Highly Efficient Batched Matrix Multiplications on SW26010 Many-core ProcessorACM Transactions on Architecture and Code Optimization10.1145/337817617:1(1-23)Online publication date: 4-Mar-2020
https://dl.acm.org/doi/10.1145/3378176
Sun QZhang CWu CZhang JLi L(2018)Bandwidth Reduced Parallel SpMV on the SW26010 Many-Core PlatformProceedings of the 47th International Conference on Parallel Processing10.1145/3225058.3225074(1-10)Online publication date: 13-Aug-2018
https://dl.acm.org/doi/10.1145/3225058.3225074
Cai YAo YYang CMa WZhao H(2018)Extreme-Scale High-Order WENO Simulations of 3-D Detonation Wave with 10 Million CoresACM Transactions on Architecture and Code Optimization10.1145/320920815:2(1-21)Online publication date: 12-Jun-2018
https://dl.acm.org/doi/10.1145/3209208
Liu CXie BLiu XXue WYang HLiu X(2018)Towards Efficient SpMV on Sunway Manycore ArchitecturesProceedings of the 2018 International Conference on Supercomputing10.1145/3205289.3205313(363-373)Online publication date: 12-Jun-2018
https://dl.acm.org/doi/10.1145/3205289.3205313
Cai YYang CMa WAo YEl-Araby EEl-Ghazawi TPanda D(2018)Extreme-scale realistic stencil computations on sunway taihulight with ten million coresProceedings of the 18th IEEE/ACM International Symposium on Cluster, Cloud and Grid Computing10.1109/CCGRID.2018.00086(566-571)Online publication date: 1-May-2018
https://dl.acm.org/doi/10.1109/CCGRID.2018.00086

View Options

Get Access

Login options

Check if you have access through your login credentials or your institution to get full access on this article.

Full Access

Get this Publication

View options

PDF

View or Download as a PDF file.

eReader

View online with eReader.

Media

Figures

Other

Tables

View Table of Contents