Changes between Version 3 and Version 4 of u/johannjc/scratchpad22
- Timestamp:
- 01/12/17 13:39:26 (8 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
u/johannjc/scratchpad22
v3 v4 1 = = Scaling ==1 = Scaling = 2 2 3 The equations in physical units look like 3 == Cross section == 4 The equations for calculating the optical depth in physical units look like 4 5 5 6 $d\tau = n \sigma dl$ 6 7 $\frac{dn}{dt} = n \sigma F$8 7 9 8 And in computational units … … 11 10 $d\tau = n_c \sigma_c dl_c$ 12 11 13 $\frac{dn_c}{dt_c} = n_c \sigma_c F_c$ 12 Dividing the first equation by the second gives 14 13 15 Dividing the first equation gives$n_c \sigma_c dl_c = n \sigma dl = n_c N_s \sigma dl_c L_s$14 $n_c \sigma_c dl_c = n \sigma dl = n_c N_s \sigma dl_c L_s$ 16 15 17 16 where $N_s$ and $L_s$ are the number density and length scales in the code. … … 21 20 $\sigma_c = \sigma N_s L_s$ 22 21 23 And dividing the second equation by the first, gives 22 == Ionizing Flux == 23 The equation for calculating the ionization rate in physical units is 24 25 $\frac{dn}{dt} = n \sigma F$ 26 27 and in computational units 28 29 $\frac{dn_c}{dt_c} = n_c \sigma_c F_c$ 30 31 32 Dividing the first equation by the second, gives 24 33 25 34 $N_s/T_s = N_s \frac{1}{N_s L_s} \frac{F}{F_c}$ … … 28 37 29 38 $F_c=F \frac{T_s}{L_s N_s}$ 39 40 == Recombination Rate == 41 42 The equation for recombination in physical units is 43 44 $\frac{d n}{dt} = \alpha n_{HII} n_e$ 45 46 and in computational units 47 48 $\frac{d n_c}{dt} = \alpha_c n_{c,HII} n_{c,e}$ 49 50 Dividing the two and solving for $\alpha_c$ gives 51 52 $\alpha_c = \alpha \frac{T_s}{N_s}$ 53 54 == Ionization Heating == 55 The equation in physical units is 56 57 $\frac{dE}{dt}=e_\Gamma \sigma n F$ 58 59 and in computational units 60 61 $\frac{dE_c}{dt_c}=e_{c,\Gamma} \sigma_c n_c F_c$ 62 63 Dividing the two and solving for $e_{c,\Gamma}$ gives 64 65 $e_{c,\Gamma} = \frac{\sigma}{\sigma_c} N_s \frac{F}{F_c}\frac{T_s}{E_s} = \frac{N_s}{E_s}$ 66 30 67 31 68