A mixed finite element for the stokes problem using quadrilateral elements

The object of this paper is to complete the results obtained in [3] by showing that the new mixed finite element that we have constructed in [3] also works for quadrilateral elements and to compare this method with the standard finite volume method. Estimates of optimal order are derived for both the new mixed finite element and an associated finite volume method.     

Multiscale methods for elliptic homogenization problems

In this article we study two families of multiscale methods for numerically solving elliptic homogenization problems. The recently developed multiscale finite element method [Hou and Wu, J Comp Phys 134 (1997), 169–189] captures the effect of microscales on macroscales through modification of finite element basis functions. Here we reformulate this method that captures the same effect through modification of bilinear forms in the finite element formulation. This new formulation is a general approach that can handle a large variety of differential problems and numerical methods. It can be easily extended to nonlinear problems and mixed finite element methods, for example. The latter extension is carried out in this article. The recently introduced heterogeneous multiscale method [Engquist and Engquist, Comm Math Sci 1 (2003), 87–132] is designed for efficient numerical solution of problems with multiscales and multiphysics. In the second part of this article, we study this method in mixed form (we call it the mixed heterogeneous multiscale method). We present a detailed analysis for stability and convergence of this new method. Estimates are obtained for the error between the homogenized and numerical multiscale solutions. Strategies for retrieving the microstructural information from the numerical solution are provided and analyzed. Relationship between the multiscale finite element and heterogeneous multiscale methods is discussed. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

A nonconforming mixed finite element for second-order elliptic problems

In a recent work, Hiptmair [Mathematisches Institut, M9404, 1994] has constructed and analyzed a family of nonconforming mixed finite elements for second-order elliptic problems. However, his analysis does not work on the lowest order elements. In this article, we show that it is possible to construct a nonconforming mixed finite element for the lowest order case. We prove the convergence and give estimates of optimal order for this finite element. Our proof is based on the use of the properties of the so-called nonconforming bubble function to control the consistency terms introduced by the nonconforming approximation. We further establish an equivalence between this mixed finite element and the nonconforming piecewise quadratic finite element of Fortin and Soulie [J. Numer. Methods Eng., 19, 505–520, 1983]. © 1997 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 13: 445–457, 1997     

A mixed finite element method for a quasi‐Newtonian fluid flow

We propose a mixed formulation for quasi‐Newtonian fluid flow obeying the power law where the stress tensor is introduced as a new variable. Based on such a formulation, a mixed finite element is constructed and analyzed. This finite element method possesses local (i.e., at element level) conservation properties (conservation of the momentum and the mass) as in the finite volume methods. We give existence and uniqueness results for the continuous problem and its approximation and we prove error bounds. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004.     

二阶椭圆问题的稳定化杂交有限元分析

0　引　　言Raviart&Thomas(1977)[13]基于Babǔska-Brezzi有限元理论[1][5]发展了二阶椭圆问题的基本杂交方法.该文指出,为确定合适的自由度,一般将杂交元刻划为非协调元.然而,对三角形偶数次杂交元和四边形杂交元而言,[13]是通过扩充手段克服有限维空间“匹配”问题的.由于扩充元的复杂性及其不再能刻划为非协调元,以致于实际计算无法选取自由度.Thomas的博士论文[15]提供了一个解决办法.即利用Gauss-Legendre数值求积分公式将扩充元近似刻划成非协调元,得到数值积分意义下的杂交方法.如此处理虽然大大简化了原杂交格式的求解过程,但数…     

Finite volumes and mixed Petrov‐Galerkin finite elements: The unidimensional problem

For Laplace operator in one space dimension, we propose to formulate the heuristic finite volume method with the help of mixed Petrov‐Galerkin finite elements. Weighting functions for gradient discretization are parameterized by some function ψ : [0, 1] → ℝ. We propose for this function ψ a compatibility interpolation condition, and we prove that such a condition is equivalent to the inf‐sup property when studying stability of the numerical scheme. In the case of stable scheme and under two distinct hypotheses concerning the regularity of the solution, we demonstrate convergence of the finite volume method in appropriate Hilbert spaces and with optimal order of accuracy. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 335–360, 2000     

A New Finite Element Space for Expanded Mixed Finite Element Method

In this article, we propose a new finite element space Λ$_h$ for the expanded mixed finite element method (EMFEM) for second-order elliptic problems to guarantee its computing capability and reduce the computation cost. The new finite element space Λ$_h$ is designed in such a way that the strong requirement V$_h\subset$Λ$_h$ in [9] is weakened to {v$_h\in$V$_h$; divv$_h$=0}$\subset$Λ$_h$ so that it needs fewer degrees of freedom than its classical counterpart. Furthermore, the new Λ$_h$ coupled with the Raviart-Thomas space satisfies the inf-sup condition, which is crucial to the computation of mixed methods for its close relation to the behavior of the smallest nonzero eigenvalue of the stiff matrix, and thus the existence, uniqueness and optimal approximate capability of the EMFEM solution are proved for rectangular partitions in $\mathbb{R}^d, d=2,3$ and for triangular partitions in $\mathbb{R}^2$. Also, the solvability of the EMFEM for triangular partition in $\mathbb{R}^3$ can be directly proved without the inf-sup condition. Numerical experiments are conducted to confirm these theoretical findings.     

