Initial Residual for p too high!
1 Attachment(s)
Hi all,
I am simulating a simple 2D flat plate flow in OpenFoam using simpleFOAM to validate a new turbulence model that I would like to implement here. This new turbulence model is a slight modification of the kOmegaSST model which includes few new terms into the nut equation. Results from kOmegaSST model are all OK! However, when I use this new model for the same grid, my initial pressure residuals oscillate around a pretty high value, at around 0.03 while the initial U,k and omega residuals are all at a reasonable convergence criterion, at around 10^5. This is not making much sense to me. Because, after some iterations(about 6~7,000) when I extract the results (even without reaching the convergence tolerance) and compare, they don't yet seem so weird or deviating highly from that of the SST's results. I have checked my code numerous times, doesn't seem to have any problems in there. Could there be a problem with the solvers I chose? I am as well wondering how the presidual is calculated in simpleFOAM. I would highly appreciate any insights into this! Thanks! Nikesh These are my settings and I've used the same for kOmegaSST (in the pic attached). 
hi
Test for one or two order lower p tolernces than U,... 1e9 or 1e10 
Thanks!
Tried, yet not much of a difference. 
Hi
Basically I am having the same problem.... I just want to ask you if you managed to do it in the end. I am simulating the flow over an aerofoil using k omega sst and simpleFoam with wall functions.......... Have you got any ideas? Thanks Andrei 
Hii Andrei,
Well, I am still stuck with the same problem. Obviously higher presidual means mass is not being conserved so well somewhere in the cells. You might want to look into your boundary conditions once more. Basically that is what I am trying to do too. And the schemes and type of mesh you are using for airfoil flow. 
Hi
Ok I understand . I am using k omega sst , wall functions and simpleFoam solver ..... I have changed my boundary conditions many times still nothing......for example how you o file looks like and something that I don't understand : how do you calculate the turbulent kinetic energy K and the rate of dissipation W(omega)....... thanks Andrei 
Nick,
There are some strategies that I would try to get these residuals down to something you would like.
There are a lot of other strategies that you can try, but this one might be sufficient. Without knowing more details like divergence schemes; mesh structure and checkmesh results; y+ values; and boundary conditions there is not much to add on my part. I would look at the thread http://www.cfdonline.com/Forums/ope...problems.html and then move on to a search of the forum. There are many threads about simpleFoam convergence...but my list of threads is not in front of me right now. good luck. 
how can decrease the initial residuals for p in unsteady problems when relaxations are not applicable?

All times are GMT 4. The time now is 02:00. 