The goal of this paper is to introduce a new multilevel solver for two-dimensional elliptic systems of nonlinear partial differential equations (PDEs), where the nonlinearity is of the type u∂v. The incompressible Navier-Stokes equations are an important representative of this class and are the target of this study. Using a first-order system least-squares (FOSLS) approach and introducing a new...