The multiphase particle-in-cell method (MP-PIC) is a numerical method for modeling particle-fluid and particle-particle interactions in a computational fluid dynamics (CFD) calculation. The MP-PIC method achieves greater stability than its particle-in-cell predecessor by simultaneously treating the solid particles as computational particles and as a continuum. In the MP-PIC approach, the particle properties are mapped from the Lagrangian coordinates to an Eulerian grid through the use of interpolation functions. After evaluation of the continuum derivative terms, the particle properties are mapped back to the individual particles. This method has proven to be stable in dense particle flows, computationally efficient, and physically accurate. This has allowed the MP-PIC method to be used as particle-flow solver for the simulation of industrial-scale chemical processes involving particle-fluid flows.
History The multiphase particle-in-cell (MP-PIC) method was originally developed for a one-dimensional case in the mid-1990s by P.J. O'Rourke (Los Alamos National Laboratory), who also coined the term MP-PIC. Subsequent extension of the method to two-dimensions was performed by D.M. Snider and O'Rourke. By 2001, D.M. Snider had extended the MP-PIC method to full three-dimensions. Currently, the MP-PIC method is used in commercial software for the simulation of particle-fluid systems and also available in MFiX suite by NETL.
Method The MP-PIC method is described by the governing equations, interpolation operators, and the particle stress model.
Governing equations
Fluid phase The multiphase particle-in-cell method assumes an incompressible fluid phase with the corresponding continuity equation,
∂ θ f ∂ t + ∇ ⋅ ( θ f u f ) = 0 , {\displaystyle {\frac {\partial \theta _{f}}{\partial t}}+\nabla \cdot (\theta _{f}\mathbf {u} _{f})=0,}
where the θ f {\displaystyle \theta _{f}\;} is the fluid volume fraction and u f {\displaystyle \mathbf {u} _{f}} is the fluid velocity. Momentum transport is given by a variation of the Navier-Stokes equations where ρ f {\displaystyle \rho _{f}} is the fluid density, p {\displaystyle p} is the fluid pressure, and g {\displaystyle \mathbf {g} } is the body force vector (gravity).
∂ θ f u f ∂ t + ∇ ⋅ ( θ f u f u f ) = − ∇ p ρ f − F ρ f + θ f g {\displaystyle {\frac {\partial \theta _{f}\mathbf {u} _{f}}{\partial t}}+\nabla \cdot (\theta _{f}\mathbf {u} _{f}\mathbf {u} _{f})=-{\frac {\nabla p}{\rho _{f}}}-{\frac {\mathbf {F} }{\rho _{f}}}+\theta _{f}\mathbf {g} }
The laminar fluid viscosity terms, not included in the fluid momentum equation, can be included if necessary but will have a negligible effect on dense particle flow. In the MP-PIC method, the fluid motion is coupled with the particle motion through F {\displaystyle \mathbf {F} } , the rate of momentum exchange per volume between the fluid and particle phases. The fluid phase equations are solved using a finite volume approach.
… excerpt ends here. Continue reading the full article.
