A symmetric mixed finite element method for nearly incompressible elasticity based on. biorthogonal system
|
|
- Harvey Roberts
- 5 years ago
- Views:
Transcription
1 A symmetric mixed finite element metod for nearly incompressible elasticity based on biortogonal systems Bisnu P. Lamicane and Ernst P. Stepan February 18, 2011 Abstract We present a symmetric version of te non-symmetric mixed finite element metod presented in [23] for nearly incompressible elasticity. Te displacement pressure formulation of linear elasticity is discretized using a Petrov Galerkin discretization for te pressure equation in [23] leading to a non-symmetric saddle point problem. A new tree-field formulation is introduced to obtain a symmetric saddle point problem wic allows us to use a biortogonal system. Working wit a biortogonal system, we can statically condense out all auxiliary variables from te saddle point problem arriving at a symmetric and positive-definite system based only on te displacement. We also derive a residual based error estimator for te mixed formulation of te problem. Keywords: Mixed finite elements, symmetric system, Petrov Galerkin discretization, biortogonal system AMS Subject Classification: 65N30, 65N15, 74B10 1 Introduction It is a well-known fact tat linear, bilinear or trilinear finite elements based on te standard displacement based formulation exibit a locking effect wen applied to nearly incompressible elasticity. A popular approac is to use a mixed formulation introducing pressure as an extra unknown. Ten one as to work wit two finite element spaces: one for discretizing te displacement field, and te oter for discretizing te pressure. Tese finite element spaces sould also be compatible in te sense tat a suitable inf- condition is satisfied. One is Scool of Matematical & Pysical Sciences, Matematics Building, University of Newcastle, University Drive, Callagan, NSW 2308, Australia, Bisnu.Lamicane@newcastle.edu.au Leibniz University Hannover, Institute for Applied Matematics, Am Welfengarten 1, Hannover, Germany, stepan@ifam.uni-annover.de 1
2 also interested in statically condensing out te pressure from te system arriving at a formulation based only on te displacement. Te pressure variable can be obtained as a post-processing step. Tere are many different metods to alleviate te locking effect in nearly incompressible elasticity using quadrilateral or exaedral meses. Some popular metods are enanced assumed strain [32, 31, 18, 1], mixed strain metod [19, 20] and strain gap metod [30]. Tese metods are based on a tree-field formulation, popularly known as te Hu-Wasizu formulation, wic is a mixed formulation based on stress, strain and displacement. Using te stability of te mini-element [2], te mixed enanced formulation is extended to simplicial meses in [34], were te pressure variable is assumed to be continuous. Te matematical analysis of te enanced assumed strain for linear elasticity can be found in [29, 7, 26]. We ave sown a uniform convergence of te tree-field formulation in te nearly incompressible elasticity using a modified Hu-Wasizu formulation. However, te matematical analysis in [26] is quite restrictive and only covers a class of quadrilateral meses. Moreover, we could not negate te presence of spurious pressure modes in te pressure or te volumetric part of te stress [7, 26]. One of te simplest mixed formulations for elasticity is a displacement pressure formulation [11, 6]. Tis formulation is obtained by introducing pressure as an extra variable, and writing a variational equation for te pressure. Starting wit tis simplest mixed formulation, a non-symmetric mixed finite element metod based on linear finite elements on simplicial meses is presented in [23] for linear elasticity and in [25] for nonlinear elasticity. Te stability of te formulation is sown by recourse to te mini-element formulation [2] of te Stokes problem. As te variational equation for te pressure is obtained by using a Petrov Galerkin formulation, were te trial and test spaces for te pressure are different, te discrete system is non-symmetric. Te trial and test spaces for te pressure form a biortogonal system, and ence te associated Gram matrix is diagonal. Tis allows an easy static condensation of te pressure from te system leading to a reduced system. As te non-symmetricity is a main drawback of tis approac, in tis paper, we present a symmetric version of tis approac. Te symmetric version is obtained by introducing a Lagrange multiplier variable as in te case of te biarmonic equation [14]. Te continuous formulation is given in Section 2, and te discretization is presented in Section 3. Te numerical analysis of te discrete formulation is presented in Section 4. Te formulation is extended to nonlinear yperelasticity in Section 5. We ave also derived residual based error estimator for our formulation in Section 6. Finally, we ave sown some numerical experiments in te last section. 2
3 2 Te boundary value problem of linear elasticity Tis section is devoted to te introduction of te boundary value problem of linear elasticity. We consider a omogeneous isotropic linear elastic material body occupying a bounded domain in R d, d {2, 3} wit Lipscitz boundary Γ. For a prescribed body force f [L 2 ()] d, te governing equilibrium equation in reads div σ = f, (1) were σ is te symmetric Caucy stress tensor. Te stress tensor σ is defined as a function of te displacement u by te Saint-Venant Kircoff constitutive law σ = 1 2 C( u + [ u]t ), (2) were C is te fourt-order elasticity tensor. Te action of te elasticity tensor C on te strain tensor ε(u) := 1 2 ( u + [ u]t ) is defined as σ = Cε(u) = (tr ε(u))1 + 2µ ε(u). (3) Here, 1 is te identity tensor, and and µ are te Lamé parameters, wic are constant in view of te assumption of a omogeneous body, and wic are assumed positive. Te displacement is assumed to satisfy te omogeneous Diriclet boundary condition u = 0 on Γ. (4) We are interested in te nearly incompressible case, wic corresponds to being very large. Here we use standard notations L 2 (), H 1 () and H 1 0 () for Sobolev spaces, see [8, 15] for details. Let V := [H 1 0 ()] d be te vector Sobolev space wit inner product (, ) 1, and norm 1, defined in te standard way: (u, v) 1, := d i=1 (u i, v i ) 1,, and te norm being induced by tis inner product. We define te bilinear form A(, ) and te linear functional l( ) by A : V V R, A(u, v) := Cε(u) : ε(v) dx, l : V R, l(v) := f v dx. Ten te standard weak form of te linear elasticity problem is as follows: given l V, find u V tat satisfies A(u, v) = l(v), v V, (5) were V is te space of continuous linear functionals on V. Te assumptions on C guarantee tat A(, ) is symmetric, continuous, and V -elliptic. Hence by using standard arguments it can be sown tat (5) as a unique solution u V. 3
4 Furtermore, we assume tat te domain is convex and polygonal or polyedral suc tat u [H 2 ()] d V, and tere exists a constant C independent of suc tat u 2, + div u 1, C f 0,. (6) Te a priori estimate (6) as been sown in [9] for te two-dimensional linear elasticity posed in a convex domain wit polygonal boundary, see [21] for te tree-dimensional case wit convex domain and polyedral boundary. As pointed out in te introduction, te linear elasticity problem can be recast into different mixed formulations. Te easiest mixed formulation is given by introducing pressure as an extra variable, wic leads to penalized Stokes equations. Defining p := div u and R := {q L 2 () : q dx = 0}, a mixed variational formulation of linear elastic problem (5) is given by: find (u, p) V R suc tat were ã(u, v) := 2µ ã(u, v) + b(v, p) = l(v), v V, b(u, q) 1 c(p, q) = 0, q R, (7) ε(u) : ε(v) dx, b(v, q) := div v q dx, and c(p, q) := p q dx. We note tat p R because of te omogeneous Diriclet boundary condition. Te existence, uniqueness and te stability of te mixed formulation (7) can be establised by using te standard saddle point teory [11, 5]. Using a pair of finite element bases forming a biortogonal system to discretize te pressure equation, an efficient numerical sceme arises [23, 25]. Since a Petrov Galerkin discretization is applied to te pressure equation, were te trial and test spaces are different, te arising linear system is non-symmetric even after eliminating te pressure. Terefore, we consider a new formulation, wic allows us to use a biortogonal system and leads to a symmetric formulation. To tis end, we write te standard weak formulation of te linearly elastic problem as a minimization problem: Since min u V 1 A(u, u) l(u). (8) 2 A(u, u) = 2µ ε(u) : ε(u) dx + div u div u dx, we introduce a pressure-like variable φ := div u R, and write te minimization problem (8) as a constrained minimization problem: ( ) 1 min 2µε(u) : ε(u) dx + φ 2 dx l(u). (9) (u,φ) V R φ= 2 div u 4
5 We write a weak variational equation for te equation φ = div u in terms of te Lagrange multiplier space R to obtain te saddle point problem of te minimization problem (9). Te saddle point formulation is to find (u, φ, ξ) V R R suc tat a((u, φ), (v, ψ)) + b((v, ψ), ξ) = l(v), (v, ψ) V R, b((u, φ), η) = 0, η R, (10) were a((u, φ), (v, ψ)) = 2µ b((u, φ), η) = ε(u) : ε(v) dx + (div u φ ) η dx. φψ dx, and Here te Lagrange multiplier plays te role of pressure of te formulation (7), i.e., ξ = φ = div u. Tis idea can be easily generalized to nonlinear elasticity if te problem can be written as a minimization problem (9). In order to sow tat te saddle point problem (10) as a unique solution, we want to apply a standard saddle point teory [11, 6, 8]. To tis end, we need to sow te following tree conditions of well-posedness. 1. Te linear form l( ), te bilinear forms a(, ) and b(, ) are continuous on te spaces on wic tere are defined. 2. Te bilinear form a(, ) is coercive on te space K defined as K = {(v, ψ) V L 2 0() : b((v, ψ), η) = 0, η L 2 0()}. 3. Te bilinear form b(, ) satisfies te inf- condition: b((v, ψ), η) β η 0,, η L 2 (v,ψ) V L 2 0 () v 1, + ψ 0() (11) 0, for a constant β > 0 independent of. It is easy to sow tat te linear form l( ), and te bilinear forms a(, ) and b(, ) are continuous. Te coercivity of te bilinear form a(, ) follows from Korn s inequality: a((u, φ), (u, φ)) = 2µ ε(u) 2 0, + φ 2 0, C( u 2 1, + φ 2 0,). (12) Tus te bilinear form a(, ) is coercive on te wole space V L 2 0(). Now we sow tat te bilinear form b(, ) satisfies te inf- condition uniformly in. 5
6 Lemma 1 Tere exists β > 0 independent of suc tat Proof b((v, ψ), η) β η 0,, η L 2 (v,ψ) V L 2 0 () v 1, + ψ 0(). (13) 0, b((v, ψ), η) (v,ψ) V L 2 0 () v 1, + ψ 0, ( ) = 1 b((v, ψ), η) b((v, ψ), η) +. 2 (v,ψ) V L 2 0 () v 1, + ψ 0, (v,ψ) V L 2 0 () v 1, + ψ 0, Now we substitute ψ = 0 in te remum for te first term, and v = 0 in te remum for te second term to get b((v, ψ), η) 1 (v,ψ) V L 2 0 () v 1, + ψ 0, 2 div v η dx + 1 v V v 1, 2 ψ η. ψ L 2 0 () ψ 0, Te result follows from te fact tat [17] v V div v η dx v 1, β 1 η 0,. Summarizing we ave proved te following teorem. Teorem 2 Te saddle point problem (10) as a unique solution (u, φ, ξ) V L 2 0() L 2 0() and u 1, + φ 0, + ξ 0, f 0,. 3 Finite element discretization We consider a quasi-uniform triangulation T of te polygonal or polyedral domain, were T consists of simplices, eiter triangles or tetraedra. Making use of te standard linear finite element space S defined on te triangulation T S := { v H 1 () : v T P 1(T ), T T } and te space of bubble functions { B := b T P d+1 (T ) : b T T = 0, and T } b T dx > 0, T T, we introduce our finite element space for te displacement as V = (S B ) d V. Te bubble function on an element T is most often defined as b T (x) = c b Π d+1 i=1 T i(x), 6
7 were T i(x) are te barycentric coordinates of te element T associated wit vertices x T i of T, i = 1,..., d + 1. Te constant c b is computed in suc a way tat te value of te bubble function at te barycenter of T is one. Let N be te number of nodes in te finite element mes, and {ϕ 1,..., ϕ N } be te finite element basis of S. Starting wit te basis of S, we construct a dual space M spanned by te basis {ξ 1,..., ξ N } so tat te basis functions of S and M satisfy a condition of biortogonality relation ξ i ϕ j dx = c j δ ij, c j 0, 1 i, j N, (14) were δ ij is te Kronecker symbol, and c j a scaling factor, wic is cosen so tat T ξ i dx = T ϕ i dx. Hence, te sets of basis functions of S and M form a biortogonal system. Te basis functions of M are constructed locally on a reference element ˆT so tat te basis functions of S and M ave te same port [23], and in eac element te sum of all te basis functions of M is one. In te following, we give te local basis functions of M on te reference triangle ˆT := {(x, y) : 0 x, 0 y, x + y 1} and on te reference tetraedron ˆT := {(x, y, z) : 0 x, 0 y, 0 z, x + y + z 1}. For te te reference triangle, we ave ˆξ 1 := 3 4x 4y, ˆξ 2 := 4x 1, and ˆξ 3 := 4y 1, associated wit its tree vertices (0, 0), (1, 0) and (0, 1), respectively, and for te reference tetraedron ˆT := {(x, y, z) : 0 x, 0 y, 0 z, x + y + z 1}, we ave ˆξ 1 := 4 5x 5y 5z, ˆξ 2 := 5x 1, and ˆξ 3 := 5y 1, ˆξ 4 := 5z 1, associated wit its four vertices (0, 0, 0), (1, 0, 0), (0, 1, 0) and (0, 0, 1), respectively. Moreover, d+1 ˆξ i=1 i = 1. Te global basis functions for te space M are constructed by glueing te local basis functions togeter. Tis procedure of constructing global basis functions for M from te local ones is te same as of constructing global basis functions for te standard finite element space S from te local ones. Tese global basis functions ten satisfy te condition of biortogonality (14) wit global finite element basis functions, and pϕ i = pξ i, 1 i N. We can see tat a local basis function of M does not assume value one at one vertex and zero at oter vertices. Since te global basis functions are extended by zero beyond teir port, tey are not continuous. Remark 3 Suc a biortogonal system as been very popular in te context of mortar finite elements [35, 22]. Construction of basis functions of M satisfying te biortogonality and an optimal approximation property for a iger order finite element space is considered in [22]. However, construction of suc a basis for iger order finite element in a simplicial mes is not so easy [27]. 7
8 We also need a subspace of S and a subspace of M aving zero average on defined as { } { } S 0 := v S : v dx = 0, M 0 := ξ M : ξ dx = 0. In [23, 25], te first equation of (7) is discretized by using a Galerkin formulation, and te second equation of (7) is discretized by using a Petrov Galerkin formulation. Te Petrov Galerkin formulation is cosen so tat te pressure solution is taken from S, wereas te test functions are taken from M. Hence te discrete formulation of variational equation (7) is a non-symmetric saddle point problem. Altoug te existence, uniqueness and stability of te solution are establised in [23] using te teory of non-symmetric saddle point problems [28, 4], many well-known linear solvers perform better for symmetric linear systems. Terefore, we consider te discretization of (10) using our biortogonal system. Te discrete formulation is to find (u, φ, ξ ) V S 0 M 0 suc tat a((u, φ ), (v, ψ )) + b((v, ψ ), ξ ) = l(v ), (v, ψ ) V S 0, b((u, φ ), η ) = 0, η M 0. (15) Te goal of coosing te different finite element bases to discretize φ and ξ is to be able to statically condense out tese variables from te system leading to a displacement-based formulation. In order to obtain an algebraic formulation of te discrete saddle point problem (15), we use te same notation for te vector representation of te solution and te solution as elements in V, S and M. Using ψ = 0 and v = 0 subsequently in (15), we ave 2µ ε(u ) : ε(v ) dx + div v ξ dx = l(v ), v V φ ψ ψ ξ dx = 0, ψ S 0. (16) Let A, M, B and D be te matrices associated wit te bilinear forms 2µ ε(u ) : ε(v ) dx, φ ψ dx, div u : η dx, and φ η dx, respectively. Ten te algebraic formulation of te saddle point problem (15) can be written as A 0 B T 0 M 1 D T B 1 D 0 u φ ξ = f 0 0, (17) were f is te vector of discretization of te linear form l( ). Note tat te first two equations in (17) correspond to te two equations in (16), and te last equation in (17) corresponds to te last equation in (15). If we look closely at te linear system (17), we find tat if te matrix D is diagonal, we can easily eliminate te degrees of freedom corresponding to φ and ξ arriving at te formulation involving only one unknown u. After statically condensing out variables φ and ξ, we arrive at te reduced system ( A + B T D 1 MD 1 B ) u = f. As te matrix D 1 is also diagonal, te system matrix is sparse. Te variational formulation of tis reduced system is given in (24). 8
9 4 A priori error estimate Now we turn our attention to sowing te well-posedness and error estimate for our discrete formulation. In order to sow te existence, uniqueness and stability of te solution of (15), we ave to prove tree conditions of well-posedness stated in Section 2 before Lemma 1. As te discretization is conforming te linear form l( ), and te bilinear forms a(, ) and b(, ) are continuous on te discrete spaces as well. An application of Korn s inequality again guarantees te coercivity of te bilinear form a(, ): a((u, φ ), (u, φ )) = 2µ ε(u ) 2 0, + φ 2 0, C( u 2 1, + φ 2 0,). (18) Hence te bilinear form a(, ) is coercive on te wole space V S 0. It remains to sow te tird condition. Tat is, te bilinear form b(, ) satisfies te discrete inf- condition: (v,ψ ) V S 0 b((v, ψ ), η ) v 1, + ψ 0, β η 0,, η M 0 (19) for a constant β > 0 independent of mes-size and. To sow tat te bilinear form b(, ) satisfies te inf- condition (19), we need te following lemma. A proof of tis lemma can be found in [23]. Lemma 4 Tere exists a constant ˆβ independent of te mes-size suc tat div v η dx v V v ˆβ η 0,, η M. 0 (20) 1, Wit te elp of tis lemma, we can sow tat te bilinear form b(, ) satisfies te inf- condition (19). Te proof of tis lemma is similar to te one in te continuous case but we give te proof for completeness. Lemma 5 Tere exists β > 0 independent of mes-size and suc tat (v,ψ ) V S 0 Proof Let η M 0. = 1 2 (v,ψ ) V S 0 ( 1 2 v V (v,ψ ) V S 0 b((v, ψ ), η ) v 1, + ψ 0, β η 0,, η M 0. (21) b((v, ψ ), η ) v 1, + ψ 0, b((v, ψ ), η ) + v 1, + ψ 0, div v η dx v 1, ψ S 0 (v,ψ ) V S 0 ψ η ψ 0,, ) b((v, ψ ), η ) v 1, + ψ 0, 9
10 were te last step is obtained by using ψ = 0 in te remum for te first term, and v = 0 in te remum for te second term from te previous step. Te result follows (wit β = ˆβ 2 ) by using Lemma 4. As all conditions of well-posedness are satisfied, te standard teory of saddle point problem yields te following two teorems [11]: Teorem 6 Te discrete problem (15) as exactly one solution (u, φ, ξ ) V S 0 M 0, wic is uniformly stable wit respect to te data f, and tere exists a constant C independent of mes-size and Lamé parameter suc tat u 1, + φ 0, + ξ 0, C f 0,. Teorem 7 Assume tat (u, φ, ξ) V R R and (u, φ, ξ ) V S 0 M 0 are te solutions of problems (10) and (15), respectively. Ten, te following error estimate olds uniformly in : u u 1, + φ φ 0, + ξ ξ 0, ( C inf u v 1, + inf φ ψ 0, + v V ψ S 0 inf η M 0 ξ η 0, ). It is well-known tat te space V and S 0 ave optimal approximation property. We now establis an approximation property of te space M. For tat purpose, we define a quasi-projection operator Q : L 2 () M by Q v ϕ dx = vϕ dx, v L 2 (), ϕ S. Since Q v M, we ave Q v = N j=1 d jξ j. We multiply bot sides of Q v = N j=1 d jξ j by ϕ i, integrate in and use te biortogonality relation between te basis functions of S and M to get d i = ϕ i v dx. c i Tis allows us to write te action of te operator Q on a function v L 2 () as N Q v = ϕ i v dx ξ i, (22) c i i=1 wic tells tat te operator Q is local in te sense tat for any v L 2 (), te value of Q v at any point in T T only depends on te value of v in S(T ), were S(T ) is te patc of an element T T. Precisely, S(T ) is defined as te interior of te closed set S(T ) = { T T : T T }. (23) 10
11 Te definition of Q allows us to eliminate te auxiliary variables φ and ξ from te discrete saddle point problem (15) arriving at a symmetric and positivedefinite problem of finding u V suc tat 2µ ε(u ) : ε(v ) dx + Q div u Q div v dx = l(v ), v V. (24) Moreover, Q is stable in te L 2 -norm [22]. Lemma 8 For eac element T T, if v L 2 (S(T )), Q v 0,T C v 0,S(T ). (25) Lemma 9 Let v H 1 (). Ten tere exists a constant C > 0 suc tat v Q v 0, C v 1,. (26) Proof We first sow tat if w L 2 () and w = 1 on S(T ), ten Q w T = w on T. (27) Let l 1,, l ns(t ) be te indices of vertices in S(T ) ordered in suc a way tat l 1,, l nt are vertices of element T, and l nt +1,, l ns(t ) are indices of remaining vertices in S(T ). Denoting te port of ϕ i by S i and using te expression of Q w from (22), we ave n T S li ϕ li w dx Q w T = ξ li. (28) c li i=1 Since w = 1 on S(T ), and te sum of all basis functions of M at eac element is one, w = n S(T ) i=1 ξ li on S(T ). We substitute tis expression of w in (28) and obtain Q w = wic concludes tat Q w = w on T. Let S(T ) be te area or te volume of te set S(T ), v H 1 (), and c v := 1 S(T ) v dx. Ten using a triangle S(T ) inequality and Lemma 8, we obtain n T i=1 v Q v 0,T v c v 0,T + Q v c v 0,T = v c v 0,T + Q (v c v ) 0,T C v c v 0,S(T ). An application of te Bramble Hilbert lemma [15] ten yields a constant C > 0 suc tat v Q v 0,T C v 1,S(T ). (29) Te estimate (26) is ten obtained by summing (29) over all elements of T and noting tat eac element is contained in only a fixed number of S(T ) for T T. Ten Teorem 7 yields te linear convergence of te energy norm of te error to zero wit respect to te mes-size if u [H 2 ()] d, φ H 1 () and ξ H 1 (). ξ li, 11
12 5 Extension to finite elasticity In tis section, we briefly outline te extension of tis symmetric formulation to te nonlinear elasticity. Te non-symmetric formulation proposed in [23] is extended to nonlinear yperelastic material law in [24]. Here we propose a symmetric mixed formulation for finite elasticity, were we can apply a biortogonal system for te discretization so tat oter auxiliary variables can be statically condensed out from te system as in te linear case. Te variational approac to te nonlinear elasticity is based on te minimization of te energy functional of te form (W (F (u)) l(u)) dx (30) over a suitable class W of te displacements [16, 13, 33], were F (u) is te deformation gradient defined as F (u) = u + 1. Let J = det(f ). We assume tat te polyconvex energy function W (F ) can be written as W (F ) = H(, J) 2 + G 2 (F ), were G 2 is independent of, and H(, J) is a function of and J. For te isotropic material te energy function W depends only on te tree principal invariants I C, II C and III C of C = F T F, were I C = tr(c), II C = 1 2 (tr2 (C) tr(c 2 )), and III C = det(c) = J 2. If te material law satisfies te two-term Mooney Rivlin law [33], we ave W (F ) = U (J) + µ 2 [(1 c m) (I C 3 2 ln(j)) + c m (II C 3 2 ln(j))] were and µ are Lamé parameters and 0 c m 1 is a material constant. Te real-valued function U is given by U(J) = 1 2 (J 1)2 or U(J) = 1 2 (ln J)2 or U(J) = 1 4 (J ln J). For tis case H(, J) = U(J), and G 2 (F ) = µ 2 [(1 c m) (I C 3 2 ln(j)) + c m (II C 3 2 ln(j))]. We recall tat te standard neo-hookean law is recovered wen c m = 0. Introducing a pressure-like variable φ = H(, J), te mixed formulation for te finite elasticity can be obtained by minimizing te functional (φ 2 + G 2 (F ) l(u)) dx (31) over te function space L 2 () W under te weak constraint (φ H(, J))η dx = 0, η L 2 (). 12
13 Ten te symmetric non-linear saddle point problem is given by te Euler- Lagrange equations of te constrained minimization problem: min (φ 2 + G 2 (F ) l(u)) dx (u,φ) W L 2 () subject to (φ H(, J))η dx = 0, η L 2 (). 6 A posteriori error estimator In tis section, using a similar approac as in [3], we derive an a posteriori error estimator of residual type. Let T be te diameter of element T, and E te size of edge E. We assume sape regularity of te mes, wic, in particular, means T T for all T, T T wit T T =. Here and in te following R S means tat te quantity R is bounded from above by C S, were C is a constant not depending on te mes size. Teorem 10 Let (u, φ, ξ) V R R and (u, φ, ξ ) V S 0 M 0 denote te solution of problems (10) and (15), respectively. Ten tere olds te a posteriori estimate + u u 1, + φ φ 0, + ξ ξ 0, [ T f + div 2µε(u ) 0,T + 1/2 E [2µε(u )n T ] 0,E T T E T φ ξ + ξ 0,T + 0,T div u φ ]. 0,T Proof As in [12], we apply Brezzi s teory [10] to obtain te following inf- condition wit a constant C > 0: Let u 1, + φ 0, + ξ 0, C (v,ψ,η) V R R {a((u, φ), (v, ψ)) + b((u, φ), η) + b((v, ψ), ξ)}. v 1,+ ψ 0,+ η 0,=1 A := a((u u, φ φ ), (v, ψ)) + b((v, ψ), ξ ξ ) I := u u 1, + φ φ 0, + ξ ξ 0, and B := div u φ. 0, Using Caucy-Scwarz inequality, we ave ( b((u u, φ φ ), η) = div (u u ) φ φ ) η dx = b((u, φ ), η) η 0, div u φ 0,, 13
14 and ence Since I C we ave (v,ψ,η) V R R v 1,+ ψ 0,+ η 0,=1 I C b((u u, φ φ ), η) η 0, B. (v,ψ,η) V R R v 1,+ ψ 0,+ η 0,=1 { a((u u, φ φ ), (v, ψ))+ b((u u, φ φ ), η) + b((v, ψ), ξ ξ ) A + b((u u, φ φ ), η). Ten wit te Galerkin ortogonality for arbitrary v A = 2µ ε(u u ) : ε(v) dx + (φ φ )ψ dx + (div v ψ ) (ξ ξ ) dx = 2µ ε(u u ) : ε(v v ) dx + (φ φ )(ψ ψ ) dx ( + div (v v ) ψ ψ ) (ξ ξ ) dx = f (v v ) dx 2µ ε(u ) : ε(v v ) dx φ (ψ ψ ) dx ( div (v v ) ψ ψ ) ξ dx = [ (f + div 2µε(u )) (v v ) dx 2µ(ε(u )n T ) (v v ) dσ T T T T ( φ ξ ) ] (ψ ψ ) dx ξ div (v v ) dx, T T were n T is te unit outward normal vector on T, and ε(u )n T is interpreted as matrix-vector product, wic produces a vector, and ence (ε(u )n T ) (v v ) is a scalar function. Tus using Caucy-Scwarz inequality and trace teorem and coosing v an interpolant to v, we obtain A C [ T f + div 2µε(u ) 0,T v 1,T + 1/2 E [2µε(u )n T ] 0,E v 1, T + T T E T φ ξ ] ψ ψ 0,T + ξ 0,T div (v v ) 0,T 0,T C [ T f + div 2µε(u ) 0,T + 1/2 E [2µε(u )n T ] 0,E T T + ξ 0,T + φ ξ ], 0,T E T }, 14
15 were [u] represents te jump of te function u across te interelement boundaries, and C now depends on v and ψ. We ave also used te fact tat v v 0,T C v 1,T. Let [ D = T f + div 2µε(u ) 0,T + 1/2 E [2µε(u )n T ] 0,E T T E T + ξ 0,T + φ ξ ]. 0,T Tis sows tat I B + D. Tus te total error I is bounded by computable error B + D. In te following, we sow tat te part B of te computable error B + D converges asymptotically to zero wen 0. To sow tis, we introduce anoter quasi-projection operator Q : L2 () S, Q v η dx = v η dx, η M. Tis operator is uniquely defined due to te biortogality relation between S and M. Since Q v S and Q w M, for v, w L 2 (), we ave Q v Q w dx = w Q v dx = v Q w dx. Moreover, we ave Q v = v for v S, and ence Q is an projection onto S. Terefore, Q is stable in te L2 -norm and as te following approximation property in te L 2 -norm: v Q v 0, v 1,, v H 1 (). Tat means Lemmas 8 and 9 can be reproduced for Q. Ten using φ = Q div u, we ave for 0 < r 1, B = div u φ = div u div u + div u Q 0, div u 0, C r div u r,, were we ave assumed tat te exact solution u H 1+r () wit r > 0. Tis sows tat B converges to zero asymptotically wen 0. Te displacement solution is computed by inverting te reduced system after statically condensing out φ and ξ. Ten te solutions φ and ξ can be efficiently computed just by inverting a diagonal matrix. Tis gives an efficient implementation of te error estimator. 15
16 7 Numerical examples In tis section, we demonstrate te performance of our new formulation for linear elasticity wit some numerical examples. Te material parameters are given in compatible units. Bot following examples do not ave pure Diriclet boundary condition, and terefore, we do not need to satisfy te zero mean condition for functions in M and S. However, in case of pure Diriclet boundary condition, bot trial and test functions in M and S sould satisfy te zero mean condition. Tat means te solutions ξ M and φ S are determined up to an additive constant. We ten searc for solutions ξ M and φ S wose values at one arbitrary vertex are prescribed a priori. Example 1: Cook s membrane Our first problem is te popular bencmark problem called Cook s membrane problem. Te Cook s membrane is a twodimensional tapered panel := conv{(0, 0), (48, 44), (48, 60), (0, 44)}, were conv{ξ} represents te convex ull of te set ξ. Te left boundary of te panel is clamped in bot directions and te rigt boundary is subjected to an inplane sear load in te positive y- direction as sown in te left picture of Figure 1. We consider te material parameters E = 250 and ν = We compute te vertical tip displacement at te top-rigt corner of te membrane using uniform refinement of te initial mes given in te rigt picture of Figure 1. In te rigt picture of Figure 1, we sow te convergence of te numerical results wit respect to te number of elements using tree finite element formulations for linear elasticity. We can see te good convergence beaviour of te two mixed formulations, wereas te standard displacement formulation sows te locking effect. Te new mixed formulation sows even better accuracy tan te mini finite element in te coarse mes. Te mini finite element proposed in [2] is applied to te linear elasticity. Example 2: Rectangular beam In tis second example, we consider a linear elastic beam of rectangular size subjected to a couple at one end, as sown in Figure 2. Along te edge x = 0, te orizontal displacement and vertical surface traction are zero. At te point (0, 0), te vertical displacement is also zero. Te exact solution is given by u(x, y) = 2f(1 ( ) ν2 ) l x El 2 y, and v(x, y) = f(1 [ ν2 ) x 2 + ν ] y(y l). El 1 ν We set L = 10, l = 2, E = 1500, ν = , and f = We ave sown te setting of te problem in Figure 2, and te discretization errors wit respect to te number of elements are presented in Figure 3. As can be seen from Figure 3, te standard approac locks completely, wereas we get very good numerical approximations wit te new approac and mini-element. 16
17 T 100N vertical tip displacement at T standard mini new mixed number of elements per edge Figure 1: Problem setting and initial mes (left) and vertical tip displacement at T versus number of elements using different formulations (rigt) l f L Figure 2: Te rectangular beam wit initial mes and problem setting 10 1 L 2 error 10 2 H 1 error standard mini new mixed O( 2 ) 10 2 standard mini new mixed O() number of elements number of elements Figure 3: Error plot versus number of elements, L 2 -norm (left) and H 1 -norm (rigt), rectangular beam Acknowledgement We are grateful to te anonymous referees for teir valuable suggestions to improve te quality of te earlier version of tis work. 17
18 References [1] F. Armero. On te locking and stability of finite elements in finite deformation plane strain problems. Computers and Structures, 75: , [2] D. N. Arnold, F. Brezzi, and M. Fortin. A stable finite element for te stokes equations. Calcolo, 21: , [3] I. Babuška and T. Strouboulis. Te Finite Element Metod and its Reliability. Oxford University Press, New York, [4] C. Bernardi, C. Canuto, and Y. Maday. Generalized Inf-Sup conditions for Cebysev spectral approximation of te Stokes problem. SIAM Journal on Numerical Analysis, 25(6): , [5] D. Braess. Stability of saddle point problems wit penalty. M 2 AN, 30: , [6] D. Braess. Finite Elements. Teory, fast solver, and applications in solid mecanics. Cambridge Univ. Press, Second Edition, [7] D. Braess, C. Carstensen, and B.D. Reddy. Uniform convergence and a posteriori estimators for te enanced strain finite element metod. Numerisce Matematik, 96(3): , [8] S.C. Brenner and L.R. Scott. Te Matematical Teory of Finite Element Metods. Springer Verlag, New York, [9] S.C. Brenner and L. Sung. Linear finite element metods for planar linear elasticity. Matematics of Computation, 59: , [10] F. Brezzi. On te existence, uniqueness and approximations of saddlepoint problems arising from Lagrange multipliers. RAIRO Anal. Numér., 8: , [11] F. Brezzi and M. Fortin. Mixed and ybrid finite element metods. Springer Verlag, New York, [12] U. Brink and E.P. Stepan. Adaptive coupling of boundary elements and mixed finite elements for incompressible elasticity. Numerical Metods for Partial Differential Equations, 17:79 92, [13] C. Carstensen and G. Dolzmann. An a priori error estimate for finite element discretizations in nonlinear elasticity for polyconvex materials under small loads. Numerisce Matematik, 97:67 80, [14] P. Ciarlet and P. Raviart. A mixed finite element metod for te biarmonic equation. In C. De Boor, editor, Symposium on Matematical Aspects of Finite Elements in Partial Differential Equations, pages , New York, Academic Press. 18
19 [15] P.G Ciarlet. Te Finite Element Metod for Elliptic Problems. Nort Holland, Amsterdam, [16] P.G. Ciarlet. Matematical Elasticity Volume I: Tree-Dimensional Elasticity. Nort-Holland, Amsterdam, [17] V. Girault and P.-A. Raviart. Finite Element Metods for Navier-Stokes Equations. Springer-Verlag, Berlin, [18] S. Glaser and F. Armero. On te formulation of enanced strain finite elements in finite deformation. Engineering Compuations, 14: , [19] E.P. Kasper and R.L. Taylor. A mixed-enanced strain metod. Part I: geometrically linear problems. Computers and Structures, 75: , [20] E.P. Kasper and R.L. Taylor. A mixed-enanced strain metod. Part II: geometrically nonlinear problems. Computers and Structures, 75: , [21] V.A. Kozlov, V.G. Maz ya, and J. Rossmann. Spectral Problems Associated wit Corner Singularities of Solutions to Elliptic Equations. Matematical Surveys and Monograps 85. American Matematical Society, Providence, RI, [22] B.P. Lamicane. Higer Order Mortar Finite Elements wit Dual Lagrange Multiplier Spaces and Applications. PD tesis, Universität Stuttgart, [23] B.P. Lamicane. A mixed finite element metod based on a biortogonal system for nearly incompressible elastic problems. In Geoffry N. Mercer and A. J. Roberts, editors, Proceedings of te 14t Biennial Computational Tecniques and Applications Conference, CTAC-2008, volume 50 of ANZIAM J., pages C324 C338, [24] B.P. Lamicane. A mixed finite element metod for nonlinear and nearly incompressible elasticity based on biortogonal systems. International Journal for Numerical Metods in Engineering, 79: , [25] B.P. Lamicane. Mortar finite elements for coupling compressible and nearly incompressible materials in elasticity. International Journal of Numerical Analysis and Modeling, 6: , [26] B.P. Lamicane, B.D. Reddy, and B.I. Wolmut. Convergence in te incompressible limit of finite element approximations based on te Hu- Wasizu formulation. Numerisce Matematik, 104: , [27] B.P. Lamicane, R.P. Stevenson, and B.I. Wolmut. Higer order mortar finite element metods in 3D wit dual Lagrange multiplier bases. Numerisce Matematik, 102:93 121,
20 [28] R.A. Nicolaides. Existence, uniqueness and approximation for generalized saddle point problems. SIAM Journal on Numerical Analysis, 19(2): , [29] B.D. Reddy and J.C. Simo. Stability and convergence of a class of enanced strain metods. SIAM Journal on Numerical Analysis, 32: , [30] G. Romano, F. Marrotti de Sciarra, and M. Diaco. Well-posedness and numerical performances of te strain gap metod. International Journal for Numerical Metods in Engineering, 51: , [31] J.C. Simo and F. Armero. Geometrically nonlinear enanced strain mixed metods and te metod of incompatible modes. International Journal for Numerical Metods in Engineering, 33: , [32] J.C. Simo and M.S. Rifai. A class of assumed strain metod and te metods of incompatible modes. International Journal for Numerical Metods in Engineering, 29: , [33] E. Stein and M. Rüter. Finite element metods for elasticity wit errorcontrolled discretization and model adaptivity. In E. Stein, R. de Borst, and T.J.R. Huges, editors, Encyclopedia of Computational Mecanics, pages Wiley, [34] R.L. Taylor. A mixed-enanced formulation for tetraedral finite elements. International Journal for Numerical Metods in Engineering, 47: , [35] B.I. Wolmut. A mortar finite element metod using dual spaces for te Lagrange multiplier. SIAM Journal on Numerical Analysis, 38: ,
Key words. Sixth order problem, higher order partial differential equations, biharmonic problem, mixed finite elements, error estimates.
A MIXED FINITE ELEMENT METHOD FOR A SIXTH ORDER ELLIPTIC PROBLEM JÉRÔME DRONIOU, MUHAMMAD ILYAS, BISHNU P. LAMICHHANE, AND GLEN E. WHEELER Abstract. We consider a saddle point formulation for a sixt order
More informationc 2004 Society for Industrial and Applied Mathematics
SIAM J NUMER ANAL Vol 4, No, pp 86 84 c 004 Society for Industrial and Applied Matematics LEAST-SQUARES METHODS FOR LINEAR ELASTICITY ZHIQIANG CAI AND GERHARD STARKE Abstract Tis paper develops least-squares
More informationarxiv: v1 [math.na] 28 Apr 2017
THE SCOTT-VOGELIUS FINITE ELEMENTS REVISITED JOHNNY GUZMÁN AND L RIDGWAY SCOTT arxiv:170500020v1 [matna] 28 Apr 2017 Abstract We prove tat te Scott-Vogelius finite elements are inf-sup stable on sape-regular
More informationNumerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems
Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for
More information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 37, pp. 166-172, 2010. Copyright 2010,. ISSN 1068-9613. ETNA A GRADIENT RECOVERY OPERATOR BASED ON AN OBLIQUE PROJECTION BISHNU P. LAMICHHANE Abstract.
More informationMIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp
MIXED DISCONTINUOUS GALERIN APPROXIMATION OF THE MAXWELL OPERATOR PAUL HOUSTON, ILARIA PERUGIA, AND DOMINI SCHÖTZAU SIAM J. Numer. Anal., Vol. 4 (004), pp. 434 459 Abstract. We introduce and analyze a
More informationarxiv: v1 [math.na] 20 Jul 2009
STABILITY OF LAGRANGE ELEMENTS FOR THE MIXED LAPLACIAN DOUGLAS N. ARNOLD AND MARIE E. ROGNES arxiv:0907.3438v1 [mat.na] 20 Jul 2009 Abstract. Te stability properties of simple element coices for te mixed
More informationDifferent Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method
WDS'10 Proceedings of Contributed Papers, Part I, 151 156, 2010. ISBN 978-80-7378-139-2 MATFYZPRESS Different Approaces to a Posteriori Error Analysis of te Discontinuous Galerkin Metod I. Šebestová Carles
More informationarxiv: v1 [math.na] 27 Jan 2014
L 2 -ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS OF BOUNDARY FLUXES MATS G. LARSON AND ANDRÉ MASSING arxiv:1401.6994v1 [mat.na] 27 Jan 2014 Abstract. We prove quasi-optimal a priori error estimates
More informationOptimization of stress modes by energy compatibility for 4-node hybrid quadrilaterals
INTERNATIONAL JOURNAL FOR NUMERIAL METHODS IN ENGINEERING Int. J. Numer. Met. Engng 2004; 59:293 33 (DOI: 0.002/nme.877) Optimization of stress modes by energy compatibility for 4-node ybrid quadrilaterals
More informationA mixed finite element method for nonlinear and nearly incompressible elasticity based on biorthogonal systems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2000; 00:1 6 [Version: 2002/09/18 v2.02] A mixed finite element method for nonlinear and nearly incompressible elasticity
More informationarxiv: v1 [math.na] 17 Jul 2014
Div First-Order System LL* FOSLL* for Second-Order Elliptic Partial Differential Equations Ziqiang Cai Rob Falgout Sun Zang arxiv:1407.4558v1 [mat.na] 17 Jul 2014 February 13, 2018 Abstract. Te first-order
More informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
More informationFinite Element Methods for Linear Elasticity
Finite Element Metods for Linear Elasticity Ricard S. Falk Department of Matematics - Hill Center Rutgers, Te State University of New Jersey 110 Frelinguysen Rd., Piscataway, NJ 08854-8019 falk@mat.rutgers.edu
More informationarxiv: v1 [math.na] 11 May 2018
Nitsce s metod for unilateral contact problems arxiv:1805.04283v1 [mat.na] 11 May 2018 Tom Gustafsson, Rolf Stenberg and Jua Videman May 14, 2018 Abstract We derive optimal a priori and a posteriori error
More informationA Two-field Finite Element Solver for Poroelasticity on Quadrilateral Meshes
A Two-field Finite Element Solver for Poroelasticity on Quadrilateral Meses Graam Harper 1, Jiangguo Liu 1, Simon Tavener 1, and Zuoran Wang 1 Colorado State University, Fort Collins, CO 80523, USA {arper,liu,tavener,wangz}@mat.colostate.edu
More informationA New Class of Zienkiewicz-Type Nonconforming Element in Any Dimensions
Numerisce Matematik manuscript No. will be inserted by te editor A New Class of Zienkiewicz-Type Nonconforming Element in Any Dimensions Wang Ming 1, Zong-ci Si 2, Jincao Xu1,3 1 LMAM, Scool of Matematical
More informationA Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Hybrid Mixed Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationAnalysis of A Continuous Finite Element Method for H(curl, div)-elliptic Interface Problem
Analysis of A Continuous inite Element Metod for Hcurl, div)-elliptic Interface Problem Huoyuan Duan, Ping Lin, and Roger C. E. Tan Abstract In tis paper, we develop a continuous finite element metod for
More informationA UNIFORM INF SUP CONDITION WITH APPLICATIONS TO PRECONDITIONING
A UNIFORM INF SUP CONDIION WIH APPLICAIONS O PRECONDIIONING KEN ANDRE MARDAL, JOACHIM SCHÖBERL, AND RAGNAR WINHER Abstract. A uniform inf sup condition related to a parameter dependent Stokes problem is
More informationA SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION
A SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION GABRIEL R. BARRENECHEA, LEOPOLDO P. FRANCA 1 2, AND FRÉDÉRIC VALENTIN Abstract. Tis work introduces and analyzes novel stable
More informationStress-based finite element methods in linear and nonlinear solid mechanics
Stress-based finite element metods in linear and nonlinear solid mecanics Benjamin Müller and Gerard Starke Fakultät für Matematik, Universität Duisburg-Essen, 45127 Essen, Germany Abstract A comparison
More informationA Mixed-Hybrid-Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Mixed-Hybrid-Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationarxiv: v1 [math.na] 12 Mar 2018
ON PRESSURE ESTIMATES FOR THE NAVIER-STOKES EQUATIONS J A FIORDILINO arxiv:180304366v1 [matna 12 Mar 2018 Abstract Tis paper presents a simple, general tecnique to prove finite element metod (FEM) pressure
More informationarxiv: v1 [math.na] 9 Sep 2015
arxiv:509.02595v [mat.na] 9 Sep 205 An Expandable Local and Parallel Two-Grid Finite Element Sceme Yanren ou, GuangZi Du Abstract An expandable local and parallel two-grid finite element sceme based on
More informationMixed Finite Element Methods for Incompressible Flow: Stationary Stokes Equations
Mixed Finite Element Metods for Incompressible Flow: Stationary Stoes Equations Ziqiang Cai, Carles Tong, 2 Panayot S. Vassilevsi, 2 Cunbo Wang Department of Matematics, Purdue University, West Lafayette,
More informationAn approximation method using approximate approximations
Applicable Analysis: An International Journal Vol. 00, No. 00, September 2005, 1 13 An approximation metod using approximate approximations FRANK MÜLLER and WERNER VARNHORN, University of Kassel, Germany,
More informationPart VIII, Chapter 39. Fluctuation-based stabilization Model problem
Part VIII, Capter 39 Fluctuation-based stabilization Tis capter presents a unified analysis of recent stabilization tecniques for te standard Galerkin approximation of first-order PDEs using H 1 - conforming
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationMIXED AND STABILIZED FINITE ELEMENT METHODS FOR THE OBSTACLE PROBLEM
MIXED AND STABILIZED FINITE ELEMENT METHODS FOR THE OBSTACLE PROBLEM TOM GUSTAFSSON, ROLF STENBERG, AND JUHA VIDEMAN Abstract. We discretize te Lagrange multiplier formulation of te obstacle problem by
More informationPosteriori Analysis of a Finite Element Discretization for a Penalized Naghdi Shell
International Journal of Difference Equations ISSN 0973-6069, Volume 8, Number 1, pp. 111 124 (2013) ttp://campus.mst.edu/ijde Posteriori Analysis of a Finite Element Discretization for a Penalized Nagdi
More informationInf sup testing of upwind methods
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Met. Engng 000; 48:745 760 Inf sup testing of upwind metods Klaus-Jurgen Bate 1; ;, Dena Hendriana 1, Franco Brezzi and Giancarlo
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
More informationAdaptive Finite Element Method
39 Capter 3 Adaptive inite Element Metod 31 Introduction As already pointed out in capter 2, singularities occur in interface problems Wen discretizing te problem (221) wit inite Elements, te singularities
More informationMass Lumping for Constant Density Acoustics
Lumping 1 Mass Lumping for Constant Density Acoustics William W. Symes ABSTRACT Mass lumping provides an avenue for efficient time-stepping of time-dependent problems wit conforming finite element spatial
More informationKey words. Finite element method; convection-diffusion-reaction; nonnegativity; boundedness
PRESERVING NONNEGATIVITY OF AN AFFINE FINITE ELEMENT APPROXIMATION FOR A CONVECTION-DIFFUSION-REACTION PROBLEM JAVIER RUIZ-RAMÍREZ Abstract. An affine finite element sceme approximation of a time dependent
More informationVolume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households
Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of
More informationA Weak Galerkin Method with an Over-Relaxed Stabilization for Low Regularity Elliptic Problems
J Sci Comput (07 7:95 8 DOI 0.007/s095-06-096-4 A Weak Galerkin Metod wit an Over-Relaxed Stabilization for Low Regularity Elliptic Problems Lunji Song, Kaifang Liu San Zao Received: April 06 / Revised:
More informationA SADDLE POINT LEAST SQUARES APPROACH TO MIXED METHODS
A SADDLE POINT LEAST SQUARES APPROACH TO MIXED METHODS CONSTANTIN BACUTA AND KLAJDI QIRKO Abstract. We investigate new PDE discretization approaces for solving variational formulations wit different types
More informationVariational Localizations of the Dual Weighted Residual Estimator
Publised in Journal for Computational and Applied Matematics, pp. 192-208, 2015 Variational Localizations of te Dual Weigted Residual Estimator Tomas Ricter Tomas Wick Te dual weigted residual metod (DWR)
More informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationEXISTENCE OF SOLUTIONS FOR A CLASS OF VARIATIONAL INEQUALITIES
Journal of Matematics and Statistics 9 (4: 35-39, 3 ISSN: 549-3644 3 doi:.3844/jmssp.3.35.39 Publised Online 9 (4 3 (ttp://www.tescipub.com/jmss.toc EXISTENCE OF SOLUTIONS FOR A CLASS OF ARIATIONAL INEQUALITIES
More informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
More informationPriority Program 1253
Deutsce Forscungsgemeinscaft Priority Program 53 Optimization wit Partial Differential Equations K. Deckelnick, A. Günter and M. Hinze Finite element approximation of elliptic control problems wit constraints
More informationAN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS
AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS BERNARDO COCKBURN, JOHNNY GUZMÁN, SEE-CHEW SOON, AND HENRYK K. STOLARSKI Abstract. Te embedded discontinuous
More informationarxiv: v1 [math.na] 19 Mar 2018
A primal discontinuous Galerkin metod wit static condensation on very general meses arxiv:1803.06846v1 [mat.na] 19 Mar 018 Alexei Lozinski Laboratoire de Matématiques de Besançon, CNRS UMR 663, Univ. Bourgogne
More informationAPPROXIMATION BY QUADRILATERAL FINITE ELEMENTS
MATHEMATICS OF COMPUTATION Volume 71, Number 239, Pages 909 922 S 0025-5718(02)01439-4 Article electronically publised on Marc 22, 2002 APPROXIMATION BY QUADRILATERAL FINITE ELEMENTS DOUGLAS N. ARNOLD,
More informationFINITE ELEMENT APPROXIMATIONS AND THE DIRICHLET PROBLEM FOR SURFACES OF PRESCRIBED MEAN CURVATURE
FINITE ELEMENT APPROXIMATIONS AND THE DIRICHLET PROBLEM FOR SURFACES OF PRESCRIBED MEAN CURVATURE GERHARD DZIUK AND JOHN E. HUTCHINSON Abstract. We give a finite element procedure for te Diriclet Problem
More informationSome Error Estimates for the Finite Volume Element Method for a Parabolic Problem
Computational Metods in Applied Matematics Vol. 13 (213), No. 3, pp. 251 279 c 213 Institute of Matematics, NAS of Belarus Doi: 1.1515/cmam-212-6 Some Error Estimates for te Finite Volume Element Metod
More informationAdvancements In Finite Element Methods For Newtonian And Non-Newtonian Flows
Clemson University TigerPrints All Dissertations Dissertations 8-3 Advancements In Finite Element Metods For Newtonian And Non-Newtonian Flows Keit Galvin Clemson University, kjgalvi@clemson.edu Follow
More informationarxiv: v2 [math.na] 5 Jul 2017
Trace Finite Element Metods for PDEs on Surfaces Maxim A. Olsanskii and Arnold Reusken arxiv:1612.00054v2 [mat.na] 5 Jul 2017 Abstract In tis paper we consider a class of unfitted finite element metods
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationEfficient, unconditionally stable, and optimally accurate FE algorithms for approximate deconvolution models of fluid flow
Efficient, unconditionally stable, and optimally accurate FE algoritms for approximate deconvolution models of fluid flow Leo G. Rebolz Abstract Tis paper addresses an open question of ow to devise numerical
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationFEM solution of the ψ-ω equations with explicit viscous diffusion 1
FEM solution of te ψ-ω equations wit explicit viscous diffusion J.-L. Guermond and L. Quartapelle 3 Abstract. Tis paper describes a variational formulation for solving te D time-dependent incompressible
More informationFINITE ELEMENT EXTERIOR CALCULUS FOR PARABOLIC EVOLUTION PROBLEMS ON RIEMANNIAN HYPERSURFACES
FINITE ELEMENT EXTERIOR CALCULUS FOR PARABOLIC EVOLUTION PROBLEMS ON RIEMANNIAN HYPERSURFACES MICHAEL HOLST AND CHRIS TIEE ABSTRACT. Over te last ten years, te Finite Element Exterior Calculus (FEEC) as
More informationA Two-Field Finite Element Solver for Poroelasticity on Quadrilateral Meshes
A Two-Field Finite Element Solver for Poroelasticity on Quadrilateral Meses Graam Harper, Jiangguo Liu, Simon Tavener, and Zuoran Wang (B) Colorado State University, Fort Collins, CO 80523, USA {arper,liu,tavener,wangz}@mat.colostate.edu
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationChapter 10. Function approximation Function approximation. The Lebesgue space L 2 (I)
Capter 1 Function approximation We ave studied metods for computing solutions to algebraic equations in te form of real numbers or finite dimensional vectors of real numbers. In contrast, solutions to
More informationCrouzeix-Velte Decompositions and the Stokes Problem
Crouzeix-Velte Decompositions and te Stokes Problem PD Tesis Strauber Györgyi Eötvös Loránd University of Sciences, Insitute of Matematics, Matematical Doctoral Scool Director of te Doctoral Scool: Dr.
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Matematik Optimal L 2 velocity error estimates for a modified pressure-robust Crouzeix-Raviart Stokes element A. Linke, C. Merdon und W. Wollner Tis preprint is also
More informationNUMERICAL DIFFERENTIATION
NUMERICAL IFFERENTIATION FIRST ERIVATIVES Te simplest difference formulas are based on using a straigt line to interpolate te given data; tey use two data pints to estimate te derivative. We assume tat
More informationTHE INF-SUP STABILITY OF THE LOWEST ORDER TAYLOR-HOOD PAIR ON ANISOTROPIC MESHES arxiv: v1 [math.na] 21 Oct 2017
THE INF-SUP STABILITY OF THE LOWEST ORDER TAYLOR-HOOD PAIR ON ANISOTROPIC MESHES arxiv:1710.07857v1 [mat.na] 21 Oct 2017 GABRIEL R. BARRENECHEA AND ANDREAS WACHTEL Abstract. Uniform LBB conditions are
More informationA finite element approximation for the quasi-static Maxwell Landau Lifshitz Gilbert equations
ANZIAM J. 54 (CTAC2012) pp.c681 C698, 2013 C681 A finite element approximation for te quasi-static Maxwell Landau Lifsitz Gilbert equations Kim-Ngan Le 1 T. Tran 2 (Received 31 October 2012; revised 29
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationEXPLICIT ERROR BOUNDS IN A CONFORMING FINITE ELEMENT METHOD
MATHEMATICS OF COMPUTATION Volume 68, Number 228, Pages 1379 1396 S 0025-5718(99)01093-5 Article electronically publised on February 24, 1999 EXPLICIT ERROR BOUNDS IN A CONFORMING FINITE ELEMENT METHOD
More informationA SPLITTING LEAST-SQUARES MIXED FINITE ELEMENT METHOD FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALSIS AND MODELING Volume 3, Number 4, Pages 6 626 c 26 Institute for Scientific Computing and Information A SPLITTING LEAST-SQUARES MIED FINITE ELEMENT METHOD FOR
More informationResearch Article Discrete Mixed Petrov-Galerkin Finite Element Method for a Fourth-Order Two-Point Boundary Value Problem
International Journal of Matematics and Matematical Sciences Volume 2012, Article ID 962070, 18 pages doi:10.1155/2012/962070 Researc Article Discrete Mixed Petrov-Galerkin Finite Element Metod for a Fourt-Order
More informationH(div) conforming and DG methods for incompressible Euler s equations
H(div) conforming and DG metods for incompressible Euler s equations Jonny Guzmán Filánder A. Sequeira Ci-Wang Su Abstract H(div) conforming and discontinuous Galerkin (DG) metods are designed for incompressible
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationA Finite Element Primer
A Finite Element Primer David J. Silvester Scool of Matematics, University of Mancester d.silvester@mancester.ac.uk. Version.3 updated 4 October Contents A Model Diffusion Problem.................... x.
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationComputers and Mathematics with Applications. A nonlinear weighted least-squares finite element method for Stokes equations
Computers Matematics wit Applications 59 () 5 4 Contents lists available at ScienceDirect Computers Matematics wit Applications journal omepage: www.elsevier.com/locate/camwa A nonlinear weigted least-squares
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
More informationSuperconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract
Superconvergence of energy-conserving discontinuous Galerkin metods for linear yperbolic equations Yong Liu, Ci-Wang Su and Mengping Zang 3 Abstract In tis paper, we study superconvergence properties of
More informationChapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction
Chapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction The a posteriori error estimation of finite element approximations of elliptic boundary value problems has reached
More informationChemnitz Scientific Computing Preprints
G. Of G. J. Rodin O. Steinbac M. Taus Coupling Metods for Interior Penalty Discontinuous Galerkin Finite Element Metods and Boundary Element Metods CSC/11-02 Cemnitz Scientific Computing Preprints Impressum:
More informationNumerical analysis of a free piston problem
MATHEMATICAL COMMUNICATIONS 573 Mat. Commun., Vol. 15, No. 2, pp. 573-585 (2010) Numerical analysis of a free piston problem Boris Mua 1 and Zvonimir Tutek 1, 1 Department of Matematics, University of
More informationLEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS
SIAM J. NUMER. ANAL. c 998 Society for Industrial Applied Matematics Vol. 35, No., pp. 393 405, February 998 00 LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS YANZHAO CAO
More informationOverlapping domain decomposition methods for elliptic quasi-variational inequalities related to impulse control problem with mixed boundary conditions
Proc. Indian Acad. Sci. (Mat. Sci.) Vol. 121, No. 4, November 2011, pp. 481 493. c Indian Academy of Sciences Overlapping domain decomposition metods for elliptic quasi-variational inequalities related
More informationCOMPUTATIONAL COMPARISON BETWEEN THE TAYLOR HOOD AND THE CONFORMING CROUZEIX RAVIART ELEMENT
Proceedings of ALGORITMY 5 pp. 369 379 COMPUTATIONAL COMPARISON BETWEEN E TAYLOR HOOD AND E CONFORMING OUZEIX RAVIART ELEMENT ROLF KRAHL AND EBERHARD BÄNSCH Abstract. Tis paper is concerned wit te computational
More information1. Introduction. Consider a semilinear parabolic equation in the form
A POSTERIORI ERROR ESTIMATION FOR PARABOLIC PROBLEMS USING ELLIPTIC RECONSTRUCTIONS. I: BACKWARD-EULER AND CRANK-NICOLSON METHODS NATALIA KOPTEVA AND TORSTEN LINSS Abstract. A semilinear second-order parabolic
More information2009 Elsevier Science. Reprinted with permission from Elsevier.
P. Hansbo and M. Juntunen. 009. Weakly imposed Diriclet boundary conditions for te Brinkman model of porous media flow. Applied Numerical Matematics, volume 59, number 6, pages 174 189. doi:10.1016/j.apnum.008.07.003.
More information3 Parabolic Differential Equations
3 Parabolic Differential Equations 3.1 Classical solutions We consider existence and uniqueness results for initial-boundary value problems for te linear eat equation in Q := Ω (, T ), were Ω is a bounded
More informationNumerical Analysis of the Double Porosity Consolidation Model
XXI Congreso de Ecuaciones Diferenciales y Aplicaciones XI Congreso de Matemática Aplicada Ciudad Real 21-25 septiembre 2009 (pp. 1 8) Numerical Analysis of te Double Porosity Consolidation Model N. Boal
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More information2 Multi-Dimensional Variational Principles y s n u = α θ pu = β n Ω Figure 3..: Two-dimensional region wit and normal vector n x were
Capter 3 Multi-Dimensional Variational Principles 3. Galerkin's Metod and Extremal Principles Te construction of Galerkin formulations presented in Capters and 2 for one-dimensional problems readily extends
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More informationApproximation of incompressible large deformation elastic problems: some unresolved issues
Comput Mec DOI 10.1007/s00466-013-0869-0 ORIGINAL PAPER Approximation of incompressible large deformation elastic problems: some unresolved issues Ferdinando Auriccio Lourenço Beirão da Veiga Carlo Lovadina
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
More informationUniform estimate of the constant in the strengthened CBS inequality for anisotropic non-conforming FEM systems
Uniform estimate of te constant in te strengtened CBS inequality for anisotropic non-conforming FEM systems R. Blaeta S. Margenov M. Neytceva Version of November 0, 00 Abstract Preconditioners based on
More information6. Non-uniform bending
. Non-uniform bending Introduction Definition A non-uniform bending is te case were te cross-section is not only bent but also seared. It is known from te statics tat in suc a case, te bending moment in
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationCONSTRUCTIVELY WELL-POSED APPROXIMATION METHODS WITH UNITY INF SUP AND CONTINUITY CONSTANTS FOR PARTIAL DIFFERENTIAL EQUATIONS
MATHEMATICS OF COMPUTATION Volume 82, Number 284, October 203, Pages 923 952 S 0025-578(203)02697-X Article electronically publised on April 23, 203 CONSTRUCTIVELY WELL-POSED APPROXIMATION METHODS WITH
More informationOPTIMAL ERROR ESTIMATES FOR THE STOKES AND NAVIER STOKES EQUATIONS WITH SLIP BOUNDARY CONDITION
Matematical Modelling and Numerical Analysis M2AN, Vol. 33, N o 5, 1999, p. 923 938 Modélisation Matématique et Analyse Numérique OPTIMAL ERROR ESTIMATES FOR THE STOKES AND NAVIER STOKES EQUATIONS WITH
More information