source: source/class/pf/Plasma.js @ 364:5443fc795688

qx.Class.define("pf.Plasma",
{
    extend: qx.core.Object,
    
    members:
    {
        R0: 0.5, // 0.3/0.7/1 ?
        
        e: null,
        i: null,
        r: null,
        P_r: 0,
        B  : 0,
        
        getMicrofieldFrequency: function(s)
        {
            var r = s.getTypicalDistance();
            var v_r = this.r.getThermalVelocity();
            var v_s = s.getThermalVelocity();
            var v = Math.sqrt(v_r*v_r + v_s*v_s);
            
            return v/r;
        },
        
        getFullDebyeLength: function(s)
        {
            var tmp = 0;
            var m = s.getM();
            
            var sa = new Array(this.e, this.i, this.r);
            for (var i = 0; i < 3; i++) {
                var os = sa[i];
                if (os.getM() <= m) {
                    var ld = os.getDebyeLength();
                    if (ld) {
                        tmp += 1/(ld*ld);
                    }
                }
                
            }
            
            return 1/Math.sqrt(tmp);
        },
        
        getPpiFactor: function(s)
        {
            var dl = this.getFullDebyeLength(s);
            var r1 = s.getTypicalDistance();
            var r;
            if (isFinite(dl) && isFinite(r1)) {
                r = r1/dl;
            } else {
                r = 0;
            }
            
            return (1 + r)*Math.exp(-r);
            // return 1 - 0.6*r
        },
        
        getRpiFactor: function(s)
        {
            var r1 = s.getTypicalDistance();
            var q_r = this.r.getQ();
            var q_s = s.getQ();
            if (q_r == 0 || q_s == 0) {
                return 1;
            }
             
            var r_m = s.getMinApproachDistance()*q_r/q_s;
            return Math.exp(-r_m/r1);
        },
        
        getDelta0: function(s)
        {
            var n_u = this.r.getN_u();
            var n_l = this.r.getN_l();
            var Z_core = this.r.getZcore();
            var ef = s.getHoltsmarkField();
            
            return 3/2*(n_u*n_u - n_l*n_l)/Z_core*ef;
        },
        
        getQsHwhm: function(s)
        {
            var n_u = this.r.getN_u();
            var n_l = this.r.getN_l();

            var ppi = this.getPpiFactor(s);
            var rpi = this.getRpiFactor(s);
            
            // 1.4385 is S_1 HWHM
            var hwhm = 1.4385*this.getDelta0(s)*ppi*rpi;
            
            // rude correction for alpha lines
            if (n_u - n_l == 1) {
                hwhm /= 2;
            }
            
            return hwhm;
        },
        
        getDynHwhm: function(s)
        {
            // neglect dynamic PP effects
            var ppi = 1.0;
            // assume dynamic RP effect == static one
            var rpi = this.getRpiFactor(s);

            var fwhm = this.getMicrofieldFrequency(s)*ppi*rpi;
            
            return fwhm/2;
        },
        
        getStarkRatio: function(s)
        {
            var R;
            var w_qs = this.getQsHwhm(s);
            if (w_qs == 0) {
                R = 0;
            } else {
                var w_dyn = this.getDynHwhm(s);
                R = w_qs/w_dyn;
            }
            
            return R;
        },
        
        getStarkQuasistaticity: function(s)
        {
            var R = this.getStarkRatio(s);
            return R/(R + this.R0);
        },
        
        getStarkHwhm: function(s)
        {
            var w_qs = this.getQsHwhm(s);
            var f  = this.getStarkQuasistaticity(s);
            
            return f*w_qs;
        },
        
        getTotalStarkHwhm: function()
        {
            var w = 0;
            
            var sa = new Array(this.e, this.i, this.r);
            for (var i = 0; i < 3; i++) {
                var s = sa[i];
                var w_s = this.getStarkHwhm(s);

                w += Math.pow(w_s, 3/2);
            }
            
            return Math.pow(w, 2/3);
        },
        
        getZeemanSplitting: function()
        {
            return pf.base.Bohr.mu_B*this.B;
        },

        getTotalFwhm: function()
        {
            var s_hwhm = this.getTotalStarkHwhm();
            var z_hwhm = this.getZeemanSplitting();
            var d_hwhm = this.r.getDopplerHwhm();
            return 2*Math.sqrt(s_hwhm*s_hwhm + z_hwhm*z_hwhm + d_hwhm*d_hwhm);
        },

        getBremsstrahlungLosses: function()
        {
            var z_i = this.i.getQ();
            var z_r = this.r.getQ();
            var v_e = this.e.getThermalVelocity();
            
            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);
        },
        
        getBremsstrahlungSpectralDensity: function()
        {
            var photon_en = this.r.getTransitionEnergy();
            return this.getBremsstrahlungLosses()/this.e.getT()*
                Math.exp(-photon_en/this.e.getT());
        },
        
        getFreeBoundLosses: function()
        {
            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);
                
                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;
                }
                
                sum += N_s*z_s*z_s*(Eb/T_e)*Math.exp(-(photon_en - Eb)/T_e);
            }

            return 8*Math.PI*pf.base.Bohr.alpha3/3*N_e/v_e*sum;
        },
        
        getPressure: function()
        {
            return this.e.getPressure() +
                   this.i.getPressure() +
                   this.r.getPressure();
        },
        
        getMagneticFieldPressure: function()
        {
            return this.B*this.B/(8*Math.PI);
        },
        
        getBeta: function()
        {
            return this.getPressure()/this.getMagneticFieldPressure();
        },
        
        getFermiEnergy: function()
        {
            return pf.base.Thermo.fermiEnergy(this.e.getN());
        },
        
        getFermiVelocity: function()
        {
            return Math.sqrt(2*this.getFermiEnergy());
        },
        
        getFermiLength: function()
        {
            var n_e = this.e.getN();
            var E_F = this.getFermiEnergy();
            return Math.sqrt(E_F/(6*Math.PI*n_e));
        },
        
        getAlfvenVelocity: function()
        {
            var v_m1;
            
            v_m1 = 1/this.e.getAlfvenVelocity(this.B) +
                   1/this.i.getAlfvenVelocity(this.B);
        
            if (this.r.getQ() != 0) {
                v_m1 += 1/this.r.getAlfvenVelocity(this.B);
            }
            
            return 1/v_m1;
        },

        getCoulombLog: function(s1, s2)
        {
            var lambda;
            
            var m1 = s1.getM();
            var m2 = s2.getM();
            var v1 = s1.getThermalVelocity();
            var v2 = s2.getThermalVelocity();
            var q1 = s1.getQ();
            var q2 = s2.getQ();
            
            // kinetic energy in the CoM frame
            var Kp = 3*m1*m2/(2*(m1 + m2))*(v1 + v2)*(v1 + v2);
            
            var b_min = Math.abs(q1*q1/Kp);
            var b_max = this.getFullDebyeLength(s1) +
                        this.getFullDebyeLength(s2);
            
            lambda = Math.log(b_max/b_min);
            if (lambda < 1.0) {
                lambda = 1.0;
            }
            
            return lambda;
        },

        getRelaxationRate: function(s1, s2)
        {
            var m1 = s1.getM();
            var m2 = s2.getM();
            var v1 = s1.getThermalVelocity();
            var v2 = s2.getThermalVelocity();
            var q1 = s1.getQ();
            var q2 = s2.getQ();
            
            var n2 = s2.getN();
            
            var lambda = this.getCoulombLog(s1, s2);
            
            // LL X (42,5) with the missing factor of 3
            return Math.sqrt(128*Math.PI)/3*q1*q1*q2*q2*n2*lambda/
                   (Math.pow(v1*v1 + v2*v2, 3/2)*m1*m2);
        },

        getDreicerField: function()
        {
            var n_e = this.e.getN();
            var T_e = this.e.getT();
            var lambda = this.getCoulombLog(this.e, this.e);
            
            return 4*Math.PI*lambda*n_e/T_e;
        },
        
        getMassDensity: function()
        {
            return this.e.getMassDensity() +
                   this.i.getMassDensity() +
                   this.r.getMassDensity();
        },
        
        // griem:1997a, Eq (7.75)
        getLteNe: function()
        {
            var Z_core = this.r.getZcore();
            var E12 = 3.0/8.0*Z_core*Z_core;
            return 5/Math.sqrt(Math.PI)*pf.base.Bohr.alpha3*
                Math.pow(E12, 3)*Math.sqrt(2*this.e.getT());
        },
        
        dE2dWl: function(dE)
        {
            var E = this.r.getTransitionEnergy();
            return (2*Math.PI)*pf.base.Bohr.c*dE/(E*E);
        },

        checkParameters: function()
        {
            if ((this.P_r < 0.0 || this.P_r > 1.0) ||
                this.i.getQ() < 1 ||
                this.P_r == 1.0 && this.r.getQ() == 0) {
                return false;
            } else {
                return true;
            }
        },

        updateParameters: function()
        {
            // Safety measures
            if (this.checkParameters() != true) {
                return false;
            }
            
            var N_e = this.e.getN();
            var N_i, N_r;
            
            if (this.P_r < 1.0) {
                N_i = N_e*(1.0 - this.P_r)/
                    (this.i.getQ()*(1.0 - this.P_r) + this.r.getQ()*this.P_r);
                
                N_r = N_i*this.P_r/(1.0 - this.P_r);
            } else {
                N_i = 0.0;
                N_r = N_e/this.r.getQ();
            }
            
            this.i.setN(N_i);
            this.r.setN(N_r);
            
            return true;
        }
    },
    
    construct: function()
    {
        var m_p = pf.base.Bohr.m_p;
        var V0  = pf.base.Bohr.V0;
        var E0  = pf.base.Bohr.E0;
        
        // Defaults: 1e16 1/cc, 1 eV
        this.e = new pf.base.Species ("e", 1.0, -1, 1e16*V0, 1.0/E0);
        this.i = new pf.base.Species ("i", m_p, +1, 1e16*V0, 1.0/E0);
        this.r = new pf.base.Radiator("r", m_p,  0,       0, 1.0/E0);

        this.P_r = 0;

        this.B   = 0;
    }
});