Numerical simulation of gas–solid flow with two fluid model in a spouted-fluid bed

The flow characteristics in a spouted-fluid bed differ from those in spouted or fluidized beds because of the injection of the spouting gas and the introduction of a fluidizing gas. The flow behavior of gas–solid phases was predicted using the Eulerian–Eulerian two-fluid model (TFM) approach with kinetic theory for granular flow to obtain the flow patterns in spouted-fluid beds. The gas flux and gas incident angle have a significant influence on the porosity and particle concentration in gas–solid spouted-fluid beds. The fluidizing gas flux affects the flow behavior of particles in the fountain. In the spouted-fluid bed, the solids volume fraction is low in the spout and high in the annulus. However, the solids volume fraction is reduced near the wall.     

Simulations of vertical jet penetration using a filtered two-fluid model in a gas–solid fluidized bed

The influence of a vertical jet located at the distributor in a cylindrical fluidized bed on the flow behavior of gas and particles was predicted using a filtered two-fluid model proposed by Sundaresan and coworkers. The distributions of volume fraction and the velocity of particles along the lateral direction were investigated for different jet velocities by analyzing the simulated results. The vertical jet penetration lengths at the different gas jet velocities have been obtained and compared with predictions derived from empirical correlations; the predicted air jet penetration length is discussed. Agreement between the numerical simulations and experimental results has been achieved.     

Magnetic resonance studies of jets in a gas–solid fluidised bed

Magnetic resonance imaging (MRI) has been used to study the behaviour of jets at the distributor of a 50 mm diameter fluidised bed of 0.5 mm diameter poppy seeds. Two perforated-plate distributors were examined, containing either 10 or 14 holes, each 1 mm diameter. Ultra-fast MR imaging was able to show the transient nature of the upper parts of the jets, where discrete bubbles are formed. Imaging in 3D showed that the central jets were the longest for flow rates below minimum fluidisation. Above minimum fluidisation, the outer jets, nearest the wall of the fluidised bed, arched inward towards the central axis. In this latter case, interpretation of the time-averaged 3D image required the use of ultra-fast MR imaging to identify the approximate height above the distributor at which discrete bubbles were formed. The apparently continuous void extending along the central axis above this height in the time-averaged 3D image was thus identified, using ultra-fast MR imaging, as representing the averaged paths of released bubbles. Time-averaged MR velocity mapping was also used to identify dead zones of stationary particles resting on the distributor between the jets. The dead zones could be observed when the superficial velocity of the gas approached minimum fluidisation, but they were smaller than those observed at lower gas superficial velocity. Comparable images of a single jet through 1.2 mm diameter poppy seeds from MRI and electrical capacitance volume tomography (ECVT) are also demonstrated.     

Two and three dimensional modeling of fluidized bed with multiple jets in a DEM–CFD framework

Fluidized beds with multiple jets have widespread industrial applications. The objective of this paper is to investigate the jet interactions and hydrodynamics of a fluidized bed with multiple jets. Discrete element modeling coupled with in-house CFD code GenIDLEST has been used to simulate a bed with nine jets. The results are compared with published experiments. Mono dispersed particles of size 550 μm are used with 1.4 times the minimum fluidization velocity of the particles. Both two and three dimensional computations have been performed. To the best of our knowledge, the results presented in this paper are the first full 3D simulations of a fluidized bed performed with multiple jets. Discrepancies between the experiment and simulations are discussed in the context of the dimensionality of the simulations. The 2D solid fraction profile compares well with the experiment close to the distributor plate. At higher heights, the 2D simulation over-predicts the solid fraction profiles near the walls. The 3D simulation on the other hand is better able to capture the solid fraction profile higher up in the bed compared to that near the distributor plate. Similarly, the normalized particle velocities and the particle fluxes compare well with the experiment closer to the distributor plate for the 2D simulation and the freeboard for the 3D simulation, respectively. A lower expanded bed height is predicted in the 2D simulation compared to the 3D simulation and the experiment. The results obtained from DEM computations show that a 2D simulation can be used to capture essential jetting trends near the distributor plate regions, whereas a full scale 3D simulation is needed to capture the bubbles near the freeboard regions. These serve as validations for the experiment and help us understand the complex jet interaction and solid circulation patterns in a multiple jet fluidized bed system.     

Statistical diagnosis of a gas–solid fluidized bed using Electrical Capacitance Tomography

