Time domain boundary elements for dynamic contact problems
|
|
- Norman Thomas
- 5 years ago
- Views:
Transcription
1 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, Germany 3: Universität Stuttgart, Germany 4: Leibniz Universität Hannover, Germany BEM on the Saar 207, Saarbrücken May 29, 207 H. Gimperlein TDBEM for dynamic contact BEM on the Saar 207 / 9
2 Wave equation U : R t Ω c x R 2 t U U = 0 in R Ω c, Ω c = R d \ Ω, U = 0 for t 0. Real problem: Lamé equation for tire dynamics tire does not penetrate road Key engineering problem, analysis challenging H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
3 Wave equation U : R t Ω c x R 2 t U U = 0 in R Ω c, Ω c = R d \ Ω, U = 0 for t 0. Real problem: Lamé equation for tire dynamics tire does not penetrate road Scalar model problem: wave equation with contact boundary conditions on G Γ := Ω { U G 0, U ν G h, U G > 0 = U G = h. ν H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
4 Wave equation U : R t Ω c x R 2 t U U = 0 in R Ω c, Ω c = R d \ Ω, U = 0 for t 0. Real problem: Lamé equation for tire dynamics tire does not penetrate road Scalar model problem: wave equation with contact boundary conditions on G Γ := Ω Recent works (computational & stabilized FEM): Renard, Hauret, Le Tallec, Wohlmuth,... H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
5 Contents U : R t Ω c x R t 2 U U = 0 in R t Ω c x, Ω c = R d \ Ω U = 0 for t 0. Analysis of dynamic contact: variational inequality on Γ := Ω, well-posedness Basics of time domain BEM Numerical analysis of dynamic contact: mixed formulation, first a priori estimates, Uzawa Numerical experiments H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
6 Dirichlet-Neumann operator for wave equation U : R t Ω x R t 2 U U = 0 in R t Ω c x, Ω c = R d \ Ω U = 0 for t 0. Dirichlet-Neumann operator on Γ := Ω S(U Γ ) = U ν contact boundary conditions on G Γ (and U = 0 on Γ \ G) { U 0, U ν h, U > 0 = U ν = h. Γ H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
7 Dirichlet-Neumann operator for wave equation U : R t Ω x R t 2 U U = 0 in R t Ω c x, Ω c = R d \ Ω U = 0 for t 0. Dirichlet-Neumann operator on Γ := Ω S(U Γ ) = U ν contact boundary conditions on G Γ (and U = 0 on Γ \ G) { { U 0, S(U Γ ) h, U 0, S(U Γ ) h, U > 0 = S(U Γ ) = h. U (S(U Γ ) h) = 0. Γ H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
8 Dirichlet-Neumann operator for wave equation U : R t Ω x R t 2 U U = 0 in R t Ω c x, Ω c = R d \ Ω U = 0 for t 0. contact boundary conditions on G Γ (and U = 0 on Γ \ G) { U 0, S(U Γ ) h, U > 0 = S(U Γ ) = h. Variational inequality: Find 0 u = U G such that Su, v u h, v u 0 v. { U 0, S(U Γ ) h, U (S(U Γ ) h) = 0. H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
9 Analytic set-up: Space-time Sobolev spaces space time anisotropic Sobolev spaces H r σ(r +, H s (Γ)), σ > 0: H r σ(r +, H s (R d )) defined using Fourier Laplace transform { ψ : supp ψ R + R d, R+iσ dω } R dξ ω 2r ( ω 2 + ξ 2 ) s Fψ(ω, ξ) 2 < d H r σ(r +, H s (R d )) supported distributions H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
10 Analytic set-up: Space-time Sobolev spaces space time anisotropic Sobolev spaces H r σ(r +, H s (Γ)), σ > 0: H r σ(r +, H s (R d )) defined using Fourier Laplace transform { ψ : supp ψ R + R d, R+iσ dω } R dξ ω 2r ( ω 2 + ξ 2 ) s Fψ(ω, ξ) 2 < d H r σ(r +, H s (R d )) supported distributions Variational inequality: Find 0 u H 2 σ (R +, H 2 (G)) such that Su, v u h, v u 0 v H 2 σ (R +, H 2 (G)). H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
11 Existence of solutions Ω = R d, G Γ bounded polygon. Existence of solutions: Lebeau Schatzmann 84, students of Eskin, including Cooper 99. Theorem (Cooper, 99): For h H 3 2 σ (R +, H 2 (Γ c )) the contact problem admits a unique solution u H 2 σ (R +, H 2 (Γ c )) +. H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
12 Existence of solutions Ω = R d, G Γ bounded polygon. Existence of solutions: Lebeau Schatzmann 84, students of Eskin, including Cooper 99. Theorem (Cooper, 99): For h H 3 2 σ (R +, H 2 (Γ c )) the contact problem admits a unique solution u H 2 σ (R +, H 2 (Γ c )) +. Basic strategy: elliptic regularization: S S ε( 2 t + R d ) solution exists for ε > 0 by standard theory control ε 0 + using improved coercivity estimate for flat Γ (from Fourier analysis): ϕ 2 2, 2,σ, σ Sϕ, ϕ C ϕ 2 0, 2,σ, H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
13 Discretization Γ = M i= Γ i (quasi-uniform) triangulation V p h piecewise polynomial functions of degree p on Γ = M i= Γ i (continuous if p ) [0, T) = L n= [t n, t n ), t n = n( t) V q t piecewise polynomial functions of degree q in time (continuous and vanishing at t = 0 if q ) tensor products in space-time: V p,q h, t = Vp h Vq t Ṽ p,q h, t subspace vanishing at boundary H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
14 A priori error analysis Variational inequality: Su, v u h, v u 0 v H 2 σ (R +, H 2 (Γc )). Discretized variational inequality in V p,q h, t : S h, t u h, t, v h, t u h, t h, v h, t u h, t 0 v h, t Ṽ p,q h, t. Crucial: not S t. H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
15 A priori error analysis Variational inequality: Su, v u h, v u 0 v H 2 σ (R +, H 2 (Γc )). Discretized variational inequality in V p,q h, t : S h, t u h, t, v h, t u h, t h, v h, t u h, t 0 v h, t Ṽ p,q h, t. Theorem Let h H 3 2 σ (R +, H 2 (G)). Then u u t,h 2 2, 2,σ, σ inf ( h S σ u 0 ϕ t,h Ṽ p,q 2,,σ u ϕ t,h 2 2, 2,σ, h, t + u ϕ t,h 2 2, 2,σ, + O(S S h, t)). Assume u H 2 +ϵ σ (R +, H 2 +ϵ (G)) for some ϵ > 0, S = S h, t. Then u u t,h 2 2, 2,σ, (hϵ + ( t) 2 +ϵ ) u 2 +ϵ, 2 +ϵ,σ + (h2ϵ + ( t) 2ϵ ) u 2 2 +ϵ, 2 +ϵ, H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
16 Mixed formulation for improved energy conservation Variational inequality: Su, v u h, v u 0 v H 2 σ (R +, H 2 (Γc )). Lagrange multiplier λ = Su h indicates contact / forces Mixed formulation: Find (u, λ) H 2 σ (R +, H 2 (G)) H 2 σ (R +, H 2 (G)) + such that { (a) Su, v λ, v = h, v (b) u, µ λ 0, for all (v, µ) H 2 σ (R +, H 2 (G)) H 2 σ (R +, H 2 (G)) +. H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
17 Mixed formulation for improved energy conservation Variational inequality: Su, v u h, v u 0 v H 2 σ (R +, H 2 (Γc )). Lagrange multiplier λ = Su h indicates contact / forces Discretized mixed formulation: Find (u h, t, λ h2, t 2 ) Ṽ, h, t V 0,0,+ h 2, t 2 such that { (a) Su h, t, v h, t λ h2, t 2, v h, t = h, v h, t (b) u h, t, µ h2, t 2 λ h2, t 2 0, for all (v h, t, µ h2, t 2 ) Ṽ, h, t V 0,0,+ h 2, t 2. Both the continuous and discretized variational inequalities admit unique solutions. H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
18 Mixed formulation for improved energy conservation Discretized mixed formulation: Find (u h, t, λ h2, t 2 ) Ṽ, h, t V 0,0,+ h 2, t 2 such that { (a) Su h, t, v h, t λ h2, t 2, v h, t = h, v h, t (b) u h, t, µ h2, t 2 λ h2, t 2 0, for all (v h, t, µ h2, t 2 ) Ṽ, h, t V 0,0,+ h 2, t 2. Error analysis based on: Theorem (space time inf sup) Let C > 0 sufficiently small, and max{h, t } min{h 2, t 2 } α > 0 such that for all λ t2,h 2 : < C. Then there exists v t,h sup, λ t2,h 2 α λ t2,h v t,h v t,h 2 0, 0, 2,σ, 2,σ. H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
19 Mixed formulation for improved energy conservation Discretized mixed formulation: Find (u h, t, λ h2, t 2 ) Ṽ, h, t V 0,0,+ h 2, t 2 such that { (a) Su h, t, v h, t λ h2, t 2, v h, t = h, v h, t (b) u h, t, µ h2, t 2 λ h2, t 2 0, for all (v h, t, µ h2, t 2 ) Ṽ, h, t V 0,0,+ h 2, t 2. The inf-sup condition implies a priori estimates similar to Brezzi Hager Raviart 78: Theorem (a priori error estimate) λ λ t2,h 2 0,,σ inf λ λ t2,h 2 2 0, λ 2,σ + ( t ) 2 u u t,h 2, 2,σ, t2,h 2 u u t,h 2, 2,σ, σ inf v t,h u v t,h 2, 2,σ, { } + inf λ t2,h 2 λ λ 2, 2,σ + λ t2,h 2 λ t2,h 2 2, 2,σ t2,h 2 H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
20 Proof of a priori estimate: λ λ t2,h 2 Inf-sup condition: For all λ t2,h 2 Algebra α λ t2 h 2 λ λ t2,h t2,h 2 0,,σ sup 2 λ t2,h 2, v t,h 2 v t,h v t,h 0, 2,σ, S(u t,h = sup u), v t,h + λ λ t2,h 2, v t,h v t,h v t,h 0, First term: duality and inverse inequality in time 2,σ,. S(u t,h u), v t,h S σ (u t,h u) 2, 2,σ v t,h 2, 2,σ, Second term: duality u t,h u 2, 2,σ, ( t ) 2 v t,h 0, 2,σ,. λ λ t2,h 2, v t,h λ λ t2,h 2 0, 2,σ v t,h 0, 2,σ,. Combine into estimate for λ λ t2,h 2 0, 2,σ. H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
21 Proof of a priori estimate: u u t,h Estimate for u without inf-sup: Use coercivity u t,h v t,h 2 2, 2,σ, σ S(u t,h v t,h ), u t,h v t,h Then Galerkin orthogonality, mapping properties, etc. standard! H. Gimperlein TDBEM for dynamic contact BEM on the Saar 207 / 9
22 Computing the Dirichlet-Neumann operator S(u Γ ) = du dν Computation through layer potentials: S = W + (K 2 )V (K 2 ) Γ Kφ(t, x) = K φ(t, x) = R + Γ R + Γ Wφ(t, x) = Here G(t, x) = δ(t x ) 4π x R + Γ G (t τ, x y) φ(τ, y) dτ ds y, n y G (t τ, x y) φ(τ, y) dτ ds y, n x 2 G n x n y (t τ, x y) φ(τ, y) dτ ds y. a fundamental solution of wave equation. H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
23 Computing the Dirichlet-Neumann operator S(u Γ ) = du dν Computation through layer potentials: S = W + (K 2 )V (K 2 ) Γ Weak formulation of Su = h H 3 2 σ (R +, H 2 (Γ)) on flat screen: Find (u, v) H 2 σ (R +, H 2 (Γ)) H 2 σ (R +, H 2 (Γ)) such that Wu + (K )v, ϕ = h, ϕ, 2 (K )u + Vv, Ψ = 0, 2 for all (ϕ, Ψ) H 2 σ (R +, H 2 (G)) H 2 σ (R +, H 2 (G)). H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
24 Computing the Dirichlet-Neumann operator S(u Γ ) = du dν Computation through layer potentials: Γ S = W + (K 2 )V (K 2 ) Discretized weak formulation of Su = h on flat screen: Find (u h, t, v h, t ) V, V, such that Wu h, t (K 2 )v h, t, ϕ h, t = h, ϕ h, t, (K 2 )u h, t + Vv h, t, t Ψ h, t = 0, for all (ϕ h, t, Ψ h, t ) V,0 V,0. Note: No t in first equation S coercive on flat screen, important later. H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
25 Uzawa algorithm Discretized mixed formulation: { (a) Su h, t, v h, t λ h2, t 2, v h, t = h, v h, t (b) u h, t, µ h2, t 2 λ h2, t 2 0, Two algorithms: Uzawa in space-time / Uzawa as time stepping scheme H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
26 Uzawa algorithm Discretized mixed formulation: { (a) Su h, t, v h, t λ h2, t 2, v h, t = h, v h, t (b) u h, t, µ h2, t 2 λ h2, t 2 0, Two algorithms: Uzawa in space-time / Uzawa as time stepping scheme Space-time Uzawa algorithm Choose ρ > 0. Set k = 0, y (0) = 0. WHILE stopping criterion not satisfied solve: Sx (k) = h + y (k) compute: y (k+) = Pr K (y (k) ρx (k) ), where (Pr K y) i = max{y i, 0} k k + ENDWHILE H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
27 Uzawa algorithm Two algorithms: Uzawa in space-time / Uzawa as time stepping scheme Space-time Uzawa algorithm Choose ρ > 0. Set k = 0, y (0) = 0. WHILE stopping criterion not satisfied solve: Sx (k) = h + y (k) compute: y (k+) = Pr K (y (k) ρx (k) ), where (Pr K y) i = max{y i, 0} k k + ENDWHILE Recall: Sx (k) = h + y (k) is a linear system with matrix ( W ((K) 2 I) ) (K 2 I) V H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
28 Uzawa algorithm Two algorithms: Uzawa in space-time / Uzawa as time stepping scheme Space-time Uzawa algorithm Choose ρ > 0. Set k = 0, y (0) = 0. WHILE stopping criterion not satisfied solve: Sx (k) = h + y (k) compute: y (k+) = Pr K (y (k) ρx (k) ), where (Pr K y) i = max{y i, 0} k k + ENDWHILE Lemma The space-time Uzawa algorithm converges, provided that 0 < ρ < 2C σ. Here C σ is the coercivity constant of S. H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
29 Uzawa algorithm Two algorithms: Uzawa in space-time / Uzawa as time stepping scheme Space-time Uzawa algorithm Choose ρ > 0. Set k = 0, y (0) = 0. WHILE stopping criterion not satisfied solve: Sx (k) = h + y (k) compute: y (k+) = Pr K (y (k) ρx (k) ), where (Pr K y) i = max{y i, 0} k k + ENDWHILE We often use a reformulation of the variational inequality as a variational inequality in every time step. efficient Uzawa in every time step. It converges in every time step. H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
30 Experiment : Dirichlet-Neumann equation Su = h on Γ = S 2 for times [0, 5], where h = Neumann data of u(t, r) = ( 3 π(4 t) 4 cos( 2 ) + 4 cos(π(4 t)))[h(4 t) H( t)]. H = Heaviside function. Approximation on uniform icosahedral meshes with t h 0.6. Plot L 2 ([0, T] Γ)-error vs. degrees of freedom 0 0 Dirichlet to Neumann equation relative L 2 error in space time DOF H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
31 Experiment 2: Contact Su h on flat screen Γ = [ 2, 2] 2 {0} G = [, ] 2 {0}, 0 < t < 5, t x = 0.7 h = e 2t t 4 cos(2πx) cos(2πy)χ [ 0.25,0.25] (x)χ [ 0.25,0.25] (y) relative L 2 ([0, T] Γ)-error against benchmark with 2800 triangles, t = Space Time Uzawa relative L 2 error in space time DOF H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
32 Experiment 2: Uzawa for contact Su h Γ = [ 2, 2] 2 {0} G = [, ] 2 {0}, 0 < t < 6, t = 0. h = e 2t t 4 cos(2πx) cos(2πy)χ [ 0.25,0.25] (x)χ [ 0.25,0.25] (y), 3200 triangles Evolution of the relative difference between solutions in L 2 (Γ) between time step Uzawa and space-time Uzawa. Difference increases sharply around the onset of contact. 0 4 Space Time vs. Time Step Uzawa on 3200 triangles 0 6 relative error Time H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
33 Experiment 3: Su h with non-flat contact Γ = [ 2, 2] 3 Γ c = 3 faces (top, left, back), 0 < t < 4, t x = 0.7 h as before, on every contact face relative L 2 ([0, T] Γ)-error against benchmark 0 Space Time Uzawa relative error in space time DOF H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
34 Experiment 3: Su h with non-flat contact H. Gimperlein (a) t=3 (b) t=5 (c) t=5.5 (d) t=6 TDBEM for dynamic contact BEM on the Saar / 9
35 Conclusions & References first a priori analysis of contact problems, using TDBEM convergence of Galerkin and mixed methods for flat contact (only case where existence of solution is known) based on inf-sup condition in space-time and refined coercivity estimates for Dirichlet-Neumann convergence of space-time Uzawa and time step Uzawa Analysis and numerics widely open: adaptivity, friction, stabilized methods,... Final remark: Theory and numerics extend to punch problems with the Neumann-Dirichlet operator. HG, Meyer, Oezdemir, Stephan, Time domain boundary elements for dynamic contact problems, preprint. H. Gimperlein TDBEM for dynamic contact BEM on the Saar / 9
Time domain boundary elements for dynamic contact problems
Time domain boundary elements for dynamic contact problems Heiko Gimperlein Fabian Meyer Ceyhun Özdemir Ernst P. Stephan (dedicated to Erwin Stein on the occasion of his 85th birthday) Abstract This article
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 informationA DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT
A DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT HUOYUAN DUAN, SHA LI, ROGER C. E. TAN, AND WEIYING ZHENG Abstract. To deal with the divergence-free
More informationDiscontinuous Petrov-Galerkin Methods
Discontinuous Petrov-Galerkin Methods Friederike Hellwig 1st CENTRAL School on Analysis and Numerics for Partial Differential Equations, November 12, 2015 Motivation discontinuous Petrov-Galerkin (dpg)
More 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 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 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 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 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 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 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 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 informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationA MHD problem on unbounded domains - Coupling of FEM and BEM
A MHD problem on unbounded domains - Coupling of FEM and BEM Wiebke Lemster and Gert Lube Abstract We consider the MHD problem on R 3 = Ω Ω E, where Ω is a bounded, conducting Lipschitz domain and Ω E
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 informationInterface and contact problems
Interface and contact problems Heiko Gimperlein Heriot Watt University and Maxwell Institute, Edinburgh and Universität Paderborn TU Wien May 30, 2016 H. Gimperlein (Edinburgh) Interface and contact problems
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 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 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 informationEnergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations of the Navier-Stokes equations
INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 19007; 1:1 [Version: 00/09/18 v1.01] nergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations
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 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 informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
More informationarxiv: v1 [math.na] 5 Jun 2018
PRIMAL-DUAL WEAK GALERKIN FINIE ELEMEN MEHODS FOR ELLIPIC CAUCHY PROBLEMS CHUNMEI WANG AND JUNPING WANG arxiv:1806.01583v1 [math.na] 5 Jun 2018 Abstract. he authors propose and analyze a well-posed numerical
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 informationAn interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes
An interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes Vincent Heuveline Friedhelm Schieweck Abstract We propose a Scott-Zhang type interpolation
More informationStabilization Techniques and a Posteriori Error Estimates for the Obstacle Problem
Applied Mathematical Sciences, Vol. 7, 2013, no. 127, 6329-6346 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2013.39504 Stabilization Techniques and a Posteriori Error Estimates for the
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 informationEquilibrated Residual Error Estimates are p-robust
Equilibrated Residual Error Estimates are p-robust Dietrich Braess, Veronika Pillwein and Joachim Schöberl April 16, 2008 Abstract Equilibrated residual error estimators applied to high order finite elements
More informationIntroduction to finite element exterior calculus
Introduction to finite element exterior calculus Ragnar Winther CMA, University of Oslo Norway Why finite element exterior calculus? Recall the de Rham complex on the form: R H 1 (Ω) grad H(curl, Ω) curl
More 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 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 informationPOISSON EQUATIONS JEONGHUN LEE (1)
POISSON EQUATIONS JEONGHUN LEE Throughout this note the symbol X À Y stands for an inequality X CY with some constant C 0 which is independent of mesh sizes h. However, we will use X CY when we want to
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
More informationA FAMILY OF MULTISCALE HYBRID-MIXED FINITE ELEMENT METHODS FOR THE DARCY EQUATION WITH ROUGH COEFFICIENTS. 1. Introduction
A FAMILY OF MULTISCALE HYBRID-MIXED FINITE ELEMENT METHODS FOR THE DARCY EQUATION WITH ROUGH COEFFICIENTS CHRISTOPHER HARDER, DIEGO PAREDES 2, AND FRÉDÉRIC VALENTIN Abstract. We aim at proposing novel
More informationDomain decomposition methods via boundary integral equations
Domain decomposition methods via boundary integral equations G. C. Hsiao a O. Steinbach b W. L. Wendland b a Department of Mathematical Sciences, University of Delaware, Newark, Delaware 19716, USA. E
More informationANALYSIS OF THE FINITE ELEMENT METHOD FOR TRANSMISSION/MIXED BOUNDARY VALUE PROBLEMS ON GENERAL POLYGONAL DOMAINS
ANALYSIS OF THE FINITE ELEMENT METHOD FOR TRANSMISSION/MIXED BOUNDARY VALUE PROBLEMS ON GENERAL POLYGONAL DOMAINS HENGGUANG LI, ANNA MAZZUCATO, AND VICTOR NISTOR Abstract. We study theoretical and practical
More informationA u + b u + cu = f in Ω, (1.1)
A WEIGHTED H(div) LEAST-SQUARES METHOD FOR SECOND-ORDER ELLIPTIC PROBLEMS Z. CAI AND C. R. WESTPHAL Abstract. This paper presents analysis of a weighted-norm least squares finite element method for elliptic
More informationOn Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1
On Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1 Jie Shen Department of Mathematics, Penn State University University Park, PA 1682 Abstract. We present some
More informationUnique continuation for the Helmholtz equation using a stabilized finite element method
Unique continuation for the Helmholtz equation using a stabilized finite element method Lauri Oksanen University College London Based on a joint work with Erik Burman and Mihai Nechita Motivation: recovering
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 informationSTOKES PROBLEM WITH SEVERAL TYPES OF BOUNDARY CONDITIONS IN AN EXTERIOR DOMAIN
Electronic Journal of Differential Equations, Vol. 2013 2013, No. 196, pp. 1 28. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu STOKES PROBLEM
More informationA Space-Time Boundary Element Method for the Wave 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 A Space-Time Boundary Element Method for the Wave Equation Marco Zank and Olaf Steinbach Institut für Numerische Mathematik Space-Time Methods for PDEs,
More informationANALYSIS OF A MODIFIED SCHRÖDINGER OPERATOR IN 2D: REGULARITY, INDEX, AND FEM
ANALYSIS OF A MODIFIED SCHRÖDINGER OPERATOR IN 2D: REGULARITY, INDEX, AND FEM HENGGUANG LI AND VICTOR NISTOR Version: 2.; Revised: April 2 by V.; Run: January 4, 2008 Abstract. Let r = (x 2 +x2 2 )/2 be
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 informationFrom Completing the Squares and Orthogonal Projection to Finite Element Methods
From Completing the Squares and Orthogonal Projection to Finite Element Methods Mo MU Background In scientific computing, it is important to start with an appropriate model in order to design effective
More informationMixed Finite Elements Method
Mixed Finite Elements Method A. Ratnani 34, E. Sonnendrücker 34 3 Max-Planck Institut für Plasmaphysik, Garching, Germany 4 Technische Universität München, Garching, Germany Contents Introduction 2. Notations.....................................
More informationGraded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons
Graded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons Anna-Margarete Sändig, Miloslav Feistauer University Stuttgart, IANS Journées
More informationChapter 5 A priori error estimates for nonconforming finite element approximations 5.1 Strang s first lemma
Chapter 5 A priori error estimates for nonconforming finite element approximations 51 Strang s first lemma We consider the variational equation (51 a(u, v = l(v, v V H 1 (Ω, and assume that the conditions
More informationAn a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element
Calcolo manuscript No. (will be inserted by the editor) An a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element Dietrich Braess Faculty of Mathematics, Ruhr-University
More informationA STABILIZED NONCONFORMING QUADRILATERAL FINITE ELEMENT METHOD FOR THE GENERALIZED STOKES EQUATIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 9, Number 2, Pages 449 457 c 2012 Institute for Scientific Computing and Information A STABILIZED NONCONFORMING QUADRILATERAL FINITE 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 informationProblems of Corner Singularities
Stuttgart 98 Problems of Corner Singularities Monique DAUGE Institut de Recherche MAthématique de Rennes Problems of Corner Singularities 1 Vertex and edge singularities For a polyhedral domain Ω and a
More informationA Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations
A Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations Songul Kaya and Béatrice Rivière Abstract We formulate a subgrid eddy viscosity method for solving the steady-state
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
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 informationAn Equal-order DG Method for the Incompressible Navier-Stokes Equations
An Equal-order DG Method for the Incompressible Navier-Stokes Equations Bernardo Cockburn Guido anschat Dominik Schötzau Journal of Scientific Computing, vol. 40, pp. 188 10, 009 Abstract We introduce
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 informationFinite Element Error Estimates in Non-Energy Norms for the Two-Dimensional Scalar Signorini Problem
Journal manuscript No. (will be inserted by the editor Finite Element Error Estimates in Non-Energy Norms for the Two-Dimensional Scalar Signorini Problem Constantin Christof Christof Haubner Received:
More informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 2198-5855 On the divergence constraint in mixed finite element methods for incompressible
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
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 informationA posteriori error estimates for non conforming approximation of eigenvalue problems
A posteriori error estimates for non conforming approximation of eigenvalue problems E. Dari a, R. G. Durán b and C. Padra c, a Centro Atómico Bariloche, Comisión Nacional de Energía Atómica and CONICE,
More informationc 2008 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 46, No. 3, pp. 640 65 c 2008 Society for Industrial and Applied Mathematics A WEIGHTED H(div) LEAST-SQUARES METHOD FOR SECOND-ORDER ELLIPTIC PROBLEMS Z. CAI AND C. R. WESTPHAL
More informationFINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
Proceedings of ALGORITMY pp. 9 3 FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION JAN STEBEL Abstract. The paper deals with the numerical simulations of steady flows
More informationConferences Fictitious domain methods for numerical realization of unilateral problems. J.H. (jointly with T. Kozubek and R. Kučera) p.
p. /2 Conferences 28 Fictitious domain methods for numerical realization of unilateral problems J.H. (jointly with T. Kozubek and R. Kučera) Lyon 23-24.6. 28 p. 2/2 Outline Fictitious domain methods -
More informationAxioms of Adaptivity (AoA) in Lecture 3 (sufficient for optimal convergence rates)
Axioms of Adaptivity (AoA) in Lecture 3 (sufficient for optimal convergence rates) Carsten Carstensen Humboldt-Universität zu Berlin 2018 International Graduate Summer School on Frontiers of Applied and
More informationMIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS
MIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS JUHO KÖNNÖ, DOMINIK SCHÖTZAU, AND ROLF STENBERG Abstract. We derive new a-priori and a-posteriori error estimates for mixed nite
More informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
More informationWEAK GALERKIN FINITE ELEMENT METHODS ON POLYTOPAL MESHES
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 12, Number 1, Pages 31 53 c 2015 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHODS ON POLYOPAL MESHES LIN
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationHigh order Galerkin appoximations for parametric second order elliptic partial differential equations
Research Collection Report High order Galerkin appoximations for parametric second order elliptic partial differential equations Author(s): Nistor, Victor; Schwab, Christoph Publication Date: 2012 Permanent
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 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 informationImplementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs
Implementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs Roman Andreev ETH ZÜRICH / 29 JAN 29 TOC of the Talk Motivation & Set-Up Model Problem Stochastic Galerkin FEM Conclusions & Outlook Motivation
More informationTomás P. Barrios 1, Gabriel N. Gatica 2,María González 3 and Norbert Heuer Introduction
ESAIM: M2AN Vol. 40, N o 5, 2006, pp. 843 869 DOI: 10.1051/m2an:2006036 ESAIM: Mathematical Modelling Numerical Analysis www.edpsciences.org/m2an A RESIDUAL BASED A POSERIORI ERROR ESIMAOR FOR AN AUGMENED
More informationLocal pointwise a posteriori gradient error bounds for the Stokes equations. Stig Larsson. Heraklion, September 19, 2011 Joint work with A.
Local pointwise a posteriori gradient error bounds for the Stokes equations Stig Larsson Department of Mathematical Sciences Chalmers University of Technology and University of Gothenburg Heraklion, September
More informationMedius analysis and comparison results for first-order finite element methods in linear elasticity
IMA Journal of Numerical Analysis Advance Access published November 7, 2014 IMA Journal of Numerical Analysis (2014) Page 1 of 31 doi:10.1093/imanum/dru048 Medius analysis and comparison results for first-order
More informationSpace time finite and boundary element methods
Space time finite and boundary element methods Olaf Steinbach Institut für Numerische Mathematik, TU Graz http://www.numerik.math.tu-graz.ac.at based on joint work with M. Neumüller, H. Yang, M. Fleischhacker,
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 informationAdaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation
Adaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA www.math.umd.edu/ rhn 7th
More informationA priori error analysis of the BEM with graded meshes for the electric eld integral equation on polyhedral surfaces
A priori error analysis of the BEM with graded meshes for the electric eld integral equation on polyhedral surfaces A. Bespalov S. Nicaise Abstract The Galerkin boundary element discretisations of the
More informationIMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1
Computational Methods in Applied Mathematics Vol. 1, No. 1(2001) 1 8 c Institute of Mathematics IMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1 P.B. BOCHEV E-mail: bochev@uta.edu
More informationarxiv: v1 [math.na] 13 Sep 2016
THE PENALTY FREE NITSCHE METHOD AND NONCONFORMING FINITE ELEMENTS FOR THE SIGNORINI PROBLEM ERIK BURMAN, PETER HANSBO, AND MATS G. LARSON arxiv:1609.03745v1 [math.na] 13 Sep 2016 Abstract. We design and
More informationDiscontinuous Galerkin Time Discretization Methods for Parabolic Problems with Linear Constraints
J A N U A R Y 0 1 8 P R E P R N T 4 7 4 Discontinuous Galerkin Time Discretization Methods for Parabolic Problems with Linear Constraints gor Voulis * and Arnold Reusken nstitut für Geometrie und Praktische
More informationMIXED FINITE ELEMENT APPROXIMATION OF THE VECTOR LAPLACIAN WITH DIRICHLET BOUNDARY CONDITIONS
MIXED FINITE ELEMENT APPROXIMATION OF THE VECTOR LAPLACIAN WITH DIRICHLET BOUNDARY CONDITIONS DOUGLAS N. ARNOLD, RICHARD S. FALK, AND JAY GOPALAKRISHNAN Abstract. We consider the finite element solution
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 informationTechnische Universität Graz
Technische Universität Graz A note on the stable coupling of finite and boundary elements O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2009/4 Technische Universität Graz A
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 informationControllability of the linear 1D wave equation with inner moving for
Controllability of the linear D wave equation with inner moving forces ARNAUD MÜNCH Université Blaise Pascal - Clermont-Ferrand - France Toulouse, May 7, 4 joint work with CARLOS CASTRO (Madrid) and NICOLAE
More informationTD 1: Hilbert Spaces and Applications
Université Paris-Dauphine Functional Analysis and PDEs Master MMD-MA 2017/2018 Generalities TD 1: Hilbert Spaces and Applications Exercise 1 (Generalized Parallelogram law). Let (H,, ) be a Hilbert space.
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 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 informationApplied/Numerical Analysis Qualifying Exam
Applied/Numerical Analysis Qualifying Exam August 9, 212 Cover Sheet Applied Analysis Part Policy on misprints: The qualifying exam committee tries to proofread exams as carefully as possible. Nevertheless,
More informationCoupling of HDG with a double-layer potential BEM
Coupling of HDG with a double-layer potential BEM Zhixing Fu, Norbert Heuer, & Francisco Javier Sayas June 29, 2018 arxiv:1212.2128v1 [math.na] 10 Dec 2012 Abstract In this paper we propose and analyze
More informationTime-dependent Dirichlet Boundary Conditions in Finite Element Discretizations
Time-dependent Dirichlet Boundary Conditions in Finite Element Discretizations Peter Benner and Jan Heiland November 5, 2015 Seminar Talk at Uni Konstanz Introduction Motivation A controlled physical processes
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 informationOrthogonality of hat functions in Sobolev spaces
1 Orthogonality of hat functions in Sobolev spaces Ulrich Reif Technische Universität Darmstadt A Strobl, September 18, 27 2 3 Outline: Recap: quasi interpolation Recap: orthogonality of uniform B-splines
More informationThe 2D Magnetohydrodynamic Equations with Partial Dissipation. Oklahoma State University
The 2D Magnetohydrodynamic Equations with Partial Dissipation Jiahong Wu Oklahoma State University IPAM Workshop Mathematical Analysis of Turbulence IPAM, UCLA, September 29-October 3, 2014 1 / 112 Outline
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 information