Elsevier

Chemical Engineering Journal

Volume 308, 15 January 2017, Pages 544-556
Chemical Engineering Journal

Computational study of the bubbling-to-slugging transition in a laboratory-scale fluidized bed

https://doi.org/10.1016/j.cej.2016.08.113Get rights and content

Highlights

  • The bubbling to slugging transition occurs over a range of gas flow.

  • Unique spatiotemporal features emerge in the bubbling-to-slugging transition.

  • High-speed pressure measurements respond to unique spatiotemporal features.

  • Bubbling-to-slugging transition spatiotemporal features vary with axial location.

  • On-line bubbling-to-slugging transition diagnostic development with a model.

Abstract

We report results from a computational study of the transition from bubbling to slugging in a laboratory-scale fluidized-bed reactor with Geldart Group B glass particles. For simulating the three-dimensional fluidized-bed hydrodynamics, we employ MFiX, a widely studied multi-phase flow simulation tool, that uses a two-fluid Eulerian-Eulerian approximation of the particle and gas dynamics over a range of gas flows. We also utilize a previously published algorithm to generate bubble statistics that can be correlated with pressure fluctuations to reveal previously unreported details about the stages through which the hydrodynamics progress during the bubbling-to-slugging transition. We expect this new information will lead to improved approaches for on-line reactor diagnostics, as well as new approaches for validating the results of computational fluidized-bed simulations with experimental measurements.

Section snippets

Introduction and background

Gas-solid fluidized-bed reactors are widely used in the chemical industry, including biomass conversion [1], [2], [3], [4], petroleum refining [5], and pharmaceutical [6], [7] and commodity chemicals production [8]. For this reason, there is widespread interest in establishing a comprehensive understanding of the gas-solid hydrodynamics to optimize processes in which fluidized-bed reactors are key components. Three of the most important hydrodynamic states or flow regimes in fluidized-bed

Fluidized bed simulation conditions

To reflect a lab-scale reactor of current relevance, we assumed the geometry of an experimental laboratory reactor used for biomass processing research at the National Renewable Energy Laboratory (NREL). A schematic of the reactor is shown in Fig. 1. The inner diameter Dr and height of the reactor Hr are 0.0508 and 1.27 m, respectively. However, for simulation purposes, the computational domain was reduced to a height of 0.4 m. Operating conditions were chosen to match baseline experiments at

Results and discussion

The visual appearance of the bubbles generated by MFiX was observed graphically in terms of void fraction iso-surfaces using Paraview [86] as illustrated in the following section below. To systematically quantify the simulated bubble patterns, we evaluated bubble statistics derived from the MFiX void fraction output using MS3DATA. We then evaluated the pressure time series features at each level to determine how they relate to the observed bubble behavior.

Conclusions

Results from the three-dimensional computational simulations of a laboratory-scale fluidized bed indicate that the transition from bubbling-to-slugging is a complex process that occurs over a range of gas flows rather than abruptly at a single flow. The transition process appears to involve a cascade of bubble coalescences that produce size and speed changes which begin near the bed surface and then progress downward toward the distributor as gas flow increases. The state of maximum slugging

Acknowledgment

This work was supported by the Bioenergy Technology Office, US Department of Energy through the Computational Pyrolysis Consortium (CPC) project. The authors would like to thank program sponsors Jeremy Leong, Cynthia Tyler, and Kevin Craig for their support and guidance. More information about the CPC project can be found at http://cpcbiomass.org/. The authors acknowledge Akhilesh Bakshi for making MS3DATA available and See Hoon Lee for discussions through email about their work related to

