7 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 regularized models for logarithmic SPDEs and their numerical approximations
Understanding the properties of the stochastic phase field models is crucial
to model processes in several practical applications, such as soft matters and
phase separation in random environments. To describe such random evolution,
this work proposes and studies two mathematical models and their numerical
approximations for parabolic stochastic partial differential equation (SPDE)
with a logarithmic Flory--Huggins energy potential. These multiscale models are
built based on a regularized energy technique and thus avoid possible
singularities of coefficients. According to the large deviation principle, we
show that the limit of the proposed models with small noise naturally recovers
the classical dynamics in deterministic case. Moreover, when the driving noise
is multiplicative, the Stampacchia maximum principle holds which indicates the
robustness of the proposed model. One of the main advantages of the proposed
models is that they can admit the energy evolution law and asymptotically
preserve the Stampacchia maximum bound of the original problem. To numerically
capture these asymptotic behaviors, we investigate the semi-implicit
discretizations for regularized logrithmic SPDEs. Several numerical results are
presented to verify our theoretical findings.Comment: 26 pages, 5 figure
A Positivity Preserving, Energy Stable Finite Difference Scheme for the Flory-Huggins-Cahn-Hilliard-Navier-Stokes System
In this paper, we propose and analyze a finite difference numerical scheme for the Cahn-Hilliard-Navier-Stokes system, with logarithmic Flory-Huggins energy potential. in the numerical approximation to the singular chemical potential, the logarithmic term and the surface diffusion term are implicitly updated, while an explicit computation is applied to the concave expansive term. Moreover, the convective term in the phase field evolutionary equation is approximated in a semi-implicit manner. Similarly, the fluid momentum equation is computed by a semi-implicit algorithm: implicit treatment for the kinematic diffusion term, explicit update for the pressure gradient, combined with semi-implicit approximations to the fluid convection and the phase field coupled term, respectively. Such a semi-implicit method gives an intermediate velocity field. Subsequently, a Helmholtz projection into the divergence-free vector field yields the velocity vector and the pressure variable at the next time step. This approach decouples the Stokes solver, which in turn drastically improves the numerical efficiency. the positivity-preserving property and the unique solvability of the proposed numerical scheme is theoretically justified, i.e., the phase variable is always between -1 and 1, following the singular nature of the logarithmic term as the phase variable approaches the singular limit values. in addition, an iteration construction technique is applied in the positivity-preserving and unique solvability analysis, motivated by the non-symmetric nature of the fluid convection term. the energy stability of the proposed numerical scheme could be derived by a careful estimate. a few numerical results are presented to validate the robustness of the proposed numerical scheme
A Modified Crank-Nicolson Numerical Scheme for the Flory-Huggins Cahn-HilliardModel
In this paper we propose and analyze a second order accurate numerical scheme for the Cahn-Hilliard equation with logarithmic Flory Huggins energy potential. a modified Crank-Nicolson approximation is applied to the logarithmic nonlinear term, while the expansive term is updated by an explicit second order Adams- Bashforth extrapolation, and an alternate temporal stencil is used for the surface diffusion term. a nonlinear artificial regularization term is added in the numerical scheme, which ensures the positivity-preserving property, i.e., the numerical value of the phase variable is always between -1 and 1 at a point-wise level. Furthermore, an unconditional energy stability of the numerical scheme is derived, leveraging the special form of the logarithmic approximation term. in addition, an optimal rate convergence estimate is provided for the proposed numerical scheme, with the help of linearized stability analysis. a few numerical results, including both the constant-mobility and solution-dependent mobility flows, are presented to validate the robustness of the proposed numerical scheme