From 8ab2148f9871920e81562a88f491456fcca4b8b5 Mon Sep 17 00:00:00 2001 From: Aaron LI Date: Fri, 25 Jan 2019 15:49:05 +0800 Subject: clusters/halo: Fold parameter "zeta_ins" into "f_acc" --- fg21sim/configs/config.spec | 10 +++------- fg21sim/extragalactic/clusters/halo.py | 5 ++--- 2 files changed, 5 insertions(+), 10 deletions(-) (limited to 'fg21sim') diff --git a/fg21sim/configs/config.spec b/fg21sim/configs/config.spec index 1ba8dd3..adc4134 100644 --- a/fg21sim/configs/config.spec +++ b/fg21sim/configs/config.spec @@ -374,19 +374,15 @@ stream = option("stderr", "stdout", "", default="stderr") # Giant radio halos # [[halos]] - # A custom parameter to tune the turbulent acceleration efficiency. - # NOTE: The smaller this parameter, the shorter the acceleration - # timescale, therefore more efficient acceleration. + # A custom factor to tune the turbulent acceleration efficiency. + # NOTE: This parameter incorporates the efficiency factor describing + # the effectiveness of the ICM plasma instabilities. f_acc = float(default=1.0, min=0.1, max=10) # The factor that is multiplied to the turbulence injection radius # to derive the radio halo radius. f_radius = float(default=1.0, min=0.5, max=10) - # An efficiency factor describing the effectiveness of plasma - # instabilities (e.g., due to spatial or temporal intermittency). - zeta_ins = float(default=0.1, min=0.1, max=1.0) - # The fraction of merger energy transferred into the turbulence. eta_turb = float(default=0.1, min=0.1, max=0.5) diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py index 81b9c0d..4e87f3f 100644 --- a/fg21sim/extragalactic/clusters/halo.py +++ b/fg21sim/extragalactic/clusters/halo.py @@ -151,7 +151,6 @@ class RadioHalo1M: self.configs = configs self.f_acc = configs.getn(comp+"/f_acc") self.f_radius = configs.getn(comp+"/f_radius") - self.zeta_ins = configs.getn(comp+"/zeta_ins") self.eta_turb = configs.getn(comp+"/eta_turb") self.eta_e = configs.getn(comp+"/eta_e") self.x_cr = configs.getn(comp+"/x_cr") @@ -334,9 +333,9 @@ class RadioHalo1M: k_L = 2 * np.pi / L_turb cs = helper.speed_sound(self.kT(t_merger)) # [km/s] v_t = self._velocity_turb(t_merger) # [km/s] - tau = self.x_cr * cs**3 / (8*k_L * self.zeta_ins * v_t**4) + tau = self.x_cr * cs**3 / (8*k_L * v_t**4) tau *= AUC.s2Gyr * AUC.kpc2km # [s kpc/km] -> [Gyr] - tau *= self.f_acc # custom tune parameter + tau /= self.f_acc # tune factor (folded with "zeta_ins") # Impose the maximum acceleration timescale if tau > tau_max: -- cgit v1.2.2