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.
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:
Provides exact resolution of contact and shock discontinuities
Preserves positivity of scalar quantities
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.
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.