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

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.

• A. Nourbakhsh

• Suspensions
• Poiseuille flow
• Effective viscosity
• Reynolds number.

[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.