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.     

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 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.     

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.     

A Two-Grid Algorithm of Fully Discrete Galerkin Finite Element Methods for a Nonlinear Hyperbolic Equation

A two-grid finite element approximation is studied in the fully discrete scheme obtained by discretizing in both space and time for a nonlinear hyperbolic equation. The main idea of two-grid methods is to use a coarse-grid space ($S_H$) to produce a rough approximation for the solution of nonlinear hyperbolic problems and then use it as the initial guess on the fine-grid space ($S_h$). Error estimates are presented in $H^1$-norm, which show that two-grid methods can achieve the optimal convergence order as long as the two different girds satisfy $h$ = $\mathcal{O}$($H^2$). With the proposed techniques, we can obtain the same accuracy as standard finite element methods, and also save lots of time in calculation. Theoretical analyses and numerical examples are presented to confirm the methods.     

An Error Analysis of Discontinuous Finite Element Methods for the Optimal Control Problems Governed by Stokes Equation

In this article, an abstract framework for the error analysis of discontinuous finite element method is developed for the distributed and Neumann boundary control problems governed by the stationary Stokes equation with control constraints. A priori error estimates of optimal order are derived for velocity and pressure in the energy norm and the L²-norm, respectively. Moreover, a reliable and efficient a posteriori error estimator is derived. The results are applicable to a variety of problems just under the minimal regularity possessed by the well-posedness of the problem. In particular, we consider the abstract results with suitable stable pairs of velocity and pressure spaces like as the lowest-order Crouzeix–Raviart finite element and piecewise constant spaces, piecewise linear and constant finite element spaces. The theoretical results are illustrated by the numerical experiments.     

Schr(o)dinger方程的时空有限元方法与守恒性

对非线性Schroedinger常微分方程，利用常微分方程连续有限元法证明了能量守恒；对非线性Schroedinger偏微分方程利用时空都连续的全离散有限元方法证明了能量积分守恒和利用空间连续、时间问断的有限元法得到电荷近似守恒，误差为高阶量．并在数值计算上探讨了守恒性和近似程度。结果与理论相吻合．     

二阶双曲方程一个新的非协调混合元超收敛性分析

研究了一类二阶双曲型方程在新混合元格式下的非协调混合有限元方法.在抛弃传统有限元分析的必要工具-Ritz投影算子的前提下,直接利用单元的插值性质,运用高精度分析和对时间t的导数转移技巧,借助于插值后处理技术,分别导出了关于原始变量u的H~1-模和通量=-▽u在L~2-模下的O(h~2)阶超逼近性质和整体超收敛结果.进一步,给出了一些数值算例验证了理论分析的正确性.     

A Priori and A Posteriori Error Estimates of Streamline Diffusion Finite Element Method for Optimal Control Problem Governed by Convection Dominated Diffusion Equation

In this paper,we investigate a streamline diffusion finite element approxi- mation scheme for the constrained optimal control problem governed by linear con- vection dominated diffusion equations.We prove the existence and uniqueness of the discretized scheme.Then a priori and a posteriori error estimates are derived for the state,the co-state and the control.Three numerical examples are presented to illustrate our theoretical results.     

Galerkin Finite Element Approximation of General Linear Second Order Hyperbolic Equations

In this article, we derive error estimates for the Galerkin approximation of a general linear second order hyperbolic equation. The results can be applied to a variety of cases, for example, vibrating systems of linked elastic bodies. The results generalize the work of Baker [1 G. A. Baker ( 1976 ). Error estimates for finite element methods for second order hyperbolic equations . SIAM J. Numer. Anal. 13 : 564 – 576 .[Crossref], [Web of Science ®] , [Google Scholar]] and also allow for viscous type damping. Splitting the proofs for the semi-discrete and fully discrete cases not only simplifies the proofs but less restrictive regularity assumptions are required.     

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

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

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.     

Error Estimates for Semi-Discrete and Fully Discrete Galerkin Finite Element Approximations of the General Linear Second-Order Hyperbolic Equation

In this article, we derive error estimates for the semi-discrete and fully discrete Galerkin approximations of a general linear second-order hyperbolic partial differential equation with general damping (which includes boundary damping). The results can be applied to a variety of cases (e.g. vibrating systems of linked elastic bodies). The results generalize pioneering work of Dupont and complement a recent article by Basson and Van Rensburg.     

求解变系数Cahn-Hilliard-Brinkman方程有限元方法的误差分析

本文研究求解变系数Cahn-Hilliard-Brinkman方程有限元方法的误差分析.在时间格式上采用能量凸分裂法以及在空间格式上采用混合有限元法进行离散,证明了全离散格式是能量衰减的.在误差分析中,利用Cauchy中值定理将含浓度和Peclet数的项分解为两项,结果表明所提出的格式在时间上是二阶精度的.     

线性Poisson-Boltzmann方程的Mortar有限元方法的数值计算

本文对分子生物物理学中产生的线性Poisson-Boltzmann方程（PBE），给出了Mortar有限元方法的计算过程，数值计算例子表明，与一般的协调有限元方法相比，用Mortar元方法求解此类有间断系数的问题有非常有效的。     

二维蜂窝结构热方程的多尺度有限元算法

研究二维蜂窝结构上的热方程,提出多尺度有限元高精度算法,并给出严格的误差估计.     

Mortar Upwind Finite Volume Element Method with Crouzeix-Raviart Element for Parabolic Convection Diffusion Problems

In this paper,we study the semi-discrete mortar upwind finite volume element method with the Crouzeix-Raviart element for the parabolic convection diffusion problems. It is proved that the semi-discrete mortar upwind finite volume element approximations derived are convergent in the H~1-and L~2-norms.     

线性对流占优扩散方程的扩展特征混合有限元法

本文针对线性对流占优扩散方程提出了一种新型数值模拟方法一扩展特征混合有限元法,即对对流部分沿特征线方向离散,而对扩散部分采用扩展混合有限元方法,同时高精度逼近未知函数,未知函数的梯度及伴随向量函数,通过严格的数值分析,得到其最优L^2模误差估计。     

求解二次Lagrangian有限元方程的代数两水平方法

基于矩阵图集的粗化算法,构造一种新的插值算子,提出了瀑布型代数两重网格法;然后结合部分几何信息,提出了求解二次Lagrangian有限元方程的代数两水平方法.数值实验表明该算法稳健性强、计算量更少.     

具有对数势能的Cahn-Hilliard方程的有限元算法

本文我们提出了具有对数势能的Cahn-Hilliard方程,在空间上采用混合有限元方法进行离散,时间上采用Crank-Nicolson格式.运用正则性,将对数势能函数F(u)的定义域的范围由(-1,1)扩展到(-∞,∞).证明该方法是能量耗散的,并计算误差估计,最后通过数值算例对理论部分进行验证.结果表明,理论与数值算...