Context Navigation

← Previous Post
Next Post →

Post common envelope simulation

I put a 1 solar mass object at the origin of the simulation box. There is no outflow from 0 to 500 day. From 500 to 1000 day, there is outflow for the boundary of the object. I list the physical quantities below:

Radius of the object that has outflow: r=46.875R_{\odot}

Mass of the object: 1M_{\odot}

Density of the outflow: \rho=10^{-11}g\cdot cm^{-3}

Temperature of the object: T_{1}=30000K

The outflow is asymmetric, in the polar direction (theta=0 and theta=pi), the outflow velocity is the maximum; in the equatorial plane, the outflow velocity is 0. The outflow velocity is a function of polar angle, there are two kinds of functions, linear and Gaussian - to vary the open angle. The Gaussian is narrower.

Linear: v_{r}(\theta)=v_{max}\times\frac{\|\pi-2\theta\|}{\pi}

Gaussian: v_{r}(\theta)=v_{max}\times\exp(-10(\frac{\|\pi-\|2\theta-\pi\|\|}{\pi})^{2})

Spherical: v_{r}(\theta)=v_{max}

Below shows the linear and Gaussian velocity v.s. angle.

The outflow temperature: 30000K

Outflow duration: 500 day - 1000 day

Mass loss: \dot{M}=\int_0^{\pi}2\pi\rho r^{2}\sin\theta v_{r}(\theta)d\theta. However, much of the outflow will be impeded by the surrounding gas, the actual mass loss will be lower than the calculated ones.

ID	v_{r}(0)(km\cdot s^{-1})	\rho (g\cdot cm^{-3})	\dot{M}(M_{\odot}\cdot yr^{-1})	simulation	approximate expansion velocity
1	300 (linear)	10^{-11}	2.31\times10^{-4}		120(km\cdot s^{-1})
				movies:density and temperature
2	400 (linear)	10^{-11}	3.07\times10^{-4}

3	400 (gaussian)	10^{-12}	1.00\times10^{-5}

4	400 (gaussian)	10^{-11}	1.00\times10^{-4}		310(km\cdot s^{-1})
				movies:density and temperature
5	300 (spherical)	10^{-11}	6.35\times10^{-4}		240(km\cdot s^{-1})
				movies:density and temperature

Simulations that has indefinite outflow duration. (it is always on)

ID	v_{r}(0)(km\cdot s^{-1})	\rho (g\cdot cm^{-3})	\dot{M}(M_{\odot}\cdot yr^{-1})	simulation	approximate expansion velocity
6	300 (spherical)	10^{-11}	6.35\times10^{-4}		movies:density

2	1000 (spherical)	10^{-16}	2.12\times10^{-8}

Some questions: where exactly does the jet emerge? What are the open angle and outflow velocity. What are the outflow density and temperature?

Posted: 7 years ago (Updated: 7 years ago)
Author: Zhuo Chen
Categories: (none)

Attachments (5)

gaussian.png (7.6 KB) - added by Zhuo Chen 7 years ago.
linear.png (8.6 KB) - added by Zhuo Chen 7 years ago.
1_137.png (402.9 KB) - added by Zhuo Chen 7 years ago.
4_105.png (462.5 KB) - added by Zhuo Chen 7 years ago.
5_177.png (357.2 KB) - added by Zhuo Chen 7 years ago.

astrobear

Context Navigation

Post common envelope simulation

Attachments (5)

Comments

Download in other formats: