Commit message (Collapse) | Author | Age | Files | Lines | |
---|---|---|---|---|---|
* | clusters/configs: Update f_acc default to 1.0; also update descriptions | Aaron LI | 2017-10-26 | 1 | -1/+1 |
| | |||||
* | clusters/halo: Calc turbulence injection scale w.r.t. sub-cluster radius | Aaron LI | 2017-10-26 | 1 | -6/+8 |
| | |||||
* | clusters/solver: Update buffer_np description and handling | Aaron LI | 2017-10-25 | 1 | -0/+2 |
| | |||||
* | clusters/halo: Fix fp_advection() w.r.t. initial electron spectrum calc. | Aaron LI | 2017-10-25 | 1 | -2/+4 |
| | |||||
* | clusters/halo: Adjust initial electron spectrum calculation | Aaron LI | 2017-10-25 | 1 | -4/+6 |
| | |||||
* | clusters/halo: Add method to derive the initial electron spectrum | Aaron LI | 2017-10-25 | 1 | -11/+51 |
| | | | | | Also add the option "time_init" to control how long a period is used to derive the initial electron spectrum. | ||||
* | clusters/halo: Fix passing the "configs" object | Aaron LI | 2017-10-25 | 1 | -3/+5 |
| | |||||
* | clusters/halo: Always pass the using configs | Aaron LI | 2017-10-24 | 1 | -3/+6 |
| | | | | | The user-using `configs` object may be different to the global `CONFIGS`, which may be forgotten to load the up-to-date user configurations. | ||||
* | clusters/halo: Change zbegin/zend to tstart/tsop (use cosmic time) | Aaron LI | 2017-10-24 | 1 | -19/+20 |
| | |||||
* | clusters/halo: Add method "_magnetic_field()" | Aaron LI | 2017-10-24 | 1 | -2/+23 |
| | |||||
* | clusters/halo: Rename property "magnetic_field" to "B_obs" | Aaron LI | 2017-10-24 | 1 | -3/+3 |
| | |||||
* | clusters/halo: Allow specify magnetic field in calc_emissivity() | Aaron LI | 2017-10-24 | 1 | -9/+14 |
| | |||||
* | clusters/halo: Use "lru_cache()", also improves some performance | Aaron LI | 2017-10-23 | 1 | -12/+22 |
| | |||||
* | clusters/halo: Add option "kT_out" to tune cluster temperature | Aaron LI | 2017-10-22 | 1 | -3/+1 |
| | |||||
* | clusters/halo: Rewrite turbulence acceleration & diffusion coefficient! | Aaron LI | 2017-10-21 | 1 | -57/+90 |
| | | | | | | | | | | | | The previous calculation of the turbulence acceleration efficiency has significant problems, which is too low, i.e., the acceleration timescale is too large (~10 Gyr). However, it is reported that turbulent acceleration has a timescale ~100 Myr (0.1 Gyr). I believe this problem is due to the over-simplification to the formula proposed by [cassano2005]. Rewrite the turbulence acceleration and diffusion coefficient calculations, adopting the method proposed by [brunetti2016]. To this end, two new options "f_lturb" and "f_acc" are introduced to tune the results. | ||||
* | clusters/halos: Add more cluster and halo properties | Aaron LI | 2017-10-21 | 1 | -7/+30 |
| | | | | | | * Rvir0, Rvir_main, Rvir_sub * kT_main, kT_sub (replace kT_merger) * radius -> Rhalo, angular_radius -> Rhalo_angular, B -> B0 | ||||
* | clusters/halo: Determine magnetic field by fraction of its energy density | Aaron LI | 2017-10-20 | 1 | -2/+2 |
| | | | | | | | | | | * Rewrite "magnetic_field()" function to calculate the mean magnetic field strength within ICM according its energy density fraction w.r.t. the ICM thermal energy density. * Remove config options "b_mean" and "b_index", which are replaced with the option "eta_b", the assumed magnetic energy density fraction w.r.t. the ICM thermal energy density. | ||||
* | clusters/halo: Add "kT_obs" and "tback_merger" properties | Aaron LI | 2017-10-19 | 1 | -1/+16 |
| | |||||
* | clusters/halo: Use "kT_cluster()" instead of "mass_to_kT()" | Aaron LI | 2017-10-19 | 1 | -1/+1 |
| | |||||
* | clusters/halo: rewrite power/flux/Tb calc with k-correction considered | Aaron LI | 2017-10-19 | 1 | -25/+94 |
| | |||||
* | clusters/halo: Add set_electron_spectrum() method; update doc a bit | Aaron LI | 2017-10-19 | 1 | -1/+26 |
| | |||||
* | clusters/halo: Improve time step determination a bit | Aaron LI | 2017-10-19 | 1 | -2/+5 |
| | |||||
* | clusters/halo: Increase the time step, but require a minimum of 20 steps | Aaron LI | 2017-10-18 | 1 | -0/+7 |
| | |||||
* | clusters: Update radio halo radius estimation | Aaron LI | 2017-10-18 | 1 | -4/+8 |
| | | | | | | | | | The simulated radio halo is assumed to have a radius of the falling sub-cluster; while previously it is assumed to be 1/4 of the virial radius of the merged cluster. The new estimation can agree better with the currently observed radio halos, which generally have a angular diameter size of ~2-7 [arcmin]. | ||||
* | clusters: Update the documents/comments a b it | Aaron LI | 2017-10-14 | 1 | -2/+2 |
| | |||||
* | clusters/halo.py: Update injection rate calculation | Aaron LI | 2017-08-26 | 1 | -10/+5 |
| | | | | | The electrons are assumed to be injected throughout the whole cluster ICM/volume. | ||||
* | clusters/halo.py: Rewrite electron acceleration coefficient calc | Aaron LI | 2017-08-14 | 1 | -26/+51 |
| | | | | | | | | | | * Adopt the electron acceleration coefficient formula from [cassano2005] * Rename method "_tau_acceleration()" to "_chi_acceleration()", and rewrite * Add property "kT_merger" * Also save "kT_merger" and "chi" into halos data Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | clusters/halo.py: Simplify two coefficients usage | Aaron LI | 2017-08-14 | 1 | -7/+5 |
| | | | | Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | clusters: Rename property "injection_rate" and save into halos data | Aaron LI | 2017-08-14 | 1 | -37/+37 |
| | | | | Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | clusters: Implement calc_{power,flux,brightness_mean} in helper | Aaron LI | 2017-08-11 | 1 | -40/+10 |
| | | | | | | Help development & debug ... Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | clusters/halo: "volume" use unit "kpc^3" instead of "cm^3" | Aaron LI | 2017-08-11 | 1 | -3/+6 |
| | | | | Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | Fix a bug in clusters/halo.py; Add a comment to utils/io/_create_dir() | Aaron LI | 2017-08-06 | 1 | -1/+1 |
| | | | | Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | clusters/halo.py: Finished brief simulation description | Aaron LI | 2017-08-01 | 1 | -3/+13 |
| | |||||
* | clusters: Fix several small issues on wrong variable names | Aaron LI | 2017-07-31 | 1 | -2/+2 |
| | | | | | | Also update logging message for RadioHalo a bit. Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | Add `fg21sim.share` module to hold globally shared "CONFIGS" & "COSMO" | Aaron LI | 2017-07-31 | 1 | -2/+1 |
| | | | | Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | clusters/halo.py: Add description for "time_crossing" property | Aaron LI | 2017-07-31 | 1 | -0/+6 |
| | | | | Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | cluster/halo.py: Add method "calc_brightness_mean()" | Aaron LI | 2017-07-30 | 1 | -0/+48 |
| | | | | Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | clusters/halo.py: Add property "angular_radius" (unit: [arcsec]) | Aaron LI | 2017-07-30 | 1 | -0/+10 |
| | | | | Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | clusters/halo.py: Add methods "calc_power()" and "calc_flux()" | Aaron LI | 2017-07-29 | 1 | -0/+50 |
| | | | | | | | Calculate the synchrotron power (a.k.a. specific luminosity) and flux density from the emissivity for the halo. Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | clusters/halo.py: Add method "calc_emissivity()" | Aaron LI | 2017-07-29 | 1 | -6/+37 |
| | | | | | | | This method invokes ``SynchrotronEmission`` to calculate the synchrotron emissivity at specified frequencies for the derived electron spectrum. Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | clusters/halo.py: Add property "volume" (unit: [cm^3]) | Aaron LI | 2017-07-29 | 1 | -0/+9 |
| | | | | Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | clusters/halo.py: Add property "magnetic_field" | Aaron LI | 2017-07-26 | 1 | -0/+11 |
| | | | | Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | clusters/halo: Change config "eta_turb" to "beta_turb" | Aaron LI | 2017-07-25 | 1 | -3/+3 |
| | | | | | | | | | | This configuration option "beta_turb" has different meaning as "eta_turb" in [Cassano & Brunetti 2005], though they can both be an indicator for the turbulence acceleration efficiency. Rename this configuration option to make the difference more clear. Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | clusters: Accept 1D numpy array and calculate values for all gamma's | Aaron LI | 2017-07-23 | 1 | -17/+16 |
| | | | | Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | clusters/halo.py: Constrain tau_acc to avoid zero diffusion | Aaron LI | 2017-07-23 | 1 | -2/+15 |
| | | | | | | | | | | Zero or negative diffusion coefficient leads to unstable or wrong results due to numerical algorithm/scheme adopted to solve the Fokker-Planck equation. Also add a NOTE to the FokkerPlanckSolver class. Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | Use all uppercase to identify the global CONFIGS and COSMO | Aaron LI | 2017-07-22 | 1 | -15/+11 |
| | | | | | | Update to use the global CONFIGS and COSMO Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | clusters: Import global "configs" to simplify parameters | Aaron LI | 2017-07-22 | 1 | -52/+54 |
| | | | | Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | clusters/halo.py: Implement "_injection_rate" property | Aaron LI | 2017-07-22 | 1 | -6/+56 |
| | | | | Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | clusters/halo.py: update parameters and methods | Aaron LI | 2017-07-22 | 1 | -46/+119 |
| | | | | | | Still WIP... Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | clusters/halo.py: More cleanups with minor updates | Aaron LI | 2017-07-22 | 1 | -70/+97 |
| | | | | Signed-off-by: Aaron LI <aly@aaronly.me> |