OutflowWind: 2.5D Tests

Test the OutflowWind module and make sure the results agree with Stone&Proga 2009.

  1. Test with no initial temperature profile

In the version of the outflow object in the code I use, some lines involving energy has E=pressure while others have gamma7*pressure. Modified them all to E=gamma7*pressure though need to double-check with Jonathan…. Added lTempProfile flag in Outflow object. The following is a 2.5D test results for high(10000K) and low(100K) temperature with point mass = 0.1 Jupiter mass. The ambient temperature are 100K and 1K correspondingly. The results shows the planet material escapes for 10000K while are bounded by the point mass for 100K.

T=10000K ​Density Movie;​Temperature Movie
T=100K ​Density Movie;​Temperature Movie
  1. Parameter comparison with Stone&Proga 2009

To compare with the results and parameters used in Stone&Proga paper, we need to double-check the planet_density together with planet_mass, planet_radius and point_mass(?) and temperatures. In the current code it specifies the density, mass & radius separately (in the current default data files, density value is ~35% off) and use the planet mass to calculate the important parameter in Stong&Proga (another important one is the planet density ) around origin. I'm modifying the code to set the planet mass & radius only then calculate the density which is probably a better way and hopefully will be more accurate to compare with the paper. No sure if I should include the point mass when calculating the density…

Comments

No comments.