The effect of numerical diffusion and the influence of computational grid over gas-solid two-phase flow in a bubbling fluidized bed

The numerical diffusion effects appear due to the discretization process of the convective terms of the transport equations. This phenomenon takes place also in the numerical simulation of gas-solid two-phase flows in bubbling fluidized beds (BFB). In the present work a comparative analysis of the numerical results obtained using two interpolation schemes for convective terms, namely FOUP (First Order UPwind) and a high order scheme (Superbee) is presented. The equations are derived by considering the Eulerian-Eulerian gas-solid two-fluid model and the kinetic theory of granular flows (KTGF) for modeling solid phase constitutive equations. For that purpose the MFIX (Multiphase Flow with Interphase eXchanges) code developed at NETL (National Energy Technology Laboratory, US Department of Energy) is used. The numerical diffusion is analyzed by considering a single bubbling detachment and its hydrodynamic process in a two-dimensional BFB. The bubble shape is used as a metric for the description of the results. The influence of the computational grid is also analyzed. It is concluded that the Superbee scheme produces better results and this scheme is recommended for discretizations of the convective terms in coarse grids. The FOUP scheme can be used only in fine grids but it requires a high computational effort. In this study it is also verified that the analysis about estimating uncertainty in grid refinement can be applied in specific points of the grid when a monotonic convergence in time and space occurs.

[1]  Richard S. Varga,et al.  Application of Oscillation Matrices to Diffusion-Convection Equations , 1966 .

[2]  R. Jackson,et al.  Frictional–collisional constitutive relations for granular materials, with application to plane shearing , 1987, Journal of Fluid Mechanics.

[3]  R. Jackson,et al.  The mechanics of fluidized beds. I. The stability of the state of uniform fluidization , 1963 .

[4]  B. V. Leer,et al.  Towards the Ultimate Conservative Difference Scheme , 1997 .

[5]  R. Bagnold Experiments on a gravity-free dispersion of large solid spheres in a Newtonian fluid under shear , 1954, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[6]  S. L. Soo,et al.  Fluid dynamics of multiphase systems , 1967 .

[7]  J. Jenkins,et al.  A theory for the rapid flow of identical, smooth, nearly elastic, spherical particles , 1983, Journal of Fluid Mechanics.

[8]  B. P. Leonard,et al.  Beyond first‐order upwinding: The ultra‐sharp alternative for non‐oscillatory steady‐state simulation of convection , 1990 .

[9]  Goodarz Ahmadi,et al.  An equation of state for dense rigid sphere gases , 1986 .

[10]  Jianmin Ding,et al.  Prediction of Radial Distribution Function of Particles in a Gas−Solid Fluidized Bed Using Discrete Hard-Sphere Model , 2009 .

[11]  He Yurong,et al.  Numerical simulations of flow behavior of gas and particles in spouted beds using frictional-kinetic stresses model , 2009 .

[12]  J. Zhu,et al.  A local oscillation-damping algorithm for higher-order convection schemes , 1988 .

[13]  George Bergeles,et al.  DEVELOPMENT AND ASSESSMENT OF A VARIABLE-ORDER NON-OSCILLATORY SCHEME FOR CONVECTION TERM DISCRETIZATION , 1998 .

[14]  Xuanhui Qu,et al.  Effect of multiple impacts on high velocity pressed iron powder , 2009 .

[15]  Guy T. Houlsby,et al.  The flow of granular materials—II Velocity distributions in slow flow , 1982 .

[16]  John Garside,et al.  Velocity-Voidage Relationships for Fluidization and Sedimentation in Solid-Liquid Systems , 1977 .

[17]  P. Sweby High Resolution Schemes Using Flux Limiters for Hyperbolic Conservation Laws , 1984 .

[18]  Yassir Makkawi,et al.  A model for gas–solid flow in a horizontal duct with a smooth merge of rapid–intermediate–dense flows , 2006 .

