diff options
| -rw-r--r-- | fg21sim/extragalactic/clusters/halo.py | 19 | ||||
| -rw-r--r-- | fg21sim/extragalactic/clusters/main.py | 2 |
2 files changed, 10 insertions, 11 deletions
diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py index 4cf6306..2b7994a 100644 --- a/fg21sim/extragalactic/clusters/halo.py +++ b/fg21sim/extragalactic/clusters/halo.py @@ -197,7 +197,7 @@ class RadioHalo1M: """ return self.age_merger - def time_turbulence(self, t=None): + def duration_turb(self, t=None): """ The duration that the turbulence persists strong enough to be able to effectively accelerate the electrons, which is @@ -209,14 +209,13 @@ class RadioHalo1M: Unit: [Gyr] """ t_merger = self._merger_time(t) - 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) + M_main = self.mass_main(t=t_merger) + M_sub = self.mass_sub(t=t_merger) L_turb = 2 * self.radius_turbulence(t_merger) + vi = helper.velocity_impact(M_main, M_sub, z_merger) uconv = AUC.kpc2km * AUC.s2Gyr # [kpc]/[km/s] => [Gyr] - time = uconv * 2*L_turb / vi # [Gyr] - return time + return uconv * 2*L_turb / vi # [Gyr] def mach_turbulence(self, t=None): """ @@ -717,8 +716,8 @@ class RadioHalo1M: 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) + tau_turb = self.duration_turb(t_merger) + return (t >= t_merger) and (t <= t_merger + tau_turb) def _energy_loss(self, gamma, t): """ @@ -884,7 +883,7 @@ class RadioHaloAM(RadioHalo1M): return (mass1 + rate * (t - t1)) @property - def time_turbulence_avg(self): + def duration_turb_avg(self): """ Calculate the time-averaged turbulence acceleration active time within the period from ``age_begin`` to ``age_obs``. @@ -893,7 +892,7 @@ class RadioHaloAM(RadioHalo1M): """ dt = self.time_step xt = np.arange(self.age_begin, self.age_obs+dt/2, step=dt) - t_turb = np.array([self.time_turbulence(t) for t in xt]) + t_turb = np.array([self.duration_turb(t) for t in xt]) return np.sum(t_turb * dt) / (len(xt) * dt) @property diff --git a/fg21sim/extragalactic/clusters/main.py b/fg21sim/extragalactic/clusters/main.py index 8affc6c..061b8ef 100644 --- a/fg21sim/extragalactic/clusters/main.py +++ b/fg21sim/extragalactic/clusters/main.py @@ -304,8 +304,8 @@ class GalaxyClusters: n_e=n_e, # [cm^-3] gamma=halo.gamma, # Lorentz factors Ke=halo.injection_rate, # [cm^-3 Gyr^-1] - time_turb=halo.time_turbulence_avg, # [Gyr] Mach_turb=halo.mach_turbulence_avg, # Mach number + tau_turb=halo.duration_turb_avg, # [Gyr] tau_acc=halo.tau_acceleration_avg, # [Gyr] tfrac_acc=halo.time_acceleration_fraction, ) |