Error estimates for mixed finite element methods for nonlinear parabolic problems

We consider a class of mixed finite element methods for nonlinear parabolic problems over a plane domain. The finite element spaces taken are Raviart-Thomas spaces of index k, k ? 0. We obtain optimal order L²- and almost optimal order L^∞-error estimates for the finite element solution and order optimal L²-error estimates for its gradient. We also derive the error estimates for the time derivatives of the solution. Our results extend those previously obtained by Johnson and Thomée for the corresponding linear problems with k ? 1.     

Solution of Maxwell equation in axisymmetric geometry by Fourier series decompostion and by use of H (rot) conforming finite element

Summary. This study deals with the mathematical and numerical solution of time-harmonic Maxwell equation in axisymmetric geometry. Using Fourier decomposition, we define weighted Sobolev spaces of solution and we prove expected regularity results. A practical contribution of this paper is the construction of a class of finite element conforming with the H (rot) space equipped with the weighted measure rdrdz. It appears as an extension of the well-known cartesian mixed finite element of Raviart-Thomas-Nédélec [11]–[15]. These elements are built from classical lagrangian and mixed finite element, therefore no special approximations functions are needed. Finally, following works of Mercier and Raugel [10], we perform an interpolation error estimate for the simplest proposed element. Received March 15, 1996 / Revised version received November 30, 1998 / Published online December 6, 1999     

SUPERCONVERGENCE OF LEAST-SQUARES MIXED FINITE ELEMENTS FOR ELLIPTIC PROBLEMS ON TRIANGULATION

In this paper, we present the least-squares mixed finite element method and investigate superconvergence phenomena for the second order elliptic boundary-value problems over triangulations. On the basis of the L~2-projection and some mixed finite element projections, we obtain the superconvergence result of least-squares mixed finite     

A family of higher order mixed finite element methods for plane elasticity

Summary The Dirichler problem for the equations of plane elasticity is approximated by a mixed finite element method using a new family of composite finite elements having properties analogous to those possessed by the Raviart-Thomas mixed finite elements for a scalar, second-order elliptic equation. Estimates of optimal order and minimal regularity are derived for the errors in the displacement vector and the stress tensor inL ²(), and optimal order negative norm estimates are obtained inH ^s () for a range ofs depending on the index of the finite element space. An optimal order estimate inL () for the displacement error is given. Also, a quasioptimal estimate is derived in an appropriate space. All estimates are valid uniformly with respect to the compressibility and apply in the incompressible case. The formulation of the elements is presented in detail.This work was performed while Professor Arnold was a NATO Postdoctoral Fellow     

On the problem of spurious eigenvalues in the approximation of linear elliptic problems in mixed form

In the approximation of linear elliptic operators in mixed form, it is well known that the so-called inf-sup and ellipticity in the kernel properties are sufficient (and, in a sense to be made precise, necessary) in order to have good approximation properties and optimal error bounds. One might think, in the spirit of Mercier-Osborn-Rappaz-Raviart and in consideration of the good behavior of commonly used mixed elements (like Raviart-Thomas or Brezzi-Douglas-Marini elements), that these conditions are also sufficient to ensure good convergence properties for eigenvalues. In this paper we show that this is not the case. In particular we present examples of mixed finite element approximations that satisfy the above properties but exhibit spurious eigenvalues. Such bad behavior is proved analytically and demonstrated in numerical experiments. We also present additional assumptions (fulfilled by the commonly used mixed methods already mentioned) which guarantee optimal error bounds for eigenvalue approximations as well.

     

Superconvergence of H 1-Galerkin Mixed Finite Element Methods for Second-Order Elliptic Equations

We derive superconvergence result for H ¹-Galerkin mixed finite element method for second-order elliptic equations over rectangular partitions. Compared to standard mixed finite element procedure, the method is not subject to the Ladyzhenskaya–Bab?ska–Brezzi (LBB) condition and the approximating finite element spaces are allowed to be of different polynomial degrees. Superconvergence estimate of order 𝒪(h ^k+3), where k ≥ 1 is the order of the approximating polynomials employed in the Raviart–Thomas elements, is established for the flux via a postprocessing technique.     

APPLICATION OF SUPERCONVERGENCE TO A MODEL FOR COMPRESSIBLE MISCIBLE DISPLACEMENT

A modification of a finite element method of Douglas and Roberts for approximating the solution of the equations describing compressible miscible displacement in a porous medium is proposed and analyzed. The pressure is treated by a parabolic mixed finite element method using a Raviart-Thomas space of index rover a quasiregular partition, An extension of the Darcy velocity along Gauss lines is used in the evaluation of the coefficients in the Galerkin procedure for the concentration. A simple computational procedure allows the superconvergence property of the fluid velocity to be retained in our total algorithm.     

