Optimal convergence analysis of mixed finite element methods for fourth‐order elliptic and parabolic problems

By using a special interpolation operator developed by Girault and Raviart (finite element methods for Navier‐Stokes Equations, Springer‐Verlag, Berlin, 1986), we prove that optimal error bounds can be obtained for a fourth‐order elliptic problem and a fourth‐order parabolic problem solved by mixed finite element methods on quasi‐uniform rectangular meshes. Optimal convergence is proved for all continuous tensor product elements of order k ≥ 1. A numerical example is provided for solving the fourth‐order elliptic problem using the bilinear element. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Estimation de l'erreur d'interpolation d'Hermite dans ℝ n

Summary This paper is devoted to study the Hermite interpolation error in an open subset of ⁿ .It follows a previous work of Arcangeli and Gout [1]. Like this one, it is based principally on the paper of Ciarlet and Raviart [7].We obtain two kinds of the Hermite interpolation error, the first from the Hermite interpolation polynomial, the other from approximation method using the Taylor polynomial.Finally in the last part we study some numerical examples concerning straight finite element methods: in the first and second examples, we use finite elements which are included in the affine theory, but it is not the case in the last example. However, in this case, it is possible to refer to the affine theory by the way of particular study (cf. Argyris et al. [2]; Ciarlet [6]; ciarlet and Raviart [7]; Raviart [11]).

     

On the generation of symmetric Lebesgue-like points in the triangle

We compute point sets on the triangle that have low Lebesgue constant, with sixfold symmetries and Gauss–Legendre–Lobatto distribution on the sides, up to interpolation degree 18. Such points have the best Lebesgue constants among the families of symmetric points used so far in the framework of triangular spectral elements.     

Spectral approximation of variationally-posed eigenvalue problems by nonconforming methods

This paper deals with the nonconforming spectral approximation of variationally posed eigenvalue problems. It is an extension to more general situations of known previous results about nonconforming methods. As an application of the present theory, convergence and optimal order error estimates are proved for the lowest order Crouzeix–Raviart approximation of the eigenpairs of two representative second-order elliptical operators.     

A numerical method based on extended Raviart–Thomas (ER‐T) mixed finite element method for solving damped Boussinesq equation

In this paper, we present the approximate solution of damped Boussinesq equation using extended Raviart–Thomas mixed finite element method. In this method, the numerical solution of this equation is obtained using triangular meshes. Also, for discretization in time direction, we use an implicit finite difference scheme. In addition, error estimation and stability analysis of both methods are shown. Finally, some numerical examples are considered to confirm the theoretical results. Copyright © 2017 John Wiley & Sons, Ltd.     

Constructing C1 triangular patch of degree six by the Boolean sum of an approximation operator and an interpolation operator

A new method to construct C¹ triangular patches which satisfy the given boundary curves and cross-boundary slopes is presented. The Boolean sum of an approximation operator and an interpolation operator is employed to construct the triangular patch. The approximation operator is used to construct a polynomial patch of degree six. The polynomial of degree six affords more freedoms, which makes the approximation operator not only approximate the given boundary interpolation conditions but also have a better approximation precision in the interior of the triangle, so that the triangular patch has a better precision on both the boundary and the interior of the triangular domain. The interpolation operator is utilized to build an interpolation patch which satisfies the given boundary conditions. The Boolean sum of the approximation and interpolation patches forms the triangular patch. Comparison results of the new method with other three methods are given.     

Mesh‐centered finite differences from unconventional mixed‐hybrid nodal finite elements

It is shown how mesh‐centered finite differences can be obtained from unconventional mixed‐hybrid nodal finite elements. The classical Raviart‐Thomas schemes of index k (RTk) are based on interpolation parameters that are cell and/or edge moments. For the unconventional form (URTk), they become point values at Gaussian points. In particular, the scheme URT1 is fully described. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

Approximation properties of lowest-order hexahedral Raviart–Thomas finite elements

Basic interpolation results are settled for lowest-order hexahedral Raviart–Thomas finite elements. Convergence in $\mathrm{H}(\mathrm{div})$ is proved for regular families of asymptotically parallelepiped meshes. The need of the asymptotically parallelepiped assumption is demonstrated with a numerical example. To cite this article: A. Bermúdez et al., C. R. Acad. Sci. Paris, Ser. I 340 (2005).     

Stability of a combined finite element ‐ finite volume discretization of convection‐diffusion equations

We consider a time‐dependent and a stationary convection‐diffusion equation. These equations are approximated by a combined finite element – finite volume method: the diffusion term is discretized by Crouzeix‐Raviart piecewise linear finite elements on a triangular grid, and the convection term by upwind barycentric finite volumes. In the nonstationary case, we use an implicit Euler approach for time discretization. This scheme is shown to be L²‐stable uniformly with respect to the diffusion coefficient. In addition, it turns out that stability is unconditional in the time‐dependent case. These results hold if the underlying grid satisfies a condition that is fulfilled, for example, by some structured meshes. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 402–424, 2012     

A PIC method for solving a paraxial model of highly relativistic beams

