اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدA posteriori error estimates by weakly symmetric stress reconstruction for the Biot problem
128
0
0.0
(
0
)
نشر من قبل
Fleurianne Bertrand
تاريخ النشر
2020
مجال البحث
الهندسة المعلوماتية
والبحث باللغة
English
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
A posteriori error estimates are constructed for the three-field variational formulation of the Biot problem involving the displacements, the total pressure and the fluid pressure. The discretization under focus is the H1({Omega})-conforming Taylor-Hood finite element combination, consisting of polynomial degrees k + 1 for the displacements and the fluid pressure and k for the total pressure. An a posteriori error estimator is derived on the basis of H(div)-conforming reconstructions of the stress and flux approximations. The symmetry of the reconstructed stress is allowed to be satisfied only weakly. The reconstructions can be performed locally on a set of vertex patches and lead to a guaranteed upper bound for the error with a constant that depends only on local constants associated with the patches and thus on the shape regularity of the triangulation. Particular emphasis is given to nearly incompressible materials and the error estimates hold uniformly in the incompressible limit. Numerical results on the L-shaped domain confirm the theory and the suitable use of the error estimator in adaptive strategies.
قيم البحث
اقرأ أيضاً
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.
The Richards equation is commonly used to model the flow of water and air through soil, and it serves as a gateway equation for multiphase flows through porous media. It is a nonlinear advection-reaction-diffusion equation that exhibits both paraboli
c-hyperbolic and parabolic-elliptic kinds of degeneracies. In this study, we provide reliable, fully computable, and locally space-time efficient a posteriori error bounds for numerical approximations of the fully degenerate Richards equation. For showing global reliability, a nonlocal-in-time error estimate is derived individually for the time-integrated $H^1(H^{-1})$, $L^2(L^2)$, and the $L^2(H^1)$ errors. A maximum principle and a degeneracy estimator are employed for the last one. Global and local space-time efficiency error bounds are then obtained in a standard $H^1(H^{-1})cap L^2(H^1)$ norm. The reliability and efficiency norms employed coincide when there is no nonlinearity. Moreover, error contributors such as flux nonconformity, time discretization, quadrature, linearization, and data oscillation are identified and separated. The estimates are also valid in a setting where iterative linearization with inexact solvers is considered. Numerical tests are conducted for nondegenerate and degenerate cases having exact solutions, as well as for a realistic case. It is shown that the estimators correctly identify the errors up to a factor of the order of unity.
Hybrid quantum/molecular mechanics models (QM/MM methods) are widely used in material and molecular simulations when MM models do not provide sufficient accuracy but pure QM models are computationally prohibitive. Adaptive QM/MM coupling methods feat
ure on-the-fly classification of atoms during the simulation, allowing the QM and MM subsystems to be updated as needed. In this work, we propose such an adaptive QM/MM method for material defect simulations based on a new residual based it a posteriori error estimator, which provides both lower and upper bounds for the true error. We validate the analysis and illustrate the effectiveness of the new scheme on numerical simulations for material defects.
In two dimensions, we propose and analyze an a posteriori error estimator for finite element approximations of the stationary Navier Stokes equations with singular sources on Lipschitz, but not necessarily convex, polygonal domains. Under a smallness
assumption on the continuous and discrete solutions, we prove that the devised error estimator is reliable and locally efficient. We illustrate the theory with numerical examples.
We present a residual-based a posteriori error estimator for the hybrid high-order (HHO) method for the Stokes model problem. Both the proposed HHO method and error estimator are valid in two and three dimensions and support arbitrary approximation o
rders on fairly general meshes. The upper bound and lower bound of the error estimator are proved, in which proof, the key ingredient is a novel stabilizer employed in the discrete scheme. By using the given estimator, adaptive algorithm of HHO method is designed to solve model problem. Finally, the expected theoretical results are numerically demonstrated on a variety of meshes for model problem.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
