Bubble Column Simulation: Different Turbulence Models different results
 

I'm simulating a Bubble Column; 0.15m by 0.9m, Single Point Sparger, hole diameter 1.96e-3m. I tried the simulation with three different turbulence models; K-eps, RSM and LES.

Superficial gas velocity 20mmps.

Turbulence Nature of the Results/Error
Model

K-eps = The bubble shoots right to the top and it appears as if jet has been sent through.

RSM = Fatal error in linear solver, the Water.turbulence eddy viscosity frequency goes out of bounds.

LES = The simulation runs well and the bubble rises as expected.

Can anyone please help me understand the probable errors in the simulation and the reason for their occurrence.

Thanks

* There are lots of multiphase parameters you should consider, like dispersion forces and things like that. That might get the simple turbulence model runs working.
* Try the SST turbulence model. It is the best all-round turbulence model these days, it is better than k-e.
* RSM models are always very hard to get to converge. Make sure your mesh is really good quality, your time step small enough and your initial condition suitable. But even then RSM is hard to get working.
* If LES is working well then why not just stick with that one?

Thanks ghorrocks....I know its a late reply...but I finally got the simulations to run without crashing. RSM was tough to run but reducing the time step did help. The RANS model hasnt given good profiles compared with the LES model.

These are the profiles I have obtained for Multipoint sparger incorporating the K eps and RANS models for simulation. Bubble Size 5mm, column dimensions 0.15m by 1m.

Any ..absolutely any comments would dearly be appreciated

Thanks

This looks typical for RSM - RSM predicts additional vorticies and lots of little features, and k-e predicts a single big vortex with something approaching solid body rotation.

What was the rationale for the simulation, if I may ask? If it is just about the traces of bubble, then did you use Langrangian particle tracking, or it was for bubble size distribution etc using population balance model (MUSIG etc)? Which forces did you enable on bubbles? From the picture it looks like Drag and Lift (since the diameter is small) are a must. You may consider dispersion force as well, if you are lucky enough to have economical convergence...

EDITED>>
Adaptive timestepping really helps, since typically it is difficult to estimate the right timescales... I am quite surprised at your selection of RSM/LES though, considering the quality of mesh these models demand. Throw in the complexity of multiphase physics and its a nightmare in terms of cost and possibility of convergence..