diff options
| author | Aaron LI <aly@aaronly.me> | 2019-01-25 15:49:05 +0800 |
|---|---|---|
| committer | Aaron LI <aly@aaronly.me> | 2019-01-25 15:49:05 +0800 |
| commit | 8ab2148f9871920e81562a88f491456fcca4b8b5 (patch) | |
| tree | 5124815c8cda1fef659527cc79ab5c9a4b4a552e | |
| parent | 427b12bb9b4f5b9923494efaf33c83cb213ab3d6 (diff) | |
| download | fg21sim-8ab2148f9871920e81562a88f491456fcca4b8b5.tar.bz2 | |
clusters/halo: Fold parameter "zeta_ins" into "f_acc"
| -rw-r--r-- | fg21sim/configs/config.spec | 10 | ||||
| -rw-r--r-- | fg21sim/extragalactic/clusters/halo.py | 5 |
2 files changed, 5 insertions, 10 deletions
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: |