该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
%1.1.2.1快群
%堆芯慢化能力
%kesi_macr=zeros(5,1);
%for i=1:5
%kesi_macr=kesi(1,i)*N(1,i)*micro_s1(1,i)
%end
kesi_macr=kesi_H*N_H*micro_s1_H+kesi_O*N_O*micro_s1_O+kesi_Zr*N_Zr*micro_s1_Zr+kesi_5*N_5*micro_s1_5+kesi_8*N_8*micro_s1_8;
%快群宏观输运截面
macr_tr1=N_5*micro_s1_5*0.9972+N_8*micro_s1_8*0.9972+N_H*micro_s1_H*0.3386+N_O*micro_s1_O*0.9583+N_Zr*micro_s1_Zr*0.9927
%快群宏观扩散系数
D1=1/(3*macr_tr1)
%1.1.2.1热群
L_sigama_th=N_5*micro_a_5+N_8*micro_a_8+N_O*micro_a_O+N_H*micro_a_H+N_Zr*micro_a_Zr
L_sigama_kTM=L_sigama_th*sqrt(293/TM) %式(2-80)
Tn=TM*(1+1.4*L_sigama_kTM/kesi_macr) %式(2-79)
ga_5=1.03888*(Tn-273)^(-0.0184)
ga_8=0.9855*(Tn-273)^(0.0043)
gf_5=1.053*(Tn-273)^(-0.0226)
micro_a2_5=micro_a_5/1.128*sqrt(293/Tn)*ga_5;
micro_a2_8=micro_a_8/1.128*sqrt(293/Tn)*ga_8;
micro_a2_Zr=micro_a_Zr/1.128*sqrt(293/Tn);
micro_a2_H=micro_a_H/1.128*sqrt(293/Tn);
micro_a2_O=micro_a_O/1.128*sqrt(293/Tn);
micro_f2_5=micro_f_5/1.128*sqrt(293/Tn)*gf_5;
%堆芯的热中子宏观吸收截面
macro_a2=N_5*micro_a2_5+N_8*micro_a2_8+N_H*micro_a2_H+N_O*micro_a2_O+N_Zr*micro_a2_Zr
%堆芯的热中子宏观裂变截面
macro_f2=N_5*micro_f2_5
%L_H2O
L_H2O=2.719/den_H20_311*((t+273)/295)^(1/2)*(1.0614/(1+0.0614*((t+273)/295)^(1/2)))^(1/2)
%水的热群扩散系数
D2_H2O=L_H2O^2*(N_H*micro_a2_H+N_H2O*V_H2O/V*micro_a2_O)
macr_tr2_H2O=1/(3*D2_H2O)
问题是:为什么从
L_sigama_th=N_5*micro_a_5+N_8*micro_a_8+N_O*micro_a_O+N_H*micro_a_H+N_Zr*micro_a_Zr
L_sigama_kTM=L_sigama_th*sqrt(293/TM) %式(2-80)
Tn=TM*(1+1.4*L_sigama_kTM/kesi_macr) %式(2-79)
开始计算的结果如下呢:
L_sigama_kTM =
.30589717806208658168495092291557e-5*17115595^(1/2)
Tn =
11683/20+.45936265372953731364559710690310e-1*17115595^(1/2)
ga_5 =
6493/6250/(6223/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(23/1250)
ga_8 =
1971/2000*(6223/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(43/10000)
gf_5 =
1053/1000/(6223/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(113/5000)
macro_a2 =
.13711453320619852340828168976573e-1*293^(1/2)/(11683/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(1/2)/(6223/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(23/1250)+.17678485256047348144979250378501e-2*293^(1/2)/(11683/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(1/2)*(6223/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(43/10000)+.84913479087977694805938554307674e-3*293^(1/2)/(11683/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(1/2)
macro_f2 =
.11909787189519441142160261181845e-1*293^(1/2)/(11683/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(1/2)/(6223/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(113/5000)
L_H2O =
0.0054
D2_H2O =
69386399805210663738033738022873015464822941615415/3138550867693340381917894711603833208051177722232017256448*293^(1/2)/(11683/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(1/2)
macr_tr2_H2O =
3138550867693340381917894711603833208051177722232017256448/60990645428780173425731655722105380593579365679949785*293^(1/2)*(11683/20+.45936265372953731364559710690310e-1*17115595^(1/2))^(1/2)