• API
• FAQ
• Tools
• Archive
SHARE
TWEET

# mat model ver 2.1

allekco Oct 29th, 2019 81 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
1. import math
2. import pylab
3. from matplotlib import mlab
4.
5.
6. t = 2700.0 + 273.0
7. T = 2973.0
8. p = 9678.1 #кг/м^3
9. Cp = 619.0
10. small_lambda = 3.64 #теплопроводность
11. big_lambda = 10196.0
12. m = 0.61 * 0.000001
13. beta = 8.555 * 0.00001 #температурный коэффициент линейного расширения
14. energy = [0.01, 0.1, 0, 10, 100, 1000, 10000, 100000, 1000000, 10000000]
15. sigma_t_o = [4, 4, 4, 4, 5, 1, 6, 6, 6, 6]
16. sigma_c_o = [3*10**(-4), 10**(-4), 3*10**(-5), 10**(-5), 5*10**(-6), 6*10**(-6),
17.              10**(-5), 3*10**(-5), 10**(-4), 10**(-4)]
18. sigma_f_238 = [3*10**(-5), 10**(-5), 3*10**(-6), 7*10**(-7), 10**(-6), 3*10**(-4),
19.                2*10**(-4), 5*10**(-5), 3*10**(-2), 2]
20. sigma_c_238 = [4, 1,  5*10**(-1), 9*10**(-1), 3*10, 4, 9*10**(-1), 2*10**(-1), 1*10**(-1), 2*10**(-4)]
21. sigma_f_235 = [10**3, 300, 80, 70, 20, 9, 3, 2, 1, 2]
22. sigma_c_235 = [200, 40, 10, 70, 10, 4, 1, 0.6, 0.1, 10**(-3)]
23. for i in range(10):
24.     sigma_t_o[i] = sigma_t_o[i] * 10**(-28)   #metres^2
25.     sigma_c_o[i] = sigma_c_o[i] * 10**(-28)
26.     sigma_c_235[i] = sigma_c_235[i] * 10**(-28)
27.     sigma_c_238[i] = sigma_c_238[i] * 10**(-28)
28.     sigma_f_235[i] = sigma_f_235[i] * 10**(-28)
29.     sigma_f_238[i] = sigma_f_238[i] * 10**(-28)
30.
31. Na = 6.022*10**23
32. Mr = ((0.007*235 + 0.993*238)*0.001)
33. Mr_UO2 = 0.27
34. N_238_nat = Na*19050/Mr
35. N_235_nat = Na*18700/Mr
36. N_o = Na*p/Mr_UO2
37. Rf_nat = []
38. Rc_nat = []
39. Rf_o = []
40. Rc_o = []
41. fi1 = 5.0*10**13/0.0001 # 1/(metr^2*sec)
42. fi2 = 2.5*10**13/0.0001
43. for i in range(10):
44.     if energy[i] < 1:
45.         fi = fi1
46.     else:
47.         fi = fi2
48.     Rf_nat.append(fi * (N_235_nat * sigma_f_235[i] + N_238_nat * sigma_f_238[i]))
49.
50. for i in range(10):
51.     if energy[i] < 1:
52.         fi = fi1
53.     else:
54.         fi = fi2
55.     Rc_nat.append(fi * (N_235_nat * sigma_c_235[i] + N_238_nat * sigma_c_238[i]))
56.
57. for i in range(10):
58.     if energy[i] < 1:
59.         fi = fi1
60.     else:
61.         fi = fi2
62.     Rf_o.append(fi * (32.4 * N_o * sigma_f_238[i] + 0.9 * N_o * sigma_f_235[i]))
63.
64. for i in range(10):
65.     if energy[i] < 1:
66.         fi = fi1
67.     else:
68.         fi = fi2
69.     Rc_o.append(fi * (32.4 * N_o * sigma_c_238[i]+ 0.9 * N_o * sigma_c_235[i])
70.                 + 33.3 * N_o * sigma_c_o[i])
71.
72. #Текущая ячейка - 1
73. #Две строки, три столбца.
74. pylab.subplot(2, 2, 1)
75. pylab.plot(energy, Rc_nat)
76. pylab.xscale('log')
77. pylab.yscale('log')
78. pylab.title("R capture for nature uran")
79.
80. pylab.subplot(2, 2, 2)
81. pylab.plot(energy, Rc_o)
82. pylab.xscale('log')
83. pylab.yscale('log')
84. pylab.title("R capture for nature UO2")
85.
86. pylab.subplot(2, 2, 3)
87. pylab.plot(energy, Rf_nat)
88. pylab.xscale('log')
89. pylab.yscale('log')
90. pylab.title("R fission for nature uran")
91.
92. print(*Rf_o)
93. pylab.subplot(2, 2, 4)
94. pylab.plot(energy, Rf_o)
95. pylab.xscale('log')
96. pylab.yscale('log')
97. pylab.title("R fission for UO2")
98.
99. pylab.show()
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.
Not a member of Pastebin yet?