From 7e3e79ccf17ac8620b1cb784b9018760ceee9733 Mon Sep 17 00:00:00 2001
From: Aaron LI <aly@aaronly.me>
Date: Wed, 31 Oct 2018 15:44:30 +0800
Subject: clusters/halo: Rework the way to disable turbulent acceleration
Introduce a new property '_acceleration_disabled' to control whether the
turbulent acceleration is disabled. With this new property, improve the
code to derive the initial electron spectrum.
---
fg21sim/extragalactic/clusters/halo.py | 29 +++++++++++++----------------
1 file changed, 13 insertions(+), 16 deletions(-)
(limited to 'fg21sim')
diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py
index c54bb73..f18e16a 100644
--- a/fg21sim/extragalactic/clusters/halo.py
+++ b/fg21sim/extragalactic/clusters/halo.py
@@ -130,6 +130,8 @@ class RadioHalo:
the final time (``zend``), which is set by the methods
``self.calc_electron_spectrum()`` or ``self.set_electron_spectrum()``.
Unit: [cm^-3]
+ _acceleration_disabled : bool
+ Whether the turbulence acceleration is intentionally disabled?
"""
# Component name
compID = "extragalactic/halos"
@@ -145,6 +147,7 @@ class RadioHalo:
self.z_merger = z_merger
self.age_merger = COSMO.age(z_merger)
+ self._acceleration_disabled = False
self._set_configs(configs)
self._set_solver()
@@ -426,14 +429,15 @@ class RadioHalo:
n0_e = n_inj * (self.age_begin - self.time_init)
logger.debug("Derive the initial electron spectrum ...")
- dt = self.time_step
- tstart = self.age_begin - self.time_init - dt
- tstop = self.age_begin - dt # avoid acceleration at the ``age_begin``
- # Use a bigger time step to save time
- self.fpsolver.tstep = 3 * dt
+ self._acceleration_disabled = True
+ dt = self.fpsolver.tstep
+ tstart = self.age_begin
+ tstop = self.age_begin + self.time_init
+ self.fpsolver.tstep = 3 * dt # Bigger step to save time
n_e = self.fpsolver.solve(u0=n0_e, tstart=tstart, tstop=tstop)
- # Restore the original time step
self.fpsolver.tstep = dt
+ self._acceleration_disabled = False
+
return n_e
def calc_electron_spectrum(self, tstart=None, tstop=None, n0_e=None):
@@ -715,18 +719,11 @@ class RadioHalo:
def _is_turb_active(self, t):
"""
- Is the turbulence acceleration is active at the given (cosmic) time?
-
- NOTE
- ----
- Considering that the turbulence acceleration is a 2nd-order Fermi
- process, it has only an effective acceleration time ~<1 Gyr.
- Therefore, only during the period that strong turbulence persists
- in the ICM that the turbulence could effectively accelerate the
- relativistic electrons.
+ Is the turbulence acceleration is active at the given time?
"""
- if t < self.age_begin:
+ if self._acceleration_disabled:
return False
+
t_merger = self._merger_time(t)
t_turb = self.time_turbulence(t_merger)
return (t >= t_merger) and (t <= t_merger + t_turb)
--
cgit v1.2.2