Paving the path towards automatic hexahedral mesh generation

Esta tesis versa sobre el desarrollo de las tecnologías para la generación de mallas de hexaedros. El proceso de generar una malla de hexaedros no es automático y su generación requiere varias horas te trabajo de un ingeniero especializado. Por lo tanto, es importante desarrollar herramientas que faciliten dicho proceso de generación. Con este fin, se presenta y desarrolla un método de proyección de mallas, una técnica de sweeping o barrido, un algoritmo para la obtención de mallas por bloques, y un entorno de generación de mallas. Las implementaciones más competitivas del método de sweeping utilizan técnicas de proyección de mallas basadas en métodos afines. Los métodos afines más habituales presentan varios problemas relacionados con la obtención de sistemas de ecuaciones normales de rango deficiente. Para solucionar dichos problemas se presenta y analiza un nuevo método afín que depende de dos parámetros vectoriales. Además, se detalla un procedimiento automático para la selección de dichos vectores. El método de proyección resultante preserva la forma de las mallas proyectadas. Esta proyección es incorporada también en una nueva herramienta de sweeping. Dicha herramienta genera capas de nodos internos que respetan la curvatura de las superficies inicial y final. La herramienta de sweeping es capaz de mallar geometrías de extrusión definidas por trayectorias curvas, secciones no constantes a lo largo del eje de sweeping, y superficies inicial y final con diferente forma y curvatura.En las últimas décadas se han propuesto varios ataques para la generación automática de mallas de hexahedros. Sin embargo, todavía no existe un algoritmo rápido y robusto que genere automáticamente mallas de hexaedros de alta calidad. Se propone un nuevo ataque para la generación de mallas por bloques mediante la representación de la geometría y la topología del dual de una malla de hexaedros. En dicho ataque, primero se genera una malla grosera de tetraedros. Después, varió polígonos planos se añaden al interior de los elementos de la malla grosera inicial. Dichos polígonos se denotan como contribuciones duales locales y representan una versión discreta del dual de una malla de hexaedros. En el último paso, la malla por bloques se obtiene como el dual de la representación del dual generada. El algoritmo de generación de mallas por bloques es aplicado a geometrías que presentan diferentes características geométricas como son superficies planas, superficies curvas, configuraciones delgadas, agujeros, y vértices con valencia mayor que tres.Las mallas se generan habitualmente con la ayuda de entornos interactivos que integran una interfaz CAD y varios algoritmos de generación de mallas. Se presenta un nuevo entorno de generación de mallas especializado en la generación de cuadriláteros y hexaedros. Este entorno proporciona la tecnología necesaria para implementar les técnicas de generación de mallas de hexaedros presentadas en esta tesis.This thesis deals with the development of hexahedral mesh generation technology. The process of generating hexahedral meshes is not fully automatic and it is a time consuming task. Therefore, it is important to develop tools that facilitate the generation of hexahedral meshes. To this end, a mesh projection method, a sweeping technique, a block-meshing algorithm, and an interactive mesh generation environment are presented and developed. Competitive implementations of the sweeping method use mesh projection techniques based on affine methods. Standard affine methods have several drawbacks related to the statement of rank deficient sets of normal equations. To overcome these drawbacks a new affine method that depends on two vector parameters is presented and analyzed. Moreover, an automatic procedure that selects these two vector parameters is detailed. The resulting projection procedure preserves the shape of projected meshes. Then, this procedure is incorporated in a new sweeping tool. This tool generates inner layers of nodes that preserve the curvature of the cap surfaces. The sweeping tool is able to mesh extrusion geometries defined by non-linear sweeping trajectories, non-constant cross sections along the sweep axis, non-parallel cap surfaces, and cap surfaces with different shape and curvature. In the last decades, several general-purpose approaches to generate automatically hexahedral meshes have been proposed. However, a fast and robust algorithm that automatically generates high-quality hexahedral meshes is not available. A novel approach for block meshing by representing the geometry and the topology of a hexahedral mesh is presented. The block-meshing algorithm first generates an initial coarse mesh of tetrahedral elements. Second, several planar polygons are added inside the elements of the initial coarse mesh. These polygons are referred as local dual contributions and represent a discrete version of the dual of a hexahedral mesh. Finally, the dual representation is dualized to obtain the final block mesh. The block-meshing algorithm is applied to mesh geometries that present different geometrical characteristics such as planar surfaces, curved surfaces, thin configurations, holes, and vertices with valence greater than three.Meshes are usually generated with the help of interactive environments that integrate a CAD interface and several meshing algorithms. An overview of a new mesh generation environment focused in quadrilateral and hexahedral mesh generation is presented. This environment provides the technology required to implement the hexahedral meshing techniques presented in this thesis.Postprint (published version

Generation of curved high-order meshes with optimal quality and geometric accuracy

We present a novel methodology to generate curved high-order meshes featuring optimal mesh quality and geometric accuracy. The proposed technique combines a distortion measure and a geometric L2-disparity measure into a single objective function. While the element distortion term takes into account the mesh quality, the L2-disparity term takes into account the geometric error introduced by the mesh approximation to the target geometry. The proposed technique has several advantages. First, we are not restricted to interpolative meshes and therefore, the resulting mesh approximates the target domain in a non-interpolative way, further increasing the geometric accuracy. Second, we are able to generate a series of meshes that converge to the actual geometry with expected rate while obtaining high-quality elements. Third, we show that the proposed technique is robust enough to handle real-case geometries that contain gaps between adjacent entities.Peer ReviewedPostprint (published version

Local bisection for conformal refinement of unstructured 4D simplicial meshes

We present a conformal bisection procedure for local refinement of 4D unstructured simplicial meshes with bounded minimum shape quality. Specifically, we propose a recursive refine-to-conformity procedure in two stages, based on marking bisection edges on different priority levels and defining specific refinement templates. Two successive applications of the first stage ensure that any 4D unstructured mesh can be conformingly refined. In the second stage, the successive refinements lead to a cycle in the number of generated similarity classes and thus, we can ensure a bound over the minimum shape quality. In the examples, we check that after successive refinement the mesh quality does not degenerate. Moreover, we refine a 4D unstructured mesh and a space-time mesh (3D + 1D) representation of a moving object

An interactive mesh generation environment for geometry-based simulations. Progress updated

We presented an overview of an interactive mesh generation environment in our previous work that it is being developed in the Laboratori de Càlcul Numèric (LaCàN). We started the development of the software in order to unify available legacy code, new developments and research algorithms in only one mesh generation package. This paper presents an update to the last environment overview after significant implementation and conceptual development. We provide a brief summary of: New features, as submapping and data output suitable for geometry-based methods as; current developments, as GUI improvement and 3D modelling tools; and future features, such as the command interpreter. Finally, we conclude that chosen mesh generation paradigm and software engineering concepts have allowed us to improve and scale the environment since last revision

A distortion measure to validate and generate curved high-order meshes on CAD surfaces with independence of parameterization

This is the accepted version of the following article: [Gargallo-Peiró, A., Roca, X., Peraire, J., and Sarrate, J. (2016) A distortion measure to validate and generate curved high-order meshes on CAD surfaces with independence of parameterization. Int. J. Numer. Meth. Engng, 106: 1100–1130. doi: 10.1002/nme.5162], which has been published in final form at http://onlinelibrary.wiley.com/doi/10.1002/nme.5162/abstractA framework to validate and generate curved nodal high-order meshes on Computer-Aided Design (CAD) surfaces is presented. The proposed framework is of major interest to generate meshes suitable for thin-shell and 3D finite element analysis with unstructured high-order methods. First, we define a distortion (quality) measure for high-order meshes on parameterized surfaces that we prove to be independent of the surface parameterization. Second, we derive a smoothing and untangling procedure based on the minimization of a regularization of the proposed distortion measure. The minimization is performed in terms of the parametric coordinates of the nodes to enforce that the nodes slide on the surfaces. Moreover, the proposed algorithm repairs invalid curved meshes (untangling), deals with arbitrary polynomial degrees (high-order), and handles with low-quality CAD parameterizations (independence of parameterization). Third, we use the optimization procedure to generate curved nodal high-order surface meshes by means of an a posteriori approach. Given a linear mesh, we increase the polynomial degree of the elements, curve them to match the geometry, and optimize the location of the nodes to ensure mesh validity. Finally, we present several examples to demonstrate the features of the optimization procedure, and to illustrate the surface mesh generation process.Peer ReviewedPostprint (author's final draft

Checking and improving the geometric accuracy of non-interpolating curved high-order meshes

Peer ReviewedPostprint (author's final draft

Unstructured and semi-structured hexahedral mesh generation methods

Discretization techniques such as the finite element method, the finite volume method or the discontinuous Galerkin method are the most used simulation techniques in ap- plied sciences and technology. These methods rely on a spatial discretization adapted to the geometry and to the prescribed distribution of element size. Several fast and robust algorithms have been developed to generate triangular and tetrahedral meshes. In these methods local connectivity modifications are a crucial step. Nevertheless, in hexahedral meshes the connectivity modifications propagate through the mesh. In this sense, hexahedral meshes are more constrained and therefore, more difficult to gener- ate. However, in many applications such as boundary layers in computational fluid dy- namics or composite material in structural analysis hexahedral meshes are preferred. In this work we present a survey of developed methods for generating structured and unstructured hexahedral meshes.Peer ReviewedPostprint (published version

Size preserving mesh generation in adaptivity processes

It is well known that the variations of the element size have to be controlled in order to generate a high-quality mesh. Hence, several techniques have been developed to limit the gradient of the element size. Although these methods allow generating high-quality meshes, the obtained discretizations do not always reproduce the prescribed size function. Specifically, small elements may not be generated in a region where small element size is prescribed. This is critical for many practical simulations, where small elements are needed to reduce the error of the numerical simulation. To solve this issue, we present the novel size-preserving technique to control the mesh size function prescribed at the vertices of a background mesh. The result is a new size function that ensures a high-quality mesh with all the elements smaller or equal to the prescribed element size. That is, we ensure that the new mesh handles at least one element of the correct size at each local minima of the size function. In addition, the gradient of the size function is limited to obtain a high-quality mesh. Two direct applications are presented. First, we show that we can reduce the number of iterations to converge an adaptive process, since we do not need additional iterations to generate a valid mesh. Second, the size-preserving approach allows to generate quadri- lateral meshes that correctly preserves the prescribed element size.Peer ReviewedPostprint (published version

A surface mesh smoothing and untangling method independent of the CAD parameterization

A method to optimize triangular and quadrilateral meshes on parameterized surfaces is proposed. The optimization procedure relocates the nodes on the surface to improve the quality (smooth) and ensures that the elements are not inverted (untangle). We detail how to express any measure for planar elements in terms of the parametric coordinates of the nodes. The extended measures can be used to check the quality and validity of a surface mesh. Then, we detail how to optimize any Jacobian-based distortion measure to obtain smoothed and untangled meshes with the nodes on the surface. We prove that this method is independent of the surface parameterization. Thus, it can optimize meshes on CAD surfaces defined by low-quality parameterizations. The examples show that the method can optimize meshes composed by a large number of inverted elements. Finally, the method can be extended to obtain high-order meshes with the nodes on the CAD surfaces.Peer ReviewedPostprint (author's final draft

Defining an2-disparity measure to check and improve the geometric accuracy of noninterpolating curved high-order meshes

We define an2-disparity measure between curved high-order meshes and parameterized manifolds in terms of an2norm. The main application of the proposed definition is to measure and improve the distance between a curved high-order mesh and a target parameterized curve or surface. The approach allows considering meshes with the nodes on top of the curve or surface (interpolative), or floating freely in the physical space (non-interpolative). To compute the disparity measure, the average of the squared point-wise differences is minimized in terms of the nodal coordinates of an auxiliary parametric high-order mesh. To improve the accuracy of approximating the target manifold with a noninterpolating curved high-order mesh, we minimize the square of the disparity measure expressed both in terms of the nodal coordinates of the physical and parametric curved high-order meshes. The proposed objective functions are continuously differentiable and thus, we are able to use minimization algorithms that require the first or the second derivatives of the objective function. Finally, we present several examples that show that the proposed methodology generates high-order approximations of the target manifold with optimal convergence rates for the geometric accuracy even when non-uniform parameterizations of the manifolds are prescribed. Accordingly, we can generate coarse curved high-order meshes significantly more accurate than finer low-order meshes that feature the same resolution.Peer ReviewedPostprint (author's final draft