Reliable Whisker Weaving via Curve Contraction

Whisker Weaving is an advancing front algorithm for all-hexahedral mesh generation. It uses global information derived from grouping the mesh dual into surfaces, the STC, to construct the connectivity of the mesh, then positions the nodes afterwards. Currently, we are able to reliably generate hexahedral meshes for complicated geometries and surface meshes. However, the surface mesh must be modified locally. Also, in large, highly-unstructured meshes, there are usually isolated regions where hex quality is poor. Reliability has been achieved by using new, provable curvecontraction algorithms to sequence the advancing front process. We have also demonstrated that sheet moving can remove certain types of invalid connectivity.     

Localized coarsening of conforming all-hexahedral meshes

Finite element mesh adaptation methods can be used to improve the efficiency and accuracy of solutions to computational modeling problems. In many applications involving hexahedral meshes, localized modifications which preserve a conforming all-hexahedral mesh are desired. Effective hexahedral refinement methods that satisfy these criteria have recently become available; however, due to hexahedral mesh topology constraints, little progress has been made in the area of hexahedral coarsening. This paper presents a new method to locally coarsen conforming all-hexahedral meshes. The method works on both structured and unstructured meshes and is not based on undoing previous refinement. Building upon recent developments in quadrilateral coarsening, the method utilizes hexahedral sheet and column operations, including pillowing, column collapsing, and sheet extraction. A general algorithm for automated coarsening is presented and examples of models that have been coarsened with this new algorithm are shown. While results are promising, further work is needed to improve the automated process.     

The receding front method applied to hexahedral mesh generation of exterior domains

Two of the most successful methods to generate unstructured hexahedral meshes are the grid-based methods and the advancing front methods. On the one hand, the grid-based methods generate high-quality hexahedra in the inner part of the domain using an inside–outside approach. On the other hand, advancing front methods generate high-quality hexahedra near the boundary using an outside–inside approach. To combine the advantages of both methodologies, we extend the receding front method: an inside–outside mesh generation approach by means of a reversed advancing front. We apply this approach to generate unstructured hexahedral meshes of exterior domains. To reproduce the shape of the boundaries, we first pre-compute the mesh fronts by combining two solutions of the Eikonal equation on a tetrahedral reference mesh. Then, to generate high-quality elements, we expand the quadrilateral surface mesh of the inner body towards the unmeshed external boundary using the pre-computed fronts as a guide.     

A dual element based geometric element transformation method for all-hexahedral mesh smoothing

Due to their increased complexity hexahedral elements are more challenging with respect to mesh generation and mesh improvement techniques than tetrahedral elements. In particular, there is a lack of geometry-based all-hexahedral smoothing methods for mesh quality improvement being easy to implement, practicable, and efficient. The recently introduced geometric element transformation method represents a new promising element oriented smoothing concept to resolve this deficiency. By giving a dual octahedron based regularizing transformation this new approach is adapted in order to smooth all-hexahedral meshes. First numerical tests indicate that the resulting smoothing method yields high quality results at least comparable to those of a state of the art global optimization-based approach while being significantly faster.     

Adaptive generation of hexahedral element mesh using an improved grid-based method

An improved grid-based algorithm for the adaptive generation of hexahedral finite element mesh is presented in this paper. It is named as the inside-out grid-based method and involves the following four steps. The first step is the generation of an initial grid structure which envelopes the analyzed solid model completely. And the elements size and density maps are constructed based on the surface curvature and local thickness of the solid model. Secondly, the core mesh is generated through removing all the undesired elements using even and odd parity rules. The third step is to magnify the core mesh in an inside-out manner through a surface node projection process using the closest position approach. To match the mesh to the characteristic boundary of the solid model, a minimal Scaled Jacobian criterion is employed. Finally, in order to handle the degenerated elements and improve the quality of the resulting mesh, two comprehensive techniques are employed: the insertion technique and collapsing technique. The present method was applied in the mesh construction of different engineering problems. Scaled Jacobian and Skew metrics are used to evaluate the hexahedral element mesh quality. The application results show that all-hexahedral element meshes which are well-shaped and capture all the geometric features of the original solid models can be generated using the inside-out grid-based method presented in this paper.     

Anisotropic quadrilateral mesh generation: An indirect approach

In this paper a new indirect approach is presented for anisotropic quadrilateral mesh generation based on discrete surfaces. The ability to generate grids automatically had a pervasive influence on many application areas in particularly in the field of Computational Fluid Dynamics. In spite of considerable advances in automatic grid generation there is still potential for better performance and higher element quality. The aim is to generate meshes with less elements which fit some anisotropy criterion to satisfy numerical accuracy while reducing processing times remarkably. The generation of high quality volume meshes using an advancing front algorithm relies heavily on a well designed surface mesh. For this reason this paper presents a new technique for the generation of high quality surface meshes containing a significantly reduced number of elements. This is achieved by creating quadrilateral meshes that include anisotropic elements along a source of anisotropy.     

