Changes between Version 28 and Version 29 of u/erica/AccretionModelingBlog
- Timestamp:
- 03/25/18 20:42:23 (7 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
u/erica/AccretionModelingBlog
v28 v29 32 32 == Research Log: == 33 33 34 === 3/18/18, $\gamma=7/5$, dynamical cases === 35 36 ==== Resolving sonic point ==== 34 === 3/18/18, $\gamma=7/5$ case=== 37 35 38 36 Am running the $\gamma=7/5$ case because the sonic point in this case is located far beyond the accretion radius ($R_s=6>>R_{acc}=.625$). This will enable us to see the generation of accretion shocks in the sim and check the behavior of the accretion algorithm under conditions of supersonic infall. We expect that the accretion algorithm should be fairly 'well-behaved' in the supersonic regime, as any pressure errors generated in the accretion kernel will be unable to propagate upstream and effect the hydrodynamical solution beyond the accretion volume. 39 37 40 For this run, I made the box bigger than the fiducial case above, but came the same effective resolution. The bigger box now contains the bondi radius ($R_{BH}=30.5$). The following table lists this sims params:38 For this run, I made the box bigger than the fiducial case above, but used the same effective resolution. The bigger box now encompasses the bondi radius ($R_{BH}=30.5$). The following table lists this sims params: 41 39 42 40 || $\gamma$ || 1.4 || … … 52 50 53 51 The final time of the simulation is ~2 sound crossing times, using the sound speed of the gas at the edge of the domain ($C_{amb}=1.283$) and the distance between this edge and the sink particle ($r=40$). 52 53 This case no longer forces an inner spherical boundary in the sim where the Bondi solution is copied. Instead the solution is sampled for each cell within the computational volume (including ghost zones) at t=0. For each subsequent cycle, the Bondi solution is pasted everywhere $r>20$. 54 55 ==== Results ==== 56 57 The non-dimensional profiles for this run line up nicely with Bondi's Fig. 3 (curves III and IV): 58 59 [[Image()]] 60 61 || Quantity || Plot || Notes || 62 || 63 54 64 55 65