[Stylianos]

Hello Everyone

I'm trying to run an axial turbine simulation via MRFSimpleFoam (Rotor-Stator with GGI interface), i've tried both kEpsilon and kOmegaSST turb-models but i get either bounding omega (in case of the kOmega) or bounding epsilon (in case of kEpsilon) with negative values dose anyone have any idea about that?

i use:
inlet:
U: flowRateInletVelocity
p: zeroGradient
k: turbulentIntensityKineticEnergyInlet
epsilon: turbulentMixingLengthDissipationRateInlet
omega: turbulentMixingLengthFrequencyInlet

outlet:
U: inletOutlet
p: fixedMeanValue
k/epsilon/omega: inletOutlet

with turbulentIntensity = 5%
and MixingLength = 0.079 (1/5*(smaller diameter of my duct))

Stelios

[Stylianos]

Managed to get kOmegaSST running,
Had to initialize omega field with omega = epsilon/(Cm*k) and not epsilon/k as mentioned in several threads.

no luck with kEpsilon thought.

Hope this is useful for someone,
Stelios

[mrshb4]

Dear Stelios
I have the same problem with my case that is a centrifugal fan with a duct, my boundary condition is almost similar to you. First, my model would diverge at about 10 iterations.(k and epsilon would blow up). I changed my fvSchemes and fvSolutions dictionaries and got better answers. At least with no divergence until about 200 iterations, but after that it would diverge again.

Could you solve your problem?
It you could, tell me what you did.

Good luck
Mohamadreza

[Stylianos]

Hello Mohamad

The only think i managed to do to get a solution was to use kOmegaSST turbulence model instead of kEpsilon. But that was some time ago and with OF 1.5-dev i don't know if there is something else going on in 1.6.

Stelios

[kaka22]

Hi Stelios,

I am using kOmegaSST turbulence model in rhoPisoFoam. I am getting bounding omega with negative values. I have also given initial value of omega, still bounding omega is appearing. My case is not crashing but it's also not showing the right flow profile. I'm guessing the problem lies with the bounding omega. Is there any suggestion that bounding omega can be eliminated by changing any of the parameters? Thanks!

Novi

[farahiam]

Hi Kaka,
Have you solved the problem?
I am also facing similar problem as yours..

Ayu

[cleoo]

Originally Posted by farahiam

Hi Kaka,
Have you solved the problem?
I am also facing similar problem as yours..

same, has anyone found a solution, my bounding omega values get really high and then the solution diverges

[hakolekar]

Any solutions ???

[ancolli]

try with non-orthogonality limiters. Control limiter based on boundedness messages on scalars

http://www.tfd.chalmers.se/~hani/kur...stPractice.pdf