ngsxfem: Geometrically unfitted discretizations with Netgen/NGSolve (
|
|
- Clemence Walton
- 5 years ago
- Views:
Transcription
1 ngsxfem: Geometrically unfitted discretizations with Netgen/NGSolve ( Christoph Lehrenfeld (with contributions from F. Heimann, J. Preuß, M. Hochsteger,...) NGSolve user meeting, Vienna, June 15, 2017 Institute for Numerical and Applied Mathematics, University of Göttingen
2 Table of contents Example problems Requirements for geometrically unfitted FEM Numerical integration on cut elements Jupyter demo Conclusion & Outlook
3 Geometry description via level sets Level set function φ = 0, x Γ(t), φ(x, t) < 0, x Ω 1(t), > 0, x Ω 2(t). Problems: PDEs on interior domain Ω 1 Ω 2 Ω 1 Γ interface problems involving Ω 1, Ω 2 and Γ PDEs on the interface Γ description only implicit Properties geometry evolution with PDE for φ meshes are unfitted C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
4 Example problem 1: Stationary fictitious domain problem Stationary geometrically unfitted domain (Γ(t) = Γ): Features: α u = f in Ω = {φ < 0}, α u n = g N on Γ N {φ = 0}, u = g D on Γ D {φ = 0}. CutFEM : Finite element space V Γ h = V h Ω Imposition of Dirichlet values = Lagrange multiplier / Nitsche techniques Weak formulations require integrals on Ω dx, Γ dx: Numerical integration on cut elements Small cuts introduce stability and conditioning problems (stabilizations?!) C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
5 Example problem 2: Stationary interface problem Stationary geometrically unfitted domain (Γ(t) = Γ): Features: α u = f in Ω 1 Ω 2, [[ α u n]] = 0 on Γ, [[u]] = 0 on Γ, u = 0 on Ω. Solution has kinks across Γ: CutFEM : Vh Γ = V h Ω1 V h Ω2 or XFEM : Vh Γ = V h Vh x Imposition of interface conditions = Lagrange multiplier / Nitsche techniques Numerical integration on cut elements Small cuts (stability / conditioning) C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) 3 20 %
6 Example problem 3: Stationary surface PDE problem (Laplace Beltrami) Stationary geometrically unfitted domain (Γ(t) = Γ): Features: Γ u = f on Γ. Solution defined only on Γ Trace/CutFEM : Vh Γ = V h Γ Numerical integration on cut elements Approximation of normal/tangential direction Ambiuigity in normal directions (near) singular matrices (conditioning) C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
7 Example problem 4: Unsteady Poisson problem Geometrically unfitted moving domain (Γ(t) = {φ(x, t) = 0}): tu u = f in Ω(t) = {φ(x, t) < 0}, α u n = g N on Γ N (t) {φ(x, t) = 0}, u = g D on Γ D (t) {φ(x, t) = 0}, u = u 0 in Ω(0) = {φ(x, 0) < 0}. t Features: Time-dependent finite element spaces Space-Time: Space-time finite element spaces Numerical integration on space-time cut elements Characteristic FEM: Characteristic back-tracking integrator + features of stationary unf. problems... Γ x C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
8 Remarks Remarks on Examples 1-4 Appear in two-phase flows Vector-valued versions relevant (Stokes interface, mean curvature calculation,..) Background discretization can be conforming or nonconforming C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) 6 40 %
9 Existing features Related existing features in NGSolve: Background finite element spaces (scalar/vector, continuous/discontinuous) Convenient integral forms (VOL/BND/BBND) Handling of additional dofs: UNUSED DOFs Easy set up of preconditioners Curved elements (isoparametric FE) Restriction of integration on some elements... bfi = SymbolicBFI( u*v, BND) bfi.setdefinedonelements(..) C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
10 The new features of ngsxfem New features in ngsxfem: CutInfo (level set cut topology): Domain type per (VOL/BND) element: NEG: completely in {φ < 0} POS: completely in {φ > 0} IF: interseceted by {φ = 0} Cut ratio T Ω / T (detect small cuts) GetDofsOfElements(Vh,neg) dof-handling for unfitted FE spaces, V h = V h Ω Extended FE spaces, V h = V h V x h Ghost penalty stabilization: Marking of relevant facets (as BitArray) Differential operator for [[ l nu]] across facets deformed mesh Θh(Th) quadrilaterals triangles Numerical integration on cut elements undeformed mesh Th deformation order: 2 Θh (Vh 2)2 integration order: 4 C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
11 Numerical integration on cut simplices Tesselation Approximation by piecewise linear level set function = piecewise linear interface Γ lin + explicit domain approx. + robust 2nd order x 3 c 3 x 4 c 1 c 2 x 1 x 2 x 3 c 3 x 1 x 2 c 3 c 2 x 1 x 2 x 1 c 3 x 4 c 1 c 2 c 3 c 1 c 2 T 1 = T a T b T c T 2 = T d C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) 9 60 %
12 Numerical integration on cut simplices Tesselation Approximation by piecewise linear level set function = piecewise linear interface Γ lin + explicit domain approx. + robust 2nd order x 3 x 4 c 3 c 3 c 3 c 2 c 1 x 2 x 1 c 2 x 3 x 2 x 2 x 1 x 1 x 1 c 3 x 4 c 3 c 2 c 1 c 2 c 1 T 1 = T a T b T c T 2 = T d fdx or uvdx in ngsxfem: Ω lin Ω lin Integrate(levelset domain={"levelset":lsetp1,"domain type":neg},cf=f) SymbolicBFI(levelset domain={"levelset":lsetp1,"domain type":neg},u*v) C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) 9 60 %
13 Parametric mapping for higher order geometrical accuracy + higher order, implicit low order, explicit Basic concept of isoparametric unfitted FEM Start from (multi-)linear level set I h φ (h) C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
14 Parametric mapping for higher order geometrical accuracy + Θ h higher order, implicit low order, explicit Basic concept of isoparametric unfitted FEM higher order, explicit Start from (multi-)linear level set I h φ (h) Construct a mapping of the underlying mesh s.t. I h φ φ Θ h C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
15 Parametric mapping for higher order geometrical accuracy + Θ h higher order, implicit low order, explicit Basic concept of isoparametric unfitted FEM higher order, explicit Start from (multi-)linear level set I h φ (h) Construct a mapping of the underlying mesh s.t. I h φ φ Θ h lsetmeshadap=levelsetmeshadaptation(mesh,order=order) deformation=lsetmeshadap.calcdeformation(levelset) Θ h mesh.setdeformation(deformation) C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
16 Examples C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
17 Numerical integration in ngsxfem Numerical integration on cut elements with φ lin + mesh transformation Θ h Integration on linear level set domains on simplices: Geometrical decomposition into simplices Integration by transformation to reference domain (Ω lin ) C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
18 Numerical integration in ngsxfem Numerical integration on cut elements with φ lin + mesh transformation Θ h Integration on linear level set domains on simplices: Geometrical decomposition into simplices Integration by transformation to reference domain (Ω lin ) Transformation factors (det(dθ h ),...) are automatically considered due to NGSolves mesh deformation handling (ALETransformation) Changes in quadrature (Ω lin i Ω i,h = Θ h(ω lin i )) Quadrature after mapping, dist( Ω i, (Ω i,h)) O(h k+1 ), ω i > 0: f dx f dx = f dx ω i det(dθ h(x i)) f (Θ h(x i)) Ω i Ω i,h Θ h (Ω lin ) i T T h i C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
19 Numerical integration in ngsxfem Numerical integration on cut elements with φ lin + mesh transformation Θ h Integration on linear level set domains on simplices: Geometrical decomposition into simplices Integration by transformation to reference domain (Ω lin ) Transformation factors (det(dθ h ),...) are automatically considered due to NGSolves mesh deformation handling (ALETransformation) Changes in quadrature (Ω lin i Ω i,h = Θ h(ω lin i )) Quadrature after mapping, dist( Ω i, (Ω i,h)) O(h k+1 ), ω i > 0: f dx f dx = f dx ω i det(dθ h(x i)) f (Θ h(x i)) Ω i Ω i,h Θ h (Ω lin ) i accuracy depends on Θ h, but cut topology of Ω lin i T T h Consequences unchanged Guaranteed stability of quadrature: Ω i,h uv dx I h(ω i,h; uv), (positiveness of quadrature weights ω i det(dθ h(x i)) ) i C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
20 Mapping of quadrature points (order = 1) quadrilaterals triangles undeformed mesh Th deformed mesh Θh(Th) deformation order: 1 Θ h (V 1 h )2 integration order: 2 C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
21 Mapping of quadrature points (order = 2) quadrilaterals triangles undeformed mesh Th deformed mesh Θh(Th) deformation order: 2 Θ h (V 2 h )2 integration order: 4 C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
22 Mapping of quadrature points (order = 3) quadrilaterals triangles undeformed mesh Th deformed mesh Θh(Th) deformation order: 3 Θ h (V 3 h )2 integration order: 6 C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
23 Mapping of quadrature points (order = 4) quadrilaterals triangles undeformed mesh Th deformed mesh Θh(Th) deformation order: 4 Θ h (V 4 h )2 integration order: 8 C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
24 Jupyter demos! C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
25 Conclusion Θ Ω h boundary approximation Ω h Ω h Demands from Ω h unfitted FEM Φ Ω h id Φ Ω h Handling cut information (cut elements...) Φ Ω h =id Numerical integration in cut elements Stabilization on parts of the mesh ch interface approximation Ω h Ω Ω 2 Ω 1 Ω Γ Θ Γ h Ω 1,h Ω 2,h Γ h Φ Γ h =id Tools in ngsxfem Ω 1,h Ω ch 1 CutInfo Φ Γ new FESpaces h Φ Γ h id (XFESpace, SpaceTimeFESpace,...) Levelset-domain integration Γ h for SymbolicBFI/LFI and Integrate Ω Γ explicit 2,h high Ω 2 Φ Γ h =id order geometry approximation with trafo Θ h Possibilities High order accurate unfitted FEM Simplex and quad meshes Fictitious domain / interface problems / surface PDEs Scalar and vector-valued problems Time dependent problem:... C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
26 Outlook: Numerical integration in space-time in ngsxfem (poster) Space time integration Space-time FE for φ lin ( order k t in time, order 1 in space) Iterated integrals with space-time prism (in regions with same cut topology) Combine with space-time finite element for Θ h : Ω [t n 1, t n ] Ω T +(t n ) t 3 = t n T (tn ) Θ h Θ h(t +(t n )) t 3 = t n Θ h(t (t n )) t 2 t 1 t 2 t 1 tfe = ScalarTimeFE(k_t) t 0 = t n 1 t 0 = t n 1 fes = H1(mesh, order=1) Q T st_fes = SpaceTimeFESpace(fes,tfe) lsetp1 = GridFunction(st_fes)... SymbolicBFI(levelset_domain={"levelset":lsetp1,..},u*v,time_order=4) Θ h (Q T ) C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
27 Thank you for your attention! Questions / Comments? C. Lehrenfeld, ngsxfem ( June 15, 2017, Vienna (NGSolve User Meeting 2017) %
Space-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
More informationHigh Order Unfitted Finite Element Methods for Interface Problems and PDEs on Surfaces
N O V E M B E R 2 0 1 6 P R E P R I N T 4 5 9 High Order Unfitted Finite Element Methods for Interface Problems and PDEs on Surfaces Christoph Lehrenfeld and Arnold Reusken Institut für Geometrie und Praktische
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 informationAnalysis of a high order trace finite element method for PDEs on level set surfaces
N O V E M B E R 2 0 1 6 P R E P R I N T 4 5 7 Analysis of a high order trace finite element method for PDEs on level set surfaces Jörg Grande *, Christoph Lehrenfeld and Arnold Reusken * Institut für Geometrie
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 informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationA 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 informationHybrid Discontinuous Galerkin methods for incompressible flow problems
Hybrid Discontinuous Galerkin methods incompressible flow problems Christoph Lehrenfeld, Joachim Schöberl MathCCES Computational Mathematics in Engineering Workshop Linz, May 31 - June 1, 2010 Contents
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 informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationA short introduction to FreeFem++
A short introduction to FreeFem++ M. Ersoy Basque Center for Applied Mathematics 12 November 2010 Outline of the talk Outline of the talk 1 Introduction 2 HOW TO solve steady pdes : e.g. Laplace equation
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 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 informationSpectral Processing. Misha Kazhdan
Spectral Processing Misha Kazhdan [Taubin, 1995] A Signal Processing Approach to Fair Surface Design [Desbrun, et al., 1999] Implicit Fairing of Arbitrary Meshes [Vallet and Levy, 2008] Spectral Geometry
More informationFEM Convergence for PDEs with Point Sources in 2-D and 3-D
FEM Convergence for PDEs with Point Sources in 2-D and 3-D Kourosh M. Kalayeh 1, Jonathan S. Graf 2 Matthias K. Gobbert 2 1 Department of Mechanical Engineering 2 Department of Mathematics and Statistics
More informationFinite Element Techniques for the Numerical Simulation of Two-Phase Flows with Mass Transport
Finite Element Techniques for the Numerical Simulation of Two-Phase Flows with Mass Transport Christoph Lehrenfeld and Arnold Reusken Preprint No. 413 December 2014 Key words: Two-phase flow, mass transport,
More informationHybrid (DG) Methods for the Helmholtz Equation
Hybrid (DG) Methods for the Helmholtz Equation Joachim Schöberl Computational Mathematics in Engineering Institute for Analysis and Scientific Computing Vienna University of Technology Contributions by
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 VOLUME MESH FINITE ELEMENT METHOD FOR PDES ON SURFACES
European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS 2012) J. Eberhardsteiner et.al. (eds.) Vienna, Austria, September 10-14, 2012 A VOLUME MESH FINIE ELEMEN MEHOD FOR
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 25: Introduction to Discontinuous Galerkin Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods
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 informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More 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 informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationNumerical analysis and a-posteriori error control for a new nonconforming quadrilateral linear finite element
Numerical analysis and a-posteriori error control for a new nonconforming quadrilateral linear finite element M. Grajewski, J. Hron, S. Turek Matthias.Grajewski@math.uni-dortmund.de, jaroslav.hron@math.uni-dortmund.de,
More informationNumerical Integration over Implicitly Defined Domains for Higher Order Unfitted Finite Element Methods
ISSN 1995-0802, Lobachevskii Journal of Mathematics, 2016, Vol. 37, No. 5, pp. 582 596. c Pleiades Publishing, Ltd., 2016. Numerical Integration over Implicitly Defined Domains for Higher Order Unfitted
More informationNitsche-type Mortaring for Maxwell s Equations
Nitsche-type Mortaring for Maxwell s Equations Joachim Schöberl Karl Hollaus, Daniel Feldengut Computational Mathematics in Engineering at the Institute for Analysis and Scientific Computing Vienna University
More informationA very short introduction to the Finite Element Method
A very short introduction to the Finite Element Method Till Mathis Wagner Technical University of Munich JASS 2004, St Petersburg May 4, 2004 1 Introduction This is a short introduction to the finite element
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 informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationMoving interface problems. for elliptic systems. John Strain Mathematics Department UC Berkeley June 2010
Moving interface problems for elliptic systems John Strain Mathematics Department UC Berkeley June 2010 1 ALGORITHMS Implicit semi-lagrangian contouring moves interfaces with arbitrary topology subject
More informationarxiv: v1 [math.na] 13 Aug 2014
THE NITSCHE XFEM-DG SPACE-TIME METHOD AND ITS IMPLEMENTATION IN THREE SPACE DIMENSIONS CHRISTOPH LEHRENFELD arxiv:1408.2941v1 [math.na] 13 Aug 2014 Abstract. In the recent paper [C. Lehrenfeld, A. Reusken,
More informationarxiv: v3 [math.na] 25 Jun 2017
ANALYSIS OF A HIGH ORDER UNFITTED FINITE ELEMENT METHOD FOR ELLIPTIC INTERFACE PROBLEMS CHRISTOPH LEHRENFELD AND ARNOLD REUSKEN arxiv:1602.02970v3 [mat.na] 25 Jun 2017 Abstract. In te context of unfitted
More informationConvergence Order Studies for Elliptic Test Problems with COMSOL Multiphysics
Convergence Order Studies for Elliptic Test Problems with COMSOL Multiphysics Shiming Yang and Matthias K. Gobbert Abstract. The convergence order of finite elements is related to the polynomial order
More informationMaximum-norm a posteriori estimates for discontinuous Galerkin methods
Maximum-norm a posteriori estimates for discontinuous Galerkin methods Emmanuil Georgoulis Department of Mathematics, University of Leicester, UK Based on joint work with Alan Demlow (Kentucky, USA) DG
More informationFinite Element Method (FEM)
Finite Element Method (FEM) The finite element method (FEM) is the oldest numerical technique applied to engineering problems. FEM itself is not rigorous, but when combined with integral equation techniques
More informationHigh-order finite volume method for curved boundaries and non-matching domain-mesh problems
High-order finite volume method for curved boundaries and non-matching domain-mesh problems May 3-7, 0, São Félix, Portugal Ricardo Costa,, Stéphane Clain,, Gaspar J. Machado, Raphaël Loubère Institut
More informationThe Discontinuous Galerkin Finite Element Method
The Discontinuous Galerkin Finite Element Method Michael A. Saum msaum@math.utk.edu Department of Mathematics University of Tennessee, Knoxville The Discontinuous Galerkin Finite Element Method p.1/41
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
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 informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
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 informationBook of Abstracts - Extract th Annual Meeting. March 23-27, 2015 Lecce, Italy. jahrestagung.gamm-ev.de
GESELLSCHAFT für ANGEWANDTE MATHEMATIK und MECHANIK e.v. INTERNATIONAL ASSOCIATION of APPLIED MATHEMATICS and MECHANICS 86 th Annual Meeting of the International Association of Applied Mathematics and
More informationSome examples in domain decomposition
Some examples in domain decomposition Frédéric Nataf nataf@ann.jussieu.fr, www.ann.jussieu.fr/ nataf Laboratoire J.L. Lions, CNRS UMR7598. Université Pierre et Marie Curie, France Joint work with V. Dolean
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationAdvanced Higher Order Finite Elements From Theory to Application with Netgen/NGSolve
Advanced Higher Order Finite Elements From Theory to Application with Netgen/NGSolve Christoph Lehrenfeld, Joachim Schöberl Center for Computational Engineering Sciences (CCES) RWTH Aachen University Vancouver
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 informationSTRUCTURAL OPTIMIZATION BY THE LEVEL SET METHOD
Shape optimization 1 G. Allaire STRUCTURAL OPTIMIZATION BY THE LEVEL SET METHOD G. ALLAIRE CMAP, Ecole Polytechnique with F. JOUVE and A.-M. TOADER 1. Introduction 2. Setting of the problem 3. Shape differentiation
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationFEM-Level Set Techniques for Multiphase Flow --- Some recent results
FEM-Level Set Techniques for Multiphase Flow --- Some recent results ENUMATH09, Uppsala Stefan Turek, Otto Mierka, Dmitri Kuzmin, Shuren Hysing Institut für Angewandte Mathematik, TU Dortmund http://www.mathematik.tu-dortmund.de/ls3
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 informationLibMesh Experience and Usage
LibMesh Experience and Usage John W. Peterson peterson@cfdlab.ae.utexas.edu and Roy H. Stogner roystgnr@cfdlab.ae.utexas.edu Univ. of Texas at Austin September 9, 2008 1 Introduction 2 Weighted Residuals
More informationDownloaded 11/21/16 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 38, No. 6, pp. A349 A3514 c 016 Society for Industrial and Applied Mathematics ROBUST PRECONDITIONING FOR XFEM APPLIED TO TIME-DEPENDENT STOKES PROBLEMS SVEN GROSS, THOMAS LUDESCHER,
More informationImplementation of the dg method
Implementation of the dg method Matteo Cicuttin CERMICS - École Nationale des Ponts et Chaussées Marne-la-Vallée March 6, 2017 Outline Model problem, fix notation Representing polynomials Computing integrals
More informationDiscrete Differential Geometry. Discrete Laplace-Beltrami operator and discrete conformal mappings.
Discrete differential geometry. Discrete Laplace-Beltrami operator and discrete conformal mappings. Technische Universität Berlin Geometric Methods in Classical and Quantum Lattice Systems, Caputh, September
More informationAn Introduction to Netgen and NGSolve Session 1 Using Netgen/NGSolve
An Introduction to Netgen and NGSolve Session 1 Using Netgen/NGSolve Christoph Lehrenfeld Center for Computational Engineering Sciences (CCES) RWTH Aachen University University of British Columbia, November
More informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
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 informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
More informationA generalized MBO diffusion generated motion for constrained harmonic maps
A generalized MBO diffusion generated motion for constrained harmonic maps Dong Wang Department of Mathematics, University of Utah Joint work with Braxton Osting (U. Utah) Workshop on Modeling and Simulation
More informationModelling of interfaces and free boundaries
University of Regensburg Regensburg, March 2009 Outline 1 Introduction 2 Obstacle problems 3 Stefan problem 4 Shape optimization Introduction What is a free boundary problem? Solve a partial differential
More informationMatrix square root and interpolation spaces
Matrix square root and interpolation spaces Mario Arioli and Daniel Loghin m.arioli@rl.ac.uk STFC-Rutherford Appleton Laboratory, and University of Birmingham Sparse Days,CERFACS, Toulouse, 2008 p.1/30
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 informationGeometric Multigrid Methods
Geometric Multigrid Methods Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University IMA Tutorial: Fast Solution Techniques November 28, 2010 Ideas
More informationNGsolve::Come to the edge
NGsolve::Come to the edge Facet spaces and hybrid methods Jay Gopalakrishnan Portland State University Winter 2015 Download code (works only on version 6.0) for these notes from here. Come to the edge,
More informationarxiv: v2 [math.na] 5 Jul 2017
Trace Finite Element Metods for PDEs on Surfaces Maxim A. Olsanskii and Arnold Reusken arxiv:1612.00054v2 [mat.na] 5 Jul 2017 Abstract In tis paper we consider a class of unfitted finite element metods
More information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More informationDISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS
DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS Sergey Korotov BCAM Basque Center for Applied Mathematics http://www.bcamath.org 1 The presentation is based on my collaboration with several
More informationMultigrid and Iterative Strategies for Optimal Control Problems
Multigrid and Iterative Strategies for Optimal Control Problems John Pearson 1, Stefan Takacs 1 1 Mathematical Institute, 24 29 St. Giles, Oxford, OX1 3LB e-mail: john.pearson@worc.ox.ac.uk, takacs@maths.ox.ac.uk
More informationOptimal preconditioners for Nitsche-XFEM discretizations of interface problems
Optimal preconditioners for Nitsche-XFEM discretizations of interface problems Christoph Lehrenfeld Arnold Reusken August 13, 2014 In the past decade, a combination of unfitted finite elements (or XFEM)
More information2. Polynomial interpolation
2. Polynomial interpolation Contents 2. POLYNOMIAL INTERPOLATION... 1 2.1 TYPES OF INTERPOLATION... 1 2.2 LAGRANGE ONE-DIMENSIONAL INTERPOLATION... 2 2.3 NATURAL COORDINATES*... 15 2.4 HERMITE ONE-DIMENSIONAL
More informationSpectral finite elements for a mixed formulation in computational acoustics taking flow effects into account
Spectral finite elements for a mixed formulation in computational acoustics taking flow effects into account Manfred Kaltenbacher in cooperation with A. Hüppe, I. Sim (University of Klagenfurt), G. Cohen
More informationDiscretization of PDEs and Tools for the Parallel Solution of the Resulting Systems
Discretization of PDEs and Tools for the Parallel Solution of the Resulting Systems Stan Tomov Innovative Computing Laboratory Computer Science Department The University of Tennessee Wednesday April 4,
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationA non-standard Finite Element Method based on boundary integral operators
A non-standard Finite Element Method based on boundary integral operators Clemens Hofreither Ulrich Langer Clemens Pechstein June 30, 2010 supported by Outline 1 Method description Motivation Variational
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationAnalysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems
Analysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems Christian Waluga 1 advised by Prof. Herbert Egger 2 Prof. Wolfgang Dahmen 3 1 Aachen Institute for Advanced Study in Computational
More informationLecture on: Numerical sparse linear algebra and interpolation spaces. June 3, 2014
Lecture on: Numerical sparse linear algebra and interpolation spaces June 3, 2014 Finite dimensional Hilbert spaces and IR N 2 / 38 (, ) : H H IR scalar product and u H = (u, u) u H norm. Finite dimensional
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 informationLibMesh Experience and Usage
LibMesh Experience and Usage John W. Peterson peterson@cfdlab.ae.utexas.edu Univ. of Texas at Austin January 12, 2007 1 Introduction 2 Weighted Residuals 3 Poisson Equation 4 Other Examples 5 Essential
More informationFEniCS Course. Lecture 0: Introduction to FEM. Contributors Anders Logg, Kent-Andre Mardal
FEniCS Course Lecture 0: Introduction to FEM Contributors Anders Logg, Kent-Andre Mardal 1 / 46 What is FEM? The finite element method is a framework and a recipe for discretization of mathematical problems
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 informationNew class of finite element methods: weak Galerkin methods
New class of finite element methods: weak Galerkin methods Xiu Ye University of Arkansas at Little Rock Second order elliptic equation Consider second order elliptic problem: a u = f, in Ω (1) u = 0, on
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
More informationApproximation of self-avoiding inextensible curves
Approximation of self-avoiding inextensible curves Sören Bartels Department of Applied Mathematics University of Freiburg, Germany Joint with Philipp Reiter (U Duisburg-Essen) & Johannes Riege (U Freiburg)
More informationJEPPIAAR ENGINEERING COLLEGE
JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III
More informationRobust preconditioning for XFEM applied to time-dependent Stokes problems
Numerical Analysis and Scientific Computing Preprint Seria Robust preconditioning for XFEM applied to time-dependent Stokes problems S. Gross T. Ludescher M.A. Olshanskii A. Reusken Preprint #40 Department
More informationAnalysis of a DG XFEM Discretization for a Class of Two Phase Mass Transport Problems
Analysis of a DG XFEM Discretization for a Class of Two Phase Mass Transport Problems Christoph Lehrenfeld and Arnold Reusken Bericht Nr. 340 April 2012 Key words: transport problem, Nitsche method, XFEM,
More informationFEM Convergence for PDEs with Point Sources in 2-D and 3-D
FEM Convergence for PDEs with Point Sources in -D and 3-D Kourosh M. Kalayeh 1, Jonathan S. Graf, and Matthias K. Gobbert 1 Department of Mechanical Engineering, University of Maryland, Baltimore County
More informationL 2 -Error Analysis of an Isoparametric Unfitted Finite Element Method for Elliptic Interface Problems
A P R I L 2 0 1 6 P R E P R I N T 4 5 0 L 2 -Error Analysis of an Isoparametric Unfitted Finite Element Metod for Elliptic Interface Problems Cristop Lerenfeld * and Arnold Reusken Institut für Geometrie
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationComputer simulation of multiscale problems
Progress in the SSF project CutFEM, Geometry, and Optimal design Computer simulation of multiscale problems Axel Målqvist and Daniel Elfverson University of Gothenburg and Uppsala University Umeå 2015-05-20
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 informationADI iterations for. general elliptic problems. John Strain Mathematics Department UC Berkeley July 2013
ADI iterations for general elliptic problems John Strain Mathematics Department UC Berkeley July 2013 1 OVERVIEW Classical alternating direction implicit (ADI) iteration Essentially optimal in simple domains
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationInstitut für Numerische und Angewandte Mathematik
Institut für Numerische und Angewandte Mathematik A note on Nitsche s method for unfitted boundary value problems de Prenter, F., Lehrenfeld, C., Massing, A. Nr. 8 Preprint-Serie des Instituts für Numerische
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
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 informationNitsche XFEM with Streamline Diffusion Stabilization for a Two Phase Mass Transport Problem
Nitsche XFEM with Streamline Diffusion Stabilization for a Two Phase Mass Transport Problem Christoph Lehrenfeld and Arnold Reusken Bericht Nr. 333 November 2011 Key words: transport problem, Nitsche method,
More information