不可压缩流问题低次元稳定有限体积数值方法研究

分析了R^d,d=2,3维不可压缩流Stokes问题低次元稳定有限体积方法,它主要利用局部压力投影方法对两种流行但不满足inf-sup条件的有限元配对(P_1-P_0和P_1-P_1)在有限体积方法的框架下进行稳定;利用有限元与有限体积方法的等价性进行有限体积方法理论分析.结果表明不可压缩流Stokes问题在f∈Hd,d=2,3维不可压缩流Stokes问题低次元稳定有限体积方法,它主要利用局部压力投影方法对两种流行但不满足inf-sup条件的有限元配对(P_1-P_0和P_1-P_1)在有限体积方法的框架下进行稳定;利用有限元与有限体积方法的等价性进行有限体积方法理论分析.结果表明不可压缩流Stokes问题在f∈H¹情况下,本文方法得到的解与稳定有限元方法解之间具有O(h1情况下,本文方法得到的解与稳定有限元方法解之间具有O(h²)阶超收敛阶结果,且稳定有限体积方法取得了与稳定有限元方法相同的收敛速度,与稳定有限元方法比较,稳定有限体积方法计算简单高效,同时保持物理守恒,因此在实际应用中具有很好的潜力。     

Development and Analysis of Higher Order Finite Volume Methods over Rectangles for Elliptic Equations

Currently used finite volume methods are essentially low order methods. In this paper, we present a systematic way to derive higher order finite volume schemes from higher order mixed finite element methods. Mostly for convenience but sometimes from necessity, our procedure starts from the hybridization of the mixed method. It then approximates the inner product of vector functions by an appropriate, critical quadrature rule; this allows the elimination of the flux and Lagrange multiplier parameters so as to obtain equations in the scalar variable, which will define the finite volume method. Following this derivation with different mixed finite element spaces leads to a variety of finite volume schemes. In particular, we restrict ourselves to finite volume methods posed over rectangular partitions and begin by studying an efficient second-order finite volume method based on the Brezzi–Douglas–Fortin–Marini space of index two. Then, we present a general global analysis of the difference between the solution of the underlying mixed finite element method and its related finite volume method. Then, we derive finite volume methods of all orders from the Raviart–Thomas two-dimensional rectangular elements; we also find finite volume methods to associate with BDFM ₂ three-dimensional rectangles. In each case, we obtain optimal error estimates for both the scalar variable and the recovered flux.     

多孔介质驱动问题的混合元最小二乘特征有限元方法及其收敛分析

１．引言多孔介质二相驱动问题的数学模型是偶合的非线性偏微分方程组的初边值问题．该问题可转化为压力方程和浓度方程［１－４］．浓度方程一般是对流占优的对流扩散方程,它的对流速度依赖于比浓度方程的扩散系数大得多的Ｆａｒｃｙ速度．因此Ｄａｒｃｙ速度的求解精度直接影响着浓度的求解精度．为了提高速度的求解精度,７０年代Ｐ．Ａ．Ｒａｖｉａｔ和Ｊ．Ｍ．Ｔｈｏｍａｓ提出混合有限元方法［５］．Ｊ．ＤｏｕｇｌａｓＪｒ,Ｔ．Ｆ．Ｒｕｓｓｅｌｌ,Ｒ．Ｅ．Ｅｗｉｎｇ,Ｍ．Ｆ．Ｗｈｅｅｌｅｒ［１］－［４］,［９］,［１２］袁…     

ERROR ESTIMATES FOR TWO-SCALE COMPOSITE FINITE ELEMENT APPROXIMATIONS OF NONLINEAR PARABOLIC EQUATIONS

We study spatially semidiscrete and fully discrete two-scale composite finite element method for approximations of the nonlinear parabolic equations with homogeneous Dirich-let boundary conditions in a convex polygonal domain in the plane.This new class of finite elements,which is called composite finite elements,was first introduced by Hackbusch and Sauter[Numer.Math.,75(1997),pp.447-472]for the approximation of partial differential equations on domains with complicated geometry.The aim of this paper is to introduce an efficient numerical method which gives a lower dimensional approach for solving par-tial differential equations by domain discretization method.The composite finite element method introduces two-scale grid for discretization of the domain,the coarse-scale and the fine-scale grid with the degrees of freedom lies on the coarse-scale grid only.While the fine-scale grid is used to resolve the Dirichlet boundary condition,the dimension of the finite element space depends only on the coarse-scale grid.As a consequence,the resulting linear system will have a fewer number of unknowns.A continuous,piecewise linear composite finite element space is employed for the space discretization whereas the time discretization is based on both the backward Euler and the Crank-Nicolson methods.We have derived the error estimates in the L∞(L2)-norm for both semidiscrete and fully discrete schemes.Moreover,numerical simulations show that the proposed method is an efficient method to provide a good approximate solution.     

Analysis of a non-standard mixed finite element method with applications to superconvergence

We show that a non-standard mixed finite element method proposed by Barrios and Gatica in 2007, is a higher order perturbation of the least-squares mixed finite element method. Therefore, it is also superconvergent whenever the least-squares mixed finite element method is superconvergent. Superconvergence of the latter was earlier investigated by Brandts, Chen and Yang between 2004 and 2006. Since the new method leads to a non-symmetric system matrix, its application seems however more expensive than applying the least-squares mixed finite element method. Dedicated to Ivan Hlaváček on the occasion of his 75th birthday     

Remarks on Multilevel Bases for Divergence-Free Finite Elements

The connection between the multilevel factorization method recently proposed by Sarin and Sameh for solving mixed discretizations of the Stokes equation using a divergence-free finite element formulation, and hierarchical basis preconditioners for the Poisson problem is established. For the 2D triangular Taylor–Hood element, a preconditioner is proposed that could be useful in fractional step methods.     

A posteriori error estimator for finite volume methods

In this paper, first we discuss a technique to compare finite volume method and some well-known finite element methods, namely the dual mixed methods and nonconforming primal methods, for elliptic equations. These both equivalences are exploited to give us a posteriori error estimator for finite volume methods. This estimator is explicitly given, easy to compute and asymptotically exact without any regularity of the solution in unstructured grids.     

Coupling of mixed finite element and stabilized boundary element methods for a fluid–solid interaction problem in 3D

We introduce and analyze the coupling of a mixed finite element and a boundary element for a three‐dimensional time‐harmonic fluid–solid interaction problem. We consider a formulation in which the Cauchy stress tensor and the rotation are the main variables in the elastic structure and use the usual pressure formulation in the acoustic fluid. The mixed variational formulation in the solid is completed with boundary integral equations relating the Cauchy data of the acoustic problem on the coupling interface. A crucial point in our formulation is the stabilization technique introduced by Hiptmair and coworkers to avoid the well‐known instability issue appearing in the boundary element method treatment of the exterior Helmholtz problem. The main novelty of this formulation, with respect to a previous approach, consists in reducing the computational domain to the solid media and providing a more accurate treatment of the far field effect. We show that the continuous problem is well‐posed and propose a conforming Galerkin method based on the lowest‐order Arnold–Falk–Winther mixed finite element. Finally, we prove that the numerical scheme is convergent with optimal order.Copyright © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1211–1233, 2014     

本征值有限元近似的一个抽象误差估计式

设T：LZ（fi）MLZ（fi）是自共轭全连续算子，SgCLZ（fi）是分片。次有限元空间，几：LZ、St是有限秩自共轭算子，IITh-Tllo、0（h、0）．考虑本征值问题：及其近似用（·，·）和DD·D【。·分别表示h（m中内积和范数·ID·卜F表示w认｝（m中范数，简记D卜队。为D卜卜·因为T自共轭全连续，所以它有可数无穷个本征值h，人，...人，...．其相应的本征函数（2丹构成完全标准直交系，所以VZELZ（m设几的本征值为A。l，汕。，...，汕n，相应的本征函数为如山，卜则。二1·不失一般性，可EitL。设tik一大干二＞，E是到Ah对应的本…     

平面线弹性纯位移边值问题低阶非协调矩形元的Robust多重网格方法

1引 言 对于各向同性,均匀介质的平面线弹性问题,当Lamé常数λ→∞(泊松率v→0.5)时,即对于几乎不可压介质,通常的协调有限元格式的解往往不再收敛到原问题的解,或者达不到最优收敛阶,这就是所谓的闭锁现象(见[3],[7],[8]及[10]).究其原因,在通常的有限元分析中,其误差估计的系数与λ有关,当λ→∞时,该系数将趋于无穷大.因此为克服闭锁现象就需要构造特殊的有限元格式,使得当λ→∞时,有限元逼近解仍然收敛到原问题的解.     

解重调和问题混合有限元方程的直接方法

§1.引言 考虑如下重调和方程的齐次边值问题: △~2w=f,在Ω中, w=?w/?v=0,在Ω上.(1.1)其中Ω是平面凸多边形区域,?Ω是Ω的边界,?/?v表示?Ω上的外法向导数.     

A mixed finite volume method for elliptic problems

We derive a novel finite volume method for the elliptic equation, using the framework of mixed finite element methods to discretize the pressure and velocities on two different grids (covolumes), triangular (tetrahedral) mesh and control volume mesh. The new discretization is defined for tensor diffusion coefficient and well suited for heterogeneous media. When the control volumes are created by connecting the center of gravity of each triangle to the midpoints of its edges, we show that the discretization is stable and first order accurate for both scalar and vector unknowns. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Semiregular Finite Elements in Solving Some Nonlinear Problems

In this paper, under the maximum angle condition, the finite element method is analyzed for nonlinear elliptic variational problem formulated in [4]. In [4] the analysis was done under the minimum angle condition.