Changes between Version 15 and Version 16 of u/erica/AccretionModelingBlog
- Timestamp:
- 03/20/18 20:04:58 (7 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
u/erica/AccretionModelingBlog
v15 v16 3 3 Am working on a 3D simulation of Bondi Flow onto a central sink particle and studying the accretion properties using the Krumholz accretion algorithm. 4 4 5 The setup has the following parameters (in computational units):5 The following are the fiducial IC's used in the simulations (in computational units), unless otherwise noted: 6 6 7 7 ||Lx=Ly=Lz|| 20 || … … 15 15 $R_{BH} = \frac{GM_*}{C_\infty^2}$ 16 16 17 $\dot{M} = 4 \pi \lambda_{crit}(G M_*)^2c_\infty^{-3}\rho_{\infty}$17 $\dot{M}_{BH} = 4 \pi \lambda_{crit}(G M_*)^2c_\infty^{-3}\rho_{\infty}$ 18 18 19 19 $R_s=\frac{5-3\gamma}{4}\frac{GM_*}{C_{\infty}^2}$ … … 25 25 In the above, $\lambda$ is a numerical coefficient which controls the strength of the accretion flow. At its critical value, $\lambda=\lambda_{crit}$, $\dot{M}$ is strongest. In the code, $\lambda_{crit}$ is solved for using Eqn. 18 in Bondi's paper. 26 26 27 The following are the fiducial IC's used in the simulations (in computational units), unless otherwise noted:27 Additionally, the following sims use (unless otherwise noted): 28 28 29 29 || $C_{\infty}$ || 1.291 || … … 32 32 || $M_*$ || .14318E+13 || 33 33 || $R_{BH}$ || 25.63 || 34 || $\dot{M}$ || 2692.17 || 35 ||Lx=Ly=Lz|| 20 || 36 ||Mx=My=Mz || 64 || 37 || !MaxLevel (AMR) || 2 || 38 || Effective resolution (dx) || .15625 || 34 || $\dot{M}_{BH}$ || 2692.17 || 39 35 40 36 == Research Log: ==