Simulation Log

[bottom]

17.02.23

  • Develop WIND Tracers
model k_sf alpha max em AGN FB Mech FB Rad FB etc
B05v_v2 0 0 0.1 o o o same as B05v but wind tracers
agnRot_ac1_v2 0 0 0.3 o o o same as agnRot_ac1 but wind tracers
agnRot_ac1_v2_em02 0 0 0.2 o o o
agnRot_ac1_v2_nofb 0 0 0.3 x x x
agnRot_ac1_v2_norad 0 0 0.3 o o x
agnRot_ac1_v2_nomech 0 0 0.3 o x o
Rotation case
agnRot_ac5_v2 0.5 0.1 0.3 o o o same as agnRot_ac5 but wind tracers
agnRot_ac5_v2_em02 0.5 0.1 0.2 o o o
agnRot_ac5_v2_alp-2 0.5 0.01 0.3 o o o
agnRot_ac5_v2_nofb 0.5 0.1 0.3 x x x
agnRot_ac5_v2_norad 0.5 0.1 0.3 o o x
agnRot_ac5_v2_nomech 0.5 0.1 0.3 o x o
agnRot_ac5_v2_k01 0.1 0.1 0.3 o o o
agnRot_ac5_v2_k03 0.3 0.1 0.3 o o o
agnRot_ac5_v2_k07 0.7 0.1 0.3 o o o
agnRot_ac5_v2_k09 0.9 0.1 0.3 o o o
17.02.17
model k_sf AGN FB SN etc
agnRot_ac6_k01 0.1 x o inherited by agnRot_ac5
agnRot_ac6_k03 0.3 x o inherited by agnRot_ac5
agnRot_ac6_k05 0.5 x o inherited by agnRot_ac5
agnRot_ac6_k07 0.7 x o inherited by agnRot_ac5
agnRot_ac7_k01 0.1 x x inherited by agnRot_ac5
agnRot_ac7_k03 0.3 x x inherited by agnRot_ac5
agnRot_ac7_k05 0.5 x x inherited by agnRot_ac5
agnRot_ac7_k07 0.7 x x inherited by agnRot_ac5
agnRot_ac8 0.5 o o agnRot_ac5 + alpha=0.01
16.12.18
Update the code with new physics including the treatment of hot accretion flow
model etc
agnRot_ac1 update for hot accretion flow otherwise same with B05v (Gan+14)
agnRot_ac2 only change for Tc (compton temperature) otherwise same with B05v
agnRot_ac3 fixed Tc ($=2.5\times 10^{7}\,K$) otherwise same with B05v
17.01.02
agnRot_ac4 agnRot_ac1 + rotating & viscosity (ssa=0.1, k_sf=0.5)
agnRot_ac5 agnRot_ac4 + vph_src * sin(x2b)
17.01.12
agnRot_ac4_k01 agnRot_ac4 + k_sf=0.1
agnRot_ac4_k03 agnRot_ac4 + k_sf=0.3
agnRot_ac4_k07 agnRot_ac4 + k_sf=0.7
agnRot_ac5_k01 agnRot_ac5 + k_sf=0.1
agnRot_ac5_k03 agnRot_ac5 + k_sf=0.3
agnRot_ac5_k07 agnRot_ac4 + k_sf=0.7
16.10.11
Clean up the code (no theta dependency of vph_src)
model ssa k_sf Cool SF StarRad SN AGN etc
agnRot_aa1 0.1 0.5 o o o o o config. is inherited from y24
agnRot_aa2 0.1 0.5 o o o o o config. is inherited from y24_4
16.10.22
agnRot_aa3 0.3 0.5 o o o o o
agnRot_aa4 0.5 0.5 o o o o o
agnRot_aa4_2 0.5 0.5 o o o o o Initial stellar age: 5 Gyr (not 2Gyr)
agnRot_aa5 0.1 0.3 o o o o o
agnRot_aa6 0.3 0.3 o o o o o
agnRot_aa7 0.5 0.3 o o o o o
16.11.09
agnRot_aa8 0.1 0.1 o o o o o
16.11.02
agnRot_ab1 0.1 0.5 o o o o x
agnRot_ab2 0.1 0.3 o o o o x
agnRot_ab3 0.1 0.1 o o o o x
16.10.07
vph_src to theta dependent
model ssa k_sf Cool SF StarRad SN AGN etc
agnRot_z01 0.1 0.5 o o o o x config. is inherited from y21
agnRot_z02 0.1 0.5 o o o o o config. is inherited from y24
16.09.29
finding appropriate k_sf (stellar rotation)
model ssa k_sf Cool SF StarRad SN AGN etc
agnRot_y41 0 0.1 o o o o x
agnRot_y42 0 0.1 o o o o x
agnRot_y43 0 0.1 o o o o x
agnRot_y44 0 0.1 o o o o x
agnRot_y45 0 0.1 o o o o x
agnRot_y46 0 0.1 o o o o x
agnRot_y47 0 0.1 o o o o x
agnRot_y48 0 0.1 o o o o x
agnRot_y49 0 0.1 o o o o x
16.09.17
model ssa Cool SF StarRad SN AGN dtvisc dtvisc2 etc
agnRot_y01 1e-1 o o o o x 0.1 1.e-3 dtvisc = max( min(dt*dtvisc, dtviscmin), dt*dtvisc2 )
v_phi > 0
agnRot_y02 1e-1 o o o o x 0.1 1.e-4
agnRot_y03 1e-1 o o o o x 0.1 1.e-5
agnRot_y04 1e-1 o o o o o 0.1 1.e-3
agnRot_y05 1e-1 o o o o o 0.1 1.e-4
agnRot_y06 1e-1 o o o o o 0.1 1.e-5
agnRot_y07 1e-1 o o o o x 0.1 1.e-3 averaged change rate
agnRot_y08 1e-1 o o o o x 0.1 1.e-4
agnRot_y09 1e-1 o o o o x 0.1 1.e-5
agnRot_y10 1e-1 o o o o o 0.1 1.e-3
agnRot_y11 1e-1 o o o o o 0.1 1.e-4
agnRot_y12 1e-1 o o o o o 0.1 1.e-5
16.09.20
modify viscosity time scale for viscosity heating term otherwise overheated
model ssa Cool SF StarRad SN AGN dtvisc dtvisc2 etc
agnRot_y21 1e-1 o o o o x 0.1 1.e-2 same as agnRot_y01
agnRot_y22 1e-1 o o o o x 0.1 1.e-3
agnRot_y23 1e-1 o o o o x 0.1 1.e-4
agnRot_y24 1e-1 o o o o o 0.1 1.e-2
agnRot_y25 1e-1 o o o o o 0.1 1.e-3
agnRot_y26 1e-1 o o o o o 0.1 1.e-4
agnRot_y24_2 1e-1 o o o o o 0.1 1.e-2 only negative chgrate
agnRot_y24_3 1e-1 o o o o o 0.1 1.e-2 no viscosity heating
agnRot_y24_4 1e-1 o o o o o 0.1 1.e-2 heating only if negative chgrate
agnRot_y31 1e-1 o o o o x 0.1 1.e-2 same as agnRot_y07
agnRot_y32 1e-1 o o o o x 0.1 1.e-3
agnRot_y33 1e-1 o o o o x 0.1 1.e-4
agnRot_y34 1e-1 o o o o o 0.1 1.e-2
agnRot_y35 1e-1 o o o o o 0.1 1.e-3
agnRot_y36 1e-1 o o o o o 0.1 1.e-4
16.09.27
model ssa Cool SF StarRad SN AGN dtvisc dtvisc2 etc
agnRot_y27 5e-2 o o o o o 0.1 1.e-3 viscosity approx. version 1. (first order)
agnRot_y28 3e-1 o o o o o 0.1 1.e-3
agnRot_y29 5e-1 o o o o o 0.1 1.e-3
16.08.24
model ssa Cool SF StarRad SN AGN etc
agnRot_s01 1e-3 o x o o x
agnRot_s02 1e-2 o x o o x
agnRot_s03 1e-1 o x o o x
agnRot_t01 1e-3 o o o o x no dstar within 500pc otherwise same as agnRot_m??
To avoid code crash, viscosity starts to work after 0.025 t_code
agnRot_t02 1e-2 o o o o x
agnRot_t03 1e-1 o o o o x
16.08.29
no dstar within 500pc
model ssa Cool SF StarRad SN AGN etc
agnRot_u01 0 x x x x x
agnRot_u02 1e-3 x x x x x
agnRot_u03 1e-2 x x x x x
agnRot_u04 5e-2 x x x x x
agnRot_v01 0 x x x x x no dark-matter only BH gravity (mbh=12 Code unit)
agnRot_v02 1e-3 x x x x x
agnRot_v03 1e-2 x x x x x
agnRot_v04 5e-2 x x x x x
agnRot_v05 1e-3 x x x x x mbh=12e2 otherwise same with agnRot_v02
agnRot_v06 1e-2 x x x x x mbh=12e2 otherwise same with agnRot_v03
agnRot_v07 5e-2 x x x x x mbh=12e2 otherwise same with agnRot_v04
agnRot_w01 1e-3 o o x x x mbh=12e2 otherwise same with agnRot_v??
agnRot_w02 1e-2 o o x x x
agnRot_w03 5e-2 o o x x x
agnRot_w04 1e-3 o o x x x mbh=12e0 & include dark-matter otherwise same with agnRot_w02
agnRot_w05 1e-2 o o x x x
agnRot_w06 5e-2 o o x x x
mbh=12e2 & no dark-matter
model ssa Cool SF StarRad SN AGN etc
agnRot_w07 1e-3 o o o o x
agnRot_w08 1e-2 o o o o x
agnRot_w09 5e-2 o o o o x
16.09.02
model ssa Cool SF StarRad SN AGN etc
agnRot_x01 5e-1 o o o o o
agnRot_x02 5e-1 o o o o x
agnRot_x03 1e-1 o o o o o
16.09.06
model ssa Cool SF StarRad SN AGN etc
agnRot_x04 2e-1 o o o o o
agnRot_x05 3e-1 o o o o o
agnRot_x06 5e-2 o o o o o
agnRot_x07 1e-2 o o o o o
16.07.22
black hole accretion test only considering black hole point mass
zoomed-in to $10^{-2} pc < r < 4 pc$
model ssa etc
bhacc_a01 0 only HSE non-rotating
bhacc_a02 0 HSE background w/ Torus-like high density-rotating ring
16.07.22
model ssa etc
bhacc_a03 1e-3 viscosity using direct calculation
bhacc_a04 1e-2
bhacc_a05 2e-2
bhacc_a06 3e-2
bhacc_a07 4e-2
bhacc_a08 5e-2
16.07.27
model ssa etc
bhacc_b01 1e-3 $\nu \propto r^{1/2}$
bhacc_b02 1e-2
bhacc_b03 2e-2
bhacc_b04 3e-2
bhacc_b05 4e-2
bhacc_b06 5e-2
16.08.10
continuous injection in the torus
model ssa etc
bhacc_c01 0 rest set is same as bhacc_a?? above
bhacc_c02 1e-3
bhacc_c03 1e-2
bhacc_c04 2e-2
bhacc_c05 3e-2
bhacc_c06 4e-2
bhacc_c06 5e-2
bhacc_d01 1e-3 rest set is same as bhacc_b?? above
bhacc_d02 1e-2
bhacc_d03 2e-2
bhacc_d04 3e-2
bhacc_d05 4e-2
bhacc_d06 5e-2
16.08.17
torus radius to $10^{-3} pc < r < 3\times10^{-3} pc$
$r_{sh} = 2.87\times 10^{-5} pc$
model ssa etc
bhacc_e01 0 rest set is same as bhacc_a?? above
bhacc_e02 1e-3
bhacc_e03 1e-2
bhacc_e04 2e-2
bhacc_e05 3e-2
bhacc_e06 4e-2
bhacc_f01 1e-3 rest set is same as bhacc_b?? above
bhacc_f02 1e-2
bhacc_f03 2e-2
bhacc_f04 3e-2
bhacc_f05 4e-2
16.08.29
gas supply at $r > 10^{-3} pc$
model ssa etc
bhacc_g01 0 rest set is same as bhacc_a?? above
bhacc_g02 1e-3
bhacc_g03 1e-2
bhacc_g04 2e-2
bhacc_g05 3e-2
bhacc_h01 0 star formation & cooling otherwise same with bhacc_g01
16.07.21
to find the reason of inconsistency bet. $alpha$ and $v_{r}$
model ssa Cool&SF StarRad SN AGN etc
agnRot_r01 0 o o x x
agnRot_r02 1e-3 o o x x
agnRot_r03 1e-2 o o x x
agnRot_r04 1e-1 o o x x
agnRot_q01 0 o x o x
agnRot_q02 1e-3 o x o x
agnRot_q03 1e-2 o x o x
agnRot_q04 1e-1 o x o x
agnRot_p01 0 o x x x
agnRot_p02 1e-3 o x x x
agnRot_p03 1e-2 o x x x
agnRot_p04 1e-1 o x x x
16.07.18
model ssa Cool&SF StarRad SN AGN etc
agnRot_o01 1e-3 o o x o
agnRot_o02 1e-2 o o x o
agnRot_o03 1e-1 o o x o
agnRot_o04 1 o o x o
agnRot_o05 10 o o x o
agnRot_o06 100 o o x o
16.07.13
model ssa Cool&SF StarRad SN AGN etc
agnRot_n01 1e-3 o o o o
agnRot_n02 1e-2 o o o o
agnRot_n03 1e-1 o o o o
agnRot_n04 1 o o o o
agnRot_n05 10 o o o o
agnRot_n06 100 o o o o
16.07.08
model ssa Cool&SF StarRad SN AGN etc
agnRot_m01 1e-3 o o o x dtviscmin: from direct change rate
agnRot_m02 1e-2 o o o x
agnRot_m03 1e-1 o o o x
agnRot_m04 1 o o o x
agnRot_m05 10 o o o x
agnRot_m06 100 o o o x
16.08.24
agnRot_m07 5e-1 o o o x
agnRot_m08 3e-1 o o o x
agnRot_m09 2e-1 o o o x
16.07.06
test run after fixing the code crash
model ssa Cool&SF StarRad SN AGN etc
agnRot_l01 1e-3 o o o x dtviscmin: diffuse time scale
before 16.07
model ssa Cool&SF StarRad SN AGN etc
agnRot_c01 0 o o o x
agnRot_b01 0 x x o x
agnRot_a01 0 x x x x

