Towards an efficient and robust high order accurate flow solver for viscous compressible flows