aboutsummaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
authorAaron LI <aly@aaronly.me>2019-01-26 22:14:45 +0800
committerAaron LI <aly@aaronly.me>2019-01-26 22:14:45 +0800
commit7eca59bde1b32b0f10cd8c2dad9c760838469709 (patch)
tree9b4f32e0b966fc2582c720d015839a93363aa752
parentc9f0dbf07992a67591881e5a092b5366b2611939 (diff)
downloadfg21sim-7eca59bde1b32b0f10cd8c2dad9c760838469709.tar.bz2
clusters/halo: Tweak {duration,velocity}_turb()
-rw-r--r--fg21sim/extragalactic/clusters/halo.py15
1 files changed, 7 insertions, 8 deletions
diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py
index 33998de..a0c49a2 100644
--- a/fg21sim/extragalactic/clusters/halo.py
+++ b/fg21sim/extragalactic/clusters/halo.py
@@ -238,10 +238,9 @@ class RadioHalo1M:
@lru_cache()
def duration_turb(self, t_merger):
"""
- The duration that the turbulence persists strong enough to be
- able to effectively accelerate the electrons, which is
- estimated as:
- τ_turb ≅ 2*L / v_impact = 4*R_turb / v_impact.
+ The duration that the turbulence persists strong enough to be able
+ to effectively accelerate the electrons, which is estimated as:
+ τ_turb ~ d / v_impact ~ 2*R_turb / v_impact.
Reference: [miniati2015],Sec.5
@@ -251,10 +250,10 @@ class RadioHalo1M:
z_merger = COSMO.redshift(t_merger)
M_main = self.mass_main(t=t_merger)
M_sub = self.mass_sub(t=t_merger)
- L_turb = 2 * self.radius_turb(t_merger)
- vi = helper.velocity_impact(M_main, M_sub, z_merger)
+ d = 2 * self.radius_turb(t_merger)
+ v_i = helper.velocity_impact(M_main, M_sub, z_merger)
uconv = AUC.kpc2km * AUC.s2Gyr # [kpc]/[km/s] => [Gyr]
- return uconv * 2*L_turb / vi # [Gyr]
+ return uconv * d / v_i # [Gyr]
@lru_cache()
def velocity_turb(self, t_merger):
@@ -296,7 +295,7 @@ class RadioHalo1M:
R_turb = self.radius_turb(t_merger) # [kpc]
rho_gas_f = helper.calc_gas_density_profile(
- M_main, z, f_rc=self.f_rc, beta=self.beta)
+ M_main+M_sub, z, f_rc=self.f_rc, beta=self.beta)
M_turb = 4*np.pi * integrate.quad(
lambda r: rho_gas_f(r) * r**2,
a=0, b=R_turb)[0] # [Msun]