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.

No abstract

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.

Article available under the terms of the CC-BY-NC-ND licence (https://creativecommons.org/licenses/by-nc-nd/4.0/) eprints@whiterose.ac.uk https://eprints.whiterose.ac.uk/ Reuse Unless indicated otherwise, fulltext items are protected by copyright with all rights reserved. The copyright exception in section 29 of the Copyright, Designs and Patents Act 1988 allows the making of a single copy solely for the purpose of non-commercial research or private study within the limits of fair dealing. The publisher or other rights-holder may allow further reproduction and re-use of this version -refer to the White Rose Research Online record for this item. Where records identify the publisher as the copyright holder, users can verify any specific terms of use on the publisher's website. Takedown AbstractIn the phase-field description of brittle fracture, the fracture-surface area can be expressed as a functional of the phase field (or damage field). In this work we study the applicability of this explicit expression as a (non-linear) pathfollowing constraint to robustly track the equilibrium path in quasi-static fracture propagation simulations, which can include snap-back phenomena. Moreover, we derive a fracture-controlled staggered solution procedure by systematic decoupling of the path-following controlled elasticity and phasefield problems. The fracture-controlled monolithic and staggered solution procedures are studied for a series of numerical test cases. The numerical results demonstrate the robustness of the new approach, and provide insight in the advantages and disadvantages of the monolithic and staggered procedures.

We present an adaptive isogeometric-analysis approach to elasto-capillary fluid-solid interaction (FSI), based on a diffuse-interface model for the binary fluid and an Arbitrary-Lagrangian-Eulerian formulation for the FSI problem.We consider approximations constructed from adaptive high-regularity truncated hierarchical splines, as employed in the isogeometric analysis (IGA) paradigm. The considered adaptive strategy comprises a two-level hierarchical a posteriori error estimate. The hierarchical a posteriori error estimate directs a support-based refinement procedure. To attain robustness of the solution procedure for the aggregated binary-fluid-solid-interaction problem, we apply a fully monolithic solution procedure and we introduce a continuation process in which the diffuse interface of the binary fluid is artificially enlarged on the coarsest levels of the adaptive-refinement procedure. To assess the capability of the presented adaptive IGA method for elasto-capillary FSI, we conduct numerical computations for a configuration pertaining to a sessile droplet on a soft solid substrate.