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).
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)
Full average (245-315)
Blowback (290-315)
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:
- "Disk" — lower density along equatorial plane
- "cold bubble" in temperature plots
COMMON ENVELOPE SIMULATIONS
Summary of recent work
- Continued to write/edit energy paper (see overleaf)
- Plots for energy paper (see en_fig.pdf, last post)
- New longer movie of AGB run (see last post)
- Energy plots for AGB run (this post)
- 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.
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.