An Interface-fitted Mesh Generator and Polytopal Element Methods for Elliptic Interface Problems
|
|
- Brett Wade
- 5 years ago
- Views:
Transcription
1 An Interface-fitted Mesh Generator and Polytopal Element Methods for Elliptic Interface Problems Long Chen University of California, Irvine Joint work with: Huayi Wei (Xiangtan University), Min Wen(UCI) Polytopal Element Methods in Mathematics and Engineering Oct 27, 2015, Atlanta, GA
2 Outline Elliptic Interface Problems Interface-fitted Mesh Generation 2D Algorithm 3D Algorithm VEM for Elliptic Interface Problems Weak Formulation Virtual Element Method Weak Galerkin Methods Numerical Results
3 ELLIPTIC INTERFACE PROBLEMS
4 A Domain with an Interface Ω Ω n Ω + Γ Figure: A square domain Ω with an interface Γ in it.
5 Elliptic Interface Problems Consider (β u) = f, in Ω\Γ (1) with prescribed jump conditions across the interface Γ: [u] Γ = u + u = q 0, [βu n ] Γ = β + u + n β u n = q 1, and boundary condition u = g on Ω.
6 Existing Work Two type of numerical methods: Numerical methods based on Cartesian meshes. Numerical methods based on interface-fitted meshes.
7 Cartesian Mesh Approach Figure: A Catesian mesh and an interface.
8 Cartesian Mesh Approach Pro: Mesh generation is very simple Con: Need to modify the stencil (FDM) or basis function (FEM) of the vertex near the interface Γ. The linear system may be non-symmetric and cause problems in fast solvers. Convergence analysis is complicated.
9 Interface-fitted Mesh Figure: An interface-fitted mesh in 2D [Wei, Chen, Huang and Zheng, SISC. 2014].
10 Interface-fitted Mesh Approach Pro: Can use standard finite element methods to discretize the interface problem. Symmetric system can be easily solved by fast solver such as algebraic multigrid solvers. Relatively easy error analysis. Con: Generate interface-fitted meshes, extremely difficult in 3D!
11 Our Goal In the past decade the mesh generation have gotten great progress. Here we aim to 1. Develop effective and robust interface-fitted (polytopal) mesh generation algorithm. 2. Solve interface problems accurately using conforming FEM (e.g. Virtual Element method) on polytopal meshes. 3. Solve the resulting algebraic system quickly using (algebraic) multigrid solver.
12 Outline Elliptic Interface Problems Interface-fitted Mesh Generation 2D Algorithm 3D Algorithm VEM for Elliptic Interface Problems Weak Formulation Virtual Element Method Weak Galerkin Methods Numerical Results
13 INTERFACE-FITTED MESH GENERATION TWO DIMENSIONAL CASE
14 Delaunay Triangulation p 6 p 7 p 3 p 5 p 1 p 4 p 2 Figure: A Delaunay triangulation of a set of points is a triangulation satisfying the empty circle condition.
15 Delaunay Triangulation (a) Non Delaunay. (b) Delaunay. Figure: Empty circle condition: a circle circumscribing any Delaunay triangle does not contain any other input points in its interior.
16 2D Interface-fitted Mesh Generation Algorithm Algorithm 1 2D Interface-fitted Mesh Generation Algorithm INPUT: Grid size: h; Level set function, ϕ(x); Square domain, Ω; OUTOUT: Interface-fitted mesh T ; 1: Find the cut points, the Cartesian mesh points near or on the interface and some auxiliary points. 2: Construct a Delaunay triangulation on these points. 3: Post processing.
17 An Example in 2D: Step 1 Figure: Step 1: Find the cut points (red), the mesh points near or on the interface (black) and the auxiliary points (magenta).
18 An Example in 2D: Step 2 Figure: Step 2: Generate a Delaunay triangulation on these points. In MATLAB, it is simply DT = delaunay(x, y).
19 An Example in 2D: Step 3 Figure: Step 3: Keep the triangles in the interface elements.
20 An Example in 2D: Step 4 Figure: An interface-fitted mesh composed by triangles and squares.
21 Features of 2D Mesh Generator Simple, efficient and semi-unstructured. Recovery the interface. The maximum angle is uniformly bounded (135 ).
22 146 C. P aum Maximum Maximal interior Angle angle: u max Condition 90 triangle Kh;2 are smaller or equal 135. To prove this, note that 45 a3 90 and a1 60. This implies a2 180 a1 a Maximal interior angle: u max 135 Subdivision of Type d) boundary subcell: Let B be a Type d) boundary subcell. Then X \ B can be approximated by four triangles Kh;1, Kh;2, Kh;3 and Kh;4 such that a2 a1. The interior angles of the triangles Kh;1, Kh;2 and Kh;4 are less than or equal to 90. The interior angles of the triangle Kh;3 are less than or equal to 135. Semi-Unstructured This can Gridsbe proven as in the Type c) 147 subcell case. (a) θ max 90 (b) θ max 135. Maximal interior angle: u max 135 Subdivision of Type b) boundary subcell: Subdivision Let Maximal B be interior aof Type angle: b) c) boundary u max subcell. subcell: 135 Let B be a Type c) boundary subcell. Then X \ B can be approximated by twothen triangles X \ BKh;1 can andbekh;2. approximated Obviously, by thethree triangles Kh;1, Kh;2 and Kh;3. The interior angles of triangle Kh;1 are less than interior equal angles to 90of. the Thetriangles interior angles Kh;1 and of Kh;3 are less than or equal to 90. The triangle Kh;2 are smaller or equal 135. To interior proveangles this, note of thethat triangle 45 Kh;2 a3 are 90less than or equal to 135. To prove this, and a1 60. This implies a2 180 a1 note a3that 135 a3. 90 and b This implies a1 180 b Analogously, we obtain a Maximal interior angle: u max 135 (c) θ max 135. (d) θ max 135. Subdivision of Type d) boundary subcell: Subdivision Let B be aof Type interior d) boundary cells: Every subcell. interior cell can be subdivided into two triangles. triangles Kh;1, Kh;2, Kh;3 and Kh;4 such Then X \ B can be approximated by four that a2 a1. The interior angles of the triangles Kh;1, Kh;2 and Kh;4 are less than or equal to 90. The interior angles of the triangle Kh;3 are less than or equal to 135. This can be proven as in the Type c) subcell case. Subdivision of Type c) boundary subcell: Let B be a Type c) boundary subcell. Then X \ B can be approximated by three triangles Kh;1, Kh;2 and Kh;3. The Figure: Maximum angles of four types of interface elements. From Semi-Unstructured Grids by C. Pflaum. Computing 2001.
23 Interface Recovery Semi-Unstructured Grids 145 Figure 1. Di erent types of boundary subcells Figure 2. Boundary subcell of a double type
24 Interface Recovery Semi-Unstructured Grids 145 Figure 1. Di erent types of boundary subcells We DO NOT code the triangulation case by case. We simply call DT = delaunay(x, y). Figure 2. Boundary subcell of a double type
25 Interface Recovery Semi-Unstructured Grids 145 Figure 1. Di erent types of boundary subcells We DO NOT code the triangulation case by case. We simply call DT = delaunay(x, y). We can PROVE the interface will be preserved in the triangulation generated by delaunay and the maximal angle is minimized. Figure 2. Boundary subcell of a double type
26 Lower Convex Hull (a) Step 1: Lift points to the paraboloid (b) Step 2: Form the lowest convex hull in R n+1. (c) Step 3: Project the lowest convex hull to R n. Figure: Delaunay triangulation is the projection of the lower convex hull of points lifted to the paraboloid f = x 2.
27 Interface Recovery Semi-Unstructured Grids 145 Figure 1. Di erent types of boundary subcells The lower convex hull when lift to R n+1 will always connect the intersection points except...
28 A Degenerate Case A B C D Figure: Add an auxiliary point in a rectangle to preserve the interface
29 Efficiency Simple, Efficient and Semi-Unstructured. The greatest advantage: localization. Only call delaunay for O( N) points near the interface. The complexity will be thus O( N) in 2-D which can be ignored comparing with the O(N) complexity for assembling the matrix and solving the matrix equation. No mesh smoothing is needed. The maximal angle is bounded by 135. Good enough for finite element methods. The mesh is only unstructured near the interface and the majority of the mesh is structured which leads to nice properties (superconvergence, multigrid etc).
30 INTERFACE-FITTED MESH GENERATION THREE DIMENSIONAL CASE
31 Challenges in Tetrahedra Mesh Generation Figure: Tetrahedra elements [Shewchuk, 2002]. Most Tetrahedra mesh generation/optimization algorithms cannot guarantee sliver-free for a general 3D domain. Advertisement: Use ODT mesh smoothing which produces fewer slivers in 3D. L. Chen. Mesh smoothing schemes based on optimal Delaunay triangulations. 2004
32 3D Interface-fitted Mesh Generation Algorithm Algorithm 2 3D Interface-fitted Mesh Generation Algorithm INPUT: Grid size: h; Level set function, ϕ(x); Cube domain, Ω; OUTOUT: Interface-fitted mesh T ; 1: Find the cut points, the Cartesian mesh points near or on the interface and the auxiliary points. 2: Construct a Delaunay triangulation on these points. 3: Post processing: use polytopal meshes.
33 Polytopal Mesh v.s. Tetrahedron Mesh Figure: Bad tetrahedron elements will be eliminated and part of their faces will become the boundary of polytopes. As a surface mesh, the triangles are more acceptable.
34 Features of 3D Mesh Generator Simple, efficient and semi-unstructured. Recovery the interface. The maximum angle of surface mesh is uniformly bounded. Boundary faces of a polytope is either triangle or square.
35 Sphere Interface Figure: Sphere surface with 3, 128 triangles: maximum angle and minimum angle
36 Sphere Interface Figure: The mesh near the interface. Use second to generate a mesh with 70, 169 points).
37 Heart interface Figure: Heart surface with 3480 triangles: maximum angle and minimum angle 4.21.
38 Heart interface Figure: The mesh near the interface. Use second to generate a mesh with 70, 521 points).
39 Quartics interface Figure: Quartics surface with 44, 512 triangles: maximum angle and minimum angle 11.2.
40 Quartics interface Figure: The mesh near the interface. Use second to generate a mesh with 553, 161 points.
41 Torus interface Figure: Torus surface with 37, 600 triangles: maximum angle and minimum angle
42 Torus interface Figure: The mesh near the interface. Use second to generate a mesh with 1, 045, 061 points.
43 Outline Elliptic Interface Problems Interface-fitted Mesh Generation 2D Algorithm 3D Algorithm VEM for Elliptic Interface Problems Weak Formulation Virtual Element Method Weak Galerkin Methods Numerical Results
44 Elliptic Interface Problems Consider (β u) = f, in Ω\Γ (2) with prescribed jump conditions across the interface Γ: [u] Γ = u + u = q 0, [βu n ] Γ = β + u + n β u n = q 1, and boundary condition u = g on Ω.
45 Weak Formulation: Case 1 For the elliptic interface problem, if the jump conditions are homogeneous, i.e., q 0 = 0 and q 1 = 0 on Γ, then the problem is equivalent to that of finding u H 1 (Ω) with u = g on Ω such that β u vdx = fv dx, v H0 1 (Ω) (3) Ω Ω
46 Weak Formulation: Case 2 If q 0 = 0, q 1 0 on Γ, then the corresponding weak formulation is : (β u, v) Ω = (f, v) Ω q 1, v Γ, v H 1 0 (Ω). (4) The jump condition [βu n ] Γ = q 1 holds in H 1/2 (Γ) sense.
47 Weak Formulation: Case 3 If q 0 0 on Γ, we can find w : Ω R with w = q 0 on Ω and w H 1 (Ω ). The zero extension of w,denote by w satisfies w H 1 (Ω Ω + ) with w = w on Ω and w = 0 on Ω +. The choice of w is not unique. The interface problem is equivalent to: find u = p w with p H0 1 (Ω) such that: (β p, v) Ω = (f, v) Ω q 1, v Γ + (β w, v) Ω, v H 1 0 (Ω).
48 Virtual Element Method Virtual element method is a generalization of the standard conforming FEM on polygon or polyhedra meshes. L. Beirão Da Veiga, F. Brezzi, A. Cangiani, G. Manzini, L. D. Marini, and A. Russo. Basic principles of virtual element methods. Mathematical Models and Methods in Applied Sciences, 23(1): , B. Ahmad, A. Alsaedi, F. Brezzi, L. D. Marini, and A. Russo. Equivalent projectors for virtual element methods. Computers and Mathematics with Applications, 66(3): , Beirão da Veiga L., Brezzi F., Marini L.D. and Russo A. The Hitchhiker s Guide to the Virtual Element Method. Math. Models Meth. Appl. Sci. 24: (2014).
49 Spaces in Virtual Element Methods We introduce the following space on a simple polygon K V k (K) := {v H 1 (K) : v K B k ( K), v P k 2 (K)}, where P k (D) is the space of polynomials of degree k on D and conventionally P 1 (D) = 0, and the boundary space B k ( K) := {v C 0 ( K) : v e P k (e) for all edges e K}. VEM vs??? (recall that N K = nk + k(k 1)/2) VEM k =1 k =2 k =3
50 Weak Galerkin Methods Given a polygon mesh T h and an integer k 1, we introduce W h = {v = {v 0, v b }, v 0 K P k (K), v b e P k 1 (e) e K, K T h }. Define a modified weak gradient K w : W k (K) P k (K) as: given v = {v 0, v b } W k (K), define K w v P k (K) such that ( K w v, p) K = (v 0, p) K + v b, n p K for all p P k (K). L. Mu, J. Wang, and X. Ye. A weak Galerkin finite element method with polynomial reduction. Journal of Computational and Applied Mathematics, 285:45 58, 2015.
51 Equivalence The WG finite element is: find u h W 0 h such that a WG (u h, v) = (f, v 0 ) v = {v 0, v b } Wh 0. where a WG (u h, v h ) := ( w u h, w v h )+ χ b (u b u 0 K ) χ b (v b v 0 K ). K T h The non-conforming VEM is: find u h Ṽ 0 h such that ā VEM h (u h, v h ) = (f, Π 0 k v h) v h Ṽ 0 h, where ā VEM h (u, v) := ( Π k u, Π k v) + χ b (u Π 0 k u) χ b(v Π 0 k v). K T h
52 Outline Elliptic Interface Problems Interface-fitted Mesh Generation 2D Algorithm 3D Algorithm VEM for Elliptic Interface Problems Weak Formulation Virtual Element Method Weak Galerkin Methods Numerical Results
53 NUMERICAL RESULTS
54 Implementation and Errors Our numerical test was based on MATLAB package ifem. L. Chen. ifem: an integrated finite element method package in MATLAB We consider the following errors ( ) u I u h A = β 1/2 h (u I u h ) 2 + β 1/2 1/2 0,Ω h (u+ I u+ h ) 2 h 0,Ω + h u I u h = max( u I u h ) 0,,Ω, u + I u+ h h ) 0,,Ω + ) h ( ) 1/2 u I u h 0 = u I u h ) 2 + u + 0,Ω I u+ h ) 2 h 0,Ω + h where u h is the numerical approximation and u I is the nodal interpolation.
55 A Cube Domain with a Sphere Interface We chose the following setting φ(x, y, z) =x 2 + y 2 + z 2 r 2 u + =10(x + y + z) u =5e x2 +y 2 +z where r = The coefficient β is piecewise constant.
56 A Cube Domain with a Sphere Surface
57 Numerical Tests: β = 1 and β + = 1 #Dof h u I u h A u I u h u I u h u I u h C 1 h Error Error 10 2 Du I Du h C 1 h u I u h C 2 h log(1/h) 10 1 log(1/h)
58 Numerical Tests: β = 1 and β + = 10 #Dof h u I u h A u I u h u I u h u I u h C 1 h Error Error 10 2 Du I Du h C 1 h u I u h C 2 h log(1/h) 10 1 log(1/h)
59 Numerical Tests:β = 1 and β + = 100 #Dof h u I u h A u I u h u I u h u I u h C 1 h Error Error 10 2 Du I Du h C 1 h u I u h C 2 h log(1/h) 10 1 log(1/h)
60 Numerical Tests:β = 1 and β + = 1000 #Dof h u I u h A u I u h u I u h u I u h Error 10 1 Error 10 2 C 1 h Du I Du h C 1 h u I u h C 2 h log(1/h) 10 1 log(1/h)
61 Numerical Tests: β = 1 and β + = #Dof h u I u h A u I u h u I u h u I u h Du I Du h 10 2 C 1 h Error 10 1 C 1 h u I u h Error C 2 h log(1/h) 10 1 log(1/h)
62 Algebraic Multigrid Solvers #Dof β + = 1 β + = 10 β + = 100 β + = 1000 β + = Table: Iteration steps of MGCG using W-cycle and AMG preconditioner with fixed β = 1 and increased β +.
63 Summary Simple, Efficient and Semi-Unstructured 2-D and 3D Interface-fitted Mesh Generator.
64 THANK YOU FOR YOUR ATTENTION!
65 THANK YOU FOR YOUR ATTENTION! Supported by NSF DMS
Tiling Three-Dimensional Space with Simplices. Shankar Krishnan AT&T Labs - Research
Tiling Three-Dimensional Space with Simplices Shankar Krishnan AT&T Labs - Research What is a Tiling? Partition of an infinite space into pieces having a finite number of distinct shapes usually Euclidean
More informationNodal Basis Functions for Serendipity Finite Elements
Nodal Basis Functions for Serendipity Finite Elements Andrew Gillette Department of Mathematics University of Arizona joint work with Michael Floater (University of Oslo) Andrew Gillette - U. Arizona Nodal
More information= β + u + n β u n = q 1,
ADAPTIVE MESH REFINEMENT AND SUPERCONVERGENCE FOR TWO DIMENSIONAL INTERFACE PROBLEMS HUAYI WEI, LONG CHEN, YUNQING HUANG, AND BIN ZHENG Abstract. Adaptive mesh refinement and the Börgers algorithm are
More informationOn a nested refinement of anisotropic tetrahedral grids under Hessian metrics
On a nested refinement of anisotropic tetrahedral grids under Hessian metrics Shangyou Zhang Abstract Anisotropic grids, having drastically different grid sizes in different directions, are efficient and
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 informationLecture 2 Unstructured Mesh Generation
Lecture 2 Unstructured Mesh Generation MIT 16.930 Advanced Topics in Numerical Methods for Partial Differential Equations Per-Olof Persson (persson@mit.edu) February 13, 2006 1 Mesh Generation Given a
More informationContents. I The Basic Framework for Stationary Problems 1
page v Preface xiii I The Basic Framework for Stationary Problems 1 1 Some model PDEs 3 1.1 Laplace s equation; elliptic BVPs... 3 1.1.1 Physical experiments modeled by Laplace s equation... 5 1.2 Other
More informationMöbius Transformations in Scientific Computing. David Eppstein
Möbius Transformations in Scientific Computing David Eppstein Univ. of California, Irvine School of Information and Computer Science (including joint work with Marshall Bern from WADS 01 and SODA 03) Outline
More informationThe Immersed Interface Method
The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial
More informationTopological Issues in Hexahedral Meshing
Topological Issues in Hexahedral Meshing David Eppstein Univ. of California, Irvine Dept. of Information and Computer Science Outline I. What is meshing? Problem statement Types of mesh Quality issues
More informationWESTFÄLISCHE WILHELMS-UNIVERSITÄT MÜNSTER
MÜNSTER Adaptive Discretization of Liftings for Curvature Regularization SIAM Conference on Imaging Science 2016, Albuquerque, NM Ulrich Hartleif, Prof. Dr. Benedikt Wirth May 23, 2016 Outline MÜNSTER
More informationarxiv: v1 [math.na] 20 Sep 2016
arxiv:1609.06236v1 [math.na] 20 Sep 2016 A Local Mesh Modification Strategy for Interface Problems with Application to Shape and Topology Optimization P. Gangl 1,2 and U. Langer 3 1 Doctoral Program Comp.
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More information13.472J/1.128J/2.158J/16.940J COMPUTATIONAL GEOMETRY
13.472J/1.128J/2.158J/16.940J COMPUTATIONAL GEOMETRY Lecture 23 Dr. W. Cho Prof. N. M. Patrikalakis Copyright c 2003 Massachusetts Institute of Technology Contents 23 F.E. and B.E. Meshing Algorithms 2
More informationSerendipity Nodal VEM spaces
Serendipity Nodal VEM spaces L. Beirão da Veiga, F. Brezzi, L.D. Marini, A. Russo Lourenço Beirão da Veiga - Dipartimento di Matematica e Applicazioni, Università di Milano Bicocca, Via Cozzi 53, I-20153,
More informationConforming Vector Interpolation Functions for Polyhedral Meshes
Conforming Vector Interpolation Functions for Polyhedral Meshes Andrew Gillette joint work with Chandrajit Bajaj and Alexander Rand Department of Mathematics Institute of Computational Engineering and
More informationGEOMETRY MODELING & GRID GENERATION
GEOMETRY MODELING & GRID GENERATION Dr.D.Prakash Senior Assistant Professor School of Mechanical Engineering SASTRA University, Thanjavur OBJECTIVE The objectives of this discussion are to relate experiences
More informationVoroCrust: Simultaneous Surface Reconstruction and Volume Meshing with Voronoi cells
VoroCrust: Simultaneous Surface Reconstruction and Volume Meshing with Voronoi cells Scott A. Mitchell (speaker), joint work with Ahmed H. Mahmoud, Ahmad A. Rushdi, Scott A. Mitchell, Ahmad Abdelkader
More informationOutline of the presentation
Surface Reconstruction Petra Surynková Charles University in Prague Faculty of Mathematics and Physics petra.surynkova@mff.cuni.cz Outline of the presentation My work up to now Surfaces of Building Practice
More informationComputational Geometry
Computational Geometry 600.658 Convexity A set S is convex if for any two points p, q S the line segment pq S. S p S q Not convex Convex? Convexity A set S is convex if it is the intersection of (possibly
More informationTHE MORTAR FINITE ELEMENT METHOD IN 2D: IMPLEMENTATION IN MATLAB
THE MORTAR FINITE ELEMENT METHOD IN D: IMPLEMENTATION IN MATLAB J. Daněk, H. Kutáková Department of Mathematics, University of West Bohemia, Pilsen MECAS ESI s.r.o., Pilsen Abstract The paper is focused
More informationSOLVING PARTIAL DIFFERENTIAL EQUATIONS ON POINT CLOUDS
SOLVING PARTIAL DIFFERENTIAL EQUATIONS ON POINT CLOUDS JIAN LIANG AND HONGKAI ZHAO Abstract. In this paper we present a general framework for solving partial differential equations on manifolds represented
More informationWe have set up our axioms to deal with the geometry of space but have not yet developed these ideas much. Let s redress that imbalance.
Solid geometry We have set up our axioms to deal with the geometry of space but have not yet developed these ideas much. Let s redress that imbalance. First, note that everything we have proven for the
More informationarxiv: v1 [math.na] 28 Oct 2015
Serendipity Nodal VEM spaces L. Beirão da Veiga, F. Brezzi, L.D. Marini, A. Russo arxiv:1510.08477v1 [math.na] 28 Oct 2015 Lourenço Beirão da Veiga - Dipartimento di Matematica, Università di Milano Statale,
More informationADAPTIVE BEM-BASED FEM ON POLYGONAL MESHES FROM VIRTUAL ELEMENT METHODS
ECCOMAS Congress 6 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 June 6 ADAPTIVE
More informationSurface Mesh Generation
Surface Mesh Generation J.-F. Remacle Université catholique de Louvain September 22, 2011 0 3D Model For the description of the mesh generation process, let us consider the CAD model of a propeller presented
More informationDomain Decomposition and hp-adaptive Finite Elements
Domain Decomposition and hp-adaptive Finite Elements Randolph E. Bank 1 and Hieu Nguyen 1 Department of Mathematics, University of California, San Diego, La Jolla, CA 9093-011, USA, rbank@ucsd.edu. Department
More informationCHARMS: A Simple Framework for Adaptive Simulation SIGGRAPH Presented by Jose Guerra
CHARMS: A Simple Framework for Adaptive Simulation SIGGRAPH 2002 Eitan Grinspun Caltech Petr Krysl UCSD Peter Schröder Caltech Presented by Jose Guerra 1 Outline Background Motivation (Element vs. Basis
More informationFrom CAD surface models to quality meshes. Patrick LAUG. Projet GAMMA. INRIA Rocquencourt. Outline
From CAD surface models to quality meshes Patrick LAUG Projet GAMMA INRIA Rocquencourt Tetrahedron II Oct. 2007 1 Outline 1. Introduction B-Rep, patches 2. CAD repair ant topology recovery 3. Discretization
More informationTwo Algorithms for Adaptive Approximation of Bivariate Functions by Piecewise Linear Polynomials on Triangulations
Two Algorithms for Adaptive Approximation of Bivariate Functions by Piecewise Linear Polynomials on Triangulations Nira Dyn School of Mathematical Sciences Tel Aviv University, Israel First algorithm from
More informationComputational Geometry. Algorithm Design (10) Computational Geometry. Convex Hull. Areas in Computational Geometry
Computational Geometry Algorithm Design (10) Computational Geometry Graduate School of Engineering Takashi Chikayama Algorithms formulated as geometry problems Broad application areas Computer Graphics,
More informationNonoscillatory Central Schemes on Unstructured Triangulations for Hyperbolic Systems of Conservation Laws
Nonoscillatory Central Schemes on Unstructured Triangulations for Hyperbolic Systems of Conservation Laws Ivan Christov Bojan Popov Department of Mathematics, Texas A&M University, College Station, Texas
More informationNumerical representation of the quality measures of triangles and. triangular meshes
Numerical representation of the quality measures of triangles and triangular meshes J. Sarrate, J. Palau and A. Huerta Departament de Matemàtica Aplicada III, E.T.S. de Ingenieros de Caminos, Canales y
More informationInvestigation of a Robust Method for Connecting Dissimilar 3D Finite Element Models. David M. Trujillo 1. December 2005
Investigation of a Robust Method for Connecting Dissimilar 3D Finite Element Models by David M. Trujillo 1 December 2005 1 Consultant, TRUCOMP, Fountain Valley, California trucomp@earthlink.net Abstract
More informationA New Smoothing Algorithm for Quadrilateral and Hexahedral Meshes
A New Smoothing Algorithm for Quadrilateral and Hexahedral Meshes Sanjay Kumar Khattri Department of Mathematics, University of Bergen, Norway sanjay@mi.uib.no http://www.mi.uib.no/ sanjay Abstract. Mesh
More informationA comparative study between some bisection based partitions in 3D
Applied Numerical Mathematics 55 (2005) 357 367 www.elsevier.com/locate/apnum A comparative study between some bisection based partitions in 3D Miguel A. Padrón a,,josép.suárez b, Ángel Plaza c a Department
More informationScientific Computing WS 2018/2019. Lecture 12. Jürgen Fuhrmann Lecture 12 Slide 1
Scientific Computing WS 2018/2019 Lecture 12 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 12 Slide 1 Recap For more discussion of mesh generation, see J.R. Shewchuk: Lecture Notes on Delaunay
More informationAdaptive numerical methods
METRO MEtallurgical TRaining On-line Adaptive numerical methods Arkadiusz Nagórka CzUT Education and Culture Introduction Common steps of finite element computations consists of preprocessing - definition
More informationINF3320 Computer Graphics and Discrete Geometry
INF3320 Computer Graphics and Discrete Geometry More smooth Curves and Surfaces Christopher Dyken, Michael Floater and Martin Reimers 10.11.2010 Page 1 More smooth Curves and Surfaces Akenine-Möller, Haines
More informationA Hybrid Geometric+Algebraic Multigrid Method with Semi-Iterative Smoothers
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 013; 00:1 18 Published online in Wiley InterScience www.interscience.wiley.com). A Hybrid Geometric+Algebraic Multigrid Method with
More informationAMS527: Numerical Analysis II
AMS527: Numerical Analysis II A Brief Overview of Finite Element Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 1 / 25 Overview Basic concepts Mathematical
More informationVoronoi diagram and Delaunay triangulation
Voronoi diagram and Delaunay triangulation Ioannis Emiris & Vissarion Fisikopoulos Dept. of Informatics & Telecommunications, University of Athens Computational Geometry, spring 2015 Outline 1 Voronoi
More informationEfficient Mesh Optimization Schemes based on Optimal Delaunay Triangulations
Efficient Mesh Optimization Schemes based on Optimal Delaunay Triangulations Long Chen a,, Michael Holst b,1 a Department of Mathematics, University of California at Irvine, Irvine, CA 92697 Email: chenlong@math.uci.edu
More information1 Appendix to notes 2, on Hyperbolic geometry:
1230, notes 3 1 Appendix to notes 2, on Hyperbolic geometry: The axioms of hyperbolic geometry are axioms 1-4 of Euclid, plus an alternative to axiom 5: Axiom 5-h: Given a line l and a point p not on l,
More informationA Data Dependent Triangulation for Vector Fields
A Data Dependent Triangulation for Vector Fields Gerik Scheuermann Hans Hagen Institut for Computer Graphics and CAGD Department of Computer Science University of Kaiserslautern, Postfach 3049, D-67653
More informationDifferential Geometry: Circle Patterns (Part 1) [Discrete Conformal Mappinngs via Circle Patterns. Kharevych, Springborn and Schröder]
Differential Geometry: Circle Patterns (Part 1) [Discrete Conformal Mappinngs via Circle Patterns. Kharevych, Springborn and Schröder] Preliminaries Recall: Given a smooth function f:r R, the function
More informationScientific Computing: Interpolation
Scientific Computing: Interpolation Aleksandar Donev Courant Institute, NYU donev@courant.nyu.edu Course MATH-GA.243 or CSCI-GA.22, Fall 25 October 22nd, 25 A. Donev (Courant Institute) Lecture VIII /22/25
More informationPreferred directions for resolving the non-uniqueness of Delaunay triangulations
Preferred directions for resolving the non-uniqueness of Delaunay triangulations Christopher Dyken and Michael S. Floater Abstract: This note proposes a simple rule to determine a unique triangulation
More informationHigh-Dimensional Computational Geometry. Jingbo Shang University of Illinois at Urbana-Champaign Mar 5, 2018
High-Dimensional Computational Geometry Jingbo Shang University of Illinois at Urbana-Champaign Mar 5, 2018 Outline 3-D vector geometry High-D hyperplane intersections Convex hull & its extension to 3
More information751 Problem Set I JWR. Due Sep 28, 2004
751 Problem Set I JWR Due Sep 28, 2004 Exercise 1. For any space X define an equivalence relation by x y iff here is a path γ : I X with γ(0) = x and γ(1) = y. The equivalence classes are called the path
More informationImage deblurring by multigrid methods. Department of Physics and Mathematics University of Insubria
Image deblurring by multigrid methods Marco Donatelli Stefano Serra-Capizzano Department of Physics and Mathematics University of Insubria Outline 1 Restoration of blurred and noisy images The model problem
More informationIMI Preprint Series INDUSTRIAL MATHEMATICS INSTITUTE 2006:02. Conforming centroidal Voronoi Delaunay triangulation for quality mesh generation. L.
INDUSTRIAL MATHEMATICS INSTITUTE 2006:02 Conforming centroidal Voronoi Delaunay triangulation for quality mesh generation L. Ju IMI Preprint Series Department of Mathematics University of South Carolina
More informationTutorial 3 Comparing Biological Shapes Patrice Koehl and Joel Hass
Tutorial 3 Comparing Biological Shapes Patrice Koehl and Joel Hass University of California, Davis, USA http://www.cs.ucdavis.edu/~koehl/ims2017/ What is a shape? A shape is a 2-manifold with a Riemannian
More informationThe goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a
The goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a coordinate system and then the measuring of the point with
More informationApproximating Polygonal Objects by Deformable Smooth Surfaces
Approximating Polygonal Objects by Deformable Smooth Surfaces Ho-lun Cheng and Tony Tan School of Computing, National University of Singapore hcheng,tantony@comp.nus.edu.sg Abstract. We propose a method
More informationKey words. weak Galerkin, finite element methods, second-order elliptic problems, a posteriori error estimate, polytopal meshes
A POSTERIORI ERROR ESTIMATES FOR THE WEAK GALERKIN FINITE ELEMENT METHODS ON POLYTOPAL MESHES HENGGUANG LI, LIN MU, AND XIU YE Abstract. In this paper, we present a simple a posteriori error estimate for
More informationConvex Hulls in Three Dimensions. Polyhedra
Convex Hulls in Three Dimensions Polyhedra Polyhedron 1.A polyhedron is the generalization of a 2- D polygon to 3-D A finite number of flat polygonal faces The boundary or surface of a polyhedron - Zero-dimensional
More informationControl Volume Finite Difference On Adaptive Meshes
Control Volume Finite Difference On Adaptive Meshes Sanjay Kumar Khattri, Gunnar E. Fladmark, Helge K. Dahle Department of Mathematics, University Bergen, Norway. sanjay@mi.uib.no Summary. In this work
More informationFinite Element Methods
Chapter 5 Finite Element Methods 5.1 Finite Element Spaces Remark 5.1 Mesh cells, faces, edges, vertices. A mesh cell is a compact polyhedron in R d, d {2,3}, whose interior is not empty. The boundary
More informationAn Agglomeration Multigrid Method for Unstructured Grids
Contemporary Mathematics Volume 218, 1998 B 0-8218-0988-1-03002-8 An Agglomeration Multigrid Method for Unstructured Grids Tony F. Chan, Jinchao Xu, and Ludmil Zikatanov 1. Introduction A new agglomeration
More informationQuadratic and cubic b-splines by generalizing higher-order voronoi diagrams
Quadratic and cubic b-splines by generalizing higher-order voronoi diagrams Yuanxin Liu and Jack Snoeyink Joshua Levine April 18, 2007 Computer Science and Engineering, The Ohio State University 1 / 24
More informationPROGRAMMING OF MULTIGRID METHODS
PROGRAMMING OF MULTIGRID METHODS LONG CHEN In this note, we explain the implementation detail of multigrid methods. We will use the approach by space decomposition and subspace correction method; see Chapter:
More informationNESTED AND FULLY AUGMENTED LINKS
NESTED AND FULLY AUGMENTED LINKS HAYLEY OLSON Abstract. This paper focuses on two subclasses of hyperbolic generalized fully augmented links: fully augmented links and nested links. The link complements
More informationPlanar Graphs. 1 Graphs and maps. 1.1 Planarity and duality
Planar Graphs In the first half of this book, we consider mostly planar graphs and their geometric representations, mostly in the plane. We start with a survey of basic results on planar graphs. This chapter
More informationACTUALLY DOING IT : an Introduction to Polyhedral Computation
ACTUALLY DOING IT : an Introduction to Polyhedral Computation Jesús A. De Loera Department of Mathematics Univ. of California, Davis http://www.math.ucdavis.edu/ deloera/ 1 What is a Convex Polytope? 2
More informationBilinear Programming
Bilinear Programming Artyom G. Nahapetyan Center for Applied Optimization Industrial and Systems Engineering Department University of Florida Gainesville, Florida 32611-6595 Email address: artyom@ufl.edu
More informationModeling Discontinuities and their Evolution within Finite Elements: Application to Material Interfaces, 3-D Cracks, and Microstructures
University of California, Davis Modeling Discontinuities and their Evolution within Finite Elements: Application to Material Interfaces, 3-D Cracks, and Microstructures N. Sukumar UC Davis Rutgers University,
More informationIsogeometric Collocation Method
Chair for Computational Analysis of Technical Systems Faculty of Mechanical Engineering, RWTH Aachen University Isogeometric Collocation Method Seminararbeit By Marko Blatzheim Supervisors: Dr. Stefanie
More informationWeek 7 Convex Hulls in 3D
1 Week 7 Convex Hulls in 3D 2 Polyhedra A polyhedron is the natural generalization of a 2D polygon to 3D 3 Closed Polyhedral Surface A closed polyhedral surface is a finite set of interior disjoint polygons
More informationCS 372: Computational Geometry Lecture 10 Linear Programming in Fixed Dimension
CS 372: Computational Geometry Lecture 10 Linear Programming in Fixed Dimension Antoine Vigneron King Abdullah University of Science and Technology November 7, 2012 Antoine Vigneron (KAUST) CS 372 Lecture
More informationSurfaces: notes on Geometry & Topology
Surfaces: notes on Geometry & Topology 1 Surfaces A 2-dimensional region of 3D space A portion of space having length and breadth but no thickness 2 Defining Surfaces Analytically... Parametric surfaces
More informationHierarchical Hybrid Grids
Hierarchical Hybrid Grids IDK Summer School 2012 Björn Gmeiner, Ulrich Rüde July, 2012 Contents Mantle convection Hierarchical Hybrid Grids Smoothers Geometric approximation Performance modeling 2 Mantle
More informationimplicit surfaces, approximate implicitization, B-splines, A- patches, surface fitting
24. KONFERENCE O GEOMETRII A POČÍTAČOVÉ GRAFICE ZBYNĚK ŠÍR FITTING OF PIECEWISE POLYNOMIAL IMPLICIT SURFACES Abstrakt In our contribution we discuss the possibility of an efficient fitting of piecewise
More information9. Three Dimensional Object Representations
9. Three Dimensional Object Representations Methods: Polygon and Quadric surfaces: For simple Euclidean objects Spline surfaces and construction: For curved surfaces Procedural methods: Eg. Fractals, Particle
More informationNonoscillatory Central Schemes on Unstructured Triangular Grids for Hyperbolic Systems of Conservation Laws
Nonoscillatory Central Schemes on Unstructured Triangular Grids for Hyperbolic Systems of Conservation Laws Ivan Christov 1,* Bojan Popov 1 Peter Popov 2 1 Department of Mathematics, 2 Institute for Scientific
More informationLecture 3. Corner Polyhedron, Intersection Cuts, Maximal Lattice-Free Convex Sets. Tepper School of Business Carnegie Mellon University, Pittsburgh
Lecture 3 Corner Polyhedron, Intersection Cuts, Maximal Lattice-Free Convex Sets Gérard Cornuéjols Tepper School of Business Carnegie Mellon University, Pittsburgh January 2016 Mixed Integer Linear Programming
More information2) For any triangle edge not on the boundary, there is exactly one neighboring
Triangulating Trimmed NURBS Surfaces Chang Shu and Pierre Boulanger Abstract. This paper describes techniques for the piecewise linear approximation of trimmed NURBS surfaces. The problem, called surface
More informationPolygonal, Polyhedral, and Serendipity Finite Element Methods
Polygonal, Polyhedral, and Serendipity Finite Element Methods Andrew Gillette Department of Mathematics University of Arizona ASU Computational and Applied Math Seminar Slides and more info at: http://math.arizona.edu/
More informationVoronoi Diagrams and Delaunay Triangulations. O Rourke, Chapter 5
Voronoi Diagrams and Delaunay Triangulations O Rourke, Chapter 5 Outline Preliminaries Properties and Applications Computing the Delaunay Triangulation Preliminaries Given a function f: R 2 R, the tangent
More information7. The Gauss-Bonnet theorem
7. The Gauss-Bonnet theorem 7.1 Hyperbolic polygons In Euclidean geometry, an n-sided polygon is a subset of the Euclidean plane bounded by n straight lines. Thus the edges of a Euclidean polygon are formed
More information1 Exercise: Heat equation in 2-D with FE
1 Exercise: Heat equation in 2-D with FE Reading Hughes (2000, sec. 2.3-2.6 Dabrowski et al. (2008, sec. 1-3, 4.1.1, 4.1.3, 4.2.1 This FE exercise and most of the following ones are based on the MILAMIN
More informationSpline Functions on Triangulations
Spline Functions on Triangulations MING-JUN LAI AND LARRY L. SCHUMAKER CAMBRIDGE UNIVERSITY PRESS Contents Preface xi Chapter 1. Bivariate Polynomials 1.1. Introduction 1 1.2. Norms of Polynomials on Triangles
More informationSimplicial Hyperbolic Surfaces
Simplicial Hyperbolic Surfaces Talk by Ken Bromberg August 21, 2007 1-Lipschitz Surfaces- In this lecture we will discuss geometrically meaningful ways of mapping a surface S into a hyperbolic manifold
More informationSurface Reconstruction. Gianpaolo Palma
Surface Reconstruction Gianpaolo Palma Surface reconstruction Input Point cloud With or without normals Examples: multi-view stereo, union of range scan vertices Range scans Each scan is a triangular mesh
More informationEmbedding Formulations, Complexity and Representability for Unions of Convex Sets
, Complexity and Representability for Unions of Convex Sets Juan Pablo Vielma Massachusetts Institute of Technology CMO-BIRS Workshop: Modern Techniques in Discrete Optimization: Mathematics, Algorithms
More informationCollocation and optimization initialization
Boundary Elements and Other Mesh Reduction Methods XXXVII 55 Collocation and optimization initialization E. J. Kansa 1 & L. Ling 2 1 Convergent Solutions, USA 2 Hong Kong Baptist University, Hong Kong
More informationParallel Adaptive Tsunami Modelling with Triangular Discontinuous Galerkin Schemes
Parallel Adaptive Tsunami Modelling with Triangular Discontinuous Galerkin Schemes Stefan Vater 1 Kaveh Rahnema 2 Jörn Behrens 1 Michael Bader 2 1 Universität Hamburg 2014 PDES Workshop 2 TU München Partial
More informationLecture 0: Reivew of some basic material
Lecture 0: Reivew of some basic material September 12, 2018 1 Background material on the homotopy category We begin with the topological category TOP, whose objects are topological spaces and whose morphisms
More informationThe THex Algorithm and a Simple Darcy Solver on Hexahedral Meshes
Available online at www.sciencedirect.com ScienceDirect Procedia Computer Science 8C (27) 93 92 International Conference on Computational Science, ICCS 27, 2-4 June 27, Zurich, Switzerland The THex Algorithm
More informationThird Order WENO Scheme on Three Dimensional Tetrahedral Meshes
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 5, No. 2-4, pp. 86-848 Commun. Comput. Phys. February 29 Third Order WENO Scheme on Three Dimensional Tetrahedral Meshes Yong-Tao Zhang 1, and Chi-Wang Shu
More informationAn introduction to mesh generation Part IV : elliptic meshing
Elliptic An introduction to mesh generation Part IV : elliptic meshing Department of Civil Engineering, Université catholique de Louvain, Belgium Elliptic Curvilinear Meshes Basic concept A curvilinear
More informationHyperbolic Structures from Ideal Triangulations
Hyperbolic Structures from Ideal Triangulations Craig Hodgson University of Melbourne Geometric structures on 3-manifolds Thurston s idea: We would like to find geometric structures (or metrics) on 3-manifolds
More informationLinear Complexity Hexahedral Mesh Generation
Linear Complexity Hexahedral Mesh Generation David Eppstein Department of Information and Computer Science University of California, Irvine, CA 92717 http://www.ics.uci.edu/ eppstein/ Tech. Report 95-51
More informationFully discrete Finite Element Approximations of Semilinear Parabolic Equations in a Nonconvex Polygon
Fully discrete Finite Element Approximations of Semilinear Parabolic Equations in a Nonconvex Polygon Tamal Pramanick 1,a) 1 Department of Mathematics, Indian Institute of Technology Guwahati, Guwahati
More informationJune 5, Institute of Structural Analysis Graz University of Technology Lessingstr. 25/II, 8010 Graz, Austria
Higher-order meshing of implicit geometries part I: Integration and interpolation in cut elements arxiv:706.00578v [cs.na] 2 Jun 207 T.P. Fries, S. Omerović, D. Schöllhammer, J. Steidl June 5, 207 Institute
More informationVoronoi Diagram. Xiao-Ming Fu
Voronoi Diagram Xiao-Ming Fu Outlines Introduction Post Office Problem Voronoi Diagram Duality: Delaunay triangulation Centroidal Voronoi tessellations (CVT) Definition Applications Algorithms Outlines
More informationThis is an author-deposited version published in: Eprints ID: 4362
This is an author-deposited version published in: http://oatao.univ-toulouse.fr/ Eprints ID: 4362 To cite this document: CHIKHAOUI Oussama, GRESSIER Jérémie, GRONDIN Gilles. Assessment of the Spectral
More informationChapter 4 Concepts from Geometry
Chapter 4 Concepts from Geometry An Introduction to Optimization Spring, 2014 Wei-Ta Chu 1 Line Segments The line segment between two points and in R n is the set of points on the straight line joining
More informationarxiv:cs/ v1 [cs.cg] 26 Feb 2002
The efficient generation of unstructured control volumes in 2D and 3D Jacek Leszczynski, Sebastian Pluta arxiv:cs/0202038v1 [cs.cg] 26 Feb 2002 Technical University of Czestochowa, Institute of Mathematics
More informationParameterization of triangular meshes
Parameterization of triangular meshes Michael S. Floater November 10, 2009 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to
More information