Example #1
0
 def P_loss(self, T_e, T_g, n_e, n_g):
     a = self.E_iz * n_e * n_g * self.K_iz(T_e)
     b = self.E_ex * n_e * n_g * self.K_ex(T_e)
     c = 3 * (m_e / self.m_i) * k * (T_e - T_g) * n_e * n_g * self.K_el(T_e)
     d = 7 * k * T_e * n_e * u_B(T_e, self.m_i) * A_eff(n_g, self.R, self.L) / self.V
     
     return a + b + c + d
Example #2
0
 def gas_heating(self, T_e, T_g, n_e, n_g):
     K_in = SIGMA_I * maxwellian_flux_speed(T_g, self.m_i)
     lambda_0 = self.R / 2.405 + self.L / pi
     # lambda_0 =np.sqrt((self.R / 2.405)**2 + (self.L / pi)**2)
     a = 3 * (m_e / self.m_i) * k * (T_e - T_g) * n_e * n_g * self.K_el(T_e)
     b = (1/4) * self.m_i * (u_B(T_e, self.m_i)**2) * n_e * n_g * K_in 
     c = self.kappa * (T_g - self.T_g_0) * self.A / (self.V * lambda_0)
     return a + b - c
Example #3
0
 def particle_balance_g(self, T_e, T_g, n_e, n_g):
     a = self.Q_g /self.V
     b = n_e * u_B(T_e, self.m_i) * A_eff_1(n_g, self.R, self.L, self.beta_i) / self.V
     c = n_e * n_g * self.K_iz(T_e)
     d = (1/4) * n_g * maxwellian_flux_speed(T_g, self.m_i) * self.A_g / self.V
     return a + b - c - d
Example #4
0
 def particle_balance_e(self, T_e, T_g, n_e, n_g):
     a = n_e * n_g * self.K_iz(T_e)
     b = n_e * u_B(T_e, self.m_i) * A_eff(n_g, self.R, self.L) / self.V
     return a - b
Example #5
0
 def flux_i(self, T_e, T_g, n_e, n_g):
     return h_L(n_g, self.L) * n_e * u_B(T_e, self.m_i)