References (98)

  • S. Nedeltchev

    New methods for flow regime identification in bubble columns and fluidized beds

    Chem. Eng. Sci.

    (2015)
  • J.P. Constantineau et al.

    Generalized bubbling–slugging fluidized bed reactor model

    Chem. Eng. Sci.

    (2007)
  • P. Lettieri et al.

    Predicting the transition from bubbling to slugging fluidization using computational fluid dynamics

    Chem. Eng. Res. Des.

    (2004)
  • J. Baeyens et al.

    An investigation into slugging fluidized beds

    Chem. Eng. Sci.

    (1974)
  • P.S.B. Stewart et al.

    Slug flow in fluidised beds

    Powder Technol.

    (1967)
  • T.C. Ho et al.

    The onset of slugging in gas solid fluidized-beds with large particles

    Powder Technol.

    (1983)
  • M. Rüdisüli et al.

    Bubble characterization in a fluidized bed by means of optical probes

    Int. J. Multiphase Flow

    (2012)
  • S.C. Saxena et al.

    Fluidization characteristics of gas-fluidized beds: air and glass beads system

    Energy

    (1989)
  • J. van der Schaaf et al.

    Non-intrusive determination of bubble and slug length scales in fluidized beds by decomposition of the power spectral density of pressure time series

    Int. J. Multiphase Flow

    (2002)
  • J.M. Weber et al.

    Bubbling fluidized bed characterization using electrical capacitance volume tomography (ECVT)

    Powder Technol.

    (2013)
  • S.C. Saxena et al.

    Pressure fluctuations in a gas fluidized bed and fluidization quality

    Energy

    (1990)
  • J.S. Halow et al.

    Observations of a fluidized bed using capacitance imaging

    Chem. Eng. Sci.

    (1993)
  • J.S. Halow et al.

    Observations of fluidized bed coalescence using capacitance imaging

    Powder Technol.

    (1992)
  • Y.T. Makkawi et al.

    Fluidization regimes in a conventional fluidized bed characterized by means of electrical capacitance tomography

    Chem. Eng. Sci.

    (2002)
  • C.C. Pain et al.

    A study of bubbling and slugging fluidised beds using the two-fluid granular temperature model

    Int. J. Multiphase Flow

    (2001)
  • S.J. Zhang et al.

    Computational investigation of slugging behaviour in gas-fluidised beds

    Powder Technol.

    (2002)
  • N. Reuge et al.

    Multifluid Eulerian modelling of a silicon fluidized bed chemical vapor deposition process: analysis of various kinetic models

    Chem. Eng. J. (Amsterdam, Neth.)

    (2009)
  • M.J.V. Goldschmidt et al.

    Hydrodynamic modelling of dense gas-fluidised beds using the kinetic theory of granular flow: effect of coefficient of restitution on bed dynamics

    Chem. Eng. Sci.

    (2001)
  • C. Loha et al.

    Euler-euler CFD modeling of fluidized bed: influence of specularity coefficient on hydrodynamic behavior

    Particuology

    (2013)
  • C. Loha et al.

    Effect of coefficient of restitution in Euler–Euler CFD simulation of fluidized-bed hydrodynamics

    Particuology

    (2014)
  • P. Fede et al.

    3D numerical simulation of a lab-scale pressurized dense fluidized bed focussing on the effect of the particle-particle restitution coefficient and particle–wall boundary conditions

    Chem. Eng. Sci.

    (2016)
  • N. Xie et al.

    Effects of using two- versus three-dimensional computational modeling of fluidized beds: Part ii, budget analysis

    Powder Technol.

    (2008)
  • A. Bakshi et al.

    Eulerian–Eulerian simulation of dense solid–gas cylindrical fluidized beds: impact of wall boundary condition and drag model on fluidization

    Powder Technol.

    (2015)
  • A. Bakshi et al.

    Study of the effect of reactor scale on fluidization hydrodynamics using fine-grid CFD simulations based on the two-fluid model

    Powder Technol.

    (2016)
  • A. Acosta-Iborra et al.

    Experimental and computational study on the bubble behavior in a 3-D fluidized bed

    Chem. Eng. Sci.

    (2011)
  • N. Xie et al.

    Effects of using two- versus three-dimensional computational modeling of fluidized beds – part i, hydrodynamics

    Powder Technol.

    (2008)
  • G. Wu et al.

    Lagrangian-Eulerian simulation of slugging fluidized bed

    Particuology

    (2012)
  • M. Ishii et al.

    Two-fluid model and hydrodynamic constitutive relations

    Nucl. Eng. Des.

    (1984)
  • N.G. Deen et al.

    Review of discrete particle modeling of fluidized beds

    Chem. Eng. Sci.

    (2007)
  • D.G. Schaeffer

    Instability in the evolution equations describing incompressible granular flow

    J. Differ. Equ.

    (1987)
  • C.M. Boyce et al.

    The origin of pressure oscillations in slugging fluidized beds: comparison of experimental results from magnetic resonance imaging with a discrete element model

    Chem. Eng. Sci.

    (2014)
  • A. Bakshi et al.

    Multiphase-flow statistics using 3D detection and tracking algorithm (MS3DATA): methodology and application to large-scale fluidized beds

    Chem. Eng. J. (Amsterdam, Neth.)

    (2016)
  • H.T. Bi

    A critical review of the complex pressure fluctuation phenomenon in gas–solids fluidized beds

    Chem. Eng. Sci.

    (2007)
  • F. Johnsson et al.

    Characterization of fluidization regimes by time-series analysis of pressure fluctuations

    Int. J. Multiphase Flow

    (2000)
  • K. Nguyen et al.

    Spatio-temporal dynamics in a train of rising bubbles

    Chem. Eng. J. Biochem. Eng. J.

    (1996)
  • A. Tufaile et al.

    Chaotic behavior in bubble formation dynamics

    Phys. A

    (2000)
  • J. van der Schaaf et al.

    Similarity between chaos analysis and frequency analysis of pressure fluctuations in fluidized beds

    Chem. Eng. Sci.

    (2004)
  • D. Kunii et al.

    Chapter 5 - Bubbles in Dense Beds, Fluidization Engineering

    (1991)
  • C. Rautenbach et al.

    A comparative study between electrical capacitance tomography and time-resolved X-ray tomography

    Flow Meas. Instrum.

    (2013)
  • Cited by (17)

    • Particle flow characteristics in a gas-solid separation fluidized bed based on machine learning

      2022, Fuel
      Citation Excerpt :

      The dissipation and interactions of the bubble phase and dense phase lead to a heterogeneous flow structure, and local bubble dynamics are characterized by chaos, which is typical of a complex system [11–13], particularly with uncertainties in the relationships among mesoscale bubbles, cluster structures, cavities and gas-solid two-phase flow characteristics and the dynamic transformation and nonuniform fluidization of different flow patterns during fluidization. Furthermore, particles in different regions have different flow forms [14]; even under the same operating conditions, the flow structure and state of the gas phase, solid phase, and gas-solid phase also change with spatial position in the bed, resulting in abnormal fluidization states such as particle back mixing [15], regional dead zone, cavitation, bridging, gully flow, slug and so on [16–19]. A study of the effect of particle flow characteristics on fluidization and separation would be beneficial to the post-industrial application of gas-solid separation fluidized bed technology.

    • Effect of bed size on the gas–solid flow characterized by pressure fluctuations in bubbling fluidized beds

      2019, Particuology
      Citation Excerpt :

      The standard deviation is usually used to represent pressure-fluctuation amplitudes and was used to determine the transition velocity from bubbling to turbulent fluidization (Bi, Ellis, Abba, & Grace, 2000; He et al., 2014) and the minimum fluidization velocity (Umf) (Felipe & Rocha, 2007), as well as to monitor the fluidization state (van Ommen, Korte, & Bleek, 2004). Skewness and kurtosis are also used to identify the transition from bubbling to slugging or turbulent fluidization (Lee & Kim, 1988; Lee, Sang, & Park, 2002; Ramirez et al., 2017). The frequency-domain analysis mainly concerns the physical interpretation of the main peak frequencies at lower frequencies (Kage, Agari, Ogura, & Matsuno, 2000) and the fall-off characteristics at higher frequencies (van der Schaaf, Johnsson, Schouten, & van den Bleek, 1999) in the power spectral density as well as the analysis of signal components in different frequency ranges corresponding to different scales by wavelet transform (Briens & Ellis, 2005; Zhao & Yang, 2003).

    View all citing articles on Scopus

    Notice: This manuscript has been authored by UT-Battelle, LLC, under Contract No. DE-AC0500OR22725 with the U.S. Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for the United States Government purposes. The Department of Energy will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doe-public-access-plan).

    View full text