Constructing Order Two Superconvergent WG Finite Elements on Rectangular Meshes

In this paper, we introduce a stabilizer free weak Galerkin (SFWG) finite element method for second order elliptic problems on rectangular meshes. With a special weak Gradient space, an order two superconvergence for the SFWG finite element solution is obtained, in both $L^2$ and $H^1$ norms. A local post-process lifts such a $P_k$ weak Galerkin solution to an optimal order $P_{k+2}$ solution. The numerical results confirm the theory.     

Energy Identities and Stability Analysis of the Yee Scheme for 3D Maxwell Equations

In this paper numerical energy identities of the Yee scheme on uniform grids for three dimensional Maxwell equations with periodic boundary conditions are proposed and expressed in terms of the $L^2$, $H^1$ and $H^2$ norms. The relations between the $H^1$ or $H^2$ semi-norms and the magnitudes of the curls or the second curls of the fields in the Yee scheme are derived. By the $L^2$ form of the identity it is shown that the solution fields of the Yee scheme is approximately energy conserved. By the $H^1$ or $H^2$ semi norm of the identities, it is proved that the curls or the second curls of the solution of the Yee scheme are approximately magnitude (or energy)-conserved. From these numerical energy identities, the Courant-Friedrichs-Lewy (CFL) stability condition is re-derived, and the stability of the Yee scheme in the $L^2$, $H^1$ and $H^2$ norms is then proved. Numerical experiments to compute the numerical energies and convergence orders in the $L^2$, $H^1$ and $H^2$ norms are carried out and the computational results confirm the analysis of the Yee scheme on energy conservation and stability analysis.     

A weak Galerkin method for second order elliptic problems with polynomial reduction

The second order elliptic equation, which is also know as the diffusion-convection equation, is of great interest in many branches of physics and industry. In this paper, we use the weak Galerkin finite element method to study the general second order elliptic equation. A weak Galerkin finite element method is proposed and analyzed. This scheme features piecewise polynomials of degree $k\geq 1$ on each element and piecewise polynomials of degree $k-1\geq 0$ on each edge or face of the element. Error estimates of optimal order of convergence rate are established in both discrete $H^1$ and standard $L^2$ norm. The paper also presents some numerical experiments to verify the efficiency of the method.     

Finite difference/$H^1$-Galerkin MFE procedure for a fractional water wave model

In this article, an $H^1$-Galerkin mixed finite element (MFE) method for solving the time fractional water wave model is presented. First-order backward Euler difference method and $L1$ formula are applied to approximate integer derivative and Caputo fractional derivative with order $1/2$, respectively, and $H^1$-Galerkin mixed finite element method is used to approximate the spatial direction. The analysis of stability for fully discrete mixed finite element scheme is made and the optimal space-time orders of convergence for two unknown variables in both $H^1$-norm and $L^2$-norm are derived. Further, some computing results for a priori analysis and numerical figures based on four changed parameters in the studied problem are given to illustrate the effectiveness of the current method     

On the Construction and Analysis of Finite Volume Element Schemes with Optimal $L^2$ Convergence Rate

We provide a general construction method for a finite volume element (FVE) scheme with the optimal $L^2$ convergence rate. The $k$-($k$-1)-order orthogonal condition (generalized) is proved to be a sufficient and necessary condition for a $k$-order FVE scheme to have the optimal $L^2$ convergence rate in 1D, in which the independent dual parameters constitute a ($k$-1)-dimension surface in the reasonable domain in $k$-dimension.In the analysis, the dual strategies in different primary elements are not necessarily to be the same, and they are allowed to be asymmetric in each primary element, which open up more possibilities of the FVE schemes to be applied to some complex problems, such as the convection-dominated problems. It worth mentioning that, the construction can be extended to the quadrilateral meshes in 2D. The stability and $H^1$ estimate are proved for completeness. All the above results are demonstrated by numerical experiments.     

A weak Galerkin finite element method for the second order elliptic problems with mixed boundary conditions

In this paper, a weak Galerkin finite element method is proposed and analyzed for the second-order elliptic equation with mixed boundary conditions. Optimal order error estimates are established in both discrete $H^1$ norm and the standard $L^2$ norm for the corresponding WG approximations. The numerical experiments are presented to verify the efficiency of the method.     

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

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

A Two-Grid Finite Element Method for Time-Dependent Incompressible Navier-Stokes Equations with Non-Smooth Initial Data

