期刊
COMMUNICATIONS IN COMPUTATIONAL PHYSICS
卷 23, 期 2, 页码 572-602出版社
GLOBAL SCIENCE PRESS
DOI: 10.4208/cicp.OA-2016-0197
关键词
Cahn-Hilliard equation; energy stable BDF; Douglas-Dupont regularization; mixed finite element; energy stability
资金
- NSFC [11671098, 11331004, 91630309]
- 111 project [B08018, NSFDMS-1418689]
- NSF [DMS-1418692]
- University of California-San Diego
- Key Laboratory of Mathematics for Nonlinear Sciences, Fudan University
- Direct For Mathematical & Physical Scien
- Division Of Mathematical Sciences [1418689] Funding Source: National Science Foundation
In this paper we present a second order accurate (in time) energy stable numerical scheme for the Cahn-Hilliard (CH) equation, with a mixed finite element approximation in space. Instead of the standard second order Crank-Nicolson methodology, we apply the implicit backward differentiation formula (BDF) concept to derive second order temporal accuracy, but modified so that the concave diffusion term is treated explicitly. This explicit treatment for the concave part of the chemical potential ensures the unique solvability of the scheme without sacrificing energy stability. An additional term A tau Delta(u(k+1) - u(k)) is added, which represents a second order Douglas-Dupont-type regularization, and a careful calculation shows that energy stability is guaranteed, provided the mild condition A >= 1/16 is enforced. In turn, a uniform in time H-1 bound of the numerical solution becomes available. As a result, we are able to establish an l(infinity)(0, T; L-2) convergence analysis for the proposed fully discrete scheme, with full O(tau(2) + h(2)) accuracy. This convergence turns out to be unconditional; no scaling law is needed between the time step size tau and the spatial grid size h. A few numerical experiments are presented to conclude the article.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据