[andrea.pasquali]

Hello,
I'd like to simulate a rotating gear system into oil. I choose the solver interFoam.
U: I compiled new BC to set a ramp of velocity of rotation (timeVaryingRotatingWallVelocity) for my rotating elements (http://www.cfd-online.com/Forums/ope...edvalue.html);
p: I used the totalPressure for the atmosphere patch and zeroGradient for all other walls and rotating element;
k-epsilon model for the turbulence.
First I tried with epsilon = 100 and my simulation died after 0.06 seconds.
Then I increased the value of epsilon = 1000 and my simulation died at 0.84 seconds.
I attached the plots of residuals and Courant number.
You can see that in the end of simulation k and epsilon start increasing very fast and bring the calculation to die in a very little time!

Now I'm thinking to run new analysis increasing the value of epsilon, but before I have some question:
1) Is it correct use k-epsilon for rotating element?
2) If yes, how can calculate the correct value of k and epsilon without find it by tentative? I try to use the formulas that are in OF user guide but I don't know how can use them with rotating elements or if are correct for rotating elements

Thanks anyone for any suggestion

Andrea