Problem with convergence of Kepsilon/LRR model
Hi Foamers,
I am running an analysis looking at the boundary layer of flow over a twodimensional bump. I am currently having trouble getting convergence for LaunderReeceRodi 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 
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 
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 
Then you ll have a problem with your BC i think :)
Can you post a screenshot of that region? Did you try Code:

2 Attachment(s)
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!

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

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

All times are GMT 4. The time now is 10:32. 