A Domain Decomposition Algorithm for Contact Problems: Analysis and Implementation
|
|
- Dustin Payne
- 5 years ago
- Views:
Transcription
1 Math. Model. Nat. Phenom. Vol. 4, No. 1, 009, pp A Domain Decomposition Algorithm for Contact Problems: Analysis and Implementation J. Haslinger a, R. Kučera b1 and T. Sassi c a Department of Numerical Mathematics, Charles University Prague, Prague, CZ b Department of Mathematics and Descriptive Geometry, VŠB-TU Ostrava, Ostrava, CZ c Department of Mathematics, University of Basse- Normandie, 1403 Caen, France Abstract. The paper deals with an iterative method for numerical solving frictionless contact problems for two elastic bodies. Each iterative step consists of a Dirichlet problem for the one body, a contact problem for the other one and two Neumann problems to coordinate contact stresses. Convergence is proved by the Banach fixed point theorem in both continuous and discrete case. Numerical experiments indicate scalability of the algorithm for some choices of the relaxation parameter. Key words: contact problem, domain decomposition method AMS subject classification: 74M15, 65Y05 1. Introduction Contact problems take an important place in the computational mechanics (see [8], [10], [5], [1], [13], [14] and references therein). Many numerical procedures have been proposed in engineering literature. The discretization of such problems leads to very large and ill-conditioned systems. Domain decomposition methods represent a possible remedy how to overcome this difficulty. Recently, a Neumann-Dirichlet algorithm for solving frictionless Signorini contact problems between two elastic bodies has been proposed and studied in the continuous and discrete setting [15], [7], []. Bayada, Sabil and Sassi proposed another approach in [3] and proved its convergence in the 1 Corresponding author. radek.kucera@vsb.cz 13
2 continuous setting. Each iterative step of their algorithm is given by a linear elasticity problem with prescribed displacements for one body and a contact problem with a prescribed gap function for the second one. Then two Neumann problems are used to ensure continuity of contact stresses along a common part of the boundary. The purpose of this contribution is to present a modified version of the previous algorithm and to compare the resulting iterative schemes. This method is an extension of the one analyzed in [16] to variational inequalities. The paper is organized as follows. Firstly in Section, we introduce the variational formulation of frictionless Signorini problems. In Section 3 we propose a modified version of the algorithm from [3] (see also [1]). Convergence in the continuous and discrete setting is established in next two sections. Finally, in Section 6 we present implementations of the modified and the original algorithm both based on dual formulations. Results of model examples will be shown in Section 7.. Formulation of the Problem Let us consider two elastic bodies occupying bounded domains Ω α R, α = 1, with sufficiently smooth boundaries, see Fig. 1. The boundary Γ α of Ω α consists of three non-empty disjoint parts Γ α u, Γ α p, Γ α c so that Γ α = Γ α u Γ α p Γ α c. The zero displacements will be prescribed on Γ α u and surface tractions of density p α (L (Γ α p )) act on Γ α p. On the contact interface represented by Γ 1 c and Γ c, we consider the contact conditions: the non-penetration of the bodies and the transmission of forces. Friction will be neglected. The contact conditions are prescribed by a predefined unit critical vector ν (L (Γ 1 c)) (e.g. the outward unit normal vector to Γ 1 c). The initial gap between the bodies in the direction of ν is given by a function g L (Γ 1 c). Let χ : Γ 1 c Γ c, χ(x) = x + g(x)ν(x), x Γ 1 c, be a one-to-one transfer mapping so that we can define ν 1 := ν on Γ 1 c and ν := ν χ 1 on Γ c. Finally, we suppose that the bodies Ω α are subject to volume forces of density f α (L (Ω α )), α = 1,. Γ 1 p Γ 1 u Ω 1 Γ 1 p Γ 1 c Γ c Γ u Ω Γ p Γ p Figure 1: Geometry of the contact problem. 14
3 Before giving the weak formulation of the contact problem, we introduce notation. Let us define the following spaces (α = 1, ): and V α = {v α H 1 (Ω α ) v α = 0 on Γ α u}, V α = V α V α, V α 0 = {v α V α v α ν α = 0 on Γ α c }, V = V 1 V, V 0 = V 1 0 V 0. The standard norm in V α, V α will be denoted by 1,Ω α in what follows. The bilinear form a : V V R representing the inner energy of the system is given by: a(v, w) = a 1 (v 1, w 1 ) + a (v, w ), v, w V, with a α (v α, w α ) = a α ijkle ij (v α )e kl (w α ) dx, α = 1,, Ω α where a α ijkl L (Ω α ) are coefficients of a symmetric, positive definite tensor of the fourth order defining the linear Hook law, v = (v 1, v ), w = (w 1, w ) and e ij (v α ) = ( i vj α + j vi α )/, i, j, k, l = 1,. By b : V R we denote the work of applied forces: b(v) = b 1 (v 1 ) + b (v ), v V, where b α (v α ) = f α v α dx + Ω α Γ α p p α v α ds, α = 1,. The weak formulation of the contact problem reads as follows: } Find u K such that a(u, v u) b(v u) v K, (P) where K is the convex set of kinematically admissible displacements, K = {v V v 1 ν v ν g on Γ 1 c} with v 1 ν := v 1 ν 1 and v ν := (v χ) (ν χ). It is well-known that (P) has a unique solution [13, 11]. The symbol σ(v α ), v α V α, α = 1, stands for the stress field in Ω α : σ ij (v α ) = a α ijkle kl (v α ), i, j = 1,. The normal and tangential components of σ(v α ) on Γ α c are given by: σ ν (v α ) := σ ij (v α )ν α i ν α j, σ t (v α ) := σ ij (v α )ν α j σ ν (v α )ν α. 15
4 3. Algorithm For the sake of simplicity of our presentation we shall suppose that there is no gap between Ω 1 and Ω in the undeformed state. Then Γ 1 c = Γ c =: Γ c, g = 0 and χ = id on Γ c. In what follows the symbol ν denotes the outward unit normal vector on Γ c with respect to Ω 1, while ν α stands for the outward unit normal vector to Ω α on the remaining parts of Γ α. Next we shall suppose that Γ c is of the class C 1,1 and Γ 1 u Γ c = Γ u Γ c, i.e., either dist(γ α u, Γ c ) > 0, α = 1, or Γ 1 u, Γ u meet Γ c at the same points. Let H 1/ (Γ c ) = {ϕ L (Γ c ) v α V α, v α = ϕ on Γ c }, H 1/ (Γ c ) = (H 1/ (Γ c )) (dual of H 1/ (Γ c )), H 1/ (Γ c ) = H 1/ (Γ c ) H 1/ (Γ c ). In view of the assumption on the mutual position of Γ α u and Γ c it is easy to see that the definition of H 1/ (Γ c ) does not depend on the choice of V α, α = 1,. It is also well-known that H 1/ (Γ c ) is the Banach space with the norm ϕ 1 := ϕ ϕ(x) ϕ(y) ds + Γ c Γc ds x y x ds y. (3.1) Γ c In addition, if v α V α then vν α H 1/ (Γ c ) as follows from the smoothness assumption on Γ c. The duality pairing between H 1/ (Γ c ) and H 1/ (Γ c ) will be denoted by, in what follows. In order to split problem (P) into two subproblems coupled through the contact interface, we introduce the mapping that builts functions in Ω α from their values on Γ c. Let P α : H 1/ (Γ c ) V α be the extension operator defined by P α ϕ = v α, ϕ H 1/ (Γ c ), where v α V α, α = 1, satisfies: } a α (v α, z α ) = 0 z α V α 0, (3.) vν α = ϕ on Γ c, or, equivalently, div σ(v α ) = 0 in Ω α, v α = 0 on Γ α u, σ(v α ) ν α = 0 on Γ α p, σ t (v α ) = 0 on Γ c, v α ν = ϕ on Γ c. With any ϕ H 1/ (Γ c ) we associate the closed convex set K (ϕ) = {v V ϕ v ν 0 on Γ c }. Our algorithm is motivated by the following observation. Let λ H 1/ (Γ c ) be given and let u 1, u be solutions of: Find u 1 := u 1 (λ) V 1 such that a 1 (u 1, v 1 ) = b 1 (v 1 ) v 1 V 1 0, (3.4) u 1 ν = λ on Γ c, 16 (3.3)
5 and Find u := u (λ) K (λ) such that a (u, v u ) b (v u ) v K (λ). Then u = (u 1, u ) solves the original contact problem if and only if the contact stresses are continuously transmitted through Γ c (the law of action and reaction). Since σ t (u 1 ) = σ t (u ) = 0 on Γ c as follows from (3.4) and (3.5), one needs only to guarantee that } (3.5) σ ν (u 1 ) = σ ν (u ) on Γ c. (3.6) Suppose that λ is chosen in such a way that w 1 ν = w ν on Γ c, where w 1 V 1, w V solve the following problems: a 1 (w 1, v 1 ) = 1 ( a1 (u 1, v 1 ) + b 1 (v 1 ) a (u, P v 1 ) + b (P v 1 )) v 1 V 1, a (w, v ) = 1 (a (u, v ) b (v ) + a 1 (u 1, P 1 v ) b 1 (P 1 v )) v V. If w 1 ν = w ν on Γ c then a 1 (w 1, w 1 ) + a (w, w ) = 0. Therefore w α = 0 in Ω α, α = 1,, implying (3.6). We propose the following algorithm. ALGORITHM 3.1 Let λ 0 H 1/ (Γ c ) and θ > 0 be given. For k 1 compute u α k, wα k, α = 1,, and λ k as follows: (Step 1) Find u 1 k V1 such that a 1 (u 1 k, v1 ) = b 1 (v 1 ) v 1 V 1 0, (3.7) u 1 k,ν = λ k 1 on Γ c. (Step ) (Step 3) Find w 1 k V1 such that Find u k K (λ k 1 ) such that a (u k, v u k ) b (v u k ) v K (λ k 1 ). a 1 (w 1 k, v1 ) = 1 ( a1 (u 1 k, v1 ) + b 1 (v 1 ) a (u k, P v 1 ) + b (P v 1 )) v 1 V 1. } } (3.8) (3.9) (Step 4) Find w k V such that a (w k, v ) = 1 (a (u k, v ) b (v ) + a 1 (u 1 k, P 1v ) b 1 (P 1 v )) v V. } (3.10) (Step 5) λ k = λ k 1 + θ(w 1 k,ν w k,ν). 17
6 Remark 1. Let us note that our algorithm is the modification of the one analyzed in [, 17], where Bayada, Sabil and Sassi used the extension operator R α : H 1/ (Γ c ) V α defined by R α ϕ = v α, ϕ H 1/ (Γ c ) with v α V α such that: } a α (v α, w α ) = 0 w α V α, w α = 0 on Γ c, v α = ϕ on Γ c. This choice of the extension operator gives a procedure in which continuity not only of the normal but also of the tangential component of the contact stress has to be ensured. 4. Convergence analysis The aim of this section is to prove convergence of ALGORITHM 3.1. We shall use the Banach fixed point theorem for a suitably defined operator. First we reformulate the steps of the previous algorithm in terms of operators acting on Γ c. Let S α : H 1/ (Γ c ) H 1/ (Γ c ) be the following Poincaré-Steklov operator: for any ϕ H 1/ (Γ c ) S α ϕ = σ ν (u α ) on Γ c, α = 1,, (4.1) where u α solves div σ(u α ) = 0 in Ω α, u α = 0 on Γ α u, σ(u α ) ν α = 0 on Γ α p, σ t (u α ) = 0 on Γ c, u α ν = ϕ on Γ c, i.e., u α = P α ϕ, α = 1,. From this and (4.1) we obtain the equivalent definition of S α : S α ϕ, ψ = a α (P α ϕ, v α ) v α V α, v α ν = ( 1) α 1 ψ on Γ c, where the sign of ψ is due to the orientation of ν on Γ c. In order to decouple the influence of the applied and contact forces, we introduce U 1 V 1 0 as the solution of a 1 (U 1, v 1 ) = b 1 (v 1 ) v 1 V 1 0, (4.3) which is the weak formulation of div σ(u 1 ) = f 1 in Ω 1, U 1 = 0 on Γ 1 u, σ(u 1 ) ν 1 = p 1 on Γ 1 p, σ t (U 1 ) = 0 on Γ c, U 1 ν = 0 on Γ c. 18 (4.) (4.4)
7 Thus u 1 k solving (3.7) can be written as: u 1 k = U 1 + P 1 λ k 1. (4.5) The function u k from (3.8) solves the following unilateral problem : div σ(u k ) = f in Ω, u k = 0 on Γ u, σ(u k ) ν = p on Γ p, σ ν (u k ) 0 on Γ c, λ k 1 u k,ν 0 on Γ c, σ ν (u k )(λ k 1 u k,ν ) = 0 on Γ c, σ t (u k ) = 0 on Γ c. (4.6) As before, we decouple u k. Let U be the solution of: div σ(u ) = f in Ω, U = 0 on Γ u, σ(u ) ν = p on Γ p, σ ν (U ) = 0 on Γ c, σ t (U ) = 0 on Γ c. (4.7) Then u k = U + ũ k, (4.8) where ũ k solves: div σ(ũ k ) = 0 in Ω, The weak formulation of (4.9) is given by: ũ k = 0 on Γ u, σ(ũ k ) ν = 0 on Γ p, σ ν (ũ k ) 0 on Γ c, λ k 1 U ν ũ k,ν 0 on Γ c, σ ν (ũ k )(λ k 1 U ν ũ k,ν ) = 0 on Γ c, σ t (ũ k ) = 0 on Γ c. Find ũ k K (λ k 1 U ν ) such that a (ũ k, v ũ k ) 0 v K (λ k 1 U ν ). } (4.9) (4.10) 19
8 Let Q : H 1/ (Γ c ) H 1/ (Γ c ) be defined by: where v K (λ U ν ) solves Then ũ k,ν = Q λ k 1 and Q λ = v ν λ H 1/ (Γ c ), a (v, z v ) 0 z K (λ U ν ). u k = U + P Q λ k 1. (4.11) The solutions to (3.9) and (3.10) can be written by means of the Green formula and the Poincaré-Steklov operators S 1 and S as follows: Therefore S 1 w 1 k,ν, v 1 ν = a 1 (w 1 k, v 1 ) = 1 σ ν(u k) σ ν (u 1 k), v 1 ν, S w k,ν, v ν = a (w k, v ) = 1 σ ν(u k) + σ ν (u 1 k), v ν. w 1 k,ν = 1 S 1 1 (σ ν (u k) σ ν (u 1 k)), w k,ν = 1 S 1 (σ ν (u k) σ ν (u 1 k)). Using these representations we can express Step 5. of the algorithm as follows: λ k = λ k 1 θ (S 1 1 S 1 )(σ ν (u 1 k) σ ν (u k)) = λ k 1 θ M(S 1u 1 k,ν S u k,ν), where M := S 1 1 S 1. Since u 1 k,ν = λ k 1 and u k,ν = U ν + Q λ k 1, it follows that λ k = λ k 1 θ M(S 1λ k 1 S U ν S Q λ k 1 ). (4.1) To simplify notation, let T θ : H 1/ (Γ c ) H 1/ (Γ c ) be defined by: T θ λ := λ θ M(S 1λ S Q λ + g 0 ), λ H 1/ (Γ c ), where g 0 = S U ν. Then (4.1) is equivalent to λ k = T θ λ k 1. (4.13) From the properties of Sα 1, α = 1, formulated in the next lemma it follows that M is a one-toone mapping. Thus if λ H 1/ (Γ c ) is a fixed point of T θ then λ is the searched u 1 ν on Γ c which guarantees the satisfaction of (3.6). To prove the existence (and uniqueness) of the fixed point and convergence of the iterative scheme (4.13), we shall need the following results. 130
9 Lemma. (i) The mappings S α : H 1/ (Γ c ) H 1/ (Γ c ), α = 1, are bounded, bijective, self-adjoint and the bilinear forms generated by S α, Sα 1 are elliptic in the corresponding spaces: c α > 0 : ( 1) α 1 S α ϕ, ϕ c α ϕ 1 d α > 0 : ( 1) α 1 µ, S 1 α µ d α µ 1 ϕ H 1/ (Γ c ), (4.14) µ H 1/ (Γ c ). (4.15) (ii) The mapping S Q : H 1/ (Γ c ) H 1/ (Γ c ) is Lipschitz continuous and monotone: L > 0 : S Q ψ S Q ϕ 1 L ψ ϕ 1 ψ, ϕ H 1/ (Γ c ), (4.16) S Q ψ S Q ϕ, ψ ϕ 0 ψ, ϕ H 1/ (Γ c ). (4.17) Proof. (i) are the standard properties of the Poincaré-Steklov operator. For the proof of (ii) we refer to [7] and also to the next section, where it will be done in the discrete case. To apply the Banach fixed point theorem, the space H 1/ (Γ c ) will be equipped with the scalar product ψ, ϕ M 1 := M 1 ψ, ϕ and the norm ϕ M 1 := ϕ, ϕ 1/ M. From the definition of M 1 and (4.14), (4.15) it follows that M 1 is an equivalent norm to 1 in H 1/ (Γ c ). Using this new norm, (4.14) and (4.16) read as follows: c α > 0 : ( 1) α 1 S α ϕ, ϕ c α ϕ M 1 ϕ H1/ (Γ c ), (4.18) L > 0 : S Q ψ S Q ϕ 1 L ψ ϕ M 1 ψ, ϕ H 1/ (Γ c ). (4.19) Finally, the norm of a linear mapping F into H 1/ (Γ c ) will be denoted by F and F if H 1/ (Γ c ) is equipped with the norm 1 and M 1, respectively. Theorem 3. There exists θ 0 > 0 such that the mapping T θ is contractive for any θ (0, θ 0 ). Proof. From the definition of T θ we have: T θ ψ T θ ϕ M = ψ θ 1 ϕ M MS 1(ψ ϕ) θ M MS Q ψ MS Q ϕ M 1 θ MS 1 (ψ ϕ), ψ ϕ M 1 + θ MS Q ψ MS Q ϕ, ψ ϕ M 1 θ MS Q ψ MS Q ϕ, MS 1 (ψ ϕ) M 1. (4.0) We shall estimate each term on the right hand-side of (4.0). It holds: MS 1 (ψ ϕ) M MS 1 ψ ϕ 1 M 1. (4.1) 131
10 From (4.19) we obtain: MS Q ψ MS Q M 1 = S Q ψ S Q ϕ, MS Q ψ MS Q ϕ Further, (4.18) yields: S Q ψ S Q ϕ 1 M(S Q ψ S Q ϕ) 1 M L ψ ϕ M 1. (4.) MS 1 (ψ ϕ), ψ ϕ M 1 = S 1 (ψ ϕ), ψ ϕ c 1 ψ ϕ M 1. (4.3) The last but one term in (4.0) will be neglected since MS Q ψ MS Q ϕ, ψ ϕ M 1 = S Q ψ S Q ϕ, ψ ϕ 0 (4.4) as follows from (4.17). Finally, MS Q ψ MS Q ϕ, MS 1 (ψ ϕ) M 1 = S Q ψ S Q ϕ, MS 1 (ψ ϕ) From (4.1)-(4.5) we see that L MS 1 ψ ϕ M 1. (4.5) T θ ψ T θ ϕ θ M 1 (1 + 4 MS 1 + θ 4 M L + θ L MS 1 θc 1 ) ψ ϕ M 1. Thus θ 0 = 4c 1 /( MS 1 + M L + L MS 1 ) has the required property. 5. Convergence in Discrete Setting In this section we shall suppose that both Ω 1, Ω are polygonal domains with triangulations Th α which are consistent with the decomposition of the boundary Γ α into Γ α u, Γ α p and Γ c, α = 1,. Moreover, let Γ c be represented by a straight segment. We restrict ourselves to the geometrical conforming situation with matching grids on Γ c. On Ω α, α = 1, we construct the finite element space Vh α of piecewise linear functions: V α h = {v α C( Ω α ) v α T P 1 (T ) for all T T α h } V α. The restriction of V α h on the interface Γ c is denoted by W α h : In view of our assumptions on T α h V α h = V h α V α 0 : W α h = V α h Γc. V α h, α = 1,, and V h = V 1 h V h. By Vα 0,h and the decomposition of Γα, Wh 1 = W h =: W h. Finally, let we denote the finite element subspace of V α 0,h = { v α h V α h v α h,ν = 0 on Γ c }. 13
11 Let P α,h : W h V α h, α = 1, be the extension operator defined by: where v α h Vα h solves P α,h ϕ h = v α h, ϕ h W h, a α (v α h, wα h ) = 0 wα h Vα 0,h, v α h,ν = ϕ h on Γ c. We define the discrete Poincaré-Steklov operator S α,h : W h (W h ) (the dual space of W h ) by S α,h ϕ h, ψ h = a α (P α,h ϕ h, v α h) v α h V α h, v α h,ν = ( 1) α 1 ψ h on Γ c. Finally, we introduce the mapping Q,h : W h W h, λ h Q,h λ h = v h,ν, where v h K h (λ h U h,ν ) solves a (v h, z h v h) 0 z h K h(λ h U h,ν). Here U h V h is the finite element solution of (4.7) and for any ϕ h W h K h(ϕ h ) = {v h V h ϕ h v h,ν 0 on Γ c }, λ h W h. To prove convergence in the discrete case, we shall suppose that the extension operators P α,h : W h V α h enjoy the following property: P α,h ϕ h 1,Ω α C ϕ h 1 ϕ h W h with C > 0 independent of h. From this we obtain the following statement. (A) Lemma 4. Let (A) be satisfied. Then S α,h q h, ϕ h C 1 q h 1 ϕ h 1, (5.1) ( 1) α 1 S α,h q h, q h C q h 1 hold for every q h, ϕ h W h, where C 1 > 0 and C > 0 are independent of h, α = 1,. Proof. Let q h, ϕ h W h be given and v α h = P α,hϕ h V α h. From the definition of S α,h and (A) we have: S α,h q h, ϕ h = a α (P α,h q h, P α,h ϕ h ) (5.) a α P α,h q h 1,Ω α P α,h ϕ h 1,Ω α proving (5.1). Further C 1 q h 1 ϕ h 1 (5.3) S α,h ϕ h, ϕ h = a α (P α,h ϕ h, ( 1) α 1 P α,h ϕ h ) ϕ h W h using that ( 1) α 1 (P α,h ϕ h ) ν = ( 1) α 1 ϕ h on Γ c. From this we arrive at (5.). 133
12 Remark 5. Analogous results hold for the operators S 1 α,h, α = 1,. Theorem 6. The discrete operator S,h Q,h : W h W h is Lipschitz and monotone: l > 0 : S,h Q,h ϕ h S,h Q,h ψ h 1 l ϕ h ψ h 1, (5.4) S,h Q,h ϕ h S,h Q,h ψ h, ϕ h ψ h 0 (5.5) hold for every ϕ h, ψ h W h. In addition, if (A) is satisfied then l in (5.4) does not depend on h. Proof. Let ϕ h, ψ h W h be given and u h K h (ϕ h Uh,ν ), v h K h (ψ h Uh,ν ) be the unique solutions of a (u h, zh u h) 0 zh K h(ϕ h Uh,ν) (5.6) and a (v h, z h v h) 0 z h K h(ψ h U h,ν). (5.7) It is easily seen that the convex sets K h (ϕ h Uh,ν ) and K h (ψ h Uh,ν ) can be written as K h(ϕ h U h,ν) = Φ h U h + K h(0), K h(ψ h U h,ν) = Ψ h U h + K h(0), respectively, where Φ h = P,h ϕ h, Ψ h = P,h ψ h. In other words: any zh K h (ϕ h Uh,ν ) can be written in the form zh = Φ h Uh + z h (5.8) for an appropriate z h K h (0) (analogously for elements of K h (ψ h Uh,ν )). Hence the solutions u h, v h of (5.6) and (5.7), respectively, can be expressed as u h = Φ h Uh + ū h, } vh = Ψ h Uh + v h, (5.9) where ū h, v h K h (0) are the solutions of a (ū h, z h ū h) a (U h Φ h, z h ū h), (5.10) a ( v h, z h v h) a (U h Ψ h, z h v h) (5.11) for every z h K h (0). Inserting z h := v h into (5.10) and z h := ū h into (5.11) and summing both inequalities, we obtain γ ū h v h 1,Ω a (ū h v h, ū h v h) a (Ψ h Φ h, ū h v h) c Ψ h Φ h 1,Ω ū h v h 1,Ω. 134
13 From this and (5.9) we have: u h v h 1,Ω c Ψ h Φ h 1,Ω, where c := const. > 0 does not depend on h. The definition of Q,h and the trace theorem yield: Q,h ϕ h Q,h ψ h 1 = u h,ν v h,ν 1 c u h v h 1,Ω with c > 0 independent of h, again. Finally, c Ψ h Φ h 1,Ω = c P,h (ψ h ϕ h ) 1,Ω S,h Q,h ϕ h S,h Q,h ψ h 1 S,h Q,h ϕ h Q,h ψ h 1 l ϕ h ψ h 1 c S 1,h P,h (ϕ h ψ h ) 1,Ω for some l > 0. If (A) is satisfied then l does not depend on h as follows from Lemma 4. To prove (5.5) let us observe first that from (5.10) it follows: a (u h, Φ h U h u h) = 0, a (u h, z h) 0 z h K h(0) (5.1) (analogously for v h ). Let ϕ h, ψ h W h be arbitrary. Then S,h Q,h ϕ h S,h Q,h ψ h, ϕ h ψ h = a (P,h Q,h ϕ h P,h Q,h ψ h, P,h ψ h P,h ϕ h ) = a (u h v h, Ψ h Φ h ), (5.13) where u h, v h are the solutions of (5.6) and (5.7), respectively. The last term in (5.13) is nonpositive. Indeed, a (u h v h, Ψ h Φ h ) = a (u h v h, Ψ h v h U h) + a (u h v h, U h Φ h + u h) + a (u h v h, v h u h) a (u h, Ψ h v h U h) a (v h, U h Φ h + u h) = a (u h, v h) a (v h, ū h) 0, where ū h, v h K h (0) are defined by (5.9). Here we used (5.1) and the ellipticity of the bilinear form a. Let T θ,h : W h W h be the mapping defined by T θ,h λ h := λ h θ M h(s 1,h λ h S,h Q,h λ h + g 0,h ), where M h := S 1 1,h S 1,h, g 0,h := S,h Uh,ν discretization of T θ introduced in Section 4. and θ > 0. This mapping can be viewed to be a 135
14 Theorem 7. Let (A) be satisfied. Then there exists θ 0 > 0 which does not depend on h such that T θ,h is contractive for any θ (0, θ 0 ). Proof. It can be done in the same way as the one of Theorem 3 using (5.4) and (5.5). Corollary 8. There exists a unique fixed point λ h of T θ,h in W h, θ (0, θ 0 ), and the method of successive approximations is convergent for any choice of the initial approximation λ h,0 W h. Let us comment on the satisfaction of (A) (see also [4]). Recall that the extension mapping P α : H 1/ (Γ c ) W α, α = 1, is defined by (3.) and P α,h : W h V α h by the finite element approximation of (3.): ϕ h W h given, the extension u α h := P α,hϕ h V α h solves a α (u α h, wα h ) = 0 wα h Vα 0,h, } u α h,ν = ϕ (5.14) h on Γ c. Next we shall suppose that the following assumptions are satisfied: (i) (regularity assumption) there exists ε (0, 1/) such that for any ϕ H 1/+ε (Γ c ) the solution u α of (3.) belongs to (H 1+ε (Ω α )) and (ii) (inverse inequality in W h ) c := c(ε) > 0 : u α 1+ε,Ω α c ϕ 1/+ε,Γc ϕ H 1/+ε (Γ c ); c > 0 : ϕ h 1/+ε ch ε ϕ h 1/ ϕ h W h h > 0. Proposition 9. Under (i) and (ii), the assumption (A) is satisfied. Proof. Let ϕ h W h be given and u α h = P α,hϕ h. Then u α h 1,Ω α u α 1,Ω α + u α h u α 1,Ω α c{ ϕ h 1/ + where c > 0 does not depend on h > 0 and Ṽ α h = {w α h V α h w α h,ν = ϕ h on Γ c }. inf w α h e V α h u α w α h 1,Ω α}, (5.15) Denote rh α : (H1+ε (Ω α )) V α : V α h the piecewise linear Lagrange interpolation operator. Since rh αuα Ṽα h, we can use wα h := rα h uα on the right of (5.15). Then u α r h u α h 1,Ω α ch ε u 1+ε,Ω α ch ε ϕ h 1/+ε. 136
15 6. Implementation In this section we describe in more details the finite element approximation of the problem and its implementation. We start with the algebraic formulation of the original non-decomposed problem. Let n α := dim V α 0,h, α = 1, and n := n 1 + n. Further, let m be the number of the contact nodes of Ω 1, i.e., the nodes of Th 1 lying on Γ 1 c \ Γ 1 u. As we consider the matching grids, the contact nodes of Ω 1 and Ω coincide when Γ 1 c = Γ c. If Γ 1 c Γ c, the contact nodes of Ω 1 are bijectively mapped onto the ones of Ω by the transfer mapping χ as mentioned in Section. By A R n n and b R n we denote the stiffness matrix and the load vector, respectively, of the whole structure. Let us note that A, b can be naturally decomposed into blocks corresponding to Ω 1 and Ω so that A = diag(a 1, A ), b = (b 1, b ), where A α R n α n α are symmetric, positive definite and b α R n α, α = 1,. We introduce the (Dirichlet) trace matrices B α R m n α, α = 1, having only one non-zero component in each row, namely 1 on the position of a contact displacement component. Moreover, we denote by N R m m the matrix projecting contact displacements to the direction of ν, i.e., each row of N has at most two non-zero entries given by the coordinates of ν in the positions of contact displacement components. Let us note that B 1, B and N have full row-ranks. Finally, suppose that g R m collects the values of the initial gap at the contact nodes. The finite element approximation of (P) leads to the following quadratic programming problem: minimize 1 u Au u b (6.1) subject to NB 1 u 1 NB u g, where u = (u 1, u ), u α R n α, α = 1,. The problem (6.1) can be solved by the following discrete counterpart of ALGORITHM 3.1. ALGORITHM 6.1 Let λ (0) R m and θ > 0 be given. For k 1 compute u (k) α, w (k) α R nα, α = 1, and λ (k) R m as follows: (Step 1) u (k) 1 := argmin 1 u 1 A 1 u 1 u 1 b 1 subject to NB 1 u 1 = λ (k 1) ; (Step ) u (k) := argmin 1 u A u u b subject to λ (k 1) NB u g; (Step 3) Solve A 1 w (k) 1 = 1 B 1 (B 1 (b 1 A 1 u (k) 1 ) B (A u (k) b )); (Step 4) Solve A w (k) = 1 B (B 1 (b 1 A 1 u (k) 1 ) B (A u (k) b )); (Step 5) λ (k) = λ (k 1) + θ(nb 1 w (k) 1 NB w (k) ). We shall show that the dual formulation simplifies considerably the implementation of the algorithm. The minimization problem in Step 1 is equivalent to solving the following saddle-point system: ( A1 B 1 N NB 1 O ) ( u1 s 1 ) = ( b1 λ (k 1) ). 137
16 Eliminating the first unknown u 1, we obtain: C 1 s 1 = d 1, (6.) where C 1 := NB 1 A 1 1 B 1 N and d 1 := NB 1 A 1 1 b 1 λ (k 1). Since C 1 is symmetric, positive definite, the system (6.) can be solved by the conjugate gradient method [9]. The constrained minimization in Step is equivalent to the saddle-point problem: Find (u, s ) R n R m + : L(u, s ) = min v R n max L(v, r ) = max r R m + r R m + min L(v, r ), v R n where L : R n R m + R is the Lagrangian defined by L(v, r ) := 1 v A v v b + r (λ (k 1) NB v g) and R m + := {r R m : r 0}. The max-min problem leads to the quadratic programming problem with the simple constraints: minimize 1 s C s s d subject to s 0, (6.3) where C := NB A 1 B N is again symmetric, positive definite and d := λ (k 1) NB A 1 b g. To solve (6.3) one can use the conjugate gradient method combined with the projected gradient technique [6]. Let s (k) 1, s (k) be the solutions to (6.), (6.3), respectively. Then the linear systems in Step 3 and Step 4 read as follows: A 1 w (k) 1 = 1 B 1 N (s (k) 1 s (k) ) and A w (k) = 1 B N (s (k) 1 s (k) ). Let us note that the actions of A 1 1 and A 1 are evaluated by the backward substitutions based on the Cholesky factorization of A 1 and A, respectively. Let us denote the relative precision of the k-th iterative step of ALGORITHM 6.1 by ε (k) λ := λ (k) λ (k 1) R m/ λ (k) R m, where R m is the Euclidean norm in R m. We terminate if ε (k) λ tol for a prescribed tolerance tol > 0. In order to increase the efficiency of the algorithm, we initialize the inner iterative solvers in Step 1 and Step by the respective results from the previous outer iterate, i.e. by s (k 1) 1 and s (k 1), and we terminate them by an adaptive (inner) terminating tolerance tol (k) in > 0. The idea is to choose tol (k) in in the outer loop: in such a way that it respects the precision ε (k 1) λ tol (k) in := r tol ε (k 1) λ, 138 achieved
17 where 0 < r tol < 1, ε (0) λ := 1. As we shall compare ALGORITHM 6.1 with the original one proposed in [3], we recall its implementation. ALGORITHM 6. Let λ (0) R m and θ > 0 be given. For k 1 compute u (k) α, w (k) α R nα, α = 1, and λ (k) R m as follows: (Step 1) u (k) 1 := argmin 1 u 1 A 1 u 1 u 1 b 1 subject to B 1 u 1 = λ (k 1) ; (Step ) u (k) := argmin 1 u A u u b subject to Nλ (k 1) NB u g; (Step 3) Solve A 1 w (k) 1 = 1 B 1 (B 1 (b 1 A 1 u (k) 1 ) B (A u (k) b )); (Step 4) Solve A w (k) = 1 B (B 1 (b 1 A 1 u (k) 1 ) B (A u (k) b )); (Step 5) λ (k) = λ (k 1) + θ(b 1 w (k) 1 B w (k) ). At the first glance ALGORITHM 6.1 and ALGORITHM 6. look similar, however, the second one uses vectors from R m. Consequently, the solution of Step 1 of ALGORITHM 6. based on the dual formulation requires to solve two times larger linear system. 7. Numerical Experiments In this section, we shall compare performances of both ALGORITHM 6.1 and ALGORITHM 6. for various values of θ and degrees of freedom n and m. In tables below we report the computational time (time), the number of the outer iterations (out:=final k), the number of the inner iteration (inn that is the total number of the matrix-vector multiplications by C 1 and C ), the final precision (final k) of the contact displacements ε λ := ε λ, the relative continuity of the normal contact stresses (n stress) and the relative continuity of the tangential contact stresses (t stress). The total efficiency of the algorithms is assessed by the ratio eff := inn/m. We set tol = 10 5, r tol = 0.5 and stop computations if the number of the outer iterations is greater than max it = 100. All computations are performed in Matlab 7.1 on Intel(R)Core(TM) Duo CPU, GHz with 1 GB RAM. The quantity inn supplies an information on the computational time as it is proportional to inn. The ratio eff represents a comparison of our algorithms with the realization of similar linear problems by the standard conjugate gradient method. It is well-known that the number of conjugate gradient iterations (i.e., the number of matrix-vector multiplications) is bounded by the size of the problem. One can say that our algorithms exhibit the complexity comparable with the conjugate gradient method when eff is less than two. Let us note that the geometry of our model problems is similar to that one depicted in Figure 1 (Γ 1 c coincides with Γ c in Example 1). All experiments use linear finite element approximation of (P). The triangulations of rectangular domains are defined by cutting rectangles obtained by an equidistant rectangulation. If Γ 1 c is a segment of a circle we define the triangulation of Ω 1 by adapting the one introduced for an appropriate rectangular domain (see Figure 3). 139
18 Example 1. Let us consider the plane elastic bodies Ω 1 = (0, 3) (1, ) and Ω = (0, 3) (0, 1) made of an isotropic, homogeneous material characterized by the Young modulus and the Poisson ratio 0.77 (steal). The decompositions of Γ 1 and Γ are as follows: Γ 1 u = {0} (1, ), Γ 1 c = (0, 3) {1}, Γ 1 p = Γ 1 \ Γ 1 u Γ 1 c, Γ u = {0} (0, 1), Γ c = (0, 3) {1}, Γ p = Γ \ Γ u Γ c. The volume forces vanish for both bodies. The non-vanishing surface tractions p 1 = (p 1 1, p 1 ) act on Γ 1 p so that p 1 1(s, ) = 0, p 1 (s, ) = p 1,L + p 1,R s, s (0, 3), p 1 1(3, s) = 0, p 1 (3, s) = p 1,B( s) + p 1,U (s 1), s (1, ), where p 1,L = 7 107, p 1,R = 1/3 107, p 1,B = and p 1,U = (a) Triangulation of bodies and applied tractions (b) Deformation Ω (c) Contact zone Ω 1 3 x Γ c (d) Contact stress Γ c Figure : Results for n = 864, m = 48. Figure shows results of computations. Comparisons of both algorithms for the best choices of θ are reported in Tables 1 and. It is seen that the performance of ALGORITHM 6.1 is considerably 140
19 higher in all cases and, moreover, the continuity of the tangential stresses is comparable to the computer precision. As inn practically does not depend on m, the algorithms exhibit the scalability property. In Tables 3 and 4 we demonstrate robustness with respect to θ. The situation is better again for ALGORITHM 6.1 (the symbol means that the terminating tolerance is not achieved in max it iterations). Table 1: ALGORITHM 6.1 for θ = 0.4. n m time out/inn ε λ n stress t stress eff /13 6.4e 6 3.7e 5 1.6e / e 6 7.e 5 3.9e / e 6 4.8e 5 4.9e / e 6 4.5e 5 5.0e / e 6 5.8e 5 6.4e /193.5e 6 1.8e 5 6.8e / e 6.4e 6 6.8e Table : ALGORITHM 6. for θ = 0.. n m time out/inn ε λ n stress t stress eff / e 6 3.8e 5.5e / e 6 4.9e 5 3.4e / e 6 6.8e 5 3.7e / e 6 5.8e 5 4.e / e 6 5.e 5 3.e /39 7.7e 6 3.4e 5.5e /43 9.4e 6 1.5e 5 1.0e Example. Let us modify the geometry of the body Ω 1 from the previous example so that Ω 1 = {(x 1, x ) R : 0 < x 1 < 3, 1 + c(x 1 ) < x < }, Γ 1 c = {(x 1, x ) R : 0 < x 1 < 3, 1 + c(x 1 ) = x }, where c is the segment of the circle defined by c(x 1 ) = x 1. The definition of Γ 1 u, Γ 1 p and the geometry of Ω are the same as in Example 1. The surface tractions p 1 = (p 1 1, p 1 ) acting on Γ 1 p are defined by p 1 1(s, ) = 0, p 1 (s, ) = p 1,L + p 1,R s, s (0, 3), p 1 1(3, s) = 0, p 1 (3, s) = p 1,B ( s) + p 1,U (s 1 c(3)), s (1 + c(3), ), 141
20 Table 3: ALGORITHM 6.1, out/inn for various θ. n m θ = 0.1 θ = 0. θ = 0.3 θ = 0.4 θ = 0.5 θ = 0.6 θ = /307 3/166 14/114 15/13 19/156 18/145 37/ /336 3/41 /178 15/135 17/155 1/187 39/ /33 31/49 3/199 15/148 19/179 3/05 51/ /381 9/47 3/08 17/167 18/176 31/79 57/ /383 9/51 4/31 16/168 18/187 4/35 51/ /391 30/76 3/31 18/193 18/186 5/57 63/ /388 30/83 5/57 16/181 17/181 5/57 61/768 Table 4: ALGORITHM 6., out/inn for various θ. n m θ = 0.1 θ = 0. θ = 0.3 θ = 0.4 θ = /550 41/301 7/04 3/184 44/ /590 44/338 33/ /60 44/350 45/ /613 44/364 58/ /649 48/414 74/ /646 45/39 6/ /645 47/43 64/570 Table 5: ALGORITHM 6.1 for θ = 0.3. n m time out/inn ε λ n stress t stress eff /10 9.9e 6 4.0e 5 1.3e /56 7.5e 6.1e 5 6.9e / e 6 7.4e 5 6.1e / e 6 5.e 5 1.1e /41 7.6e 6.4e 5 9.3e / e 6.e 5 8.e / e 6 9.9e 6 1.0e
21 where p 1,L = 7 108, p 1,R = 1/3 107, p 1,B = 4/(1 c(3)) 106 and p 1,U = 4/(1 c(3)) The results of computations are shown in Figure 3. In Tables 5-8 we present analogous comparisons as in Example 1 (with max it=00). The conclusions are similar but the superiority of ALGORITHM 6.1 is now less expressive. The small number of iterations for n = 4, 7 and some θ corresponds to the wrong (zero) approximation of the contact stress on the coarsest meshes. (a) Triangulation of bodies and applied tractions (b) Deformation 1.01 (c) Contact zone 1.5 x 108 (d) Contact stress Ω Ω Γ c 1 Γ c Figure 3: Results for n = 864, m = Conclusions The paper analyzes the new variant of the domain decomposition algorithm for solving frictionless contact problems. Convergence is proved in the continuous as well as discrete setting using the Banach fixed point theorem and, in addition, it is shown its independence on the discretization parameter h. The algorithm can be conceived as a modification of the one proposed in [3]. The main difference consists in the fact that the tangential contact stress is not treated iteratively as it vanishes a-priori. Consequently, the smaller subproblems have to be solved during the iterative process and the algorithm exhibits a better behavior than the original one. The analysis of the paper can be extended to contact problems with friction. 143
22 Table 6: ALGORITHM 6. for θ = 0.3. n m time out/inn ε λ n stress t stress eff / e 6.6e 5.5e /46 9.7e 6 1.1e 5 1.1e / e 6 5.3e 6 4.8e / e 6 1.7e 5 1.1e / e 6 3.8e 5 5.7e / e 6.7e 5 5.e / e 6 5.3e 5 1.7e Table 7: ALGORITHM 6.1, out/inn for various θ. n m θ = 0.1 θ = 0. θ = 0.3 θ = 0.4 θ = 0.5 θ = 0.6 θ = /614 48/3 31/10 4/164 18/13 15/10 1/ /693 47/37 31/56 3/190 18/15 14/119 1/ /685 50/394 41/375 49/469 89/ /751 46/404 38/361 48/49 81/ / /776 47/479 40/41 48/510 73/ / / /837 46/489 38/433 66/778 79/910 1/ / /890 44/505 51/634 9/ / /1618 1/1499 Table 8: ALGORITHM 6., out/inn for various θ. n m θ = 0.1 θ = 0. θ = 0.3 θ = 0.4 θ = / /58 45/396 34/300 7/ /117 74/653 50/46 38/360 31/ /1136 7/571 50/441 37/ / /701 51/513 44/ /139 75/796 5/583 63/ / /899 5/635 66/ /1579 7/99 50/680 56/
23 Acknowledgements Supported by the grant 01/07/094 of the Grant Agency of the Czech Republic and by the Research Project MSM of the Czech Ministry of Education. References [1] L. Baillet, T. Sassi. Simulations numériques de différentes méthodes d éments finis pour les problémes contact avec frottement. C. R. Acad. Sci, Paris, Ser. IIB, 331 (003), [] G. Bayada, J. Sabil, T. Sassi. Algorithme de décomposition de domaine pour un probléme de Signorini sans frottement. C. R. Acad. Sci. Paris, Ser. I335 (00), [3] G. Bayada, J. Sabil, T. Sassi. A Neumann-Neumann domain decomposition algorithm for the Signorini problem. Appl. Math. Letters, 17 (004), [4] P. E. Bjorstad, O. B. Widlund. Iterative methods for the solution of elliptic problems on regions partitioned into substructures. SIAM J. Numerical Analysis, 3 (1986), No. 6, [5] P. W. Christensen, A. Klarbring, J. S. Pang, N. Strömberg. Formulation and comparison of algorithms for frictional contact problems. Internat. J. Numer. Methods Engrg., 4 (1998), No. 1, [6] Z. Dostál, J. Schöberl. Minimizing quadratic functions over non-negative cone with the rate of convergence and finite termination. Comput. Optim. Appl., 30 (005), No. 1, [7] C. Eck, B. Wohlmuth. Convergence of a Contact-Neumann iteration for the solution of twobody contact problems. Mathematical Models and Methods in Applied Sciences, 13 (003), No. 8, [8] R. Glowinski, J. L. Lions, R. Trémoliére. Numerical analysis of variational inequalities. Studies in Mathematics and its Applications, Volume VIII, North-Holland, Amsterdam, [9] G. H. Golub, C. F. Van Loan. Matrix computation. The Johns Hopkins University Press, Baltimore, [10] J. Haslinger, Z. Dostál, R. Kučera. On a splitting type algorithm for the numerical realization of contact problems with Coulomb friction. Comput. Methods Appl. Mech. Engrg., 191 (00), No. 1-, [11] J. Haslinger, I. Hlaváček, J. Nečas. Numerical methods for unilateral problems in solid mechanics. Handbook of Numerical Analysis, Volume IV, Part, North Holland, Amsterdam,
24 [1] M. A. Ipopa. Algorithmes de Décomposition de Domaine pour les problémes de Contact: Convergence et simulations numériques. Thesis, Universit de Caen, 008. [13] N. Kikuchi, J. T. Oden. Contact problems in elasticity: A study of variational inequalities and finite element methods. SIAM, Philadelphia, [14] R. Kornhuber, R. Krause. Adaptive multigrid methods for Signorini s problem in linear elasticity. Comput. Vis. Sci., 4 (001), No. 1, 9 0., [15] R. Krause, B. Wohlmuth. A Dirichlet-Neumann type algorithm for contact problems with friction. Comput. Vis. Sci., 5 (00), No. 3, [16] P. Le Tallec. Domain decomposition methods in computational mechanics. Comput. Mech. Adv., 1 (1994), No., [17] J. Sabil. Modélisation et méthodes de décomposition de domaine pour des problémes de contact. Thesis, INSA de Lyon,
A domain decomposition algorithm for contact problems with Coulomb s friction
A domain decomposition algorithm for contact problems with Coulomb s friction J. Haslinger 1,R.Kučera 2, and T. Sassi 1 1 Introduction Contact problems of elasticity are used in many fields of science
More informationA new algorithm for solving 3D contact problems with Coulomb friction
A new algorithm for solving 3D contact problems with Coulomb friction Radek Kučera VŠB-TU Ostrava, Department of Mathematics and Descriptive Geometry, name@email.address Summary. The paper deals with solving
More informationFETI domain decomposition method to solution of contact problems with large displacements
FETI domain decomposition method to solution of contact problems with large displacements Vít Vondrák 1, Zdeněk Dostál 1, Jiří Dobiáš 2, and Svatopluk Pták 2 1 Dept. of Appl. Math., Technical University
More informationShort title: Total FETI. Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ Ostrava, Czech Republic
Short title: Total FETI Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ-70833 Ostrava, Czech Republic mail: zdenek.dostal@vsb.cz fax +420 596 919 597 phone
More informationAn example of nonuniqueness for the continuous static unilateral contact model with Coulomb friction
An example of nonuniqueness for the continuous static unilateral contact model with Coulomb friction Un exemple de non-unicité pour le modèle continu statique de contact unilatéral avec frottement de Coulomb
More informationScalable BETI for Variational Inequalities
Scalable BETI for Variational Inequalities Jiří Bouchala, Zdeněk Dostál and Marie Sadowská Department of Applied Mathematics, Faculty of Electrical Engineering and Computer Science, VŠB-Technical University
More informationA smooth variant of the fictitious domain approach
A smooth variant of the fictitious domain approach J. Haslinger, T. Kozubek, R. Kučera Department of umerical Mathematics, Charles University, Prague Department of Applied Mathematics, VŠB-TU, Ostrava
More informationOn Mixed Methods for Signorini Problems
Annals of University of Craiova, Math. Comp. Sci. Ser. Volume 30, 2003, Pages 45 52 ISS: 1223-6934 On Mixed Methods for Signorini Problems Faker Ben Belgacem, Yves Renard, and Leila Slimane Abstract. We
More informationMATHEMATICAL MODELLING OF HUMAN LIMB IN MATLAB
MATHEMATICAL MODELLING OF HUMAN LIMB IN MATLAB J. Daněk Department of Mathematics, University of West Bohemia Abstract The contribution deals with mathematical modelling of human limb and numerical simulation
More informationConstruction of a New Domain Decomposition Method for the Stokes Equations
Construction of a New Domain Decomposition Method for the Stokes Equations Frédéric Nataf 1 and Gerd Rapin 2 1 CMAP, CNRS; UMR7641, Ecole Polytechnique, 91128 Palaiseau Cedex, France 2 Math. Dep., NAM,
More informationOn finite element uniqueness studies for Coulomb s frictional contact model
On finite element uniqueness studies for Coulomb s frictional contact model Patrick HILD We are interested in the finite element approximation of Coulomb frictional unilateral contact problem in linear
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 informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationCONVERGENCE OF A NEUMANN-DIRICHLET ALGORITHM FOR TWO-BODY CONTACT PROBLEMS WITH NON LOCAL COULOMB S FRICTION LAW
ESAIM: MAN 4 008 43 6 DOI: 10.1051/man:008003 ESAIM: Mathematical Modelling and Numerical Analysis www.esaim-man.org CONVERGENCE OF A NEUMANN-DIRICHLET ALGORITHM FOR TWO-BODY CONTACT PROBLEMS WITH NON
More informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationA uniqueness criterion for the Signorini problem with Coulomb friction
A uniqueness criterion for the Signorini problem with Coulomb friction Yves REARD 1 Abstract he purpose of this paper is to study the solutions to the Signorini problem with Coulomb friction (the so-called
More informationThe Mortar Finite Element Method for Contact Problems
The Mortar Finite Element Method for Contact Problems F. Ben Belgacem, P. Hild, P. Laborde Mathématiques pour l Industrie et la Physique, Unité Mixte de Recherche CNRS UPS INSAT (UMR 5640), Université
More informationAn additive average Schwarz method for the plate bending problem
J. Numer. Math., Vol. 10, No. 2, pp. 109 125 (2002) c VSP 2002 Prepared using jnm.sty [Version: 02.02.2002 v1.2] An additive average Schwarz method for the plate bending problem X. Feng and T. Rahman Abstract
More informationDomain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions
Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions Ernst P. Stephan 1, Matthias Maischak 2, and Thanh Tran 3 1 Institut für Angewandte Mathematik, Leibniz
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY 11, USA Dan Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
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 informationThe All-floating BETI Method: Numerical Results
The All-floating BETI Method: Numerical Results Günther Of Institute of Computational Mathematics, Graz University of Technology, Steyrergasse 30, A-8010 Graz, Austria, of@tugraz.at Summary. The all-floating
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationTwo-scale Dirichlet-Neumann preconditioners for boundary refinements
Two-scale Dirichlet-Neumann preconditioners for boundary refinements Patrice Hauret 1 and Patrick Le Tallec 2 1 Graduate Aeronautical Laboratories, MS 25-45, California Institute of Technology Pasadena,
More informationA Locking-Free MHM Method for Elasticity
Trabalho apresentado no CNMAC, Gramado - RS, 2016. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics A Locking-Free MHM Method for Elasticity Weslley S. Pereira 1 Frédéric
More information1.The anisotropic plate model
Pré-Publicações do Departamento de Matemática Universidade de Coimbra Preprint Number 6 OPTIMAL CONTROL OF PIEZOELECTRIC ANISOTROPIC PLATES ISABEL NARRA FIGUEIREDO AND GEORG STADLER ABSTRACT: This paper
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
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 informationConstrained Minimization and Multigrid
Constrained Minimization and Multigrid C. Gräser (FU Berlin), R. Kornhuber (FU Berlin), and O. Sander (FU Berlin) Workshop on PDE Constrained Optimization Hamburg, March 27-29, 2008 Matheon Outline Successive
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 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 informationIterative Methods for Linear Systems
Iterative Methods for Linear Systems 1. Introduction: Direct solvers versus iterative solvers In many applications we have to solve a linear system Ax = b with A R n n and b R n given. If n is large the
More informationOn Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities
On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
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 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 informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
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 informationEXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 147, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE AND REGULARITY OF SOLUTIONS FOR
More informationWeighted Regularization of Maxwell Equations Computations in Curvilinear Polygons
Weighted Regularization of Maxwell Equations Computations in Curvilinear Polygons Martin Costabel, Monique Dauge, Daniel Martin and Gregory Vial IRMAR, Université de Rennes, Campus de Beaulieu, Rennes,
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 informationMultigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids
Multigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids Long Chen 1, Ricardo H. Nochetto 2, and Chen-Song Zhang 3 1 Department of Mathematics, University of California at Irvine. chenlong@math.uci.edu
More informationVariational Formulations
Chapter 2 Variational Formulations In this chapter we will derive a variational (or weak) formulation of the elliptic boundary value problem (1.4). We will discuss all fundamental theoretical results that
More informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
More informationCIMPA Summer School on Current Research in Finite Element Methods
CIMPA Summer School on Current Research in Finite Element Methods Local zooming techniques for the finite element method Alexei Lozinski Laboratoire de mathématiques de Besançon Université de Franche-Comté,
More informationChapter 3 Conforming Finite Element Methods 3.1 Foundations Ritz-Galerkin Method
Chapter 3 Conforming Finite Element Methods 3.1 Foundations 3.1.1 Ritz-Galerkin Method Let V be a Hilbert space, a(, ) : V V lr a bounded, V-elliptic bilinear form and l : V lr a bounded linear functional.
More informationPenalty Robin-Robin Domain Decomposition Schemes 2 for Contact Problems of Nonlinear Elastic Bodies 3 UNCORRECTED PROOF
1 Penalty Robin-Robin Domain Decomposition Schemes 2 for Contact Problems of Nonlinear Elastic Bodies 3 Ihor I. Prokopyshyn 1, Ivan I. Dyyak 2, Rostyslav M. Martynyak 1,and 4 Ivan A. Prokopyshyn 2 5 1
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 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 informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
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 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 informationTHE solution of the absolute value equation (AVE) of
The nonlinear HSS-like iterative method for absolute value equations Mu-Zheng Zhu Member, IAENG, and Ya-E Qi arxiv:1403.7013v4 [math.na] 2 Jan 2018 Abstract Salkuyeh proposed the Picard-HSS iteration method
More informationON THE ANALYSIS OF A VISCOPLASTIC CONTACT PROBLEM WITH TIME DEPENDENT TRESCA S FRIC- TION LAW
Electron. J. M ath. Phys. Sci. 22, 1, 1,47 71 Electronic Journal of Mathematical and Physical Sciences EJMAPS ISSN: 1538-263X www.ejmaps.org ON THE ANALYSIS OF A VISCOPLASTIC CONTACT PROBLEM WITH TIME
More informationInteraction of Incompressible Fluid and Moving Bodies
WDS'06 Proceedings of Contributed Papers, Part I, 53 58, 2006. ISBN 80-86732-84-3 MATFYZPRESS Interaction of Incompressible Fluid and Moving Bodies M. Růžička Charles University, Faculty of Mathematics
More informationFrom the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes
From the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes Dylan Copeland 1, Ulrich Langer 2, and David Pusch 3 1 Institute of Computational Mathematics,
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY, USA. Dan_Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
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 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 informationTheoretically supported scalable FETI for numerical solution of variational inequalities
Theoretically supported scalable FETI for numerical solution of variational inequalities Zdeněk Dostál and David Horák Abstract The FETI method with a natural coarse grid is combined with recently proposed
More informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationPDE & FEM TERMINOLOGY. BASIC PRINCIPLES OF FEM.
PDE & FEM TERMINOLOGY. BASIC PRINCIPLES OF FEM. Sergey Korotov Basque Center for Applied Mathematics / IKERBASQUE http://www.bcamath.org & http://www.ikerbasque.net 1 Introduction The analytical solution
More informationPhD dissertation defense
Isogeometric mortar methods with applications in contact mechanics PhD dissertation defense Brivadis Ericka Supervisor: Annalisa Buffa Doctor of Philosophy in Computational Mechanics and Advanced Materials,
More informationDiscrete Maximum Principle for a 1D Problem with Piecewise-Constant Coefficients Solved by hp-fem
The University of Texas at El Paso Department of Mathematical Sciences Research Reports Series El Paso, Texas Research Report No. 2006-10 Discrete Maximum Principle for a 1D Problem with Piecewise-Constant
More informationThe Plane Stress Problem
The Plane Stress Problem Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) February 2, 2010 Martin Kronbichler (TDB) The Plane Stress Problem February 2, 2010 1 / 24 Outline
More informationGAKUTO International Series
GAKUTO International Series Mathematical Sciences and Applications, Vol.28(2008) Proceedings of Fourth JSIAM-SIMMAI Seminar on Industrial and Applied Mathematics, pp.139-148 A COMPUTATIONAL APPROACH TO
More informationGeometric Multigrid Methods
Geometric Multigrid Methods Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University IMA Tutorial: Fast Solution Techniques November 28, 2010 Ideas
More informationNUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL
Proceedings of ALGORITMY 212 pp. 29 218 NUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL MARTIN HADRAVA, MILOSLAV FEISTAUER, AND PETR SVÁČEK Abstract. The present paper
More informationTermination criteria for inexact fixed point methods
Termination criteria for inexact fixed point methods Philipp Birken 1 October 1, 2013 1 Institute of Mathematics, University of Kassel, Heinrich-Plett-Str. 40, D-34132 Kassel, Germany Department of Mathematics/Computer
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 informationON SOME ELLIPTIC PROBLEMS IN UNBOUNDED DOMAINS
Chin. Ann. Math.??B(?), 200?, 1 20 DOI: 10.1007/s11401-007-0001-x ON SOME ELLIPTIC PROBLEMS IN UNBOUNDED DOMAINS Michel CHIPOT Abstract We present a method allowing to obtain existence of a solution for
More informationDuality method in limit analysis problem of non-linear elasticity
Duality method in limit analysis problem of non-linear elasticity Igor A. Brigadnov Department of Computer Science North-Western State Technical University Millionnaya 5, St. Petersburg, 191186, Russia
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 informationNUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE
Proceedings of the Czech Japanese Seminar in Applied Mathematics 2005 Kuju Training Center, Oita, Japan, September 15-18, 2005 pp. 69 76 NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING
More informationA POSTERIORI ERROR ESTIMATES OF TRIANGULAR MIXED FINITE ELEMENT METHODS FOR QUADRATIC CONVECTION DIFFUSION OPTIMAL CONTROL PROBLEMS
A POSTERIORI ERROR ESTIMATES OF TRIANGULAR MIXED FINITE ELEMENT METHODS FOR QUADRATIC CONVECTION DIFFUSION OPTIMAL CONTROL PROBLEMS Z. LU Communicated by Gabriela Marinoschi In this paper, we discuss a
More informationEfficient smoothers for all-at-once multigrid methods for Poisson and Stokes control problems
Efficient smoothers for all-at-once multigrid methods for Poisson and Stoes control problems Stefan Taacs stefan.taacs@numa.uni-linz.ac.at, WWW home page: http://www.numa.uni-linz.ac.at/~stefant/j3362/
More informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
More informationA Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract
A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin
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 WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATIONS
A WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATIONS LIN MU, JUNPING WANG, YANQIU WANG, AND XIU YE Abstract. This article introduces and analyzes a weak Galerkin mixed finite element method
More informationProjected Schur Complement Method for Solving Non-Symmetric Saddle-Point Systems (Arising from Fictitious Domain Approaches)
Projected Schur Complement Method for Solving Non-Symmetric Saddle-Point Systems (Arising from Fictitious Domain Approaches) Jaroslav Haslinger, Charles University, Prague Tomáš Kozubek, VŠB TU Ostrava
More informationDomain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions
Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions Bernhard Hientzsch Courant Institute of Mathematical Sciences, New York University, 51 Mercer Street, New
More informationNew constructions of domain decomposition methods for systems of PDEs
New constructions of domain decomposition methods for systems of PDEs Nouvelles constructions de méthodes de décomposition de domaine pour des systèmes d équations aux dérivées partielles V. Dolean?? F.
More informationDomain Decomposition Methods for Mortar Finite Elements
Domain Decomposition Methods for Mortar Finite Elements Dan Stefanica Courant Institute of Mathematical Sciences New York University September 1999 A dissertation in the Department of Mathematics Submitted
More informationA Posteriori Estimates for Cost Functionals of Optimal Control Problems
A Posteriori Estimates for Cost Functionals of Optimal Control Problems Alexandra Gaevskaya, Ronald H.W. Hoppe,2 and Sergey Repin 3 Institute of Mathematics, Universität Augsburg, D-8659 Augsburg, Germany
More informationBoundary Value Problems and Iterative Methods for Linear Systems
Boundary Value Problems and Iterative Methods for Linear Systems 1. Equilibrium Problems 1.1. Abstract setting We want to find a displacement u V. Here V is a complete vector space with a norm v V. In
More informationPreconditioned space-time boundary element methods for the heat equation
W I S S E N T E C H N I K L E I D E N S C H A F T Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
More informationFrom the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes
www.oeaw.ac.at From the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes D. Copeland, U. Langer, D. Pusch RICAM-Report 2008-10 www.ricam.oeaw.ac.at From the Boundary Element
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 informationA Nonoverlapping Subdomain Algorithm with Lagrange Multipliers and its Object Oriented Implementation for Interface Problems
Contemporary Mathematics Volume 8, 998 B 0-88-0988--03030- A Nonoverlapping Subdomain Algorithm with Lagrange Multipliers and its Object Oriented Implementation for Interface Problems Daoqi Yang. Introduction
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 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 informationSome 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 informationThe mortar element method for quasilinear elliptic boundary value problems
The mortar element method for quasilinear elliptic boundary value problems Leszek Marcinkowski 1 Abstract We consider a discretization of quasilinear elliptic boundary value problems by the mortar version
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
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 informationExistence of minimizers for the pure displacement problem in nonlinear elasticity
Existence of minimizers for the pure displacement problem in nonlinear elasticity Cristinel Mardare Université Pierre et Marie Curie - Paris 6, Laboratoire Jacques-Louis Lions, Paris, F-75005 France Abstract
More information