Update 10/2

MHD

The bug was tracked down last week. The low-field test case now works (known working up to a stellar beta of 10-2), so I'm experimenting with some values of the ramping to get the higher-field cases working.

Charge Exchange

Here's a movie up to frame 30:

https://www.pas.rochester.edu/~adebrech/ChargeExchange/HD209458b_charge_exchange0000.png

https://www.pas.rochester.edu/~adebrech/ChargeExchange/HD209458b_charge_exchange_species0000.png

And here's the latest frame:

https://www.pas.rochester.edu/~adebrech/ChargeExchange/HD209458b_charge_exchange0034.png

I've implemented a check to make sure the stellar wind has a greater total pressure than the ambient, and if not to replace the ambient (or material at a density threshold approximately equal to the ambient density) with the stellar wind (which will be a much more significant shock to the simulation, I think, but probably can't be helped). Still testing this.

Scaling Test

Still fighting with the KNL nodes. I'm talking with Baowei to try to figure out why.

Comments

No comments.