def test_cross_vs_ji_held_K2(self): """confirm cross agrees with ji-held model when K=2""" kappa_i_hat_brag = _nondim_tc_i_braginskii(self.big_hall, "cross") kappa_i_hat_jh = _nondim_tc_i_ji_held(self.big_hall, 1, 0, 100, "cross", K=2) assert np.isclose(kappa_i_hat_brag, kappa_i_hat_jh, rtol=2e-2)
def test_known_values_perp(self): """check some known values""" kappa_i_hat = _nondim_tc_i_braginskii(self.big_hall, field_orientation="perp") expected = 2.0 # Braginskii '65 eq (2.16) assert np.isclose(kappa_i_hat * self.big_hall**2, expected, atol=1e-1)
def test_unmagnetized(self): """confirm perp -> par as B -> 0""" kappa_i_hat_par = _nondim_tc_i_braginskii(self.small_hall, "par") kappa_i_hat_perp = _nondim_tc_i_braginskii(self.small_hall, "perp") assert np.isclose(kappa_i_hat_par, kappa_i_hat_perp, rtol=1e-3)
def test_known_values_par(self): """check some known values""" kappa_i_hat = _nondim_tc_i_braginskii(self.big_hall, field_orientation="par") expected = 3.9 # Braginskii '65 eq (2.15) assert np.isclose(kappa_i_hat, expected, atol=1e-1)