def test_cross_vs_ji_held(self, Z):
"""cross should roughly agree with ji-held"""
kappa_e_hat_cross_brag = _nondim_tc_e_braginskii(self.big_hall, Z, "cross")
kappa_e_hat_cross_jh = _nondim_tc_e_ji_held(self.big_hall, Z, "cross")
testTrue = np.isclose(kappa_e_hat_cross_brag, kappa_e_hat_cross_jh, rtol=2e-2)
errStr = (
f"braginskii cross {kappa_e_hat_cross_brag} should agree "
f"with ji-held {kappa_e_hat_cross_jh}."
)
assert testTrue, errStr
def test__nondim_tc_e_spitzer(Z):
"""test _nondim_tc_e_spitzer function"""
kappa = _nondim_tc_e_spitzer(Z)
if Z == 1:
kappa_check = 3.203
rtol = 1e-3
elif Z in [2, 4]:
kappa_check = _nondim_tc_e_braginskii(0, Z, "par")
rtol = 2e-2
elif Z == 16:
kappa_check = _nondim_tc_e_ji_held(0, Z, "par")
rtol = 2e-2
elif Z == np.inf:
kappa_check = _nondim_tc_e_ji_held(0, 1e6, "par")
rtol = 2e-2
assert np.isclose(kappa, kappa_check, rtol=rtol)
def test_unmagnetized(self, Z):
"""confirm perp -> par as B -> 0"""
kappa_e_hat_par = _nondim_tc_e_braginskii(self.small_hall, Z, "par")
kappa_e_hat_perp = _nondim_tc_e_braginskii(self.small_hall, Z, "perp")
assert np.isclose(kappa_e_hat_par, kappa_e_hat_perp, rtol=1e-3)
def test_known_values_perp(self, Z, field_orientation, expected):
"""check some known values"""
kappa_e_hat = _nondim_tc_e_braginskii(self.big_hall, Z,
field_orientation)
assert np.isclose(kappa_e_hat * self.big_hall**2, expected, atol=1e-1)
