RSM convergence
Hi all,
I need some advice... I'm running RANS (with the LRR Reynold Stress Turbulence Model) simulations for atmospheric flow over a regular array of simple buildings. I'm aware that RSM models are notoriously tricky to get converged but, unfortunately, I must use it (or I can run LES, but I'd prefer not to step in that field... yet). The problem is that it seems like I'm not able to reach complete convergence, at least basing this on the (initial) residuals only. While I can consider Ux residuals as ~ converged (they go flats at ~ 1e-04), I cannot say the same for other variables (notably p, whose residuals remains ~0.005). I tried with kEpsilon and I was able to reach convergence (residuals <1e-04). Searching in this forum I've discovered how to, at least, start with a converging RSM simulation by initialise it with a partially converged kEpsilon simulation. Unfortunately this doesn't help after a good number of iterations, when the residuals start becoming flat. Initially I tried with simpleFoam (steady state), but soon I discovered that the RSM model is able to catch some unsteady features of the flow (in particular vortex shedding) and so I couldn't expect it to converge with steady state simulations. I'm now running pisoFoam and I'm using second order discretisation schemes. The mesh I'm using is a fine Cartesian mesh with low expansion ratios. I'm going to try further refinements in some parts of the domain, but I'm running out of tricks, now. Is anyone able to suggest me anything else? Is there anything I'm doing wrong in your opinion? Please ask me if you need more information on my case study. Thanks in advance for your help! Matteo |
small improvements
2 Attachment(s)
Hello,
after some useful advice and more tweaking of the boundary conditions I was able to improve the results. I'm now able to reach values of residuals for almost all variables lower than 1e-05 with first order discretisation (I need second order, though...). The only exception is the p residual which does not go below 0.002. Does anyone have any idea? For info, these are my 'system' files: controlDict: Code:
application pisoFoam; Code:
ddtSchemes Code:
solvers Matteo |
Hello mcarpe,
I'd like to know: How have you determine the relaxation factors. I'd like to learn that. Cheers |
Quote:
They are in my fvSolution file because of some previous attempts with simpleFoam. Anyway I'm afraid there isn't any simple rule for determining the UR factors. Best Matteo |
Hello mcarpe,
in interFoam they work. I had problems i my case, but I used your RF's. and I have convergence now. However, I had to modify them to 0.4 0.4 0.4 and so on. kinda regards! |
Convergence problems Fluent RSM
Hello!
I am doing a convergence analysis of a single tube with heat flux as boundary condition in ANSYS Fluent with RSM solver and scalable as wall function. I'm studing the surface temperature for this analysis. When I refine my mesh, and scalable wall function limit the y+ value into 11.25 for over refine mesh, the temperature change and it trend doesn't converge. I don't know why if the temperature is in the same position change when during the refinament and y* and the boundary condition are equal for both case. Any of you have any suggestion to achieve a convergence result?. Thank so much! Rafa |
All times are GMT -4. The time now is 19:48. |