The (Isogeometric) Finite Cell Method -in which a domain is immersed in a structured background mesh -suffers from conditioning problems when cells with small volume fractions occur. In this contribution, we establish a rigorous scaling relation between the condition number of (I)FCM system matrices and the smallest cell volume fraction. Ill-conditioning stems either from basis functions being small on cells with small volume fractions, or from basis functions being nearly linearly dependent on such cells. Based on these two sources of ill-conditioning, an algebraic preconditioning technique is developed, which is referred to as Symmetric Incomplete Permuted Inverse Cholesky (SIPIC). A detailed numerical investigation of the effectivity of the SIPIC preconditioner in improving (I)FCM condition numbers and in improving the convergence speed and accuracy of iterative solvers is presented for the Poisson problem and for two-and three-dimensional problems in linear elasticity, in which Nitche's method is applied in either the normal or tangential direction. The accuracy of the preconditioned iterative solver enables mesh convergence studies of the finite cell method.

Immersed finite element methods generally suffer from conditioning problems when cut elements intersect the physical domain only on a small fraction of their volume. De Prenter et al. [Computer Methods in Applied Mechanics and Engineering, 316 (2017) pp. 297-327] present an analysis for symmetric positive definite (SPD) immersed problems, and for this class of problems an algebraic preconditioner is developed. In this contribution the conditioning analysis is extended to immersed finite element methods for systems that are not SPD and the preconditioning technique is generalized to a connectivity-based preconditioner inspired by Additive-Schwarz preconditioning. This Connectivity-based Additive-Schwarz (CbAS) preconditioner is applicable to problems that are not SPD and to mixed problems, such as the Stokes and Navier-Stokes equations. A detailed numerical investigation of the effectivity of the CbAS preconditioner to a range of flow problems is presented.

Ill-conditioning of the system matrix is a well-known complication in immersed finite element methods and trimmed isogeometric analysis. Elements with small intersections with the physical domain yield problematic eigenvalues in the system matrix, which generally degrades efficiency and robustness of iterative solvers. In this contribution we investigate the spectral properties of immersed finite element systems treated by Schwarz-type methods, to establish the suitability of these as smoothers in a multigrid method. Based on this investigation we develop a geometric multigrid preconditioner for immersed finite element methods, which provides mesh-independent and cut-element-independent convergence rates. This preconditioning technique is applicable to higher-order discretizations, and enables solving large-scale immersed systems in parallel, at a computational cost that scales linearly with the number of degrees of freedom. The performance of the preconditioner is demonstrated for conventional Lagrange basis functions and for isogeometric discretizations with both uniform B-splines and locally refined approximations based on truncated hierarchical B-splines.

The finite cell method is a highly flexible discretization technique for numerical analysis on domains with complex geometries. By using a non-boundary conforming computational domain that can be easily meshed, automatized computations on a wide range of geometrical models can be performed. Application of the finite cell method, and other immersed methods, to large real-life and industrial problems is often limited due to the conditioning problems associated with these methods. These conditioning problems have caused researchers to resort to direct solution methods, which significantly limit the maximum size of solvable systems. Iterative solvers are better suited for large-scale computations than their direct counterparts due to their lower memory requirements and suitability for parallel computing. These benefits can, however, only be exploited when systems are properly conditioned. In this contribution we present an Additive-Schwarz type preconditioner that enables efficient and parallel scalable iterative solutions of large-scale multi-level hp-refined finite cell analyses.

Nitsche's method is a popular approach to implement Dirichlet-type boundary conditions in situations where a strong imposition is either inconvenient or simply not feasible. The method is widely applied in the context of unfitted finite element methods. From the classical (symmetric) Nitsche's method it is well-known that the stabilization parameter in the method has to be chosen sufficiently large to obtain unique solvability of discrete systems. In this short note we discuss an often used strategy to set the stabilization parameter and describe a possible problem that can arise from this. We show that in specific situations error bounds can deteriorate and give examples of computations where Nitsche's method yields large and even diverging discretization errors.

No abstract

We propose an adaptive mesh refinement strategy for immersed isogeometric analysis, with application to steady heat conduction and viscous flow problems. The proposed strategy is based on residual-based error estimation, which has been tailored to the immersed setting by the incorporation of appropriately scaled stabilization and boundary terms. Element-wise error indicators are elaborated for the Laplace and Stokes problems, and a THB-spline-based local mesh refinement strategy is proposed. The error estimation and adaptivity procedure are applied to a series of benchmark problems, demonstrating the suitability of the technique for a range of smooth and non-smooth problems. The adaptivity strategy is also integrated into a scan-based analysis workflow, capable of generating error-controlled results from scan data without the need for extensive user interactions or interventions.

This review paper discusses the developments in immersed or unfitted finite element methods over the past decade. The main focus is the analysis and the treatment of the adverse effects of small cut elements. We distinguish between adverse effects regarding the stability and adverse effects regarding the conditioning of the system, and we present an overview of the developed remedies. In particular, we provide a detailed explanation of Schwarz preconditioning, element aggregation, and the ghost penalty formulation. Furthermore, we outline the methodologies developed for quadrature and weak enforcement of Dirichlet conditions, and we discuss open questions and future research directions.