diff options
Diffstat (limited to 'fg21sim')
| -rw-r--r-- | fg21sim/extragalactic/clusters/halo.py | 17 | ||||
| -rw-r--r-- | fg21sim/extragalactic/clusters/main.py | 2 |
2 files changed, 15 insertions, 4 deletions
diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py index 6605763..6b5a683 100644 --- a/fg21sim/extragalactic/clusters/halo.py +++ b/fg21sim/extragalactic/clusters/halo.py @@ -210,6 +210,15 @@ class RadioHalo: z=self.z_merger, configs=self.configs) @property + def mach_turbulence(self): + """ + The turbulence Mach number determined from its velocity dispersion. + """ + cs = helper.speed_sound(self.kT_main) # [km/s] + v_turb = self._velocity_turb() # [km/s] + return v_turb / cs + + @property def radius_virial_obs(self): """ The virial radius of the "current" cluster (``M_obs``) at @@ -342,7 +351,7 @@ class RadioHalo: R_vir = helper.radius_virial(mass=self.M_main, z=self.z_merger) L = self.f_lturb * R_vir # [kpc] cs = helper.speed_sound(self.kT_main) # [km/s] - v_turb = self._velocity_turb(t=self.age_merger) # [km/s] + v_turb = self._velocity_turb() # [km/s] tau = (self.x_cr * cs**3 * L / (16*np.pi * self.zeta_ins * v_turb**4)) # [s kpc/km] tau *= AUC.s2Gyr * AUC.kpc2km # [Gyr] @@ -805,7 +814,7 @@ class RadioHalo: mass = rate * (t - t_merger) + self.M_main return mass - def _velocity_turb(self, t): + def _velocity_turb(self, t=None): """ Calculate the turbulence velocity dispersion (i.e., turbulence Mach number). @@ -834,8 +843,10 @@ class RadioHalo: The turbulence velocity dispersion Unit: [km/s] """ - mass = self.M_main + self.M_sub + if t is None: + t = self.age_merger z = COSMO.redshift(t) + mass = self.M_main + self.M_sub R_vir = helper.radius_virial(mass=mass, z=z) * AUC.kpc2cm # [cm] v2 = np.sqrt(AC.G * self.M_sub*AUC.Msun2g / R_vir) # [cm/s] v2 *= AUC.cm2km # [km/s] diff --git a/fg21sim/extragalactic/clusters/main.py b/fg21sim/extragalactic/clusters/main.py index c441977..57811af 100644 --- a/fg21sim/extragalactic/clusters/main.py +++ b/fg21sim/extragalactic/clusters/main.py @@ -293,7 +293,7 @@ class GalaxyClusters: ("Rhalo", halo.radius), # [kpc] ("Rhalo_angular", halo.angular_radius), # [arcsec] ("volume", halo.volume), # [kpc^3] - ("Mach_turb", halo.Mach_turbulence), # turbulence Mach number + ("Mach_turb", halo.mach_turbulence), # turbulence Mach number ("tau_acc", halo.tau_acceleration), # [Gyr] ("Ke", halo.injection_rate), # [cm^-3 Gyr^-1] ("gamma", halo.gamma), # Lorentz factors |