COMMON ENVELOPE SIMULATIONS

New Work

  1. Explored the case of binding energy relative to the companion, and plotted several figures.
  2. Snapshots of density at various times.

New Plots of run 143
The plots below are similar to those from the last post but with one main difference: the binding energy contours are for the binding energy relative to the companion. The contours are -0.5 (grey), -0.25 (light grey), 0 (white), 0.25 (thin line, light grey). That is, (i) the kinetic energy density is calculated using the velocity in the frame of the companion, and (ii) the RG core is not included in the calculation of the gravitational potential energy density.

  • The normalized binding energy, with extrema at -1 and 1, is given by:

http://www.pas.rochester.edu/~lchamandy/Graphics/RGB/Post-sink_particle/Post-modified_Lane_Emden/Run143/v_phi_over_v_Kepler_rel_2_v_xyVec_rel_2_E_bind_normalized_rel_2_sliceP2_faceon_lview_60.0Rsun_P1dir_left_0_1_rainbow_0173.png http://www.pas.rochester.edu/~lchamandy/Graphics/RGB/Post-sink_particle/Post-modified_Lane_Emden/Run143/v_phi_over_c_s_rel_2_v_xyVec_rel_2_E_bind_normalized_rel_2_sliceP2_faceon_lview_60.0Rsun_P1dir_left_0_1_rainbow_0173.png

  • Below are the same plots but with an unsaturated color scheme:

http://www.pas.rochester.edu/~lchamandy/Graphics/RGB/Post-sink_particle/Post-modified_Lane_Emden/Run143/v_phi_over_v_Kepler_rel_2_v_xyVec_rel_2_E_bind_normalized_rel_2_sliceP2_faceon_lview_60.0Rsun_P1dir_left_-0.25_1.5_rainbow_0173.png http://www.pas.rochester.edu/~lchamandy/Graphics/RGB/Post-sink_particle/Post-modified_Lane_Emden/Run143/v_phi_over_c_s_rel_2_v_xyVec_rel_2_E_bind_normalized_rel_2_sliceP2_faceon_lview_60.0Rsun_P1dir_left_-0.25_1.5_rainbow_0173.png

  • Below are the same plots but zoomed out, with slightly different scaling for vectors:

http://www.pas.rochester.edu/~lchamandy/Graphics/RGB/Post-sink_particle/Post-modified_Lane_Emden/Run143/v_phi_over_v_Kepler_rel_2_v_xyVec_rel_2_E_bind_normalized_rel_2_sliceP2_faceon_lview_120.0Rsun_P1dir_left_-0.25_1.5_rainbow_0173.png http://www.pas.rochester.edu/~lchamandy/Graphics/RGB/Post-sink_particle/Post-modified_Lane_Emden/Run143/v_phi_over_c_s_rel_2_v_xyVec_rel_2_E_bind_normalized_rel_2_sliceP2_faceon_lview_120.0Rsun_P1dir_left_-0.25_1.5_rainbow_0173.png

  • Below are color plots of the binding energy relative to the companion, since color is easier to understand than contours, at different zoom levels:

http://www.pas.rochester.edu/~lchamandy/Graphics/RGB/Post-sink_particle/Post-modified_Lane_Emden/Run143/E_bind_normalized_rel_2___sliceP2_faceon_cmlview_60.0Rsun_P1dir_left_-1.0_1.0_difference_0173.png http://www.pas.rochester.edu/~lchamandy/Graphics/RGB/Post-sink_particle/Post-modified_Lane_Emden/Run143/E_bind_normalized_rel_2___sliceP2_faceon_cmlview_200.0Rsun_P1dir_left_-1.0_1.0_difference_0173.png http://www.pas.rochester.edu/~lchamandy/Graphics/RGB/Post-sink_particle/Post-modified_Lane_Emden/Run143/E_bind_normalized_rel_2___sliceP2_faceon_cmlview_400.0Rsun_P1dir_left_-1.0_1.0_difference_0173.png

  • Below: Very zoomed out version, with left plot showing binding energy relative to the companion and right plot showing binding energy relative to lab frame (copied from the last post):

http://www.pas.rochester.edu/~lchamandy/Graphics/RGB/Post-sink_particle/Post-modified_Lane_Emden/Run143/E_bind_normalized_rel_2___sliceP2_faceon_cmlview_1200.0Rsun_P1dir_left_-1.0_1.0_difference_0173.png http://www.pas.rochester.edu/~lchamandy/Graphics/RGB/Post-sink_particle/Post-modified_Lane_Emden/Run143/E_bind_normalized___sliceP2_faceon_cmlview_600.0Rsun_P1dir_left_-1.00_1.00_difference_0173.png

  • Below: Same view as above but showing density (left) and temperature (right), for comparison

http://www.pas.rochester.edu/~lchamandy/Graphics/RGB/Post-sink_particle/Post-modified_Lane_Emden/Run143/rho___sliceP2_faceon_lview_1200.0Rsun_P1dir_left_1e-9_1e-2_orangehot_0173.png http://www.pas.rochester.edu/~lchamandy/Graphics/RGB/Post-sink_particle/Post-modified_Lane_Emden/Run143/Temperature___sliceP2_faceon_lview_1200.0Rsun_P1dir_left_100.0_100000000.0_orangehot_0173.png

  • Below: Density at t = 0, 10, 20, 30 and 40 days

http://www.pas.rochester.edu/~lchamandy/Graphics/RGB/Post-sink_particle/Post-modified_Lane_Emden/Run143/rho___sliceP2_faceon_lview_250.0Rsun_P1dir_left_1e-08_0.01_orangehot_0000.png http://www.pas.rochester.edu/~lchamandy/Graphics/RGB/Post-sink_particle/Post-modified_Lane_Emden/Run143/rho___sliceP2_faceon_lview_250.0Rsun_P1dir_left_1e-08_0.01_orangehot_0043.png
http://www.pas.rochester.edu/~lchamandy/Graphics/RGB/Post-sink_particle/Post-modified_Lane_Emden/Run143/rho___sliceP2_faceon_lview_250.0Rsun_P1dir_left_1e-08_0.01_orangehot_0086.png http://www.pas.rochester.edu/~lchamandy/Graphics/RGB/Post-sink_particle/Post-modified_Lane_Emden/Run143/rho___sliceP2_faceon_lview_250.0Rsun_P1dir_left_1e-08_0.01_orangehot_0130.png
http://www.pas.rochester.edu/~lchamandy/Graphics/RGB/Post-sink_particle/Post-modified_Lane_Emden/Run143/rho___sliceP2_faceon_lview_250.0Rsun_P1dir_left_1e-08_0.01_orangehot_0173.png

Next steps

  • Produce the same figures but for run 132 which had Krumholz accretion.
  • Separation vs time plot for both runs on the same graph.
  • Plot(s?) showing orbits of both runs.
  • Polish accretion rate plots for both runs and redo plot for run 132 with higher time resolution.

Comments

No comments.