I think you have missed the point all together…
I thought I presented the issue very clearly in a step by step method.
So, in an even more concise presentation:
- The integration of the X component of the Wall Shear Stress vector over the geometry surface should result in the value of viscous drag force (if that is incorrect, please show me why it is incorrect…)
- Visually the surface mapping of the X component of the Wall Shear Stress vector, seems to indicate that the result of that integration will be a -ve number (do you agree with that?)
- Regardless of what we visually can say, what is the sign of the actual mathematical integration? My ParaView mathematical integration of that X component results in a -ve value while OpenFoam shows a +ve number. (Please show me that I have incorrectly mathematically integrated the X component in ParaView…)
- Right now, the only issue is… Why does OpenFoam show +ve Viscous X Force and ParaView mathematical integration results in a -ve value, of the same magnitude as OpenFoam?
Forget everything about what is logical here, forget whether Bryces y+=60 Wall Function mesh is relevant, forget about even trying to see which result matches experimental values better and forget about divergence schemes…
Concentrate on finding out why ParaView integrates to a -ve value and OpenFoam integrates to a +ve number…
I hope I made a mistake 