Mixed finite volume methods for semiconductor device simulation

We deal with the two-dimensional numerical solution of the Van Roosbroeck system, widely employed in modern semiconductor device simulation. Using the well-known Gummel's decoupled algorithm leads to the iterative solution of a nonlinear Poisson equation for the electric potential and two linearized continuity equations for the electron and hole current densities. The numerical approximation is based on the dual mixed formulation for a self-adjoint second-order elliptic operator by using the Raviart-Thomas (RT) finite elements of lowest degree on a triangular partition of the device domain. In this article, we propose a suitable variant of the RT method, based on the diagonalization of the element mass matrix. This is achieved by use of an appropriate numerical integration that eliminates the fluxes and gives rise to a cell-centered finite volume scheme for the scalar unknown with the same approximation properties of the mixed approach, but at a reduced computational cost. The above procedure suggests also a natural way to introduce in the frame of the classical Box Method (BM) suitable vector basis functions (edge elements) to represent the current field over each mesh triangle. This issue may be profitably employed both as a postprocessing tool, as well as a technique for solving the current continuity equations when source terms depending on the current itself are included in the mathematical model. Simulations of realistic semiconductor devices are then included to demonstrate the accuracy and stability of the new method. © 1997 John Wiley & Sons, Inc.     

L ∞-error estimates of triangular mixed finite element methods for optimal control problems governed by semilinear elliptic equations

We look at L ^∞-error estimates for convex quadratic optimal control problems governed by nonlinear elliptic partial differential equations. In so doing, use is made of mixed finite element methods. The state and costate are approximated by the lowest order Raviart-Thomas mixed finite element spaces, and the control, by piecewise constant functions. L ^∞-error estimates of optimal order are derived for a mixed finite element approximation of a semilinear elliptic optimal control problem. Finally, numerical tests are presented which confirm our theoretical results.     

Comparison of a mixed least-squares formulation using different approximation spaces

The main goal of the present work is the comparison of the performance of a least-squares mixed finite element formulation where the solution variables (displacements and stresses) are interpolated using different approximation spaces. Basis for the formulation is a weak form resulting from the minimization of a least-squares functional, compare e.g. [1]. As suitable functions for standard interpolation polynomials of Lagrangian type are chosen. For the conforming discretization of the Sobolev space vector-valued Raviart-Thomas interpolation functions, see also [2], are used. The resulting elements are named as P_mP_k and RT_mP_k. Here m (stresses) and k (displacements) denote the approximation order of the particular interpolation function. For the comparison we consider a two-dimensional cantilever beam under plain strain conditions and small strain assumptions. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Superconvergence and a posteriori error estimates of RT1 mixed methods for elliptic control problems with an integral constraint

In this paper, we investigate the superconvergence property and a posteriori error estimates of mixed finite element methods for a linear elliptic control problem with an integral constraint. The state and co-state are approximated by the order k = 1 Raviart-Thomas mixed finite element spaces and the control variable is approximated by piecewise constant functions. Approximations of the optimal control of the continuous optimal control problem will be constructed by a projection of the discrete adjoint state. It is proved that these approximations have convergence order h ². Moreover, we derive a posteriori error estimates both for the control variable and the state variables. Finally, a numerical example is given to demonstrate the theoretical results.     

Superconvergence of triangular mixed finite elements for optimal control problems with an integral constraint

In this paper, we shall investigate the superconvergence property of quadratic elliptical optimal control problems by triangular mixed finite element methods. The state and co-state are approximated by the order k = 1 Raviart-Thomas mixed finite elements and the control is discretized by piecewise constant functions. We prove the superconvergence error estimate of h² in L²-norm between the approximated solution and the interpolation of the exact control variable. Moreover, by postprocessing technique, we find that the projection of the discrete adjoint state is superclose (in order h²) to the exact control variable.     

鱼骨形网格上二阶方程混合元的超收敛

In this paper, superconvergence of the lowest order Raviart-Thomas mixed finite element approximation for second order Neumann boundary value problem on fishbone shape meshes is analyzed. The main term of the error between the exact solution and the finite element interpolating function is determined by Bramble-Hilbert lemma on the individual finite element. A part of the main term of the error on two adjacent finite elements can be cancelled along the special direction, and thus the higher order error estimate is obtained on the whole domain by summation. Compared with the general finite element error estimate,the convergence rate can be increased from order one to order two in L2-norm by postprocessing superconvergence technique.     