We analyze here, a two-grid finite element method for the two dimensional time-dependent incompressible Navier-Stokes equations with non-smooth initial data. It involves solving the non-linear Navier-Stokes problem on a coarse grid of size $H$ and solving a Stokes problem on a fine grid of size $h, h <相似文献   

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.     

Optimal and Pressure-Independent $L^2$ Velocity Error Estimates for a Modified Crouzeix-Raviart Stokes Element with BDM Reconstructions

Nearly all inf-sup stable mixed finite elements for the incompressible Stokes equations relax the divergence constraint. The price to pay is that a priori estimates for the velocity error become pressure-dependent, while divergence-free mixed finite elements deliver pressure-independent estimates. A recently introduced new variational crime using lowest-order Raviart-Thomas velocity reconstructions delivers a much more robust modified Crouzeix-Raviart element, obeying an optimal pressure-independent discrete $H^1$ velocity estimate. Refining this approach, a more sophisticated variational crime employing the lowest-order BDM element is proposed, which also allows proving an optimal pressure-independent $L^2$ velocity error. Numerical examples confirm the analysis and demonstrate the improved robustness in the Navier-Stokes case.     

Error Analysis of SAV Finite Element Method to Phase Field Crystal Model

In this paper, we construct and analyze an energy stable scheme by combining the latest developed scalar auxiliary variable (SAV) approach and linear finite element method (FEM) for phase field crystal (PFC) model, and show rigorously that the scheme is first-order in time and second-order in space for the $L^2$ and $H^{-1}$ gradient flow equations. To reduce efficiently computational cost and capture accurately the phase interface, we give a simple adaptive strategy, equipped with a posteriori gradient estimator, i.e., $L^2$ norm of the recovered gradient. Extensive numerical experiments are presented to verify our theoretical results and to demonstrate the effectiveness and accuracy of our proposed method.     

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.     

PPIFE Method with Non-Homogeneous Flux Jump Conditions and Its Efficient Numerical Solver for Elliptic Optimal Control Problems with Interfaces

In this paper, we design a partially penalized immersed finite element method for solving elliptic interface problems with non-homogeneous flux jump conditions. The method presented here has the same global degrees of freedom as classic immersed finite element method. The non-homogeneous flux jump conditions can be handled accurately by additional immersed finite element functions. Four numerical examples are provided to demonstrate the optimal convergence rates of the method in $L^{\infty}$, $L^{2}$ and $H^{1}$ norms. Furthermore, the method is combined with post-processing technique to solve elliptic optimal control problems with interfaces. To solve the resulting large-scale system, block diagonal preconditioners are introduced. These preconditioners can lead to fast convergence of the Krylov subspace methods such as GMRES and are independent of the mesh size. Four numerical examples are presented to illustrate the efficiency of the numerical schemes and preconditioners.     

Two-Grid Finite Element Method with Crank-Nicolson Fully Discrete Scheme for the Time-Dependent Schrödinger Equation

In this paper, we study the Crank-Nicolson Galerkin finite element method and construct a two-grid algorithm for the general two-dimensional time-dependent Schrödinger equation. Firstly, we analyze the superconvergence error estimate of the finite element solution in $H^1$ norm by use of the elliptic projection operator. Secondly, we propose a fully discrete two-grid finite element algorithm with Crank-Nicolson scheme in time. With this method, the solution of the Schrödinger equation on a fine grid is reduced to the solution of original problem on a much coarser grid together with the solution of two Poisson equations on the fine grid. Finally, we also derive error estimates of the two-grid finite element solution with the exact solution in $H^1$ norm. It is shown that the solution of two-grid algorithm can achieve asymptotically optimal accuracy as long as mesh sizes satisfy $H = \mathcal{O}(h^{\frac{1}{2}})$.     

An Adaptive Mesh-Refining Algorithm Allowing for an <Emphasis Type="Italic">H</Emphasis><Superscript>1</Superscript> Stable <Emphasis Type="Italic">L</Emphasis><Superscript>2</Superscript> Projection onto Courant Finite Element Spaces

