UNIFORM STABILITY AND ERROR ANALYSIS FOR SOME DISCONTINUOUS GALERKIN METHODS

In this paper,we provide a number of new estimates on the stability and convergence of both hybrid discontinuous Galerkin(HDG)and weak Galerkin(WG)methods.By using the standard Brezzi theory on mixed methods,we carefully define appropriate norms for the various discretization variables and then establish that the stability and error estimates hold uniformly with respect to stabilization and discretization parameters.As a result,by taking appropriate limit of the stabilization parameters,we show that the HDG method converges to a primal conforming method and the WG method converges to a mixed conforming method.     

A weak Galerkin mixed finite element method for the Helmholtz equation with large wave numbers

In this article, a new weak Galerkin mixed finite element method is introduced and analyzed for the Helmholtz equation with large wave numbers. The stability and well‐posedness of the method are established for any wave number k without mesh size constraint. Allowing the use of discontinuous approximating functions makes weak Galerkin mixed method highly flexible in term of little restrictions on approximations and meshes. In the weak Galerkin mixed finite element formulation, approximation functions can be piecewise polynomials with different degrees on different elements and meshes can consist elements with different shapes. Suboptimal order error estimates in both discrete H¹ and L² norms are established for the weak Galerkin mixed finite element solutions. Numerical examples are tested to support the theory.     

A unified study of continuous and discontinuous Galerkin methods

A unified study is presented in this paper for the design and analysis of different finite element methods(FEMs), including conforming and nonconforming FEMs, mixed FEMs, hybrid FEMs, discontinuous Galerkin(DG) methods, hybrid discontinuous Galerkin(HDG) methods and weak Galerkin(WG) methods.Both HDG and WG are shown to admit inf-sup conditions that hold uniformly with respect to both mesh and penalization parameters. In addition, by taking the limit of the stabilization parameters, a WG method is shown to converge to a mixed method whereas an HDG method is shown to converge to a primal method. Furthermore,a special class of DG methods, known as the mixed DG methods, is presented to fill a gap revealed in the unified framework.     

Weak Galerkin finite element method for a class of quasilinear elliptic problems

The purpose of this paper is to study the weak Galerkin finite element method for a class of quasilinear elliptic problems. The weak Galerkin finite element scheme is proved to have a unique solution with the assumption that guarantees the corresponding operator to be strongly monotone and Lipschitz-continuous. An optimal error estimate in a mesh-dependent energy norm is established. Some numerical results are presented to confirm the theoretical analysis.     

Weak Galerkin finite element methods for Parabolic equations

A newly developed weak Galerkin method is proposed to solve parabolic equations. This method allows the usage of totally discontinuous functions in approximation space and preserves the energy conservation law. Both continuous and discontinuous time weak Galerkin finite element schemes are developed and analyzed. Optimal‐order error estimates in both H¹ and L² norms are established. Numerical tests are performed and reported. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A Modified Weak Galerkin Finite Element Method for Sobolev Equation

For Sobolev equation, we present a new numerical scheme based on a modified weak Galerkin finite element method, in which differential operators are approximated by weak forms through the usual integration by parts. In particular, the numerical method allows the use of discontinuous finite element functions and arbitrary shape of element. Optimal order error estimates in discrete $H^1$ and $L^2$ norms are established for the corresponding modified weak Galerkin finite element solutions. Finally, some numerical results are given to verify theoretical results.     

A weak Galerkin method for diffraction gratings

This paper concerns a weak Galerkin method (WGM) for the diffraction of a time-harmonic incident wave impinging upon a one-dimensional periodic grating structure. The existence and uniqueness of the weak Galerkin solution to the grating problem are established using a variational approach. The convergence rate of the proposed WGM is systematically analyzed. Numerical simulations are presented to verify the efficiency of the WGM for solving grating problems.     

A POSTERIORI ERROR ESTIMATES FOR A MODIFIED WEAK GALERKIN FINITE ELEMENT APPROXIMATION OF SECOND ORDER ELLIPTIC PROBLEMS WITH DG NORM

