Simulation Update

  • Ionization-only test: Matches exactly the results from our original ionization tests (way back when) up to the final frame at 4ish hours (~0.205 time units). Since this simulation is running in 3D, it doesn't go as far as the potentially troublesome simulation, but I think it's safe to go back to 2D.

http://www.pas.rochester.edu/~adebrech/PlanetIonization/ionization_only0000.png

The long term behavior still appears to not be what we expect — given the behavior of the 1D simulation (see below), it may be due to a higher flux than expected being used.

1D test: With the HotJupiter branch line transfer module, we get a stable ionization equilibrium (did rerun to make sure none of the parameters had been changed). With the radiation pressure implemented, the ionization front moves across the simulation domain to the right, indicating something incorrect happening. However, the momentum is a constant 0, so it doesn't appear to be radiation pressure-related.

http://www.pas.rochester.edu/~adebrech/PlanetIonization/bad_front0000.png

  • Investigated neutral tail coming from the first planet. Velocity vectors show it's originating almost entirely well outside of the planet boundary.

http://www.pas.rochester.edu/~adebrech/PlanetIonization/neutral_tail_e_lim0000.png

  • Have simulations queued for parameter space. For production runs, using total of 150 frames for 15 computational times (5.25 days, or ~30 planetary crossing times).
Rotating Frame
Run # MP (MJ) Flux (phot/cm2/s) Status
1 0.07 2x1013 117 frames complete, running
2 0.263435 2x1013 107 frames complete, queued to restart
3 0.263435 2x1017 Unqueued
4 0.07 2x1017 Unqueued


Non-rotating Frame
Run # MP (MJ) Flux (phot/cm2/s) Status
1 0.07 2x1013 Unqueued
2 0.263435 2x1013 Unqueued
3 0.263435 2x1017 Unqueued
4 0.07 2x1017 Unqueued

Was missing the TABLES from the non-rotating version, not sure why I wasn't getting any error output.

Comments

No comments.