wiki:u/erica/AccretionModelingBlog

Version 6 (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 addition, the Bondi solution requires a numerical constant, , which controls the strength of the accretion flow. At its critical value, , is strongest. The case of flow is an extreme limit for the solution space. This value of has a few interesting features:

In the code, is solved for numerically for given value of . In this simulation, these parameters are:

.2526
1.66

Additionally, we have:

.2526
1.66

Library

Attachments (29)

Note: See TracWiki for help on using the wiki.