def rho(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 85
'''
T = KtoC(T)
return 8940 * (1 - 3e-6 * (T - 20) *
(-2.49e-9 * T**3 + 8.18e-6 * T**2 + 3.16e-3 * T + 16.8))
def Syavg(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 44
'''
T = KtoC(T)
return 1.28 * (225.75 - 0.73683 * T + 2.5036e-3 * T**2 -
5.4546e-6 * T**3 + 6.4366e-9 * T**4 - 3.029e-12 * T**5)
def rho(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 97
'''
T = KtoC(T)
return 8900 * (1 - 3e-6 * (T - 20) *
(7.2e-9 * T**3 - 9.05e-6 * T**2 + 6.24e-3 * T + 16.6))
def Suavg(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 47
'''
T = KtoC(T)
return (529.75 - 0.95180 * T + 2.5732e-3 * T**2 - 2.3316e-6 * T**3 -
4.3944e-10 * T**4 + 3.4942e-13 * T**5) * 1.112
def CTE(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 40
T_ref = 293 K
'''
T = KtoC(T)
return 15.13 + 7.93e-3 * T - 3.33e-6 * T**2
def k(T: 'Kelvin'):
'''
DEMO_D_2MKKGB v0 Table 20.1
'''
T = KtoC(T)
return T * (0.19706 - 4.3053e-4 * T + 3.817e-7 * T**2 -
1.158e-0 * T**3)
def Suavg(T: 'Kelvin'):
'''
DEMO_D_2MKKGB v0 Table
'''
T = KtoC(T)
return (670.1 - 0.904 * T + 0.00401 * T**2 - 1.091e-5 * T**3 +
1.115e-8 * T**4 - 4.75e-12 * T**5)
def Sy(T: 'Kelvin'):
'''
DEMO_D_2MKKGB v0 Table 6.4
'''
T = KtoC(T)
return (491.5 - 0.627 * T + 0.00464 * T**2 - 1.744e-5 * T**3 +
2.68e-8 * T**4 - 1.59e-11 * T**5)
def rho(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 129
'''
T = KtoC(T)
return 1830 * (1 - 3e-6 *
(4.1e-9 * T**3 - 1.2e-5 * T**2 + 1.33E-2 * T + 11.3) *
(T - 20))
def Cp(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Table A.A01.2.6-1
'''
T = KtoC(T)
t = [20, 100, 200, 300, 400, 500, 600, 700, 800]
cp = [1807, 2046, 2295, 2497, 2660, 2791, 2998, 2989, 3071]
return interp1d(t, cp)(T)
def CTE(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Table A.S31.2.1-1
'''
T = KtoC(T)
t = [20, 50, 100, 150, 200, 250, 300, 400, 450, 500, 550, 600]
a = [
16.7, 17, 17.3, 17.5, 17.7, 17.8, 18, 18.1, 18.2, 18.4, 18.5, 18.6
]
return interp1d(CtoK(t), a)(T)
def E(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 128
'''
T = KtoC(T)
return 306.78 - 2.9281e-2 * T - 4.5069e-5 * T**2
def CTE(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 83
'''
T = KtoC(T)
return 16.73 + 5.26e-3 * T - 5.53e-6 * T**2 + 4.2e-9 * T**3
def Cp(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 139
'''
T = KtoC(T)
return 128.308 + 3.2797e-2 * T - 3.4097e-6 * T**2
def CTE(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 127
'''
T = KtoC(T)
return 4.1e-9 * T**3 - 1.2e-5 * T**2 + 1.33e-2 * T + 11.3
def Sy(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 91
'''
T = KtoC(T)
return 584 - 5.04e-2 * T - 5e-6 * T**2
def k(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 138
'''
T = KtoC(T)
return 174.9274 - 0.1067 * T + 5.0067e-5 * T**2 - 7.8349e-9 * T**3
def CTE(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 135
'''
T = KtoC(T)
return 3.9225 + 5.8352e-4 * T + 5.7054e-11 * T**2 - 2.0463e-14 * T**3
def Sy(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 131
'''
T = KtoC(T)
return 252.2 - 1.524e-1 * T - 1.45E-4 * T**2
def E(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 84
'''
T = KtoC(T)
return 117 - 2.17e-2 * T - 6.29e-5 * T**2
def k(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 98
'''
T = KtoC(T)
return 2.11e-7 * T**3 - 2.83e-4 * T**2 - 1.38e-1 * T + 323
def k(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 86
'''
T = KtoC(T)
return 402 + 2.46e-6 * T**2 - 7.08e-2 * T
def E(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 96
'''
T = KtoC(T)
return 128 - 2.59e-2 * T - 4.87e-5 * T**2
def Su(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 92
'''
T = KtoC(T)
return 212 - 2.94e-1 * T + 8.17e-5 * T**2
def k(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 130
'''
T = KtoC(T)
return 209 - 4.8e-1 * T + 9.38E-4 * T**2 - 9.37e-7 * T**3 + 3.34e-10 * T**4
def E(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 136
'''
T = KtoC(T)
return 397.903 - 2.3066e-3 * T - 2.7162e-5 * T**2
def Cp(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 43
'''
T = KtoC(T)
return 462.69 + 0.520265 * T - 1.7117e-3 * T**2 + 3.3658e-6 * T**3 - 2.1958e-9 * T**4
def rho(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 137
'''
T = KtoC(T)
return 1000 * (19.3027 - 2.3786e-4 * T - 2.2448e-8 * T**2)
def Su(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 132
'''
T = KtoC(T)
return 392.8 - 3.023e-1 * T - 2.011e-4 * T**2
def Cp(T: 'Kelvin'):
'''
ITER_D_222RLN v3.3 Equation 87
'''
T = KtoC(T)
return 387 + 4.36e-5 * T**2 + 6.35e-2 * T