[bgraves]

Hi Foamers,
I am running an analysis looking at the boundary layer of flow over a two-dimensional bump. I am currently having trouble getting convergence for Launder-Reece-Rodi model, after about a day of running the residuals blow up. I am currently using the following Boundary Conditions:

Inlet
U fixedValue
P fixedValue
k/epsilon fixedValue
nut calculated

Outlet
U zeroGradient
P zeroGradient
k/epsilon zeroGradient
nut zeroGradient

Wall
U fixed Value 0 0 0
P zeroGradient
k kqRWallFunction
epsilon epsilonWallFunction
nut nutkWallFunction

I tried to use the method found on this site for calculating the value for Reynolds stress, and got that far. If this helps my colleague ran the same analysis using StarCCM+ and his simulation worked, and his results were good.

If anyone has any suggestions of how to improve my simulations it would be greatly appreciated.

Thanks,
Ben

[Tobi]

Hi,

did you have a look at the last timesteps. It should be possible to see the Problem zones there. Maybe you have a Problem in a cell or at a BC (outlet for example).

Tobi

[bgraves]

Hi Tobi,

Thanks for the speedy reply! I did check, by looking in paraFoam and the velocity goes to infinity by the 20th iteration.

Ben

[Tobi]

Then you ll have a problem with your BC i think

Can you post a screenshot of that region?

Did you try

inletOutlet 
pressureInletOutletVelocity

Tobi

[bgraves]

That did get the code to run however now my residuals are terrible, and the results are bogus. The nice thing is that now the code runs the problem must be in the solving methods right? Since I know the grid is good. Attached is the system and 0 file for my problem if you wouldnt mind looking at it when you have a chance I would greatly appreciate it. I would like to say thank you for helping me!

[Tobi]

would it be possible to upload or send me the whole case?

[bgraves]

I could try to send you the file the constants folder is to big to upload to this website. where should I send it