A boundary integral method for computing forces on particles in unsteady Stokes and linear viscoelastic fluids. (15th January 2016)

• Record Type: Journal Article
• Authors: Feng, H.; Córdoba, A.; Hernandez, F.; Indei, T.; Li, S.; Li, X.; Schieber, J. D.
• Is Part Of: International journal for numerical methods in fluids. Volume 82:Number 4(2016)
• Page Start: 198
• Page End: 217
• Publication Date: 2016-01-15
• Subjects: particulate flow -- linear viscoelastic fluid -- unsteady stokes -- boundary integral method; Fluid dynamics -- Mathematics -- Periodicals; 532
• DOI: 10.1002/fld.4216
• Languages: English
• ISSNs: 0271-2091

Abstract

Summary: Accurately characterizing the forces acting on particles in fluids is of fundamental importance for understanding particle dynamics and binding kinetics. Conventional asymptotic solutions may lead to poor accuracy for neighboring particles. In this paper, we develop an accurate boundary integral method to calculate forces exerted on particles for a given velocity field. We focus our study on the fundamental two‐bead oscillating problem in an axisymmetric frame. The idea is to exploit a correspondence principle between the unsteady Stokes and linear viscoelasticity in the Fourier domain such that a unifying boundary integral formulation can be established for the resulting Brinkman equation. In addition to the dimension reduction vested in a boundary integral method, our formulation only requires the evaluation of single‐layer integrals, which can be carried out efficiently and accurately by a hybrid numerical integration scheme based on kernel decompositions. Comparison with known analytic solutions and existing asymptotic solutions confirms the uniform third‐order accuracy in space of our numerical scheme. Copyright © 2016 John Wiley & Sons, Ltd. Abstract : We develop an accurate boundary integral method to calculate forces exerted on particles in unsteady Stokes flow and linear viscoelastic fluids. Our numerical method is third‐order accurate uniformly in space and corrects the error due to the poles at the axis of symmetry.