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
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
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
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
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)