Hi Martin,
Thank you for your help.
I plot the X-axis component of “Shear Stress” from the postprocessing output of Code_Satrune on the Boundary surface, which is the solver’s standard output file. The result plotted is circumferentially averaged. Yes, the x direction is the flow direction.
In answer to your second question, I use the k-omega SST turbulence model without wall function.
I have a little update since yesterday night. I gave a blending factor 0.8 to the spatial scheme for momentum equaiton (say velocity in the GUI) with “Automatic” scheme (I guess it is “Centred”). The Drag force collected with the subroutine in my initial post gives an error with the measurement for about 10%, a bit better. This also mitigates the shear stress oscillation from my previous simulation (which is almost the same issue as seen in shear stress is checkboarding - #5 by romoli).
I read another recent post (Skin friction coefficient calculation) saying the drag force prediction is around 4% with the version 8.1. So I am compiling this version now, and will give an update later.
Best regards,
Feng