Continuous-discontinuous finite element method for convection-diffusion problems with characteristic layers

We study convergence properties of a numerical method for convection-diffusion problems with characteristic layers on a layer-adapted mesh. The method couples standard Galerkin with an h-version of the nonsymmetric discontinuous Galerkin finite element method with bilinear elements. In an associated norm, we derive the error estimate as well as the supercloseness result that are uniform in the perturbation parameter. Applying a post-processing operator for the discontinuous Galerkin method, we construct a new numerical solution with enhanced convergence properties.     

分层网格上奇异摄动问题的一致NIPG分析

本文采用非对称内罚间断有限元方法(以下简称NIPG方法)求解一维对流扩散型奇异摄动问题.理论上证明了采用拉格朗日线性元的NIPG方法在分层网格上至多相差一个关于摄动参数对数因子的拟最优阶的一致收敛性,即在能量范数度量下其误差估计为O((log~2(1/e))/N),其中N为网格剖分中单元个数.数值算例验证了理论分析的正确性.     

A TWO-LEVEL FINITE ELEMENT GALERKIN METHOD FOR THE NONSTATIONARY NAVIER-STOKES EQUATIONS II： TIME DISCRETIZATION

In this article we consider the fully discrete two-level finite element Galerkin method for the two-dimensional nonstationary incompressible Navier-Stokes equations. This method consists in dealing with the fully discrete nonlinear Navier-Stokes problem on a coarse mesh with width $H$ and the fully discrete linear generalized Stokes problem on a fine mesh with width $h << H$. Our results show that if we choose $H=O(h^{1/2}$) this method is as the same stability and convergence as the fully discrete standard finite element Galerkin method which needs dealing with the fully discrete nonlinear Navier-Stokes problem on a fine mesh with width $h$. However, our method is cheaper than the standard fully discrete finite element Galerkin method.     

NUMERICAL METHODS FOR MISCIBLE DISPLACEMENT IN POROUS MEDIA INFLUENCED BY MOBILE WATER AND IMMOBILE WATER

1　ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｒｅｃｅｎｔｙｅａｒｓ,ｔｈｅｉｎｔｅｎｔｉｏｎａｌｏｒａｃｃｉｄｅｎｔａｌｒｅｌｅａｓｅｏｆｔｈｅｃｈｅｍｉｃａｌｗａｓｔｅｓｏｎｓｏｉｌｓｈａｓｆｕｒｔｈｅｒｓｔｉｍｕｌａｔｅｄｃｕｒｒｅｎｔｉｎｔｅｒｅｓｔｓｉｎｔｈｅｍｏｖｅｍｅｎｔｏｆｃｈｅｍｉｃａｌｓ.Ｄｉｓｐｌａｃｅｍｅｎｔｓｔｕｄｉｅｓｈａｖｅｂｅｃｏｍｅｉｍｐｏｒｔａｎｔｔｏｏｌｓｉｎｓｏｉｌｐｈｙｓｉｃｓ,ｐａｒｔｉｃｕｌａｒｌｙｆｏｒｐｒｅｄｉｃｔｉｎｇｔｈｅｍｏｖｅｍｅｎｔｏｆｐｅｓｔｃｉｄｅｓ…     

Convergence of Weak Galerkin Finite Element Method for Second Order Linear Wave Equation in Heterogeneous Media

Weak Galerkin finite element method is introduced for solving wave equation with interface on weak Galerkin finite element space $(\mathcal{P}_k(K), \mathcal{P}_{k−1}(∂K), [\mathcal{P}_{k−1}(K)]^2).$ Optimal order a priori error estimates for both space-discrete scheme and implicit fully discrete scheme are derived in $L^∞(L^2)$ norm. This method uses totally discontinuous functions in approximation space and allows the usage of finite element partitions consisting of general polygonal meshes. Finite element algorithm presented here can contribute to a variety of hyperbolic problems where physical domain consists of heterogeneous media.     

输运方程谱间断流线扩散耦合方法

本文讨论了谱有限元方法,构造了求解Ｂｏｌｔｚｍ ａｎｎ 方程球谐函数谱展开和间断流线扩散有限元耦合格式．建立了这种耦合方法的稳定性及最优阶收敛性误差估计．得到了比标准有限元更高的精度．     

A supercloseness result for the discontinuous Galerkin stabilization of convection–diffusion problems on Shishkin meshes

We consider a convection–diffusion problem with Dirichlet boundary conditions posed on a unit square. The problem is discretized using a combination of the standard Galerkin FEM and an h–version of the nonsymmetric discontinuous Galerkin FEM with interior penalties on a layer–adapted mesh with linear/bilinear elements. With specially chosen penalty parameters for edges from the coarse part of the mesh, we prove uniform convergence (in the perturbation parameter) in an associated norm. In the same norm we also establish a supercloseness result. Numerical tests support our theoretical estimates.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

On the coupling of finite volume and discontinuous Galerkin method for elliptic problems

The coupling of cell-centered finite volume method with primal discontinuous Galerkin method is introduced in this paper for elliptic problems. Convergence of the method with respect to the mesh size is proved. Numerical examples confirm the theoretical rates of convergence. Advantages of the coupled scheme are shown for problems with discontinuous coefficients or anisotropic diffusion matrix.     

A uniformly convergent continuous-discontinuous Galerkin method for singularly perturbed problems of convection-diffusion type

