Stabilized mixed finite element methods for linear elasticity on simplicial grids in Rn\mathbb{R}^{n}

In this paper, we design two classes of stabilized mixed finite element methods for linear elasticity on simplicial grids. In the first class of elements, we use $\boldsymbol{H}(\mathbf{div}, \Omega; \mathbb{S})$-$P_k$ and $\boldsymbol{L}^2(\Omega; \mathbb{R}^n)$-$P_{k-1}$ to approximate the stress and displacement spaces, respectively, for $1\leq k\leq n$, and employ a stabilization technique in terms of the jump of the discrete displacement over the faces of the triangulation under consideration; in the second class of elements, we use $\boldsymbol{H}_0^1(\Omega; \mathbb{R}^n)$-$P_{k}$ to approximate the displacement space for $1\leq k\leq n$, and adopt the stabilization technique suggested by Brezzi, Fortin, and Marini. We establish the discrete inf-sup conditions, and consequently present the a priori error analysis for them. The main ingredient for the analysis is two special interpolation operators, which can be constructed using a crucial $\boldsymbol{H}(\mathbf{div})$ bubble function space of polynomials on each element. The feature of these methods is the low number of global degrees of freedom in the lowest order case. We present some numerical results to demonstrate the theoretical estimates.Comment: 16 pages, 1 figur

A priori and a posteriori analysis of non-conforming finite elements with face penalty for advection-diffusion equations

We analyse a non-conforming finite-element method to approximate advection-diffusion-reaction equations. The method is stabilized by penalizing the jumps of the solution and those of its advective derivative across mesh interfaces. The a priori error analysis leads to (quasi-)optimal estimates in the mesh size (sub-optimal by order ½ in the L2-norm and optimal in the broken graph norm for quasi-uniform meshes) keeping the Péclet number fixed. Then, we investigate a residual a posteriori error estimator for the method. The estimator is semi-robust in the sense that it yields lower and upper bounds of the error which differ by a factor equal at most to the square root of the Péclet number. Finally, to illustrate the theory we present numerical results including adaptively generated meshe

Analysis of the discontinuous Galerkin method for elliptic problems on surfaces

We extend the discontinuous Galerkin (DG) framework to a linear second-order elliptic problem on a compact smooth connected and oriented surface. An interior penalty (IP) method is introduced on a discrete surface and we derive a-priori error estimates by relating the latter to the original surface via the lift introduced in Dziuk (1988). The estimates suggest that the geometric error terms arising from the surface discretisation do not affect the overall convergence rate of the IP method when using linear ansatz functions. This is then verified numerically for a number of test problems. An intricate issue is the approximation of the surface conormal required in the IP formulation, choices of which are investigated numerically. Furthermore, we present a generic implementation of test problems on surfaces.Comment: 21 pages, 4 figures. IMA Journal of Numerical Analysis 2013, Link to publication: http://imajna.oxfordjournals.org/cgi/content/abstract/drs033? ijkey=45b23qZl5oJslZQ&keytype=re