A fully upwind implicit PNS solver for the computation of 3D supersonic flows with crossflow separation