I will try it tomorrow, now I don’t have access to the files. However for the moment I can say that that is a behaviour that has replicated itself across more than 20 runs, disregarding the mesh, the freestream values and even the solver used (tried with both pressure based and density based).

—-

Moreover today I have tried for the first time to reproduce the same kind of plot with an open source software, exporting the solution data of a case affected by the problem in the Ensight gold format, and they are displayed correctly and qualitatively in line with the expected physics. I honestly don’t know the details of the implementations, but I have tried two approaches:

1. Construct the Wall Shear vector from the exported field values of x_wall_shear,…,z_wall_shear, and plotted the “streamlines” following that vector field, constricted at the surface.

2. Plotted a Line Integral Convolution of the Wall Shear vector field (computed in the same way as point 1).

Both the approaches works. This should rule out a problem in the solver, do you agree?