A high order Galerkin discretization scheme is used for solving steady compressible Navier-Stokes equations. The pointwise numerical fluxes are separated into convective fluxes, acoustic fluxes, and viscous fluxes. The separation of convective flux and vi
