Accurate determination of surface normal stress in viscous flow from a consistent boundary flux method