From c9b0744c1056ebcaf425c4ff7d6eb98786b40699 Mon Sep 17 00:00:00 2001 From: Aaron LI Date: Mon, 1 Jan 2018 16:57:16 +0800 Subject: clusters/halo: impose the maximum acceleration timescale --- fg21sim/extragalactic/clusters/halo.py | 12 +++++++++--- 1 file changed, 9 insertions(+), 3 deletions(-) diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py index a8f218a..1da89b3 100644 --- a/fg21sim/extragalactic/clusters/halo.py +++ b/fg21sim/extragalactic/clusters/halo.py @@ -1,4 +1,4 @@ -# Copyright (c) 2017 Weitian LI +# Copyright (c) 2017-2018 Weitian LI # MIT license """ @@ -746,12 +746,18 @@ class RadioHalo: ---------- Ref.[donnert2013],Eq.(15) """ + # Maximum acceleration timescale when no turbulence acceleration + # NOTE: see the above WARNING! + tau_max = 10.0 # [Gyr] if (t < self.age_merger) or (t > self.age_merger+self.time_turbulence): - # NO acceleration (see also the above NOTE and WARNING!) - tau_acc = 10 # [Gyr] + # NO active turbulence acceleration + tau_acc = tau_max else: # Turbulence acceleration tau_acc = self.tau_acceleration # [Gyr] + # Impose the maximum acceleration timescale + if tau_acc > tau_max: + tau_acc = tau_max gamma = np.asarray(gamma) diffusion = gamma**2 / 4 / tau_acc -- cgit v1.2.2