Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- 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
- Rf_o.append(fi * (32.4 * N_o * sigma_f_238[i][0] + 0.9 * N_o * sigma_f_235[i][0])
- + 33.3 * N_o * sigma_c_o[i][0])
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement