Fluidized Bed Simulation
Theory and Physics
Overview
Professor, what exactly does fluidized bed simulation do?
A fluidized bed is a device where gas is blown from below into a particle-packed bed, causing the particles to suspend and mix. It's a core chemical engineering technology used in petroleum refining FCC (Fluid Catalytic Cracking), coal gasification, biomass combustion, pharmaceutical granulation coating, and more. CFD is used to predict internal particle behavior and gas mixing.
What methods are there for CFD of fluidized beds?
Governing Equations
Please explain the equations for TFM (Two-Fluid Model).
It solves the continuity and momentum equations for both the gas phase and the solid phase. KTGF (Kinetic Theory of Granular Flow) is used for solid phase stress.
Solid phase pressure is derived from Granular Temperature $\Theta_s$.
What is the key parameter for fluidization?
It's the minimum fluidization velocity $U_{mf}$. Particles begin to suspend when gas velocity exceeds this. It can be estimated from the Ergun equation.
In the fluidized state, pressure loss balances with bed weight. $\Delta p = (1-\varepsilon_{mf})(\rho_s - \rho_g) g L$ is the criterion for fluidization.
Geldart Classification
Does the fluidization behavior differ depending on particle type?
The basic classification is Geldart's (1973).
| Group | Particle Size | Fluidization Characteristics | Example |
|---|---|---|---|
| A | 20~100 μm | Uniform expansion followed by bubble formation | FCC catalyst |
| B | 100~1000 μm | Direct bubble fluidization | Sand, glass beads |
| C | < 20 μm | Strong cohesiveness, difficult to fluidize | Wheat flour, talc |
| D | > 1000 μm | Spout formation | Grains, coal lumps |
The Discovery of Fluidization—The Eve of the FCC Process and the Fluidized Bed Revolution
The industrial application of Fluidized Bed technology expanded rapidly in the 1940s, triggered by Standard Oil (now ExxonMobil) developing the Fluid Catalytic Cracking (FCC) process. This phenomenon, where sand particles levitate with air and behave "like a liquid," was said to look like magic to chemical engineers in the early 20th century. The Ergun equation (1952), the fundamental theory of fluidization, remains at the core of fluidized bed design today, semi-empirically linking ε (void fraction) and ΔP (pressure loss). CFD simulation of fluidized beds heavily depends on how this Ergun model represents inter-particle forces.
Physical Meaning of Each Term
- Temporal Term $\partial(\rho\phi)/\partial t$: Imagine the moment you turn on a faucet. At first, water comes out spluttering and unstable, but after a while, the flow becomes steady, right? This "period of change" is described by the temporal term. The pulsation of blood flow with a heartbeat, or the flow fluctuation each time an engine valve opens and closes—all are unsteady phenomena. So what is steady-state analysis? It looks only at "after sufficient time has passed and the flow has settled down"—meaning setting this term to zero. Since computational cost drops significantly, solving first with a steady-state approach is a basic CFD strategy.
- Convection Term $\nabla \cdot (\rho \mathbf{u} \phi)$: What happens if you drop a leaf into a river? It gets carried downstream by the flow, right? This is "convection"—the effect where fluid motion transports things. Warm air from a heater reaching the far corner of a room is also because the "carrier," air, transports heat via convection. Here's the interesting part—this term contains "velocity × velocity," making it nonlinear. That is, as flow speed increases, this term rapidly strengthens, making control difficult. This is the root cause of turbulence. A common misconception: "Convection and conduction are similar" → They're completely different! Convection is carried by flow, conduction is transmitted by molecules. There's an order of magnitude difference in efficiency.
- Diffusion Term $\nabla \cdot (\Gamma \nabla \phi)$: Have you ever put milk in coffee and left it? Even without stirring, after a while it naturally mixes. That's molecular diffusion. Now a question—honey and water, which flows more easily? Obviously water, right? Honey has high viscosity ($\mu$), so it flows poorly. Higher viscosity strengthens the diffusion term, making the fluid move in a "sluggish" manner. In low Reynolds number flows (slow, viscous), diffusion is dominant. Conversely, in high Re number flows, convection overwhelms and diffusion plays a supporting role.
- Pressure Term $-\nabla p$: When you push a syringe plunger, liquid shoots out forcefully from the needle tip, right? Why? Because the plunger side is high pressure, the needle tip is low pressure—this pressure difference becomes the force pushing the fluid. Dam discharge works on the same principle. On a weather map, where isobars are tightly packed? That's right, strong winds blow. "Flow is generated where there is a pressure difference"—this is the physical meaning of the pressure term in the Navier-Stokes equations. A point of confusion here: "Pressure" in CFD is often gauge pressure, not absolute pressure. If results go wrong immediately after switching to compressible analysis, it might be due to confusing absolute/gauge pressure.
- Source Term $S_\phi$: Heated air rises—why? Because it becomes lighter (lower density) than its surroundings, so buoyancy pushes it up. This buoyancy is added to the equation as a source term. Other examples: chemical reaction heat generated by a gas stove flame, Lorentz force applied to molten metal by an electromagnetic pump in a factory... These are all actions that "inject energy or force into the fluid from the outside," expressed by the source term. What happens if you forget the source term? In natural convection analysis, forgetting buoyancy means the fluid doesn't move at all—a physically impossible result where warm air doesn't rise in a heated room in winter.
Assumptions and Applicability Limits
- Continuum Assumption: Valid for Knudsen number Kn < 0.01 (mean free path ≪ characteristic length)
- Newtonian Fluid Assumption: Linear relationship between shear stress and strain rate (viscosity model needed for non-Newtonian fluids)
- Incompressibility Assumption (for Ma < 0.3): Treat density as constant. Consider compressibility effects for Mach number ≥ 0.3
- Boussinesq Approximation (Natural Convection): Consider density change only in the buoyancy term, use constant density in other terms
- Non-applicable Cases: Rarefied gas (Kn > 0.1), supersonic/hypersonic flow (shock capturing required), free surface flow (VOF/Level Set etc. required)
Dimensional Analysis and Unit Systems
| Variable | SI Unit | Notes / Conversion Memo |
|---|---|---|
| Velocity $u$ | m/s | When converting from volumetric flow rate for inlet conditions, pay attention to cross-sectional area units |
| Pressure $p$ | Pa | Distinguish between gauge and absolute pressure. Use absolute pressure for compressible analysis |
| Density $\rho$ | kg/m³ | Air: ~1.225 kg/m³@20°C, Water: ~998 kg/m³@20°C |
| Viscosity Coefficient $\mu$ | Pa·s | Be careful not to confuse with kinematic viscosity $\nu = \mu/\rho$ [m²/s] |
| Reynolds Number $Re$ | Dimensionless | $Re = \rho u L / \mu$. Indicator for Laminar/turbulent transition |
| CFL Number | Dimensionless | $CFL = u \Delta t / \Delta x$. Directly related to time step stability |
Numerical Methods and Implementation
Details of Numerical Methods
Please tell me the key numerical points for fluidized bed CFD.
Fluidized bed simulation with TFM (Eulerian Granular Model) has several unique challenges.
Mesh and Mesoscale Structure
In fluidized beds, dense particle structures called "clusters" are important. Cluster size is about 10~100 times the particle diameter, so the mesh needs to be sufficiently fine to resolve this.
| Mesh | Resolution | Computational Cost | Accuracy |
|---|---|---|---|
| Fine | $\Delta x \approx 5 d_p$ | Very High | High |
| Standard | $\Delta x \approx 10$~$20 d_p$ | Moderate | Good |
| Coarse + Filter | $\Delta x > 50 d_p$ | Low | Requires filter model |
What happens if clusters cannot be resolved with a coarse mesh?
It overestimates bed expansion and underestimates gas bypass. In other words, it appears more uniformly fluidized than reality. You need to correct with Filtered TFM (Igci et al., 2008; Ozel et al., 2013) or use a sufficiently fine mesh.
Drag Model Selection
The most important closure model in fluidized beds is the gas-solid interphase drag.
| Model | Characteristics | Recommended Use |
|---|---|---|
| Gidaspow | Switches between Ergun + Wen-Yu | BFB (Bubbling Fluidized Bed) standard |
| Syamlal-O'Brien | Continuous, adjustable parameters | General purpose |
| EMMS | Considers mesoscale structure | CFB (Circulating Fluidized Bed) |
| Koch-Hill | LBM database | High accuracy |
What is the EMMS model?
The Energy Minimization Multi-Scale (EMMS) model, proposed by Li & Kwauk (Chinese Academy of Sciences), reflects the gas bypass effect due to cluster structure in the drag force. It can incorporate some cluster influence even with coarse meshes, so it's widely used in industrial-scale circulating fluidized beds.
Time Step and Computation Time
Fluidized bed TFM requires unsteady calculation, needing computation for several to tens of seconds of physical time.
| Parameter | Recommended Value | Remarks |
|---|---|---|
| $\Delta t$ | $10^{-4}$~$10^{-3}$ s | Courant Number < 0.5 |
| Physical Time | 5~30 s | Until statistical steady state is reached |
| Averaging Start | After 2~5 s | Exclude initial transient |
TFM vs DEM-CFD—The Two Major Trends in Fluidized Bed Simulation
Fluidized bed CFD broadly has two approaches: Two-Fluid Model (TFM/Euler-Euler) and DEM-CFD (Euler-Lagrange). TFM treats particles as a continuum, making it scalable to systems with over a million particles, but individual particle contacts are averaged out and lost. DEM-CFD tracks individual particles, surpassing TFM in physical accuracy, but computational cost increases sharply when particle count exceeds 100,000. Full-scale CFD of industrial-scale fluidized beds (diameter 3 m × height 10 m) is still realistically handled by TFM even in the 2020s, with DEM-CFD playing a role in validation and closure model development.
Upwind Scheme (Upwind)
1st Order Upwind: Large numerical diffusion but stable. 2nd Order Upwind: Improved accuracy but risk of oscillation. Essential for high Reynolds number flows.
Central Differencing (Central Differencing)
2nd order accuracy, but numerical oscillation occurs for Pe number > 2. Suitable for low Reynolds number diffusion-dominated flows.
TVD Scheme (MUSCL, QUICK, etc.)
Suppresses numerical oscillation while maintaining high accuracy using limiter functions. Effective for capturing shocks and steep gradients.
Finite Volume Method vs Finite Element Method
FVM: Naturally satisfies conservation laws. Mainstream in CFD.FEM: Advantageous for complex shapes and multiphysics. Mesh-free methods like SPH are also developing.
CFL Condition (Courant Number)
Explicit method: CFL ≤ 1 is the stability condition. Implicit method: Stable even for CFL > 1, but affects accuracy and iteration count.LES: CFL ≈ 1 recommended. Physical meaning: Information should not travel more than one cell per time step.
Residual Monitoring
Convergence is judged when residuals for Continuity, momentum, and energy equations drop by 3~4 orders of magnitude. Mass conservation residual is particularly important.
Relaxation Factor
Pressure: 0.2~0.3, Velocity: 0.5~0.7 are typical initial values. Reduce if diverging. Increase after convergence to accelerate.
Internal Iteration for Unsteady Calculation
Iterate within each time step until steady solution converges. Internal iteration count: 5~20 times as a guideline. If residuals fluctuate between time steps, review the time step size.
Analogy for SIMPLE Method
The SIMPLE method is an "alternating adjustment" technique. First, velocity is tentatively determined (predictor step), then pressure is corrected so that mass conservation is satisfied with that velocity (corrector step), and velocity is revised with the corrected pressure—this back-and-forth is repeated to approach the correct solution. It resembles two people leveling a shelf: one adjusts the height, the other balances it, and they repeat this alternately.
Analogy for Upwind Scheme
The upwind scheme is a method that "stands in the river flow and prioritizes upstream information." A person in the river cannot tell the source of the water by looking downstream—it's a discretization method reflecting the physics that upstream information determines downstream. Accuracy is 1st order, but it's highly stable because it correctly captures flow direction.
Practical Guide
Practical Guide
Related Topics
なった
詳しく
報告