We consider the steady-state equations for a Compressible fluid. For low-speed flow, the system is stiff because the ratio of the convective speed to the speed of sound is quite small. To overcome this difficulty, we alter the time evolution of the equations but retain the same steady-state analytic equations. To achieve high numerical resolution, we also alter the artificial viscosity of the numerical scheme, which is implemented conveniently by using other sets of variables in addition to the conservative variables. We investigate the effect of the artificial dissipation within this preconditioned system. We consider both the nonconservative and conservative formulations for artificial viscosity and examine their effect on the accuracy and convergence of the numerical solutions. The numerical results for viscous three-dimensional wing flows and two-dimensional multi-element airfoil flows indicate that efficient multigrid computations of flows with arbitrarily low Mach numbers are now possible with only modifications to existing compressible Navier-Stokes codes. The conservative formulation for artificial viscosity, coupled with the preconditioning, offers a viable computational fluid dynamics (CFD) tool for analyzing problems that contain both incompressible and compressible flow regimes.