A mixed finite element method for solving the nonstationary stokes equation

Summary This paper deals with a mixed finite element method for approximating a fourth order initial value problem arising from the nonstationary Stokes problem. For piecewise linear shape functions error estimates are given with convergence rates similar to the elliptic case. Some numerical computations will illustrate the theoretical results.     

解三维抛物型方程的一个ADI格式

构造了一个解三维抛物型方程的高精度ADI格式,格式绝对稳定,截断误差为O（△t^2＋△x^4）;然后应用Richerdson外推法,外推一次得到了具有O（△t^3＋△x^6）阶精度的近似解.     

解高维热传导方程的一个高精度ADI格式

构造了一个解四维热传导方程的一个高精度ADI格式,格式绝对稳定,截断误差阶达到O(△t~2 △x~4).可用追赶法求解.     

A new mixed finite element method for Burgers’ equation

In this paper, anH ¹-Galerkin mixed finite element method is used to approximate the solution as well as the flux of Burgers’ equation. Error estimates have been derived. The results of the numerical experiment show the efficacy of the mixed method and justifies the theoretical results obtained in the paper.     

Analysis of mixed finite element methods on locally refined grids

Summary We consider the mixed finite element method for locally refined triangulations. A local projection operator is defined to satisfy the commutation property that is required in the general theory of mixed methods. Our results can be applied to every known space of arbitrary order over rectangles or triangles. Optimal-order error estimates and superconvergence for the flux along the Gauss lines are established.     

二次Lagrangian有限元方程的几何多重网格法

为了构造快速求解二次Lagrangian有限元方程的几何多重网格法,在选择二次Lagrangian有限元空间和一系列线性Lagrangian有限元空间分别作为最细网格层和其余粗网格层以及构造一种新限制算子的基础上,提出了一种新的几何多重网格法,并对它的计算量进行了估计.数值实验结果,与通常的几何多重网格法和AMG01法相比,表明了新算法计算量少且稳健性强.     

A finite difference scheme for solving the heat transport equation at the microscale

Heat transport at the microscale is of vital importance in microtechnology applications. In this study, we develop a finite difference scheme of the Crank‐Nicholson type by introducing an intermediate function for the heat transport equation at the microscale. It is shown by the discrete energy method that the scheme is unconditionally stable. Numerical results show that the solution is accurate. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 697–708, 1999     

A mixed hybrid formulation based on oscillated finite element polynomials for solving Helmholtz problems

A mixed-hybrid-type formulation is proposed for solving Helmholtz problems. This method is based on (a) a local approximation of the solution by oscillated finite element polynomials and (b) the use of Lagrange multipliers to “weakly” enforce the continuity across element boundaries. The computational complexity of the proposed discretization method is determined mainly by the total number of Lagrange multiplier degrees of freedom introduced at the interior edges of the finite element mesh, and the sparsity pattern of the corresponding system matrix. Preliminary numerical results are reported to illustrate the potential of the proposed solution methodology for solving efficiently Helmholtz problems in the mid- and high-frequency regimes.     

A finite element lumped mass scheme for solving eigenvalue problems of circular arches

Summary In this paper, we present a finite element lumped mass scheme for eigenvalue problems of circular arch structures, and give error estimates for the approximation. They assert that approximate eigenvalues and eigenfuctions converge to the exact ones. Some numerical examples are also given to illustrate our results.     

A finite element lumped mass scheme for solving eigenvalue problems of circular arches

Summary This study is a continuation of a previous paper [4] in which the numerical results are given by using single precision arithmetic. In this paper, we show the numerical results which experess the sharper convergence properties than those of [4], by using double precision arithmetic.Dedicated to Prof. Masaya Yamaguti on the occasion of his 60th birthday     

A finite element scheme for domains with corners

This paper is concerned with the rate of convergence of the finite element method on polygonal domains in weighted Sobolev spaces. It is shown that the use of different spaces of trial and test functions will restrict the usual low rate of convergence to a neighborhood of each vertex of the polygonal domain.L ₂-convergence and lower bounds on the error are also studied.This research was supported in part by the Atomic Energy Commission under contract no. AEC AT-(40-1)-3443/4.This research was supported in part by the U.S. Naval Academy Research Council.     

Sobolev方程一个新的H~1-Galerkin混合有限元分析

研究了Sobolev方程的H~1-Galerkin混合有限元方法.利用不完全双二次元Q_2~-和一阶BDFM元,建立了一个新的混合元模式,通过Bramble-Hilbert引理,证明了单元对应的插值算子具有的高精度结果.进一步,对于半离散和向后欧拉全离散格式,分别导出了原始变量u在H~1-模和中间变量p在H(div)-模意义下的超逼近性质.     

