Update 8/3

  • ATHENA and AstroBEAR comparisons are working nicely. There are still a lot of differences between the simulations, however. Dashed lines are ATHENA, solid are AstroBEAR.

​Initial profile:

Initial conditions

Flattening radius - ATHENA flattens at r_ib rather than r_mask
Electron densities - ATHENA goes to near zero inside planet, AstroBEAR is just 10-10 of the H density (this is probably not significant in the end)
r_ob - ATHENA profile extends to r_zero, but uses values at r_ob for density and pressure at edge
Pressure profiles differ

These last two seem likely to cause many of the differences seen once things actually start happening.

Frame 1

Looks like this ATHENA run was actually run with ramping - John seemed to believe it was not, so either the flag he's using isn't successful or he grabbed the wrong simulation accidentally
dE/dt for PdV, Lyman-alpha, and ionization heating all differ fairly significantly
AstroBEAR recombination is incredibly small - likely not attributable to differences in profiles

Other differences at this point are almost certainly attributable to the difference in flux or the difference in outer boundaries.

​Adjusted outward to r_zero:

With the profile extended to rzero (and the resolution increased by one level), the location of the outer boundary becomes closer, but the values no longer match. So something is going wrong either in the calculation or import of the profile, still.

  • Gave Diyora edge-centered VTK output from VisIt volume plot to start with. I almost have the layout working, but can't seem to get the refined data (only the coarsest is showing in the chombo). Is only level 0 data loaded BeforeGlobalStep?

http://www.pas.rochester.edu/~adebrech/HD209458b/density_layout.png

  • Did testing for OpenMP version of code. Problem appears sometime between 5 tasks/4 threads and 10 tasks/12 threads (haven't determined exact cutoff yet).
  • Working on presentation for eclipse at Ogden Farmers' Library.

Attachments (1)

Comments

No comments.