四阶抛物方程H1-Galerkin混合有限元方法的超逼近及最优误差估计

 本文基于双线性元及零阶Raviart-Thomas元 (R-T)对四阶抛物方程建立了半离散和向后欧拉全离散H¹-Galerkin混合有限元格式. 利用积分恒等式技巧和单元的特殊构造, 证明了关于上述两元的两个新的重要性质. 进而导出了这两种格式下相关变量的最优误差估计和超逼近性质.     

Minimization- and Saddle-Point-Based Modeling of Diffusion-Deformation-Processes in Hydrogels

Hydrogels have gained importance during the last years due to their wide range of synthetically fabricable elastic properties as well their increasing meaning in biomedical applications. Future exploitation of the vast prospects of hydrogels is however only feasible by establishing reliable material models that precisely capture their behavior in different environments. To this end, we propose a consistent variational framework for deformation-diffusion processes, offering a canonically compact approach to the chemo-mechanical coupling of hydrogels via a saddle-point as well as a new minimization formulation. The work depicts the construction of rate-type potentials for the chemo-mechanical evolution problem and their transformation into time-discrete incremental potentials. In terms of spatial discretization, the finite element method is employed, benefiting from the intrinsic symmetric structure of the variational foundation. While the saddle-point formulation yields the well-known LBB condition as a constraint for finite element interpolations, on the part of its minimizing counterpart H(Div, ℬ︁)-conforming elements have to be chosen. We illustrate appropriate solutions to both challenges, using mixed Taylor-Hood for the saddle-point and Raviart-Thomas elements for the minimization formulation and discuss advantages of the new approach. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

TRANSCENDENTAL AND ALGEBRAIC EXTENSIONS OF COMMUTATIVE RINGS

Abstract

Transcendental and algebraic elements over commutative rings are defined. Rings with zero nil radical are considered. For a transcendental over R, necessary and sufficient conditions are derived for elements of R[α] to be algebraic or transcendental over R. For R a ring with identity and a finite number of minimal prime ideals, necessary and sufficient conditions are given for any element in a unitary overring of R to be algebraic or transcendental over R. It is proved that if α is algebraic Over R, so is every element of R[α]. It is show that if R is Noetherian, β is algebraic over R[α] and α is algebraic over R, then, under certain conditions, β is algebraic over R. If R has a finite number of minimal prime ideals, P₁,…,P_k, which are pairwise comaximal, then if t is transcendental over R, R[t] can be obtained by adjoining k algebraic elements a_i over R to R whose defining polynomials are in P_i [x], and conversely, if such elements are adjoined to R, they generate an element transcendental over R.     

二阶椭圆方程Dirichlet边值问题混合元的超收敛

关于二阶椭圆方程Dirichlet边值问题混合元的超收敛,在正则矩形网格上,林群和林甲富在文[1]中,采用一阶Raviart-Thomas混合元空间,对有限元解经后处理后,其收敛于精确解的速度从二阶提高到四阶.本文拟将这一结果进行推广,讨论二阶椭圆方程Dirichlet边值问题的k阶Raviart-Thomas混合元的超收敛,得到了以k 3阶速度收敛于精确解的有限元解.     

A posteriori error estimates of mixed DG finite element methods for linear parabolic equations

In this article, we analyse a posteriori error estimates of mixed finite element discretizations for linear parabolic equations. The space discretization is done using the order λ?≥?1 Raviart–Thomas mixed finite elements, whereas the time discretization is based on discontinuous Galerkin (DG) methods (r?≥?1). Using the duality argument, we derive a posteriori l ^∞(L ²) error estimates for the scalar function, assuming that only the underlying mesh is static.     

Error estimates and superconvergence of semidiscrete mixed methods for optimal control problems governed by hyperbolic equations

In this paper, we investigate the L ^??(L ₂)-error estimates and superconvergence of the semidiscrete mixed finite elementmethods for quadratic optimal control problems governed by linear hyperbolic equations. The state and the co-state are discretized by the order k Raviart-Thomas mixed finite element spaces and the control is approximated by piecewise polynomials of order k(k ?? 0). We derive error estimates for approximation of both state and control. Moreover, we present the superconvergence analysis for mixed finite element approximation of the optimal control problems.     

Error estimates of mixed finite element methods for quadratic optimal control problems

In this paper, we investigate the error estimates for the solutions of optimal control problems by mixed finite element methods. The state and costate are approximated by Raviart-Thomas mixed finite element spaces of order k and the control is approximated by piecewise polynomials of order k. Under the special constraint set, we will show that the control variable can be smooth in the whole domain. We derive error estimates of optimal order both for the state variables and the control variable.