Changes between Version 25 and Version 26 of u/erica/AccretionModelingBlog


Ignore:
Timestamp:
03/20/18 20:26:03 (7 years ago)
Author:
Erica Kaminski
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • u/erica/AccretionModelingBlog

    v25 v26  
    5353|| $R_{s}$ || .0128 ||
    5454
    55 Note, the sonic radius is not resolved by ~ factor of 20 in this simulation, as seen when checking dx quoted above. In this sim, the Bondi solution was initialized everywhere in the grid at t=0 (except within a small inner radius, described below). In subsequent timesteps, the Bondi solution is re-pasted into the ghost zones and a small inner spherical region at the origin. A sink particle is initialized at the origin with mass given above ($M_{sink}=M_*$). The particle can accrete gas following the Krumholz prescription. The only source of gravity in the sims is the point gravity object associated with the particle (self-gravity is turned off).
     55Note, the sonic radius is not resolved by ~ factor of 20 in this simulation, as seen when checking dx quoted above.
     56
     57In this sim, the Bondi solution was initialized everywhere in the grid at t=0 (except within a small inner radius, described below). In subsequent timesteps, the Bondi solution is re-pasted into the ghost zones and a small inner spherical region at the origin. A sink particle is initialized at the origin with mass given above ($M_{sink}=M_*$). The particle can accrete gas following the Krumholz prescription. The only source of gravity in the sims is the point gravity object associated with the particle (self-gravity is turned off).
    5658
    5759The simulation setup reproduces the correct nondimensional profiles. Here's a comparison of the t=0 nondimensional density and velocity profiles astrobear calculates for the bondi module (left), compared to the Bondi solutions for $\gamma=5/3$ (curves II and III for the case of $\lambda=\lambda_{crit}$).