The Least Squares Finite Element Method for Elasticity Interface Problem on Unfitted Mesh

Smoothness and regularity of solutions to PDEs (35B65) Solutions to PDEs in closed form (35C05) PDEs in connection with mechanics of deformable solids (35Q74) Preconditioners for iterative methods (65F08) Iterative numerical methods for linear systems (65F10) Numerical computation of matrix norms, conditioning, scaling (65F35) Numerical optimization and variational techniques (65K10) Stability and convergence of numerical methods for boundary value problems involving PDEs (65N12) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30) Linear elasticity with initial stresses (74B10)

Abstract: In this paper, we propose and analyze the least squares finite element methods for the linear elasticity interface problem in the stress-displacement system on unfitted meshes. We consider the cases that the interface is

C^{2}

or polygonal, and the exact solution

(s i g m a, u)

belongs to

H^{s} (d i v; O m e g a_{0} c u p O m e g a_{1}) i m e s

H^{1+s}(Omega_0 cup Omega_1)

w i t h

s > 1/2

. T w o t y p e s o f l e a s t s q u a r e s f u n c t i o n a l s a r e d e f i n e d t o s e e k t h e n u m e r i c a l s o l u t i o n . T h e f i r s t i s d e f i n e d b y s i m p l y a p p l y i n g t h e

L^2

n o r m l e a s t s q u a r e s p r i n c i p l e, a n d r e q u i r e s t h e c o n d i t i o n

s geq 1

. T h e s e c o n d i s d e f i n e d w i t h a d i s c r e t e m i n u s n o r m, w h i c h i s r e l a t e d t o t h e i n n e r p r o d u c t i n

H^{-1/2}(Gamma)

. T h e u s e o f t h i s d i s c r e t e m i n u s n o r m r e s u l t s i n a m e t h o d o f o p t i m a l c o n v e r g e n c e r a t e s a n d a l l o w s t h e e x a c t s o l u t i o n h a s t h e r e g u l a r i t y o f a n y

s > 1/2

. T h e s t a b i l i t y n e a r t h e i n t e r f a c e f o r b o t h m e t h o d s i s g u a r a n t e e d b y t h e g h o s t p e n a l t y b i l i n e a r f o r m s a n d w e c a n d e r i v e t h e r o b u s t c o n d i t i o n n u m b e r e s t i m a t e s . T h e c o n v e r g e n c e r a t e s u n d e r

L^2$ norm and the energy norm are derived for both methods. We illustrate the accuracy and the robustness of the proposed methods by a series of numerical experiments for test problems in two and three dimensions.

Recommendations

Cited in

(4)

This page was built for publication: The Least Squares Finite Element Method for Elasticity Interface Problem on Unfitted Mesh

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6440347)