Fluidization experiments were performed using several particle size distributions of spherical glass particles, ranging from Geldart B to D. An Electrical Capacitance Tomography (ECT) tomograph was utilized in the present study and its usefulness as a diagnostic tool is illustrated. During the experiments a 10.4 cm diameter column was utilized and the column was operated at atmospheric pressure and room temperature (cold fluidized bed). Statistical analyses were performed on the average solid fraction data obtained using the ECT tomograph. Using the time domain skewness and kurtosis the time series could be characterised and the quality of fluidization is determined at different superficial gas velocities (Azizpour, H., Sotudeh-Gharebagh, R., Zarghami, R., Abbasi, M., Mostoufi, N., Mahjoob, M., 2011. Characterization of gas–solid fluidized bed hydrodynamics by vibration signal analysis. International Journal of Multiphase Flow, 37, 788–793). Statistical analysis is also used to characterise the influence of small particles on the bed hydrodynamics.     

Orientation of cylindrical particles in gas–solid circulating fluidized bed

The orientation of cylindrical particles in a gas–solid circulating fluidized bed was investigated by establishing a three-dimensional Euler–Lagrange model on the basis of rigid kinetics, impact kinetics and gas–solid two-phase flow theory. The resulting simulation indicated that the model could well illustrate the orientation of cylindrical particles in a riser during fluidization. The influences of bed structure and operation parameters on orientation of cylindrical particles were then studied and compared with related experimental results. The simulation results showed that the majority of cylindrical particles move with small nutation angles in the riser, the orientation of cylindrical particles is affected more obviously by their positions than by their slenderness and local gas velocities. The simulation results well agree with experiments, thus validating the proposed model and computation.     

Particle resolved direct numerical simulation of a liquid–solid fluidized bed: Comparison with experimental data

Particle-resolved direct numerical simulations of a 3-D liquid–solid fluidized bed experimentally investigated by Aguilar-Corona (2008) have been performed at different fluidization velocities (corresponding to a range of bed solid volume fraction between 0.1 and 0.4), using Implicit Tensorial Penalty Method. Particle Reynolds number and Stokes number are O(100) and O(10), respectively. In this paper, we compare the statistical quantities computed from numerical results with the experimental data obtained with 3-D trajectography and High Frequency PIV. Fluidization law predicted by the numerical simulations is in very good agreement with the experimental curve and the main features of trajectories and Lagrangian velocity signal of the particles are well reproduced by the simulations. The evolution of particle and flow velocity variances as a function of bed solid volume fraction is also well captured by the simulations. In particular, the numerical simulations predict the right level of anisotropy of the dispersed phase fluctuations and its independence of bed solid volume fraction. They also confirm the high value of the ratio between the fluid and the particle phase fluctuating kinetic energy. A quick analysis suggests that the fluid velocity fluctuations are mainly driven by fluid–particle wake interactions (pseudo-turbulence) whereas the particle velocity fluctuations derive essentially from the large scale flow motion (recirculation). Lagrangian autocorrelation function of particle fluctuating velocity exhibits large-scale oscillations, which are not observed in the corresponding experimental curves, a difference probably due to a statistical averaging effect. Evolution as a function of the bed solid volume fraction and the collision frequency based upon transverse component of particle kinetic energy correctly matches the experimental trend and is well fitted by a theoretical expression derived from Kinetic Theory of Granular Flows.     

Characterization of gas–solid fluidized bed hydrodynamics by vibration signature analysis

Vibration measurement, as a non-intrusive technique, was used to characterize the hydrodynamics of fluidized beds. A series of experiments were performed in a lab-scale fluidized bed using two accelerometers for measuring the vibration of the bed and a pressure probe for measuring pressure fluctuations. The output signals were analyzed by statistical methods. The results show that the vibration technique can predict transition velocities at high velocities and indicate that analyzing the vibration signals can be an effective non-intrusive technique to characterize the hydrodynamics of fluidized beds. It was shown that transition from bubbling to turbulent velocity can be determined from the variation of standard deviation and kurtosis of vibration signals against superficial gas velocity of the bed. However, this point could be determined only from standard deviation of pressure fluctuations, and not from skewness or kurtosis of pressure fluctuations.     

Anti-wear beam effects on gas–solid hydrodynamics in a circulating fluidized bed

