Changeset 302:a6d1f859b52e for source/class

Timestamp:

Sep 19, 2011, 7:20:22 PM (13 years ago)

Author:

Evgeny Stambulchik <Evgeny.Stambulchik@…>

Branch:

default

Phase:

public

Message:

Simplified expressions for radiative losses.

File:

• source/class/pf/Plasma.js (modified) (4 diffs)

source/class/pf/Plasma.js

@@ -167 +167 @@
             var v_e = this.e.getThermalVelocity();
             
-            return 8*Math.PI/(3*Math.pow(pf.base.Bohr.c, 3))*v_e*this.e.getN()*
+            return 8*Math.PI*pf.base.Bohr.alpha3/3*v_e*this.e.getN()*
                 (this.i.getN()*z_i*z_i + this.r.getN()*z_r*z_r);
         },
@@ -180 +180 @@
         getFreeBoundLosses: function()
         {
-            var z_i = this.i.getQ();
-            var m_i = this.i.getM();
-            var z_r = this.r.getQ();
-            var m_r = this.r.getM();
-            var v_e = this.e.getThermalVelocity();
-
-            return -8*Math.PI/(3*Math.pow(pf.base.Bohr.c, 3))*v_e*this.e.getN()*
-                (this.i.getN()*z_i*z_i*pf.base.Bohr.bindingEnergy(z_i, m_i, 1) +
-                 this.r.getN()*z_r*z_r*pf.base.Bohr.bindingEnergy(z_r, m_r, 1))/
-                 this.e.getT();
-        },
-        
-        getFreeBoundSpectralDensity: function()
-        {
-            var photon_en = this.r.getTransitionEnergy();
             var T_e = this.e.getT();
             var v_e = this.e.getThermalVelocity();
             var N_e = this.e.getN();
-            
+
             var sum = 0;
             
@@ -214 +199 @@
                 var Eb = -pf.base.Bohr.bindingEnergy(z_s, m_s, 1);
                 
+                sum += N_s*z_s*z_s*(Eb/T_e);
+            }
+            
+            return 8*Math.PI*pf.base.Bohr.alpha3/3*N_e*v_e*sum;
+        },
+        
+        getFreeBoundSpectralDensity: function()
+        {
+            var photon_en = this.r.getTransitionEnergy();
+            var T_e = this.e.getT();
+            var v_e = this.e.getThermalVelocity();
+            var N_e = this.e.getN();
+            
+            var sum = 0;
+            
+            var sa = [this.i, this.r];
+            for (var i =0; i < sa.length; i++) {
+                var s = sa[i];
+                var N_s = s.getN();
+                if (!N_s) {
+                    continue;
+                }
+                
+                var z_s = s.getQ();
+                var m_s = s.getM();
+                
+                var Eb = -pf.base.Bohr.bindingEnergy(z_s, m_s, 1);
+                
                 if (photon_en < Eb) {
                     continue;
@@ -221 +234 @@
             }
 
-            return 8*Math.PI/(3*Math.pow(pf.base.Bohr.c, 3))*v_e*N_e/T_e*sum;
+            return 8*Math.PI*pf.base.Bohr.alpha3/3*N_e/v_e*sum;
         },