Experimental study and modeling of particle drying in a continuously-operated horizontal fluidized bed

Multistage fluidized beds are frequently used for product drying in industry. One advantage of these fluidized beds is that they can achieve a high throughput, when operated continuously. In this study, γ-Al₂O₃ particles were dried in a pilot-scale horizontal fluidized bed, without considering any comminution effects. For each experiment, the particle moisture content distribution and residence time distribution were determined. To take into account particle back mixing in our experiments, a one-dimensional population balance model that considers particle residence time was introduced into a fluidized bed-drying model. Experimental particle residence time distributions were reproduced using a tank-in-series model. Subsequently, the moisture content distribution was implemented, as a second dimension to the population balance in this model. These two-dimensional simulations were able to describe the experimental data, especially the spread in the residual particle moisture distribution, much more accurately than one-dimensional simulations. Using this novel two-dimensional model, the effects of different operating parameters (process gas temperature, solid feed rate, superficial air velocity) on the particle moisture content distribution were systematically studied.     

Optimization of Fine Solid Drying in Bubble Fluidized Bed

A commonly used method to dry fine solid particles is drying in a fluidized bed. This paper presents the optimization problem of fluidized drying of fine solids. A drying process proceeding in a three-stage cascade of fluidized cross-current dryers was considered. Solid flows from stage to stage, and fresh gas is introduced to each stage of the cascade. The hydrodynamics of bubble fluidized bed and kinetics of heat and mass transfer are taken into account. The bed hydrodynamics is described by a two-phase model. The drying process considered proceeds in the second period of drying. To optimize this problem a generalized version of a discrete algorithm with constant Hamiltonian was used. The optimization procedure is presented in the paper. In optimization calculations, gas parameters (temperature, humidity and flow rate) minimizing total process cost are sought. The results of calculation are presented as graphs. The results obtained and the conclusions drawn are discussed.     

Simulation of mass balance behavior in a large-scale circulating fluidized bed reactor

We determine using a compound model the influence of the mass of granular matter on the behavior of a supercritical circulating fluidized bed (CFB) reactor. Population balance enables a stationary-regime modeling of the mass flow of granular matter inside a CFB unit in a large-scale. The simulation includes some important dynamic processes of gas-particle flows in fluidized bed such as attrition, fragmentation, elutriation, and fuel combustion. Numerical calculations with full boiler loading were performed of operational parameters such as furnace temperature, furnace pressure, feeding materials mass flows, and excess air ratio. Furthermore, three bed inventory masses were adopted as experimental variables in the simulation model of mass balance. This approach enables a sensitivity study of mass flows of granular matter inside a CFB facility. Some computational results from this population balance model obtained for a supercritical CFB reactor are presented that show consistency with the operational data for large-scale CFB units.     

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 GenlDLEST 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.     

A model for expansion ratio in liquid-solid fluidized beds

Liquid-solid fluidized beds are used in mineral processing industries to separate particles based on parti- cle size, density, and shape. Understanding the expanded fluidized bed is vital for accurately assessing its performance. Expansion characteristics of the fluidized bed were studied by performing several experi- ments with iron ore, chromite, quartz, and coal samples. Using water as liquid medium, experiments were conducted to study the effects of particle size, particle density, and superficial velocity on fluidized bed expansion. The experimental data were utilized to develop an empirical mathematical model based on dimensional analysis to estimate the expansion ratio of the fluidized bed in terms of particle character- istics, operating and design parameters. The predicted expansion ratio obtained from the mathematical model is in good agreement with the experimental data.     

Insights in hydrodynamics of bubbling fluidized beds at elevated pressure by DEM–CFD approach

