Hua Tong

Hua Tong, Yongjie Jessica Zhang

Conforming hex meshes are widely regarded as an effective computational domain for simulation because of their nice numerical properties, yet automatically decomposing a general 3D volume into a conforming hex mesh remains a formidable challenge. Among existing approaches, methods that construct an adaptive Cartesian grid and subsequently convert it into a conforming mesh stand out for their robustness. However, topological conversion schemes require strict compatibility conditions that inevitably increase element count. State-of-the-art 2-refinement octree methods employ weakly-balanced and generalized pairing conditions to yield low element counts, but suffer from critical limitations: primal cell information is lost after dualization, and resulting dual cells often exhibit non-planar quad faces. Alternatively, 3-refinement 27-tree methods directly generate conforming hex meshes through template-based replacement, producing higher-quality elements with planar faces, but previous techniques impose far stricter conditions, severely over-refining grids by factors of ten to one hundred. This article introduces a novel 3-refinement approach using a moderately-balanced condition, slightly stronger than weakly-balanced but substantially more relaxed than prior 3-refinement requirements. The key insight is that recursively applying local refinements can isolate and reduce complex configurations to simpler cases covered by a fundamental template set. Two open-sourced variants are provided: one optimized for speed, and another trading some computational cost for marginally reduced element counts. Compared to previous 3-refinement methods, they significantly reduce final hex element counts while preserving min SJ values and guaranteeing convex polyhedral cells; relative to 2-refinement state-of-the-art, they also achieve a lower Hausdorff ratio using slightly fewer elements.

Hua Tong, Kuanren Qian, Eni Halilaj et al.

High-quality mesh generation is the foundation of accurate finite element analysis. Due to the vast interior vertices search space and complex initial boundaries, mesh generation for complicated domains requires substantial manual processing and has long been considered the most challenging and time-consuming bottleneck of the entire modeling and analysis process. In this paper, we present a novel computational framework named ``SRL-assisted AFM" for meshing planar geometries by combining the advancing front method with neural networks that select reference vertices and update the front boundary using ``policy networks." These deep neural networks are trained using a unique pipeline that combines supervised learning with reinforcement learning to iteratively improve mesh quality. First, we generate different initial boundaries by randomly sampling points in a square domain and connecting them sequentially. These boundaries are used for obtaining input meshes and extracting training datasets in the supervised learning module. We then iteratively improve the reinforcement learning model performance with reward functions designed for special requirements, such as improving the mesh quality and controlling the number and distribution of extraordinary points. Our proposed supervised learning neural networks achieve an accuracy higher than 98% on predicting commercial software. The final reinforcement learning neural networks automatically generate high-quality quadrilateral meshes for complex planar domains with sharp features and boundary layers.

Yuxuan Yu, Jiashuo Liu, Hua Tong et al.

Polycube structures provide parametric domains for all-hexahedral (all-hex) mesh generation and analysis-suitable volumetric spline construction in isogeometric analysis (IGA). Recent learning-based polycube pipelines have improved automation, yet several challenges remain when handling complex CAD geometries. These challenges include the limited diversity of primitive geometries, restricted grid configurations, and the increasing cost of genus-guided context search during inference as both the primitive set and the grid size grow. In this paper, we present {Scalable DDPM-Polycube}, an extended diffusion-based polycube construction method that addresses these limitations. First, we expand the primitive set from two primitive geometries to three by introducing a blind-hole cube primitive, thereby improving the representation of local hole-like features that do not change the global genus. Second, we extend the grid configuration from the previous $2\times 1$ setting to an enlarged three-dimensional grid configuration, which increases representational capacity and reduces mapping distortion for complex geometries. Third, we develop a genus-guided context generation strategy together with a hierarchical verification procedure, enabling robust context generation in both user-guided and automated modes. Once a valid polycube structure is generated, it is used for parametric mapping, all-hex control mesh generation, and volumetric spline construction. Experimental results demonstrate that scalable DDPM-Polycube improves the generality, scalability, and automation of diffusion-based polycube generation, and supports hex mesh generation and volumetric spline construction for IGA applications on complex geometries.

Yuxuan Yu, Yuzhuo Fang, Hua Tong et al.

In this paper, we present a novel algorithm that integrates deep learning with the polycube method (DL-Polycube) to generate high-quality hexahedral (hex) meshes, which are then used to construct volumetric splines for isogeometric analysis. Our DL-Polycube algorithm begins by establishing a connection between surface triangular meshes and polycube structures. We employ deep neural network to classify surface triangular meshes into their corresponding polycube structures. Following this, we combine the acquired polycube structural information with unsupervised learning to perform surface segmentation of triangular meshes. This step addresses the issue of segmentation not corresponding to a polycube while reducing manual intervention. Quality hex meshes are then generated from the polycube structures, with employing octree subdivision, parametric mapping and quality improvement techniques. The incorporation of deep learning for creating polycube structures, combined with unsupervised learning for segmentation of surface triangular meshes, substantially accelerates hex mesh generation. Finally, truncated hierarchical B-splines are constructed on the generated hex meshes. We extract trivariate Bézier elements from these splines and apply them directly in isogeometric analysis. We offer several examples to demonstrate the robustness of our DL-Polycube algorithm.

Hua Tong

The quality of mesh generation has long been considered a vital aspect in providing engineers with reliable simulation results throughout the history of the Finite Element Method (FEM). The element extraction method, which is currently the most robust method, is used in business software. However, in order to speed up extraction, the approach is done by finding the next element that optimizes a target function, which can result in local mesh of bad quality after many time steps. We provide TreeMesh, a method that uses this method in conjunction with reinforcement learning (also possible with supervised learning) and a novel Monte-Carlo tree search (MCTS) (Coulom(2006), Kocsis and Szepesvári(2006), Browne et~al.(2012)). The algorithm is based on a previously proposed approach (Pan et~al.(2021)). After making many improvements on DRL (algorithm, state-action-reward setting) and adding a MCTS, it outperforms the former work on the same boundary. Furthermore, using tree search, our program reveals much preponderance on seed-density-changing boundaries, which is common on thin-film materials.