Changeset 301:7bef42595c00 for source
- Timestamp:
- Sep 19, 2011, 6:58:03 PM (13 years ago)
- Branch:
- default
- Phase:
- public
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
source/class/pf/Plasma.js
r295 r301 181 181 { 182 182 var z_i = this.i.getQ(); 183 var m_i = this.i.getM(); 183 184 var z_r = this.r.getQ(); 185 var m_r = this.r.getM(); 184 186 var v_e = this.e.getThermalVelocity(); 185 var m_r = this.r.getM();186 187 187 188 return -8*Math.PI/(3*Math.pow(pf.base.Bohr.c, 3))*v_e*this.e.getN()* 188 (this.i.getN()*z_i*z_i*pf.base.Bohr.bindingEnergy(z_i, m_ r, 1) +189 (this.i.getN()*z_i*z_i*pf.base.Bohr.bindingEnergy(z_i, m_i, 1) + 189 190 this.r.getN()*z_r*z_r*pf.base.Bohr.bindingEnergy(z_r, m_r, 1))/ 190 191 this.e.getT(); … … 194 195 { 195 196 var photon_en = this.r.getTransitionEnergy(); 196 return this.getFreeBoundLosses()/this.e.getT()* 197 Math.exp(-photon_en/this.e.getT()); 197 var T_e = this.e.getT(); 198 var v_e = this.e.getThermalVelocity(); 199 var N_e = this.e.getN(); 200 201 var sum = 0; 202 203 var sa = [this.i, this.r]; 204 for (var i =0; i < sa.length; i++) { 205 var s = sa[i]; 206 var N_s = s.getN(); 207 if (!N_s) { 208 continue; 209 } 210 211 var z_s = s.getQ(); 212 var m_s = s.getM(); 213 214 var Eb = -pf.base.Bohr.bindingEnergy(z_s, m_s, 1); 215 216 if (photon_en < Eb) { 217 continue; 218 } 219 220 sum += N_s*z_s*z_s*(Eb/T_e)*Math.exp(-(photon_en - Eb)/T_e); 221 } 222 223 return 8*Math.PI/(3*Math.pow(pf.base.Bohr.c, 3))*v_e*N_e/T_e*sum; 198 224 }, 199 225
Note: See TracChangeset
for help on using the changeset viewer.