|
[Sponsors] |
k-omega-SST breakdown with high resolution mesh |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 27, 2011, 01:58 |
k-omega-SST breakdown with high resolution mesh
|
#1 |
Senior Member
Roman Thiele
Join Date: Aug 2009
Location: Eindhoven, NL
Posts: 374
Rep Power: 21 |
Lately I ran into problems using buoyantBoussinesqSimpleFoam together with the k-omega-SST model by menter when I increase the mesh resolution in order to find mesh independency.
when the cell count gets too large the whole simulation breaks down, it does not converge, and there is neither heat transfer in the annulus nor does the velocity field fully develop. the curious thing about it is, that with higher inlet velocities, this cell count can be higher, has anyone an idea, why this could be and/or if there are limitation in the usage of kOmegaSST?
__________________
~roman |
|
May 5, 2011, 05:20 |
updated bug report for this
|
#2 |
Senior Member
Roman Thiele
Join Date: Aug 2009
Location: Eindhoven, NL
Posts: 374
Rep Power: 21 |
maybe somebody has an idea. here is also the updated bug report:
http://www.openfoam.com/mantisbt/view.php?id=179 the solution is in the bug report. but for clarity, if someone wants to run the kOmegaSST in OpenFOAM in low Re mode, the following should be considered: http://www.openfoam.com/mantisbt/view.php?id=179#c351 quote from henry: Note that the k-omega SST model we provide is in high-Re form and does not include the wall-damping terms often included in the k-omega model for near-wall and low-Re flow. However, you can still use the k-omega SST model for low-Re and near wall flow for a range of resolutions if you use a continuous wall-function (which in OpenFOAM-1.7.x is named nutSpalartAllmarasWallFunction for historical reasons) and this should be used as the wall BC in nut. The BC of k for the continuous wall-function should be kqRWallFunction. If these changes do not help it may be worth investigating the viscosity averaging in omegaWallFunctionFvPatchScalarField: scalar omegaVis = 6.0*nuw[faceI]/(beta1_*sqr(y[faceI])); scalar omegaLog = sqrt(k[faceCellI])/(Cmu25*kappa_*y[faceI]); omega[faceCellI] = sqrt(sqr(omegaVis) + sqr(omegaLog)); we have found cases for which this causes a sudden change in the viscosity near the wall if the mesh is sufficiently fine and that just using the logarithmic part give more continuous behavior: omega[faceCellI] = omegaLog;
__________________
~roman Last edited by romant; May 5, 2011 at 06:56. Reason: change in bug status |
|
January 13, 2012, 14:20 |
|
#3 |
Member
Timo K.
Join Date: Feb 2010
Location: University of Stuttgart
Posts: 66
Rep Power: 16 |
Hallo Roman,
do you still have this problem? I also encountered problems with convergence with finer meshes for the kOmegaSST model. The problem comes, since the modification of nut with sqrt(2) was introduced about one year ago. nut_ = a1_*k_/max(a1_*omega_, F2()*sqrt(2.0)*mag(symm(fvc::grad(U_)))); Without this factor sqrt(2) there is no problem with convergence, but the results are poorer compared to CFX. An unevaluated thesis of me is, that this factor is additionally not used in some areas. And/or there should be a slightly modification in the omega-equation (only in some areas), instead of: 2*gamma*S2 --> gamma*omega/k*P (or similar, sorry it is friday evening) as you can see in: Development and Application of SST-SAS Turbulence Model in the DESIDER Project, of Egorov and Menter Are there any comments on that? |
|
February 19, 2012, 08:14 |
|
#4 |
Senior Member
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 368
Rep Power: 20 |
Hi all,
I have opened a tread about this issue here: http://www.cfd-online.com/Forums/ope...estigated.html Regards V. |
|
July 3, 2012, 16:21 |
|
#5 |
Member
Jubayer
Join Date: Oct 2009
Location: The University of Western Ontario, London, Ontario
Posts: 42
Blog Entries: 1
Rep Power: 16 |
Hi Roman,
Have you resolved the issue with this high resolution mesh with buoyantBoussinesqSimpleFoam solver and k-omega SST model? I am facing a similar problem. For my wind loading on solar panel problem, I have used pisoFoam solver k-omega SST without any problem. However, for my heat transfer problem, I increased the mesh resolution close to the surface of the panel (y+<5) and used buoyantBoussinesqPimpleFoam solver. My simulation is blowing up (large negative bounding values for k and omega). I have tried various schemes in fvschemes but still could not find a solution. What else should I consider for making my simulation converge? Jubayer |
|
July 4, 2012, 02:44 |
|
#6 |
Senior Member
Roman Thiele
Join Date: Aug 2009
Location: Eindhoven, NL
Posts: 374
Rep Power: 21 |
Hej Jubayer,
I have found a solution that is more or less connected to the smoothness of the mesh. If the cell growth is too larger close to the wall, you will end up with a non-physical solution and/or blowing up solution. I made sure that the cell growth is only 1.1 close to the wall and going towards the core one could increase the cell growth. In snappyHexMesh this is solved by actually using the right amount of boundary layers with a slow growth. In blockMesh, one would need to use different block towards the core of the flow and use a block with very smooth grading towards the wall.
__________________
~roman |
|
Tags |
heat transfer, komegasst, turbulence |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[ICEM] surface mesh merging problem | everest | ANSYS Meshing & Geometry | 44 | April 14, 2016 06:41 |
[ICEM] Negative volume error in hybrid mesh | siw | ANSYS Meshing & Geometry | 4 | September 3, 2014 05:25 |
[ICEM] Unstructure Meshing Around Imported Plot3D Structured Mesh ICEM | kawamatt2 | ANSYS Meshing & Geometry | 17 | December 20, 2011 11:45 |
larger blend factor or high resolution ? | amine | CFX | 2 | March 4, 2008 20:28 |
upwind vs high resolution | Richard | CFX | 3 | February 20, 2006 06:35 |