22 research outputs found
Positivity-preserving, energy stable numerical schemes for the Cahn-Hilliard equation with logarithmic potential
We present and analyze finite difference numerical schemes for the Allen
Cahn/Cahn-Hilliard equation with a logarithmic Flory Huggins energy potential.
Both the first order and second order accurate temporal algorithms are
considered. In the first order scheme, we treat the nonlinear logarithmic terms
and the surface diffusion term implicitly, and update the linear expansive term
and the mobility explicitly. We provide a theoretical justification that, this
numerical algorithm has a unique solution such that the positivity is always
preserved for the logarithmic arguments. In particular, our analysis reveals a
subtle fact: the singular nature of the logarithmic term around the values of
and 1 prevents the numerical solution reaching these singular values, so
that the numerical scheme is always well-defined as long as the numerical
solution stays similarly bounded at the previous time step. Furthermore, an
unconditional energy stability of the numerical scheme is derived, without any
restriction for the time step size. The unique solvability and the
positivity-preserving property for the second order scheme are proved using
similar ideas, in which the singular nature of the logarithmic term plays an
essential role. For both the first and second order accurate schemes, we are
able to derive an optimal rate convergence analysis, which gives the full order
error estimate. The case with a non-constant mobility is analyzed as well. We
also describe a practical and efficient multigrid solver for the proposed
numerical schemes, and present some numerical results, which demonstrate the
robustness of the numerical schemes
Energy Stable Numerical Schemes for Ternary Cahn-Hilliard System
We present and analyze a uniquely solvable and unconditionally energy stable numerical scheme for the ternary Cahn-Hilliard system, with a polynomial pattern nonlinear free energy expansion. One key difficulty is associated with presence of the three mass components, though a total mass constraint reduces this to two components. Another numerical challenge is to ensure the energy stability for the nonlinear energy functional in the mixed product form, which turns out to be non-convex, non-concave in the three-phase space. to overcome this subtle difficulty, we add a few auxiliary terms to make the combined energy functional convex in the three-phase space, and this, in turn, yields a convex-concave decomposition of the physical energy in the ternary system. Consequently, both the unique solvability and the unconditional energy stability of the proposed numerical scheme are established at a theoretical level. in addition, an optimal rate convergence analysis in the ℓ∞(0,T;HN-1)∩ℓ2(0,T;HN1) norm is provided, with Fourier pseudo-spectral discretization in space, which is the first such result in this field. to deal with the nonlinear implicit equations at each time step, we apply an efficient preconditioned steepest descent (PSD) algorithm. a second order accurate, modified BDF scheme is also discussed. a few numerical results are presented, which confirm the stability and accuracy of the proposed numerical scheme