A parallel domain decompositionbased implicit method for the Cahn–Hilliard–Cook phasefield equation in 3D
Abstract
We present a numerical algorithm for simulating the spinodal decomposition described by the three dimensional Cahn–Hilliard–Cook (CHC) equation, which is a fourthorder stochastic partial differential equation with a noise term. The equation is discretized in space and time based on a fully implicit, cellcentered finite difference scheme, with an adaptive timestepping strategy designed to accelerate the progress to equilibrium. At each time step, a parallel Newton–Krylov–Schwarz algorithm is used to solve the nonlinear system. We discuss various numerical and computational challenges associated with the method. The numerical scheme is validated by a comparison with an explicit scheme of high accuracy (and unreasonably high cost). We present steady state solutions of the CHC equation in two and three dimensions. The effect of the thermal fluctuation on the spinodal decomposition process is studied. We show that the existence of the thermal fluctuation accelerates the spinodal decomposition process and that the final steady morphology is sensitive to the stochastic noise. We also show the evolution of the energies and statistical moments. In terms of the parallel performance, it is found that the implicit domain decomposition approach scales well on supercomputers with a large number of processors.
 Authors:
 Department of Applied Physics and Applied Mathematics, Columbia University, New York, NY 10027 (United States)
 Institute of Software, Chinese Academy of Sciences, Beijing 100190 (China)
 (China)
 Department of Computer Science, University of Colorado Boulder, Boulder, CO 80309 (United States)
 Computer, Electrical and Mathematical Sciences and Engineering, King Abdullah University of Science and Technology, Thuwal 23955 (Saudi Arabia)
 (United States)
 Publication Date:
 OSTI Identifier:
 22465608
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Computational Physics; Journal Volume: 285; Other Information: Copyright (c) 2015 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; 97 MATHEMATICAL METHODS AND COMPUTING; ACCURACY; ALGORITHMS; FIELD EQUATIONS; FLUCTUATIONS; MATHEMATICAL SOLUTIONS; NONLINEAR PROBLEMS; PARTIAL DIFFERENTIAL EQUATIONS; STEADYSTATE CONDITIONS; STOCHASTIC PROCESSES; SUPERCOMPUTERS; THREEDIMENSIONAL CALCULATIONS
Citation Formats
Zheng, Xiang, Yang, Chao, State Key Laboratory of Computer Science, Chinese Academy of Sciences, Beijing 100190, Cai, XiaoChuan, Email: cai@cs.colorado.edu, Keyes, David, and Department of Applied Physics and Applied Mathematics, Columbia University, New York, NY 10027. A parallel domain decompositionbased implicit method for the Cahn–Hilliard–Cook phasefield equation in 3D. United States: N. p., 2015.
Web. doi:10.1016/J.JCP.2015.01.016.
Zheng, Xiang, Yang, Chao, State Key Laboratory of Computer Science, Chinese Academy of Sciences, Beijing 100190, Cai, XiaoChuan, Email: cai@cs.colorado.edu, Keyes, David, & Department of Applied Physics and Applied Mathematics, Columbia University, New York, NY 10027. A parallel domain decompositionbased implicit method for the Cahn–Hilliard–Cook phasefield equation in 3D. United States. doi:10.1016/J.JCP.2015.01.016.
Zheng, Xiang, Yang, Chao, State Key Laboratory of Computer Science, Chinese Academy of Sciences, Beijing 100190, Cai, XiaoChuan, Email: cai@cs.colorado.edu, Keyes, David, and Department of Applied Physics and Applied Mathematics, Columbia University, New York, NY 10027. 2015.
"A parallel domain decompositionbased implicit method for the Cahn–Hilliard–Cook phasefield equation in 3D". United States.
doi:10.1016/J.JCP.2015.01.016.
@article{osti_22465608,
title = {A parallel domain decompositionbased implicit method for the Cahn–Hilliard–Cook phasefield equation in 3D},
author = {Zheng, Xiang and Yang, Chao and State Key Laboratory of Computer Science, Chinese Academy of Sciences, Beijing 100190 and Cai, XiaoChuan, Email: cai@cs.colorado.edu and Keyes, David and Department of Applied Physics and Applied Mathematics, Columbia University, New York, NY 10027},
abstractNote = {We present a numerical algorithm for simulating the spinodal decomposition described by the three dimensional Cahn–Hilliard–Cook (CHC) equation, which is a fourthorder stochastic partial differential equation with a noise term. The equation is discretized in space and time based on a fully implicit, cellcentered finite difference scheme, with an adaptive timestepping strategy designed to accelerate the progress to equilibrium. At each time step, a parallel Newton–Krylov–Schwarz algorithm is used to solve the nonlinear system. We discuss various numerical and computational challenges associated with the method. The numerical scheme is validated by a comparison with an explicit scheme of high accuracy (and unreasonably high cost). We present steady state solutions of the CHC equation in two and three dimensions. The effect of the thermal fluctuation on the spinodal decomposition process is studied. We show that the existence of the thermal fluctuation accelerates the spinodal decomposition process and that the final steady morphology is sensitive to the stochastic noise. We also show the evolution of the energies and statistical moments. In terms of the parallel performance, it is found that the implicit domain decomposition approach scales well on supercomputers with a large number of processors.},
doi = {10.1016/J.JCP.2015.01.016},
journal = {Journal of Computational Physics},
number = ,
volume = 285,
place = {United States},
year = 2015,
month = 3
}

Phase transitions of macromolecular microsphere composite hydrogels based on the stochastic Cahn–Hilliard equation
We use the stochastic Cahn–Hilliard equation to simulate the phase transitions of the macromolecular microsphere composite (MMC) hydrogels under a random disturbance. Based on the Flory–Huggins lattice model and the Boltzmann entropy theorem, we develop a reticular free energy suit for the network structure of MMC hydrogels. Taking the random factor into account, with the timedependent GinzburgLandau (TDGL) mesoscopic simulation method, we set up a stochastic Cahn–Hilliard equation, designated herein as the MMCTDGL equation. The stochastic term in the equation is constructed appropriately to satisfy the fluctuationdissipation theorem and is discretized on a spatial grid for the simulation. A semiimplicitmore » 
A Fourier Spectral Moving Mesh Method for the CahnHilliard Equation with Elasticity
In recent years, Fourier spectral methods have emerged as competitive numerical methods for largescale phase field simulations of microstructures in computational materials sciences. To further improve their effectiveness, we recently developed a new adaptive Fourierspectral semiimplicit method (AFSIM) for solving the phase field equation by combining an adaptive moving mesh method and the semiimplicit Fourier spectral algorithm. In this paper, we present the application of AFSIM to the CahnHilliard equation with inhomogeneous, anisotropic elasticity. Numerical implementations and test examples in both two and three dimensions are considered with a particular illustration using the wellstudied example of misfitting particles in amore » 
Supersymmetry in the CahnHilliardCook theory
We demonstrate that the early stage of spinodal decomposition and continuous ordering is supersymmetric in systems with longrange interactions. Indeed, the CahnHilliardCook linear theory can be described equivalently as an evolution of the order parameter (bosons) or objects that obey an exclusion principle (fermions). Coupled with a cluster representation of the CahnHilliardCook theory, supersymmetry suggests a new physical picture of the early stage of spinodal decomposition and continuous ordering that provides a greater understanding of earlystage morphology, and could provide a basis for an improved description of the nonlinear regime. 
Existence of equilibria for the CahnHilliard equation via local minimizers of the perimeter
In this paper we shall be interested in establishing existence of certain equilibria of the CahnHiliard equation. 20 refs.