Roe Flux-Difference Splitting Scheme
The inviscid flux vector appearing in Equation 25.5-6 is evaluated by a standard upwind, flux-difference splitting [ 308]. This approach acknowledges that the flux vector contains characteristic information propagating through the domain with speed and direction according to the eigenvalues of the system. By splitting into parts, where each part contains information traveling in a particular direction (i.e., characteristic information), and upwind differencing the split fluxes in a manner consistent with their corresponding eigenvalues, we obtain the following expression for the discrete flux at each face:
Here is the spatial difference . The fluxes and are computed using the (reconstructed) solution vectors and on the "right'' and "left'' side of the face. The matrix is defined by
where is the diagonal matrix of eigenvalues and M is the modal matrix that diagonalizes A, where A is the inviscid flux Jacobian .
For the non-preconditioned system (and an ideal gas) Equation 25.5-10 reduces to Roe's flux-difference splitting [ 308] when Roe-averaged values are used to evaluate . At present, arithmetic averaging of states and is used.
In its current form, Equation 25.5-10 can be viewed as a second-order central difference plus an added matrix dissipation. The added matrix dissipation term is not only responsible for producing an upwinding of the convected variables, and of pressure and flux velocity in supersonic flow, but it also provides the pressure-velocity coupling required for stability and efficient convergence of low-speed and incompressible flows.
An alternative way to compute the flux vector appearing in Equation 25.5-6 is by using a flux-vector splitting scheme [ 57]. The scheme, called Advection Upstream Splitting Method (AUSM), was first introduced by Liou and Steffen in 1993 [ 215]. The AUSM scheme first computes a cell interface Mach number based on the characteristic speeds from the neighboring cells. The interface Mach number is then used to determine the upwind extrapolation for the convection part of the inviscid fluxes. A separate Mach number splitting is used for the pressure terms. Generalized Mach number based convection and pressure splitting functions were proposed by Liou [ 214] and the new scheme was termed AUSM+. The AUSM+ scheme has several desirable properties:
The AUSM+ scheme avoids using an explicit artificial dissipation, by proposing a numerical flux of the form:
Here is the mass flux through the interface, which is computed using the fourth order polynomial functions of the left and right side (of the interface) Mach numbers.
FLUENT utilizes an all-speed version of the AUSM+ scheme based on the low-Mach preconditioning.
Low Diffusion Roe Flux Difference Splitting Scheme
In order to reduce dissipation in LES calculations, FLUENT uses a modified version of the Roe Flux Difference Splitting scheme, called the Low Diffusion Roe Flux Difference Splitting scheme. The scheme includes low Mach number preconditioning, in which the artificial dissipation term has been reduced [ 46] through the use of a hybrid scheme that combines a central scheme and a second-order upwind scheme (Roe's Flux Difference scheme).
| The low diffusion discretization must be used only for subsonic flows. For high Mach number flows, you should switch to the second-order upwind scheme.
| The low diffusion discretization is only available with the implicit-time formulation (dual-time-stepping). When running LES with the explicit time formulation, you will need to use the second-order upwind scheme.