Benchmark 3D: the VAG scheme
|
|
- Derick Sutton
- 5 years ago
- Views:
Transcription
1 Benchmark D: the VAG scheme Robert Eymard, Cindy Guichard and Raphaèle Herbin Presentation of the scheme Let Ω be a bounded open domain of R, let f L 2 Ω) and let Λ be a measurable function from Ω to the set M R) of matrices, such that for a.e. x Ω, Λx) is symmetric, and such that the set of its eigenvalues is included in [λ,λ], where 0 < λ λ. We wish to approximate the function u solution of u H0 Ω) and v H0 Ω), Λx) ux) vx)dx = f x)vx)dx, ) Ω Ω by the approximate gradient scheme [4, 2] which reads: U X D, V X D, Λx) D Ux) D V x)dx = f x)π D V x)dx, 2) Ω Ω where Π D is a reconstruction operator and D a discrete gradient operator which act on the discrete functional space X D, where the index D denotes the discretization; these operators are defined as defined as follows:. M is the set of control volumes, that are disjoint open subsets of Ω such that K = Ω, V = V int V ext is the set of vertices of the mesh; any element K of M is defined by its vertices s V K, its faces σ F K ; each face is also defined by the set of its vertices s V σ, using a suitable geometric definition for R. Eymard Université Paris-Est, robert.eymard@univ-mlv.fr C. Guichard IFP Energies nouvelles and Université Paris-Est, guichard@ifpenergiesnouvelles.fr R. Herbin Université Aix-Marseille, Raphaele.Herbin@latp.univ-mrs.fr Work supported by Groupement MOMAS and ANR VFSitCom
2 2 R. Eymard, C. Guichard and R. Herbin the resulting surface in the case of non-planar faces; we assume that Λ is constant on all K M, and we denote by Λ K its value in K; 2. the set of discrete unknowns X D is the finite dimensional vector space on R, containing all real families U = u K ),u s ) s V ), such that u s = 0 if s V ext ;. the mapping Π D : X D L 2 Ω) maps U = u K ),u s ) s V ) X D to the piecewise constant function u D L 2 Ω) equal to u K on each cell K M ; 4. the mapping D : X D L 2 Ω) is the reconstruction of a gradient from the values U = u K ),u s ) s V ) X D ; different expressions for this reconstruction are proposed below, which all lead to convergent gradient schemes in the sense of [4, 2]. Their theoretical analysis is related to that of the SUSHI scheme [, ]. Detailed numerical results are given in this paper only using the method described in subsection.2; the differences obtained using the other expressions are commented in the last section. The exterior faces are those of the form K Ω for any boundary control volume K, and the interior faces are those of the form K L for two neighbouring control volumes K and L. For any face σ, we define a point x σ, which is a barycentre with non-negative weights β σ,s of the elements of the set V σ including all the vertices of the face, and the value u σ is defined by u σ = s V σ β σ,s u s with s V σ β σ,s =. In the next three subsections, we describe three ways of defining a gradient operator which satisfies the VAG requirements. The first gradient is constructed from the Stokes formula on the cells of the mesh we call it the primal cell to distinguish it from further constructed cells), and requires a stabilization. The second gradient and third gradients are constructed on tetrahedral or octahedral sub cells of the primal mesh, and are natively stable.. Stabilised gradient on the primal mesh cells For a face σ F, we denote by τ any triangular sub-face with vertices x σ, s and s, where s and s are two consecutive vertices of σ. The barycentre x τ of each sub-face τ may thus be expressed by the following barycentric combination: x τ = s V σ β τ,s s with s V σ β τ,s =, where β τ,s 0 for all s V σ. We then define β τ,s = 0 for all s V \ V σ. Next, we reconstruct a value u τ at the point x τ, by u τ = s Vσ β τ,s u s. Let K M be an element of the mesh. We denote by T K the set of all sub-faces of the faces of K. We first define, for U = u K ),u s ) s V ), an approximation of the gradient on cell K:
3 Benchmark D: the VAG scheme K U = K where we denote by τ u τ u K )n K,τ = u s u K )v K,s, ) τ T K s V K v K,s = K β τ,s τ n K,τ, τ T K where n K,τ is the unit normal vector to τ, outward to K, and τ, K are respectively the area and the volume of τ and K. We then define a partition M K,s of K there is no need to define this partition precisely), such that M K,s = K /N K, where N K is the number of vertices of K and we introduce R K,s U = u s u K K U s x K ). We then define, for a given γ > 0, the constant value K,s U in M K,s : K,s U = K u + γr K,s Uv K,s. We finally define a piecewise constant gradient by D Ux) = K,s U for a.e. x M K,s. This scheme is denoted by VAG in [5]..2 Piecewise constant gradient on octahedral sub-cells For a given face σ of a control volume K and for any vertex s of σ, we respectively denote by s and s + the preceding and the following vertices of s in the face σ defining any orientation on σ), and we consider the degenerate) octahedron, denoted by V K,σ,s and depicted in Fig., whose vertices are A = x K, A 2 = x σ, A = 2 s +s), A 5 = s, A 6 = 2 s+ +s) and A 4 = 2 x σ +s) note that all these octahedra are disjoint, and that the union of their closure is Ω). The approximate values of U at the vertices of V K,σ,s are respectively u = u K, u 2 = u σ, u = 2 u s + u s), u 5 = u s, u 6 = 2 u s + + u s) and u 4 = 2 u σ + u s ) the main diagonals of V K,σ,s are therefore A,A 4 ), A 2,A 5 ) and A,A 6. We then define the following approximate gradient: K,σ,s U = i= u i+ u i ) A i+ A i+4 A i+2 A i+5 Det A i+ A i+4, A i+2 A i+5, A i A i+ ), 4) setting A 7 = A and A 8 = A 2. We finally define a piecewise constant gradient by D Ux) = K,σ,s U for a.e. x V K,σ,s. Remark that, denoting for simplicity V instead of V K,σ,s, defining F V as the set of the 8 triangular faces of V and V τ as the set of the vertices of each triangular face τ of V, one may check that K,σ,s U = V τ n V,τ τ F ) u s. 5) V s V τ
4 4 R. Eymard, C. Guichard and R. Herbin A 5 = s A = s A 6 A 4 s + A = x K A A 0 = x K A 2 = x σ A = x σ s A 2 = s Fig. The octahedral left) and tetrahedral right) cells for the definition of the gradient We then set define K U, used in the tables below, by K K U = σ F K This scheme is denoted by VAGR in [5]. s V σ V K,σ,s K,σ,s U.. Piecewise constant gradient on tetrahedral sub cells For a given face σ of a control volume K and for any pair of consecutive vertices s,s ) of σ, we consider the tetrahedron, denoted by V K,σ,s,s and depicted in Fig., whose vertices are A 0 = x K, A = x σ, A 2 = s and A = s note that all these tetrahedra are disjoint, and that the union of their closure is Ω). The approximate values of U at the vertices of V K,σ,s,s are respectively u 0 = u K, u = u σ, u 2 = u s and u = u s. We then define the following approximate gradient: K,σ,s,s U = i= u i u 0 ) A 0 A i+ A 0 A i+2 Det A 0 A i+, A 0 A i+2, A 0 A i ), 6) where A 4 = A and A 5 = A 2. We finally define a piecewise constant gradient by D Ux) = K,σ,s,s U for a.e. x V K,σ,s,s. Remark that 5) also holds in this case, denoting V instead of V K,σ,s,s.
5 Benchmark D: the VAG scheme 5 2 Numerical results We provide the detailed numerical results obtained, using the scheme VAGR for computing the discrete gradient. In the numerical implementation, the values u K are locally eliminated, and the unknowns of the linear solver are the values u s. Denoting by denotes the Euclidean norm, the norms used in the bench tables have been computed using the following formulae: erl2 = normg = K K U, K u K ux K )/ 2 /2 K ux K ) 2, /2 ergrad = K K U ux K ) )/ 2 K ux K ) 2, ) /2 ener = K K U ux K ) Λ 2 / K ux K ) Λ 2, setting, for any K M, ξ 2 Λ = Λ Kξ ξ for all ξ R. Test Mild anisotropy, ux,y,z) = + sinπx)sin π y + 2 min = 0, max = 2, Tetrahedral meshes sin π z E E-02.95E+00.99E+00.77E E E-0.97E+00.99E+00.78E E-02 9.E-0.98E+00.99E+00.79E E E-0.99E E+00.79E E-0.49E-0.99E E+00.79E E-0.8E-0.99E E+00.80E E E E E E+00.79E-0.5E+00.77E-0.5E E-0.94E+00.44E-0.05E+00.42E-0.08E E-0 2.0E+00.E-0.4E+00.E-0.7E E-0 2.8E E-02.06E E-02.0E E-0 2.E E-02.E E-02.5E+00
6 6 R. Eymard, C. Guichard and R. Herbin Test Mild anisotropy, ux,y,z) = + sinπx)sin π y + 2 sin π z + min = 0, max = 2, Voronoi meshes E-02.56E-0 2.5E+00.86E+00.4E E-0.79E-0.95E+00.8E+00.4E E E E+00.9E+00.60E E-02.20E E+00.9E+00.66E E-02.85E E+00.97E+00.70E E E E E-0-8.4E E-0.65E E-0.68E E E+00.55E E+00.62E-0.90E E-02.68E+00.9E-0.5E+00.2E-0.42E E E E-02.E E-02.29E+00 Test Mild anisotropy, ux,y,z) = + sinπx)sin π y + 2 min = 0, max = 2, Kershaw meshes sin π z E-02.0E-02.96E+00.96E+00.56E E-02.06E-02.98E+00.99E+00.68E E-04.75E-0.99E E+00.74E E E E E+00.78E E E E E E-0.09E E E E E E-0.74E E-0.54E E E-02.E E-02.26E E-02.29E+00 Test Mild anisotropy, ux,y,z) = + sinπx)sin π y + 2 min = 0, max = 2, Checkerboard meshes sin π z E-02.54E E+00.85E+00.4E E-0 4.0E E+00.96E+00.70E E-02.0E E+00.99E+00.78E E E E E+00.79E E-0 6.6E E E+00.80E+00
7 Benchmark D: the VAG scheme E-0-4.7E E E-0.7E+00.50E-0.72E+00.52E-0.54E E E E-02.47E E-02.4E E-0.99E E-02.2E E-02.8E E E+00.8E-02.E+00.2E-02.08E+00 Test 2 Heterogeneous anisotropy, ux,y,z) = x y 2 z + xsin2πxz)sin2πxy)sin2πz), min = 0.862, max =.048, Prism meshes E-0-8.4E-0.0E E-0.5E E-0-8.9E-0.04E+00.0E+00.66E E E-0.05E+00.0E+00.69E E E-0.04E+00.0E+00.70E E E-0 -.8E E E+00.7E-02 2.E+00.64E-02 2.E E E+00.67E E+00.6E E E E E E E E+00 Test Flow on random meshes, ux, y, z) = sin2πx) sin2πy) sin2πz), min =, max =, Random meshes E E-0.68E E-0.5E E E-0.2E E E E E-0.06E E-0.44E E E-0.0E E-0.56E E E E E E E-0.90E+00.8E-0 2.E E E E-02.45E E-02.9E E E+00.45E-02.28E+00.74E-02.68E+00
8 8 R. Eymard, C. Guichard and R. Herbin Test 4 Flow around a well, Well meshes E E-0 5.2E E+00.68E E E-0 5.E+00 5.E+00.65E E-0.57E-0 5.E+00 5.E+00.64E E-0.20E-0 5.E+00 5.E+00.6E E E E E+00.6E E E E E+00.6E E E E E+00.6E E E E E-0.2E E-02.05E E-02.08E E-0.E+00.2E E+00.5E E E-04.7E E-0.9E+00 7.E-0.79E E E E-0.7E E-0.68E E E+00.55E-0.8E+00.47E-0.80E E E-0.26E E-0.9E E-0 Test 5 Discontinuous permeability, ux,y,z) = α i sin2πx)sin2πy)sin2πz), min = 0, max =, Locally refined meshes E E E+02.00E+02.24E E E E+0.54E E E E E E E E+0-9.4E E+0 9.4E E E E+0.00E E E E E E E+00.0E E-0.82E E E E E+00.86E-0.92E+00.80E-0.99E E E E-02.76E E E E E E-02.46E+00.26E E+00 Comments on the results The results obtained using ) VAG) instead of 4) VAGR) are systematically less precise, except in the test5 case, where we obtained the following tables:
9 Benchmark D: the VAG scheme E E E+02.00E E E E+0 7.7E+0.54E E E E E E+0 8.9E E+0-9.4E E+0 9.4E+0 9.4E E E+0 9.9E E E E E E E-0.60E E-0 4.6E E-0 6.9E E-0.05E E-0.85E E-0.85E E E E E E E E-0 2.E+00.7E-02.90E+00.54E E+00 The results using 6) are very similar to those obtained using 4) VAGR). For both ) VAG) and 4) VAGR), we have chosen the conjugate gradient solver of the ISTL library with ILU0) preconditioning with tolerance or reduction factor) set to 0 0. The following observations have been made on the computing times, using ) VAG) we may expect that similar observations could be done with VAGR).. On the fourth Kershaw mesh and test, we obtain the following CPU times using the conjugate gradient solver of the PETSC library: with ILU2), s, with ILU), 7s, with ILU0), 0s, and with Jacobi, s, which shows that the ILU0) preconditioning seems the fastest one on this case. Note that this computing time is depending on the unknown orderings. For the bench computations, we used the recursive domain decomposition ordering, which is the most efficient for direct solvers, and the respective computing times with PETSC CG+ILU0) and with ISTL CG+ILU0) are 0. and.2 s. Using the reverse Cuthill - McKee ordering, we respectively obtain 4.4 s and 5. s with PETSC CG+ILU0) and ISTL CG+ILU0). 2. The computing times, for the conjugate gradient solver of the PETSC library with ILU) preconditioning, in the test case on tetrahedral meshes 2 to 5, have been approximately equal to 0.0, 0.0, 0.04, 0.08, and 0.6 s, showing the possibility to apply this method on much larger meshes. Finally, we may not exclude that the systematic choice of computing the L 2 error with respect to the values in the control volumes instead of the vertex values, makes all these results somewhat pessimistic. References. R. Eymard, T. Gallouët, and R. Herbin. Discretisation of heterogeneous and anisotropic diffusion problems on general non-conforming meshes, sushi: a scheme using stabilisation and hybrid interfaces. IMA J. Numer. Anal., 04):009 04, 200. see also
10 0 R. Eymard, C. Guichard and R. Herbin 2. R. Eymard, C. Guichard, and R. Herbin. Small-stencil D schemes for diffusive flows in porous media. submitted, 200. see also R. Eymard, T. Gallouët and R. Herbin. Benchmark D: the SUSHI scheme these proceedings, R. Eymard and R. Herbin. Gradient schemes for diffusion problem. these proceedings, R. Eymard, G. Henry, R. Herbin, F. Hubert, R. Klöfkorn and G. Manzini. D Benchmark on Discretization Schemes for Anisotropic Diffusion Problems on General Grids. these proceedings, 20. The paper is in final form and no similar paper has been or is being submitted elsewhere.
3D Benchmark on Discretization Schemes for Anisotropic Diffusion Problems on General Grids
3D Benchmark on Discretization Schemes for Anisotropic Diffusion Problems on General Grids Robert Eymard, Gérard Henry, Rapahele Herbin, Florence Hubert, Robert Klofkorn, Gianmarco Manzini To cite this
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 informationStability Analysis of the Muscl Method on General Unstructured Grids for Applications to Compressible Fluid Flow
Stability Analysis of the Muscl Method on General Unstructured Grids for Applications to Compressible Fluid Flow F. Haider 1, B. Courbet 1, J.P. Croisille 2 1 Département de Simulation Numérique des Ecoulements
More informationWeno Scheme for Transport Equation on Unstructured Grids with a DDFV Approach
Weno Scheme for Transport Equation on Unstructured Grids with a DDFV Approach Florence Hubert and Rémi Tesson Abstract In this paper we develop a DDFV approach for WENO scheme on unstructred grids for
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 informationLecture 15: More Iterative Ideas
Lecture 15: More Iterative Ideas David Bindel 15 Mar 2010 Logistics HW 2 due! Some notes on HW 2. Where we are / where we re going More iterative ideas. Intro to HW 3. More HW 2 notes See solution code!
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 informationarxiv: v1 [math.na] 7 Jul 2014
arxiv:1407.1567v1 [math.na] 7 Jul 2014 FINITE VOLUME SCHEMES FOR DIFFUSION EQUATIONS: INTRODUCTION TO AND REVIEW OF MODERN METHODS JEROME DRONIOU School of Mathematical Sciences, Monash University Victoria
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 informationFrom Hyperbolic Diffusion Scheme to Gradient Method: Implicit Green-Gauss Gradients for Unstructured Grids
Preprint accepted in Journal of Computational Physics. https://doi.org/10.1016/j.jcp.2018.06.019 From Hyperbolic Diffusion Scheme to Gradient Method: Implicit Green-Gauss Gradients for Unstructured Grids
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 informationGraph Partitioning for High-Performance Scientific Simulations. Advanced Topics Spring 2008 Prof. Robert van Engelen
Graph Partitioning for High-Performance Scientific Simulations Advanced Topics Spring 2008 Prof. Robert van Engelen Overview Challenges for irregular meshes Modeling mesh-based computations as graphs Static
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 informationALF USER GUIDE. Date: September
ALF USER GUIDE R. VERFÜRTH Contents 1. Introduction 1 2. User interface 2 3. Domain and initial grid 3 4. Differential equation 9 5. Discretization 12 6. Solver 12 7. Mesh refinement 14 8. Number of nodes
More informationParameterization. Michael S. Floater. November 10, 2011
Parameterization Michael S. Floater November 10, 2011 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to generate from point
More informationA face-centred finite volume method for second-order elliptic problems
A face-centred finite volume method for second-order elliptic problems arxiv:1712.6173v1 [math.na] 17 Dec 217 Ruben Sevilla Zienkiewicz Centre for Computational Engineering, College of Engineering, Swansea
More informationpα i + q, where (n, m, p and q depend on i). 6. GROMOV S INVARIANT AND THE VOLUME OF A HYPERBOLIC MANIFOLD
6. GROMOV S INVARIANT AND THE VOLUME OF A HYPERBOLIC MANIFOLD of π 1 (M 2 )onπ 1 (M 4 ) by conjugation. π 1 (M 4 ) has a trivial center, so in other words the action of π 1 (M 4 ) on itself is effective.
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 informationA NEW LEVEL SET METHOD FOR MOTION IN NORMAL DIRECTION BASED ON A FORWARD-BACKWARD DIFFUSION FORMULATION
A NEW LEVEL SET METHOD FOR MOTION IN NORMAL DIRECTION BASED ON A FORWARD-BACKWARD DIFFUSION FORMULATION KAROL MIKULA AND MARIO OHLBERGER Abstract. We introduce a new level set method for motion in normal
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 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 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 informationKeywords: Block ILU preconditioner, Krylov subspace methods, Additive Schwarz, Domain decomposition
BLOCK ILU PRECONDITIONERS FOR PARALLEL AMR/C SIMULATIONS Jose J. Camata Alvaro L. G. A. Coutinho Federal University of Rio de Janeiro, NACAD, COPPE Department of Civil Engineering, Rio de Janeiro, Brazil
More informationDual Interpolants for Finite Element Methods
Dual Interpolants for Finite Element Methods Andrew Gillette joint work with Chandrajit Bajaj and Alexander Rand Department of Mathematics Institute of Computational Engineering and Sciences University
More informationOn the Construction, Comparison, and Local Characteristic Decomposition for High-Order Central WENO Schemes
Journal of Computational Physics 8, 87 09 (00) doi:0.006/jcph.00.79 On the Construction, Comparison, and Local Characteristic Decomposition for High-Order Central WENO Schemes Jianxian Qiu, and Chi-Wang
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 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 informationSELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND
Student Submission for the 5 th OpenFOAM User Conference 2017, Wiesbaden - Germany: SELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND TESSA UROIĆ Faculty of Mechanical Engineering and Naval Architecture, Ivana
More informationArtificial Intelligence for Robotics: A Brief Summary
Artificial Intelligence for Robotics: A Brief Summary This document provides a summary of the course, Artificial Intelligence for Robotics, and highlights main concepts. Lesson 1: Localization (using Histogram
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 information3D NURBS-ENHANCED FINITE ELEMENT METHOD
7th Workshop on Numerical Methods in Applied Science and Engineering (NMASE 8) Vall de Núria, 9 a 11 de enero de 28 c LaCàN, www.lacan-upc.es 3D NURBS-ENHANCED FINITE ELEMENT METHOD R. Sevilla, S. Fernández-Méndez
More informationHigh Order Nédélec Elements with local complete sequence properties
High Order Nédélec Elements with local complete sequence properties Joachim Schöberl and Sabine Zaglmayr Institute for Computational Mathematics, Johannes Kepler University Linz, Austria E-mail: {js,sz}@jku.at
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 informationINVESTIGATING the numerical schemes with high accuracy
March 18-20 2015 Hong ong Monotone Finite Volume Schemes for Diffusion Equations on Distorted Quadrilateral Meshes Jing-yan Yue Guang-wei Yuan and Zhi-qiang Sheng Abstract We construct a new monotone finite
More informationThe 3D DSC in Fluid Simulation
The 3D DSC in Fluid Simulation Marek K. Misztal Informatics and Mathematical Modelling, Technical University of Denmark mkm@imm.dtu.dk DSC 2011 Workshop Kgs. Lyngby, 26th August 2011 Governing Equations
More informationSome Computational Results for Dual-Primal FETI Methods for Elliptic Problems in 3D
Some Computational Results for Dual-Primal FETI Methods for Elliptic Problems in 3D Axel Klawonn 1, Oliver Rheinbach 1, and Olof B. Widlund 2 1 Universität Duisburg-Essen, Campus Essen, Fachbereich Mathematik
More informationPARALLEL MULTILEVEL TETRAHEDRAL GRID REFINEMENT
PARALLEL MULTILEVEL TETRAHEDRAL GRID REFINEMENT SVEN GROSS AND ARNOLD REUSKEN Abstract. In this paper we analyze a parallel version of a multilevel red/green local refinement algorithm for tetrahedral
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 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 informationTowards a complete FEM-based simulation toolkit on GPUs: Geometric Multigrid solvers
Towards a complete FEM-based simulation toolkit on GPUs: Geometric Multigrid solvers Markus Geveler, Dirk Ribbrock, Dominik Göddeke, Peter Zajac, Stefan Turek Institut für Angewandte Mathematik TU Dortmund,
More informationLecture 5: Simplicial Complex
Lecture 5: Simplicial Complex 2-Manifolds, Simplex and Simplicial Complex Scribed by: Lei Wang First part of this lecture finishes 2-Manifolds. Rest part of this lecture talks about simplicial complex.
More informationScalable, Hybrid-Parallel Multiscale Methods using DUNE
MÜNSTER Scalable Hybrid-Parallel Multiscale Methods using DUNE R. Milk S. Kaulmann M. Ohlberger December 1st 2014 Outline MÜNSTER Scalable Hybrid-Parallel Multiscale Methods using DUNE 2 /28 Abstraction
More informationA GEOMETRIC DESCRIPTION OF THE PSL 4 (R)-HITCHIN COMPONENT
A GEOMETRIC DESCRIPTION OF THE PSL 4 (R)-HITCHIN COMPONENT SAMUEL A. BALLAS Abstract. These notes form a rough outline of the correspondence between the PSL 4 (R)- Hitchin component and convex foliated
More informationOpen and Closed Sets
Open and Closed Sets Definition: A subset S of a metric space (X, d) is open if it contains an open ball about each of its points i.e., if x S : ɛ > 0 : B(x, ɛ) S. (1) Theorem: (O1) and X are open sets.
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 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 informationBasics of Combinatorial Topology
Chapter 7 Basics of Combinatorial Topology 7.1 Simplicial and Polyhedral Complexes In order to study and manipulate complex shapes it is convenient to discretize these shapes and to view them as the union
More informationA Hash Data Structure for Adaptive PDE Solvers Based on Discontinuous Galerkin Discretizations
A Hash Data Structure for Adaptive PDE Solvers Based on Discontinuous Galerkin Discretizations Kolja Brix, Ralf Massjung, and Alexander Voß Bericht Nr. 302 Juli 2009 Key words: Discontinuous Galerkin method,
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 informationHPC Algorithms and Applications
HPC Algorithms and Applications Dwarf #5 Structured Grids Michael Bader Winter 2012/2013 Dwarf #5 Structured Grids, Winter 2012/2013 1 Dwarf #5 Structured Grids 1. dense linear algebra 2. sparse linear
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 informationREINER HORST. License or copyright restrictions may apply to redistribution; see
MATHEMATICS OF COMPUTATION Volume 66, Number 218, April 1997, Pages 691 698 S 0025-5718(97)00809-0 ON GENERALIZED BISECTION OF n SIMPLICES REINER HORST Abstract. A generalized procedure of bisection of
More informationSpace Filling Curves and Hierarchical Basis. Klaus Speer
Space Filling Curves and Hierarchical Basis Klaus Speer Abstract Real world phenomena can be best described using differential equations. After linearisation we have to deal with huge linear systems of
More informationA Voronoi-Based Hybrid Meshing Method
A Voronoi-Based Hybrid Meshing Method Jeanne Pellerin, Lévy Bruno, Guillaume Caumon To cite this version: Jeanne Pellerin, Lévy Bruno, Guillaume Caumon. A Voronoi-Based Hybrid Meshing Method. 2012. hal-00770939
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 informationSurface Parameterization
Surface Parameterization A Tutorial and Survey Michael Floater and Kai Hormann Presented by Afra Zomorodian CS 468 10/19/5 1 Problem 1-1 mapping from domain to surface Original application: Texture mapping
More informationAn Interface-fitted Mesh Generator and Polytopal Element Methods for Elliptic Interface Problems
An Interface-fitted Mesh Generator and Polytopal Element Methods for Elliptic Interface Problems Long Chen University of California, Irvine chenlong@math.uci.edu Joint work with: Huayi Wei (Xiangtan University),
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 informationAn easy treatment of hanging nodes in hp-finite elements
An easy treatment of hanging nodes in hp-finite elements Paolo Di Stolfo a,, Andreas Schröder a, Nils Zander b, Stefan Kollmannsberger b a Department of Mathematics, Universität Salzburg, Hellbrunner Str.
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 informationC 1 Quintic Spline Interpolation Over Tetrahedral Partitions
C 1 Quintic Spline Interpolation Over Tetrahedral Partitions Gerard Awanou and Ming-Jun Lai Abstract. We discuss the implementation of a C 1 quintic superspline method for interpolating scattered data
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 informationHigh Order Weighted Essentially Non-Oscillatory Schemes for Convection. Dominated Problems. Chi-Wang Shu 1
High Order Weighted Essentially Non-Oscillatory Schemes for Convection Dominated Problems Chi-Wang Shu Division of Applied Mathematics, Brown University, Providence, Rhode Island 09 ABSTRACT High order
More informationMath 462: Review questions
Math 462: Review questions Paul Hacking 4/22/10 (1) What is the angle between two interior diagonals of a cube joining opposite vertices? [Hint: It is probably quickest to use a description of the cube
More informationObtaining the H and T Honeycomb from a Cross-Section of the 16-cell Honeycomb
Bridges 2017 Conference Proceedings Obtaining the H and T Honeycomb from a Cross-Section of the 16-cell Honeycomb Hideki Tsuiki Graduate School of Human and Environmental Studies, Kyoto University Yoshida-Nihonmatsu,
More informationScalable Algorithms in Optimization: Computational Experiments
Scalable Algorithms in Optimization: Computational Experiments Steven J. Benson, Lois McInnes, Jorge J. Moré, and Jason Sarich Mathematics and Computer Science Division, Argonne National Laboratory, Argonne,
More informationMATHEMATICAL ANALYSIS, MODELING AND OPTIMIZATION OF COMPLEX HEAT TRANSFER PROCESSES
MATHEMATICAL ANALYSIS, MODELING AND OPTIMIZATION OF COMPLEX HEAT TRANSFER PROCESSES Goals of research Dr. Uldis Raitums, Dr. Kārlis Birģelis To develop and investigate mathematical properties of algorithms
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 informationLinear Finite Element Methods
Chapter 3 Linear Finite Element Methods The finite element methods provide spaces V n of functions that are piecewise smooth and simple, and locally supported basis function of these spaces to achieve
More informationChapter 13. Boundary Value Problems for Partial Differential Equations* Linz 2002/ page
Chapter 13 Boundary Value Problems for Partial Differential Equations* E lliptic equations constitute the third category of partial differential equations. As a prototype, we take the Poisson equation
More informationNumerical Experiments
77 Chapter 4 Numerical Experiments 4.1 Error estimators and adaptive refinement Due to singularities the convergence of finite element solutions on uniform grids can be arbitrarily low. Adaptivity based
More informationKatholieke Universiteit Leuven Department of Computer Science
Numerical solution of partial differential equations with Powell-Sabin splines Hendrik Speleers Paul Dierckx Stefan Vandewalle Report TW 48, October 24 Katholieke Universiteit Leuven Department of Computer
More informationParameterization II Some slides from the Mesh Parameterization Course from Siggraph Asia
Parameterization II Some slides from the Mesh Parameterization Course from Siggraph Asia 2008 1 Non-Convex Non Convex Boundary Convex boundary creates significant distortion Free boundary is better 2 Fixed
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 informationPersistent Homology and Nested Dissection
Persistent Homology and Nested Dissection Don Sheehy University of Connecticut joint work with Michael Kerber and Primoz Skraba A Topological Data Analysis Pipeline A Topological Data Analysis Pipeline
More informationComputing NodeTrix Representations of Clustered Graphs
Journal of Graph Algorithms and Applications http://jgaa.info/ vol. 22, no. 2, pp. 139 176 (2018) DOI: 10.7155/jgaa.00461 Computing NodeTrix Representations of Clustered Graphs Giordano Da Lozzo Giuseppe
More informationCAT(0)-spaces. Münster, June 22, 2004
CAT(0)-spaces Münster, June 22, 2004 CAT(0)-space is a term invented by Gromov. Also, called Hadamard space. Roughly, a space which is nonpositively curved and simply connected. C = Comparison or Cartan
More informationOn the graphical display of Powell-Sabin splines: a comparison of three piecewise linear approximations
On the graphical display of Powell-Sabin splines: a comparison of three piecewise linear approximations Hendrik Speleers Paul Dierckx Stefan Vandewalle Report TW515, January 008 Ò Katholieke Universiteit
More informationReducing the memory footprint of an Eikonal solver
SpezialForschungsBereich F 32 Karl Franzens Universita t Graz Technische Universita t Graz Medizinische Universita t Graz Reducing the memory footprint of an Eikonal solver D. Ganellari G. Haase SFB-Report
More informationCOMPARING SEMI-NORMS ON HOMOLOGY
COMPARING SEMI-NORMS ON HOMOLOGY J.-F. LAFONT AND CH. PITTET Abstract. We compare the l 1 -seminorm 1 and the manifold seminorm man on n-dimensional integral homology classes. We explain how it easily
More informationClustering algorithms and introduction to persistent homology
Foundations of Geometric Methods in Data Analysis 2017-18 Clustering algorithms and introduction to persistent homology Frédéric Chazal INRIA Saclay - Ile-de-France frederic.chazal@inria.fr Introduction
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 informationMorse Theory. Investigates the topology of a surface by looking at critical points of a function on that surface.
Morse-SmaleComplex Morse Theory Investigates the topology of a surface by looking at critical points of a function on that surface. = () () =0 A function is a Morse function if is smooth All critical points
More informationTopology Changes for Surface Meshes in Active Contours
Topology Changes for Surface Meshes in Active Contours Jochen Abhau University of Innsbruck, Department of Computer Science FSP Conference Vorau, November 2006 Outline 1 Introduction 2 Topological Adaptation
More informationHigh Performance Computing: Tools and Applications
High Performance Computing: Tools and Applications Edmond Chow School of Computational Science and Engineering Georgia Institute of Technology Lecture 15 Numerically solve a 2D boundary value problem Example:
More informationγ 2 γ 3 γ 1 R 2 (b) a bounded Yin set (a) an unbounded Yin set
γ 1 γ 3 γ γ 3 γ γ 1 R (a) an unbounded Yin set (b) a bounded Yin set Fig..1: Jordan curve representation of a connected Yin set M R. A shaded region represents M and the dashed curves its boundary M that
More informationIMAGE ANALYSIS DEDICATED TO POLYMER INJECTION MOLDING
Image Anal Stereol 2001;20:143-148 Original Research Paper IMAGE ANALYSIS DEDICATED TO POLYMER INJECTION MOLDING DAVID GARCIA 1, GUY COURBEBAISSE 2 AND MICHEL JOURLIN 3 1 European Polymer Institute (PEP),
More informationGeometric Modeling Mortenson Chapter 11. Complex Model Construction
Geometric Modeling 91.580.201 Mortenson Chapter 11 Complex Model Construction Topics Topology of Models Connectivity and other intrinsic properties Graph-Based Models Emphasize topological structure Boolean
More informationStabilized Finite Element Method for 3D Navier-Stokes Equations with Physical Boundary Conditions
Stabilized Finite Element Method for 3D Navier-Stokes Equations with Physical Boundary Conditions Mohamed Amara, Daniela Capatina, David Trujillo Université de Pau et des Pays de l Adour Laboratoire de
More informationGRAPH CENTERS USED FOR STABILIZATION OF MATRIX FACTORIZATIONS
Discussiones Mathematicae Graph Theory 30 (2010 ) 349 359 GRAPH CENTERS USED FOR STABILIZATION OF MATRIX FACTORIZATIONS Pavla Kabelíková Department of Applied Mathematics FEI, VSB Technical University
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 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 information1 2 (3 + x 3) x 2 = 1 3 (3 + x 1 2x 3 ) 1. 3 ( 1 x 2) (3 + x(0) 3 ) = 1 2 (3 + 0) = 3. 2 (3 + x(0) 1 2x (0) ( ) = 1 ( 1 x(0) 2 ) = 1 3 ) = 1 3
6 Iterative Solvers Lab Objective: Many real-world problems of the form Ax = b have tens of thousands of parameters Solving such systems with Gaussian elimination or matrix factorizations could require
More informationAPPROXIMATING PDE s IN L 1
APPROXIMATING PDE s IN L 1 Veselin Dobrev Jean-Luc Guermond Bojan Popov Department of Mathematics Texas A&M University NONLINEAR APPROXIMATION TECHNIQUES USING L 1 Texas A&M May 16-18, 2008 Outline 1 Outline
More informationImplicit schemes for wave models
Implicit schemes for wave models Mathieu Dutour Sikirić Rudjer Bo sković Institute, Croatia and Universität Rostock April 17, 2013 I. Wave models Stochastic wave modelling Oceanic models are using grids
More informationExact Volume of Hyperbolic 2-bridge Links
Exact Volume of Hyperbolic 2-bridge Links Anastasiia Tsvietkova Louisiana State University Parts of this work are joint with M. Thislethwaite, O. Dasbach Anastasiia Tsvietkova (Louisiana State University
More informationElement Quality Metrics for Higher-Order Bernstein Bézier Elements
Element Quality Metrics for Higher-Order Bernstein Bézier Elements Luke Engvall and John A. Evans Abstract In this note, we review the interpolation theory for curvilinear finite elements originally derived
More informationRight-angled hyperbolic polyhedra
Right-angled hyperbolic polyhedra Andrei Vesnin Sobolev Institute of Mathematics, Novosibirsk, Russia The Fourth Geometry Meeting dedicated to the centenary of A. D. Alexandrov August 20 24, 2012 A.D.Alexandrov:
More informationTiling 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 information1 Unstable periodic discrete minimal surfaces
1 Unstable periodic discrete minimal surfaces Konrad Polthier Technische Universität Berlin, Institut für Mathematik polthier@math.tu-berlin.de Summary. In this paper we define the new alignment energy
More information