In this paper,we derive a residual based a posteriori error estimator for a modified weak Galerkin formulation of second order elliptic problems.We prove that the error estimator used for interior penalty discontinuous Galerkin methods still gives both upper and lower bounds for the modified weak Galerkin method,though they have essentially different bilinear forms.More precisely,we prove its reliability and efficiency for the actual error measured in the standard DG norm.We further provide an improved a priori error estimate under minimal regularity assumptions on the exact solution.Numerical results are presented to verify the theoretical analysis.     

Superconvergence of a combined mixed finite element and discontinuous Galerkin approximation for an incompressible miscible displacement problem

A combined mixed finite element and discontinuous Galerkin approximation for an incompressible miscible displacement problem which includes molecular diffusion and dispersion in porous media is studied. That is to say, the mixed finite element method is applied to the flow equation, and the transport equation is solved by an interior penalty discontinuous Galerkin method. Convolution of the Darcy velocity approximation with the Bramble-Schatz kernel function and averaging are applied in the evaluation of the coefficients in the Galerkin procedure for the concentration. A superconvergence estimate is obtained. Numerical experimental results are presented to verify the theoretical analysis.     

Low Regularity Error Analysis for Weak Galerkin Finite Element Methods for Second Order Elliptic Problems

This paper presents error estimates in both an energy norm and the $L^2$-norm for the weak Galerkin (WG) finite element methods for elliptic problems with low regularity solutions. The error analysis for the continuous Galerkin finite element remains same regardless of regularity. A totally different analysis is needed for discontinuous finite element methods if the elliptic regularity is lower than H-1.5. Numerical results confirm the theoretical analysis.     

求解Brinkman方程的修正弱有限元方法

本文针对Brinkman方程引入了一种修正弱Galerkin(MWG)有限元方法.我们通过具有两个离散弱梯度算子的变分形式来逼近模型. 在MWG方法中, 分别用次数为$k$和$k-1$的不连续分段多项式来近似速度函数$u$和压力函数$p$. MWG方法的主要思想是用内部函数的平均值代替边界函数. 因此, 与WG方法相比, MWG方法在不降低准确性的同时, 具有更少的自由度, 对于任意次数不超过$k-1$ 的多项式,MWG方法均可以满足稳定性条件. MWG 方法具有高度的灵活性, 它允许在具有一定形状正则性的任意多边形或多面体上使用不连续函数. 针对$H^1$和$L^22$范数下的速度和压力近似解, 建立了最优阶误差估计. 数值算例表明了该方法的准确性, 收敛性和稳定性.     

Superconvergence of a combined mixed finite element and discontinuous Galerkin method for a compressible miscible displacement problem

A combined mixed finite element and discontinuous Galerkin method for a compressible miscible displacement problem which includes molecular diffusion and dispersion in porous media is investigated. That is to say, the mixed finite element method with Raviart-Thomas space is applied to the flow equation, and the transport one is solved by the symmetric interior penalty discontinuous Galerkin (SIPG) approximation. Based on projection interpolations and induction hypotheses, a superconvergence estimate is obtained. During the analysis, an extension of the Darcy velocity along the Gauss line is also used in the evaluation of the coefficients in the Galerkin procedure for the concentration.     

A Fourier spectral-discontinuous Galerkin method for time-dependent 3-D Schrödinger–Poisson equations with discontinuous potentials

In this paper, we propose a high order Fourier spectral-discontinuous Galerkin method for time-dependent Schrödinger–Poisson equations in 3-D spaces. The Fourier spectral Galerkin method is used for the two periodic transverse directions and a high order discontinuous Galerkin method for the longitudinal propagation direction. Such a combination results in a diagonal form for the differential operators along the transverse directions and a flexible method to handle the discontinuous potentials present in quantum heterojunction and supperlattice structures. As the derivative matrices are required for various time integration schemes such as the exponential time differencing and Crank Nicholson methods, explicit derivative matrices of the discontinuous Galerkin method of various orders are derived. Numerical results, using the proposed method with various time integration schemes, are provided to validate the method.     

A convergent FEM-DG method for the compressible Navier–Stokes equations

