An upwind difference scheme for the double-adiabatic equations (Q676345)

An upstream differencing numerical scheme is devised for the double-adiabatic equations which describe a plasma with anisotropic pressure. The Godunov-type method uses an approximate Riemann solver, which is constructed with an explicit representation of the complete eigensystem of the Jacobian matrix of the flux vector. A correction term is introduced into the numerical flux which allows switching off the instabilities. First results of calculations for the interaction of the solar wind with a comet show that the pressure anisotropies due to the ion pickup modify the plasma flow considerably.