Solving the Vlasov–Maxwell problem can lead to very expensive computations. To construct a simpler model, Laval et al. [G. Laval, S. Mas-Gallic, P.A. Raviart, Paraxial approximation of ultrarelativistic intense beams, Numer. Math. 69 (1) (1994) 33–60] proposed to exploit the paraxial property of the charged particle beams, i.e the property that the particles of the beam remain close to an optical axis. They so constructed a paraxial model and performed its mathematical analysis. In this paper, we investigate how their framework can be adapted to handle the axisymmetric geometry, and its coupling with the Vlasov equation. First, one constructs numerical schemes and error estimates results for this discretization are reported. Then, a Particle In Cell (PIC) method, in the case of highly relativistic beams is proposed. Finally, numerical results are given. In particular, numerical comparisons with the Vlasov–Poisson model illustrate the possibilities of this approach.     

Convergence Conditions for Interpolation Fractions at the Nodes Distinct from the Singular Points of a Function

We consider an interpolation process for the class of functions with finitely many singular points by means of the rational functions whose poles coincide with the singular points of the function under interpolation. The interpolation nodes constitute a triangular matrix and are distinct from the singular points of the function. We find a necessary and sufficient condition for uniform convergence of sequences of interpolation fractions to the function under interpolation on every compact set disjoint from the singular points of the function and other conditions for convergence.Original Russian Text Copyright © 2005 Lipchinskii A. G.__________Translated from Sibirskii Matematicheskii Zhurnal, Vol. 46, No. 4, pp. 822–833, July–August, 2005.     

On Synge-type angle condition for <Emphasis Type="Italic">d</Emphasis>-simplices

The maximum angle condition of J. L. Synge was originally introduced in interpolation theory and further used in finite element analysis and applications for triangular and later also for tetrahedral finite element meshes. In this paper we present some of its generalizations to higher-dimensional simplicial elements. In particular, we prove optimal interpolation properties of linear simplicial elements in ? ^d that degenerate in some way.     

Uniform error estimates for triangular finite element solutions of advection-diffusion equations

In this paper, the authors use the integral identities of triangular linear elements to prove a uniform optimal-order error estimate for the triangular element solution of two-dimensional time-dependent advection-diffusion equations. Also the authors introduce an interpolation postprocessing operator to get the superconvergence estimate under the ε weighted energy norm. The estimates above depend only on the initial and right data but not on the scaling parameter ε.     

Interpolation of analytic functions with finitely many singularities

We consider an interpolation process for the class of functions with finitely many singular points by means of rational functions whose poles coincide with the singular points of the function under interpolation. The interpolation nodes form a triangular matrix. We find necessary and sufficient conditions for the uniform convergence of sequences of interpolation fractions to the function under interpolation on every compact set disjoint from the singular points of the function and other conditions for convergence. We generalize and improve the familiar results on the interpolation of functions with finitely many singular points by rational fractions and of entire functions by polynomials.     

三角形单元上二次Lagrange型插值多项式与被插函数的误差估计

冯康在文［２］中证明了三角形单元Ｃ上一次Ｌａｇｒａｎｇｅ型插值函数Ｕ与被插函数ｕ的误差估计为：这里θ为三角单元Ｃ的最大内角，ｈ为最大边长．本文将该结果推广至二次Ｌａｇｒａｎｇｅ型插值多项式。并得到了相应的误差估计。     

任意三角形单元上的三次Lagrange型插值逼近

本文推广了文［２］中的结果，对于任意三角形单元的三次Ｌａｇｒａｎｇｅ型插值多项式给出了原函数ｕ与被插函数Ｕ之间的误差估计     

多项时间分数阶扩散方程各向异性线性三角元的高精度分析

在各向异性网格下,针对具有Caputo导数的二维多项时间分数阶扩散方程,给出了线性三角形元的高精度分析.首先,基于线性三角形元和改进的L1格式,建立了一个全离散逼近格式,并证明了其无条件稳定性;其次,利用有限元插值算子与Riesz投影算子之间的关系及相关的高精度结果,导出了超逼近性质.进而,借助于插值后处理技术得到了超收敛估计.值得指出的是,单独利用插值算子或Riesz投影都无法得到上述超逼近和超收敛结果.最后,利用数值算例验证了理论分析的正确性.此外,对一些常见的有限单元在该方程的数值逼近方面,作了进一步探讨.     

跳行范德蒙矩阵的三角分解

跳行范德蒙矩阵是一种重要的矩阵,在函数插值等方面有着重要的应用.根据跳行范德蒙矩阵的特殊结构,将跳行范德蒙矩阵分解为一系列下三角矩阵与一系列上三角矩阵的乘积.进一步给出了其逆矩阵分解为一系列上三角矩阵与一系列下三角矩阵的乘积的表达式.     

插值三角多项式在Orlicz空间逼近的渐进等式

对于具有等距分布插值结点的三角多项式,借助广义的Minkowski不等式在Orlicz空间内建立了由三角多项式逼近的渐近等式.并对于Orlicz空间内不同的函数类给出不同的结果.