Suppose $\cal{S}^1({\cal T})\subset H^1(\Omega)$ is the $P_1$-finite element space of $\cal{T}$-piecewise affine functions based on a regular triangulation $\cal{T}$ of a two-dimensional surface $\Omega$ into triangles. The $L^2$ projection $\Pi$ onto $\cal{S}^1(\cal{T})$ is $H^1$ stable if $\norm{\Pi v}{H^1(\Omega)}\le C\norm{v}{H^1(\Omega)}$ for all $v$ in the Sobolev space $H^1(\Omega)$ and if the bound $C$ does not depend on the mesh-size in $\cal{T}$ or on the dimension of $\cal{S}^1(\cal{T})$. \hskip 1em A red–green–blue refining adaptive algorithm is designed which refines a coarse mesh $\cal{T}_0$ successively such that each triangle is divided into one, two, three, or four subtriangles. This is the newest vertex bisection supplemented with possible red refinements based on a careful initialization. The resulting finite element space allows for an $H^1$ stable $L^2$ projection. The stability bound $C$ depends only on the coarse mesh $\cal{T}_0$ through the number of unknowns, the shapes of the triangles in $\cal{T}_0$, and possible Dirichlet boundary conditions. Our arguments also provide a discrete version $\norm{h_\cal{T}^{-1}\,\Pi v}{L^2(\Omega)}\le C\norm{h_\cal{T}^{-1}\,v}{L^2(\Omega)}$ in $L^2$ norms weighted with the mesh-size $h_\T$.     

A C0-WEAK GALERKIN FINITE ELEMENT METHOD FOR THE TWO-DIMENSIONAL NAVIER-STOKES EQUATIONS IN STREAM-FUNCTION FORMULATION

We propose and analyze a $C^0$-weak Galerkin (WG) finite element method for the numerical solution of the Navier-Stokes equations governing 2D stationary incompressible flows. Using a stream-function formulation, the system of Navier-Stokes equations is reduced to a single fourth-order nonlinear partial differential equation and the incompressibility constraint is automatically satisfied. The proposed method uses continuous piecewise-polynomial approximations of degree $k+2$ for the stream-function $\psi$ and discontinuous piecewise-polynomial approximations of degree $k+1$ for the trace of $\nabla\psi$ on the interelement boundaries. The existence of a discrete solution is proved by means of a topological degree argument, while the uniqueness is obtained under a data smallness condition. An optimal error estimate is obtained in $L^2$-norm, $H^1$-norm and broken $H^2$-norm. Numerical tests are presented to demonstrate the theoretical results.     

A New Discontinuous Galerkin Method for Parabolic Equations with Discontinuous Coefficient

In this paper, a new discontinuous Galerkin method is developed for the parabolic equation with jump coefficients satisfying the continuous flow condition. Theoretical analysis shows that this method is $L^2$ stable. When the finite element space consists of interpolative polynomials of degrees $k$, the convergent rate of the semi-discrete discontinuous Galerkin scheme has an order of$\mathcal{O}(h^k)$. Numerical examples for both 1-dimensional and 2-dimensional problems demonstrate the validity of the new method.     

Stabilized Nonconforming Mixed Finite Element Method for Linear Elasticity on Rectangular or Cubic Meshes

Based on the primal mixed variational formulation, a stabilized nonconforming mixed finite element method is proposed for the linear elasticity on rectangular and cubic meshes. Two kinds of penalty terms are introduced in the stabilized mixed formulation, which are the jump penalty term for the displacement and the divergence penalty term for the stress. We use the classical nonconforming rectangular and cubic elements for the displacement and the discontinuous piecewise polynomial space for the stress, where the discrete space for stress are carefully chosen to guarantee the well-posedness of discrete formulation. The stabilized mixed method is locking-free. The optimal convergence order is derived in the $L^2$-norm for stress and in the broken $H^1$-norm and $L^2$-norm for displacement. A numerical test is carried out to verify the optimal convergence of the stabilized method.     

AN ERROR ANALYSIS METHOD SPP-BEAM AND A CONSTRUCTION GUIDELINE OF NONCONFORMING FINITE ELEMENTS FOR FOURTH ORDER ELLIPTIC PROBLEMS

Under two hypotheses of nonconforming finite elements of fourth order elliptic problems,we present a side-patchwise projection based error analysis method(SPP-BEAM for short).Such a method is able to avoid both the regularity condition of exact solutions in the classical error analysis method and the complicated bubble function technique in the recent medius error analysis method.In addition,it is universal enough to admit generalizations.Then,we propose a sufficient condition for these hypotheses by imposing a set of in some sense necessary degrees of freedom of the shape function spaces.As an application,we use the theory to design a P3 second order triangular H2 non-conforming element by enriching two P4 bubble functions and,another P4 second order triangular H2 nonconforming finite element,and a P3 second order tetrahedral H2 non-conforming element by enriching eight P4 bubble functions,adding some more degrees of freedom.     

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.