Implementation of Polygonal Mesh Refinement in MATLAB

Yue Yu

DOI: https://doi.org/10.48550/arXiv.2101.03456

2021-01-10

Abstract:We present a simple and efficient MATLAB implementation of the local refinement for polygonal meshes. The purpose of this implementation is primarily educational, especially on the teaching of adaptive virtual element methods.

Numerical Analysis

What problem does this paper attempt to address?

A new quality preserving polygonal mesh refinement algorithm for Virtual Element Methods

Stefano Berrone,Alessandro D'Auria

DOI: https://doi.org/10.48550/arXiv.2112.09888

2021-12-18

Numerical Analysis

Abstract:Mesh adaptivity is a useful tool for efficient solution to partial differential equations in very complex geometries. In the present paper we discuss the use of polygonal mesh refinement in order to tackle two common issues: first, adaptively refine a provided good quality polygonal mesh preserving quality, second, improve the quality of a coarse poor quality polygonal mesh during the refinement process on very complex domains. For finite element methods and triangular meshes, convergence of a posteriori mesh refinement algorithms and optimality properties have been widely investigated, whereas convergence and optimality are still open problems for polygonal adaptive methods. In this article, we propose a new refinement method for convex cells with the aim of introducing some properties useful to tackle convergence and optimality for adaptive methods. The key issues in refining convex general polygons are: a refinement dependent only on the marked cells for refinement at each refinement step; a partial quality improvement, or, at least, a non degenerate quality of the mesh during the refinement iterations; a bound of the number of unknowns of the discrete problem with respect to the number of the cells in the mesh. Although these properties are quite common for refinement algorithms of triangular meshes, these issues are still open problems for polygonal meshes
A unified approach for fairing arbitrary polygonal meshes

Guiqing Li,Hujun Bao,Weiyin Ma

DOI: https://doi.org/10.1016/j.gmod.2004.03.001

IF: 1.094

2004-01-01

Graphical Models

Abstract:Arbitrary polygonal meshes play an important role in 3D model representation and surface modeling. Most of the previous smoothing approaches are however, mainly designed for triangular meshes. This paper proposes a unified approach to address fairing problems of arbitrary polygonal meshes. The approach is built upon an assumption that all face centroids of a smoothed mesh locate on its noise-free surface and remain unchanged during the smoothing process in order to prevent the resulted mesh from shrinkage and preserve surface features. Boundary smoothing is considered in a similar way as well but the process is independent to mesh processing such that two, noisy meshes sharing a boundary still adjoin continuously after smoothing. Finally, a Taubin-like inflation operator is employed to further improve the ability of anti-shrinkage. Although very simple iterative equation is employed in the smoothing procedure, experiments demonstrate that the approach is able to not only effectively address the problem of shrinkage and distortion but also preserve features very well.
A Parallel Local Reconnection Approach for Tetrahedral Mesh Improvement

Mengmeng Shang,Chaoyan Zhu,Jianjun Chen,Zhoufang Xiao,Yao Zheng

DOI: https://doi.org/10.1016/j.proeng.2016.11.062

2016-01-01

Procedia Engineering

Abstract:A multi-threaded parallel local reconnection algorithm is proposed for tetrahedral meshes. It defines a feature point within the region involved in each operation, and sorts the features points along a Hilbert curve. The decomposition of this Hilbert curve results in a load-balanced distribution of local operations. Meanwhile, the regions of concurrently executed local operations are separated far away, such that the possibility of interference is reduced to a very low level. Finally, a parallel mesh improver is developed by combining the proposed algorithm with a parallel mesh smoothing algorithm, and its effectiveness and efficiency is verified in various numerical experiments.
Tetrahedral mesh improvement by shell transformation

Jianjun Chen,Jianjing Zheng,Yao Zheng,Zhoufang Xiao,Hang Si,Yufeng Yao

DOI: https://doi.org/10.1007/s00366-016-0480-z

2016-01-01

Abstract:Existing flips for tetrahedral meshes simply make a selection from a few possible configurations within a single shell (i.e., a polyhedron that can be filled up with a mesh composed of a set of elements that meet each other at one edge), and their effectiveness is usually confined. A new topological operation for tetrahedral meshes named shell transformation is proposed. Its recursive callings execute a sequence of shell transformations on neighboring shells, acting like composite edge removal transformations. Such topological transformations are able to perform on a much larger element set than that of a single flip, thereby leading the way towards a better local optimum solution. Hence, a new mesh improvement algorithm is developed by combining this recursive scheme with other schemes, including smoothing, point insertion, and point suppression. Numerical experiments reveal that the proposed algorithm can well balance some stringent and yet sometimes even conflict requirements of mesh improvement, i.e., resulting in high-quality meshes and reducing computing time at the same time. Therefore, it can be used for mesh quality improvement tasks involving millions of elements, in which it is essential not only to generate high-quality meshes, but also to reduce total computational time for mesh improvement.
Machine learning based refinement strategies for polyhedral grids with applications to virtual element and polyhedral discontinuous Galerkin methods

