Note: the equations below require MathML
to be displayed.
Since central differencing is used for Aero Troll CFD, artificial dissipation is added to control the odd-even uncoupling of the grid and to control nonlinear effects such as shocks (1). The nonlinear scalar model is adapted from reference (2) and from information provided in reference (3).
The numerical implementation of the spectral radius,
, in Aero Troll CFD deviates from what seems to be the standard implementation and instead uses the default implementation as stated in the OVERFLOW manual (3), i.e. the spectral radius in the direction of differentiation.
Note(*) a "2" has been added in front of the spectral radius equation above. The value exists because of historical reasons related to testing Aero Troll CFD. The original idea was to keep the dissipation constants close to OVERFLOW to try to parallel some published OVERFLOW results. However, it turns out that some cases go unstable using the published dissipation constants, therefore the dissipation constants were internally doubled by inserting the "2." The decision was made not to remove the "2" and not redefine the dissipation constants for the public release of Aero Troll CFD so old test cases would not break.
The formulation for
for Aero Troll CFD varies from what is provided in literature. Literature provides the following formulation:
Instead, Aero Troll CFD uses the following formulation for
The reason the different formulation is used is that the formulation provided in literature has problems converging for some supersonic cases due to discrete switching between the
values. On such case for a NACA 0012 airfoil at Mach 1.4 can be seen here
The values for the dissipation constants
shown above come from the OVERFLOW manual (3).
Axisymmetric Artificial Dissipation
are constant in the η direction, therefore:
Artificial Dissipation Jacobian
The max eigenvalue is
Diagonal Artificial Dissipation Jacobian
It can be seen that the last two equations are uncoupled from the first three, therefore the eigenvalues for the last two equations are:
The spectral radius, i.e. the maximum eigenvalue, for the last two equations is γ
The eigenvalues for the first three equations are:
Therefore the spectral radius for the first three equations is 1.
The above diagonalized artificial dissipation Jacobian was compared to results for one were the diagonal had all γs. In general, the two forms of the diagonalized Jacobian behave similar with the one with all γs converging slightly faster on some results. Therefore the diagonalized artificial dissipation Jacobian in Aero Troll contains all γs on the diagonal.
1) Pulliam, T. H., "Artificial Dissipation Models for the Euler Equations," AIAA Journal, Vol. 24, No. 12, Dec. 1986.
2) Slooff, J. W, and Schmidt, W., "Computational Aerodynamics Based on the Euler Equations," AGARD-AG-325, AGARD, Sept. 1984.
3) Nichols, R. H., and Buning P. G., "User's Manual for OVERFLOW 2.1, Version 2.1t Aug 2008," Aug 2008.