[19]  M. Syamlal,et al.  The effect of numerical diffusion on simulation of isolated bubbles in a gas-solid fluidized bed , 2001 .

[20]  T. B. Anderson,et al.  Fluid Mechanical Description of Fluidized Beds. Equations of Motion , 1967 .

[21]  B. P. Leonard,et al.  A stable and accurate convective modelling procedure based on quadratic upstream interpolation , 1990 .

[22]  A. W. Jenike,et al.  A theory of flow of particulate solids in converging and diverging channels based on a conical yield function , 1987 .

[23]  D. Drew Mathematical Modeling of Two-Phase Flow , 1983 .

[24]  R. Courant,et al.  On the solution of nonlinear hyperbolic differential equations by finite differences , 1952 .

[25]  P. Gaskell,et al.  Curvature‐compensated convective transport: SMART, A new boundedness‐ preserving transport algorithm , 1988 .

[26]  B. P. Leonard,et al.  Simple high-accuracy resolution program for convective modelling of discontinuities , 1988 .

[27]  Ulrich Renz,et al.  Eulerian simulation of bubble formation at a jet in a two-dimensional fluidized bed , 1997 .

[28]  James D. Murray,et al.  On the mathematics of fluidization Part 1. Fundamental equations and wave propagation , 1965, Journal of Fluid Mechanics.

[29]  Jam Hans Kuipers,et al.  Critical comparison of hydrodynamic models for gas-solid fluidized beds - Part II: freely bubbling gas-solid fluidized beds , 2005 .

[30]  van den Cm Bleek,et al.  Eulerian simulations of bubbling behaviour in gas-solid fluidised beds , 1998 .

[31]  R. Collins,et al.  An extension of Davidson's theory of bubbles in fluidized beds , 1965 .

[32]  Jam Hans Kuipers,et al.  Critical comparison of hydrodynamic models for gas-solid fluidized beds - Part I: bubbling gas-solid fludized beds operated with a jet , 2005 .

[33]  Prabhu R. Nott,et al.  Frictional–collisional equations of motion for participate flows and their application to chutes , 1990, Journal of Fluid Mechanics.

[34]  P. Roache Perspective: A Method for Uniform Reporting of Grid Refinement Studies , 1994 .

[35]  Dimitri Gidaspow,et al.  Fluidization in Two-Dimensional Beds with a Jet. 2. Hydrodynamic Modeling , 1983 .

[36]  Jam Hans Kuipers,et al.  Computer simulation of the hydrodynamics of a two-dimensional gas-fluidized bed , 1993 .

[37]  K. E. Starling,et al.  Equation of State for Nonattracting Rigid Spheres , 1969 .

[38]  D. Jeffrey,et al.  Kinetic theories for granular flow: inelastic particles in Couette flow and slightly inelastic particles in a general flowfield , 1984, Journal of Fluid Mechanics.

[39]  David G. Schaeffer,et al.  Instability in the evolution equations describing incompressible granular flow , 1987 .

[40]  J. X. Bouillard,et al.  Hydrodynamics of fluidization: Fast-bubble simulation in a two-dimensional fluidized bed , 1991 .

[41]  Rajamani Krishna,et al.  Validation of the Eulerian simulated dynamic behaviour of gas-solid fluidised beds , 1999 .

[42]  He Yurong,et al.  Computer simulations of gas–solid flow in spouted beds using kinetic–frictional stress model of granular flow , 2004 .

[43]  H. Enwald,et al.  Eulerian two-phase flow theory applied to fluidization , 1996 .

[44]  Chia-Jung Hsu Numerical Heat Transfer and Fluid Flow , 1981 .

[45]  M. Syamlal,et al.  MFIX documentation theory guide , 1993 .

[46]  D. Gidaspow Multiphase Flow and Fluidization , 1994 .

[47]  D. Spalding A novel finite difference formulation for differential expressions involving both first and second derivatives , 1972 .