Hierarchical multigrid approaches for the finite cell method on uniform and multi-level hp-refined grids

Abstract

This contribution presents a hierarchical multigrid approach for the solution of large-scale finite cell problems on both uniform grids and multi-level hp-discretizations. The proposed scheme leverages the hierarchical nature of the basis functions utilized in the finite cell method and the multi-level hp-method, which is attributed to the use of high-order integrated Legendre basis functions and overlay meshes, to yield a simple and elegant multigrid scheme. This simplicity is reflected in the fact that all restriction and prolongation operators reduce to binary matrices that do not need to be explicitly constructed. The coarse spaces are constructed over the different polynomial orders and refinement levels of the immersed discretization. Elementwise and patchwise additive Schwarz smoothing techniques are used to mitigate the influence of the cut cells leading to convergence rates that are independent of the cut configuration, mesh size and in certain scenarios even the polynomial order. The multigrid approach is applied to second-order problems arising from the Poisson equation and linear elasticity. A series of numerical examples demonstrate the applicability of the scheme for solving large immersed systems with multiple millions and even billions of unknowns on massively parallel machines.

References in corpus (11)

• Condition number analysis and preconditioning of the finite cell method
• The aggregated unfitted finite element method for elliptic problems
• Preconditioning immersed isogeometric finite element methods with application to flow problems
• Robust and parallel scalable iterative solutions for large-scale finite cell analyses
• CutIGA with Basis Function Removal
• Robust and scalable domain decomposition solvers for unfitted finite element methods
• Distributed-memory parallelization of the aggregated unfitted finite element method
• A 3D benchmark problem for crack propagation in brittle fracture
• The aggregated unfitted finite element method on parallel tree-based adaptive meshes
• A p-multigrid method enhanced with an ILUT smoother and its comparison to h-multigrid methods within Isogeometric Analysis
• Parallelization of the multi-level hp-adaptive finite cell method

Cited by in corpus (4)

• Industrial scale large eddy simulations (LES) with adaptive octree meshes using immersogeometric analysis
• An EigenValue Stabilization Technique for Immersed Boundary Finite Element Methods in Explicit Dynamics
• Algorithms for Parallel Generic $hp$-adaptive Finite Element Software
• Efficient Numerical Integration for Finite Element Trunk Spaces in 2D and 3D using Machine Learning: A new Optimisation Paradigm to Construct Application-Specific Quadrature Rules