wiki:u/erica/AccretionModelingBlog

Version 8 (modified by Erica Kaminski, 7 years ago) ( diff )

Bondi Accretion

3/18/2018

Am working on a 3D simulation of Bondi Flow onto a central sink particle and studying the accretion properties using the Krumholz accretion algorithm.

The setup has the following parameters (in computational units):

Lx=Ly=Lz 20
MaxLevel (AMR) 2
Mx=My=Mz 64
Effective resolution (dx) .15625
Sink dx = 4dx .625

From Bondi (1952) we have the following quantities (in order: Bondi radius, accretion rate, sonic radius, nondimensional density and velocity):

In the above, is a numerical coefficient which controls the strength of the accretion flow. At its critical value, , is strongest. In the code, is solved for using Eqn. 18 in Bondi's paper. In this simulation, these parameters are:

.2526
1.66

Additionally, we have:

25

The case of flow is an extreme limit for the solution space. This value of has a few interesting features:

Results

The simulation setup reproduces the correct nondimensional profiles:

No image "Bondi_comparison.png" attached to u/erica/AccretionModelingBlog

Library

Attachments (29)

Note: See TracWiki for help on using the wiki.