Changes between Version 30 and Version 31 of u/erica/AccretionModelingBlog


Ignore:
Timestamp:
03/25/18 22:12:51 (7 years ago)
Author:
Erica Kaminski
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • u/erica/AccretionModelingBlog

    v30 v31  
    6161Note -- I did this lineout from the center of the first zone in the octant surrounding the sink, to the zone center at the edge of the boundary. Aside from doing a spherical average over shells, I think this will provide the most accurate results from Visit (rather than along cell boundaries, which I am not sure how Visit interpolates..)
    6262
    63 Calculating the mass flux through spherical shells surrounding the sink seems to agree well with the predicted value of $\dot{M}_{BH}=8646.9$. I tried at 3 different radii: $r=5<R_s$, $r=7>R_s$, $r=25>OBS$, where $OBS$ is the spherical boundary in the grid beyond which the Bondi solution is pasted into cells.
     63Calculating the mass flux through spherical shells surrounding the sink seems to agree well with the predicted value of $\dot{M}_{BH}=8646.9$. I tried at 3 different radii: $r=5<R_s$, $r=7>R_s$, $r=25>OBS$, where $OBS$ is the spherical boundary in the grid beyond which the Bondi solution is pasted into cells:
    6464
    65 [[Image()]]
     65[[Image(mdot_comp_gamma1pt4.png, 35%)]]
     66
     67Despite these agreements with theory, and the fact that the sonic point lies well beyond the accretion kernel, there is a spurious outflow being generated within the accretion kernel. Contrary to intuition, this outflow propagates beyond the accretion kernel into the grid, disrupting the upstream hydrodynamic solution.. This happens ''very'' quick in the simulation ($\sim1/30 ~t_{sc, amb}$). Here is a snapshot of density illustrating this effect:
     68
     69[[Image(density_blowup_gamma1pt4.png, 35%)]]
    6670
    6771[[Image()]]