Changes between Version 9 and Version 10 of RadiationDriving


Ignore:
Timestamp:
02/22/17 13:30:39 (8 years ago)
Author:
Zhuo Chen
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • RadiationDriving

    v9 v10  
    55[http://adsabs.harvard.edu/abs/2017arXiv170206160C Mass Transfer and Disk Formation in AGB Binary Systems]
    66
    7 We apply 2D ray tracing algorithm in radiative transfer calculation. Radiation will transfer momentum (exert radiation force on fluid) without energy transfer. Therefore, we are solving the following equation.
     7We 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.
    88
    99$\frac{\partial\rho}{\partial t}+\nabla\cdot\left(\rho\mathbf{u}\right)=0$
     
    1212
    1313$\frac{\partial E}{\partial t}+\nabla\cdot\left[\left(E+P\right)\mathbf{u}\right]=\mathbf{f}_{rad}\cdot\mathbf{u}$
     14
     15We 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$