Help with Kepsilon for 2D turbulence modelling
1 Attachment(s)
Hi,
I'm simulating a flow though an volute casing. I'm in the initial stages of showing grid independence. the kepsilon model was working fine until the no. of elements was around 2.5 lakhs. in the 4 lakh element model im not getting a converged solution. the residuals for continuity are not stable now. I'm using the simple algorithm, with 2nd order upwind diff(used 1st order too; but same problem) 0.5 relaxation. its a singlephase flow. Please suggest ways to converge this solution. 
Re : Help with Kepsilon for 2D turbulence modelling
Hi aditya.pandare
It seems to me that there are two possibilities : numerical diffusion problem or wrong nearwall treatment. For the case of nearwall treatment, as you probably know, when you fine you mesh near wall, you decrease the y_plus value of your closest cell to wall. If y_plus values had decreased under 5 or 6, you need to use enhanced wall functions to have solution for viscous sublayer. I would think that there would be numerical diffusion problem if you had not the same issue with a higher order discretization scheme. Orkun 
another possiblity:
change from single precesion to double precesion 
Hi orkun and ztdep! thnks for your reply!
I am getting a y+ value of around 3 to 4. So I am running a simulation using enhanced wall treatment. I will also run a double precision solution after this run! another question, how do we find the value of velocity gradient at the wall? since the velocity profile equations themselves are functions of friction velocity (which in turn is viscosity multiplied by velocity gradient at wall) thanks again 
You can use wall shear stress to calculate velocity gradient, or you may calculate it by a UDF.

even After using the enhanced wall treatment, the solution isnt converging. I am checking convergence on the basis of areaweighted avg along a line.
Is it because i'm using an "outflow" boundary condition? would a pressure outlet be more suitable? 
All times are GMT 4. The time now is 17:25. 