Lagrangian numerical simulation of particulate flows

The Lagrangian numerical simulation (LNS) scheme presented in this paper is motivated by the multiphase particle-in-cell (MP-PIC). In this numerical s cheme we solve the fluid phase continuity and momentum equations on an Eule rian grid. The particle motion is governed by Newton's law thus following t he Lagrangian approach. Momentum exchange from the particle-to-fluid is mod eled in the fluid phase momentum equation. Forces acting on the particle in clude drag from the fluid, body force and force due to interparticle stress . There is a freedom to use different models for these forces and to introd uce other forces. The effect of viscous stresses are included in the fluid phase equations. The volume fraction of the particles appear in the fluid p hase continuity and momentum equations. A finite volume method is used to s olve for the fluid phase equations on an Eulerian grid. Particle positions are updated using the Runge-Kutta scheme. This numerical scheme can handle a range of particle loadings and particle types. The LNS scheme is implemen ted using an efficient three-dimensional time-dependent finite volume algor ithm. We use a Chorin-type pressure-correction based fractional-step scheme on a non-staggered cartesian grid. In this paper, we consider only incompr essible Newtonian suspending fluid. However, the average velocity field of the fluid phase is not divergence-free because its effective density is not constant. Our pressure-correction based fractional-step scheme accounts fo r varying properties in the fluid phase equations. This method can also acc ount for suspending fluids with non-constant properties. The numerical sche me is verified by comparing results with test cases and experiments. (C) 20 01 Elsevier Science Ltd. All rights reserved.