Anti-wear beams installed on water walls of circulating fluidized bed (CFB) boilers are one of the most effective ways to protect against water-wall erosion. Beam effects from, for example, beam size and superficial gas velocity were investigated on gas–solid hydrodynamics in a CFB test rig using CFD simulations and experimental methods. The downward flow of the wall layer solids is observed to be disrupted by the beam but is then restored some distance further downstream. When falling solids from the wall layer hit the anti-wear beam, the velocity of the falling solids decreases rapidly. A fraction of the solids accumulates on the beam. Below the beams, the falling solids have reduced velocities but upward-moving solids were observed on the wall. The effect of the beam increases with width and superficial gas velocity. Wear occurs mainly above the beam and its variation with width is different above to below the beam. There is an optimum width that, when combined with beam height, results in less erosion.     

Lagrangian–Eulerian simulation of slugging fluidized bed

This work studies gas–solid slugging fluidized beds with Type-D particles, using two-dimensional simulations based on discrete element model (DEM). DEM performance is quantitatively validated by two commonly accepted correlations for determining slugging behavior. The voidage profiles simulated with bed height corresponding to Baeyens and Geldart (1974) correlation for onset of slugging demonstrate a transitional flow pattern from free bubbling to slugging. The present calculated values for the maximum slugging bed height are in good agreement with the correlation from Matsen et al. (1969). Simulations show that fluidized beds with Type-D particles can operate in the round-nosed slugging regime and also shows that wall slugs and square-nosed slugs tend to be formed with increase in superficial gas velocity and in bed height, respectively.     

Validation of CFD-DEM simulation of a liquid–solid fluidized bed by dynamic analysis of time series

Liquid–solid fluidized beds (LSFB) modeling validation is crucial for establishing design rules and monitoring tools. However, it generally relies on comparing global variables, which overlook dynamic features that influence reaction outputs. This work aims to implement time series analysis tools to compare Radioactive Particle Tracking data with a simulation consisting of Computational Fluid Dynamics coupled with Discrete-Element Method. Experiments have been performed in a pilot-scale LSFB of calcium alginate spheres fluidized with a calcium chloride solution. The Diks’ test indicates that the simulation can capture the LSFB behavior. It also allows diagnosing flow regime transitions from the simulation. Trends of solid dispersion coefficients and mixing times predicted by the simulation are in good agreement with the experiments.     

Improvement,validation and application of CFD/DEM model to dense gas–solid flow in a fluidized bed

Dense gas–solid flow with solid volume fraction greater than 10% and at moderate Reynolds number is important in many industrial facilities such as fluidized beds. In this work, the Euler–Lagrange approach in combination with a deterministic collision model is applied to a laboratory-scale fluidized bed. The fluid–particle interaction is studied using a new procedure called the offset method, which results in several numbers of spatial displacements of the fluid grid. The proposed method is highly precise in determining porosity and momentum transfer, thus improving simulation accuracy. A validation study was carried out to assess the results using this in-house CFD/DEM code against 5-s operation of a Plexiglas spouted-fluidized bed, showing good qualitative correlation of solid distribution in the bed and acceptable quantitative agreement of pressure drops at different positions in the bed. In view of high computing cost, special emphasis is placed on effective program design, such as application of advanced detection algorithm for particle–particle/wall collisions, the multi-grid method and parallel calculation. In this context, the influence of increasing the processor number, up to 36, on calculation efficiency was investigated.     

Numerical simulation of counter-current flow field in the downcomer of a liquid–solid circulating fluidized bed

A comprehensive study on the hydrodynamics in the downcomer of a liquid–solid circulating fluidized bed (LSCFB) is crucial in the control and optimization of the extraction process using an ion exchange LSCFB. A computational fluid dynamics model is proposed in this study to simulate the counter-current two-phase flow in the downcomer of the LSCFB. The model is based on the Eulerian–Eulerian approach incorporating the kinetic theory of granular flow. The predicted results agree well with our earlier experimental data. Furthermore, it is shown that the bed expansion of the particles in the downcomer is directly affected by the superficial liquid velocity in downcomer and solids circulation rate. The model also predicts the residence time of solid particles in the downcomer using a pulse technique. It is demonstrated that the increase in the superficial liquid velocity decreases the solids dispersion in the downcomer of the LSCFB.     

The characteristics of gas–solid flow and wall heat transfer in a fluidized bed reactor