Automatic unstructured all-hexahedral mesh generation from B-Reps for non-manifold CAD assemblies

This paper describes an automatic and robust approach to convert non-manifold CAD assemblies into unstructured all-hexahedral meshes conformal to the given B-Reps (boundary-representations) and with sharp feature preservation. In previous works, we developed an octree-based isocontouring method to construct unstructured hexahedral meshes for arbitrary non-manifold and manifold domains. However, sharp feature preservation still remains a challenge, especially for non-manifold CAD assemblies. In this paper, boundary features such as NURBS (non-uniform rational B-Splines) curves and surface patches are first extracted from the given B-Reps. Features shared by multiple components are identified and distinguished. To preserve these non-manifold features, one given surface patch may need to be split into several small ones. An octree-based algorithm is then carried out to create an unstructured all-hexahedral base mesh, detecting and preserving all the sharp features via a curve and surface parametrization. Two sets of local refinement templates are provided for adaptive mesh generation, along with a novel 2-refinement implementation. Vertices in the base mesh are categorized into four groups based on the given non-manifold topology, and each group is relocated using various methods with all sharp features preserved. After this stage, a novel two-step pillowing technique is developed for such complicated non-manifold domains to eliminate triangle-shaped quadrilateral elements along the curves and “doublets”, handling non-manifold and manifold features in different ways. Finally, a combination of smoothing and optimization is used to further improve the mesh quality. Our algorithm is automatic and robust for non-manifold and manifold domains. We have applied our algorithm to several complicated CAD assemblies.     

Facial hexahedral mesh transferring by volumetric mapping based on harmonic fields

Hexahedral mesh has obvious mechanical advantages over tetrahedral mesh, but it is no trivial task to generate hexahedral mesh for complex object shapes such as individual faces. This paper presents a novel method to generate patient-specific hexahedral meshes of facial soft tissue models, based on a volumetric cross-parameterization mapping from a standard hexahedral mesh to the individual model. The volumetric parameterization is constructed based on triple of the volumetric harmonic fields, which are adapted to be as close to mutually orthogonal as possible, to achieve some quasi-conformal effect. In addition, some piecewise constraints on the harmonic fields are added to ensure anatomical feature correspondence. Experimental results show that our approach works efficiently for facial soft tissue modeling, avoids element flipping and preserves mesh element angles to a significant extent.     

A methodology for quadrilateral finite element mesh coarsening

High fidelity finite element modeling of continuum mechanics problems often requires using all quadrilateral or all hexahedral meshes. The efficiency of such models is often dependent upon the ability to adapt a mesh to the physics of the phenomena. Adapting a mesh requires the ability to both refine and/or coarsen the mesh. The algorithms available to refine and coarsen triangular and tetrahedral meshes are very robust and efficient. However, the ability to locally and conformally refine or coarsen all quadrilateral and all hexahedral meshes presents many difficulties. Some research has been done on localized conformal refinement of quadrilateral and hexahedral meshes. However, little work has been done on localized conformal coarsening of quadrilateral and hexahedral meshes. A general method which provides both localized conformal coarsening and refinement for quadrilateral meshes is presented in this paper. This method is based on restructuring the mesh with simplex manipulations to the dual of the mesh. In addition, this method appears to be extensible to hexahedral meshes in three dimensions. Sandia National Laboratories is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, for the United States Department of Energy under Contract DE-AC04-94AL85000.     

A control approach for pore size distribution in the bone scaffold based on the hexahedral mesh refinement

Tissue engineering is the application of that knowledge to the building or repairing of tissues. Generally, engineered tissue is a combination of living cells and a support structure called scaffolds. Modeling, design and fabrication of tissue scaffold with intricate architecture, porosity and pore size for desired tissue properties presents a challenge in tissue engineering. In this paper, a control approach for pore size distribution in the bone scaffold based on the hexahedral mesh refinement is presented. Firstly, the bone scaffold modeling approach based on the shape function in the finite element method is provided. The resulting various macroporous morphologies can be obtained. Then conformal refinement algorithm for all-hexahedral element mesh is illustrated. Finally, a modeling approach for constructing tissue engineering (TE) bone scaffold with defined pore size distribution is presented. Before the conformal refinement of all-hexahedral element mesh, a 3D mesh with various hexahedral elements must be provided. If all the pores in the bone scaffold need to be reduced, that means that the whole hexahedral mesh needs to be refined. Then the solid entity can be re-divided with altered subdivision parameters. If the pores in the local regions of bone need to be reduced, that means that 3D hexahedral mesh in the local regions needs to be refined. Based on SEM images, the pore size distribution in the normal bone can be obtained. Then, according to the conformal refinement of all-hexahedral element meshes, defined hexahedral size distribution can be gained, which leads to generate defined pore size distribution in the bone scaffold, for the pore morphology and size are controlled by various subdivided hexahedral elements. Compared to other methods such as varying processing parameters in supercritical fluid processing and multi-interior architecture design, the method proposed in this paper enjoys easy-controllability and higher accuracy.     

