All-Hex Meshing Strategies For Densely Packed Spheres

Yu-Hsiang Lan,Paul Fischer,Elia Merzari,Misun Min

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

2021-06-02

Abstract:We develop an all-hex meshing strategy for the interstitial space in beds of densely packed spheres that is tailored to turbulent flow simulations based on the spectral element method (SEM). The SEM achieves resolution through elevated polynomial order N and requires two to three orders of magnitude fewer elements than standard finite element approaches do. These reduced element counts place stringent requirements on mesh quality and conformity. Our meshing algorithm is based on a Voronoi decomposition of the sphere centers. Facets of the Voronoi cells are tessellated into quads that are swept to the sphere surface to generate a high-quality base mesh. Refinements to the algorithm include edge collapse to remove slivers, node insertion to balance resolution, localized refinement in the radial direction about each sphere, and mesh optimization. We demonstrate geometries with 10^2-10^5 spheres using approximately 300 elements per sphere (for three radial layers), along with mesh quality metrics, timings, flow simulations, and solver performance.

Computational Engineering, Finance, and Science

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to generate high - quality all - hexahedral meshes in a closely - packed sphere bed for turbulence simulations based on the Spectral Element Method (SEM). Specifically, the paper focuses on how to create an all - hexahedral mesh with strong adaptability and a small number of elements in the void space between spheres, which is crucial for achieving high - resolution turbulence simulations. Since SEM achieves resolution by increasing the polynomial order \(N\) and requires two to three orders of magnitude fewer elements than the standard finite - element method, strict requirements are imposed on mesh quality and consistency. The paper proposes a mesh - generation algorithm based on Voronoi decomposition centered on the sphere center. High - quality base meshes are generated through tessellation, quadrangulation of Voronoi cells and sweeping towards the sphere surface, and a series of improvements including edge collapse, node insertion, local refinement and mesh optimization are carried out.