How to create faceSet/cellSet with wall faces exceeding a user-defined y+ value
my issue is the following. I have a quite complex 3D geometry and I am using simpleFoam and the kOmegaSST model. I am calculating the y+ value of the wall adjacent cells with my plusPostRANS utility, see http://www.cfd-online.com/Forums/ope...-testcase.html
Suppose I have a grid with: y+_max = 2, y+_min = 0.002 and y+_avg = 0.5. In order to locate the region of all the wall faces (or wall adjacent cell) with y+>1.0 I would like to store either the wall faces or the wall adjacent cells in a faceSet or a cellSet respectively.
The corresponding code snippet of the plusPostRANS utility is
const fvPatch& currPatch = patches[patchi];
if (typeid(currPatch) == typeid(wallFvPatch))
const scalarField& YpTemp = yPlusTemp.boundaryField()[patchi];
So my question is the following: Within the loop, how can I store the wall faces or wall adjacent cells, which have y+>1, into a faceSet or cellSet respectively. With this I would be able to nicely locate them in Paraview using the include sets option. As you can see, I have the y+ values stored in a scalarField YpTemp on which I perform for instance the min and max operations.
As always, any help is appreciated!
Did you consider the option to write the value of y+ to a volScalarField. I think that is much more easy to post-process in paraview.
indeed, thats possible, although I would prefer to create a seperated faceSet or cellSet as it is done when failing a checkMesh check. The reason is, that I would have to play around with contour scaling (max. = user-defined y+ limit) and colors. For me the much more elegant way is to create a faceSet/cellSet, although its maybe slightly more complicated to code. Thatswhy I asked for it and not the volScalarField :)
Anyway for the time being I will create a volScalarField.
Thanks and best regards,
In paraview, you can use the treshold value, and select only the cells with y+ in a certain range of specified values. This has the same effect as specifying a cellSet. This is the route that I would follow probably.
you are right, that should give the same... thanks for the hint! I will give it a try.
|All times are GMT -4. The time now is 13:25.|