An interior surface generation method for all-hexahedral meshing

This paper describes an all-hexahedral generation method focusing on how to create interior surfaces. It is well known that a solid homeomorphic to a ball with even number of bounding quadrilaterals can be partitioned into a compatible hexahedral mesh where each associated hexahedron corresponds to the intersection of three interior surfaces that are dual to the original hexahedral mesh. However, no such method for creating dual interior surfaces has been developed for generating all-hexahedral meshes of volumes covered with simply connected quadrilaterals. We generate an interior surface as an orientable regular homotopy (or more definitively a sweep) by splitting a dual cycle into several pieces at self-intersecting points and joining the three connected pieces, if the self-intersecting point-types are identical, while we generate a non-orientable surface (containing Möbius bands) if the self-intersecting point-types are distinct. Stitching these simple interior surfaces together allows us to compose more complex interior surfaces. Thus, we propose a generalized method of generating a hexahedral mesh topology by directly creating the interior surface arrangement. We apply the present framework to Schneiders’ open pyramid problem and show an arrangement of interior surfaces that decompose Schneiders’ pyramid into 146 hexahedra.     

A Case Study in Hexahedral Mesh Generation: Simulation of the Human Mandible

We provide a case study for the generation of pure hexahedral meshes for the numerical simulation of physiological stress scenarios of the human mandible. Du to its complex and very detailed free-form geometry, the mandible model is very demanding. This test case is used as a running example to demonstrate the applicability of a combinatorial approach for the generation of hexahedral meshes by means of successive dual cycle eliminations, which has been proposed by the second author in previous work. We report on the progress and recent advances of the cycle elimination scheme. The given input data, a surface triangulation obtained from computed tomography data, requires a substantial mesh reduction and a suitable conversion into a quadrilateral surface mesh as a first step, for which we use mesh clustering and b-matching techniques. Several strategies for improved cycle elimination orders are proposed. They lead to a significant reduction in the mesh size and a better structural quality. Based on the resulting combinatorial meshes, gradient-based optimized smoothing with the condition number of the Jacobian matrix as objective together with mesh untangling techniques yielded embeddings of a satisfactory quality. To test our hexahedral meshes for the mandible model within an FEM simulation we used the scenario of a bite on a ‘hard nut.’ Our simulation results are in good agreement with observations from biomechanical experiments.     

Hexahedral mesh smoothing via local element regularization and global mesh optimization

The quality of finite element meshes is one of the key factors that affect the accuracy and reliability of finite element analysis results. In order to improve the quality of hexahedral meshes, we present a novel hexahedral mesh smoothing algorithm which combines a local regularization for each hexahedral mesh, using dual element based geometric transformation, with a global optimization operator for all hexahedral meshes. The global optimization operator is composed of three main terms, including the volumetric Laplacian operator of hexahedral meshes and the geometric constraints of surface meshes which keep the volumetric details and the surface details, and another is the transformed node displacements condition which maintains the regularity of all elements. The global optimization operator is formulated as a quadratic optimization problem, which is easily solved by solving a sparse linear system. Several experimental results are presented to demonstrate that our method obtains higher quality results than other state-of-the-art approaches.     

扫掠法有限元网格生成方法

为了提高有限元网格的生成质量,扫掠法生成六面体网格过程中内部节点定位成为关键一步,在研究复杂扫掠体六面体有限元网格生成算法过程中,提出了一种基于扫掠法的六面体网格生成算法,算法利用源曲面已经划分好的网格和连接曲面的结构化网格,用仿射映射逐层投影,生成目标曲面,提出基于Roca算法的内部节点定位的新算法,运用由外向内推进的波前法思想,生成全部的六面体网格。通过实例表明,该算法快速,稳定,可靠,可处理大量复杂2.5维实体六面体网格生成问题。     

Simultaneous untangling and smoothing of quadrilateral and hexahedral meshes using an object-oriented framework

In this work, we present a simultaneous untangling and smoothing technique for quadrilateral and hexahedral meshes. The algorithm iteratively improves a quadrilateral or hexahedral mesh by minimizing an objective function defined in terms of a regularized algebraic distortion measure of the elements. We propose several techniques to improve the robustness and the computational efficiency of the optimization algorithm. In addition, we have adopted an object-oriented paradigm to create a common framework to smooth meshes composed by any type of elements, and using different minimization techniques. Finally, we present several examples to show that the proposed technique obtains valid meshes composed by high-quality quadrilaterals and hexahedra, even when the initial meshes contain a large number of tangled elements.     

