3.2. Fixed Bed (Ergun Correlation)

3.2.1. Case Description

In this verification case, the pressure drop along a fixed bed of spherical particles is evaluated using the Ergun[5] dense flow drag law, which is modeled using the 2-Way Fluent (unresolved) CFD coupling method. For more information about how CFD models and methods are implemented in Rocky, refer to the Rocky DEM-CFD Coupling Technical Manual. The case illustration is shown in Figure 3.9: Fixed bed with 0.5 m height of spherical particles within a container..

Figure 3.9: Fixed bed with 0.5 m height of spherical particles within a container.

Fixed bed with 0.5 m height of spherical particles within a container.

3.2.2. Analytical Formulation

The Ergun [5] correlation for the pressure losses accompanying the flow of fluids through a fixed bed of granular material can be expressed as follows:

(3–3)

where:

  • is the pressure drop suffered by the fluid that flows though the fixed bed of granular material.

  • L is the fixed bed's height.

  • is the solid volume fraction within the fixed bed.

  • is the fluid volume fraction within the fixed bed.

  • us is the superficial relative velocity of the fluid in the direction of the flow.

  • is the fluid dynamic viscosity.

  • is the fluid density.

  • dp is the diameter of the spherical particles within the bed.

3.2.3. Input Data & Setup

The equations shown in the last section can be resolved and equivalent results can be calculated by Rocky considering the same input data and boundary conditions. The input parameters for this verification case setup are presented in Table 3.5: Verification case input parameters..

Table 3.5: Verification case input parameters.

Parameter

Value

Unit

Physical Model:

Gravity Y-9.81

Geometry (Container):

Length (X)0.1m
Length (Y)0.8m
Length (Z)0.1m

Solid Properties (Geometry):

Material Density7850
Material Young's Modulus 1011
Material Poisson's Ratio0.3-

Solid Properties (Particle):

Bed Height 0.5m
ShapeSphere-
Diameter3.95mm
RotationsEnabled-
Input TypeVolume Fill-
Input Mass10kg
Material Bulk Density500
Material Young's Modulus108
Material Possions's Ratio0.3-

Fluid Properties (Fluent):

Domain Length (X)0.1m
Domain Length (Y)0.8m
Domain Length (Z)0.1m
Fluent Mesh Size10mm
Inlet Velocity0.275m/s
Outlet Gauge Pressure0Pa
Walls Shear ConditionFree Slip-
Density1000
Dynamic Viscosity0.001
2-Way Fluent Properties:
Drag LawErgun[5]-
Solver Parameters:
Simulation Duration1.0s

3.2.4. Results

After running the Rocky case as specified, the results can then be compared to the analytical values. Here, the static pressure loss along the bed height is shown in Figure 3.10: Rocky simulation results using 2-Way Fluent coupling showing the static pressure along the fixed bed at the last timestep. and Figure 3.11: Comparison of static pressure [kPa] for the Rocky simulation and reference results.. The numerical solution given by Rocky presents strongly correlated values to those obtained by the analytical expression.

Figure 3.10: Rocky simulation results using 2-Way Fluent coupling showing the static pressure along the fixed bed at the last timestep.

Rocky simulation results using 2-Way Fluent coupling showing the static pressure along the fixed bed at the last timestep.

Figure 3.11: Comparison of static pressure [kPa] for the Rocky simulation and reference results.

Comparison of static pressure [kPa] for the Rocky simulation and reference results.

The absolute and relative errors for the pressure drop are compared in Figure 3.12: Absolute and relative errors.. The maximum absolute error for the pressure drop error is around 1.0 kPa. The maximum relative error is around 2.3%.

Figure 3.12: Absolute and relative errors.

Absolute and relative errors.

Table 3.6: Pressure target lists the value for the total pressure drop on the fixed bed. This includes the target value calculated by the analytical expression as compared to the value calculated by Rocky. A Ratio of 1.00 shows a strong correlation between the results.

Table 3.6: Pressure target

Target [kPa]Rocky [kPa]Ratio[-]
117.42116.900.99