P.F. Antonietti,F. Dassi,E. Manuzzi

DOI: https://doi.org/10.1016/j.jcp.2022.111531

IF: 4.645

2022-11-15

Journal of Computational Physics

Abstract:We propose two new strategies based on Machine Learning techniques to handle polyhedral grid refinement, to be possibly employed within an adaptive framework. The first one employs the k-means clustering algorithm to partition the points of the polyhedron to be refined. This strategy is a variation of the well known Centroidal Voronoi Tessellation. The second one employs Convolutional Neural Networks to classify the “shape” of an element so that “ad-hoc” refinement criteria can be defined. This strategy can be used to enhance existing refinement strategies, including the k-means strategy, at a low online computational cost. We test the proposed algorithms considering two families of finite element methods that support arbitrarily shaped polyhedral elements, namely the Virtual Element Method (VEM) and the Polygonal Discontinuous Galerkin (PolyDG) method. We demonstrate that these strategies do preserve the structure and the quality of the underlaying grids, reducing the overall computational cost and mesh complexity.

physics, mathematical,computer science, interdisciplinary applications
Flattenable Mesh Processing by Controllable Laplacian Evolution

Hongwei Lin,Yunbo Zhang,Charlie C. L. Wang,Shuming Gao

DOI: https://doi.org/10.1115/detc2010-28301

2010-01-01

Abstract:Models represented by polygonal meshes have been more and more widely used in CAD/CAM systems. In sheet manufacturing industries, the flattenability of a model is very important. Prior methods for processing the flattenability of a mesh surface usually employ a constrained optimization framework, which takes the positions of all its non-boundary vertices as variables in computation. For a mesh surface with hundred thousands of vertices, solving such an optimization is very time-consuming, and may exceed the capacity of main memory. In this paper, we develop a controllable evolution method to process the flattenability of a given mesh patch. It decouples the global optimization problem into a sequence of local controllable evolution steps, each of which has only one variable. Therefore, mesh surfaces with large number of vertices can be processed.
Refinement Computations of Electromagnetic Fields Using FE and Meshless Methods

SL Ho,SY Yang,HC Wong,EWC Lo,GZ Ni

DOI: https://doi.org/10.1109/tmag.2005.844553

IF: 1.848

2005-01-01

IEEE Transactions on Magnetics

Abstract:A refinement algorithm for electromagnetic field computations using a combination of finite element and meshless methods is introduced. Bridging scales are used to separate the finite element and meshless shape functions to make the refinement hierarchical and to uphold the mathematical properties such as consistency and linear independence for all the bases. To facilitate the application of the proposed algorithm, details about the node addition, requirements for the node distribution, and relationships between the finite element and meshless shape functions, as well as the determination of the stop criterion are also fully addressed. Primary numerical results are reported to demonstrate and validate the applicability and advantages of the proposed algorithm over traditional ones.
Noise Removal Algorithm for Polygonal Meshes

XG Liu,HJ Bao,QS Peng

IF: 4.269

2000-01-01

Progress in Natural Science Materials International

Abstract:An effective noise removal algorithm which does not cause the shrinkage of surface is presented. By introducing the constraints of keeping all triangle barycenters invariant at each smoothing step, the algorithm turns the problem of mesh fairing into a constrained minimization of the energy functional, which is then solved efficiently by our iterative method. Experimental results demonstrate that our algorithm not only can preserve the good shape of the original surface while quickly removing the noise, hut also has the advantages of fast convergence, low computational cost and stable performance.
Refinement of polygonal grids using Convolutional Neural Networks with applications to polygonal Discontinuous Galerkin and Virtual Element methods

P. F. Antonietti,E. Manuzzi

DOI: https://doi.org/10.1016/j.jcp.2021.110900

2021-10-02

Abstract:We propose new strategies to handle polygonal grids refinement based on Convolutional Neural Networks (CNNs). We show that CNNs can be successfully employed to identify correctly the "shape" of a polygonal element so as to design suitable refinement criteria to be possibly employed within adaptive refinement strategies. We propose two refinement strategies that exploit the use of CNNs to classify elements' shape, at a low computational cost. We test the proposed idea considering two families of finite element methods that support arbitrarily shaped polygonal elements, namely Polygonal Discontinuous Galerkin (PolyDG) methods and Virtual Element Methods (VEMs). We demonstrate that the proposed algorithms can greatly improve the performance of the discretization schemes both in terms of accuracy and quality of the underlying grids. Moreover, since the training phase is performed off-line and is independent of the differential model the overall computational costs are kept low.