TUM.GeoFrame: automated high-order hexahedral mesh generation for shell-like structures

This paper presents a fully automated high-order hexahedral mesh generation algorithm for shell-like structures based on enhanced sweeping methods. Traditional sweeping techniques create all-hexahedral element meshes for solid structures by projecting an initial single surface mesh along a specified trajectory to a specified target surface. The work reported here enhances the traditional method for thin solids by creating conforming high-order all-hexahedral finite element meshes on an enhanced surface model with surfaces intersecting in parallel, perpendicular and skew-angled directions. The new algorithm is based on cheap projection rules separating the original surface model into a set of disjoint single surfaces and a so-called interface skeleton. The core of this process is reshaping the boundary representations of the initial surfaces, generating new sweeping templates along the intersection curves and joining the single swept hex meshes in an independently generated interface mesh.     

A two-dimensional paving mesh generator for triangles with controllable aspect ratio and quadrilaterals with high quality

Certain classes of problems result in solution fields of which the characteristic length scales vary with the orientation. Often the orientation of these length scales is related to the orientation of the boundaries. Such solution fields can be captured by the finite element method, using a mesh that is refined in the direction of the short length scales and coarse in the other directions. These meshes contain elements with high aspect ratios in a predefined pattern.The mesh generator presented here can render triangles with high aspect ratios through a paving algorithm. The paving algorithm that is employed applies both triangles and quadrilaterals, combining the advantages of both to render a qualitatively good, oriented triangular mesh, with a concentration of elements in the direction where the internal length scales of the solution field are the shortest.The mesh generator produces triangles with one (almost) orthogonal corner. When low aspect ratio triangles are generated, these are well suited for conversion to quadrilateral elements. Test results indicate that quadrilateral meshes converted from the mesh generator introduced here have a considerably better quality than those converted from several other triangular mesh generators.     

四边形网格生成中的前沿边生长改进算法

为提高B样条曲面重构中点云四边形网格的生成效率和质量,对现有的四边形网格Q-Morph前沿边生长算法进行改进,提出面向四边形网格生成的三角网格拓扑优化方法,通过设定生长限制条件和调整网格顶点度,保证全局四边形网格质量,实现适合复杂曲面重构的规则四边形网格获取。实例结果表明,该算法效率高、适应性好,生成的四边形网格具有分布均匀、不规则网格数量少的优点。     

A strategy of automatic hexahedral mesh generation by using an improved whisker-weaving method with a surface mesh modification procedure

One of the demands for three dimensional (3D) finite element analyses is the development of an automatic hexahedral mesh generator. For this problem, several methods have been proposed by many researchers. However, reliable automatic hexahedral mesh generation has not been developed at present. In this paper, a new strategy of fully automatic hexahedral mesh generation is proposed. In this strategy, the prerequisite for generating a hexahedral mesh is a quadrilateral surface mesh. From the given surface mesh, combinatorial dual cycles (sheet loops for the whisker-weaving algorithm) are generated to produce a hexahedral mesh. Since generating a good quality hexahedral mesh does not depend only on the quality of quadrilaterals of the surface mesh but also on the quality of the sheet loops generated from it, a surface mesh modification method to remove self-intersections from sheet loops is developed. Next, an automatic hexahedral mesh generator by the improved whisker-weaving algorithm is developed in this paper. By creating elements and nodes on 3D real space during the weaving process, it becomes possible to generate a hexahedral mesh with fewer bad-quality elements. Several examples will be presented to show the validity of the proposed mesh generation strategy.     

Hexahedral Mesh Generation by Successive Dual Cycle Elimination

We propose a new method for constructing all-hexahedral finite element meshes. The core of our method is to build up a compatible combinatorial cell complex of hexahedra for a solid body which is topologically a ball, and for which a quadrilateral surface mesh of a certain structure is prescribed. The step-wise creation of the hex complex is guided by the cycle structure of the combinatorial dual of the surface mesh. Our method transforms the graph of the surface mesh iteratively by changing the dual cycle structure until we get the surface mesh of a single hexahedron. Starting with a single hexahedron and reversing the order of the graph transformations, each transformation step can be interpreted as adding one or more hexahedra to the so far created hex complex. Given an arbitrary solid body, we first decompose it into simpler subdomains equivalent to topological balls by adding virtual 2-manifolds. Secondly, we determine a compatible quadrilateral surface mesh for all subdomains created. Then, in the main part we can use the core routine to build up a hex complex for each subdomain independently. The embedding and smoothing of the combinatorial mesh(es) finishes the mesh generation process. First results obtained for complex geometries are encouraging.