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. OK, I Understand
Not a member of Pastebin yet?
Sign Up, it unlocks many cool features!
 
Top