[Spidermantuba]

I presently have written a code that evolves Euler's equations with Newtonian self-gravity on an Eulerian Cartesian grid. I use the Kurganov-Tadmor central difference scheme for the advection.

I would like to be able to numerically conserve both linear and angular momenta. Does anyone know of any way to do this? I have figured out the fundamental reason the K-T method does not conserve angular momentum when advecting linear momenta - it does not preserve the symmetry of the stress tensor. But I can't figure out how to alter the scheme to account for this effect.


Any suggestions?