PhD dissertation defense
|
|
- Evelyn Payne
- 5 years ago
- Views:
Transcription
1 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, XXVIII Cycle ( ) Pavia, January 19th, 2017
2 Context of the work This research work is part of a project between the of Pavia and the company which aims to treat complex contact problems with isogeometric analysis. It has been done in collaboration with A. Buffa ( Pavia), G. Elber ( ), M. Martinelli ( Pavia), F. Massarwi ( ), L. Wunderlich ( ) and B. Wohlmuth ( ). Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
3 Outline 1 IGA basics 2 Isogeometric mortar methods 3 Applications in contact mechanics Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
4 IGA basics Univariate B-Splines The parametric unit interval 0 1 Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
5 IGA basics Univariate B-Splines The parametric unit interval is split by the breakpoints ζ j (j = 1,..., N) leading to a parametric mesh U = {0, 1 5, 2 5, 3 5, 4 5, 1} Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
6 IGA basics Univariate B-Splines The parametric unit interval is split by the breakpoints ζ j (j = 1,..., N) leading to a parametric mesh U = {0, 1 5, 2 5, 3 5, 4 5, 1} Each breakpoint ζ j has a multiplicity m j. The breakpoints and their multiplicities define the knot vector Ξ = {ζ 1,..., ζ }{{} 1, ζ 2,..., ζ 2,..., ζ }{{} N,..., ζ N }. }{{} m 1 times m 2 times m N times Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
7 IGA basics Univariate B-Splines The parametric unit interval is split by the breakpoints ζ j (j = 1,..., N) leading to a parametric mesh U = {0, 1 5, 2 5, 3 5, 4 5, 1} Each breakpoint ζ j has a multiplicity m j. The breakpoints and their multiplicities define the knot vector Ξ = {ζ 1,..., ζ }{{} 1, ζ 2,..., ζ 2,..., ζ }{{} N,..., ζ N }. }{{} m 1 times m 2 times m N times Given a degree p and a knot vector Ξ, the univariate (p + 1)-order B-Splines can be defined recursively. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
8 IGA basics Univariate B-Splines Important properties Univariate quadratique B-Spline functions defined on the knot vector Ξ = {0, 0, 0, 1 5, 2 5, 2 5, 3 5, 4 5, 1, 1, 1} We note the C 0 continuity in ζ j = 2 5. Each B p i is a piecewise positive polynomial of degree p and has a local support. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
9 IGA basics Univariate B-Splines Important properties Univariate quadratique B-Spline functions defined on the knot vector Ξ = {0, 0, 0, 1 5, 2 5, 2 5, 3 5, 4 5, 1, 1, 1} We note the C 0 continuity in ζ j = 2 5. Each B p i is a piecewise positive polynomial of degree p and has a local support. On each element, at most (p + 1) functions have non-zero values. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
10 IGA basics Univariate B-Splines Important properties Univariate quadratique B-Spline functions defined on the knot vector Ξ = {0, 0, 0, 1 5, 2 5, 2 5, 3 5, 4 5, 1, 1, 1} We note the C 0 continuity in ζ j = 2 5. Each B p i is a piecewise positive polynomial of degree p and has a local support. On each element, at most (p + 1) functions have non-zero values. The inter-element continuity is defined by the breakpoint multiplicity: the basis is C p m j at ζ j. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
11 IGA basics Multivariate B-Splines and NURBS For any direction δ, we have p δ, U δ, Ξ δ and n δ. Then we introduce U = (U 1 U 2... U d ), Ξ = (Ξ 1 Ξ 2... Ξ d ) and I = {i = (i 1,..., i d ) : 1 i δ n δ } to define multivariate B-Splines by a tensor product of univariate B-Splines. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
12 IGA basics Multivariate B-Splines and NURBS For any direction δ, we have p δ, U δ, Ξ δ and n δ. Then we introduce U = (U 1 U 2... U d ), Ξ = (Ξ 1 Ξ 2... Ξ d ) and I = {i = (i 1,..., i d ) : 1 i δ n δ } to define multivariate B-Splines by a tensor product of univariate B-Splines. Given a set of positive weights {ω i, i I}, the Non-Uniform Rational B-Splines functions are defined as rational functions of multivariate B-Spline functions by N p i (ζ) = ω i B p i (ζ) i I ω i B p i (ζ). Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
13 IGA basics Multivariate B-Splines and NURBS For any direction δ, we have p δ, U δ, Ξ δ and n δ. Then we introduce U = (U 1 U 2... U d ), Ξ = (Ξ 1 Ξ 2... Ξ d ) and I = {i = (i 1,..., i d ) : 1 i δ n δ } to define multivariate B-Splines by a tensor product of univariate B-Splines. Given a set of positive weights {ω i, i I}, the Non-Uniform Rational B-Splines functions are defined as rational functions of multivariate B-Spline functions by N p i (ζ) = ω i B p i (ζ) i I ω i B p i (ζ). S p (Ξ) is the spline space spanned by the functions Ŝ p i (ζ), N p (Ξ) the NURBS space spanned by the functions N p i (ζ). Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
14 IGA basics Geometrical parametrisation In CAGD, NURBS and splines are widely spread as they are capable of accurately describing various geometries. Given a set of Control Points C i R d, i I, a NURBS parametrization of a curve (d=1), a surface (d = 2) or a solid (d = 3) is given by F : Ω Ω F(ζ) = i I C i Np i (ζ). F Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
15 IGA basics Geometrical approximation The computational domain Ω is generally decomposed into K non-overlapping subdomains, i.e., K Ω = Ω k, Ω i Ω j =, i j. k=1 Some of the reasons: Geometrical definition purposes Different materials Figure: A tire kindly given by Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
16 IGA basics Some model problems Strong formulations Static mechanical balance problem div(σ) = f in Ω, u = u D on Γ D, σ n = l on Γ N. Laplace problem u = f in Ω, u = u D on Γ D, u n = l on Γ N. under the following assumptions: linear elasticity ε = 1 2 ( u + ( u)t ) = s u, small displacement-deformation σ = λ Lamé tr(ε) I + 2 µ Lamé ε. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
17 IGA basics Some model problems Weak formulations Find u V D such that a(u, v) = f (v), v V 0. Laplace problem a(u, v) = u v dx, Ω f (v) = f v dx + l v ds. Ω Γ N Static mechanical balance problem a(u, v) = λ Lamé div(u) div(v) + Ω 2 µ Lamé s u : s v dx, f (v) = f v dx + l v ds. Ω Γ N Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
18 IGA basics Motivations Each Ω k is parametrised by a NURBS parametrization F k. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
19 IGA basics Motivations Each Ω k is parametrised by a NURBS parametrization F k. Each subdomain is discretized independently non-conforming meshes at the subdomain interfaces! Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
20 IGA basics Motivations Each Ω k is parametrised by a NURBS parametrization F k. Each subdomain is discretized independently non-conforming meshes at the subdomain interfaces! The approximation space on Ω k is V k,h = {v k = v k F 1 k, v k (N p k (Ξ k )) d P }, i.e., a push-forward of a NURBS space. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
21 IGA basics Motivations Each Ω k is parametrised by a NURBS parametrization F k. Each subdomain is discretized independently non-conforming meshes at the subdomain interfaces! The approximation space on Ω k is V k,h = {v k = v k F 1 k, v k (N p k (Ξ k )) d P }, i.e., a push-forward of a NURBS space. at the interfaces? Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
22 Where we are 1 IGA basics 2 Isogeometric mortar methods 3 Applications in contact mechanics Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
23 Isogeometric mortar methods The discrete mortar problem We define the interface γ k1 k 2 by γ k1 k 2 = Ω k1 Ω k2 k 1, k 2 such that 1 k 1 < k 2 K. For each interface, one of the adjacent subdomains is chosen as the master side and the other one as the slave side. The skeleton Γ is the union of all the interfaces. Ω 1 γ 2 Ω 1 γ1 Ω M(2) Ω S(2) Ω S(1) n Ω 1 3 Ω 2 γ 1 Ω S(1) n 2 Ω M(1) n 1 Ω 2 γ 2 Ω S(2) n 2 Ω M(1) Ω M(2) Ω 3 Figure: Geometrical conforming case (left) and slave conforming case (right) Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
24 Isogeometric mortar methods The discrete mortar problem Mortar method weak continuity imposition of u at each interface. Discrete mortar saddle point formulation Find (u h, λ h ) V h M h, such that with V h = a(u h, v h ) + b(v h, λ h ) = f (v h ), v h V h, b(u h, φ h ) = 0, φ h M h, K V k,h and b(u h, φ h ) = k=1 Γ u h φ h ds. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
25 Isogeometric mortar methods The discrete mortar problem Mortar method weak continuity imposition of u at each interface. Discrete mortar saddle point formulation Find (u h, λ h ) V h M h, such that with V h = a(u h, v h ) + b(v h, λ h ) = f (v h ), v h V h, b(u h, φ h ) = 0, φ h M h, K V k,h and b(u h, φ h ) = k=1 Γ u h φ h ds. Choose properly M h Two requirements on each interface to have an optimal method an appropriate approximation order for the multiplier space a uniform inf-sup stability between the primal and multiplier spaces Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
26 Where we are 1 IGA basics 2 Isogeometric mortar methods The discrete mortar problem Lagrange multiplier spaces A key issue: the evaluation of the mortar integrals 3 Applications in contact mechanics Contact problem definition Example of treatment of deformable-deformable contact in a Lagrange multiplier formulation Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
27 Isogeometric mortar methods Lagrange multiplier spaces M h = L M l,h where M l,h is the multiplier space on an interface. l=1 From now on, we focus on the one interface case. M h is generated from a parametric multiplier space M h. We consider the following choices for this latter one: choice 1: Mh S p S ( Γ), choice 2: Mh S p S 1 ( Γ), choice 3: Mh S p S 2 ( Γ). Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
28 Isogeometric mortar methods Lagrange multiplier spaces Choice 1 M h S p S ( Γ) A boundary modification can be necessary. 2 1 B 2 i h 2h 3h ζ Figure: Left boundary modification of B-Spline functions of degree p = 3 A degree reduction in the boundary elements. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
29 Isogeometric mortar methods Lagrange multiplier spaces Choice 2 M h S p S 1 ( Γ) Chapelle-Bathe 1 tests to estimate the stability constant 10 0 Stability constant P8 P7 P7 P6 P6 P5 P5 P4 P4 P3 P3 P2 P2 P1 P1 P0 O(E 1 ) Inverse of the number of elements Unstable pairings 1 The inf-sup test. D. Chapelle, K. J. Bathe. Computers & Structures, Vol. 47 (1993), No. 4/5, pp Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
30 Isogeometric mortar methods Lagrange multiplier spaces Choice 3 M h S p S 2 ( Γ) Proof for the p S /p S 2 pairing in 2D inf φ S p S 2 sup û S p S 0 Γ φ û dŝ û L 2 φ L 2 α Regularity of the mapping Quasi-uniformity of the meshes Choose û S p S 0 such that 2 xxû = φ D : S p S 0 S p S 1 \ R and D : S p S 1 \ R S p S 2 are bijections Work with Sobolev norms Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
31 Isogeometric mortar methods Lagrange multiplier spaces inf φ S p S 2 sup û S p S 0 Γ φ û dŝ û L 2 φ L 2 α inf sup φ M h u W h From parametric to physical space Γ φ u ds u L 2 φ L 2 α W h = {u Γ, u V S,h } H 1 0 (Γ) Variable change: Γ φ u ds = Γ φ û ρ dŝ, Super-convergence results a Π : L 2 ( Γ) M h φρ Π( φρ) L 2 (ˆΓ) Ch φ L 2 (ˆΓ) a Superconvergence in Galerkin Finite Element Methods. L. Wahlbin. Lecture Notes in Mathematics (1995), Vol. 1605, Springer, Berlin Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
32 Isogeometric mortar methods Lagrange multiplier spaces inf φ S p S 2 sup û S p S 0 Γ φ û dŝ û L 2 φ L 2 α inf sup φ M h u W h From parametric to physical space Γ φ u ds u L 2 φ L 2 α W h = {u Γ, u V S,h } H0 1(Γ) Variable change: Γ φ u ds = Γ φ û ρ dŝ, Super-convergence results a For h small enough, the stability holds! a Superconvergence in Galerkin Finite Element Methods. L. Wahlbin. Lecture Notes in Mathematics (1995), Vol. 1605, Springer, Berlin Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
33 Isogeometric mortar methods Lagrange multiplier spaces Dependence of the stability constant on the degree Stability constant Pp Pp Pp P(p 2) p Exponential dependence in p Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
34 Isogeometric mortar methods Lagrange multiplier spaces - Summary For V S,h a push-forward of a NURBS space of degree p S, and : M h a push-forward of a spline space of degree p S = problems at cross points. M h a push-forward of a spline space of degree p S 1 = stability is NOK. M h a push-forward of a spline space of degree p S 2 = stability is OK. And more generally, for a push-forward of a spline space of degree p S 2k (k N and k > 1). Pairing p S p S Pairing p S (p S 2) u u h L 2 (Ω) p S + 1 p S λ λ h L 2 (Γ) p S 1 2 p S 1 Table: Optimal order of convergence of isogeometric mortar methods Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
35 Where we are 1 IGA basics 2 Isogeometric mortar methods The discrete mortar problem Lagrange multiplier spaces A key issue: the evaluation of the mortar integrals 3 Applications in contact mechanics Contact problem definition Example of treatment of deformable-deformable contact in a Lagrange multiplier formulation Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
36 Isogeometric mortar methods A key issue: the evaluation of the mortar integrals We recall the bilinear form b(u, φ) = φ u S ΓS Γ S ds φ u M ΓM Γ S ds. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
37 Isogeometric mortar methods A key issue: the evaluation of the mortar integrals We recall the bilinear form b(u, φ) = φ u Γ S ds φ u + Γ S ds. A challenge: the evaluation of the second integral due to the product of functions which are defined on different meshes. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
38 Isogeometric mortar methods A key issue: the evaluation of the mortar integrals We recall the bilinear form b(u, φ) = φ u Γ S ds φ u + Γ S ds. A challenge: the evaluation of the second integral due to the product of functions which are defined on different meshes. Finite element litterature Dedicated quadrature rules have been studied. a,b Approximating these integrals taking the quadrature of one side is inducing large (consistency or approximation) errors. a Numerical quadrature and mortar methods. L. Cazabeau, C. Lacour, Y. Maday. Computational Science for the 21st Century, John Wiley and Sons (1997), pp b The influence of quadrature formulas in 2D and 3D mortar element methods. Y. Maday, F. Rapetti, B.I. Wohlmuth. Recent developments in domain decomposition methods. Some papers of the workshop on domain decomposition. ETH Zürich, Switzerland, June , pp Springer (2002) Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
39 Isogeometric mortar methods A key issue: the evaluation of the mortar integrals Slave integration approach One could expect the sensitivity with respect to the slave quadrature to be less than for finite element methods. a(ũ h, v h ) + (v h v + h ) λ h = f (v h ), v h V h, (ũ h ũ+ h ) φ h = 0, φ h M h. Mixed integration approach We also consider another approach which leads to a non-symmetric saddle point problem. It is motivated by different requirements for the integration of the primal and multiplier test functions. a(ũ h, v h ) + v h λ h + v + h λ h = f (v h ), v h V h, (ũ h ũ+ h ) φ h = 0, φ h M h. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
40 λ Isogeometric mortar methods A key issue: the evaluation of the mortar integrals Numerical results Problem settings Standard Poisson problem u = f Analytical solution slave y 1 0 x x Boundary conditions imposed such that no need to handle a cross point master meshes at refinement level 1 Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
41 Isogeometric mortar methods A key issue: the evaluation of the mortar integrals Numerical results - Slave integration approach L 2 primal and multiplier errors as a function of the number of additional quadrature points 10 4 u uh L 2 (Ω) Exact Refinement level number u n λ h L 2 (Γ) Exact Refinement level number P5 P5 Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
42 Isogeometric mortar methods A key issue: the evaluation of the mortar integrals Numerical results - Mixed integration approach L 2 primal error as a function of the number of additional quadrature points u uh L 2 (Ω) Exact O(h 4.5 ) Refinement level number P4 P2 Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
43 Isogeometric mortar methods A key issue: the evaluation of the mortar integrals Numerical results - Mixed integration approach L 2 primal error as a function of the number of additional quadrature points u uh L 2 (Ω) Exact O(h 2.5 ) Refinement level number P4 P0 Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
44 Where we are 1 IGA basics 2 Isogeometric mortar methods 3 Applications in contact mechanics Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
45 Applications in contact mechanics Contact problem definition - Strong formulation Assumption: frictionless contact the contact variables are thus the gap g N and the contact pressure σ N defined as: 2 types g N = (x S x M ) n M = [x], σ N = (σ S n S ) n M. Rigid-Deformable contact Deformable-Deformable contact g N = u S n M + (X S X M ) n M g N = (u S u M ) n M +(X S X M ) n M Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
46 Applications in contact mechanics Contact problem definition - Strong formulation Assumption: frictionless contact the contact variables are thus the gap g N and the contact pressure σ N defined as: 2 types g N = (x S x M ) n M = [x], σ N = (σ S n S ) n M. Rigid-Deformable contact Deformable-Deformable contact g N = u S n M + g g N = (u S u M ) n M +g Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
47 Applications in contact mechanics Contact problem definition - Strong formulation Assumption: frictionless contact the contact variables are thus the gap g N and the contact pressure σ N defined as: 2 types g N = (x S x M ) n M = [x], σ N = (σ S n S ) n M. Rigid-Deformable contact Deformable-Deformable contact On Ω = k Ω k with either k = {S}, or k = {S, M}, the mechanical equilibrium with contact conditions is written div(σ k ) = f k in Ω k, u k = u D,k on Γ D,k, σ k n k = l k on Γ N,k, g N 0, σ N 0, g N σ N = 0 on Γ C,S. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
48 Applications in contact mechanics Contact problem definition - Weak formulations Different formulations to enforce the contact constraints: penalty, Lagrange multiplier, augmented lagrangian, Nitsche s. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
49 Applications in contact mechanics Contact problem definition - Weak formulations Different formulations to enforce the contact constraints: penalty, Lagrange multiplier, augmented lagrangian, Nitsche s. In the following, we focus on a Lagrange multiplier formulation. To do so, we introduce a multiplier space M defined by M = {φ H 1/2 (Γ C,S )}, as its subspace M = {φ M, Γ C,S φ [u] ds 0, [u] H 1/2 (Γ C,S ), [u] 0}. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
50 Applications in contact mechanics Contact problem definition - Weak formulations Different formulations to enforce the contact constraints: penalty, Lagrange multiplier, augmented lagrangian, Nitsche s. In the following, we focus on a Lagrange multiplier formulation. We also introduce the following bilinear and linear forms: b : V M R, b(u, φ) = φ [u] ds, Γ C,S g : M R, g(φ) = φ g ds. Γ C,S Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
51 Applications in contact mechanics Contact problem definition - Weak formulations Different formulations to enforce the contact constraints: penalty, Lagrange multiplier, augmented lagrangian, Nitsche s. In the following, we focus on a Lagrange multiplier formulation. Contact problem - Lagrange multiplier formulation Find (u, λ) V D M such that { a(u, v) + b(v, λ) = f (v), v V0, b(u, λ φ) g(λ φ), φ M. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
52 Applications in contact mechanics Contact problem definition - Semi-discrete formulations On each subdomain Ω k with either k = {S} (R-D contact) or k = {S, M} (D-D contact), based on its relative NURBS geometrical parametrization, we introduce the displacement approximation space V k,h = {v k = v k F 1 k, v k (N p k (Ξ k )) du }. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
53 Applications in contact mechanics Contact problem definition - Semi-discrete formulations On each subdomain Ω k with either k = {S} (R-D contact) or k = {S, M} (D-D contact), based on its relative NURBS geometrical parametrization, we introduce the displacement approximation space V k,h = {v k = v k F 1 k, v k (N p k (Ξ k )) du }. We note that in the D-D contact case, the discrete product space V h = V k,h V forms a (H 1 (Ω S Ω M )) du non-conforming k={s, M} space discontinuous over the potential region of contact Γ C,S. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
54 Applications in contact mechanics Contact problem definition - Semi-discrete formulations On each subdomain Ω k with either k = {S} (R-D contact) or k = {S, M} (D-D contact), based on its relative NURBS geometrical parametrization, we introduce the displacement approximation space V k,h = {v k = v k F 1 k, v k (N p k (Ξ k )) du }. We note that in the D-D contact case, the discrete product space V h = V k,h V forms a (H 1 (Ω S Ω M )) du non-conforming k={s, M} space discontinuous over the potential region of contact Γ C,S. Need to complete the definition of the discrete formulation to define remaining discrete datas. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
55 Applications in contact mechanics Contact problem definition - To discrete formulations Discrete multiplier space For V S,h a push-forward of a NURBS space of degree p S and : M h a push-forward of a spline space of degree p S = cross points M h a push-forward of a spline space of degree p S 2 = stability OK And more generally, for a push-forward of a spline space of degree p S 2k (k N and k > 1). Discrete gap Use of a lumped L 2 -projection Π into M h defined as Γ Π : R M h, ( Πg N ) j = C,S g N B φ j ds, j = 1,..., n Mh. ds Γ C,S B φ j Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
56 Where we are 1 IGA basics 2 Isogeometric mortar methods The discrete mortar problem Lagrange multiplier spaces A key issue: the evaluation of the mortar integrals 3 Applications in contact mechanics Contact problem definition Example of treatment of deformable-deformable contact in a Lagrange multiplier formulation Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
57 Applications in contact mechanics Example of treatment of deformable-deformable contact in a Lagrange multiplier formulation Numerical resolution strategy The Lagrange multiplier formulation of the contact problem is solved with an active set strategy written as: Find (u h, λ h ) V D,h M h such that a(u h, v h ) + λ h [v h ] ds = f (v h ), v h in V 0,h, ACT ACT φ h [u h ] ds = ACT φ h g ds, φ h in M h. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
58 Applications in contact mechanics Example of treatment of deformable-deformable contact in a Lagrange multiplier formulation Numerical resolution strategy The Lagrange multiplier formulation of the contact problem is solved with an active set strategy written as: Find (u h, λ h ) V D,h M h such that a(u h, v h ) + λ h [v h ] ds = f (v h ), v h in V 0,h, ACT Need to focus on ACT φ h [u h ] ds = ACT the discrete active set region definition, φ h g ds, φ h in M h. suitable integration methods to approximate the mixed term contributions. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
59 Applications in contact mechanics Example of treatment of deformable-deformable contact in a Lagrange multiplier formulation Giving a discrete definition of the dual cone M h, e.g., M h,e = {φ h M h, φ h 0}, M h,w = {φ h M h, φ h B j ds 0, j = 1,..., n Mh }, Γ C,S M h,p = {φ h = α k B k, α k 0, k = 1,..., n Mh }, k consists in giving a discrete definition of the active vs inactive contact region. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
60 Applications in contact mechanics Example of treatment of deformable-deformable contact in a Lagrange multiplier formulation Giving a discrete definition of the dual cone M h, e.g., M h,e = {φ h M h, φ h 0}, M h,w = {φ h M h, φ h B j ds 0, j = 1,..., n Mh }, Γ C,S M h,p = {φ h = α k B k, α k 0, k = 1,..., n Mh }, k consists in giving a discrete definition of the active vs inactive contact region. According to the R-D contact case study, we choose to define it as the support of active multiplier functions. This set of functions is denoted CPλ ACT while the set of inactive ones CPλ INA. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
61 Applications in contact mechanics Example of treatment of deformable-deformable contact in a Lagrange multiplier formulation Active set strategy (i) Initialise CPλ ACT (ii) Compute [ and CP INA λ K B (,CP ACT λ ) B 0 (CPλ ACT, ) (CP ACT λ, CPλ ACT ) ] [ U Λ (CP ACT λ ) ] [ = F G (CP ACT λ ) ] (iii) (iv) Check convergence, i.e.,: CP ACT λ Update CP ACT λ and CP INA λ and CPλ INA stable and go to (ii) until convergence is reached The constraints are checked on the multiplier control values λ j discrete gap control values ( Πg N ) j. The 9 possible cases are: and the λ j < 0, λ j > 0, λ j = 0, ( Πg N ) j = 0 Actif (a) ( Πg N ) j = 0 Inactif (b) ( Πg N ) j = 0 Inactif (c) ( Πg N ) j > 0 Actif (d) ( Πg N ) j > 0 Inactif (e) ( Πg N ) j > 0 Inactif (f) ( Πg N ) j < 0 Actif (g) ( Πg N ) j < 0 Actif (h) ( Πg N ) j < 0 Actif (i) Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
62 Applications in contact mechanics Example of treatment of deformable-deformable contact in a Lagrange multiplier formulation Active set strategy (i) Initialise CPλ ACT (ii) Compute [ and CP INA λ K B (,CP ACT λ ) B 0 (CPλ ACT, ) (CP ACT λ, CPλ ACT ) ] [ U Λ (CP ACT λ ) ] [ = F G (CP ACT λ ) ] (iii) (iv) Check convergence, i.e.,: CP ACT λ Update CP ACT λ [ ] B S B = B M [ ] B B S = B M and CP INA λ and CPλ INA stable and go to (ii) until convergence is reached B S ij = B λ j (x S ) n M (x M ) B S i (x S ) ds, ACT B S ij = B λ i (x S ) n M (x M ) B S j (x S ) ds, ACT B M ij = B λ j (x S ) n M (x M ) B M i (x M ) ds, ACT B M ij = B λ i (x S ) n M (x M ) B M j (x M ) ds. ACT Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
63 Applications in contact mechanics Example of treatment of deformable-deformable contact in a Lagrange multiplier formulation Active set strategy (i) Initialise CPλ ACT (ii) Compute [ and CP INA λ K B (,CP ACT λ ) B 0 (CPλ ACT, ) (CP ACT λ, CPλ ACT ) ] [ U Λ (CP ACT λ ) ] [ = F G (CP ACT λ ) ] (iii) (iv) Check convergence, i.e.,: CP ACT λ Update CP ACT λ and CP INA λ and CPλ INA stable and go to (ii) until convergence is reached B S ij = ACT λ B S j (x S ) n M (x M ) B S i (x S ), B S ij = ACT λ B S i (x S ) n M (x M ) B S j (x S ), B M ij = ACT λ B M j ( x S ) n M (x M ) B M i (x M ), B M ij = ACT λ B S i (x S ) n M (x M ) B M j (x M ). Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
64 Applications in contact mechanics Example of treatment of deformable-deformable contact in a Lagrange multiplier formulation Numerical results - Test 1 - Transmission of a constant pressure Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
65 Applications in contact mechanics Example of treatment of deformable-deformable contact in a Lagrange multiplier formulation Numerical results - Test 1 - Transmission of a constant pressure P2 P2, non-conforming meshes, 4 2 elements for each subdomain Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
66 Applications in contact mechanics Example of treatment of deformable-deformable contact in a Lagrange multiplier formulation Numerical results - Test 1 - Transmission of a constant pressure P2 P2, non-conforming meshes, elements for each subdomain Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
67 Applications in contact mechanics Example of treatment of deformable-deformable contact in a Lagrange multiplier formulation Numerical results - Test 2 conforming meshes non-conforming meshes Full contact vs Partial contact at refinement level 2 Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
68 Applications in contact mechanics Example of treatment of deformable-deformable contact in a Lagrange multiplier formulation Numerical results - Test 2 - Full contact L 2 multiplier errors 10 2 λ λh L 2 (ΓC,S) Conf. meshes Mix. int. app. Sla. int. app. O(h 0.5 ), O(h 1 ),O(h 1.5 ) primal dof number P2 P2 Ref: P4 P4, elements for each subdomain Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
69 Applications in contact mechanics Example of treatment of deformable-deformable contact in a Lagrange multiplier formulation Numerical results - Test 2 - Full contact P2 P2 Ref. level nb. 3 Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
70 Applications in contact mechanics Example of treatment of deformable-deformable contact in a Lagrange multiplier formulation Numerical results - Test 2 - Partial contact L 2 multiplier errors 10 2 λ λh L 2 (ΓC,S) Conf. meshes Mix. Int. App. Slav. int. App. O(h 0.5 ), O(h 1 ),O(h 1.5 ) primal dof number P2 P2 Ref: P4 P4, elements for each subdomain Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
71 Where we are 1 IGA basics 2 Isogeometric mortar methods The discrete mortar problem Lagrange multiplier spaces A key issue: the evaluation of the mortar integrals 3 Applications in contact mechanics Contact problem definition Example of treatment of deformable-deformable contact in a Lagrange multiplier formulation Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
72 Conclusion & perspectives Theoretical and numerical study of isogeometric mortar methods. Study of alternative integration methods to alleviate the construction of the merged mesh. From these former results, proposition of contact methods for rigid-deformable and deformable-deformable contact. Additionally, work on a segmentation process for IGA mortar applications. It has been thought to be suitable for mortar-like contact methods. Perspectives: generalisation of the proposed segmentation process, extension of the methods to the large deformation cases, consideration of contact problems with friction. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
73 Conclusion & perspectives Isogeometric mortar methods. E. Brivadis, A. Buffa, B. Wohlmuth and L. Wunderlich. Comput. Methods Appl. Mech. Eng. 284 (2015), pp The Influence of Quadrature Errors on Isogeometric Mortar Methods. E. Brivadis, A. Buffa, B. Wohlmuth and L. Wunderlich. Isogeometric Analysis and Applications Ed. by B. Jttler and B. Simeon. Vol , pp Isogeometric contact mechanics applications. E. Brivadis, A. Buffa. In preparation. Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
74 Thanks for your attention Ericka Brivadis Isogeometric mortar & contact January 19th, / 45
Isogeometric mortaring
Isogeometric mortaring E. Brivadis, A. Buffa, B. Wohlmuth, L. Wunderlich IMATI E. Magenes - Pavia Technical University of Munich A. Buffa (IMATI-CNR Italy) IGA mortaring 1 / 29 1 Introduction Splines Approximation
More informationIsogeometric Analysis:
Isogeometric Analysis: some approximation estimates for NURBS L. Beirao da Veiga, A. Buffa, Judith Rivas, G. Sangalli Euskadi-Kyushu 2011 Workshop on Applied Mathematics BCAM, March t0th, 2011 Outline
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 informationDual-Primal Isogeometric Tearing and Interconnecting Solvers for Continuous and Discontinuous Galerkin IgA Equations
Dual-Primal Isogeometric Tearing and Interconnecting Solvers for Continuous and Discontinuous Galerkin IgA Equations Christoph Hofer and Ulrich Langer Doctoral Program Computational Mathematics Numerical
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 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 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 informationAn introduction to Isogeometric Analysis
An introduction to Isogeometric Analysis A. Buffa 1, G. Sangalli 1,2, R. Vázquez 1 1 Istituto di Matematica Applicata e Tecnologie Informatiche - Pavia Consiglio Nazionale delle Ricerche 2 Dipartimento
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 Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements
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 A Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements Matthias Gsell and Olaf Steinbach Institute of Computational Mathematics
More informationAND BARBARA I. WOHLMUTH
A QUASI-DUAL LAGRANGE MULTIPLIER SPACE FOR SERENDIPITY MORTAR FINITE ELEMENTS IN 3D BISHNU P. LAMICHHANE AND BARBARA I. WOHLMUTH Abstract. Domain decomposition techniques provide a flexible tool for the
More informationBasics and some applications of the mortar element method
GAMM-Mitt. 28, No. 2, 97 123 (2005) Basics and some applications of the mortar element method Christine Bernardi 1, Yvon Maday 1, and Francesca Rapetti 2 1 Laboratoire Jacques-Louis Lions, C.N.R.S. & université
More informationThe Mortar Boundary Element Method
The Mortar Boundary Element Method A Thesis submitted for the degree of Doctor of Philosophy by Martin Healey School of Information Systems, Computing and Mathematics Brunel University March 2010 Abstract
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 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 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 informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
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 informationMatrix assembly by low rank tensor approximation
Matrix assembly by low rank tensor approximation Felix Scholz 13.02.2017 References Angelos Mantzaflaris, Bert Juettler, Boris Khoromskij, and Ulrich Langer. Matrix generation in isogeometric analysis
More informationIsogeometric modeling of Lorentz detuning in linear particle accelerator cavities
Isogeometric modeling of Lorentz detuning in linear particle accelerator cavities Mauro Bonafini 1,2, Marcella Bonazzoli 3, Elena Gaburro 1,2, Chiara Venturini 1,3, instructor: Carlo de Falco 4 1 Department
More informationSpline Element Method for Partial Differential Equations
for Partial Differential Equations Department of Mathematical Sciences Northern Illinois University 2009 Multivariate Splines Summer School, Summer 2009 Outline 1 Why multivariate splines for PDEs? Motivation
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 informationTime domain boundary elements for dynamic contact problems
Time domain boundary elements for dynamic contact problems Heiko Gimperlein (joint with F. Meyer 3, C. Özdemir 4, D. Stark, E. P. Stephan 4 ) : Heriot Watt University, Edinburgh, UK 2: Universität Paderborn,
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 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 space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh Rome, 10th Apr 2017 Joint work with: Emmanuil
More informationThe variational collocation method
The variational collocation method Hector Gomez a,, Laura De Lorenzis b a Departamento de Métodos Matemáticos, Universidade da Coruña, Campus de A Coruña, 15071, A Coruña, Spain. b Institut für Angewandte
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 informationIsogeometric Analysis with Geometrically Continuous Functions on Two-Patch Geometries
Isogeometric Analysis with Geometrically Continuous Functions on Two-Patch Geometries Mario Kapl a Vito Vitrih b Bert Jüttler a Katharina Birner a a Institute of Applied Geometry Johannes Kepler University
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 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 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 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 informationConstruction of wavelets. Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam
Construction of wavelets Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam Contents Stability of biorthogonal wavelets. Examples on IR, (0, 1), and (0, 1) n. General domains
More informationUniversität Stuttgart
Universität Stuttgart Second Order Lagrange Multiplier Spaces for Mortar Finite Elements in 3D Bishnu P. Lamichhane, Barbara I. Wohlmuth Berichte aus dem Institut für Angewandte Analysis und Numerische
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
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 informationA FETI-DP Method for Mortar Finite Element Discretization of a Fourth Order Problem
A FETI-DP Method for Mortar Finite Element Discretization of a Fourth Order Problem Leszek Marcinkowski 1 and Nina Dokeva 2 1 Department of Mathematics, Warsaw University, Banacha 2, 02 097 Warszawa, Poland,
More informationBenchmarking high order finite element approximations for one-dimensional boundary layer problems
Benchmarking high order finite element approximations for one-dimensional boundary layer problems Marcello Malagù,, Elena Benvenuti, Angelo Simone Department of Engineering, University of Ferrara, Italy
More informationTrefftz-discontinuous Galerkin methods for time-harmonic wave problems
Trefftz-discontinuous Galerkin methods for time-harmonic wave problems Ilaria Perugia Dipartimento di Matematica - Università di Pavia (Italy) http://www-dimat.unipv.it/perugia Joint work with Ralf Hiptmair,
More informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationCondition number estimates for matrices arising in the isogeometric discretizations
www.oeaw.ac.at Condition number estimates for matrices arising in the isogeometric discretizations K. Gahalaut, S. Tomar RICAM-Report -3 www.ricam.oeaw.ac.at Condition number estimates for matrices arising
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 informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationAn Extended Finite Element Method for a Two-Phase Stokes problem
XFEM project An Extended Finite Element Method for a Two-Phase Stokes problem P. Lederer, C. Pfeiler, C. Wintersteiger Advisor: Dr. C. Lehrenfeld August 5, 2015 Contents 1 Problem description 2 1.1 Physics.........................................
More informationTrefftz-Discontinuous Galerkin Methods for Acoustic Scattering - Recent Advances
Trefftz-Discontinuous Galerkin Methods for Acoustic Scattering - Recent Advances Ilaria Perugia Dipartimento di Matematica Università di Pavia (Italy) In collaboration with Ralf Hiptmair, Christoph Schwab,
More informationStochastic multiscale modeling of subsurface and surface flows. Part III: Multiscale mortar finite elements for coupled Stokes-Darcy flows
Stochastic multiscale modeling of subsurface and surface flows. Part III: Multiscale mortar finite elements for coupled Stokes-Darcy flows Ivan otov Department of Mathematics, University of Pittsburgh
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh & Department of Mathematics, University of
More informationInstitut de Recherche MAthématique de Rennes
LMS Durham Symposium: Computational methods for wave propagation in direct scattering. - July, Durham, UK The hp version of the Weighted Regularization Method for Maxwell Equations Martin COSTABEL & Monique
More informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
More informationA Priori Error Analysis for Space-Time Finite Element Discretization of Parabolic Optimal Control Problems
A Priori Error Analysis for Space-Time Finite Element Discretization of Parabolic Optimal Control Problems D. Meidner and B. Vexler Abstract In this article we discuss a priori error estimates for Galerkin
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 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 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 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 information20. A Dual-Primal FETI Method for solving Stokes/Navier-Stokes Equations
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 23 DDM.org 2. A Dual-Primal FEI Method for solving Stokes/Navier-Stokes
More informationGeoPDEs. An Octave/Matlab software for research on IGA. R. Vázquez. IMATI Enrico Magenes, Pavia Consiglio Nazionale della Ricerca
GeoPDEs An Octave/Matlab software for research on IGA R. Vázquez IMATI Enrico Magenes, Pavia Consiglio Nazionale della Ricerca Joint work with C. de Falco and A. Reali Supported by the ERC Starting Grant:
More informationHybridized DG methods
Hybridized DG methods University of Florida (Banff International Research Station, November 2007.) Collaborators: Bernardo Cockburn University of Minnesota Raytcho Lazarov Texas A&M University Thanks:
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 informationA BIVARIATE SPLINE METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS IN NON-DIVERGENCE FORM
A BIVARIAE SPLINE MEHOD FOR SECOND ORDER ELLIPIC EQUAIONS IN NON-DIVERGENCE FORM MING-JUN LAI AND CHUNMEI WANG Abstract. A bivariate spline method is developed to numerically solve second order elliptic
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 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 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 informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationarxiv: v1 [math.na] 1 May 2013
arxiv:3050089v [mathna] May 03 Approximation Properties of a Gradient Recovery Operator Using a Biorthogonal System Bishnu P Lamichhane and Adam McNeilly May, 03 Abstract A gradient recovery operator based
More informationSpace-time sparse discretization of linear parabolic equations
Space-time sparse discretization of linear parabolic equations Roman Andreev August 2, 200 Seminar for Applied Mathematics, ETH Zürich, Switzerland Support by SNF Grant No. PDFMP2-27034/ Part of PhD thesis
More informationBUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM
BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM ERIK BURMAN AND BENJAMIN STAMM Abstract. We propose a low order discontinuous Galerkin method for incompressible flows. Stability of the
More informationIso-P2 P1/P1/P1 Domain-Decomposition/Finite-Element Method for the Navier-Stokes Equations
Contemporary Mathematics Volume 28, 998 B -828-988--35-6 Iso-P2 P/P/P Domain-Decomposition/Finite-Element Method for the Navier-Stokes Equations Shoichi Fujima. Introduction In recent years, parallel computers
More informationLocal Mesh Refinement with the PCD Method
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 1, pp. 125 136 (2013) http://campus.mst.edu/adsa Local Mesh Refinement with the PCD Method Ahmed Tahiri Université Med Premier
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 informationUnified A Posteriori Error Control for all Nonstandard Finite Elements 1
Unified A Posteriori Error Control for all Nonstandard Finite Elements 1 Martin Eigel C. Carstensen, C. Löbhard, R.H.W. Hoppe Humboldt-Universität zu Berlin 19.05.2010 1 we know of Guidelines for Applicants
More informationSTABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS
STABILIZED DISCONTINUOUS FINITE ELEMENT APPROXIMATIONS FOR STOKES EQUATIONS RAYTCHO LAZAROV AND XIU YE Abstract. In this paper, we derive two stabilized discontinuous finite element formulations, symmetric
More informationNumerical Solutions of Laplacian Problems over L-Shaped Domains and Calculations of the Generalized Stress Intensity Factors
WCCM V Fifth World Congress on Computational Mechanics July 7-2, 2002, Vienna, Austria Eds.: H.A. Mang, F.G. Rammerstorfer, J. Eberhardsteiner Numerical Solutions of Laplacian Problems over L-Shaped Domains
More informationVirtual Element Methods for general second order elliptic problems
Virtual Element Methods for general second order elliptic problems Alessandro Russo Department of Mathematics and Applications University of Milano-Bicocca, Milan, Italy and IMATI-CNR, Pavia, Italy Workshop
More informationDiscontinuous Galerkin Method for interface problem of coupling different order elliptic equations
Discontinuous Galerkin Method for interface problem of coupling different order elliptic equations Igor Mozolevski, Endre Süli Federal University of Santa Catarina, Brazil Oxford University Computing Laboratory,
More informationA Stabilized Finite Element Method for the Darcy Problem on Surfaces
arxiv:1511.03747v1 [math.na] 12 Nov 2015 A Stabilized Finite Element Method for the Darcy Problem on Surfaces Peter Hansbo Mats G. Larson October 8, 2018 Abstract We consider a stabilized finite element
More informationOn a Discontinuous Galerkin Method for Surface PDEs
On a Discontinuous Galerkin Method for Surface PDEs Pravin Madhavan (joint work with Andreas Dedner and Bjo rn Stinner) Mathematics and Statistics Centre for Doctoral Training University of Warwick Applied
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationICES REPORT March R. Kruse, N. Nguyen-Thanh, L. De Lorenzis, and T.J.R. Hughes
ICES REPORT 15-08 March 2015 Isogeometric collocation for large deformation elasticity and frictional contact problems by R. Kruse,. guyen-thanh, L. De Lorenzis, and T.J.R. Hughes The Institute for Computational
More informationApproximation of Geometric Data
Supervised by: Philipp Grohs, ETH Zürich August 19, 2013 Outline 1 Motivation Outline 1 Motivation 2 Outline 1 Motivation 2 3 Goal: Solving PDE s an optimization problems where we seek a function with
More informationReduction of Finite Element Models of Complex Mechanical Components
Reduction of Finite Element Models of Complex Mechanical Components Håkan Jakobsson Research Assistant hakan.jakobsson@math.umu.se Mats G. Larson Professor Applied Mathematics mats.larson@math.umu.se Department
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 informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
More informationAn hp-adaptive Mixed Discontinuous Galerkin FEM for Nearly Incompressible Linear Elasticity
An hp-adaptive Mixed Discontinuous Galerkin FEM for Nearly Incompressible Linear Elasticity Paul Houston 1 Department of Mathematics, University of Leicester, Leicester LE1 7RH, UK (email: Paul.Houston@mcs.le.ac.uk)
More informationNumerical Analysis of Electromagnetic Fields
Pei-bai Zhou Numerical Analysis of Electromagnetic Fields With 157 Figures Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Hong Kong Barcelona Budapest Contents Part 1 Universal Concepts
More informationSupraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives
Supraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives Etienne Emmrich Technische Universität Berlin, Institut für Mathematik, Straße
More informationFinite Elements. Colin Cotter. January 15, Colin Cotter FEM
Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation
More informationPseudo-divergence-free element free Galerkin method for incompressible fluid flow
International Workshop on MeshFree Methods 23 Pseudo-divergence-free element free Galerkin method for incompressible fluid flow Y. Vidal () and A. Huerta (2) Abstract: Incompressible modelling in finite
More informationTHE PATCH RECOVERY FOR FINITE ELEMENT APPROXIMATION OF ELASTICITY PROBLEMS UNDER QUADRILATERAL MESHES. Zhong-Ci Shi and Xuejun Xu.
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 9, Number, January 2008 pp. 63 82 THE PATCH RECOVERY FOR FINITE ELEMENT APPROXIMATION OF ELASTICITY PROBLEMS UNDER
More informationOVERLAPPING SCHWARZ ALGORITHMS FOR ALMOST INCOMPRESSIBLE LINEAR ELASTICITY TR
OVERLAPPING SCHWARZ ALGORITHMS FOR ALMOST INCOMPRESSIBLE LINEAR ELASTICITY MINGCHAO CAI, LUCA F. PAVARINO, AND OLOF B. WIDLUND TR2014-969 Abstract. Low order finite element discretizations of the linear
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 informationNon-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 17, 2017, Lesson 5
Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, attilio.frangi@polimi.it Politecnico di Milano, February 17, 2017, Lesson 5 1 Politecnico di Milano, February 17, 2017, Lesson 5 2 Outline
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 informationGoal. Robust A Posteriori Error Estimates for Stabilized Finite Element Discretizations of Non-Stationary Convection-Diffusion Problems.
Robust A Posteriori Error Estimates for Stabilized Finite Element s of Non-Stationary Convection-Diffusion Problems L. Tobiska and R. Verfürth Universität Magdeburg Ruhr-Universität Bochum www.ruhr-uni-bochum.de/num
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 information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
More informationThe Mortar Wavelet Method Silvia Bertoluzza Valerie Perrier y October 29, 1999 Abstract This paper deals with the construction of wavelet approximatio
The Mortar Wavelet Method Silvia Bertoluzza Valerie Perrier y October 9, 1999 Abstract This paper deals with the construction of wavelet approximation spaces, in the framework of the Mortar method. We
More informationFinite volume method for nonlinear transmission problems
Finite volume method for nonlinear transmission problems Franck Boyer, Florence Hubert To cite this version: Franck Boyer, Florence Hubert. Finite volume method for nonlinear transmission problems. Proceedings
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationngsxfem: Geometrically unfitted discretizations with Netgen/NGSolve (https://github.com/ngsxfem/ngsxfem)
ngsxfem: Geometrically unfitted discretizations with Netgen/NGSolve (https://github.com/ngsxfem/ngsxfem) Christoph Lehrenfeld (with contributions from F. Heimann, J. Preuß, M. Hochsteger,...) lehrenfeld@math.uni-goettingen.de
More information