Proteus: A direct forcing method in the simulations of particulate flows. (English) Zbl 1076.76568

Feng, Zhi-Gang; Michaelides, Efstathios E.

J. Comput. Phys. 202, No. 1, 20-51 (2005).

Summary: A new and efficient direct numerical method for the simulation of particulate flows is introduced. The method combines desired elements of the immersed boundary method, the direct forcing method and the lattice Boltzmann method. Adding a forcing term in the momentum equation enforces the no-slip condition on the boundary of a moving particle. By applying the direct forcing scheme, Proteus eliminates the need for the determination of free parameters, such as the stiffness coefficient in the penalty scheme or the two relaxation parameters in the adaptive-forcing scheme. The method presents a significant improvement over the previously introduced immersed-boundary-lattice-Boltzmann method (IB-LBM) where the forcing term was computed using a penalty method and a user-defined parameter. The method allows the enforcement of the rigid body motion of a particle in a more efficient way. Compared to the ”bounce-back” scheme used in the conventional LBM, the direct-forcing method provides a smoother computational boundary for particles and is capable of achieving results at higher Reynolds number flows. By using a set of Lagrangian points to track the boundary of a particle, Proteus eliminates any need for the determination of the boundary nodes that are prescribed by the ”bounce-back” scheme at every time step. It also makes computations for particles of irregular shapes simpler and more efficient. Proteus has been developed in two- as well as three-dimensions. This new method has been validated by comparing its results with those from experimental measurements for a single sphere settling in an enclosure under gravity. As a demonstration of the efficiency and capabilities of the present method, the settling of a large number (1232) of spherical particles is simulated in a narrow box under two different boundary conditions. It is found that when the no-slip boundary condition is imposed at the front and rear sides of the box the particles motion is significantly hindered. Under the periodic boundary conditions, the particles move faster. The simulations show that the sedimentation characteristics in a box with periodic boundary conditions at the two sides are very close to those found in the sedimentation of two-dimensional circular particles.

MSC:

76M25	Other numerical methods (fluid mechanics) (MSC2010)
76M28	Particle methods and lattice-gas methods
76T20	Suspensions

Keywords:

immersed boundary method; lattice Boltzmann method; spherical particles; sedimentation

Software:

Proteus

References:

[1]	Abraham, F., Functional dependence of drag coefficient of a sphere on Reynolds number, Phys. fluids, 13, 2194, (1970)
[2]	Aidun, C.K.; Lu, Y.-N.; Ding, E.-J., Direct analysis of particulate suspensions with inertia using the discrete Boltzmann equation, J. fluid mech., 373, 287, (1998) · Zbl 0933.76092
[3]	T. Bohlin, Terminal velocities of solid spheres in cylindrical enclosures, Transactions of the Royal Institute of Technology, Stockholm, 1960, Report# 155
[4]	Brady, J.F.; Bossis, G., Stokesian dynamics, Ann. rev. fluid. mech., 20, 111-157, (1988)
[5]	Chen, S.; Doolen, G.D., Lattice Boltzmann method for fluid flows, Annu. rev. fluid mech., 30, 329-364, (1998) · Zbl 1398.76180
[6]	Faxen, H., Der widerstand gegen die bewegung einer starren kugel in einer zum den flussigkeit, die zwischen zwei parallelen ebenen winden eingeschlossen ist, Annalen der physik, 68, 89-119, (1922) · JFM 48.1319.03
[7]	Feng, J.; Hu, H.H.; Joseph, D.D., Direct simulation of initial value problems for the motion of solid bodies in a Newtonian fluid. part 1: sedimentation, J. fluid mech., 261, 95-134, (1994) · Zbl 0800.76114
[8]	Feng, J.; Hu, H.H.; Joseph, D.D., Direct simulation of initial value problems for the motion of solid bodies in a Newtonian fluid. part 2: Couette and Poiseuille flows, J. fluid mech., 277, 271-301, (1994) · Zbl 0876.76040
[9]	Feng, Z.-G.; Michaelides, E.E., An immersed boundary method combined with lattice Boltzmann method for solving fluid and particles interaction problems, J. comput. phys., 195, 602-628, (2004) · Zbl 1115.76395
[10]	Feng, Z.G.; Michaelides, E.E., Hydrodynamic force on spheres in cylindrical and prismatic enclosures, Int. J. multiphase flow, 28, 479-496, (2002) · Zbl 1136.76503
[11]	Feng, Z.G.; Michaelides, E.E., Interparticle forces and lift on a particle attached to a solid boundary in suspension flow, Phys. fluids, 14, 49-60, (2002) · Zbl 1184.76158
[12]	Fogelson, A.L.; Peskin, C.S., A fast numerical method for solving the three-dimensional Stokes equation in the presence of suspended particles, J. comput. phys., 79, 50-69, (1988) · Zbl 0652.76025
[13]	Glowinski, R.; Pan, T.-W.; Hesla, T.I.; Joseph, D.D., A distributed Lagrange multiplier/fictitious domain method for particulate flows, Int. J. multiphase flow, 25, 755-794, (1999) · Zbl 1137.76592
[14]	Glowinski, R.; Pan, T.-W.; Hesla, T.I.; Joseph, D.D.; Periaux, J., A fictitious domain approach to the direct numerical simulation of incompressible viscous flow past moving rigid bodies: application to particulate flow, J. comput. phys., 169, 363-426, (2001) · Zbl 1047.76097
[15]	Goldstein, D.; Handler, R.; Sirovich, L., Modeling a no-slip flow boundary with an external force field, J. comput. phys., 105, 354-366, (1993) · Zbl 0768.76049
[16]	J. Happer, H. Brenner, Low Reynolds number hydroynamics, Martinus Nijhoff Publ. (1986) reprint
[17]	Höfler, K.; Schwarzer, S., Navier-Stokes simulation with constraint forces: finite-difference method for particle-laden flows and complex geometries, Phys. rev. E, 61, 7146-7160, (2000)
[18]	Kalthoff, W.; Schwarzer, S.; Herrmann, H.J., An algorithm for the simulation of particulate suspensions with inertia effects, Phys. rev. E, 56, 2234-2242, (1996)
[19]	Ladd, A.J.C., Numerical simulations of particulate suspensions via a discretized Boltzmann equation. part I. theoretical foundation, J. fluid mech., 271, 285-310, (1994) · Zbl 0815.76085
[20]	Ladd, A.J.C., Numerical simulations of particulate suspensions via a discretized Boltzmann equation. part II: numerical results, J. fluid mech., 271, 311-339, (1994) · Zbl 0815.76085
[21]	Ladd, A.J.C.; Verberg, R., Lattice-Boltzmann simulation of particle-fluid suspensions, J. statistical phys., 104, 1191-1251, (2001) · Zbl 1046.76037
[22]	Michaelides, E.E., Hydrodynamic force and heat/mass transfer from particles, bubbles and drops - the freeman scholar lecture, J. fluids eng., 125, 209-238, (2003)
[23]	Mohd-Yusof, J., Combined immersed boundaries/B-splines methods for simulations of flows in complex geometries, (1997), Annual Research Briefs, Center for Turbulence Research Stanford University
[24]	Peskin, C.S., Numerical analysis of blood flow in the heart, J. comput. phys., 25, 220-252, (1977) · Zbl 0403.76100
[25]	Peskin, C.S., The immersed boundary method, Acta numerica, 11, 479-517, (2002) · Zbl 1123.74309
[26]	ten Cate, A.; Nieuwstad, C.H; Derksen, J.J.; Van den Akker, H.E.A., Particle imaging velocimetry experiments and lattice-Boltzmann simulations on a single sphere settling under gravity, Phys. fluids, 14, 4012-4025, (2002) · Zbl 1185.76073
[27]	Tryggvason, T.; Bunner, B.; Esmaeeli, A.; Juric, D.; Al-Rawahi, N.; Tauber, W.; Han, J.; Nas, S.; Jan, Y.-J., A front-tracking method for the computations of multiphase flow, J. comput. phys., 169, 708-759, (2001) · Zbl 1047.76574
[28]	Zhang, Z.; Prosperetti, A., A method for particle simulation, J. appl. mech.-trans. ASME, 70, 64-74, (2003) · Zbl 1110.74798