View Issue Details
ID | Project | Category | View Status | Date Submitted | Last Update |
---|---|---|---|---|---|
0000179 | OpenFOAM | Bug | public | 2011-04-08 15:48 | 2014-02-10 13:16 |
Reporter | Assigned To | henry | |||
Priority | normal | Severity | major | Reproducibility | always |
Status | closed | Resolution | no change required | ||
Platform | AMD64 | OS | Linux Ubuntu | OS Version | 10.10 |
Summary | 0000179: k-omega sst model breaks down with high resolution grid | ||||
Description | The k-omega model does break down and does not produce predictions at all if the grid resolution goes up. | ||||
Steps To Reproduce | increase the resolution for a k-omega low-Reynolds number simulation, but keep the cell closest to the wall the same. all other cells should not be smaller than the cell closest to the wall. | ||||
Additional Information | the cases in the attached file only need to be run with buoyantBoussinesqSimpleFoam | ||||
Tags | No tags attached. | ||||
|
the case files for this can be found under http://www.reactor.sci.kth.se/files/breakdownKOmega.zip it is a simulation of a heated rod submersed in liquid metal. the crossectional results for k, omega, p_rgh and p can be found in the figures within the zip file. the velocity distribution and temperature distribution are found in radialUdist and radialTempDist within the cases respectively. they can be plotted with any plotting program that supports text file input. the first column is the radial position and the 4 column is the velocity. case 36 is the original case with 36000 cells, sens01 with 60000 cells and sens02 with 105000 cells. when increasing the velocity one can see that the case will not break down. the velocity can be changed in the file, input values. for example setting 1.55 m/s the case will produce reasonable results can converges after about 24000 iterations |
|
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; |
Date Modified | Username | Field | Change |
---|---|---|---|
2011-04-08 15:48 |
|
New Issue | |
2011-05-05 10:11 |
|
Note Added: 0000350 | |
2011-05-05 10:15 |
|
Note Edited: 0000350 | |
2011-05-05 10:18 |
|
Note Edited: 0000350 | |
2011-05-05 11:03 | henry | Note Added: 0000351 | |
2011-05-23 10:46 |
|
Status | new => closed |
2011-05-23 10:46 |
|
Assigned To | => henry |
2011-05-23 10:46 |
|
Resolution | open => no change required |