Changes between Version 9 and Version 10 of RadiationDriving
- Timestamp:
- 02/22/17 13:30:39 (8 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
RadiationDriving
v9 v10 5 5 [http://adsabs.harvard.edu/abs/2017arXiv170206160C Mass Transfer and Disk Formation in AGB Binary Systems] 6 6 7 We apply 2D ray tracing algorithm in radiative transfer calculation. Radiation will transfer momentum (exert radiation force on fluid) withoutenergy transfer. Therefore, we are solving the following equation.7 We apply 2D ray tracing algorithm in radiative transfer calculation. Radiation will transfer momentum (exert radiation force on fluid). However, we do not consider radiation energy transfer. Therefore, we are solving the following equation. 8 8 9 9 $\frac{\partial\rho}{\partial t}+\nabla\cdot\left(\rho\mathbf{u}\right)=0$ … … 12 12 13 13 $\frac{\partial E}{\partial t}+\nabla\cdot\left[\left(E+P\right)\mathbf{u}\right]=\mathbf{f}_{rad}\cdot\mathbf{u}$ 14 15 We average the radiation flux over the azimutal angle ($2\pi$) and divide the polar angle ($\pi$) into n small subdivision. The default is $n=30$ and you can change it in physics/physics_declarations.f90 DustSublimationRadius=30. The radiation will be aborsbed by the gas and dust by opacity. In our paper, opacity can be a function of distance to the luminous object. You can checkout source/radforce.f90 Kappa_try(pos). 16 17 $\frac{d I}{d x}=-\kappa I$