23.5.3. Convective Fluxes

23.5.3.1. Roe Flux-Difference Splitting Scheme

The inviscid flux vector appearing in Equation 23–81 is evaluated by a standard upwind, flux-difference splitting  [559]. 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 (that is, 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:

(23–85)

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

(23–86)

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 23–85 reduces to Roe’s flux-difference splitting  [559] when Roe-averaged values are used to evaluate . At present, arithmetic averaging of states and is used.

In its current form, Equation 23–85 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.

23.5.3.2. AUSM+ Scheme

An alternative way to compute the flux vector appearing in Equation 23–81 is by using a flux-vector splitting scheme  [107]. The scheme, called Advection Upstream Splitting Method (AUSM), was first introduced by Liou and Steffen in 1993 [381]. 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 [380] and the new scheme was termed AUSM+. The AUSM+ scheme has several desirable properties:

  1. Provides exact resolution of contact and shock discontinuities

  2. Preserves positivity of scalar quantities

  3. Free of oscillations at stationary and moving shocks

The AUSM+ scheme avoids using an explicit artificial dissipation, by proposing a numerical flux of the form:

(23–87)

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.

Ansys Fluent utilizes an all-speed version of the AUSM+ scheme based on the low-Mach preconditioning.

23.5.3.3. Low Diffusion Roe Flux Difference Splitting Scheme

In order to reduce dissipation in LES calculations, Ansys 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  [87] 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 Roe-FDS is available in special circumstances when the LES viscous model is selected and when using the implicit-time formulation (dual-time-stepping). It is recommended that you use the second-order upwind scheme for maximum accuracy when selecting this flux scheme.


Important:  The low diffusion discretization must be used only for subsonic flows. For high Mach number flows, it is recommended that you switch to standard Roe-FDS or AUSM+ flux schemes and second-order upwind discretization.



Important:  When running LES with the explicit-time formulation, Low Diffusion Roe-FDS is not available and you will need to either use the standard Roe-FDS or the AUSM+ flux scheme and second-order upwind discretization.


To learn how to apply the different convective fluxes, see Convective Flux Types in the User's Guide.