Abstract:In this paper, we provide a theoretical analysis for a preconditioned steepest descent (PSD) iterative solver that improves the computational time of a finite difference numerical scheme for the Cahn-Hilliard equation with Flory-Huggins energy potential. In the numerical design, a convex splitting approach is applied to the chemical potential such that the logarithmic and the surface diffusion terms are treated implicitly while the expansive concave term is treated with an explicit update. The nonlinear and singular nature of the logarithmic energy potential makes the numerical implementation very challenging. However, the positivity-preserving property for the logarithmic arguments, unconditional energy stability, and optimal rate error estimates have been established in a recent work and it has been shown that successful solvers ensure a similar positivity-preserving property at each iteration stage. Therefore, in this work, we will show that the PSD solver ensures a positivity-preserving property at each iteration stage. The PSD solver consists of first computing a search direction (involved with solving a Poisson-like equation) and then takes a one-parameter optimization step over the search direction in which the Newton iteration becomes very powerful. A theoretical analysis is applied to the PSD iteration solver and a geometric convergence rate is proved for the iteration. In particular, the strict separation property of the numerical solution, which indicates a uniform distance between the numerical solution and the singular limit values of $\pm 1$ for the phase variable, plays an essential role in the iteration convergence analysis. A few numerical results are presented to demonstrate the robustness and efficiency of the PSD solver.

Convergence analysis of a block preconditioned steepest descent eigensolver with implicit deflation

Convergence Analysis of Projected Gradient Descent for Schatten-p Nonconvex Matrix Recovery

Convergence Analysis of a Preconditioned Steepest Descent Solver for the Cahn-Hilliard Equation with Logarithmic Potential

A New Projected Variant of the Deflated Block Conjugate Gradient Method

Riemannian Acceleration with Preconditioning for symmetric eigenvalue problems

Geodesic Convexity of the Symmetric Eigenvalue Problem and Convergence of Steepest Descent

A Relaxed Deteriorated Pss Preconditioner for Nonsymmetric Saddle Point Problems from the Steady Navier-Stokes Equation

EPIC: a provable accelerated Eigensolver based on Preconditioning and Implicit Convexity

A projected preconditioned conjugate gradient algorithm for computing many extreme eigenpairs of a Hermitian matrix

A class of new extended generalized shift-splitting preconditioners for block saddle point problems

A Note on Parallel Preconditioning for the All-at-Once Solution of Riesz Fractional Diffusion Equations

On Preconditioned Euler-extrapolated Single-Step Hermitian and Skew-Hermitian Splitting Method for Complex Symmetric Linear Systems

Deflation for the off-diagonal block in symmetric saddle point systems

A Two-Level Block Preconditioned Jacobi–Davidson Method for Multiple and Clustered Eigenvalues of Elliptic Operators

A Refined Block Conjugate Gradient Method for Symmetric Eigenvalue Problems

Optimal Convergence Rate of Self-Consistent Field Iteration for Solving Eigenvector-dependent Nonlinear Eigenvalue Problems

A Note on the Block Alternating Splitting Implicit Iteration Method for Complex Saddle‐point Problems

Global Convergence of Non-Convex Gradient Descent for Computing Matrix Squareroot

Spectral Properties of the Preconditioned Ahss Iteration Method for Generalized Saddle Point Problems

A spectrally preconditioned and initially deflated variant of the restarted block GMRES method for solving multiple right-hand sides linear systems

Convergence analysis of iterative solution with inexact block preconditioning for weak Galerkin finite element approximation of Stokes flow