Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAn Extended Galerkin Analysis for Linear Elasticity with Strongly Symmetric Stress Tensor
96
0
0.0
(
0
)
Added by
Qingguo Hong
Publication date
2020
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
This paper presents an extended Galerkin analysis for various Galerkin methods of the linear elasticity problem. The analysis is based on a unified Galerkin discretization formulation for the linear elasticity problem consisting of four discretization variables: strong symmetric stress tensor $boldsymbol sigma_h$, displacement $u_h$ inside each element and the modifications of these two variables $check boldsymbol sigma_h$ and $check u_h$ on elementary boundaries. Motivated by many relevant methods in literature, this formulation can be used to derive most existing discontinuous, nonconforming and conforming Galerkin methods for linear elasticity problem and especially to develop a number of new discontinuous Galerkin methods. Many special cases of this four-field formulation are proved to be hybridizable and can be reduced to some known hybridizable discontinuous Galerkin, weak Galerkin and local discontinuous Galerkin methods by eliminating one or two of the four fields. As certain stabilization parameter tends to infinity, this four-field formulation is proved to converge to some conforming and nonconforming mixed methods for linear elasticity problem. Two families of inf-sup conditions, one known as $H^1$-philic and another known as H(div)-phillic, are proved to be uniformly valid with respect to different choices of discrete spaces and parameters. These inf-sup conditions guarantee the well-posedness of the new proposed formulations and also offer a new and unified analysis for many existing methods in literature as a by-product.
rate research
Read More
A stress equilibration procedure for linear elasticity is proposed and analyzed in this paper with emphasis on the behavior for (nearly) incompressible materials. Based on the displacement-pressure approximation computed with a stable finite element pair, it constructs an $H (text{div})$-conforming, weakly symmetric stress reconstruction. Our focus is on the Taylor-Hood combination of continuous finite element spaces of polynomial degrees $k+1$ and $k$ for the displacement and the pressure, respectively. Our construction leads then to reconstructed stresses by Raviart-Thomas elements of degree $k$ which are weakly symmetric in the sense that its anti-symmetric part is zero tested against continuous piecewise polynomial functions of degree $k$. The computation is performed locally on a set of vertex patches covering the computational domain in the spirit of equilibration cite{BraSch:08}. Due to the weak symmetry constraint, the local problems need to satisfy consistency conditions associated with all rigid body modes, in contrast to the case of Poissons equation where only the constant modes are involved. The resulting error estimator is shown to constitute a guaranteed upper bound for the error with a constant that depends only on the shape regularity of the triangulation. Local efficiency, uniformly in the incompressible limit, is deduced from the upper bound by the residual error estimator.
A general analysis framework is presented in this paper for many different types of finite element methods (including various discontinuous Galerkin methods). For second order elliptic equation, this framework employs $4$ different discretization variables, $u_h, bm{p}_h, check u_h$ and $check p_h$, where $u_h$ and $bm{p}_h$ are for approximation of $u$ and $bm{p}=-alpha abla u$ inside each element, and $ check u_h$ and $check p_h$ are for approximation of residual of $u$ and $bm{p} cdot bm{n}$ on the boundary of each element. The resulting 4-field discretization is proved to satisfy inf-sup conditions that are uniform with respect to all discretization and penalization parameters. As a result, most existing finite element and discontinuous Galerkin methods can be analyzed using this general framework by making appropriate choices of discretization spaces and penalization parameters.
An interface/boundary-unfitted eXtended hybridizable discontinuous Galerkin (X-HDG) method of arbitrary order is proposed for linear elasticity interface problems on unfitted meshes with respect to the interface and domain boundary. The method uses piecewise polynomials of degrees $k (>= 1)$ and $k-1$ respectively for the displacement and stress approximations in the interior of elements inside the subdomains separated by the interface, and piecewise polynomials of degree $k$ for the numerical traces of the displacement on the inter-element boundaries inside the subdomains and on the interface/boundary of the domain. Optimal error estimates in $L^2$-norm for the stress and displacement are derived. Finally, numerical experiments confirm the theoretical results and show that the method also applies to the case of crack-tip domain.
For the Hodge--Laplace equation in finite element exterior calculus, we introduce several families of discontinuous Galerkin methods in the extended Galerkin framework. For contractible domains, this framework utilizes seven fields and provides a unifying inf-sup analysis with respect to all discretization and penalty parameters. It is shown that the proposed methods can be hybridized as a reduced two-field formulation.
Linear poroelasticity models have a number of important applications in biology and geophysics. In particular, Biots consolidation model is a well-known model that describes the coupled interaction between the linear response of a porous elastic medium and a diffusive fluid flow within it, assuming small deformations. Although deterministic linear poroelasticity models and finite element methods for solving them numerically have been well studied, there is little work to date on robust algorithms for solving poroelasticity models with uncertain inputs and for performing uncertainty quantification (UQ). The Biot model has a number of important physical parameters and inputs whose precise values are often uncertain in real world scenarios. In this work, we introduce and analyse the well-posedness of a new five-field model with uncertain and spatially varying Youngs modulus and hydraulic conductivity field. By working with a properly weighted norm, we establish that the weak solution is stable with respect to variations in key physical parameters, including the Poisson ratio. We then introduce a novel locking-free stochastic Galerkin mixed finite element method that is robust in the incompressible limit. Armed with the `right norm, we construct a parameter-robust preconditioner for the associated discrete systems. Our new method facilitates forward UQ, allowing efficient calculation of statistical quantities of interest and is provably robust with respect to variations in the Poisson ratio, the Biot--Willis constant and the storage coefficient, as well as the discretization parameters.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