Numerical Analysis,Machine Learning
A Nyström-based finite element method on polygonal elements

Akash Anand,Jeffrey S. Ovall,Steffen Weisser

DOI: https://doi.org/10.1016/j.camwa.2018.03.007

2017-08-24

Abstract:We consider families of finite elements on polygonal meshes, that are defined implicitly on each mesh cell as solutions of local Poisson problems with polynomial data. Functions in the local space on each mesh cell are evaluated via Nyström discretizations of associated integral equations, allowing for curvilinear polygons and non-polynomial boundary data. Several experiments demonstrate the approximation quality of interpolated functions in these spaces.

Numerical Analysis
Effective polygonal mesh generation and refinement for VEM

Stefano Berrone,Fabio Vicini

2024-03-15

Abstract:In the present work we introduce a novel refinement algorithm for two-dimensional elliptic partial differential equations discretized with Virtual Element Method (VEM). The algorithm improves the numerical solution accuracy and the mesh quality through a controlled refinement strategy applied to the generic polygonal elements of the domain tessellation. The numerical results show that the outlined strategy proves to be versatile and applicable to any two-dimensional problem where polygonal meshes offer advantages. In particular, we focus on the simulation of flow in fractured media, specifically using the Discrete Fracture Network (DFN) model. A residual a-posteriori error estimator tailored for the DFN case is employed. We chose this particular application to emphasize the effectiveness of the algorithm in handling complex geometries. All the numerical tests demonstrate optimal convergence rates for all the tested VEM orders.

Numerical Analysis
Local Refinement and Multilevel Preconditioning: Implementation and Numerical Experiments

Burak Aksoylu,Stephen Bond,Michael Holst

DOI: https://doi.org/10.48550/arXiv.1001.1370

2010-01-08

Numerical Analysis

Abstract:In this paper, we examine a number of additive and multiplicative multilevel iterative methods and preconditioners in the setting of two-dimensional local mesh refinement. While standard multilevel methods are effective for uniform refinement-based discretizations of elliptic equations, they tend to be less effective for algebraic systems which arise from discretizations on locally refined meshes, losing their optimal behavior in both storage and computational complexity. Our primary focus here is on BPX-style additive and multiplicative multilevel preconditioners, and on various stabilizations of the additive and multiplicative hierarchical basis method (HB), and their use in the local mesh refinement setting. In this article, we describe in detail the implementation of these types of algorithms, including detailed discussions of the datastructures and traversal algorithms we employ for obtaining optimal storage and computational complexity in our implementations. We show how each of the algorithms can be implemented using standard datatypes available in languages such as C and FORTRAN, so that the resulting algorithms have optimal (linear) storage requirements, and so that the resulting multilevel method or preconditioner can be applied with optimal (linear) computational costs. Our implementations are performed in both C and MATLAB using the Finite Element ToolKit (FETK), an open source finite element software package. We finish the paper with a sequence of numerical experiments illustrating the effectiveness of a number of BPX and stabilized HB variants for several examples requiring local refinement.
A coarsening algorithm on adaptive red-green-blue refined meshes

Stefan A. Funken,Anja Schmidt

DOI: https://doi.org/10.1007/s11075-020-01003-7

IF: 2.37

2020-10-01

Numerical Algorithms

Abstract:Abstract Adaptive meshing is a fundamental component of adaptive finite element methods. This includes refining and coarsening meshes locally. In this work, we are concerned with the red-green-blue refinement strategy in two dimensions and its counterpart-coarsening. In general, coarsening algorithms are mostly based on an explicitly given refinement history. In this work, we present a coarsening algorithm on adaptive red-green-blue meshes in two dimensions without explicitly knowing the refinement history. To this end, we examine the local structure of these meshes, find an easy-to-verify criterion to adaptively coarsen red-green-blue meshes, and prove that this criterion generates meshes with the desired properties. We present a MATLAB implementation built on the red-green-blue refinement routine of the -package (Funken and Schmidt 2018, 2019).

mathematics, applied
Improvement of boundary aberration and computation efficiency in meshless method

Li MA,Rong-hao BAO,Shuang-lian WANG,Yi-mu GUO

DOI: https://doi.org/10.3785/j.issn.1008-973X.2007.06.019

2007-01-01

