Formulation fonction-courant et tourbillon du problème de Stokes dans un domaine bidimensionnel multiplement connexe

We consider the formulation of the Stokes problem in a multiply connected two-dimensional domain where the unknowns are the stream-function and the vorticity. We derive its equivalence with a finite system of several variational problems. This leads to the construction of a finite element discretization of this problem. The analysis and a numerical experiment prove the convergence of the method. To cite this article: M. Amara et al., C. R. Acad. Sci. Paris, Ser. I 342 (2006).     

A parameter-free stabilized finite element method for scalar advection-diffusion problems

We formulate and study numerically a new, parameter-free stabilized finite element method for advection-diffusion problems. Using properties of compatible finite element spaces we establish connection between nodal diffusive fluxes and one-dimensional diffusion equations on the edges of the mesh. To define the stabilized method we extend this relationship to the advection-diffusion case by solving simplified one-dimensional versions of the governing equations on the edges. Then we use H(curl)-conforming edge elements to expand the resulting edge fluxes into an exponentially fitted flux field inside each element. Substitution of the nodal flux by this new flux completes the formulation of the method. Utilization of edge elements to define the numerical flux and the lack of stabilization parameters differentiate our approach from other stabilized methods. Numerical studies with representative advection-diffusion test problems confirm the excellent stability and robustness of the new method. In particular, the results show minimal overshoots and undershoots for both internal and boundary layers on uniform and non-uniform grids.     

An augmented mixed finite element method for 3D linear elasticity problems

In this paper we introduce and analyze a new augmented mixed finite element method for linear elasticity problems in 3D. Our approach is an extension of a technique developed recently for plane elasticity, which is based on the introduction of consistent terms of Galerkin least-squares type. We consider non-homogeneous and homogeneous Dirichlet boundary conditions and prove that the resulting augmented variational formulations lead to strongly coercive bilinear forms. In this way, the associated Galerkin schemes become well posed for arbitrary choices of the corresponding finite element subspaces. In particular, Raviart-Thomas spaces of order 0 for the stress tensor, continuous piecewise linear elements for the displacement, and piecewise constants for the rotation can be utilized. Moreover, we show that in this case the number of unknowns behaves approximately as 9.5 times the number of elements (tetrahedrons) of the triangulation, which is cheaper, by a factor of 3, than the classical PEERS in 3D. Several numerical results illustrating the good performance of the augmented schemes are provided.     

Stable hp mixed finite elements based on the Hellinger–Reissner principle

In the Hellinger–Reissner formulation for linear elasticity, both the displacement u and the stress σ are taken as unknowns, giving rise to a saddle point problem. We present new pairings of quadrilateral ‘trunk’ finite element spaces for this method and prove stability (and optimality) in terms of both h and p. The effect of mesh shape regularity on the stability constant is explicitly tracked. Our results provide a theoretical basis for recent numerical experiments (in the context of a mixed p formulation for viscoelasticity) that showed these spaces worked well computationally.     

Reuse,Recycle, Reduce (3R) – strategies for the calculation of transient magnetic fields

The discretization of transient magneto-dynamic field problems with geometric discretization schemes such as the Finite Integration Technique or the Finite-Element Method based on Whitney form functions results in nonlinear differential-algebraic systems of equations of index 1. Their time integration with embedded s-stage singly diagonal implicit Runge–Kutta methods requires the solution of s nonlinear systems within one time step. Accelerated solution of these schemes is achieved with techniques following so-called 3R-strategies (“reuse, recycle, reduce”). This involves e.g. the solution of the linear(-ized) equations in each time step where the solution process of the iterative preconditioned conjugate gradient method reuses and recycles spectral information of linear systems from previous stages. Additionally, a combination of an error controlled spatial adaptivity and an error controlled implicit Runge–Kutta scheme is used to reduce the number of unknowns for the algebraic problems effectively and to avoid unnecessary fine grid resolutions both in space and time. First numerical results for 2D nonlinear magneto-dynamic problems validate the presented approach and its implementation. The space discretization in the numerical examples is done by Lagrangian nodal finite elements but the presented algorithms also work in combination with other discretization schemes for the Maxwell equations such as the Whitney vector finite elements.     

A finite volume scheme for anisotropic diffusion problems

