AstroBEAR/ATHENA debugging

Initial Conditions

Initial conditions for the planet are finally correct. The ambient media still differ - the planet ends slightly earlier, and the profile is very steep near the edge.

Density

http://www.pas.rochester.edu/~adebrech/PlanetIonization/comparisons/test4/ramp_densities0000.png

Temperature

http://www.pas.rochester.edu/~adebrech/PlanetIonization/comparisons/test4/ramp_temp0000.png

Pressure

http://www.pas.rochester.edu/~adebrech/PlanetIonization/comparisons/test4/ramp_pressure0000.png

Ionization Fraction

http://www.pas.rochester.edu/~adebrech/PlanetIonization/comparisons/test4/ramp_ionization0000.png

Radiative Transfer

  • The optical depths are close, but not quite identical. At ~300 s,

http://www.pas.rochester.edu/~adebrech/PlanetIonization/comparisons/test4/ramp_tau0001.png

Working backwards, I find \sigmaH = ~6.28x10-18 for AstroBEAR, and ~6.14x10-14 for ATHENA - tested, and the difference isn't significant for the optical depth:

http://www.pas.rochester.edu/~adebrech/PlanetIonization/comparisons/test4/tau_new_sigma.png

http://www.pas.rochester.edu/~adebrech/PlanetIonization/comparisons/test4/sigma_diff.png (ATHENA in blue, AstroBEAR in orange)

  • Ionization/recombination rates and heating/cooling rates are qualitatively similar, but still differ by 1-2 orders of magnitude.

Ionization

http://www.pas.rochester.edu/~adebrech/PlanetIonization/comparisons/test4/ramp_ionization0001.png

(Without advective terms)

http://www.pas.rochester.edu/~adebrech/PlanetIonization/comparisons/test4/ion_recom_no_adv.png

Theory Comparisons

It appears that the actual data from ATHENA doesn't match the theoretical calculation (using the data for n, X, and \tau in the equations on the LineTransfer page).

http://www.pas.rochester.edu/~adebrech/PlanetIonization/comparisons/test4/ionization.png

Heating

http://www.pas.rochester.edu/~adebrech/PlanetIonization/comparisons/test4/ramp_heating-cooling0001.png

(Without advective terms)

http://www.pas.rochester.edu/~adebrech/PlanetIonization/comparisons/test4/heating_cooling_no_adv.png

  • Electron densities are correct:

http://www.pas.rochester.edu/~adebrech/PlanetIonization/comparisons/test4/ramp_densities0001.png

Dynamics

For the first time, the dynamics are beginning to look sort of similar. Radial velocities are too large outside of the planet for AstroBEAR:

http://www.pas.rochester.edu/~adebrech/PlanetIonization/comparisons/test4/ramp_vx0001.png

And the evolution of the temperature on the day side is correct:

http://www.pas.rochester.edu/~adebrech/PlanetIonization/comparisons/test4/ramp_temp0001.png

Comments

No comments.