[FLASH-USERS] Fwd: Re: 2.5D split hydro & centrifugal force

[Chris Lindner]

Hi Chris,

My group has implemented angular momentum in 2D and 1D simulations by - as
Sean mentioned - adding a mass scalar which represents angular momentum.  If
you just use velz in cylindrical or spherical coordinates, you end up
creating some ambiguity about how total and kinetic energy should be
calculated.  We incorporate the centrifugal force term in the
gravity_accelOneRow .  We also include an alpha-viscosity treatment of
angular momentum transport.  We edited the diffuse_species module and a few
other spots to do this.  We describe our implementation and some tests of
this we did in some of our recent papers:
http://arxiv.org/PS_cache/arxiv/pdf/1108/1108.1415v1.pdf
http://arxiv.org/PS_cache/arxiv/pdf/0910/0910.4989v1.pdf .  Let me know if
you run into any trouble or have some more questions about it.

Chris Lindner
lindner at astro.as.utexas.edu
NSF Graduate Fellow
Department of Astronomy
The University of Texas at Austin
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20110811/8919a02e/attachment.htm>