Posts for the month of December 2018

Update 12/17

Radiation Pressure Paper

Figures pretty much done. Need to put the medium flux and high flux average results side by side, and most could use some minor tweaking for presentation. Currently only using the streamline image for the high flux case, but may want to add at least radiation pressure (especially if it turns out that a lot of the absorption is happening in a fast-moving shell).

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

High flux results are as expected on average, but the results for the blowback phase are worrisome - could use a second opinion.

Bubbleout (245-290)

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

Full average (245-315)

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

Blowback (290-315)

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

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

Will have a first draft shared with everyone by end of day Friday.

Ideas for improving simulation

  • Adding a tracer to see the exact layers structure of the envelope

Updates 12/03

PNe in big box

Use the updated problem module, 4 level AMR, and a cubic box with 128 kR_sun radius (64000 in computational unit, base grid 1000 R_sun). I did two runs over the past week. The first one has minTemp=10K, and mindensity=1d-10 for all. In the second one, I raised those two value to doubled for frame 0, and then lower them back for the rest of the simulation.

  • Links to movies and graphs:

Nov29 Dec2

  • "Disk" — lower density along equatorial plane
  • "cold bubble" in temperature plots

http://www.pas.rochester.edu/~yzou5/PNe_simulation_movies/Dec2_bigbox/rho0000.png http://www.pas.rochester.edu/~yzou5/PNe_simulation_movies/Dec2_bigbox/rho0001.png http://www.pas.rochester.edu/~yzou5/PNe_simulation_movies/Dec2_bigbox/rho0054.png

http://www.pas.rochester.edu/~yzou5/PNe_simulation_movies/Dec2_bigbox/Tmep0000.png http://www.pas.rochester.edu/~yzou5/PNe_simulation_movies/Dec2_bigbox/Tmep0001.png http://www.pas.rochester.edu/~yzou5/PNe_simulation_movies/Dec2_bigbox/Tmep0054.png

http://www.pas.rochester.edu/~yzou5/PNe_simulation_movies/Dec2_bigbox/M0000.png http://www.pas.rochester.edu/~yzou5/PNe_simulation_movies/Dec2_bigbox/M0001.png http://www.pas.rochester.edu/~yzou5/PNe_simulation_movies/Dec2_bigbox/M0054.png

COMMON ENVELOPE SIMULATIONS

Summary of recent work

  1. Continued to write/edit energy paper (see overleaf)
  2. Plots for energy paper (see en_fig.pdf, last post)
  3. New longer movie of AGB run (see last post)
  4. Energy plots for AGB run (this post)
  5. Plan for other papers

Energy paper

Overleaf: https://www.overleaf.com/project/5bd20d99379a851a13bfd103

2D plots for AGB run 164: http://www.pas.rochester.edu/~lchamandy/CE_papers/Energy/en_fig_run164.pdf

Planned papers

  • Parameter space—dependence on secondary mass (~15 pages):
    • Compare energy/mass unbinding for different runs.
    • Analysis of orbital dynamics including damping force. Compare with theory and compare the three different runs.
    • Convergence tests.
    • Suggest to extend the runs in time somewhat (but stop before material exits box).
    • ALTERNATIVELY: could study orbital dynamics for the fiducial run in a shorter separate paper.
  • Jet (Letter): See previous posts
  • AGB (Letter)
    • Present the run since not really done by anyone before
    • Compare with RGB run
      • Orbit
      • Energy
      • Orbital dynamics
    • QUESTION: Do we perform a new run or use the present run? Material has just started to exit the box!

Update 12/3

Radiation Pressure

Another 3 frames done. I have it running still, but if this is sufficient, I can stop it and start postprocessing so we can get the remaining plots for the paper.

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

Parameter space

Spent a good portion of last week retrieving data and making movies. Checking the proofs and sending back today - should be completely done by tomorrow evening.