Result

  1. Jonathan helped in debugging. The code can restart on bluehive without producing error.
  1. I ran a 4 AU simulation this week. This simulation has a pulsating AGB star which described before. The difference is: I injected an outflow during 16.83 year and 26.83 year. So the total outflow time window is 10 years. Then the AGB star is set back to normal pulsating state. The result is very exciting. Below are three density plot and contour plot from the simulation. The burst of the outflow has ended and the AGB star is back to pulsating state. The speed of the outflow is 20km/s (remember the pulsating amplitude is only 5km/s), and the escape velocity (with respect to the AGB star only) where the outflow is ejected is about 40km/s. However, we can not treat the gas as ballistic matter. The pressure and temperature will push much of the gas out.

A large amount of the gas is gravitationally bound. I mistakenly overwrite the result of later time. A new simulation is running. I saw the late time evolution, a smooth circumbinary disk will form around the binary. (click the image to enlarge)

A disk forms around the secondary. (click the image to enlarge)

  1. Baowei helped me to set up the code on Bluestreak. I can use 8192 cores on bluestreak. The 3 AU, q=0.1 simulation is schedules to run on bluestreak.

Attachments (3)

Download all attachments as: .zip

Comments

No comments.