diff options
| author | Aaron LI <aly@aaronly.me> | 2019-01-26 22:14:45 +0800 |
|---|---|---|
| committer | Aaron LI <aly@aaronly.me> | 2019-01-26 22:14:45 +0800 |
| commit | 7eca59bde1b32b0f10cd8c2dad9c760838469709 (patch) | |
| tree | 9b4f32e0b966fc2582c720d015839a93363aa752 /fg21sim/extragalactic/clusters | |
| parent | c9f0dbf07992a67591881e5a092b5366b2611939 (diff) | |
| download | fg21sim-7eca59bde1b32b0f10cd8c2dad9c760838469709.tar.bz2 | |
clusters/halo: Tweak {duration,velocity}_turb()
Diffstat (limited to 'fg21sim/extragalactic/clusters')
| -rw-r--r-- | fg21sim/extragalactic/clusters/halo.py | 15 |
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] |