Changes between Version 30 and Version 31 of u/erica/AccretionModelingBlog
- Timestamp:
- 03/25/18 22:12:51 (7 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
u/erica/AccretionModelingBlog
v30 v31 61 61 Note -- 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..) 62 62 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 .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: 64 64 65 [[Image()]] 65 [[Image(mdot_comp_gamma1pt4.png, 35%)]] 66 67 Despite 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%)]] 66 70 67 71 [[Image()]]