This paper presents a new numerical method for the compressible Navier–Stokes equations governing the flow of an ideal isentropic gas. To approximate the continuity equation, the method utilizes a discontinuous Galerkin discretization on piecewise constants and a basic upwind flux. For the momentum equation, the method is a new combined discontinuous Galerkin and finite element method approximating the velocity in the Crouzeix–Raviart finite element space. While the diffusion operator is discretized in a standard fashion, the convection and time-derivative are discretized using discontinuous Galerkin on the element average velocity and a Lax–Friedrich type flux. Our main result is convergence of the method to a global weak solution as discretization parameters go to zero. The convergence analysis constitutes a numerical version of the existence analysis of Lions and Feireisl.     

Weak Galerkin finite element methods combined with Crank-Nicolson scheme for parabolic interface problems

This article is devoted to the a priori error estimates of the fully discrete Crank-Nicolson approximation for the linear parabolic interface problem via weak Galerkin finite element methods (WG-FEM). All the finite element functions are discontinuous for which the usual gradient operator is implemented as distributions in properly defined spaces. Optimal order error estimates in both $L^{\infty}(H^1)$ and $L^{\infty}(L^2)$ norms are established for lowest order WG finite element space $({\cal P}_{k}(K),\;{\cal P}_{k-1}(\partial K),\;\big[{\cal P}_{k-1}(K)\big]^2)$. Finally, we give numerical examples to verify the theoretical results.     

Fully discrete dynamic mesh discontinuous Galerkin methods for the Cahn-Hilliard equation of phase transition

Fully discrete discontinuous Galerkin methods with variable mesh- es in time are developed for the fourth order Cahn-Hilliard equation arising from phase transition in materials science. The methods are formulated and analyzed in both two and three dimensions, and are proved to give optimal order error bounds. This coupled with the flexibility of the methods demonstrates that the proposed discontinuous Galerkin methods indeed provide an efficient and viable alternative to the mixed finite element methods and nonconforming (plate) finite element methods for solving fourth order partial differential equations.

     

Strong convergence of the discontinuous Galerkin scheme for the low regularity miscible displacement equations

Strong convergence of the numerical solution to a weak solution is proved for a nonlinear coupled flow and transport problem arising in porous media. The method combines a mixed finite element method for the pressure and velocity with an interior penalty discontinuous Galerkin method in space for the concentration. Using functional tools specific to broken Sobolev spaces, the convergence of the broken gradient of the numerical concentration to the weak solution is obtained in the L² norm. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 489–513, 2017     

弱有限元方法简论

本文简述弱有限元方法(weak Galerkin finite element met,hods)的数学基本原理和计算机实现.弱有限元方法对间断函数引入广义弱微分,并将其应用于偏微分方程相应的变分形式进行数值求解,而数值解的弱连续性则通过稳定子或光滑子来实现.弱有限元方法针对广义函数而构建,是经典有限元方法的一种自然拓广,且能够弥补经典有限元方法的某些缺憾,也因此在科学与工程计算领域具有广泛的应用前景.     

Space-Time adaptive algorithm for the mixed parabolic problem

In this paper we present an a-posteriori error estimator for the mixed formulation of a linear parabolic problem, used for designing an efficient adaptive algorithm. Our space-time discretization consists of lowest order Raviart-Thomas finite element over graded meshes and discontinuous Galerkin method with variable time step. Finally, several examples show that the proposed method is efficient and reliable.     

Mixed Discontinuous Galerkin Time-Stepping Method for Linear Parabolic Optimal Control Problems

In this paper, we discuss the mixed discontinuous Galerkin (DG) finite element approximation to linear parabolic optimal control problems. For the state variables and the co-state variables, the discontinuous finite element method is used for the time discretization and the Raviart-Thomas mixed finite element method is used for the space discretization. We do not discretize the space of admissible control but implicitly utilize the relation between co-state and control for the discretization of the control. We derive a priori error estimates for the lowest order mixed DG finite element approximation. Moveover, for the element of arbitrary order in space and time, we derive a posteriori $L^2(0, T ;L^2(Ω))$ error estimates for the scalar functions, assuming that only the underlying mesh is static. Finally, we present an example to confirm the theoretical result on a priori error estimates.