A NEW multigrid relaxation scheme is developed for the steady-state solution of the Euler and Navier-Stokes equations. The lower-upper Symmetric-Gauss-Seidel method (LUSGS) does not require flux splitting for approximate Newton iteration. The present method, which is vectorizable and unconditionally stable, needs only scalar diagonal inversions. Application to transonic flow shows that the new ...