Numerical simulation of rigid particles in Stokes flow: lubrication correction for general shapes of particles

Aline Lefebvre-Lepot¹ ; Flore Nabet¹

¹ CMAP, CNRS, Ecole Polytechnique, Institut Polytechnique de Paris, 91128 Palaiseau, France.

Mathematical modelling of natural phenomena, Tome 16 (2021), article no. 45

Résumé

We address the problem of numerical simulation of suspensions of rigid particles in a Stokes flow. We focus on the inclusion of the singular short range interaction effects (lubrication effects) in the simulations when the particles come close one to another. The problem is solved without introducing new hypothesis nor model. As in Lefebvre-Lepot et al. [J. Fluid Mech. 769 (2015) 369–386], the key idea is to decompose the velocity and pressure flows in a sum of a singular and a regular part. In this article, the singular part is computed using an explicit asymptotic expansion of the solution when the distance goes to zero. This expansion is similar to the asymptotic expansion proposed in Hillairet and Kelai [Asymptotic Anal. 95 (2015) 187–241] but is more appropriate for numerical simulations of suspensions. It can be computed for any locally convex (that is the particles have to be convex close to the contact point) and regular shape of particles. Using Hillairet and Kelai [Asymptotic Anal. 95 (2015) 187–241] as an intermediate result, we prove that the remaining part is regular in the sense that it is bounded independently of the distance. As a consequence, only a small number of degrees of freedom are necessary to obtain accurate results. The method is tested in dimension 2 for clusters of two or three aligned particles with general rigid velocities. We show that, as expected, the convergence is independent of the distance.

Détail

DOI : 10.1051/mmnp/2021037

Bibliographie

[1] B. Cichocki, B.U. Felderhof, K. Hinsen, E. Wajnryb, J. Bławzdziewicz Friction and mobility of many spheres in stokes flow J. Chem. Phys 1994 3780 3790

[2] R.G. Cox The motion of suspended particles almost in contact Int. J. Multiphase Flow 1974 343 371

[3] L. Durlofsky, J.F. Brady, G. Bossis Dynamic simulation of hydrodynamically interacting particles J. Fluid Mech 1987 21 49

[4] FreeFem webpage, https://freefem.org/. Accessed: 2019-09-30.

[5] S. Gallier, E. Lemaire, L. Lobry, P. François A fictitious domain approach for the simulation of dense suspensions J. Comput. Phys 2014 367 387

[6] B.D. Goddard, R.D. Mills-Williams, J. Sun The singular hydrodynamic interactions between two spheres in stokes flow Phys. Fluids 2020 062001

[7] M Hillairet, T Kelai Justification of lubrication approximation: an application to fluid/solid interactions Asymptotic Anal 2015 187 241

[8] A.J.C. Ladd Numerical simulations of particulate suspensions via a discretized Boltzmann equation. Part 1. Theoretical foundation J. Fluid Mech. 1994 285 309

[9] A.J.C. Ladd Numerical simulations of particulate suspensions via a discretized Boltzmann equation. Part 2. Numerical results J. Fluid Mech. 1994 311 339

[10] B. Lambert, L. Weynans, M. Bergmann Local lubrication model for spherical particles within incompressible Navier-Stokes flows Phys. Rev. E 2018 033313

[11] A. Lefebvre-Lepot Numerical simulation of suspensions: lubrication correction, including fluid correction Equations aux dérivées partielles et leurs applications Actes du colloque Edp-Normandie. Caen 2017 2017 87 100

[12] A. Lefebvre-Lepot, B. Merlet, T.N. Nguyen An accurate method to include lubrication forces in numerical simulations of dense Stokesian suspensions J. Fluid Mech 2015 369 386

[13] A.D. Maude End effects in a falling-sphere viscometer Br. J. Appl. Phys 1961 293

[14] N.A. Patankar, P. Singh, D.D. Joseph, R. Glowinski, T.-W. Pan A new formulation of the distributed Lagrange multiplier/fictitious domain method for particulate flows Int. J. Multiphase Flow 2000 1509 1524

[15] A.S. Sangani, G. Mo Inclusion of lubrication forces in dynamic simulations Phys. Fluids 1994 1653 1662

[16] K. Yeo, M.R. Maxey Simulation of concentrated suspensions using the force-coupling method J. Comput. Phys 2010 2401 2421