A numerical simulation was conducted to study the effect of pressure on bubble dynamics in a gas–solid fluidized bed. The gas flow was modeled using the continuum theory and the solid phase, by the discrete element method (DEM). To validate the simulation results, calculated local pressure fluctuations were compared with corresponding experimental data of 1-mm polyethylene particles. It was shown that the model successfully predicts the hydrodynamic features of the fluidized bed as observed in the experiments. Influence of pressure on bubble rise characteristics such as bubble rise path, bubble stability, average bubbles diameter and bubble velocity through the bed was investigated. The simulation results are in conformity with current hydrodynamic theories and concepts for fluidized beds at high pressures. The results show further that elevated pressure reduces bubble growth, velocity and stability and enhances bubble gyration through the bed, leading to change in bed flow structure.     

Influence of rolling friction on single spout fluidized bed simulation

In this paper we study the effect of rolling friction on the dynamics in a single spout fluidized bed using Discrete Element Method (DEM) coupled to Computational Fluid Dynamics (CFD). In a first step we neglect rolling friction and show that the results delivered by the open source CFD–DEM framework applied in this study agree with previous simulations documented in literature. In a second step we include a rolling friction sub-model in order to investigate the effect of particle non-sphericity. The influence of particle–particle as well as particle–wall rolling friction on the flow in single spout fluidized bed is studied separately. Adequate rolling friction model parameters are obtained using first principle DEM simulations and data from literature. Finally, we demonstrate the importance of correct modelling of rolling friction for coupled CFD–DEM simulations of spout fluidized beds. We show that simulation results can be improved significantly when applying a rolling friction model, and that experimental data from literature obtained with Positron Emission Particle Tracking (PEPT) technique can be satisfactorily reproduced.     

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.     

Experimental characterization and CFD simulation of gas maldistribution in turbulent fluidization of Group A particles

The study presented hereby investigates experimentally and with CFD simulations the gas distribution effect on the hydrodynamic of a Geldart Group A turbulent fluidized bed. Experiments were carried out on a cold flow fluidized bed column with an even and uneven gas distribution. Local solid volume fraction profiles were measured using optical probes at different bed heights and along two radial directions. Optical probe measurements allow catching a clear hydrodynamic difference between both even and uneven gas distributions. These results were then used to assess CFD simulations with the code Barracuda™ (MP-PIC approach). It is noteworthy that the choice of drag correlation and boundary conditions strongly influences the agreement between the experimental and CFD results. Once the correct parameters are chosen, CFD simulations captured the effect of gas distribution changes.     

Flow structure and bubble dynamics in supercritical water fluidized bed and gas fluidized bed: A comparative study

Supercritical water (SCW) fluidized bed is a new reactor concept for hydrogen production from biomass or coal gasification. In this paper, a comparative study on flow structure and bubble dynamics in a supercritical water fluidized bed and a gas fluidized bed was carried out using the discrete element method (DEM). The results show that supercritical water condition reduces the incipient fluidization velocity, changes regime transitions, i.e. a homogeneous fluidization was observed when the superficial velocity is in the range of the minimum fluidization velocity and minimum bubbling velocity even the solids behave as Geldart B powders in the gas fluidized bed. Bubbling fluidization in the supercritical water fluidized bed was formed after superficial velocity exceeds the minimum bubbling velocity, as in the gas fluidized bed. Bubble is one of the most important features in fluidized bed, which is also the emphasis in this paper. Bubble growth was effectively suppressed in the supercritical water fluidized bed, which resulted in a more uniform flow structure. By analyzing a large number of bubbles, bubble dynamic characteristics such as diameter distribution, frequency, rising path and so on, were obtained. It is found that bubble dynamic characteristics in the supercritical water fluidized bed differ a lot from that in the gas fluidized bed, and there is a better fluidization quality induced by the bubble dynamics in the supercritical water fluidized bed.     

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.     

Projective synchronization of nonlinear-coupled spatiotemporal chaotic systems

The projective synchronization of one-dimensional discrete spatiotemporal chaotic systems is discussed in this paper. The coupling equation is determined by suitably separating the linear term and the nonlinear term of the dynamic function, and two coupled map lattices reach projective synchronization by the nonlinear-coupling method. Besides, this method is expanded to the projective synchronization of the complex network composed by coupled map lattices. Numerical simulations show the effectiveness of the scheme.     

