ECCOMAS 2024

Large-Scale Parallel Mesh Adaptation for High-Order Discontinuous Galerkin Methods Applied to Scale-Resolving Simulations

  • Coulaud, Olivier (Cenaero)
  • Rocca, Andrea (Cenaero)
  • Toulorge, Thomas (Cenaero)

Please login to view abstract download link

Large-Eddy Simulation (LES) and Direct Numerical Simulation (DNS) are becoming applicable to problems of increasing industrial relevance by virtue of two main recent advances, namely i) the availability and efficient exploitation of very large-scale computing resources, and ii) the maturation of high-order numerical methods capable of dealing with complex geometries without accuracy degradation (particularly on unstructured meshes). Despite this progress, computational cost remains a challenge for such scale-resolving simulation methods to be applied to flows of higher Reynolds number. Moreover, the pre-processing phase still requires significant efforts and specific skills from experienced practitioners, as the cost and turnaround time of such computations does not favour manual mesh improvement by the user. Automatic mesh adaptation techniques have the potential to alleviate both of these problems. In this work, we integrate a metric-based anisotropic mesh adaptation procedure in the workflow of massively parallel scale-resolving simulations based on a high-order Discontinuous Galerkin spatial scheme [1]. The adaptation is driven by an error estimator that relies on the log-simplex method to extend Hessian-based methods to high-order Finite Element solutions [2]. In the case of LES, the estimator is supplemented with physics-based adaptation criteria. Most importantly, the adaptation procedure is parallelized for distributed-memory systems and includes a load balancing capability which provides very good scalability. In addition to academic test cases demonstrating the performance of the mesh adaptation method, we will show its benefits on a practical turbomachinery application, namely the LES of flow around a low-pressure compressor blade in cascade configuration. [1] C. Carton de Wiart, K. Hillewaert. Development and validation of a massively parallel high-order solver for DNS and LES of industrial flows. In IDIHOM: Industrialization of High-Order Methods-A Top-Down Approach, pp. 251–292, Springer, 2015. [2] O. Coulaud, P. Schrooyen, A. Loseille. Anisotropic mesh adaptation for high-order Finite Element spaces with the log-simplex method. Application to Discontinuous Galerkin methods. Submitted to J. Comput. Phys., 2023.