A boundary integral formulation for elastically deformable particles in a viscous fluid

This paper reports a formulation and implementation of a mixed (both direct and indirect) boundary element method using the double layer and its adjoint in a form suitable for solving Stokes flow problems involving elastically deformable particles. The formulation is essentially the Completed Double Layer Boundary Element Method for solving an exterior traction problem for the surrounding fluid or solid phase, followed by an interior displacement, and a mobility problem (if required) for the elastic particles. At the heart of the method is a deflation procedure that allows iterative solution strategies to be adopted, effectively opens the way for large-scale simulations of suspensions of deformable particles to be performed. Several problems are considered, to illustrate and benchmark the method. In particular, an analytical solution for an elastic sphere in an elongational flow is derived. The stresslet calculations for an elastic sphere in shear and elongational flows indicate that elasticity of the inclusions can potentially lead to positive second normal stress difference in shear flow, and an increase in the tensile resistance in elongational flow.