Update 10/21

MHD

I think a CFL issue has been successfully ruled out (Jonathan, if you think a direct wave speed comparison would be useful, I can do that - but based on the tests already done, don't think they're necessary).

Is the prolongation of the aux fields still the most likely place for a bug? It may be time to sit down for a bug hunt.

Charge exchange

Still moving right along. Hoping we reach a steady state soon. Here's the latest frame. (At some point I'll get a movie, but waiting to hear back about the paper so we can hopefully archive the rad pressure run and make room….)

http://www.pas.rochester.edu/~adebrech/ChargeExchange/charge_exchange_species0017.png

http://www.pas.rochester.edu/~adebrech/ChargeExchange/charge_exchange0017.png

I meant to turn off the lower y boundary, turned off the upper one instead - been corrected, and shouldn't affect the simulation as it progresses (but I've queued the restart from an earlier frame anyway).

AMR line transfer

Set up a stack, still need to change ray tracing to use it. Think it should be a stack of rays rather than ray positions, which means the main ray tracing loop needs to be rejiggered again.

Once the MHD bug is worked out I should have time to go a little deeper into this, see if there's a better way to approximate the graph that represents the dependencies.

Comments

No comments.