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)
