Efficient local energy dissipation preserving algorithms for the Cahn–Hilliard equation
Introduction
There has been an increasing emphasis on the free interface problem of multi-phase incompressible flows which is often modeled by the diffuse interface method with a long history [39], [44]. Its main idea is to apply one or more continuous functions to describe the volume fraction in the multi-phase immiscible system and represent the fluid interfaces by a thin and smooth transition layers, see, e.g., [1], [3], [7], [20], [22], [33], [34], [37], [47] and [38], [43] for details. The Cahn–Hilliard equation induced by the variational approach, is a widely used phase-field model.
In this paper, we consider the following Cahn–Hilliard (CH) equation [6] where ϕ is an unknown real-valued phase function, which stands for the relative concentration of one phase, μ is the chemical potential which is the variation of the free energy, κ is a parameter measuring the strength of the conformational entropy. Here the free energy is chosen as where Ω is the domain that the fluid occupies, and is the bulk energy density. For the choice of nonlinear potential , Ginzburg–Landau double-well type potential, and Logarithmic Flory–Huggins potential [41], are two widely used nonlinear bulk potentials in the literature, where γ measures the strength of the repulsive potential. As is well known, the CH equation, with suitable boundary conditions, possesses the mass conservation law (MCL) and the energy dissipation law (EDL)
To date, extensive mathematical studies have been carried out for the CH equation, see, e.g., [7], [8], [9], [10], [19], [21], [31], [32], [37], [40], [42], [47] for details. For the analytical front, the steady state solution, the existence and the uniqueness of the system, we refer to [13], [24] and the references therein. For the numerical front, various numerical methods have been considered in the literature. Elliott et al. [12] discussed the physical background and its numerical solution by the Galerkin finite element method. Later, Feng and Prohl [17] developed the mixed finite element method which was convergent with quasi-optimal order in time and optimal order in space. Sun [40] proposed the finite difference method by the method of reduction of order, which is conditionally stable and not mass conserving. Choo et al. [5] gave a conservative nonlinear finite difference scheme, which was proven to preserve the total mass, yet the energy-based stability was not discussed. Subsequently, Furihata [16] gave a conservative finite difference scheme, which was proven to be unconditionally stable in the sense of energy decay. Chen and Shen [8] developed a semi-implicit Fourier spectral scheme, in which the linear part was treated implicitly and the nonlinear part was calculated explicitly. Later, Zhu et al. [49] proposed the semi-implicit Fourier spectral method for the CH equation with a variable mobility. Shen and Yang [42] obtained energy stability method without any a priori assumption on the numerical solution. Zhang and Qiao [50] discussed a finite difference scheme and proposed an adaptive time-stepping technique to quickly solve the 2D CH equation. He et al. [31] proposed a large time-stepping methods. Guo et al. [27] used the convex splitting skills, whose idea was proposed more than two decades earlier by Elliott et al. and Eyre in Ref. [14], [15]. Li et al. [35], [36] developed the semi-implicit Fourier spectral method, which was proved unconditional energy stability for modified energy functionals by the new stabilization techniques. Recently, the new technique called “Invariant Energy Quadratization” (IEQ) which was successfully applied to different phase-field models by authors [30], [46], [48], was extended to handle the CH equation in Ref. [25], [26]. The efficient iterative solution of linear system was investigated for the coupled Navier–Stokes Cahn–Hilliard system, in which the block-triangular preconditioners were introduced in Ref. [2]. Gomez et al. gave a good reference book for stable time integration of phase-field equations in Ref. [23]. We notice that most existing methods and their stabilities which are discussed in the literature are based on the global energy dissipation property which is defined on the global space region and depends on the suitable boundary conditions. Here in this paper, we explore to construct novel schemes basing on our observation that CH equation has a local energy dissipation law which is valid in any time–space region and independent of the boundary conditions. It is known that boundary condition of the system is very important in many practical problems.
Proposition 1 The system (1.1) possesses the local energy dissipation law (LEDL) which can be derived by multiplying μ and to the first line and the second line of (1.1) respectively, and adding the two resluting formulas together. Compared with (1.4), the LEDL (1.5) has an additional flux term which can supply much more information of the original system. Moreover, the LEDL is independent of boundary conditions. Therefore, it is valid in any time space domain for the CH equation with any boundary condition. Under suitable boundary conditions, the LEDL (1.5) implies the total EDL (1.4). This observation motivates us to construct numerical schemes preserving the LEDL (1.5) for the CH equation.
Actually, how to construct numerical schemes preserving certain local invariant quantities for the continuous dynamical systems is brought into sharp focus in recent years. The numerical methods preserving some structural properties of the system are sometime referred as geometric integrators or structure-preserving numerical methods. So far, on construction of the local structure-preserving schemes for some classical PDEs, many results have been achieved [28], [45] which further show that numerical methods with local property supply richer information of the original system than that of schemes with global property.
In this paper, we aim to construct three novel numerical schemes to preserve discrete analogues of the local energy dissipation law (1.5) for the 2D Cahn–Hilliard equation which are achieved through the concatenation of basic algorithms suggested by Wang et al. [45], including the implicit midpoint method, the leap-frog method and the discrete variational derivative method proposed by Furihata et al. [11], [18] for the nonlinear term, on the associated first-order subsystems. The corresponding discrete local energy dissipation laws are rigorously proved by the Leibnitz rules, with no dependence of the boundary conditions. While imposed of the periodic boundary condition or the homogeneous boundary condition, the proposed schemes thereby possess the global EDL as well as MCL.
The remainder of this paper is arranged as follows. Some operator definitions and theirs properties are given in Section 2. In Section 3, we construct three novel numerical schemes preserving the LEDL for solving the CH equation. When the equation is imposed on appropriate boundary conditions, the schemes will conserve the global MCL and EDL. Numerical experiments are presented in Section 4 and some concluding remarks are given in the final section.
Section snippets
Operator definitions and properties
To simplify the presentation, we introduce some notations firstly. For a positive integer , we define the time-step as . Let be two positive integers. The spatial domain is uniformly partitioned with mesh size and The approximation of the value of the function at the node is denoted by . A grid function is called periodic if
Local energy dissipation preserving algorithms
In this section, we apply the local structure-preserving algorithms proposed in [28], [45] for the CH equation (1.1) to develop a series of local energy dissipation preserving algorithms.
Numerical experiments
In this section, we present some numerical results to illustrate the efficiency and accuracy of the proposed methods. We notice that the proposed methods lead to nonlinear algebraic systems, which can be solved using simple fixed-point iterative methods (see [29]). Due to periodic boundary conditions, the fast Fourier transform can be applied for every iteration step. In the following experiments, we take .
Example 4.1 In order to test the convergence rate, we first consider the system (1.1)
Concluding remarks
In this paper, we observe that the CH equation possesses a local energy dissipation law. Based on this observation, we derive the three LEDL numerical schemes for the 2D CH equation by the concatenating method. In particular, when the nonlinear bulk potential function is chosen as the double-well potential, the method [50] proposed by Zhang and Qiao is a special case of our scheme II, which means the method [50] possesses a local property. Utilizing the discrete Leibnitz rule and the
Acknowledgements
This work is supported by the National Natural Science Foundation of China (under Grant Nos. 11771213 and 41504078), the National Key Research and Development Project of China (under Grant No. 2016YFC0600310), the Natural Science Foundation of Jiangsu Province (Grant No. SBK2018040413), the Major Projects of Natural Sciences of University in Jiangsu Province of China (Grant No. 15KJA110002) and the Priority Academic Program Development of Jiangsu Higher Education Institutions, China.
References (50)
- et al.
Local energy-preserving and momentum-preserving algorithms for coupled nonlinear Schrödinger system
J. Comput. Phys.
(2013) - et al.
Conservative nonlinear difference scheme for the Cahn–Hilliard equation-II
Comput. Math. Appl.
(2000) - et al.
High accuracy solutions to energy gradient flows from material science models
J. Comput. Phys.
(2014) - et al.
Applications of semi-implicit Fourier-spectral method to phase field equations
Comput. Phys. Commun.
(1998) Finite difference schemes for that inherit energy conservation or dissipation property
J. Comput. Phys.
(1999)- et al.
Provably unconditionally stable, second-order time-accurate, mixed variational methods for phase-field models
J. Comput. Phys.
(2011) - et al.
Some new structure-preserving algorithms for general multi-symplectic formulations of Hamiltonian PDEs
J. Comput. Phys.
(2014) - et al.
A conservative Fourier pseudo-spectral method for the nonlinear Schrödinger equation
J. Comput. Phys.
(2017) - et al.
On large time-stepping methods for the Cahn–Hilliard equation
Appl. Numer. Math.
(2007) - et al.
Fast and accurate algorithms for simulating coarsening dynamics of Cahn–Hilliard equations
Comput. Mater. Sci.
(2015)
Liquid-vapor phase transition: thermomechanical theory, entropy stable numerical formulation, and boiling simulations
Comput. Methods Appl. Mech. Eng.
A phase field model for the mixture of two incompressible fluids and its approximation by a Fourier-spectral method
Physica D
Linear, first and second-order, unconditionally energy stable numerical schemes for the phase field model of homopolymer blends
J. Comput. Phys.
A novel linear second order unconditionally energy stable scheme for a hydrodynamic Q-tensor model of liquid crystals
Comput. Methods Appl. Mech. Eng.
Diffuse-interface methods in fluid mechanics
Annu. Rev. Fluid Mech.
Preconditioning of a coupled Cahn–Hilliard Navier–Stokes system
Commun. Comput. Phys.
Free energy of a nonuniform system. II. Thermodynamic basis
J. Chem. Phys.
Free energy of a nonuniform system. I. Interfacial energy free energy
J. Chem. Phys.
Efficient energy stable schemes with spectral discretization in space for anisotropic Cahn–Hilliard systems
Commun. Comput. Phys.
A second-order, weakly energy-stable pseudo-spectral scheme for the Cahn–Hilliard equation and its solution by the homogeneous linear iteration method
J. Sci. Comput.
A general framework for deriving integral preserving numerical methods for PDEs
SIAM J. Sci. Comput.
Numerical studies of the Cahn–Hilliard equation for phase separation
IMA J. Appl. Math.
On the Cahn–Hilliard equation
Arch. Ration. Mech. Anal.
The global dynamics of discrete semilinear parabolic equations
SIAM J. Numer. Anal.
Unconditionally gradient stable time marching the Cahn–Hilliard equation
Comput. Math. Models Microstruct. Evol.
Cited by (8)
Two classes of linearly implicit local energy-preserving approach for general multi-symplectic Hamiltonian PDEs
2020, Journal of Computational PhysicsLinearly implicit local energy-preserving algorithm for a class of multi-symplectic Hamiltonian PDEs
2022, Computational and Applied MathematicsLocal Energy Dissipation Rate Preserving Approximations to Driven Gradient Flows with Applications to Graphene Growth
2022, Journal of Scientific ComputingLocal structure-preserving algorithms for the molecular beam epitaxy model with slope selection
2021, Discrete and Continuous Dynamical Systems - Series BBenchmark Problems for the Numerical Discretization of the Cahn-Hilliard Equation with a Source Term
2021, Discrete Dynamics in Nature and SocietyThe Navier–Stokes–Cahn–Hilliard model with a high-order polynomial free energy
2020, Acta Mechanica