 Multiphase particleincell method

The multiphase particleincell method (MPPIC) is a numerical method for modeling particlefluid and particleparticle interactions in a computational fluid dynamics (CFD) calculation. The MPPIC method achieves greater stability than its particleincell predecessor by simultaneously treating the solid particles as computational particles and as a continuum. In the MPPIC 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.^{[1]} This method has proven to be stable in dense particle flows (>5% by volume), computationally efficient^{[2]}, and physically accurate^{[3]}. This has allowed the MPPIC method to be used as particleflow solver for the simulation of industrialscale chemical processes involving particlefluid flows.
Contents
History
The multiphase particleincell (MPPIC) method was originally developed for a onedimensional case in the mid1990s by P.J. O'Rourke (Los Alamos National Laboratory)^{[1]}, who also coined the term MPPIC. Subsequent extension of the method to twodimensions was performed by D.M. Snider and O'Rourke^{[4]}. By 2001, D.M. Snider had extended the MPPIC method to full threedimensions^{[2]}. Currently, the MPPIC method is used in commercial software for the simulation of particlefluid systems.
Method
The MPPIC method is described by the governing equations, interpolation operators, and the particle stress model.
Governing equations
Fluid phase
The multiphase particleincell method assumes an incompressible fluid phase with the corresponding continuity equation,
where the is the fluid volume fraction and is the fluid velocity. Momentum transport is given by a variation of the NavierStokes equations where is the fluid density, is the fluid pressure, and is the body force vector (gravity).
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 MPPIC method, the fluid motion is coupled with the particle motion through , the rate of momentum exchange per volume between the fluid and particle phases. The fluid phase equations are solved using a finite volume approach.
Particle phase
The particle phase is described by a probability distribution function (PDF), which indicates the likelihood of finding a particle with a velocity , particle density , particle volume at location and time . The particle PDF changes in time as described by
where is the particle acceleration.
A numerical solution of the particle phase is obtained by dividing the distribution into a finite number of "computational particles" that each represent a number of real particles with identical mass density, volume, velocity and location. At each time step, the velocity and location of each computational particle are updated using a discretized form of the above equations. The use of computational particles allows for a significant reduction in computational requirements with a negligible impact on accuracy under many conditions. The use of the computational particle in the Multiphase ParticleinCell method allows a full particle size distribution (PSD) to be modeled within the system as well as the modeling of polydisperse solids.^{[5]}
Identities of the particle probability distribution function
The following local particle properties are determined from integrating the particle probability distribution function:
 Particle volume fraction:
 Average particle density:
 Mean particle velocity:
Interphase coupling
The particle phase is coupled to the fluid phase through the particle acceleration term, , defined as
In the acceleration term, is determined from the particle drag model and is determined from the interparticle stress model.
The momentum of the fluid phase is coupled to the particle phase through the rate of momentum exchange, . This is defined from the particle population distribution as
Interpolation operators
The transfer of particle properties between the Lagrangian particle space and the Eulerian grid is performed using linear interpolation functions. Assuming a rectilinear grid consisting of rectangular cuboid cells, the scalar particle properties are interpolated to the cell centers while the vector properties are interpolated to cell faces. In three dimensions, trilinear interpolation functions and definitions for the products and gradients of interpolated properties are provided by Snider for three dimensional models^{[2]}.
Particle stress model
The effects of particle packing are modeled in the MPPIC method with the use of a function of particle stress. Snider (2001) has suggested calculating the particle stress , as
where is the closepack volume fraction and , , and are constants.
Limitations of the multiphase particleincell method
 Particle Shape  In the MPPIC method, all particles are assumed to be spherical. Corrections for nonspherical particles can be included in particle drag model but for highly nonspherical particles, the true interactions may not be well represented.
 Particle size with respect to grid size  The size of particles must be small compared to the Eularian grid in the MPPIC approach for accurate interpolation.
Extensions
 Chemical reactions – Coupling the local Eulerian values for fluid velocity in the MPPIC method with equations for diffusional mass transfer allows the transport of a chemical species within the fluidparticle system to be modeled. Reaction kinetics dependent on particle density, surface area, or volume can be included as well for applications in catalysis^{[6]}, gasification^{[7]}, or solid deposition.
 Liquid Injection  MPPIC method was extended by Zhao, O'Rourke, and Snider to model the coating of particle with a liquid^{[8]}.
 Thermal Modeling  Conductive and convective heat transfer can be included by coupling MPPIC variables with equations for heat transfer. Commercial implementations of MPPIC method include radiative heat transfer as well^{[9]}.
Applications
 Biomass gasifiers ^{[10]}
 Chemical looping combustion (CLC)^{[11]}^{[12]}^{[13]}^{[14]}^{[15]}
 Circulating Fluidized Bed Combustion ^{[16]}
 Coal gasifiers^{[7]}^{[17]}
 Cyclones^{[18]}
 Fluid Catalytic Cracking Reactors and Regenerators
 Fluidized bed dryers^{[19]}
 Fluidized Bed Reactors^{[20]}
 Liquidsolid settlers^{[5]}
 Metal casting ^{[13]}^{[21]}^{[22]}
 Particle jets^{[23]}
 Polysilicon deposition ^{[24]}
 Spray Coating^{[8]}
Software
 Barracuda by CPFD Software
