From 211cad9fc0db7b14908fda113a5e3df7cd31676f Mon Sep 17 00:00:00 2001 From: Aaron LI Date: Fri, 18 Jan 2019 15:47:55 +0800 Subject: clusters/halo: Update time_turbulence() calculation --- fg21sim/extragalactic/clusters/halo.py | 17 +++++++++-------- 1 file changed, 9 insertions(+), 8 deletions(-) (limited to 'fg21sim/extragalactic/clusters') diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py index 3984821..99f2306 100644 --- a/fg21sim/extragalactic/clusters/halo.py +++ b/fg21sim/extragalactic/clusters/halo.py @@ -199,11 +199,12 @@ class RadioHalo1M: def time_turbulence(self, t=None): """ - The duration that the compressive turbulence persists, which is + The duration that the turbulence persists strong enough to be + able to effectively accelerate the electrons, which is estimated as: - τ_turb ≅ L / v_impact = 2*R_turb / v_impact. - During this period, the merger-induced turbulence is regarded - to accelerate the relativistic electrons effectively. + τ_turb ≅ 2*L / v_impact = 4*R_turb / v_impact. + + Reference: [miniati2015],Sec.5 Unit: [Gyr] """ @@ -211,10 +212,10 @@ class RadioHalo1M: mass_main = self.mass_main(t=t_merger) mass_sub = self.mass_sub(t=t_merger) z_merger = COSMO.redshift(t_merger) - vi = helper.velocity_impact(mass_main, mass_sub, z_merger) # [km/s] - distance = 2 * self.radius_turbulence(t_merger) - uconv = AUC.kpc2km * AUC.s2Gyr # [s kpc/km] => [Gyr] - time = uconv * distance / vi # [Gyr] + vi = helper.velocity_impact(mass_main, mass_sub, z_merger) + L_turb = 2 * self.radius_turbulence(t_merger) + uconv = AUC.kpc2km * AUC.s2Gyr # [kpc]/[km/s] => [Gyr] + time = uconv * 2*L_turb / vi # [Gyr] return time def mach_turbulence(self, t=None): -- cgit v1.2.2