Parallel graph-grammar-based algorithm for the longest-edge refinement of triangular meshes and the pollution simulations in Lesser Poland area

In this paper, we propose parallel graph-grammar-based algorithm for the longest-edge refinements and the pollution simulations in Lesser Poland area. We introduce graph-grammar productions for Rivara’s longest-edged algorithm for the local refinement of unstructured triangular meshes. We utilize the hyper-graph to represent the computational mesh and the graph-grammar productions to express the longest-edge mesh refinement algorithm. The parallelism in the original Rivara’s longest edge refinement algorithm is obtained by processing different longest edge refinement paths in different three ads. Our graph-grammar-based algorithm allows for additional parallelization within a single longest-edge refinement path. The graph-grammar-based algorithm automatically guarantees the validity and conformity of the generated mesh; it prevents the generation of duplicated nodes and edges, elongated elements with Jacobians converging to zero, and removes all the hanging nodes automatically from the mesh. We test the algorithm on generating a surface mesh based on a topographic data of Lesser Poland area. The graph-grammar productions also generate the layers of prismatic three-dimensional elements on top of the triangular mesh, and they break each prismatic element into three tetrahedral elements. Next, we propose graph-grammar productions generating element matrices and right-hand-side vectors for each tetrahedral element. We utilize the Streamline Upwind Petrov–Galerkin (SUPG) stabilization for the pollution propagation simulations in Lesser Poland area. We use the advection–diffusion-reaction model, the Crank–Nicolson time integration scheme, and the graph-grammar-based interface to the GMRES solver.

[1]  Annegret Habel,et al.  May we introduce to you: hyperedge replacement , 1986, Graph-Grammars and Their Application to Computer Science.

[2]  Ewa Grabska,et al.  Graph Transformations for Modeling hp-Adaptive Finite Element Method with Triangular Elements , 2008, ICCS.

[3]  Keshav Pingali,et al.  Graph Grammar based Multi-thread Multi-frontal Direct Solver with Galois Scheduler , 2014, ICCS.

[4]  Keshav Pingali,et al.  Priority Queues Are Not Good Concurrent Priority Schedulers , 2015, Euro-Par.

[5]  M. Rivara Lepp-bisection algorithms, applications and mathematical properties , 2009 .

[6]  Keshav Pingali,et al.  The tao of parallelism in algorithms , 2011, PLDI '11.

[7]  Peter Sanders,et al.  Network Flow-Based Refinement for Multilevel Hypergraph Partitioning , 2019, ACM J. Exp. Algorithmics.

[8]  D. N. Herting,et al.  Finite elements: Computational aspects: Vol. III, by G.F. Carey and J. Tinsley Oden, Prentice-Hall, Englewood Cliffs, NJ, 1984 , 1985 .

[9]  M. Rivara NEW LONGEST-EDGE ALGORITHMS FOR THE REFINEMENT AND/OR IMPROVEMENT OF UNSTRUCTURED TRIANGULATIONS , 1997 .

[10]  M. Rivara,et al.  Multithread parallelization of Lepp-bisection algorithms , 2012 .

[11]  Keshav Pingali,et al.  Ordered vs. unordered: a comparison of parallelism and work-efficiency in irregular algorithms , 2011, PPoPP '11.

[12]  Igor Kossaczký A recursive approach to local mesh refinement in two and three dimensions , 1994 .

[13]  Keshav Pingali,et al.  Optimistic parallelism requires abstractions , 2007, PLDI '07.

[14]  Leszek Siwik,et al.  A Massively Parallel Algorithm for the Three-Dimensional Navier-Stokes-Boussinesq Simulations of the Atmospheric Phenomena , 2020, ICCS.

[15]  Christophe Geuzaine,et al.  Gmsh: A 3‐D finite element mesh generator with built‐in pre‐ and post‐processing facilities , 2009 .

[16]  Victor M. Calo,et al.  Residual-based multiscale turbulence modeling: Finite volume simulations of bypass transition , 2005 .

[17]  M. Rivara Algorithms for refining triangular grids suitable for adaptive and multigrid techniques , 1984 .

[18]  M. Rivara Mesh Refinement Processes Based on the Generalized Bisection of Simplices , 1984 .

[19]  Annegret Habel,et al.  Some Structural Aspects of Hypergraph Languages Generated by Hyperedge Replacement , 1987, STACS.

[20]  L. Demkowicz One and two dimensional elliptic and Maxwell problems , 2006 .

[21]  Agustí Pérez-Foguet,et al.  Adaptive finite element simulation of stack pollutant emissions over complex terrains , 2013 .

[22]  T. Hughes,et al.  A new finite element formulation for computational fluid dynamics: VI. Convergence analysis of the generalized SUPG formulation for linear time-dependent multi-dimensional advective-diffusive systems , 1987 .

[23]  Anna Paszynska,et al.  Graph Transformations for Modeling Parallel hp-Adaptive Finite Element Method , 2007, PPAM.

[24]  Ewa Grabska,et al.  Graph Transformations for Modeling hp-Adaptive Finite Element Method with Mixed Triangular and Rectangular Elements , 2009, ICCS.

[25]  H. Bruyninckx The European Environment Agency , 2017 .