CONVERGENCE ANALYSIS OF A FULLY DISCRETE FINITE DIFFERENCE SCHEME FOR THE CAHN-HILLIARD-HELE-SHAW EQUATION

被引:118
作者
Chen, Wenbin [1 ]
Liu, Yuan [2 ]
Wang, Cheng [3 ]
Wise, Steven M. [4 ]
机构
[1] Fudan Univ, Sch Math Sci, Shanghai Key Lab Contemporary Appl Math, Shanghai 200433, Peoples R China
[2] Fudan Univ, Sch Math Sci, Shanghai 200433, Peoples R China
[3] Univ Massachusetts, Dept Math, N Dartmouth, MA 02747 USA
[4] Univ Tennessee, Dept Math, Knoxville, TN 37996 USA
基金
美国国家科学基金会;
关键词
Cahn-Hilliard-Hele-Shaw; Darcy's law; convex splitting; finite difference method; unconditional energy stability; discrete Gagliardo-Nirenberg inequality; discrete Gronwall inequality; DIFFUSE INTERFACE MODEL; ENERGY STABLE SCHEME; THIN-FILM MODEL; ELEMENT-METHOD; ALLEN-CAHN; APPROXIMATIONS; FLUID; RECONNECTION; SIMULATION; STABILITY;
D O I
10.1090/mcom3052
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
We present an error analysis for an unconditionally energy stable, fully discrete finite difference scheme for the Cahn-Hilliard-Hele-Shaw equation, a modified Cahn-Hilliard equation coupled with the Darcy flow law. The scheme, proposed by S. M. Wise, is based on the idea of convex splitting. In this paper, we rigorously prove first order convergence in time and second order convergence in space. Instead of the (discrete) L-s(infinity) (0, T; L-h(2))boolean AND L-s(2) (0, T; H-h(2)) error estimate, which would represent the typical approach, we provide a discrete L-s(infinity) (0, T; H-h(1))boolean AND L-s(2) (0, T; H-h(3)) error estimate for the phase variable, which allows us to treat the nonlinear convection term in a straightforward way. Our convergence is unconditional in the sense that the time step s is in no way constrained by the mesh spacing h. This is accomplished with the help of an L-s(2) (0, T; H-h(3)) bound of the numerical approximation of the phase variable. To facilitate both the stability and convergence analyses, we establish a finite difference analog of a Gagliardo-Nirenberg type inequality.
引用
收藏
页码:2231 / 2257
页数:27
相关论文
共 48 条
[21]   ERROR-ESTIMATES WITH SMOOTH AND NONSMOOTH DATA FOR A FINITE-ELEMENT METHOD FOR THE CAHN-HILLIARD EQUATION [J].
ELLIOTT, CM ;
LARSSON, S .
MATHEMATICS OF COMPUTATION, 1992, 58 (198) :603-630
[22]   A NONCONFORMING FINITE-ELEMENT METHOD FOR THE TWO-DIMENSIONAL CAHN-HILLIARD EQUATION [J].
ELLIOTT, CM ;
FRENCH, DA .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 1989, 26 (04) :884-903
[23]   Unconditionally gradient stable time marching the Cahn-Hilliard equation [J].
Eyre, DJ .
COMPUTATIONAL AND MATHEMATICAL MODELS OF MICROSTRUCTURAL EVOLUTION, 1998, 529 :39-46
[24]   Error analysis of a mixed finite element method for the Cahn-Hilliard equation [J].
Feng, XB ;
Prohl, A .
NUMERISCHE MATHEMATIK, 2004, 99 (01) :47-84
[25]   Fully discrete dynamic mesh discontinuous Galerkin methods for the Cahn-Hilliard equation of phase transition [J].
Feng, Xiaobing ;
Karakashian, Ohannes A. .
MATHEMATICS OF COMPUTATION, 2007, 76 (259) :1093-1117
[26]   Analysis of finite element approximations of a phase field model for two-phase fluids [J].
Feng, Xiaobing ;
He, Yinnian ;
Liu, Chun .
MATHEMATICS OF COMPUTATION, 2007, 76 (258) :539-571
[27]   Fully discrete finite element approximations of the Navier-Stokes-Cahn-Hilliard diffuse interface model for two-phase fluid flows [J].
Feng, Xiaobing .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2006, 44 (03) :1049-1072
[28]   ANALYSIS OF A DARCY-CAHN-HILLIARD DIFFUSE INTERFACE MODEL FOR THE HELE-SHAW FLOW AND ITS FULLY DISCRETE FINITE ELEMENT APPROXIMATION [J].
Feng, Xiaobing ;
Wise, Steven .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2012, 50 (03) :1320-1343
[29]  
Feng XB, 2008, J COMPUT MATH, V26, P767
[30]   LONG TIME STABILITY OF A CLASSICAL EFFICIENT SCHEME FOR TWO-DIMENSIONAL NAVIER-STOKES EQUATIONS [J].
Gottlieb, S. ;
Tone, F. ;
Wang, C. ;
Wang, X. ;
Wirosoetisno, D. .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2012, 50 (01) :126-150