A stabilized mixed finite element method for the biharmonic equation based on biorthogonal systems

We propose a stabilized finite element method for the approximation of the biharmonic equation with a clamped boundary condition. The mixed formulation of the biharmonic equation is obtained by introducing the gradient of the solution and a Lagrange multiplier as new unknowns. Working with a pair of bases forming a biorthogonal system, we can easily eliminate the gradient of the solution and the Lagrange multiplier from the saddle point system leading to a positive definite formulation. Using a superconvergence property of a gradient recovery operator, we prove an optimal a priori estimate for the finite element discretization for a class of meshes.     

对流占优扩散方程的分裂特征混合有限元方法

利用修正的特征线方法,构建一类求解对流占优扩散方程的分裂特征混合有限元算法.在新的算法中,混合系统的系数矩阵对称正定,且原未知函数u与流函数σ=-ε▽u可分离求解.推导了加权能量模意义下的最优阶误差估计,并给出数值算例验证理论上的分析结果.     

A new efficient energy-preserving finite volume element scheme for the improved Boussinesq equation

In this paper, a new linearized energy-preserving Crank-Nicolson finite volume element scheme is derived for the improved Boussinesq equation. The fully discrete scheme can be shown to conserve both mass and energy in the discrete setting. It is proved that the scheme is uniquely solvable and convergent with the rate of order two in a discrete L₂ norm. At last, a series of numerical experiments on typical improved Boussinesq and Sine–Gordon equations are provided to verify our theoretical results and to show the efficiency and accuracy of the proposed scheme.     

A combined finite volume–finite element scheme for the discretization of strongly nonlinear convection–diffusion–reaction problems on nonmatching grids

We propose and analyze in this paper a numerical scheme for nonlinear degenerate parabolic convection–diffusion–reaction equations in two or three space dimensions. We discretize the time evolution, convection, reaction, and source terms on a given grid, which can be nonmatching and can contain nonconvex elements, by means of the cell‐centered finite volume method. To discretize the diffusion term, we construct a conforming simplicial mesh with the vertices given by the original grid and use the conforming piecewise linear finite element method. In this way, the scheme is fully consistent and the discrete solution is naturally continuous across the interfaces between the subdomains with nonmatching grids, without introducing any supplementary equations and unknowns or using any interpolation at the interfaces. We allow for general inhomogeneous and anisotropic diffusion–dispersion tensors, propose two variants corresponding respectively to arithmetic and harmonic averaging, and use the local Péclet upstream weighting in order to only add the minimal numerical diffusion necessary to avoid spurious oscillations in the convection‐dominated case. The scheme is robust, efficient since it leads to positive definite matrices and one unknown per element, locally conservative, and satisfies the discrete maximum principle under the conditions on the simplicial mesh and the diffusion tensor usual in the finite element method. We prove its convergence using a priori estimates and the Kolmogorov relative compactness theorem and illustrate its behavior on a numerical experiment. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

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.     

The Nitsche mortar method for matching grids in a mixed finite element method

It is well-known that nonmatching grids are often used in finite element methods. Usually, grids are being matched along lines or surfaces that divide a domain into subdomains. Such lines or surfaces are called interfaces. The interface matching means the satisfaction of some continuity conditions when crossing the interface. The direct matching procedures fall into three groups: Methods that use Lagrange multipliers, mortar methods based on the Nitsche technique, and penalty methods.     

A new minimization protocol for solving nonlinear Poisson–Boltzmann mortar finite element equation

The nonlinear Poisson–Boltzmann equation (PBE) is a widely-used implicit solvent model in biomolecular simulations. This paper formulates a new PBE nonlinear algebraic system from a mortar finite element approximation, and proposes a new minimization protocol to solve it efficiently. In particular, the PBE mortar nonlinear algebraic system is proved to have a unique solution, and is equivalent to a unconstrained minimization problem. It is then solved as the unconstrained minimization problem by the subspace trust region Newton method. Numerical results show that the new minimization protocol is more efficient than the traditional merit least squares approach in solving the nonlinear system. At least 80 percent of the total CPU time was saved for a PBE model problem. AMS subject classification (2000)  65N30, 65H10, 65K10, 92-08     

A combined finite element and Marker and cell method for solving Navier-Stokes equations

Summary A new method is developed to approximate Navier-Stokes equations for incompressible fluids on polygonal domains. This method is based on a simple finite element approximation adapted to the Marker and Cell technique. Its main originality lies in the fact that the components of the velocity are approximated on distinct interlaced networks. The error analysis shows that this method is of order one.