Changeset 100:e7d27289faaf

Timestamp:

Jan 11, 2009, 12:03:25 PM (16 years ago)

Author:

fnevgeny

Branch:

default

Phase:

public

Message:

Added Coulomb logarithm.

File:

• source/class/pf/base/Species.js (modified) (1 diff)

source/class/pf/base/Species.js

@@ -180 +180 @@
         {
             return 1.602164e-12*this._n*this._t;
+        },
+        
+        getCoulombLog: function(s2)
+        {
+            var lambda;
+            if (this.id == "e" && s2.id == "e") {
+                lambda = 24 - Math.log(Math.sqrt(this._n)/this._t)
+            } else
+            if (this.id != "e" && s2.id != "e") {
+                lambda = 23 - Math.log(
+                    this._q*s2._q*(this._m + s2._m)/(this._m*s2._t + s2._m*this._t)*
+                    Math.sqrt(this._n*this._q*this._q/this._t + s2._n*s2._q*s2._q/s2._t));
+            } else {
+                var se, si;
+                if (this.id == "e") {
+                    se = this;
+                    si = s2;
+                } else {
+                    si = this;
+                    se = s2;
+                }
+                
+                var tmm = si._t*se._m/si._m;
+                
+                if (se._t < si._q*tmm) {
+                    lambda = 30 - Math.log(Math.sqrt(si._n)/Math.pow(si._t, 3/2)*
+                        si._q*si._q/si._m);
+                } else
+                if (tmm   < 10*si._q*si._q && se._t > 10*si._q*si._q) {
+                    lambda = 24 - Math.log(Math.sqrt(se._n)/se._t);
+                } else
+                if (tmm < se._t && se._t < 10*si._q*si._q) {
+                    lambda = 23 -
+                        Math.log(Math.sqrt(se._n)*si._q/Math.pow(se._t, 3/2));
+                } else {
+                    // no way???
+                    lambda = 0;
+                }
+            }
+            
+            return lambda;
         }
     },