What is the best way to solve elasticity problems where the geometry is variable? Explicit meshing can be costly and cumbersome. The diffuse boundary method is an attractive alternative that replaces explicit boundaries with source terms, and allows for a geometry-agnostic meshing scheme with regularized boundaries. However, this can result in solver issues because the problem becomes singular in the region that does not contain material. In this work we show that it is possible to efficiently solve these kinds of near-singular problems using BSAMR and maintain excellent solver convergence.
Full article: https://doi.org/10.1007/s00466-023-02325-8
Abstract: Many solid mechanics problems on complex geometries are conventionally solved using discrete boundary methods. However, such an approach can be cumbersome for problems involving evolving domain boundaries due to the need to track boundaries and constant remeshing. The purpose of this work is to present a comprehensive strategy for efficiently solving such problems on an adaptive structured grid, while expositing some of the basic yet important nuances associated with solving near-singular problems in strong form. We employ a robust smooth boundary method (SBM) that represents complex geometry implicitly, in a larger and simpler computational domain, as the support of a smooth indicator function. We present the resulting semidefinite equations for mechanical equilibrium, in which inhomogeneous boundary conditions are replaced by source terms. In this work, we present a computational strategy for efficiently solving near-singular SBM-based solid mechanics problems. We use the block-structured adaptive mesh refinement method, coupled with a geometric multigrid solver for an efficient solution of mechanical equilibrium. We discuss some of the practical numerical strategies for implementing this method, notably including the importance of grid versus node-centered fields. We demonstrate the solver’s accuracy and performance for three representative examples: (a) plastic strain evolution around a void, (b) crack nucleation and propagation in brittle materials, and (c) structural topology optimization. In each case, we show that very good convergence of the solver is achieved, even with large near-singular areas, and that any convergence issues arise from other complexities, such as stress concentrations.