A new finite volume for the discretization of anisotropic diffusion problems on general unstructured meshes in any space dimension is presented. The convergence of the approximate solution and its discrete gradient is proven. The efficiency of the scheme is illustrated by numerical results. To cite this article: R. Eymard et al., C. R. Acad. Sci. Paris, Ser. I 339 (2004).     

Local flux mimetic finite difference methods

We develop a local flux mimetic finite difference method for second order elliptic equations with full tensor coefficients on polyhedral meshes. To approximate the velocity (vector variable), the method uses two degrees of freedom per element edge in two dimensions and n degrees of freedom per n-gonal mesh face in three dimensions. To approximate the pressure (scalar variable), the method uses one degree of freedom per element. A specially chosen quadrature rule for the L ²-product of vector-functions allows for a local flux elimination and reduction of the method to a cell-centered finite difference scheme for the pressure unknowns. Under certain assumptions, first-order convergence is proved for both variables and second-order convergence is proved for the pressure. The assumptions are verified on simplicial meshes for a particular quadrature rule that leads to a symmetric method. For general polyhedral meshes, non-symmetric methods are constructed based on quadrature rules that are shown to satisfy some of the assumptions. Numerical results confirm the theory.     

A three field stabilized finite element method for the Stokes equationsUne méthode d'éléments finis stabilisée à trois champs pour les équations de Stokes

We consider in this work the boundary value problem for Stokes equations on a two dimensional domain in cases where non-standard boundary conditions are given. We study the cases where pressure and normal or tangential components of the velocity are given in different parts of the boundary and solve the problem with a minimal regularity. We introduce the problem and its variational formulation which is a mixed one. The principal unknowns are the pressure and the vorticity, the multiplier is the velocity. We present the numerical discretization which needs some stabilization. We prove the convergence and the behavior of the a priori error estimates. Some numerical tests are also presented. To cite this article: M. Amara et al., C. R. Acad. Sci. Paris, Ser. I 334 (2002) 603–608.     

A posteriori error estimation in the conforming finite element method based on its local conservativity and using local minimization

We present in this Note fully computable a posteriori error estimates allowing for accurate error control in the conforming finite element discretization of pure diffusion problems. The derived estimates are based on the local conservativity of the conforming finite element method on a dual grid associated with simplex vertices rather than directly on the Galerkin orthogonality. To cite this article: M. Vohralík, C. R. Acad. Sci. Paris, Ser. I 346 (2008).     

Superconvergence analysis of splitting positive definite nonconforming mixed finite element method for pseudo-hyperbolic equations

In this paper, a new splitting positive definite nonconforming mixed finite element method is proposed for pseudo-hyperbolic equations, in which a quasi-Wilson quadrilateral element is used for the flux p, and the bilinear element is used for u. Superconvergence results in ｜｜·｜｜div,h norm for p and optimal error estimates in L2 norm for u are derived for both semi-discrete and fully discrete schemes under almost uniform meshes.     

A mortar element method for hyperbolic problems

A non-conforming finite element method based on non-overlapping domain decomposition is extended to linear hyperbolic problems. The method is based on streamline-diffusion/discontinuous Galerkin methods and the mortar element method. A weak flux continuity condition at the inflow interface is enforced by means of Lagrange multipliers. This weak flux continuity condition replaces the usual mortar condition for elliptic problems, and allows non-matching grids at the subdomain interfaces. To cite this article: Y. Bourgault, A. El Boukili, C. R. Acad. Sci. Paris, Ser. I 338 (2004).     

Admissible solutions of k-extended finite state set and sequence compound decision problems

We treat situations in which independent structurally identical decision problems are to be faced either simultaneously or serially. Recent work on such compound decision problems has centered on finding procedures that satisfy the strengthened asymptotic optimality property of Gilliland and Hannan (Ann. Math. Statist.40 (1969), 1536–1541) and on providing decision rules that at once satisfy an admissibility property and the classical asymptotic optimality property. We suggest a simplified method of accomplishing the first of these goals in general situations and then provide decision rules for finite state components simultaneously admissible and satisfying the strengthened optimality property.     

