A new defect‐correction method for the stationary Navier–Stokes equations based on local Gauss integration