Simultaneous solution algorithms for Eulerian-Eulerian gas-solid flow models: Stability analysis and convergence behaviour of a point and a plane solver

  • Authors:

  • Juray De Wilde;Jan Vierendeels;Geraldine J. Heynderickx;Guy B. Marin

  • Affiliations:

  • Laboratorium voor Petrochemische Techniek, Ghent University, Krijgslaan 281, Blok S5, B-9000 Ghent, Belgium and Fluid Mechanics Laboratory, Department of Fluid, Heat and Combustion Mechanics, Ghen ...;Fluid Mechanics Laboratory, Department of Fluid, Heat and Combustion Mechanics, Ghent University, St.-Pietersnieuwstraat 41, B-9000 Ghent, Belgium;Laboratorium voor Petrochemische Techniek, Ghent University, Krijgslaan 281, Blok S5, B-9000 Ghent, Belgium;Laboratorium voor Petrochemische Techniek, Ghent University, Krijgslaan 281, Blok S5, B-9000 Ghent, Belgium

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2005

Quantified Score

Hi-index 31.45

Visualization

Abstract

Simultaneous solution algorithms for Eulerian-Eulerian gas-solid flow models are presented and their stability analyzed. The integration algorithms are based on dual-time stepping with fourth-order Runge-Kutta in pseudo-time. The domain is solved point or plane wise. The discretization of the inviscid terms is based on a low-Mach limit of the multi-phase preconditioned advection upstream splitting method (MP-AUSMP). The numerical stability of the simultaneous solution algorithms is analyzed in 2D with the Fourier method. Stability results are compared with the convergence behaviour of 3D riser simulations. The impact of the grid aspect ratio, preconditioning, artificial dissipation, and the treatment of the source terms is investigated. A particular advantage of the simultaneous solution algorithms is that they allow a fully implicit treatment of the source terms which are of crucial importance for the Eulerian-Eulerian gas-solid flow models and their solution. The numerical stability of the optimal simultaneous solution algorithm is analyzed for different solids volume fractions and gas-solid slip velocities. Furthermore, the effect of the grid resolution on the convergence behaviour and the simulation results is investigated. Finally, simulations of the bottom zone of a pilot-scale riser with a side solids inlet are experimentally validated.