[mwaqas]

Hello Everyone

I am working on twoPhaseEulerFoam fluidised bed 2D OpenFoam 2.3.0. I am trying to validate my setting for twoPhaseEulerFoam (fluidised bed) by reproducing already existing results . I am using same parameters as mentioned in papers (refrence papers are mentioned at the end) (didn't change fvScheme).
viscosityModel Gidaspow;
conductivityModel Gidaspow;
granularPressureModel Lun;
frictionalStressModel JohnsonJackson;
radialModel SinclairJackson.

Bed width, m 0.28
Bed height, m 1.0
Bed depth, m 0.025
Initial bed height, m 0.4
Initial solid packing 0.6
Superficial gas velocity, m/s 0.38, 0.46
Min. fluidization velocity, m/s 0.065
Gas density, kg/m 3 1.2
Gas kinematic viscosity, m 2 /s 1.4*10-5
Particle diameter, m 2.8*10-4
Particle density, kg/m 3 2500
Particle–particle restitution coefficient 0.9
Specularity coefficient 0.5
alpha max 0.65
alpha min 0.50

I tried both equilibrium on and true. I am facing following problems in my simulation

1) I am getting same results for pressure drop and bed expansion ratio, but in my simulation I am getting max. solid volume fraction up to 0.48, while in these papers they are getting up to 0.55

2) In my simulation time-averaged particles velocity is up to 0.12 m/s (max. value, I am using temporal statistics filter in paraview for time-averaged calculation) while in reference papers velocity is of the order of 0.8 m/s

3)when I use "patchIntegrate phi inlet/outlet" command in terminal window, both mass flow rate show a difference of 5-6 %. Although patchIntegrate phi.Particles shows no solid mass loss, as I am using 0/U.particles outlet= fixedValue 0. Then why it is showing this mass difference.

4)In my simulation, the U.particles velocity in the grids where no solid exists is not equal to zero.

It would be really nice and helpful If could guide me regarding these problems.

1) http://www.sciencedirect.com/science...09250905004653
2) http://www.sciencedirect.com/science...09250905004653

Best Regards
mwaqas