Three-Dimensional Numerical Simulation of Drops Suspended in Poiseuille Flow: Effect of Reynolds Number

Authors: A. Nourbakhsh

Abstract:

A finite difference/front tracking method is used to study the motion of three-dimensional deformable drops suspended in plane Poiseuille flow at non-zero Reynolds numbers. A parallel version of the code was used to study the behavior of suspension on a reasonable grid resolution (grids). The viscosity and density of drops are assumed to be equal to that of the suspending medium. The effect of the Reynolds number is studied in detail. It is found that drops with small deformation behave like rigid particles and migrate to an equilibrium position about half way between the wall and the centerline (the Segre-Silberberg effect). However, for highly deformable drops there is a tendency for drops to migrate to the middle of the channel, and the maximum concentration occurs at the centerline. The effective viscosity of suspension and the fluctuation energy of the flow across the channel increases with the Reynolds number of the flow.

Keywords: Suspensions, Poiseuille flow, Effective viscosity, Reynolds number.

Digital Object Identifier (DOI): doi.org/10.5281/zenodo.1097088

References:

[1] G. Segre, and A. Silberberg, “Behaviour of macroscopic rigid spheres in poiseuille flow, part1. determination of local concentration by statistical analysis of particle passages through crossed light beams,” J. Fluid Mech. vol. 14, 1962a, pp. 115.
[2] G. Segre, and A. Silberberg, “Behaviour of macroscopic rigid spheres in poiseuille flow, part2. experimental results and interpretation,” J. Fluid Mech. vol. 14, 1962b, pp. 136.
[3] H. L. Goldsmith, and S. G. Mason, “The flow of suspensions through tubes. I. Single spheres, rods and discs,” J. Colloid Sci. vol. 17, 1962, pp. 448.
[4] T. A. Kowalewski, “Concentration and velocity measurement in the flow of droplet suspensions through a tube,” Experiments in Fluids. vol. 2, 1984, pp. 213.
[5] J. P. Matas, J. F. Morris, and E. Guazzelli, “Inertial migration of rigid spherical particles in Poiseuille flow,” J. Fluid Mech. vol. 515, 2004, pp. 171.
[6] H. Zhou, and C. Pozrikidis, “Pressure-driven flow of suspensions of liquid drops,” Phys. Fluids. vol. 6, 1994, pp. 80.
[7] P. Nott, and J. Brady, “Pressure-driven flow of suspensions: simulation and theory,” J. Fluid Mech. vol. 275, 1994, pp. 157.
[8] M. Loewenberg, and E. J. Hinch, “Numerical simulation of a concentrated emulsion in shear flow,” J. Fluid Mech. vol. 321, 1996, pp. 395.
[9] S. K. Doddi, and P. Bagchi, “Three-dimensional computational modeling of multiple deformable cells flowing in micro vessels,” Physical Review E. vol. 79, 2009, pp. 046318-1.
[10] X. Li, and C. Pozrikidis, “Wall-bounded shear flow and channel flow of suspensions of liquid drops,” Int. J. Multiphase Flow. vol. 26, 2000, pp. 1247.
[11] D. J. Feng, H. H. Hu, and D. D. Joseph, “Direct simulation of initial value problems for the motion of solid bodies in a Newtonian fluid. Part1. Sedimentation,” J. Fluid Mech. vol. 261, 1994a, pp. 95.
[12] D. J. Feng, H. H. Hu, and D. D. Joseph, “Direct simulation of initial value problems for the motion of solid bodies in a Newtonian fluid. Part2. Couette and Poiseuille flows,” J. Fluid Mech. vol. 277, 1994b, pp. 271.
[13] S. S. Mortazavi, and G. Tryggvason, “A numerical study of the motion of drops in poiseuille flow, part1: Lateral migration of one drop,” J. Fluid Mech. vol. 411, 2000, pp. 325.
[14] A. Nourbakhsh, S. S. Mortazavi, “A three-dimensional study of the motion of a drop in plane Poiseuille flow at finite Reynolds numbers,” Iranian J. Science and Technology, Transaction B: Engineering. vol. 34, No. B2, 2010, pp. 179.
[15] M. Bayareh, and S. S. Mortazavi, “Numerical simulation of the motion of a single drop in a shear flow at finite Reynolds numbers,” Iranian J. Science and Technology, Transaction B: Engineering. vol. 33, No. B5, 2009, pp. 441.
[16] S. S. Mortazavi, Y. Afshar, and H. Abbuspour, “Numerical simulation of two-dimensional drops suspended in simple shear flow at nonzero Reynolds numbers,” J. Fluids. Eng. vol. 133, 2011, pp. 031303-1.
[17] M. Bayareh, and S. S. Mortazavi, “Binary collision of drops in simple shear flow at finite Reynolds numbers: Geometry and viscosity ratio effects,” Advances in Eng. Software. vol. 42, 2011, pp. 604.
[18] J. Lovick, and P. Angeli, “Droplet size and velocity profiles in liquidliquid horizontal flows,” Chem. Eng. Sci. vol. 59, 2004, pp. 3105.
[19] J. Adams, “MUDPACK: Multigrid FORTRAN software for the efficient solution of linear elliptic partial differential equations,” Appl. Math. Comput. vol. 34, 1989, pp. 113.
[20] M. Gorokhovski, and M. Herrmann, “Modeling primary atomization,” Annual Review of Fluid Mech. vol. 40, 2008, pp. 343.
[21] S. O. Unverdi, and G. Tryggvason, “A front-tracking method for viscous incompressible multi-fluid flows,” J. Comput. Phys. vol. 100, 1962a, pp. 25.
[22] S. O. Unverdi, and G. Tryggvason, “Computations of multi-fluid flows.” Physics. vol. 60(D), 1962b, pp. 70.
[23] G. Tryggvason, B. Bunner, A. Esmaeeli, D. Juric, N. Al-Rawahi, W. Tauber, J. Han, S. Nas, and Y. J. Jan, “A front-tracking method for the Computations of Multiphase Flow,” J. Computational Physics. vol. 169, 2001, pp. 708.