I now checked Aspect ratio, Non-orthogonality, Volume ratio, Edge ratio, Cell volume and non of these show issues in the Mesh Quality at the locations mentioned in the Runs:
Run 3: at Gauge pressure = -1.04349e+15 at position: (-54.24 m, -76.35 m, 0.5804 m).
and Run 4 at: Gauge pressure = 1.41853e+15 at position: (-86.87 m, -46.04 m, 0.5278 m).
I was not able to find ‘Geometry primitives’ option to locate these points, so at the end I looked at an earlier successful run with the same model and with x, y, z (in Solution Fields) I found the locations. As said no issues were seen there in the Mesh quality view.
At the end I reduced the Roughness height of Custom (from 10m to 1m) and now I can Run it successfully. If I try Roughness heights 3, 5 or10m they all cause errors in run (I did not try 2 yet).
So I did not yet find a useful workflow to get it working.