Implicit residual smoothing operators for accelerating the convergence of explicit, multistage, artificial compressibility algorithms are developed using ideas from pressure-based methods. The velocity derivatives in the continuity equation and the pressure gradient terms in the momentum equations are discretized in time implicitly. The discrete system of equations is linearized in time producing a block implicit operator which is approximately factorized and diagonalized via a similarity transformation. The so-derived diagonal operator depends only on the metrics of the geometric transformation and can, thus, be implemented in an efficient and straightforward manner. It is combined with the standard implicit residual smoothing operator and incorporated in a four-stage Runge–Kutta algorithm also enhanced with local time-stepping and multigrid acceleration. Linear stability analysis for the three-dimensional Navier–Stokes equations and calculations for laminar flows through curved square ducts and pipes demonstrate the damping properties and efficiency of the proposed approach particularly on large-aspect ratio, highly skewed meshes.