blob: 6ed557a6115c1338d907fce1038d54010da0945e (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
// Scilab code Exa4.8.3: To calculate the energy released and temperature required for fusion of given gases : Page 194 (2011)
// Firstly calculate for B-10
Z_B = 5; // Atomic number of B-10
r_B = 5.17; // Seperation of two nuclei, fm
K = 1.38e-023; // Boltzmann's constant
F = 1/137; // Fine structure constant
E = 197.5*1.6e-013; // Energy, J
V_c_B = F*Z_B^2*E/r_B; // Coulomb barrier for B-10, J
T_B = 2/3*V_c_B/K; // Temperature required to overcome the barrier for B-10, K
// Now calculate for Mg-24
Z_Mg = 12; // Atomic number of Mg-24
r_Mg = 6.92; // Seperation of two nuclei, fm
K = 1.38e-023; // Boltzmann's constant
F = 1/137; // Fine structure constant
E = 197.5*1.6e-013; // Energy, J
V_c_Mg = F*Z_Mg^2*E/r_Mg; // Coulomb barrier for Mg-24, J
T_Mg = 2/3*V_c_Mg/K; // Temperature required to overcome the barrier for Mg-24, K
printf("\nFor B-10 \n Energy released = %4.2e J \n Temperature required = %4.1e K \nFor Mg-24 \n Energy released = %4.2e J \n Temperature required = %4.2e K", V_c_B,T_B, V_c_Mg, T_Mg)
// Result
// For B-10
// Energy released = 1.12e-012 J
// Temperature required = 5.4e+010 K
// For Mg-24
// Energy released = 4.80e-012 J
// Temperature required = 2.32e+011 K
|
