Changes between Version 28 and Version 29 of u/erica/AccretionModelingBlog


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

Legend:

Unmodified
Added
Removed
Modified
  • u/erica/AccretionModelingBlog

    v28 v29  
    3232== Research Log: ==
    3333
    34 === 3/18/18, $\gamma=7/5$, dynamical cases ===
    35 
    36 ==== Resolving sonic point ====
     34=== 3/18/18, $\gamma=7/5$ case===
    3735
    3836Am 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.
    3937
    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:
     38For 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:
    4139
    4240|| $\gamma$ || 1.4 ||
     
    5250
    5351The 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
     53This 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
     57The 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
    5464
    5565