Unconditionally energy stable second-order numerical schemes for the Functionalized Cahn–Hilliard gradient flow equation based on the SAV approach
Introduction
Adding hydrophilic groups to the backbones of hydrophobic polymers, which is called functionalization, can make amphiphilic polymer materials. The Functionalized Cahn–Hilliard (FCH) free energy was proposed in [1] and improved in [2], [3], [4] to describe a binary mixture composed of the amphiphilic polymer and the polar solvent (such as water). In general, the FCH energy is non-increasing when the phase separation occurs in the mixture. Meanwhile, many different kinds of nanostructures gradually emerge, for example micelles, bilayers, pores, and so on. To study these interesting structures, the so-called FCH gradient flow equation was proposed, see [4]. This equation can be obtained by applying a dissipation mechanism on the FCH free energy. The FCH equation is a special case of the FCH gradient flow equation and can be formulated as where are constants, is the phase-field variable and varies between and , and the double-well potential function is a nonlinear function. The FCH equation is actually an gradient flow of the FCH free energy.
Generally speaking, it is almost impossible to obtain an analytic solution to the FCH equation, although a recent work [5] has established the global-in-time Gevrey regularity solutions for it. Hence, the approximate solution obtained by the numerical method becomes the focus. However, it is not easy to solve the FCH equation (1.1) numerically due to its high nonlinearity and high derivative order. In general, the explicit Euler method is one of the simplest methods to obtain the numerical solution. But for the FCH equation, a very small time step is necessary to ensure the numerical stability in calculations when this method is adopted. Hence, the explicit Euler method is not popular although it is very easy to implement. Therefore, designing some highly efficient numerical schemes for the FCH equation or the FCH gradient flow equation has become very important. In recent years, based on the ideas of highly efficient algorithms of Cahn–Hilliard-type equations, see [6], [7], [8], [9], [10], [11], [12], [13], [14], [15], [16], [17], [18], [19], a number of numerical schemes have sprung up for approximating the solution to the FCH equation.
The fully implicit scheme was constructed and used in [20] to solve the FCH equation. Lots of numerical experiments have demonstrated that large time steps indeed can be used in this scheme, see [20], [21], [22]. However, a nonlinear system has to be solved, which may make the computation less efficient if an efficient iterative solver is not available. Hence, in order to improve the computational efficiency, both an adaptive time strategy and a highly efficient preconditioned iterative solver were designed by the authors. Additionally, higher order schemes of the FCH equation were also constructed in [23] by using the fully implicit method. However, the unconditional energy stability and unique solvability of all the above fully implicit schemes cannot be guaranteed, which greatly limits their applications.
Another classic method is Eyre’s convex-splitting method [24], which was introduced by Eyre for solving the Cahn–Hilliard equation. This method is widely used to design unconditionally energy stable schemes for phase field models, see [9], [12], [19]. A first-order scheme was proposed in [23] for the FCH equation based on the convex-splitting approach. This scheme satisfies the unconditional energy stability for Galerkin-type spatial discretization methods. However, its application is limited due to the lack of efficient iterative solvers. Another slightly modified version of this scheme, which we call implicit–explicit (IMEX) scheme in this paper, was also proposed in [23]. Unfortunately, the IMEX scheme is neither unconditionally energy stable nor uniquely solvable, although the decay in energy was verified by many numerical experiments. Based on the theory of convex-splitting, Feng and his co-workers established a novel first-order convex-splitting scheme by adding two auxiliary terms to the FCH energy [25]. The unconditional energy stability and unconditional unique solvability of the scheme were proved at the theoretical level. In addition, a second-order scheme based on the BDF2 formula was also presented in [25], but without any corresponding theoretical results. The main advantage of the convex-splitting method is that one can use it to design unconditionally energy stable schemes as long as a correct or reasonable splitting of the energy can be found. However, it should be noted that a nonlinear equation, in general, will be generated when the convex-splitting method is considered, which may make the calculation slow if an efficient iterative solver is unavailable.
The third approach is the so-called stabilization method that achieves great success in solving the classic Cahn–Hilliard equation [16] or other similar phase field equations [6], [11], [15], [18]. By adding stabilization terms to the semi-implicit scheme (i.e., evaluating the linear terms implicitly and the nonlinear terms explicitly), a linear, first-order, unconditionally energy stable scheme can be constructed. For the FCH equation, Chen and his co-workers put up a first-order scheme by using this method in [26]. Then, Guo and his co-workers adopted it and the local discontinuous Galerkin method to study the FCH equation in [27]. Unfortunately, neither of the two works gave the theoretical analysis of the energy stability (in time) of this time-semi-discrete scheme.
The Scalar Auxiliary Variable (SAV) method is very effective and easy to implement for solving the gradient flow equations [28], [29], [30], [31], [32]. One only needs to solve several linear, constant-coefficient equations in a SAV scheme. In [22], the SAV method was used for solving the FCH equation. However, large time steps are unallowable in the second-order scheme, although it is unconditionally energy stable. To solve this problem, a stabilization term was added to the second-order SAV scheme. It was shown in [22] that the new stabilized second-order scheme was unconditionally energy stable and unconditionally uniquely solvable. This work is remarkable and worthy of recognition. However, an extra linear transformation has to be considered for using this scheme when the potential function in (1.1a) is asymmetric. In addition, due to the introduction of the stabilization term, it is so hard, or almost impossible, to apply variable time steps to it without affecting the unconditional energy stability.
In this paper, we study the FCH gradient flow equation with various dissipation mechanisms, not just the FCH equation (1.1). Note that most of the above works focus on the study of the FCH equation. Based on the SAV approach, we successfully design three new second-order linear schemes for the FCH gradient flow equation (2.5) where an asymmetric double-well potential function is considered. The success is mainly due to the fact that a new auxiliary variable is introduced in a subtle way which is different from the one adopted in [22]. Specially, in order to greatly improve the computational efficiency, we construct a linear, second-order, unconditionally energy stable numerical scheme with variable time steps, which is almost impossible by using the splitting manner in [22]. We use the Fourier pseudo-spectral method [33] to discrete spatial variables in this paper when the periodic boundary conditions are considered. The Fourier pseudo-spectral method has a very low phase error and allows us to use fewer spatial nodes to achieve higher accuracy compared to other methods, for instance the finite element method [34], the finite volume method [35], [36] and the finite difference method [8], [25]. Moreover, by means of the Fast Fourier Transform (FFT) algorithm, the computational efficiency can be greatly improved. In addition, the Legendre spectral method [37] or the Chebyshev spectral method [38] can be used for non-periodic boundary conditions. However, it should be noted that the numerical solutions in this case may be polluted by the rounding-off errors resulting from approximations of high order derivatives explicitly.
The rest of the paper is organized as follows. In Section 2, we introduce the FCH free energy and the FCH gradient flow equation, respectively. In Section 3, we introduce the auxiliary variable so that the standard SAV method can be used to rewrite the original FCH gradient flow equation. In Section 4, we present the constructed second-order SAV schemes and establish the corresponding unconditional energy stability and unique solvability for each scheme. In Section 5, we perform some numerical experiments to verify the theoretical results. Finally, some conclusions follow in Section 6.
Section snippets
FCH gradient flow equation
We firstly introduce some necessary and useful notations which will be used throughout this paper. We denote the inner product and norm associated with the space by and , respectively. And is used to represent the norm of the space .
The FCH free energy [2], [4] with general form can be written as where represents the width of the polymer–solvent interface, is a constant and usually set to be or , and are
The SAV approach for the FCH gradient flow equation
To construct the unconditionally energy stable schemes for Eq. (2.5) based on the SAV approach, the first thing we should do is to introduce an auxiliary variable. The simplest approach is to set [43] with a positive constant such that . Note that the energy in (2.1) is bounded from below, which has been proven theoretically in [22], [25], [44]. The above manner of introducing the auxiliary variable is reasonable and feasible from a mathematical point of view. Then we
Numerical schemes for the FCH gradient flow equation
In this section, we design three second-order accurate (in time) numerical schemes to solve the system (3.27) and present the theoretical analysis results. The unconditional energy stability and unique solvability are our main concerns. In addition, an effective adaptive time strategy is introduced in the end. We emphasize that all the schemes designed below for also apply to other types of the FCH gradient flows.
Numerical experiments
In this section, we carry out some numerical experiments to verify the accuracy and efficiency of the numerical schemes designed in this paper. Additionally, we also perform some numerical comparisons between the BDF2 version of the fully implicit scheme and our SAV schemes. These comparisons mainly focus on the accuracy and computational efficiency.
We use the periodic boundary conditions for all the following computations. The Fourier pseudo-spectral method is employed to discretize the space
Conclusions
In this work we consider the FCH gradient flow equation (2.5) where an asymmetric double-well potential function is considered. The difficulty in solving this equation lies mainly in how to deal with nonlinear terms. With the help of the SAV method, we successfully construct three second-order accurate (in time) schemes for the FCH gradient flow equation. All the schemes are linear and easy to solve. More specifically, only two linear constant-coefficient equations of the same order as the
Acknowledgements
C.Zhang and J.Ouyang are grateful to the support from the National Natural Science Foundation of China (Grant Nos. 11671321, 11971387, 11901051, 91434201). C.Zhang also acknowledges partial support from the Natural Science Foundation of Shaanxi Province, China (Grant No. 2020JQ-338).
References (51)
- et al.
Decoupled energy stable schemes for phase-field vesicle membrane model
J. Comput. Phys.
(2015) - et al.
Efficient energy stable schemes for isotropic and strongly anisotropic Cahn–Hilliard systems with the Willmore regularization
J. Comput. Phys.
(2018) - et al.
An energy stable fourth order finite difference scheme for the Cahn–Hilliard equation
J. Comput. Appl. Math.
(2019) - et al.
Stabilized linear semi-implicit schemes for the nonlocal Cahn–Hilliard equation
J. Comput. Phys.
(2018) - et al.
Analysis of a novel finite element method for a modified Cahn–Hilliard–Hele–Shaw system
J. Comput. Appl. Math.
(2020) - et al.
An efficient two-grid method for the Cahn–Hilliard equation with the concentration-dependent mobility and the logarithmic Flory-Huggins bulk potential
Appl. Math. Comput.
(2020) - et al.
High accuracy solutions to energy gradient flows from material science models
J. Comput. Phys.
(2014) - et al.
Numerical comparison of modified-energy stable SAV-type schemes and classical BDF methods on benchmark problems for the functionalized Cahn–Hilliard equation
J. Comput. Phys.
(2020) - et al.
Efficient spectral-Galerkin methods for systems of coupled second-order equations and their applications
J. Comput. Phys.
(2012) - et al.
A variant of scalar auxiliary variable approaches for gradient flows
J. Comput. Phys.
(2019)
The scalar auxiliary variable (SAV) approach for gradient flows
J. Comput. Phys.
Efficient SAV approach for imaginary time gradient flows with applications to one- and multi-component Bose–Einstein Condensates
J. Comput. Phys.
Appropriate stabilized Galerkin approaches for solving two-dimensional coupled Burgers’ equations at high Reynolds numbers
Comput. Math. Appl.
Numerical simulation of viscoelastic fluid flows using a robust FVM framework on triangular grid
J. Non-Newton. Fluid Mech.
Solving the regularized, strongly anisotropic Cahn–Hilliard equation by an adaptive nonlinear multigrid method
J. Comput. Phys.
Curvature driven flow of bi-layer interfaces
Physica D
Efficient, adaptive energy stable schemes for the incompressible Cahn–Hilliard Navier–Stokes phase-field models
J. Comput. Phys.
An adaptive time-stepping strategy for solving the phase field crystal model
J. Comput. Phys.
Interfacial tension of evaporating emulsion droplets containing amphiphilic block copolymers: effects of solvent and polymer composition
J. Colloid Interface Sci.
Origin of the scattering peak in microemulsions
J. Chem. Phys.
Variational models of network formation and ion transport: applications to perfluorosulfonate ionomer membranes
Polymers
Correlation between structural and interfacial properties of amphiphilic systems
Phys. Rev. Lett.
Pem fuel cells: a mathematical overview
SIAM J. Appl. Math.
Global-in-time gevrey regularity solutions for the functionalized Cahn–Hilliard equation
Discrete Contin. Dyn. Syst. Ser. S
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.
Cited by (10)
SAV Fourier-spectral method for diffuse-interface tumor-growth model
2023, Computers and Mathematics with ApplicationsHigh-order time-accurate, efficient, and structure-preserving numerical methods for the conservative Swift–Hohenberg model
2021, Computers and Mathematics with ApplicationsNumerical analysis of an unconditionally energy-stable reduced-order finite element method for the Allen-Cahn phase field model
2021, Computers and Mathematics with ApplicationsHighly accurate, linear, and unconditionally energy stable large time-stepping schemes for the Functionalized Cahn–Hilliard gradient flow equation
2021, Journal of Computational and Applied MathematicsNumerical analysis of a second-order IPDGFE method for the Allen–Cahn equation and the curvature-driven geometric flow
2021, Computers and Mathematics with Applications