[Rafal_Hemzaczek]

Hello all,
I am running a 2D compressible k-omega SST simulation of the internal flow inside a C-D nozzle using shockFluid OpenFOAM V13. I have a good understanding of the flow as it can be estimated analytically, there is real test data and I had run an inviscid simulation of it previously.


I used this information to calculate the cell size required for y+ = 1.
Having then ran the simulation I used the command:
`foamPostProcess -solver shockFluid -func yPlus`


The results are suprising as they are saying I was off by 2-3 orders of magnitude and that my mesh is very over resolved (y+ max: 0.0032).
However I have my suspicions:


a) The cmd line logs first say: `transport sutherland;`in the thermodynamics package description, but then follows with:
`Selecting default RAS thermophysical transport model unityLewisEddyDiffusivity`
I think this is what is causing the results to be incorrect as the transport model should be Sutherland. [see screenshot 1]


b) I did a manual check of y+ at the throat using the values of wall shear stress that were generated using:
`foamPostProcess -solver shockFluid -func wallShearStress`
This results in the same contradicting cmd message (different thermophysical transport models) however the values it generates are belivable (14 vs my original estimate of 19).
Using this and all other values from the simulation output the manually calculated y+ for the first cell is 0.58 [screenshot 2] which is not too far off my estimates (meshed for target y+ at throat of 0.8).


I am teaching myself CFD from textbooks and quite new so I am not sure if I have understood something incorrectly or the issue is with the way I have called up the yPlus post processing utility?


I would just like to validate my mesh is correctly sized for resolving the flow near the walls everywhere (not just the nozzle throat), before I carry on experimenting


Thank you in advance!
Raf