0000179: k-omega sst model breaks down with high resolution grid - OpenFOAM Issue Tracking

ID	Project	Category	View Status	Date Submitted	Last Update

0000179	OpenFOAM	Bug	public	2011-04-08 15:48	2014-02-10 13:16

Reporter	~~user171~~	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.

~~user171~~ 2011-05-05 10:11 ~0000350 Last edited: 2011-05-05 10:18	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

henry 2011-05-05 11:03 manager ~0000351	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.0nuw[faceI]/(beta1_sqr(y[faceI])); scalar omegaLog = sqrt(k[faceCellI])/(Cmu25kappa_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	~~user171~~	New Issue
2011-05-05 10:11	~~user171~~	Note Added: 0000350
2011-05-05 10:15	~~user171~~	Note Edited: 0000350
2011-05-05 10:18	~~user171~~	Note Edited: 0000350
2011-05-05 11:03	henry	Note Added: 0000351
2011-05-23 10:46	~~user2~~	Status	new => closed
2011-05-23 10:46	~~user2~~	Assigned To	=> henry
2011-05-23 10:46	~~user2~~	Resolution	open => no change required