Analysis of Mixed Interior Penalty Discontinuous Galerkin Methods for the Cahn-Hilliard Equation and the Hele-Shaw Flow
نویسندگان
چکیده
This paper proposes and analyzes two fully discrete mixed interior penalty discontinuous Galerkin (DG) methods for the fourth order nonlinear Cahn-Hilliard equation. Both methods use the backward Euler method for time discretization and interior penalty discontinuous Galerkin methods for spatial discretization. They differ from each other on how the nonlinear term is treated, one of them is based on fully implicit time-stepping and the other uses the energy-splitting timestepping. The primary goal of the paper is to prove the convergence of the numerical interfaces of the DG methods to the interface of the Hele-Shaw flow. This is achieved by establishing error estimates that depend on −1 only in some low polynomial orders, instead of exponential orders. Similar to [15], the crux is to prove a discrete spectrum estimate in the discontinuous Galerkin finite element space. However, the validity of such a result is not obvious because the DG space is not a subspace of the (energy) space H1(Ω) and it is larger than the finite element space. This difficult is overcome by a delicate perturbation argument which relies on the discrete spectrum estimate in the finite element space proved in [15]. Numerical experiment results are also presented to gauge the theoretical results and the performance of the proposed fully discrete mixed DG methods.
منابع مشابه
Numerical Analysis and Scientific Computing Preprint Seria C0 Interior Penalty Discontinuous Galerkin approximation of a sixth order Cahn-Hilliard equation modeling microemulsification processes
Microemulsions can be modeled by an initial-boundary value problem for a sixth order Cahn-Hilliard equation. Introducing the chemical potential as a dual variable, a Ciarlet-Raviart type mixed formulation yields a system consisting of a linear second order evolutionary equation and a nonlinear fourth order equation. The spatial discretization is done by a C0 Interior Penalty Discontinuous Galer...
متن کاملA Posteriori Error Estimates for Finite Element Approximations of the Cahn-hilliard Equation and the Hele-shaw Flow
This paper develops a posteriori error estimates of residual type for conforming and mixed finite element approximations of the fourth order Cahn-Hilliard equation ut + ∆ ` ε∆u− ε−1f(u) ́ = 0. It is shown that the a posteriori error bounds depends on ε−1 only in some low polynomial order, instead of exponential order. Using these a posteriori error estimates, we construct an adaptive algorithm f...
متن کاملAn efficient fully-discrete local discontinuous Galerkin method for the Cahn-Hilliard-Hele-Shaw system
متن کامل
A Hybridized Crouziex-Raviart Nonconforming Finite Element and Discontinuous Galerkin Method for a Two-Phase Flow in the Porous Media
In this study, we present a numerical solution for the two-phase incompressible flow in the porous media under isothermal condition using a hybrid of the linear lower-order nonconforming finite element and the interior penalty discontinuous Galerkin (DG) method. This hybridization is developed for the first time in the two-phase modeling and considered as the main novelty of this research.The p...
متن کاملA discontinuous Galerkin method for the Cahn-Hilliard equation
A discontinuous Galerkin finite element method has been developed to treat the high-order spatial derivatives appearing in the Cahn–Hilliard equation. The Cahn–Hilliard equation is a fourth-order nonlinear parabolic partial differential equation, originally proposed to model phase segregation of binary alloys. The developed discontinuous Galerkin approach avoids the need for mixed finite elemen...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- SIAM J. Numerical Analysis
دوره 54 شماره
صفحات -
تاریخ انتشار 2016