Abstract:Two usual methods for constructing the shape functions of meshless methods, smooth particles hydrodynamics (SPH) and moving least square (MLS), were analyzed and compared. The reason of boundary aberration in SPH was indicated, and the virtual points method which assigned extra points outside the boundary was presented to control the boundary aberration. Allocation method and element-free Galerkin method were used to resolve the one-dimensional and two-dimensional differential equations. Computation results showed that Lagrange multiplier method has better accuracy than point interpolation method in dealing with boundary conditions. Cantilever beam example illustrated the influence of support radius and Gauss integral point on accuracy and cost. The result indicates that single Gauss integral point is affected by the variation of support radius, however higher order Gauss integral points reduce the sensitivity to the variation of support radius and improve the computational accuracy, while the computation cost increases.
An Ultralight Geometry Processing Library for Parallel Mesh Refinement

Wang Bohan,Chen Bo,Yu Kaixin,Xie Lijun,Chen Jianjun

DOI: https://doi.org/10.1186/s42774-021-00093-8

2022-01-01

Advances in Aerodynamics

Abstract:In applications such as parallel mesh refinement, it remains a challenging issue to ensure the refined surface respects the original Computer-Aided Design (CAD) model accurately. In this paper, an ultralight geometry processing library is developed to resolve this issue effectively and efficiently. Here, we say the kernel is ultralight because it has a very small set of data-structures and algorithms by comparison with industrial-level geometry kernels. Within the library, a simplified surface boundary representation (B-rep) and a radial edge structure are developed respectively to depict the geometry model and the surface mesh, plus hash tables that record the connections between the geometry model and the surface mesh. Based on these data structures, a set of efficient algorithms are developed, which initializes the connection tables, projects a point back to the original geometry, etc. With these data-structure and algorithmic infrastructures set up, the callings of eight well-designed Application Programming Interfaces (APIs) are powerful enough to enable the parallel mesh refinement algorithm outputs a mesh respecting the input CAD model accurately. Numerical experiments will be finally presented to evaluate the performance of the overall parallel mesh refinement algorithm and the algorithms in relation with the developed library.
The lowest-order Neural Approximated Virtual Element Method on polygonal elements

Stefano Berrone,Moreno Pintore,Gioana Teora

2024-09-24

Abstract:The lowest-order Neural Approximated Virtual Element Method on polygonal elements is proposed here. This method employs a neural network to locally approximate the Virtual Element basis functions, thereby eliminating issues concerning stabilization and projection operators, which are the key components of the standard Virtual Element Method. We propose different training strategies for the neural network training, each correlated by the theoretical justification and with a different level of accuracy. Several numerical experiments are proposed to validate our procedure on general polygonal meshes and demonstrate the advantages of the proposed method across different problem formulations, particularly in cases where the heavy usage of projection and stabilization terms may represent challenges for the standard version of the method. Particular attention is reserved to triangular meshes with hanging nodes which assume a central role in many virtual element applications.

Numerical Analysis
Quasi-optimal mesh generation for the virtual element method: A fully adaptive remeshing procedure

Daniel van Huyssteen,Felipe Lopez Rivarola,Guillermo Etse,Paul Steinmann

2024-07-19

