In this paper, we construct and analyze a uniquely solvable, positivity preserving and unconditionally energy stable finite-difference scheme for the periodic three-component Macromolecular Microsphere Composite (MMC) hydrogels system, a ternary Cahn
-Hilliard system with a Flory-Huggins-deGennes free energy potential. The proposed scheme is based on a convex-concave decomposition of the given energy functional with two variables, and the centered difference method is adopted in space. We provide a theoretical justification that this numerical scheme has a pair of unique solutions, such that the positivity is always preserved for all the singular terms, i.e., not only two phase variables are always between $0$ and $1$, but also the sum of two phase variables is between $0$ and $1$, at a point-wise level. In addition, we use the local Newton approximation and multigrid method to solve this nonlinear numerical scheme, and various numerical results are presented, including the numerical convergence test, positivity-preserving property test, energy dissipation and mass conservation properties.
In this paper, two finite difference numerical schemes are proposed and analyzed for the droplet liquid film model, with a singular Leonard-Jones energy potential involved. Both first and second order accurate temporal algorithms are considered. In t
he first order scheme, the convex potential and the surface diffusion terms are implicitly, while the concave potential term is updated explicitly. Furthermore, we provide a theoretical justification that this numerical algorithm has a unique solution, such that the positivity is always preserved for the phase variable at a point-wise level, so that a singularity is avoided in the scheme. In fact, the singular nature of the Leonard-Jones potential term around the value of 0 prevents the numerical solution reaching such singular value, so that the positivity structure is always preserved. Moreover, an unconditional energy stability of the numerical scheme is derived, without any restriction for the time step size. In the second order numerical scheme, the BDF temporal stencil is applied, and an alternate convex-concave decomposition is derived, so that the concave part corresponds to a quadratic energy. In turn, the combined Leonard-Jones potential term is treated implicitly, and the concave part the is approximated by a second order Adams-Bashforth explicit extrapolation, and an artificial Douglas-Dupont regularization term is added to ensure the energy stability. The unique solvability and the positivity-preserving property for the second order scheme could be similarly established. In addition, optimal rate convergence analysis is provided for both the first and second order accurate schemes. A few numerical simulation results are also presented, which demonstrate the robustness of the numerical schemes.
We develop a theoretical foundation for the application of Nesterovs accelerated gradient descent method (AGD) to the approximation of solutions of a wide class of partial differential equations (PDEs). This is achieved by proving the existence of an
invariant set and exponential convergence rates when its preconditioned version (PAGD) is applied to minimize locally Lipschitz smooth, strongly convex objective functionals. We introduce a second-order ordinary differential equation (ODE) with a preconditioner built-in and show that PAGD is an explicit time-discretization of this ODE, which requires a natural time step restriction for energy stability. At the continuous time level, we show an exponential convergence of the ODE solution to its steady state using a simple energy argument. At the discrete level, assuming the aforementioned step size restriction, the existence of an invariant set is proved and a matching exponential rate of convergence of the PAGD scheme is derived by mimicking the energy argument and the convergence at the continuous level. Applications of the PAGD method to numerical PDEs are demonstrated with certain nonlinear elliptic PDEs using pseudo-spectral methods for spatial discretization, and several numerical experiments are conducted. The results confirm the global geometric and mesh size-independent convergence of the PAGD method, with an accelerated rate that is improved over the preconditioned gradient descent (PGD) method.
We discuss two doubly degenerate Cahn-Hilliard (DDCH) models for isotropic surface diffusion. Degeneracy is introduced in both the mobility function and a restriction function associated to the chemical potential. Our computational results suggest th
at the restriction functions yield more accurate approximations of surface diffusion. We consider a slight generalization of a model that has appeared before, which is non-variational, meaning there is no clear energy that is dissipated along the solution trajectories. We also introduce a new variational and, more precisely, energy dissipative model, which can be related to the generalized non-variational model. For both models we use formal matched asymptotics to show the convergence to the sharp interface limit of surface diffusion.
In this paper, we study a novel second-order energy stable Backward Differentiation Formula (BDF) finite difference scheme for the epitaxial thin film equation with slope selection (SS). One major challenge for the higher oder in time temporal discre
tization is how to ensure an unconditional energy stability and an efficient numerical implementation. We propose a general framework for designing the higher order in time numerical scheme with unconditional energy stability by using the BDF method with constant coefficient stabilized terms. Based on the unconditional energy stability property, we derive an $L^infty_h (0,T; H_{h}^2)$ stability for the numerical solution and provide an optimal the convergence analysis. To deal with the 4-Laplacian solver in an $L^{2}$ gradient flow at each time step, we apply an efficient preconditioned steepest descent algorithm and preconditioned nonlinear conjugate gradient algorithm to solve the corresponding nonlinear system. Various numerical simulations are present to demonstrate the stability and efficiency of the proposed schemes and slovers.
