A NEW CLASS OF MIXED FINITE ELEMENT METHODS FOR REISSNER MINDLIN PLATES
|
|
- Blaze Weaver
- 5 years ago
- Views:
Transcription
1 A NEW CLASS OF IXED FINITE ELEENT ETHODS FOR REISSNER INDLIN PLATES C. LOVADINA Abstract. A new class of finite elements for Reissner indlin plate problem is presented. The family is based on a modified mixed formulation recently introduced by Arnold and Brezzi. A result of stability and convergence uniformly in the thickness is provided. Key words. Plate, shear locking, finite element. AS subject classification. 65N30 1. Introduction. The Reissner indlin theory is widely used to describe the bending behaviour of an elastic plate loaded by a transverse force. Despite its simple formulation, it was suddenly noticed that discretization by means of finite elements was not at all straightforward, since in most cases numerical experiments exhibited a bad lack of convergence whenever the plate thickness was too small. Nowadays this undesirable phenomenon, known as shear locking (cf. [11], is well understood: as the thickness tends to zero, the Reissner indlin model enforces the Kirchhoff constraint, which is usually too severe on the discrete level when low order elements are employed. ore precisely, Kirchhoff constraint arises from the shear stress term which has to be taken into account in Reissner indlin formulation. The most popular way to overcome shear locking phenomenon is to reduce the influence of the shear energy by considering a suitable mixed formulation. A number of locking free elements based on this philosophy have been proposed and studied (cf. [6], [8]. Unfortunately, most of the mixed formulations require a clever choice of spaces of discretization in order to success the approximation and they do not generalize to shells. Although many methods for both Koiter and Naghdi shells are used in engineering practice (cf [4], [15], the problem of their convergence and stability has remained completely unsolved. Arnold and Brezzi (cf. [] have recently proposed a modification to the standard mixed formulation of Reissner indlin plate problem which seems to be capable to allow a wider choice of elements. ore importantly, its generalization to the Naghdi shell problem is the only one for which a first theoretical analysis has been developped so far (cf. [3]. Having this in mind, we face the study of elements for Reissner indlin plates based on the modified formulation as a first step of research, firmly believing that they can lead to competitive methods for shells. In this paper we present a new class of finite elements based on the formulation given in []. For our family we prove stability and convegence uniformly in the thickness. An outline of the paper is as follows. In Section we first recall the Reissner indlin model. Then we briefly introduce both the standard mixed formulation and the modified one. Section 3 is the core of the paper. By using a macroelement technique we prove a general stability result. The key idea is the possibility, under suitable hypotheses, to construct an interpolation operator à la Fortin (cf. [7] not componentwise but mixing the components of rotations and vertical displacements. Istituto di Analisi Numerica del C.N.R., Via Abbiategrasso 09, 7100 Pavia, Italia. 1
2 c. lovadina In Section 4 we will apply the results of Section 3 to the case of the family P k+ P k+ P k, showing that conditions of Theorem 3.1 are fulfilled for this choice. We conclude the paper by stating a uniform convergence result which follows immediately from our analysis and standard mixed finite element theory. An appendix in which we briefly discuss the implementation features of our methods is also provided. Our notations agree with those commonly used in the literature (cf. [7] and [10], for instance.. The Reissner indlin model and the modified mixed formulation. Let us denote with A = Ω ( t/, t/ the region in R 3 occupied by an undeformed elastic plate of thickness t > 0. It s well known that the Reissner indlin plate model describes the bending behaviour of the plate in terms of the transverse displacements and of the fiber rotations normal to the midplane Ω. If we restrict our attention to the case of a clamped plate, the problem consists in finding the couple (θ(t, w(t that minimizes the following functional Π t (θ(t, w(t = 1 λt a(θ(t, θ(t + θ(t w(t 0,Ω fw(t dx dy (.1 Ω over the space V = Θ W = (H0 1(Ω H0 1 (Ω, where (i θ = (θ 1, θ and w are the fiber rotations and the transverse displacements, respectively; (ii a(, : Θ Θ R is a bilinear form defined by a(θ, η = {( E θ1 1(1 ν Ω x + ν θ ( η1 y x + ν θ 1 x + θ η y y ν ( θ1 y + θ ( η1 x y + η dx dy x where E is the Young s modulus and ν is the Poisson s ratio (0 < ν < 1/; (iii λ = Ek/(1 + ν with k shear correction factor (usually taken as 5/6; (iv ft 3 is the transverse load force density per unit area. By Korn s inequality, a(, is a coercive form over Θ and this implies that there exists a unique solution (θ(t, w(t in V of the Problem P t :. For t > 0 fixed, find (θ(t, w(t in V such that a(θ(t, η + λt (θ(t w(t, η v = Ω fv dx dy (η, v V (. The idea of discretizing by means of finite elements based on formulation (. is typically a very naive one, because of the well known shear locking phenomenon (cf. [11]. A common way to escape this undesirable lack of convergence is to discretize a mixed formulation derived by introducing the scaled shear stress ξ = λt (θ w
3 a new class of finite elements for plates 3 as independent unknown. functional This leads to consider the saddle point problem for the Π t (θ, w, ξ = 1 a(θ, θ λ 1 t ξ 0,Ω + (ξ, θ w (f, w (.3 on V (L (Ω. Hence the variational plate problem now reads as follows: Problem P t :. For t > 0 fixed, find (θ(t, w(t, ξ(t in V (L (Ω such that a(θ(t, η + (ξ, η v = (f, v (η, v V (.4 (τ, θ(t w(t λ 1 t (ξ, τ = 0 τ (L (Ω In agreement with the notations of [7], let us now introduce the differential operators { ϕ rot : ϕ rot ϕ = y, ϕ x rot : χ rot χ = χ 1 y + χ x. In what follows, we also need the Hilbert space Γ = H 0 (rot; Ω defined by H 0 (rot; Ω = { χ : χ (L (Ω, rot χ L (Ω, χ t = 0 on Ω χ H 0(rot;Ω := χ 0,Ω + rot χ 0,Ω (here t is the unit tangent to Ω and its dual space { Γ = H 1 (div ; Ω = γ : γ (H 1 (Ω div γ H 1 (Ω. The space Γ will be equipped with the norm γ Γ := γ 1,Ω + div γ 1,Ω which is equivalent to the natural dual norm induced by H 0 (rot ; Ω. We are now ready to claim that the following proposition holds (for a proof see [7], for instance. Proposition.1. Given t > 0, there exists a unique triple (θ(t, w(t, ξ(t in V (L (Ω satisfying equations (.4. Furthermore θ(t θ 0 in Θ w(t w 0 in W ξ(t ξ 0 in Γ and (θ 0, w 0, ξ 0 in V Γ solves a(θ 0, η < ξ 0, η v > Γ Γ= (f, v (η, v V (.5 < τ, θ 0 w 0 > Γ Γ= 0 τ Γ
4 4 c. lovadina Therefore, from the general theory of mixed methods, a choice V h Γ h V Γ of finite element spaces will{ be stable if it satisfies (cf. [7]: (S.1 (η h, v h Vh = (η h, v h V h : (η h v h, τ h = 0 τ h Γ h : a(η h, η h α 0 ( η h 1,Ω + v h 1,Ω (S. τ h Γ h : (η sup h v h, τ h β 0 τ (η,v h V h (η h h, v h h Γ 1,Ω with α 0 and β 0 independent of h. Remark. Even though we only require the inclusion Γ h Γ, it s obvious that every natural choice for Γ h will always satisfy Γ h (L (Ω. Note that finding finite element spaces that fulfill both condition (S.1 and (S. is not so trivial since a(, considered as a bilinear form on V is coercive only on V = { (η, v V : η = v and not on the whole V. Hence, if we try to enrich V h in order to get (S., then (S.1 may fail; on the other hand if we let V h be small with respect to Γ h, condition (S.1 will become easier to verify but (S. may be false. Recently, Arnold and Brezzi (cf. [] have proposed a modification to the mixed formulation (.4 by taking as Lagrange multiplier and setting γ = λ(t α(θ w ã(θ, w; η, v = a(θ, η + λα(θ w, η v where α is a parameter such that 0 < α < t (cf. Appendix. Therefore the weak formulation for the plate problem turns out to be ã(θ, w; η, v + (γ, η v = (f, v (η, v V t (s, θ w (.6 λ(1 αt (γ, s = 0 s (L (Ω By Poincaré s inequality and Korn s inequality it s easy to show that ã(, is indeed bounded and coercive on the whole V. This means that in choosing finite element spaces for the discretization of (.6 we do not have to take care of condition (S.1 anymore: any choice satisfying (S. will be stable and the following proposition is true (cf. [7]. Proposition.. Let us choose Θ h Θ, W h W, Γ h Γ and consider the discretized problem of (.6: Find (θ h, w h, γ h Θ h W h Γ h = V h Γ h, solution of equations ã(θ h, w h ; η h, v h + (γ h, η h v h = (f, v h (η h, v h V h t (s h, θ h w h (.7 λ(1 αt (γ, s h h = 0 s h Γ h Problem (.7 has a unique solution; moreover, if condition (S. is satisfied and (θ, w, γ is the solution of (.6, then there exists a constant C independent of t and h such that
5 a new class of finite elements for plates 5 θ θ h 1,Ω + w w h 1,Ω + γ γ h t C ( inf η η η Θ h h 1,Ω + inf w v h 1,Ω + inf γ s h t v h W h s h h Γ h (.8 where γ t := γ 1,Ω + div γ 1,Ω + t γ 0,Ω 3. A stability result. Let us start by recalling some standard definitions and notations useful to the subsequent analysis. First of all, we say that a family {T h h>0 of triangulations of Ω is regular if there exists a constant σ > 0 such that h T σρ T T h>0 T h where h T is the diameter of triangle T and ρ T is the maximum diameter of the circles contained in T. Furthermore, a macroelement is the union of a fixed number of neighboring triangles along a well-defined pattern. A macroelement = m i=1 T i is said to be equivalent to a reference macroelement = m i=1 T i if there is a mapping F : for which the following conditions are fulfilled (cf. [13]: (1 F is a continuous bijection. ( T i = F ( T i i 1 i m (3 F T i = F Ti F 1 T i where F Ti and F Ti are the affine functions from the standard reference triangle onto T i and T i, respectively. From a given triangulation T h of Ω it s always possible to derive (obviously not in a unique fashion a macroelement mesh h in such a way that each T T h is covered by some macroelement in h and each macroelement is equivalent to a certain reference macroelement. Let us now introduce finite element spaces V h Γ h V Γ for the discretization of (.6. For every macroelement in h we will denote with V 0, the space V 0, = { (η h, v h V h : (η h, v h = (0, 0 in Ω \. We further define the space Γ as Finally, let us set Γ = { s h Γ h : s h = 0 in Ω \. N = { s Γ : (s, η v = 0 (η, v V 0,. In what follows, we will always assume that the family {T h h>0 arising from finite element discretization will be regular; moreover, there will be only a fixed finite number of reference macroelements { 1,..., r to which each macroelement h>0 h is equivalent. We are now in the position to state and prove our main result of this section. Theorem 3.1. Given Θ h Θ, W h W, Γ h Γ subspaces for discretization of problem (.6, suppose that the following conditions are met:
6 6 c. lovadina (C1 There exists a linear operator π h : W W h such that (C1.1 (C1. π h v 1,Ω c v 1,Ω, c independent of h Ω (π hv v s h = 0 s h Γ h (C From T h it is possible to derive h in such a way that (C.1 Γ h = h Γ (direct sum in L (C. each N is the trivial space. Then condition (S. of the previous section is fulfilled. Proof. Let (η, v V be given. Fix an arbitrary macroelement h ; let us denote with i the index 1 i r such that is equivalent to i and with h = max i h Ti if = m i T i. Consider the problem to find (η B, vb V 0, solution of (η B vb s = where Π 1 η is the Clément interpolated of η. (η Π 1 η s s Γ, (3.1 By (C. it follows that system (3.1 is solvable. Let us take the solution of minimum norm. An easy scaling argument and the features of the Clément interpolating operator (cf. [7] show that there exists c( i > 0 such that η B 1, + h vb 1, c( i η 1, (3. Since h is obviously bounded by mes(ω, inequality (3. implies that c 1 ( i > 0 : η B 1, + vb 1, c 1( i η 1,. (3.3 Now let us set η I = Π 1 η + v I = π h v + η B v B (3.4 Thanks to (C.1, every s h Γ h can be written in a unique fashion as s h = s h, where s h, Γ ; hence we have
7 a new class of finite elements for plates 7 (η I v I s h = (by (3.4 Ω = = [ [ ] Π 1 η + η B (vb + π hv s h, = ] (Π 1 η + η B vb s h, π h v s h, = = (by (C1. and (3.1 = ( η s h, v s h, = = (again (C.1 = Ω (η v s h Hence for every (η, v V we have found Π h (η, v = (η I, v I V h such that (η I v I s h = (η v s h s h Γ h (3.5 Ω Ω Let us now estimate η I 1,Ω + vi 1,Ω : η I 1,Ω + vi 1,Ω = Π1 η + η B + π h v + 1,Ω v B 1,Ω ( Π 1 η 1,Ω + π hv 1,Ω + ( η B 1, + vb 1, ( ( (by 3.3 and (C1.1 c v 1,Ω + η 1,Ω + c 1 ( i η 1,. Note that, above, we have also used the continuity of the operator Π 1. Since there s only a finite number of reference macroelements, we obtain η I 1,Ω + vi 1,Ω C 1( v 1,Ω + η 1,Ω (3.6 with C 1 = max {c, c 1 ( i,..., c 1 ( i. Therefore, by Poincaré s inequality, we finally have Π h (η, v 1,Ω C (η, v 1,Ω (3.7 with C independent of h. It s well-known (cf. [7], for instance that (3.5 together with (3.7 implies condition (S. of the previous section and the proof is complete. Remark. Note that in Theorem 3.1 we have done nothing but employing Fortin s criterion (cf. [7]. So far every interpolation operator built in order to meet this
8 8 c. lovadina criterion split into Π h = (π 1, π where π 1 : Θ Θ h and π : W W h ; hence the operator acted componentwise. Here we have presented a construction of Π h for which this splitting is not necessary (cf. equation (3.1. There are several cases where the difficulty consists in building π 1 : in such a situation our result may be of great help (cf. Section Analysis of the class (P k+ P k+ P k. In this section we will apply our results to show convergence and stability of a new class of mixed finite element methods based on formulation (.6. For every k 0 let us introduce the following spaces of discretization: Θ h = { η h Θ : η h T (P k+ (T T T h W h = { v h W : v h T P k+ (T T T h (4.1 Γ h = { s h (L (Ω : s h T (P k (T T T h The goal of the next proposition is to verify, under an extremely weak hypothesis on T h, that conditions of Theorem 3.1 are fulfilled. Proposition 4.1. Let us choose the approximation spaces as in (4.1 and let us suppose that every h>0 h is the union of at least three triangles. Then conditions of Theorem (3.1 are met. Proof. It s easily seen that we can build π h : W W h such that (C1.1 and (C1. are true. Indeed, set first πh 1 : W W h as the interpolating operator of Clément; we have (cf. [7] h s T v πh 1v s,t c v 1,Ω s = 0, 1 (4. T T h Secondly, for every v W take πh v W h satisfying: πh v(a = 0 a vertex of T T h (πh v vp k = 0 p k P k (e e edge of T T h e (πh v vq k 1 = 0 q k 1 P k 1 (T for k > 0 Standard arguments show that If we now set T π h v 1,T c ( h 1 T v 0,T + v 1,T (4.3 (4.4 π h v = π 1 h v + π h( v π 1 h v (4.5 it s easy to see from (4., (4.3 and (4.4 that π h : W W h is a linear bounded operator for which (C1.1 and (C1. are fulfilled. As far as condition (C.1 is concerned, note that s h Γ h is fully discontinuous along adjacent triangles, thus (C.1 is trivially satisfied. Therefore it remains to verify (C..
9 a new class of finite elements for plates 9 Fix a macroelement = m i=1 T i and take s N. By choosing (0, v V 0,, integrating by parts yields 0 = v s = m i=1 { v div s + v s n Ti T i T i (4.6 Since equation (4.6 is true for every (0, v V 0, vanishing on the boundary of macroelement, it follows that s N has to satisfy div s = 0 in each T i (4.7 s n is continuous along the internal interfaces of T i Hence there exists a unique ϕ such that { ϕ L 1 k+1 (/R s = rot ϕ (4.8 where L 1 k+1 ( is the usual space of continuous functions which are polynomials of degree k+1 on each T i of. Let us now take (η, 0 V 0,. Since s N, we have η s = 0 (η, 0 V 0,. (4.9 By (4.8 we obtain η rot ϕ = 0 (η, 0 V 0, (4.10 from which we wish to deduce rot ϕ = 0. To perform this step, consider on the Stokes like problem given χ (L (, find (u, p ( L 1 k+ ( L 1 k+1 (/R such that u : v p rot v = χ v v ( L 1 k+ ( (4.11 u rot q = 0 q L 1 k+1 ( where L 1 k+ ( = L1 k+ ( H1 0 (. D. Boffi (cf. [5] has shown that discretization of the Stokes problem by means of a standard P k+ element for velocity and a standard continuous P k+1 element for pressure is stable for every k 0, whenever in the mesh there are at least three elements. It s obvious that this is also true for problem (4.11. In particular, stability implies that there are not spurious pressure modes, i.e. { rot q v = 0 v ( L 1 k+ ( = rot q = 0 (4.1
10 10 c. lovadina Compared to (4.10, implication (4.1 is exactly what we are looking for. s N s = 0 and the proof is now complete. Hence Remark. Note that it s always possible to derive, from a given regular family {T h h>0, a macroelement family { h h>0 which fulfills the assumption of Proposition 4.1 and for which there are only a finite number of reference macroelements, provided in each T h there are at least three triangles. Recalling Proposition. and standard interpolation theory, it s straightforward to obtain the following error bound Proposition 4.. Let (θ, w, γ be the solution of problem (.6 and (θ h, w h, γ h be the solution of the discretized problem by means of the spaces (4.1, then there exists a constant C independent of h and t such that θ θ h 1,Ω + w w h 1,Ω + γ γ h t Ch k+1 ( θ k+,ω + w k+,ω + γ k+1,ω ( Appendix. A brief discussion about implementation of the methods. Let us start by recalling the modified mixed formulation (.6: given f L (Ω and t > 0, find (θ, w, γ V (L (Ω such that ã(θ, w; η, v + (γ, η v = (f, v (η, v V t (s, θ w (A.1 λ(1 αt (γ, s = 0 s (L (Ω and its discrete counterpart find (θ h, w h, γ h V h Γ h such that ã(θ h, w h ; η h, v h + (γ h, η h v h = (f, v h (η h, v h V h t (s h, θ h w h λ(1 αt (γ, s h h = 0 s h Γ h (A. We have already noticed that α, which has the units of reciprocal of the square of length, is a parameter with 0 < α < t. If we now choose the spaces Θ h, W h, Γ h as in (4.1, from the second equation of (A. we deduce that γ h = λ(t αp h (θ h w h (A.3 where P h denotes the L orthogonal projection operator onto Γ h. The variable γ h can be eliminated from equations (A. and it follows that (θ h, w h V h is the minimizer over V h of the functional E h (η, v = 1 a(η, η + λα η v 0,Ω + λ (t α P h (η v 0,Ω (f, v (A.4 From a computational point of view is much more convenient to implement the Euler Lagrange equations arising from functional (A.4 instead of solving equations (A..
11 a new class of finite elements for plates 11 Let us now turn our attention on the parameter α. Looking at (A.4, we find that choosing a very small α means essentially that we have decided to treat the shear energy term by a reduced integration procedure. On the other hand, a choice of α close to t means that the shear energy term is exactly integrated. Therefore, we expect that a very small α may lead to spurious modes, while the discretization may exhibit a severe locking whenever α is too large. Hence a suitable choice of the parameter has to be done in order to achieve acceptable convergence results. Our guess is that the parameter α may be chosen so that the first two terms in the right hand side of (A.4 are well balanced. Numerical experiments for the case P P P 0, i.e. the lowest order scheme of the family described in section 4, show that such a behaviour happens, indeed. Just to give an idea, we report here some numerical results concerning the discretization of a clamped square plate of unit length side and of thickness t = We used a uniform decomposition into 64 equal subsquares each of one divided into two triangles by a diagonal, thus leading to a triangular mesh of 18 elements. Due to symmetry, only one quadrant has been discretized in actual computations. Young s modulus and Poisson s ratio have been chosen as E = and ν = 0.3, respectively. Finally, two different loads have been selected. The first one constists in choosing a uniform load f = 1. It turns out that exact solution of vertical displacement variable at the center of the plate is (cf. [14]: w c = We have obtained, for different values of α, the following results: α = 10 5 α = 10 4 α = 10 3 α = 10 α = 10 α = 1 α = 10 5 w c = w c = w c = w c = w c = w c = w c = The second load is of purely academic interest: it has been designed in order to have the exact solution explicitly known in the whole square domain. This choice allows to compute the discrete L -norm relative error for both rotations and deflections (cf. [9], as displayed below:
12 1 c. lovadina α = 10 5 E θ = E w = α = 10 4 E θ = E w = α = 10 3 E θ = E w = α = 10 E θ = E w = α = 10 E θ = E w = α = 1 E θ = E w = α = 10 5 E θ = E w = Note that, within a reasonable range of values, the choice of α is not so crucial, showing that the method promises to be robust. For further numerical results about schemes based on the modified formulation used in the present work, we refer to [9]; in particular more detailed recipes for a good choice of the parameter α have been therein suggested. REFERENCES [1] D.N. Arnold, Innovative finite element methods for plates, ath. Applic. Comp., V.10 (1991, pp [] D.N. Arnold and F. Brezzi, Some new elements for the Reissner indlin plate model, in Boundary value problems for partial differential equations and applications, J.L. Lions and C. Baiocchi, ed.. [3] D.N. Arnold and F. Brezzi, Locking free finite elements for shells, to appear. [4] K.J. Bathe and E. Wilson, Numerical methods in finite elements analysis, Englewood Cliffs, NJ, [5] D. Boffi, Stability of higher order triangular Hood Taylor methods for the stationary Stokes equations, ath. odels and ethods in Appl. Sci., to appear. [6] F. Brezzi, K.J. Bathe and. Fortin, ixed interpolated elements for Reissner indlin plates, Internat. J. Numer. ethods Engrg., 8 (1989, pp [7] F. Brezzi and. Fortin, ixed and Hybrid Finite Element ethod, Springer-Verlag, [8] F. Brezzi,. Fortin and R. Stenberg, Error analysis of mixed interpolated elements for Reissner indlin plates, ath. odels and ethods in Appl. Sci., (1991, pp [9] C. Chinosi and C. Lovadina, Numerical analysis of some mixed finite element methods for Reissner indlin plates, Comput. echanics, to appear. [10] P.G. Ciarlet, The Finite Element ethod for Elliptic Equations, North Holland, [11] T.J.R. Hughes, The Finite Element ethod, Englewood Cliffs NY, [1] J. Pitkäranta, Analysis of some low order finite element schemes for indlin Reissner and Kirchhoff plates, Numer. ath., 53 (1988, pp [13] R. Stenberg, Analysis of mixed finite element methods for the Stokes problem: a unified approach, ath. Comp., 4 (1984, pp [14] S. Timoshenko and S. Woinowki Kriegel, Theory of plates and shells, cgraw Hill, Singapore, [15] O.C. Zienkiewicz and R.L. Taylor, The finite element methods, cgraw Hill, New York, NY, 1989.
Some New Elements for the Reissner Mindlin Plate Model
Boundary Value Problems for Partial Differential Equations and Applications, J.-L. Lions and C. Baiocchi, eds., Masson, 1993, pp. 287 292. Some New Elements for the Reissner Mindlin Plate Model Douglas
More informationANALYSIS OF MIXED FINITE ELEMENTS FOR LAMINATED COMPOSITE PLATES
ANALYSIS OF MIXED FINITE ELEMENTS FOR LAMINATED COMPOSITE PLATES F. Auricchio Dipartimento di Meccanica Strutturale Università di Pavia, Italy C. Lovadina Dipartimento di Ingegneria Meccanica e Strutturale
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationINNOVATIVE FINITE ELEMENT METHODS FOR PLATES* DOUGLAS N. ARNOLD
INNOVATIVE FINITE ELEMENT METHODS FOR PLATES* DOUGLAS N. ARNOLD Abstract. Finite element methods for the Reissner Mindlin plate theory are discussed. Methods in which both the tranverse displacement and
More informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
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 informationLocking phenomena in Computational Mechanics: nearly incompressible materials and plate problems
Locking phenomena in Computational Mechanics: nearly incompressible materials and plate problems C. Lovadina Dipartimento di Matematica Univ. di Pavia IMATI-CNR, Pavia Bologna, September, the 18th 2006
More informationMIXED FINITE ELEMENTS FOR PLATES. Ricardo G. Durán Universidad de Buenos Aires
MIXED FINITE ELEMENTS FOR PLATES Ricardo G. Durán Universidad de Buenos Aires - Necessity of 2D models. - Reissner-Mindlin Equations. - Finite Element Approximations. - Locking. - Mixed interpolation or
More informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationA-priori and a-posteriori error estimates for a family of Reissner-Mindlin plate elements
A-priori and a-posteriori error estimates for a family of Reissner-Mindlin plate elements Joint work with Lourenco Beirão da Veiga (Milano), Claudia Chinosi (Alessandria) and Carlo Lovadina (Pavia) Previous
More informationA UNIFORMLY ACCURATE FINITE ELEMENT METHOD FOR THE REISSNER MINDLIN PLATE*
SIAM J. NUMER. ANAL. c 989 Society for Industrial and Applied Mathematics Vol. 26, No. 6, pp. 5, December 989 000 A UNIFORMLY ACCURATE FINITE ELEMENT METHOD FOR THE REISSNER MINDLIN PLATE* DOUGLAS N. ARNOLD
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 informationANALYSIS OF A LINEAR LINEAR FINITE ELEMENT FOR THE REISSNER MINDLIN PLATE MODEL
Mathematical Models and Methods in Applied Sciences c World Scientific Publishing Company ANALYSIS OF A LINEAR LINEAR FINITE ELEMENT FOR THE REISSNER MINDLIN PLATE MODEL DOUGLAS N. ARNOLD Department of
More informationFind (u,p;λ), with u 0 and λ R, such that u + p = λu in Ω, (2.1) div u = 0 in Ω, u = 0 on Γ.
A POSTERIORI ESTIMATES FOR THE STOKES EIGENVALUE PROBLEM CARLO LOVADINA, MIKKO LYLY, AND ROLF STENBERG Abstract. We consider the Stokes eigenvalue problem. For the eigenvalues we derive both upper and
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationENHANCED STRAIN METHODS FOR ELASTICITY PROBLEMS
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2004 P. Neittaanmäki, T. Rossi, S. Korotov, E. Oñate, J. Périaux, and D. Knörzer (eds.) Jyväskylä, 24 28 July 2004
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
MATHEMATICS OF COMPUTATION Volume 75, Number 256, October 2006, Pages 1659 1674 S 0025-57180601872-2 Article electronically published on June 26, 2006 ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED
More informationDominique Chapelle 1, Anca Ferent 1 and Patrick Le Tallec 2
Mathematical Modelling and Numerical Analysis ESAIM: M2AN Modélisation Mathématique et Analyse Numérique M2AN, Vol. 37, N o, 2003, pp. 43 58 DOI: 0.05/m2an:200305 THE TREATMENT OF PINCHING LOCKING IN 3D-SHELL
More informationMULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract
More informationb i (x) u + c(x)u = f in Ω,
SIAM J. NUMER. ANAL. Vol. 39, No. 6, pp. 1938 1953 c 2002 Society for Industrial and Applied Mathematics SUBOPTIMAL AND OPTIMAL CONVERGENCE IN MIXED FINITE ELEMENT METHODS ALAN DEMLOW Abstract. An elliptic
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 informationA NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION
A NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION JOHNNY GUZMÁN, ABNER J. SALGADO, AND FRANCISCO-JAVIER SAYAS Abstract. The analysis of finite-element-like Galerkin discretization techniques for the
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
More informationANALYSIS OF A MIXED-SHEAR-PROJECTED QUADRILATERAL ELEMENT METHOD FOR REISSNER-MINDLIN PLATES
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 14, Number 1, Pages 48 62 c 2017 Institute for Scientific Computing and Information ANALYSIS OF A MIXED-SHEAR-PROJECTED QUADRILATERAL ELEMENT
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationTwo Nonconforming Quadrilateral Elements for the Reissner-Mindlin Plate
Two Nonconforming Quadrilateral Elements for the Reissner-Mindlin Plate Pingbing Ming and Zhong-ci Shi Institute of Computational Mathematics & Scientific/Engineering Computing, AMSS, Chinese Academy of
More informationFinite Elements for Elastic Shell Models in
Elastic s in Advisor: Matthias Heinkenschloss Computational and Applied Mathematics Rice University 13 April 2007 Outline Elasticity in Differential Geometry of Shell Geometry and Equations The Plate Model
More informationNONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS
NONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS J.-L. GUERMOND 1, Abstract. This paper analyzes a nonstandard form of the Stokes problem where the mass conservation
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationMathematical modelling and numerical analysis of thin elastic shells
Mathematical modelling and numerical analysis of thin elastic shells Zhang Sheng Wayne State University, Michigan Zhang Sheng (szhang@wayne.edu) Thin elastic shell May, 202 / 9 Models for mechanics of
More informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural
More informationNumerische Mathematik
Numer. Math. (2012) 122:61 99 DOI 10.1007/s00211-012-0456-x Numerische Mathematik C 0 elements for generalized indefinite Maxwell equations Huoyuan Duan Ping Lin Roger C. E. Tan Received: 31 July 2010
More informationERROR ESTIMATES FOR LOW-ORDER ISOPARAMETRIC QUADRILATERAL FINITE ELEMENTS FOR PLATES
ERROR ESTIATES FOR LOW-ORDER ISOPARAETRIC QUADRILATERAL FINITE ELEENTS FOR PLATES RICARDO G DURÁN, ERWIN HERNÁNDEZ, LUIS HERVELLA-NIETO, ELSA LIBERAN, AND RODOLFO RODRíGUEZ Abstract This paper deals with
More informationA Finite Element Method for the Surface Stokes Problem
J A N U A R Y 2 0 1 8 P R E P R I N T 4 7 5 A Finite Element Method for the Surface Stokes Problem Maxim A. Olshanskii *, Annalisa Quaini, Arnold Reusken and Vladimir Yushutin Institut für Geometrie und
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationA LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES
A LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES P. HANSBO Department of Applied Mechanics, Chalmers University of Technology, S-4 96 Göteborg, Sweden E-mail: hansbo@solid.chalmers.se
More informationANALYSIS OF SOME LOW ORDER QUADRILATERAL REISSNER-MINDLIN PLATE ELEMENTS
MATHEMATICS OF COMPUTATION Volume 00, Number 0, Pages 000 000 S 0025-5718(XX)0000-0 ANALYSIS OF SOME LOW ORDER QUADRILATERAL REISSNER-MINDLIN PLATE ELEMENTS PINGBING MING AND ZHONG-CI SHI Abstract. Four
More informationNONCONFORMING MIXED ELEMENTS FOR ELASTICITY
Mathematical Models and Methods in Applied Sciences Vol. 13, No. 3 (2003) 295 307 c World Scientific Publishing Company NONCONFORMING MIXED ELEMENTS FOR ELASTICITY DOUGLAS N. ARNOLD Institute for Mathematics
More informationA primer on Numerical methods for elasticity
A primer on Numerical methods for elasticity Douglas N. Arnold, University of Minnesota Complex materials: Mathematical models and numerical methods Oslo, June 10 12, 2015 One has to resort to the indignity
More informationRemarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable?
Remarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable? Thomas Apel, Hans-G. Roos 22.7.2008 Abstract In the first part of the paper we discuss minimal
More informationDiscontinuous Petrov-Galerkin Methods
Discontinuous Petrov-Galerkin Methods Friederike Hellwig 1st CENTRAL School on Analysis and Numerics for Partial Differential Equations, November 12, 2015 Motivation discontinuous Petrov-Galerkin (dpg)
More informationMIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS
MIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS JUHO KÖNNÖ, DOMINIK SCHÖTZAU, AND ROLF STENBERG Abstract. We derive new a-priori and a-posteriori error estimates for mixed nite
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationLeast-squares Finite Element Approximations for the Reissner Mindlin Plate
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl., 6, 479 496 999 Least-squares Finite Element Approximations for the Reissner Mindlin Plate Zhiqiang Cai, Xiu Ye 2 and Huilong Zhang
More informationFinite Element Approximation of Piezoelectric Plates
Finite Element Approximation of Piezoelectric Plates F. Auricchio Department of Structural Mechanics, University of Pavia, 27100 Pavia, Italy P. Bisegna Department of Civil Engineering, University of Roma
More informationPREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO
PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationLOCKING-FREE FINITE ELEMENT METHODS FOR SHELLS
MATHEMATICS OF COMPUTATION Volume 66, Number 217, January 1997, Pages 1 14 S 0025-5718(97)00785-0 LOCKING-FREE FINITE ELEMENT METHODS FOR SHELLS DOUGLAS N. ARNOLD AND FRANCO BREZZI Abstract. We propose
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationMultigrid Methods for Saddle Point Problems
Multigrid Methods for Saddle Point Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University Advances in Mathematics of Finite Elements (In
More informationQuestions (And some Answers) (And lots of Opinions) On Shell Theory
1 Questions (And some Answers) (And lots of Opinions) On Shell Theory 2 Shells Q: What is a shell? A: A three-dimensional elastic body occupying a thin neighborhood of a two-dimensional submanifold of
More informationINF-SUP CONDITION FOR OPERATOR EQUATIONS
INF-SUP CONDITION FOR OPERATOR EQUATIONS LONG CHEN We study the well-posedness of the operator equation (1) T u = f. where T is a linear and bounded operator between two linear vector spaces. We give equivalent
More informationIMPLEMENTATION OF FUNCTIONAL-TYPE A POSTERIORI ERROR ESTIMATES FOR LINEAR PROBLEMS IN SOLID MECHANICS. Maksim Frolov
IMPLEMENTATION OF FUNCTIONAL-TYPE A POSTERIORI ERROR ESTIMATES FOR LINEAR PROBLEMS IN SOLID MECHANICS Maksim Frolov Peter the Great St.Petersburg Polytechnic University AANMPDE11, FINLAND 06 10 August
More informationMathematics FINITE ELEMENT ANALYSIS AS COMPUTATION. What the textbooks don't teach you about finite element analysis. Chapter 3
Mathematics FINITE ELEMENT ANALYSIS AS COMPUTATION What the textbooks don't teach you about finite element analysis Chapter 3 Completeness and continuity: How to choose shape functions? Gangan Prathap
More informationNumerical Methods for the Navier-Stokes equations
Arnold Reusken Numerical Methods for the Navier-Stokes equations January 6, 212 Chair for Numerical Mathematics RWTH Aachen Contents 1 The Navier-Stokes equations.............................................
More informationA DECOMPOSITION RESULT FOR BIHARMONIC PROBLEMS AND THE HELLAN-HERRMANN-JOHNSON METHOD
Electronic ransactions on Numerical Analysis. Volume 45, pp. 257 282, 2016. Copyright c 2016,. ISSN 1068 9613. ENA A DECOMPOSIION RESUL FOR BIHARMONIC PROBLEMS AND HE HELLAN-HERRMANN-JOHNSON MEHOD WOLFGANG
More informationA posteriori error estimates for non conforming approximation of eigenvalue problems
A posteriori error estimates for non conforming approximation of eigenvalue problems E. Dari a, R. G. Durán b and C. Padra c, a Centro Atómico Bariloche, Comisión Nacional de Energía Atómica and CONICE,
More informationFinite Elements for the Stokes Problem
Finite Elements for the Stokes Problem Daniele Boffi 1, Franco Brezzi 2, and Michel Fortin 3 1 Dipartimento di Matematica F. Casorati, Università degli studi di Pavia, Via Ferrata 1, 27100 Pavia, Italy
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationNumerische Mathematik
Numer. Math. 69: 59 522 (1995) Numerische Mathematik c Springer-Verlag 1995 Locking and robustness in the finite element method for circular arch problem Zhimin Zhang Department of Mathematics, Texas Tech
More informationKonstantinos Chrysafinos 1 and L. Steven Hou Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ANALYSIS AND APPROXIMATIONS OF THE EVOLUTIONARY STOKES EQUATIONS WITH INHOMOGENEOUS
More informationOn the Numerical Modelling of Orthotropic Large Strain Elastoplasticity
63 Advances in 63 On the Numerical Modelling of Orthotropic Large Strain Elastoplasticity I. Karsaj, C. Sansour and J. Soric Summary A constitutive model for orthotropic yield function at large strain
More informationFinite Elements. Colin Cotter. January 18, Colin Cotter FEM
Finite Elements January 18, 2019 The finite element Given a triangulation T of a domain Ω, finite element spaces are defined according to 1. the form the functions take (usually polynomial) when restricted
More informationA posteriori error estimates applied to flow in a channel with corners
Mathematics and Computers in Simulation 61 (2003) 375 383 A posteriori error estimates applied to flow in a channel with corners Pavel Burda a,, Jaroslav Novotný b, Bedřich Sousedík a a Department of Mathematics,
More informationGlowinski Pironneau method for the 3D ω-ψ equations
280 GUERMOND AND QUARTAPELLE Glowinski Pironneau method for the 3D ω-ψ equations Jean-Luc Guermond and Luigi Quartapelle 1 LIMSI CNRS, Orsay, France, and Dipartimento di Fisica, Politecnico di Milano,
More informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
More informationNumerische Mathematik
Numer. Math. 2015 130:395 423 DOI 10.1007/s00211-014-0672-7 Numerische Mathematik Weakly over-penalized discontinuous Galerkin schemes for Reissner Mindlin plates without the shear variable Paulo Rafael
More informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
More informationChapter 2 Finite Element Spaces for Linear Saddle Point Problems
Chapter 2 Finite Element Spaces for Linear Saddle Point Problems Remark 2.1. Motivation. This chapter deals with the first difficulty inherent to the incompressible Navier Stokes equations, see Remark
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationOctober 25, 2013 INNER PRODUCT SPACES
October 25, 2013 INNER PRODUCT SPACES RODICA D. COSTIN Contents 1. Inner product 2 1.1. Inner product 2 1.2. Inner product spaces 4 2. Orthogonal bases 5 2.1. Existence of an orthogonal basis 7 2.2. Orthogonal
More informationSTOKES PROBLEM WITH SEVERAL TYPES OF BOUNDARY CONDITIONS IN AN EXTERIOR DOMAIN
Electronic Journal of Differential Equations, Vol. 2013 2013, No. 196, pp. 1 28. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu STOKES PROBLEM
More informationUNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES
UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
More informationMATHEMATICAL AND NUMERICAL ANALYSIS OF SOME PARAMETER-DEPENDENT MODELS
MATHEMATICAL AND NUMERICAL ANALYSIS OF SOME PARAMETER-DEPENDENT MODELS A SUMMARY OF MY RECENT WORK SHENG ZHANG 1. A general formulation of the problems Much of my recent work is on mathematical models
More informationFINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
Proceedings of ALGORITMY pp. 9 3 FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION JAN STEBEL Abstract. The paper deals with the numerical simulations of steady flows
More informationMORTAR MULTISCALE FINITE ELEMENT METHODS FOR STOKES-DARCY FLOWS
MORTAR MULTISCALE FINITE ELEMENT METHODS FOR STOKES-DARCY FLOWS VIVETTE GIRAULT, DANAIL VASSILEV, AND IVAN YOTOV Abstract. We investigate mortar multiscale numerical methods for coupled Stokes and Darcy
More informationCOMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Computational Methods in Elasticity - Michel SALAUN
COMPUAIONAL MEHODS IN ELASICIY Michel SALAUN Départment de Mathématiques et Informatique, ENSICA, oulouse, France Keywords Elastodynamics, Energy minimization, Finite element method, Kirchhoff- Love model,
More informationMIXED FINITE ELEMENTS, COMPATIBILITY CONDITIONS AND APPLICATIONS
D.Boffi F.Brezzi L.Demkowicz R.G.Durán R.S.Falk M.Fortin MIXED FINITE ELEMENTS, COMPATIBILITY CONDITIONS AND APPLICATIONS Lectures given at the C.I.M.E. Summer School held in Cetraro, Italy, June 26-July
More informationA DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT
A DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT HUOYUAN DUAN, SHA LI, ROGER C. E. TAN, AND WEIYING ZHENG Abstract. To deal with the divergence-free
More informationON SINGULAR PERTURBATION OF THE STOKES PROBLEM
NUMERICAL ANALYSIS AND MATHEMATICAL MODELLING BANACH CENTER PUBLICATIONS, VOLUME 9 INSTITUTE OF MATHEMATICS POLISH ACADEMY OF SCIENCES WARSZAWA 994 ON SINGULAR PERTURBATION OF THE STOKES PROBLEM G. M.
More informationSome Remarks on the Reissner Mindlin Plate Model
Some Remarks on the Reissner Mindlin Plate Model Alexandre L. Madureira LNCC Brazil Talk at the Workshop of Numerical Analysis of PDEs LNCC February 10, 2003 1 Outline The 3D Problem and its Modeling Full
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationA NEW FAMILY OF STABLE MIXED FINITE ELEMENTS FOR THE 3D STOKES EQUATIONS
MATHEMATICS OF COMPUTATION Volume 74, Number 250, Pages 543 554 S 0025-5718(04)01711-9 Article electronically published on August 31, 2004 A NEW FAMILY OF STABLE MIXED FINITE ELEMENTS FOR THE 3D STOKES
More informationFinite element approximation on quadrilateral meshes
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2001; 17:805 812 (DOI: 10.1002/cnm.450) Finite element approximation on quadrilateral meshes Douglas N. Arnold 1;, Daniele
More informationWeak Galerkin Finite Element Scheme and Its Applications
Weak Galerkin Finite Element Scheme and Its Applications Ran Zhang Department of Mathematics Jilin University, China IMS, Singapore February 6, 2015 Talk Outline Motivation WG FEMs: Weak Operators + Stabilizer
More informationHybridized Discontinuous Galerkin Methods
Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, 2010 1
More informationON NONCONFORMING LINEAR CONSTANT ELEMENTS FOR SOME VARIANTS OF THE STOKES EQUATIONS
ON NONCONFORMING LINEAR CONSTANT ELEMENTS FOR SOME VARIANTS OF THE STOKES EQUATIONS DOUGLAS N ARNOLD (*) Nota (**) presenta dal s c Franco Brezzi nella seduta del 24 6 93 Sunto L uso di elementi finiti
More informationCheckerboard instabilities in topological shape optimization algorithms
Checkerboard instabilities in topological shape optimization algorithms Eric Bonnetier, François Jouve Abstract Checkerboards instabilities for an algorithm of topological design are studied on a simple
More informationIntroduction to finite element exterior calculus
Introduction to finite element exterior calculus Ragnar Winther CMA, University of Oslo Norway Why finite element exterior calculus? Recall the de Rham complex on the form: R H 1 (Ω) grad H(curl, Ω) curl
More informationChapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems
Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationOn convergence rate of the augmented Lagrangian algorithm for nonsymmetric saddle point problems
Applied Numerical athematics 54 2005) 122 134 wwwelseviercom/locate/apnum On convergence rate of the augmented Lagrangian algorithm for nonsymmetric saddle point problems G Awanou a,,jlai b a Institute
More informationAn introduction to the mathematical theory of finite elements
Master in Seismic Engineering E.T.S.I. Industriales (U.P.M.) Discretization Methods in Engineering An introduction to the mathematical theory of finite elements Ignacio Romero ignacio.romero@upm.es October
More informationVirtual Elements for plates and shells
Virtual Elements for plates and shells Claudia Chinosi Dipartimento di Scienze e Innovazione Tecnologica, Università del Piemonte Orientale, Alessandria, Italia E-mail: claudia.chinosi@unipmn.it eywords:
More informationSCALE INVARIANT FOURIER RESTRICTION TO A HYPERBOLIC SURFACE
SCALE INVARIANT FOURIER RESTRICTION TO A HYPERBOLIC SURFACE BETSY STOVALL Abstract. This result sharpens the bilinear to linear deduction of Lee and Vargas for extension estimates on the hyperbolic paraboloid
More informationACCURATE MODELLING OF STRAIN DISCONTINUITIES IN BEAMS USING AN XFEM APPROACH
VI International Conference on Adaptive Modeling and Simulation ADMOS 213 J. P. Moitinho de Almeida, P. Díez, C. Tiago and N. Parés (Eds) ACCURATE MODELLING OF STRAIN DISCONTINUITIES IN BEAMS USING AN
More informationError estimates for the Raviart-Thomas interpolation under the maximum angle condition
Error estimates for the Raviart-Thomas interpolation under the maximum angle condition Ricardo G. Durán and Ariel L. Lombardi Abstract. The classical error analysis for the Raviart-Thomas interpolation
More information