Skew Symmetric Schemes

Skew symmetric finite difference schemes for compressible and incompressible flows are subject of academic research since they where proposed in the early eighties (Feiereisen, et al., 1981;
Tadmor, 1984). They allow to strictly avoid numerical damping. They however find little application in simulations of engineering practice. This might be partly due to the different, uncommon concepts
compared to finite volume (FV) schemes, due to the sometimes complex implementation, but also due to the lack of easy to use concepts for handling central issues like transformed grids, time integration,
boundary conditions and block-block interfaces.


The aim of our research is to develop the skew symmetric scheme further to a valuable scheme for both research and engineering applications.

The skew symmetric discretization builds on rewriting the Navier-Stokes or, shown here, Euler equations in the skew symmetric form, in one dimension given by

Here ρ, p, γ are density, pressure, the adiabatic index and ui (i =1,2,3) are the velocity components. The dots in the momentum equation (2) mark, that the derivatives act also on the u right of the parenthesis.
The unusual time operator in the momentum equation is a consequence of the skew symmetric rewriting of the equation.

When using skew symmetric derivative matrices it can be shown that the nonlinear transport conserves the kinetic energy in the discrete case, as in the analytical theory. This is in contrast to the majority
of numerical schemes with the exception of spectral methodes. Further the schemes used by us also conserve mass, total momentum and total energy, even though equations (2) and (3) are clearly not in
divergence form, i.e. are not equivalent to FV. The time has to be discretized carefully to keep the exact conservation. One way to do it is due to Morinishi (Morinishi, 1987) by rewriting it in (2) as

This can be discretized in time by an adopted, implicit mid-point rule, which gives a time integration scheme similar to one obtained by (Subbareddy, 2009), being fully conservative for finite time steps.
Other time integration schemes developed by us, fully and partly conservative can be derived.

Grafik2 Kontakt

Julius Reiss

Dr. Rer. Nat. Julius Reiss

FG Numerische Fluiddynamik