diff options
Diffstat (limited to '1970/CH15/EX15.1')
-rwxr-xr-x | 1970/CH15/EX15.1/CH15Exa1.sce | 20 |
1 files changed, 20 insertions, 0 deletions
diff --git a/1970/CH15/EX15.1/CH15Exa1.sce b/1970/CH15/EX15.1/CH15Exa1.sce new file mode 100755 index 000000000..0a58da351 --- /dev/null +++ b/1970/CH15/EX15.1/CH15Exa1.sce @@ -0,0 +1,20 @@ +// Scilab code Exa15.1 : : Page-652 (2011) +clc; clear; +N_0_235 = 1; // Number of uranium atom +N_0_c = 10^5; // Number of graphite atoms per uranium atom +sigma_a_235 = 698; // Absorption cross section for uranium, barns +sigma_a_c = 0.003; // Absorption cross section for graphite, barns +f = N_0_235*sigma_a_235/(N_0_235*sigma_a_235+N_0_c*sigma_a_c ); // Thermal utilization factor +eta = 2.08; // Number of fast fission neutron produced +k_inf = eta*f; // Multiplication factor +L_m = 0.54; // Material length, metre +L_sqr = ((L_m)^2*(1-f)); // diffusion length, metre +tau = 0.0364; // Age of the neutron +B_sqr = 3.27; // Geometrical buckling +k_eff = round (k_inf*exp(-tau*B_sqr)/(1+L_sqr*B_sqr)); // Effective multiplication factor +N_lf = k_eff/k_inf; // Non leakage factor +lf = (1-N_lf)*100; // Leakage factor, percent +printf("\n Total leakage factor = %4.1f percent",lf) + +// Result +// Total leakage factor = 31.3 percent
\ No newline at end of file |