Numerical study using computational fluid dynamics has been carried out to investigate the heat transfer characteristics of a laboratory fluidized bed reactor. The fluidized bed reactor of vTI (Johann Heinrich von Thünen-Institute)-Institute of Wood Technology and Wood Biology is modeled. For the simulation of multiphase flow and thermal fields, an Eulerian–Eulerian approach is applied. The flow and thermal characteristics of the reactor are fully investigated for the wide range of superficial gas velocities and two different particle diameters. In particular, the contributions of the gas bubble and emulsion phase flows on the wall heat transfer are scrutinized. From the predicted results, it is fully elucidated that particular near-wall bubble motions mainly govern the wall heat transfer.     

Apparent viscosity of mixed particle system in pulsed gas–solid separation fluidized bed

The apparent viscosity reflects the resistance of the fluidized medium in the bed to the beneficiation particles, which directly affects the separation time and mismatch content. So, the falling-ball method was used to measure the apparent viscosity of a binary medium in a pulsed fluidized bed by varying the gas velocity, pulsation frequency, and fine particle content. The results show that with increasing gas velocity and fine particle content, the apparent viscosity of the bed gradually decreased, whereas it first decreased and then increased with pulsation frequency increasing and achieved a minimum value in the range of 4–6 Hz. Within limits, the adjustment of gas velocity and fine content can effectively reduce the apparent viscosity and improve the separation process. A model for predicting the apparent viscosity in a pulsation separation fluidized bed was established with good accuracy.     

A new drag model for TFM simulation of gas–solid bubbling fluidized beds with Geldart-B particles

In this work, a new drag model for TFM simulation in gas–solid bubbling fluidized beds was proposed, and a set of equations was derived to determine the meso-scale structural parameters to calculate the drag characteristics of Geldart-B particles under low gas velocities. In the new model, the meso-scale structure was characterized while accounting for the bubble and meso-scale structure effects on the drag coefficient. The Fluent software, incorporating the new drag model, was used to simulate the fluidization behavior. Experiments were performed in a Plexiglas cylindrical fluidized bed consisting of quartz sand as the solid phase and ambient air as the gas phase. Comparisons based on the solids hold-up inside the fluidized bed at different superficial gas velocities, were made between the 2D Cartesian simulations, and the experimental data, showing that the results of the new drag model reached much better agreement with experimental data than those of the Gidaspow drag model did.     

Numerical study on hydrodynamics of gas–solids circulating fluidized bed with L-valve

L-valve is often used as a non-mechanical valve for the circulation of solids in gas–solids fluidized bed (GSFB) due to its advantages in simple construction and easy control. The information on solids circulation rate as well as the hydrodynamics performance of the CFB with L-valve is of great importance for its better control and design. This paper proposes a Eulerian-Eulerian approach based numerical model integrating the computational fluid dynamics (CFD) with turbulent model, the kinetic theory of granular flow (KTGF) and the drag model, thus the solids circulation rate and the local phase velocity as well as solids volume fraction can be predicted simultaneously. With this model, the hydrodynamics performance of the full loop GSCFB with a L-valve is analyzed in detail. It is found that the drag model affects the simulation significantly and the (energy minimization multiscale) EMMS method shows good performance in the full-loop simulation of GSCFB.     

Influence of vertical heat exchanger tubes,their arrangement and the column diameter on the hydrodynamics in a gas–solid bubbling fluidized bed

The hydrodynamic behavior of a cold-flow gas–solid fluidized bed with an inner diameter of 22 cm is investigated by means of an ultra-fast X-ray tomographic setup. In the case of an exothermal reaction, heat exchanger tubes are required to remove the reaction heat out of the bubbling fluidized bed reactor. For the examined cold-flow model, the heat exchanger tubes are replaced by vertical internals that serve as placeholder. The influence of vertical internals on the bubble properties for different spatial configurations (square and circular arrangements) is investigated in addition to measurements without internals. Furthermore, the hydrodynamic results of the Ø 22 cm column are compared with an available data set which is based on measurements that were conducted in a column with an inner diameter of 14 cm. The objective of this paper is to provide measurement data for the scale-up process as well as for various computer models simulating a bubbling fluidized bed with vertical internals. It was found that the scale-up process from pilot plants to an industrial scale may be simplified if vertical internals are present, independently of the geometric arrangement.