16.07.01

  • Fixing the MPI miscommunication problem in alpha_old.F

16.06.30

  • Finding the proper alpha parameter showing the reasonable inflow which is consistent with Hopkins et al. 2011
  • Running w/ Single core
model ssa etc
agnRot_i01 1e-3 $\rho_{0} = 10^{-10}$ and cool&viscosity on
agnRot_i02 1e-2
agnRot_i03 1e-1
agnRot_i04 1

16.06.29

model ssa etc crash
agnRot_e01 1e-3 everything but AGN Feedback & d0=1d-10 62
agnRot_e02 1e-2 everything but AGN Feedback 180
agnRot_e03 1e-1 everything but AGN Feedback
agnRot_e04 1 everything but AGN Feedback
agnRot_g01 0 everything but AGN Feedback no crash
agnRot_h01 1e-3 same w/ agnRot_e01 but d0 becomes 1/10 (d0=1d-11) of it to test of possibility
that the code crashed by thermal instability
crash
  • fixing the code-crash
model x1min nr / nth etc
test 2.5e-3 120 / 30 $r1 \equiv x1min$
test2 2.5e-2 96 / 26 10 r1
test3 5e-3 112 / 24 2 r1
test4 1e-2 104 / 22 4 r1
  • try w/ single core
model etc
test7 $\rho_{0} = 10^{-11}$