aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim/extragalactic/clusters/main.py
diff options
context:
space:
mode:
authorAaron LI <aly@aaronly.me>2017-10-21 22:13:42 +0800
committerAaron LI <aly@aaronly.me>2017-10-21 22:13:42 +0800
commit2ee32e7196f83a9350a2b83c03846f6282d7ffe9 (patch)
tree7171acc49e70b31ef796c4a58f262c0367a6972f /fg21sim/extragalactic/clusters/main.py
parenta8a02a96a78a85ffacf52aa43c158ecc24391af3 (diff)
downloadfg21sim-2ee32e7196f83a9350a2b83c03846f6282d7ffe9.tar.bz2
clusters/halo: Rewrite turbulence acceleration & diffusion coefficient!
The previous calculation of the turbulence acceleration efficiency has significant problems, which is too low, i.e., the acceleration timescale is too large (~10 Gyr). However, it is reported that turbulent acceleration has a timescale ~100 Myr (0.1 Gyr). I believe this problem is due to the over-simplification to the formula proposed by [cassano2005]. Rewrite the turbulence acceleration and diffusion coefficient calculations, adopting the method proposed by [brunetti2016]. To this end, two new options "f_lturb" and "f_acc" are introduced to tune the results.
Diffstat (limited to 'fg21sim/extragalactic/clusters/main.py')
-rw-r--r--fg21sim/extragalactic/clusters/main.py3
1 files changed, 2 insertions, 1 deletions
diff --git a/fg21sim/extragalactic/clusters/main.py b/fg21sim/extragalactic/clusters/main.py
index 506e77c..dfc8e6e 100644
--- a/fg21sim/extragalactic/clusters/main.py
+++ b/fg21sim/extragalactic/clusters/main.py
@@ -281,8 +281,9 @@ class GalaxyClusters:
("Rhalo", halo.radius), # [kpc]
("Rhalo_angular", halo.angular_radius), # [arcsec]
("volume", halo.volume), # [kpc^3]
+ ("Mach_turb", halo.Mach_turbulence), # turbulence Mach number
+ ("tau_acc", halo.tau_acceleration), # [Gyr]
("Ke", halo.injection_rate), # [cm^-3 Gyr^-1]
- ("chi", halo._chi_acceleration()), # [Gyr^-1]
("gamma", halo.gamma), # Lorentz factors
("n_e", n_e), # [cm^-3]
])