References
 ^ ^{a} ^{b} Andrews, M.J. and O'Rourke, P.J. (1996). The Multiphase ParticleinCell (MPPIC) Method for Dense Particle Flows. International Journal of Multiphase Flow, 22(2):379–402.
 ^ ^{a} ^{b} ^{c} Snider, D.M. (2001). An Incompressible ThreeDimensional Multiphase ParticleinCell Model for Dense Particle Flows. Journal of Computational Physics, 170:523–549.
 ^ Snider, D. (2007). Three fundamental granular flow experiments and CPFD predictions. Powder Technology 176: 3646.
 ^ Snider, D.M., O'Rourke, P.J., and Andrews, M.J. (1997). An Incompressible TwoDimensional Multiphase ParticleInCell Model for Dense Particle Flows, NM, LA17280MS (Los Alamos National Laboratories, Los Alamos, NM)
 ^ ^{a} ^{b} Sundaresan, S. (2010). Challenges in the Analysis of HighVelocity GasParticle Flows in Large Devices, University of Houston Neal Amundson Memorial Lecture Series, 2010.
 ^ Snider, D. and Banerjee, S. (2010). Heterogeneous gas chemistry in the CPFD Eulerian–Lagrangian numerical scheme (ozone decomposition). Powder Technology 199(1):100–106
 ^ ^{a} ^{b} Snider, D.M., Clark, S.M., O'Rourke, P.J. (2011). Eulerian–Lagrangian method for threedimensional thermal reacting flow with application to coal gasifiers. Chemical Engineering Science 66:1285–1295.
 ^ ^{a} ^{b} Zhao, P., O'Rourke, P.J., Snider, D. Threedimensional simulation of liquid injection, film formation and transport, in fluidized beds. Particuology 7:337346
 ^ CPFD Software, LLC. Barracuda 14.4 Released. http://www.cpfdsoftware.com/news/barracuda_14.4_released Retrieved Feb 8, 2011
 ^ Blaser, P. and Chandran, R. (2009). Computational Simulation of Fluidization Dynamics Inside a Commercial Biomass Gasifier. AIChE 2009 Annual Meeting.
 ^ Williams, K., Snider, D., Guenther, C. (2010) CFD Simulations of the NETL Chemical Looping Experiment, AIChE 2010 National Meeting, http://www.aicheproceedings.org/2010/Fall/data/papers/Paper202402.html Retrieved Feb 8, 2011
 ^ Snider, D., Guenther, C., Dalton J., Williams, K. (2010) CPFD EulerianLagrangian Numerical Scheme Applied to the NETL Benchtop Chemical Looping Experiment. Proceedings of the 1st International Conference on Chemical Looping
 ^ ^{a} ^{b} Yeomans, N., and Blaser, P. (2006). Predicting the Process, Foundry Management & Technology, January 2006, pp 48–49.
 ^ Blaser, P., and Yeomans, N. (2006). Sand Core Engineering & Process Modeling, Japan Foundry Society, Vol. 2, No. 2, February 2006, pp. 420–427.
 ^ Schleg, P. (2003). Technology of Metalcasting, American Foundry Society, Des Plaines, IL, pp. 1 and 39.
 ^ Weng, M., Nies, M., and Plackmeyer, J. (2010). Comparison between Measurements and Numerical Simulation of Particle Flow and Combustion at the CFBC Plant Duisburg. 5. Internationaler VGBWorkshop "Betriebserfahrungen mit Wirbelschichtfeuerungen 2010"
 ^ Snider, D., Clark, S.(2009). CPFD EulerianLagrangian Method for Three Dimensional Thermal Reacting Flow. 2009 AIChE National Meeting, http://www.aicheproceedings.org/2009/Fall/data/papers/Paper149130.html Retrieved Feb 19, 2011
 ^ Williams, K., Snider, D., Badalassi, V., Reddy Karri, S.B., Knowlton, T.M., and Cocco, R.A. (2006). Computational Particle Fluid Dynamics Simulations and Validation for Cyclones: High and Low Loadings. AIChE 2006 National Meeting http://aiche.confex.com/aiche/2006/preliminaryprogram/abstract_76001.htm Retrieved Feb. 19, 2011
 ^ Cocco, R. and Williams, K. (2004). Optimization of Particle Residence Time Inside Commercial Dryers with Arenaflow. AIChE 2004 National Meeting
 ^ Karimipour, S. and Pugsley, T. (2009). Application of the ParticleinCell Approach for the Simulation of Bubbling Fluidized Beds of Geldhart A Particles, Seventh International Conference on CFD in the Minerals and Process Industries.
 ^ Lefebvre, D., Mackenbrock, A., Vidal, V., and Haigh, P. (2005). Development and use of simulation in the design of blown cores and moulds, Foundry Trade Journal, February 2005.
 ^ Winartomo, B., Vroomen, U., and BuhrigPolaczek, A., Pelzer, M. (2005). Multiphase modeling of core shooting processes, International Journal of Cast Metals Research, Vol. 18, No. 1.
 ^ O'Rourke, P.J., Snider, D.M. (2010). An improved collision damping time for MPPIC calculations of dense particle flows with applications to polydisperse sedimenting beds and colliding particle jets. Chemical Engineering Science, 65:6014–6028.
 ^ Parker, J. (2011). Validation of CFD Model for Polysilicon Deposition and Production of Silicon Fines in a Silane Deposition FBR, International Journal of Chemical Reactor Engineering, Vol. 9, A40
Categories: Computational physics
 Numerical differential equations
 Computational fluid dynamics
Wikimedia Foundation. 2010.