We outline a new class of robust and efficient methods for solving the Navier– Stokes equations. We describe a general solution strategy that has two basic building blocks: an implicit time integrator using a stabilized trapezoid rule with an explicit Adams–Bashforth method for error control, and a robust Krylov subspace solver for the spatially discretized system. We present numerical experiments illustrating the potential of our approach. Key words. time-stepping, adaptivity, Navier–Stokes, preconditioning, fast solvers AMS subject classifications. 65M12, 65M15, 65M20 DOI. 10.1137/080728032
David A. Kay, Philip M. Gresho, David F. Griffiths