[selig5576]

Hi all,

I am writing a solver for the 2D compressible Navier-Stokes equations in conservative form. When it is in conservative form it can be written as a 2D hyperbolic conservation law, but with a viscous flux. To start off simple, I would like to apply the Rusanov Flux (Local Lax Friedrichs) for my flux discretization, how I am unsure how to handle the viscous terms. A part of me wants to separate the fluxes and perform 2nd-order central differencing on the viscous terms and use a standard Rusanov flux for the advection terms, but that would not maintain conservative form.

I have 2 books on solving the compressible Navier-Stokes, and both use MacCormack's method. Furthermore, it uses backward/forward differencing to achieve 2nd-order accuracy. Would I apply a similar approach when using the Local Lax Friedrichs on the viscous flux?