VERIFICATION OF SMOM AND QMOM POPULATION BALANCE MODELING IN CFD CODE USING ANALYTICAL SOLUTIONS FOR BATCH PARTICULATE PROCESSES

For many processes of industrial significance, due to the strong coupling between particle interactions and fluid dynamics, the population balance must be solved as part of a computational fluid dynamics （CFD） simulation. In this work, a CFD based population balance model is tested using a batch crystallization reactor. In this CFD model, the population balance is solved by the standard method of moments （SMOM） and the quadrature method of moments （QMOM）. The results of these simulations are compared to analytical solutions for the population balance in a batch tank where 1） nucleation, 2） growth, 3） aggregation, and 4） breakage are taking place separately. The results of these comparisons show that the first 6 moments of the population balance are accurately predicted for nucleation, growth, aggregation and breakage at all times.     

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.     

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.     

A LBM–DEM solver for fast discrete particle simulation of particle–fluid flows

The lattice Boltzmann method (LBM) for simulating fluid phases was coupled with the discrete element method (DEM) for studying solid phases to formulate a novel solver for fast discrete particle simulation (DPS) of particle–fluid flows. The fluid hydrodynamics was obtained by solving LBM equations instead of solving the Navier–Stokes equation by the finite volume method (FVM). Interparticle and particle–wall collisions were determined by DEM. The new DPS solver was validated by simulating a three-dimensional gas–solid bubbling fluidized bed. The new solver was found to yield results faster than its FVM–DEM counterpart, with the increase in the domain-averaged gas volume fraction. Additionally, the scalability of the LBM–DEM DPS solver was superior to that of the FVM–DEM DPS solver in parallel computing. Thus, the LBM–DEM DPS solver is highly suitable for use in simulating dilute and large-scale particle–fluid flows.     

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 exoerimental data than those of the Gidasoow dra~ model did.     

Fourier analysis of the Eulerian–Lagrangian least squares collocation method

A Fourier analysis was performed in order to study the numerical characteristics of the effective Eulerian–Lagrangian least squares collocation (ELLESCO) method. As applied to the transport equation, ELLESCO requires a C¹-continuous trial space and has two degrees of freedom per node. Two coupled discrete equations are generated for a typical interior node for a one-dimensional problem. Each degree of freedom is expanded separately in a Fourier series and is substituted into the discrete equations to form a homogeneous matrix equation. The required singularity of the system matrix leads to a ‘physical’ amplification factor that characterizes the numerical propagation of the initial conditions and a ‘computational’ one that can affect stability. Unconditional stability for time-stepping weights greater than or equal to 0-5 is demonstrated. With advection only, ELLESCO accurately propagates spatial wavelengths down to 2Δx. As the dimensionless dispersion number becomes large, implicit formulations accurately propagate the phase, but the higher-wave-number components are underdamped. At large dispersion numbers, phase errors combined with underdamping cause oscillations in Crank–Nicolson solutions. These effects lead to limits on the temporal discretization when dispersion is present. Increases in the number of collocation points per element improve the spectral behaviour of ELLESCO.     

Application of CFD-DEM to the study of solid exchange in a dual-leg fluidized bed

The CFD-DEM model was developed to simulate solid exchange behavior between two half beds in a bench-scale two-dimensional dual-leg fluidized bed （DL-FB）. Power spectrum density （PSD） analysis was applied to obtain the dominant frequency （F） of the simulated differential particle number （APLR） between the two half beds. Effects of fluidization velocity （u） and bed material inventory （H） on the solid exchange behavior were studied using the CFD-DEM model. Not only snapshots of the simulated particle flow patterns using the OpenGL code but also the dominant frequency of APLR was similar to the experimental results. The simulation results show that higher fluidization velocity assists the exchange of more particles between the two half beds, but the dispersion of clusters on the bed surface into single particles decreases the cluster exchange frequency. A greater bed material inventory results in more intense cluster exchange. The cluster exchange frequency decreases with an increase of the bed material inventory.