SUMMARYWe consider higher-order mixed ÿnite elements with continuous pressures for the computation of stationary compressible ows at low Mach number. The proposed approach is based on a fully coupled treatment of the governing equations and therefore, for steady-state calculations, does not rely on timestepping techniques. The non-linear problem is solved by means of a quasi-Newton iteration. The strongly coupled system resulting from higher-order discretization of the linearized equations requires adequate solvers. We propose a new scheme based on multigrid methods with varying FEM ansatz orders on the grid hierarchy as well as multiplicative smoothers based on blocking techniques. Computational results are described for a benchmark conÿguration including a ow with heat transfer in the low Mach number regime. Furthermore, the issue of anisotropic grids is addressed in that context.