mat model ver 2.0

import math
import matplotlib.pyplot as plt


t = 2700.0 + 273.0
T = 2973.0
p = 9678.1
Cp = 619.0
small_lambda = 3.64 #теплопроводность
big_lambda = 10196.0
m = 0.61 * 0.000001
beta = 8.555 * 0.00001 #температурный коэффициент линейного расширения
energy = [0.01, 0.1, 0, 10, 100, 1000, 10000, 100000, 1000000, 10000000]
sigma_t_o = [(4, 0.01), (4, 0.1), (4, 0), (4, 10), (5, 100), (1, 1000),
             (6, 10000), (6, 100000), (6, 1000000), (6, 10000000)]
sigma_c_o = [(3*10**(-4), 0.01), (10**(-4), 0.1), (3*10**(-5), 0), (10**(-5), 10), (5*10**(-6), 100), (6*10**(-6), 1000),
             (10**(-5), 10000), (3*10**(-5), 100000), (10**(-4), 1000000), (10**(-4), 10000000)]
sigma_f_238 = [(3*10**(-5), 0.01), (10**(-5), 0.1), (3*10**(-6), 0), (7*10**(-7), 10), (10**(-6), 100), (3*10**(-4), 1000),
               (2*10**(-4), 10000), (5*10**(-5), 100000), (3*10**(-2), 1000000), (2, 10000000)]
sigma_c_238 = [(4, 0.01), (1, 0.1), (5*10**(-1), 0), (9*10**(-1), 10), (3*10, 100), (4, 1000),
               (9*10**(-1), 10000), (2*10**(-1), 100000), (1*10**(-1), 1000000), (2*10**(-4), 10000000)]
sigma_f_235 = [(10**3, 0.01), (300, 0.1), (80, 0), (70, 10), (20, 100), (9, 1000),
               (3, 10000), (2, 100000), (1, 1000000), (2, 10000000)]
sigma_c_235 = [(200, 0.01), (40, 0.1), (10, 0), (70, 10), (10, 100), (4, 1000),
               (1, 10000), (0.6, 100000), (0.1, 1000000), (10**(-3), 10000000)]
Na = 6.022*10**23
N_238_nat = Na*19.05*10**6/(0.007*235 + 0.993*238) #metr^(-3)
N_235_nat = Na*2.1*10**14/(0.007*235 + 0.993*238)
N_o = Na*p/270.03
Rf_nat = []
Rc_nat = []
Rf_o = []
Rc_o = []
for i in range(10):
    if energy[i] < 1:
        fi = 5.0*10**13/0.0001
    else:
        fi = 2.5*10**13/0.0001
    Rf_nat.append(fi * (N_235_nat * sigma_f_235[i][0] + N_238_nat * sigma_f_238[i][0]))

for i in range(10):
    if energy[i] < 1:
        fi = 5.0*10**13/0.0001
    else:
        fi = 2.5*10**13/0.0001
    Rc_nat.append(fi * (N_235_nat * sigma_c_235[i][0] + N_238_nat * sigma_c_238[i][0]))

for i in range(10):
    if energy[i] < 1:
        fi = 5.0*10**13/0.0001
    else:
        fi = 2.5*10**13/0.0001
    Rf_o.append(fi * (32.4 * N_o * sigma_f_238[i][0] + 0.9 * N_o * sigma_f_235[i][0]))

for i in range(10):
    if energy[i] < 1:
        fi = 5.0*10**13/0.0001
    else:
        fi = 2.5*10**13/0.0001
    Rc_o.append(fi * (32.4 * N_o * sigma_c_238[i][0] + 0.9 * N_o * sigma_c_235[i][0])
                + 33.3 * N_o * sigma_c_o[i][0])
#for i in range(10):
 #   Rf_o[i] = math.log10(Rf_o[i])
  #  energy[i] = math.log10(energy[i])
plt.plot(energy, Rf_nat)
plt.title('R fission for nature uran')
plt.yscale('log')
plt.xscale('log')
plt.show()

plt.plot(energy, Rf_o)
plt.title('R fission for UO2')
plt.yscale('log')
plt.xscale('log')
plt.show()

plt.plot(energy, Rc_nat)
plt.title('R capture for nature uran')
plt.yscale('log')
plt.xscale('log')
plt.show()

plt.plot(energy, Rc_o)
plt.title('R capture for nature UO2')
plt.yscale('log')
plt.xscale('log')
plt.show()