Abstract:The mesh flexibility offered by the virtual element method through the permission of arbitrary element geometries, and the seamless incorporation of `hanging' nodes, has made the method increasingly attractive in the context of adaptive remeshing. There exists a healthy literature concerning error estimation and adaptive refinement techniques for virtual elements while the topic of adaptive coarsening (i.e. de-refinement) is in its infancy. The creation of a quasi-optimal mesh is based on the principle of quasi-even error distribution over the elements which inherently relies on localized refinement and coarsening techniques. Thus, necessitating a fully adaptive remeshing procedure. In this work a novel fully adaptive remeshing procedure for the virtual element method is presented. Additionally, novel procedures are proposed for the identification of elements qualifying for refinement or coarsening based on user-defined targets. Specifically, a target percentage error, target number of elements, or target number of nodes can be selected. Numerical results demonstrate that the adaptive remeshing procedure can meet any prescribed target while creating a quasi-optimal mesh. The proposed fully adaptive procedure is of particular interest in engineering applications requiring an efficient simulation of a given accuracy, or desiring a simulation with the maximum possible accuracy for a given computational constraint.

Numerical Analysis
An Edge-Based Anisotropic Mesh Refinement Algorithm and Its Application to Interface Problems

Duan Wang,Ruo Li,Ningning Yan

DOI: https://doi.org/10.4208/cicp.210709.121109a

IF: 3.791

2010-01-01

Communications in Computational Physics

Abstract:Based on an error estimate in terms of element edge vectors on arbitrary unstructured simplex meshes, we propose a new edge-based anisotropic mesh refinement algorithm. As the mesh adaptation indicator, the error estimate involves only the gradient of error rather than higher order derivatives. The preferred refinement edge is chosen to reduce the maximal term in the error estimate. The algorithm is implemented in both two- and three-dimensional cases, and applied to the singular function interpolation and the elliptic interface problem. The numerical results demonstrate that the convergence order obtained by using the proposed anisotropic mesh refinement algorithm can be higher than that given by the isotropic one.
On mesh coarsening procedures for the virtual element method

Daniel van Huyssteen,Felipe Lopez Rivarola,Guillermo Etse,Paul Steinmann

2023-08-15

Abstract:In the context of adaptive remeshing, the virtual element method provides significant advantages over the finite element method. The attractive features of the virtual element method, such as the permission of arbitrary element geometries, and the seamless permission of 'hanging' nodes, have inspired many works concerning error estimation and adaptivity. However, these works have primarily focused on adaptive refinement techniques with little attention paid to adaptive coarsening (i.e. de-refinement) techniques that are required for the development of fully adaptive remeshing procedures. In this work novel indicators are proposed for the identification of patches/clusters of elements to be coarsened, along with a novel procedure to perform the coarsening. The indicators are computed over prospective patches of elements rather than on individual elements to identify the most suitable combinations of elements to coarsen. The coarsening procedure is robust and suitable for meshes of structured and unstructured/Voronoi elements. Numerical results demonstrate the high degree of efficacy of the proposed coarsening procedures and sensible mesh evolution during the coarsening process. It is demonstrated that critical mesh geometries, such as non-convex corners and holes, are preserved during coarsening, and that meshes remain fine in regions of interest to engineers, such as near singularities.

Numerical Analysis,Analysis of PDEs
Local Element Operations for Curved Simplex Meshes

Andrew Shi,Per-Olof Persson

DOI: https://doi.org/10.1002/nme.7379

2023-10-17

Abstract:Mesh optimization procedures are generally a combination of node smoothing and discrete operations which affect a small number of elements to improve the quality of the overall mesh. These procedures are useful as a post-processing step in mesh generation procedures and in applications such as fluid simulations with severely deforming domains. In order to perform high-order mesh optimization, these ingredients must also be extended to high-order (curved) meshes. In this work, we present a method to perform local element operations on curved meshes. The mesh operations discussed in this work are edge/face swaps, edge collapses, and edge splitting (more generally refinement) for triangular and tetrahedral meshes. These local operations are performed by first identifying the patch of elements which contain the edge/face being acted on, performing the operation as a straight-sided one by placing the high-order nodes via an isoparametric mapping from the master element, and smoothing the high-order nodes on the elements in the patch by minimizing a Jacobian-based high-order mesh distortion measure. Since the initial straight-sided guess from the placement of the nodes via the isoparametric mapping frequently results in invalid elements, the distortion measure must be regularized which allows for mesh untangling for the optimization to succeed. We present several examples in 2D and 3D to demonstrate these local operations and how they can be combined with a high-order node smoothing procedure to maintain mesh quality when faced with severe deformations.

Numerical Analysis

Implementation of Polygonal Mesh Refinement in MATLAB

A new quality preserving polygonal mesh refinement algorithm for Virtual Element Methods

A unified approach for fairing arbitrary polygonal meshes

A Parallel Local Reconnection Approach for Tetrahedral Mesh Improvement

Tetrahedral mesh improvement by shell transformation

Machine learning based refinement strategies for polyhedral grids with applications to virtual element and polyhedral discontinuous Galerkin methods

Flattenable Mesh Processing by Controllable Laplacian Evolution

Refinement Computations of Electromagnetic Fields Using FE and Meshless Methods

Noise Removal Algorithm for Polygonal Meshes

Refinement of polygonal grids using Convolutional Neural Networks with applications to polygonal Discontinuous Galerkin and Virtual Element methods

A Nyström-based finite element method on polygonal elements

Effective polygonal mesh generation and refinement for VEM

Local Refinement and Multilevel Preconditioning: Implementation and Numerical Experiments

A coarsening algorithm on adaptive red-green-blue refined meshes

Improvement of boundary aberration and computation efficiency in meshless method

An Ultralight Geometry Processing Library for Parallel Mesh Refinement

The lowest-order Neural Approximated Virtual Element Method on polygonal elements

Quasi-optimal mesh generation for the virtual element method: A fully adaptive remeshing procedure

An Edge-Based Anisotropic Mesh Refinement Algorithm and Its Application to Interface Problems

On mesh coarsening procedures for the virtual element method

Local Element Operations for Curved Simplex Meshes