A NEW NONCONFORMING MIXED FINITE ELEMENT METHOD FOR LINEAR ELASTICITY

2006 ◽  
Vol 16 (07) ◽  
pp. 979-999 ◽  
Author(s):  
SON-YOUNG YI

We have developed new nonconforming mixed finite element methods for linear elasticity with a pure traction (displacement) boundary condition based on the Hellinger–Reissner variational principle using rectangular elements. Convergence analysis yields an optimal (suboptimal) convergence rate of [Formula: see text] for the L2-error of the stress and [Formula: see text] for the displacement in the pure traction (displacement) boundary problem. However, numerical experiments have yielded optimal-order convergence rates for both stress and displacement in both problems and have shown superconvergence for the displacement at the midpoint of each element. Moreover, we observed that the optimal convergence rates are still valid for large λ.

2020 ◽  
Vol 2020 ◽  
pp. 1-8
Author(s):  
Gwanghyun Jo ◽  
J. H. Kim

In this paper, we suggest a new patch condition for nonconforming mixed finite elements (MFEs) on parallelepiped and provide a framework for the convergence. Also, we introduce a new family of nonconforming MFE space satisfying the new patch condition. The numerical experiments show that the new MFE shows optimal order convergence in Hdiv and L2-norm for various problems with discontinuous coefficient case.


2019 ◽  
Vol 53 (6) ◽  
pp. 2081-2108
Author(s):  
Eldar Khattatov ◽  
Ivan Yotov

Two non-overlapping domain decomposition methods are presented for the mixed finite element formulation of linear elasticity with weakly enforced stress symmetry. The methods utilize either displacement or normal stress Lagrange multiplier to impose interface continuity of normal stress or displacement, respectively. By eliminating the interior subdomain variables, the global problem is reduced to an interface problem, which is then solved by an iterative procedure. The condition number of the resulting algebraic interface problem is analyzed for both methods. A multiscale mortar mixed finite element method for the problem of interest on non-matching multiblock grids is also studied. It uses a coarse scale mortar finite element space on the non-matching interfaces to approximate the trace of the displacement and impose weakly the continuity of normal stress. A priori error analysis is performed. It is shown that, with appropriate choice of the mortar space, optimal convergence on the fine scale is obtained for the stress, displacement, and rotation, as well as some superconvergence for the displacement. Computational results are presented in confirmation of the theory of all proposed methods.


2020 ◽  
Vol 2020 ◽  
pp. 1-11
Author(s):  
Chenguang Zhou ◽  
Yongkui Zou ◽  
Shimin Chai ◽  
Fengshan Zhang

This paper is devoted to the numerical analysis of weak Galerkin mixed finite element method (WGMFEM) for solving a heat equation with random initial condition. To set up the finite element spaces, we choose piecewise continuous polynomial functions of degree j+1 with j≥0 for the primary variables and piecewise discontinuous vector-valued polynomial functions of degree j for the flux ones. We further establish the stability analysis of both semidiscrete and fully discrete WGMFE schemes. In addition, we prove the optimal order convergence estimates in L2 norm for scalar solutions and triple-bar norm for vector solutions and statistical variance-type convergence estimates. Ultimately, we provide a few numerical experiments to illustrate the efficiency of the proposed schemes and theoretical analysis.


2014 ◽  
Vol 2014 ◽  
pp. 1-14 ◽  
Author(s):  
Jinfeng Wang ◽  
Meng Zhao ◽  
Min Zhang ◽  
Yang Liu ◽  
Hong Li

We discuss and analyze anH1-Galerkin mixed finite element (H1-GMFE) method to look for the numerical solution of time fractional telegraph equation. We introduce an auxiliary variable to reduce the original equation into lower-order coupled equations and then formulate anH1-GMFE scheme with two important variables. We discretize the Caputo time fractional derivatives using the finite difference methods and approximate the spatial direction by applying theH1-GMFE method. Based on the discussion on the theoretical error analysis inL2-norm for the scalar unknown and its gradient in one dimensional case, we obtain the optimal order of convergence in space-time direction. Further, we also derive the optimal error results for the scalar unknown inH1-norm. Moreover, we derive and analyze the stability ofH1-GMFE scheme and give the results of a priori error estimates in two- or three-dimensional cases. In order to verify our theoretical analysis, we give some results of numerical calculation by using the Matlab procedure.


Mathematics ◽  
2019 ◽  
Vol 7 (4) ◽  
pp. 325 ◽  
Author(s):  
Muhammad Arshad ◽  
Madiha Sana ◽  
Muhammad Mustahsan

In this paper, the multiblock mortar mixed approximation of second order parabolic partial differential equations is considered. In this method, the simulation domain is decomposed into the non-overlapping subdomains (blocks), and a physically-meaningful boundary condition is set on the mortar interface between the blocks. The governing equations hold locally on each subdomain region. The local problems on blocks are coupled by introducing a special approximation space on the interfaces of neighboring subdomains. Each block is locally covered by an independent grid and the standard mixed finite element method is applied to solve the local problem. The unique solvability of the discrete problem is shown, and optimal order convergence rates are established for the approximate velocity and pressure on the subdomain. Furthermore, an error estimate for the interface pressure in mortar space is presented. The numerical experiments are presented to validate the efficiency of the method.


Sign in / Sign up

Export Citation Format

Share Document