三维多面体网格上扩散方程的保正格式

 针对三维任意(星形)多面体网格, 本文构造了扩散方程的一种单元中心型非线性有限体积格式, 证明了该格式具有保正性. 在该格式设计中, 除引入网格中心量外, 还引入网格节点量和网格面中心量作为中间未知量, 它们将用网格中心未知量线性组合表示, 使得格式仅有网格中心未知量作为基本未知量. 在节点量计算中, 利用网格面上的调和平均点, 设计了一种适用于三维多面体网格的局部显式加权方法. 该格式适用于求解非平面的网格表面和间断扩散系数的问题. 数值例子验证了它对光滑解具有二阶精度和保正性.     

Error estimates of the lumped mass finite element method for semilinear elliptic problems

We are concerned with the semilinear elliptic problems. We first investigate the L²-error estimate for the lumped mass finite element method. We then use the cascadic multigrid method to solve the corresponding discrete problem. On the basis of the finite element error estimates, we prove the optimality of the proposed multigrid method. We also report some numerical results to support the theory.     

Résolution numérique d'un problème elliptique fortement anisotrope en deux dimensions par une méthode de paramétrisation

We introduce a numerical method for solving an anisotropic elliptic problem. We address the case where the direction of the anisotropy varies, and the anisotropy is high. A finite volume scheme is implemented to solve the problem for small anisotropy ratio, then the parameterization method consists in devising an extrapolation of the solution of the anisotropic problem by combining solutions of a sequence of isotropic problems. To cite this article: P. Guillaume, V. Latocha, C. R. Acad. Sci. Paris, Ser. I 337 (2003).     

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.     

A posteriori error analysis of multipoint flux mixed finite element methods for interface problems

In this paper, the multipoint flux mixed finite element method is used to approximate the flux of two-dimensional elliptic interface problems. Within the class of modified quasi-monotonically distributed coefficients, we derive uniformly robust residual-type a posteriori error estimators for the flux error. Based on the residual-type estimator, we further develop robust implicit and explicit recovery-type estimators through gradient recovery in H(curl) conforming finite element spaces. Numerical experiments are presented to support the theoretical results.     

Adaptive p-version of the finite element method for solving boundary value problems for ordinary second-order differential equations

We suggest an adaptive strategy for constructing a hierarchical basis for a p-version of the finite element method used to solve boundary value problems for second-order ordinary differential equations. The choice of the order of an element on each grid interval is based on estimates of the change, in the norm of C, of the approximate solution or the value of the functional to be minimized when increasing the degree of the basis function added on this interval. The results of numerical experiments estimating the method efficiency are given for sample problems whose solutions have singularities of the boundary layer type. We make a comparison with the p-version of the finite element method, which uses a uniform growth of the degree of the basis functions, and with the h-version, which uses uniform grid refinement along with an adaptive grid refinement and coarsening strategy.     

Gradient WEB-spline finite element method for solving two-dimensional quasilinear elliptic problems

In this paper, a two-dimensional quasilinear elliptic problem of the form -divF(x,▽u)=g(x)$-\mathit{divF}(\mathbf{x}\text{,}▽u)=g(\mathbf{x})$ is considered. This problem is ill-conditioned and we therefore propose a modified iterative algorithm based on coupling of the Sobolev space gradient method and WEB-spline finite element method. Applying the preconditioned iterative method, which has been already provided by Farago and Karatson (2001) [1] reduces the our considered problem to a sequence of linear Poisson’s problems. Then the WEB-spline finite element method is applied to the approximate solution of these Poisson’s problems. In this sense, a convergence theorem is proved and the advantages of this technique than the gradient finite element method (GFEM) is also described. Finally, the presented method is tested on some examples and compared with GFEM. It is shown that the gradient WEB-spline finite element method gives better test results.     

Relationship between multiscale enrichment and stabilized finite element methods for the generalized Stokes problem

We derive a new stabilized finite element method for the generalized Stokes problem starting from the non-stable continuous ${\mathbb{P}}_1/{\mathbb{P}}_1$ finite element space enriched with multiscale functions. The stabilization parameter is related with the enrichment functions which are analytically computed from a boundary value problem at the element level leading to a method which is free of constants. Optimal error estimates are obtained in natural norms and numerical tests validate the method. To cite this article: G.R. Barrenechea, F. Valentin, C. R. Acad. Sci. Paris, Ser. I 341 (2005).