SHARE
TWEET

mat model ver 1.1

allekco Oct 29th, 2019 82 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1.  
  2. t = 2700.0 + 273.0
  3. T = 2973.0
  4. p = 9678.1
  5. Cp = 619.0
  6. small_lambda = 3.64 #теплопроводность
  7. big_lambda = 10196.0
  8. m = 0.61 * 0.000001
  9. beta = 8.555 * 0.00001 #температурный коэффициент линейного расширения
  10. energy = [0.01, 0.1, 0, 10, 100, 1000, 10000, 100000, 1000000, 10000000]
  11. sigma_t_o = [(4, 0.01), (4, 0.1), (4, 0), (4, 10), (5, 100), (1, 1000),
  12.              (6, 10000), (6, 100000), (6, 1000000), (6, 10000000)]
  13. 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),
  14.              (10**(-5), 10000), (3*10**(-5), 100000), (10**(-4), 1000000), (10**(-4), 10000000)]
  15. 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),
  16.                (2*10**(-4), 10000), (5*10**(-5), 100000), (3*10**(-2), 1000000), (2, 10000000)]
  17. sigma_c_238 = [(4, 0.01), (1, 0.1), (5*10**(-1), 0), (9*10**(-1), 10), (3*10, 100), (4, 1000),
  18.                (9*10**(-1), 10000), (2*10**(-1), 100000), (1*10**(-1), 1000000), (2*10**(-4), 10000000)]
  19. sigma_f_235 = [(10**3, 0.01), (300, 0.1), (80, 0), (70, 10), (20, 100), (9, 1000),
  20.                (3, 10000), (2, 100000), (1, 1000000), (2, 10000000)]
  21. sigma_c_235 = [(200, 0.01), (40, 0.1), (10, 0), (70, 10), (10, 100), (4, 1000),
  22.                (1, 10000), (0.6, 100000), (0.1, 1000000), (10**(-3), 10000000)]
  23. Na = 6.022*10**23
  24. N_238_nat = Na*19.05*10**6/(0.007*235 + 0.993*238) #metr^(-3)
  25. N_235_nat = Na*2.1*10**14/(0.007*235 + 0.993*238)
  26. N_o = Na*p/270.03
  27. Rf_nat = []
  28. Rc_nat = []
  29. Rf_o = []
  30. Rc_o = []
  31. for i in range(10):
  32.     if energy[i] < 1:
  33.         fi = 5.0*10**13/0.0001
  34.     else:
  35.         fi = 2.5*10**13/0.0001
  36.     Rf_nat.append(fi * (N_235_nat * sigma_f_235[i][0] + N_238_nat * sigma_f_238[i][0]))
  37.  
  38. for i in range(10):
  39.     if energy[i] < 1:
  40.         fi = 5.0*10**13/0.0001
  41.     else:
  42.         fi = 2.5*10**13/0.0001
  43.     Rc_nat.append(fi * (N_235_nat * sigma_c_235[i][0] + N_238_nat * sigma_c_238[i][0]))
  44.  
  45. for i in range(10):
  46.     if energy[i] < 1:
  47.         fi = 5.0*10**13/0.0001
  48.     else:
  49.         fi = 2.5*10**13/0.0001
  50.     Rf_o.append(fi * (32.4 * N_o * sigma_f_238[i][0] + 0.9 * N_o * sigma_f_235[i][0]))
  51.  
  52. for i in range(10):
  53.     if energy[i] < 1:
  54.         fi = 5.0*10**13/0.0001
  55.     else:
  56.         fi = 2.5*10**13/0.0001
  57.     Rf_o.append(fi * (32.4 * N_o * sigma_f_238[i][0] + 0.9 * N_o * sigma_f_235[i][0])
  58.                 + 33.3 * N_o * sigma_c_o[i][0])
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
 
Top