Update 11/20

  • From discussion on Wednesday, see this page.
  • For planet density in the high-flux case (assuming recombination-limited flows), see previous blog post.

Runs

Currently running small planet, low flux case (still recombination-limited, but looks like thanks to tidal forces the surface is staying well above the reset radius). Modified parameters (from first link above):

Mp 0.07 MJ
TP 1000 K

Also queued high flux, same size planet case. Python script​ was modified so that the planet profile is calculated assuming the totally-neutral surface is at the planet radius (though this is absolutely not the case for the new density). Modified parameters:

Ionizing flux 2x1017 phot/cm2/s
5.3x10-9 g/cm3

Code stuff

Got another soft crash on the small planet - need to investigate (race condition seems likely - Jonathan suggested layout_comms).

Running ionization-only test with radiation pressure implemented again to confirm and investigate results (second link above).

Further thoughts

Do we want to think about including a stellar wind at all? My thought: maybe for later simulations, but not for the parameter study, and probably not for the radiation pressure.

We had discussed possibly increasing the distance in y to capture more of the up- and downwind arms — may be nice for production runs at least (I wouldn't expect it to break anything, but could be worth testing first).

Attachments (1)

Download all attachments as: .zip

Comments

No comments.