In this paper, we introduce a coupled approach of local discontinuous Galerkin and standard finite element method for solving singularly perturbed convection-diffusion problems. On Shishkin mesh with linear elements, a rate O(N^-1lnN) in an associated norm is established, where N is the number of elements. Numerical experiments complement the theoretical results. Moreover, a rate O(N^-2ln²N) in a discrete L^∞ norm, and O(N^-2) in L² norm, are observed numerically on the Shishkin mesh.     

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     

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.     

奇异摄动问题最优阶一致收敛的间断有限元分析

采用非对称内罚间断有限元方法(以下简称NIPG方法)求解一维对流扩散型奇异摄动问题.理论上证明了采用拉格朗日线性元的NIPG方法在Bakhvalov-Shishkin网格上具有最优阶的一致收敛性,即在能量范数度量下其误差估计为O(N~(-1)),其中N为网格剖分中单元个数.数值算例验证了理论分析的正确性.     

Analysis of a DG method for singularly perturbed convection-diffusion problems

In this article, we studied a discontinuous Galerkin finite element method for convection-diffusion-reaction problems with singular perturbation. Our approach is highly flexible by allowing the use of discontinuous approximating function on polytopal mesh without imposing extra conditions on the convection coefficient. A priori error estimate is devised in a suitable energy norm on general polytopal mesh. Numerical examples are provided.     

Construction of locally conservative fluxes for the SUPG method

We consider the construction of locally conservative fluxes by means of a simple postprocessing technique obtained from the finite element solutions of advection diffusion equations. It is known that a naive calculation of fluxes from these solutions yields nonconservative fluxes. We consider two finite element methods: the usual continuous Galerkin finite element method for solving nondominating advection diffusion equations and the streamline upwind/Petrov‐Galerkin method for solving advection dominated problems. We then describe the postprocessing technique for constructing conservative fluxes from the numerical solutions of the general variational formulation. The postprocessing technique requires solving an auxiliary Neumann boundary value problem on each element independently and it produces a locally conservative flux on a vertex centered dual mesh relative to the finite element mesh. We provide a convergence analysis for the postprocessing technique. Performance of the technique and the convergence behavior are demonstrated through numerical examples including a set of test problems for advection diffusion equations, advection dominated equations, and drift‐diffusion equations. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1971–1994, 2015     

Alternating Direction Implicit Galerkin Finite Element Method for the Two-Dimensional Time Fractional Evolution Equation

New numerical techniques are presented for the solution of the two-dimensional time fractional evolution equation in the unit square. In these methods, Galerkin finite element is used for the spatial discretization, and, for the time stepping, new alternating direction implicit (ADI) method based on the backward Euler method combined with the first order convolution quadrature approximating the integral term are considered. The ADI Galerkin finite element method is proved to be convergent in time and in the $L^2$ norm in space. The convergence order is$\mathcal{O}$($k$|ln $k$|+$h^r$), where $k$ is the temporal grid size and $h$ is spatial grid size in the $x$ and $y$ directions, respectively. Numerical results are presented to support our theoretical analysis.     

裂缝-孔隙介质中地下水污染问题的Galerkin交替方向有限元方法

考虑裂缝 孔隙介质中地下水污染问题均匀化模型的数值模拟．对压力方程采用混合元方法,对浓度方程采用Ｇａｌｅｒｋｉｎ交替方向有限元方法,对吸附浓度方程采用标准Ｇａｌｅｒｋｉｎ方法,证明了交替方向有限元格式具有最优犔２ 和犎１ 模误差估计．     

A Characteristic Finite Volume Element Method for the Air Pollution Model

In this article, a characteristic finite volume element method is presented for solving air pollution models. The convection term is discretized using the characteristic method and diffusion term is approximated by finite volume element method. Compared with standard finite volume element method, our proposed method is more accurate and efficient, especially suitable to solve convection-dominated problems. The proposed numerical schemes are analyzed for convergence in L ₂ norm. Some numerical results are presented to demonstrate the efficiency and accuracy of the method.     

A parameter robust Petrov–Galerkin scheme for advection–diffusion–reaction equations

In this paper, a singularly perturbed convection diffusion boundary value problem, with discontinuous diffusion coefficient is examined. In addition to the presence of boundary layers, strong and weak interior layers can also be present due to the discontinuities in the diffusion coefficient. A priori layer adapted piecewise uniform meshes are used to resolve any layers present in the solution. Using a Petrov–Galerkin finite element formulation, a fitted finite difference operator is shown to produce numerical approximations on this fitted mesh, which are uniformly second order (up to logarithmic terms) globally convergent in the pointwise maximum norm.     

A New Hybridized Mixed Weak Galerkin Method for Second-Order Elliptic Problems

In this paper, a new hybridized mixed formulation of weak Galerkin method is studied for a second order elliptic problem. This method is designed by approximate some operators with discontinuous piecewise polynomials in a shape regular finite element partition. Some discrete inequalities are presented on discontinuous spaces and optimal order error estimations are established. Some numerical results are reported to show super convergence and confirm the theory of the mixed weak Galerkin method.     

Interior penalty discontinuous approximations of convection–diffusion problems with parabolic layers

Summary A nonsymmetric discontinuous Galerkin finite element method with interior penalties is considered for two–dimensional convection–diffusion problems with regular and parabolic layers. On an anisotropic Shishkin–type mesh with bilinear elements we prove error estimates (uniformly in the perturbation parameter) in an integral norm associated with this method. On different types of interelement edges we derive the values of discontinuity–penalization parameters. Numerical experiments complement the theoretical results.