核反应堆物理分析课后习题参考答案

核反应堆物理分析答案第一章1-1.某压水堆采用UO2作燃料，其富集度为2.43%（质量），密度为10000kg/m3。试计算：当中子能量为0.0253eV时，UO2的宏观吸收截面和宏观裂变截面。解：由18页表1-3查得，0.0253eV时：由289页附录3查得，0.0253eV时：以c5表示富集铀内U-235与U的核子数之比，表示富集度，则有：所以，1-2.某反应堆堆芯由U-235,H2O和Al组成，各元素所占体积比分别为0.002,0.6和0.398，计算堆芯的总吸收截面(E=0.0253eV)。解：由18页表1-3查得，0.0253eV时：由289页附录3查得，0.0253eV时：可得天然U核子数密度则纯U-235的宏观吸收截面：总的宏观吸收截面：1-61-12题每秒钟发出的热量：每秒钟裂变的U235：运行一年的裂变的U235：消耗的u235质量：需消耗的煤：1-10.为使铀的η＝1.7，试求铀中U-235富集度应为多少(E=0.0253eV)。解：由18页表1-3查得，0.0253eV时：由定义易得：为使铀的η＝1.7，富集度.一核电站以富集度20%的U-235为燃料，热功率900MW,年负荷因子(实际年发电量/额定年发电量)为0.85,U-235的俘获－裂变比取0.169,试计算其一年消耗的核燃料质量。解：该电站一年释放出的总能量=对应总的裂变反应数=因为对核燃料而言：核燃料总的核反应次数=消耗的U-235质量=消耗的核燃料质量=第二章.某裂变堆，快中子增殖因数1.05，逃脱共振俘获概率0.9，慢化不泄漏概率0.952，扩散不泄漏概率0.94，有效裂变中子数1.335，热中子利用系数0.882，试计算其有效增殖因数和无限介质增殖因数。解：无限介质增殖因数：有效增殖因数：不泄漏概率：2-1.H和O在1000eV到1eV能量范围内的散射截面近似为常数，分别为20b和38b。计算H2O的ξ以及在H2O中中子从1000eV慢化到1eV所需的平均碰撞次数。解：不难得出，H2O的散射截面与平均对数能降应有下述关系：σH2O∙ξH2O=2σH∙ξH+σO∙ξO即：(2σH+σO)∙ξH2O=2σH∙ξH+σO∙ξOξH2O=（2σH∙ξH+σO∙ξO）/(2σH+σO)查附录3，可知平均对数能降：ξH=1.000，ξO=0.120，代入计算得：ξH2O=(2×20×1.000+38×0.120)/(2×20+38)=0.571可得平均碰撞次数：Nc=ln(E2/E1)/ξH2O=ln(1000/1)/0.571=12.09≈12.12-6.在讨论中子热化时，认为热中子源项Q(E)是从某给定分界能Ec以上能区的中子，经过弹性散射慢化而来的。设慢化能谱服从Ф(E)=Ф/E分布，试求在氢介质内每秒每单位体积内由Ec以上能区，（1）散射到能量E（E<Ec）的单位能量间隔内之中子数Q(E)；（2）散射到能量区间ΔEg=Eg-1-Eg内的中子数Qg。解：（1）由题意可知：对于氢介质而言，一次碰撞就足以使中子越过中能区，可以认为宏观截面为常数：在质心系下，利用各向同性散射函数：。已知，有：（这里隐含一个前提：E/α>E’）（2）利用上一问的结论：2-8.计算温度为535.5K，密度为0.802×103kg/m3的H2O的热中子平均宏观吸收截面。解：已知H2O的相关参数，M=18.015g/mol，ρ=0.802×103kg/m3，可得：m-3已知玻尔兹曼常数k=1.38×10-23J•K-1，则：kTM=1.38×10-23×535.5=739.0(J)=0.4619(eV)查附录3，得热中子对应能量下，σa=0.664b，ξ=0.948，σs=103b，σa=0.664b，由“1/v”律：0.4914(b)由56页（2-81）式，中子温度：577.8(K)对于这种”1/v”介质，有：n0.4192(b)所以：1.123(m-1)三章3.1有两束方向相反的平行热中子束射到235U薄片上，设其上某点自左面入射的中子束强度为1012cm-2·s-1。自右面入射的中子束强度2×1012cm-2·s-1。计算：（1）该点的中子通量密度；（2）该点的中子流密度；（3）设Σa=19.2×102m-1，求该点的吸收率。解：（1）由定义可知：（2）若以向右为正方向：3×1012(cm-2·s-1)-1×1012(cm-2·s-1)可见其方向垂直于薄片表面向左。（3）19.2•3×1012=5.76×1013(cm-3·s-1)3.2设在x处中子密度的分布函数是其中：λ，ɑ为常数，μ是与x轴的夹角。求：（1）中子总密度n(x)；（2）与能量相关的中子通量密度φ(x,E)；（3）中子流密度J(x,E)。解：由于此处中子密度只与夹角φ为方向角，则有：与x轴的夹角有关，不妨视μ为极角，定义在Y-Z平面的投影上与Z轴的（1）根据定义：可见，上式可积的前提应保证ɑ<0，则有：（2）令mn为中子质量，则（等价性证明：如果不作坐标变换，则依据投影关系可得：则涉及角通量的、关于空间角的积分：对比：可知两种方法的等价性。）（3）根据定义式：利用不定积分：（其中n为正整数），则：3.7设一立方体反应堆，边长ɑ=9m。中子通量密度分布为已知D=0.84×10-2m，L=0.175m。试求：（1）表达式；（2）从两端及侧面每秒泄漏的中子数；（3）每秒被吸收的中子数（设外推距离很小可略去）。解：有必要将坐标原点取在立方体的几何中心，以保证中子通量始终为正。为简化表达式起见，不妨设φ0=3×1013cm-2•s-1。（1）利用Fick’sLaw：（2）先计算上端面的泄漏率：同理可得，六个面上总的泄漏率为：L=1.7×1017(s-1)其中，两端面的泄漏率为L/3=5.8×1016-1(如果有同学把问题理解成‘六个面’上总的泄漏，也不算错)(s)；侧面的泄漏率为L-L/3=1.2×1017(s-1)（3）由可得由于外推距离可忽略，只考虑堆体积内的吸收反应率：1.24×1020(s-1)3.8圆柱体裸堆内中子通量密度分布为其中，H，R为反应堆的高度和半径（假定外推距离可略去不计）。试求：（1）径向和轴向的平均中子通量密度与最大中子通量密度之比；（2）每秒从堆侧表面和两个端面泄漏的中子数；（3）设H=7m，R=3m，反应堆功率为10MW，σf,5=410b，求反应堆内235U的装载量。解：有必要将坐标原点取在圆柱体的几何中心，以保证中子通量始终为正。为简化表达式起见，不妨设φ0=1012cm-2•s-1。且借用上一题的D值。（1）先考虑轴向：且在整个堆内只在z=0时为0，故有：径向：且在整个堆内只在r=0时为0，故有：已知，所以：0.611（2）先计算上端面的泄漏率：易知，两端面总泄漏率为侧面泄漏率：2.93×1014(s-1)利用Bessel函数微分关系式：，且已知J1(2.405)=0.5191，可得：所以：4.68×1014(s-1)（3）已知每次裂变释能(J)所以：其中：利用Bessel函数的积分关系式：，可得已知：J1(0)=0，J1(2.405)=0.5191，所以：=5.44×1017(m•s-1)所以：106/(3.2×10-11×410×10-28×5.44×1017)=1.40×1024(m-3)所需235U装载量：10-3×1.40×1024×3.14×32×7×235/(6.02×1023)=108(kg)3.9试计算E=0.025eV时的铍和石墨的扩散系数。解：查附录3可得，对于E=0.025eV的中子：/m-1BeC8.653.850.92590.9444对于Be：同理可得，对于C：D=0.0917(m)3-12试计算T=535K，ρ=802kg/m3时水的热中子扩散系数和扩散长度。0.0416(m)解：查79页表3-2可得，294K时：m，由定义可知：所以：0.00195(m)（另一种方法：如果近似认为水的微观散射截面在热能区为常数，且不受温度影响，查附表3可得：在T=535K，ρ=802kg/m3时，水的分子数密度：103×802×6.02×1023/18=2.68×1028(m-3)所以：276(m-1)1/（3×2.68×103×0.676）=0.00179(m)这一结果只能作为近似值）中子温度利用56页（2-81）式计算：其中，介质吸收截面在中子能量等于kTM=7.28×1021J=0.0461eV再利用“1/v”律：0.4920(b)Tn=535×(1+0.46×36×0.4920/103)=577(K)（若认为其值与在0.0253eV时的值相差不大，直接用0.0253eV热中子数据计算：Tn=535×(1+0.46×36×0.664/103)=592(K)这是一种近似结果）（另一种方法：查79页表3-2，利用293K时的平均宏观吸收截面与平均散射截面：(m-1)1/（3×0.0016×0.676）=308(m-1)进而可得到Tn=592K）利用57页（2-88）式0.414×10-28(m2)1.11(m-1)802/(3×1000×0.0016×0.676)=247(m-1)0.0424(m)（此题如果利用79页（3-77）式来计算：由于水是“1/v”介质，非1/v修正因子为1：代入中子温度可得：0.0340(m)这是错误的！因为（3-74）式是在（3-76）式基础上导出的，而（3-76）式是栅格的计算公式，其前提是核子数密度不随温度变化）3.13如图3-15所示，在无限介质内有两个源强为Ss-1的点源，试求P1和P2点的中子通量密度和中子流密度。解：按图示定义平面坐标。O假设该介质无吸收、无散射，则在P2点，来自左右两个点源的中子束流强度均为I+=I-=S/4πa2，可知：在P1点，来自左右两个点源的中子束流强度均为，且其水平方向的投影分量恰好大小相等、方向相反，可得：其方向沿Y轴正向。若考虑介质对中子的吸收及散射，设总反应截面为，则上述结果变为：（注意：如果有同学用解扩散方程的方法，在有限远处的通量密度同时与x、y、z有关。）3-16设有一强度为I（m-2•s-1）的平行中子束入射到厚度为a的无限平板层上。试求：（1）中子不遭受碰撞而穿过平板的概率；（2）平板内中子通量密度的分布；（3）中子最终扩散穿过平板的概率。解：（1）（2）此情况相当于一侧有强度为I的源，建立以该侧所在横坐标为x原点的一维坐标系，则扩散方程为：边界条件：i.ii.方程普遍解为：由边界条件i可得：由边界条件ii可得：所以：（也可使用双曲函数形式：方程普遍解为：由边界条件i可得：由边界条件ii可得：所以：可以证明这两种解的形式是等价的）（3）此问相当于求x=a处单位面积的泄漏率与源强之比：（或用双曲函数形式：）3-17设有如图3-16所示的单位平板状“燃料栅元”，燃料厚度为2a，栅元厚度为2b，假定热中子在慢化剂内以均匀分布源（源强为S）出现。在栅元边界上的中子流为零（即假定栅元之间没有中子的净转移）。试求：（1）屏蔽因子Q，其定义为燃料表面上的中子通量密度与燃料内平均中子通量密度之比；（2）中子被燃料吸收的份额。解：（1）以栅元几何中线对应的横坐标点为原点，建立一维横坐标系。在这样对称的几何条件下，对于所要解决的问题，我们只需对x>0的区域进行讨论。燃料内的单能中子扩散方程：边界条件：i.ii.通解形式为：利用Fick’sLaw：代入边界条件i：代入边界条件ii：所以（2）把该问题理解为“燃料内中子吸收率/燃料和慢化剂内总的中子吸收率”，设燃料和慢化剂的宏观吸收截面分别为和，则有：回顾扩散长度的定义，可知：，所以上式化为：（这里是将慢化剂中的通量视为处处相同，大小为S，其在b处的流密度自然为0，但在a处情况特殊：如果认为其流密度也为0，就会导致没有向燃料内的净流动、进而燃料内通量为0这一结论！所以对于这一极度简化的模型，应理解其求解的目的，不要严格追究每个细节。）3-21解：（1）建立以无限介质内任一点为原点的球坐标系（对此问题表达式较简单），建立扩散方程：即：边界条件：i.,ii.设存在连续函数满足：可见，函数满足方程，其通解形式：由条件i可知：C=0，由方程（2）可得：再由条件ii可知：A=0，所以：（实际上，可直接由物理模型的特点看出通量处处相等这一结论，进而其梯度为0）（2）此时须以吸收片中线上任一点为原点建立一维直角坐标系，先考虑正半轴，建立扩散方程：即：，x>0边界条件：i.,ii.,iii.对于此“薄”吸收片，可以忽略其厚度内通量的畸变。参考上一问中间过程，可得通解形式：由条件ii可得：由条件iii可得：C=0所以：对于整个坐标轴，只须将式中坐标加上绝对值号，证毕。3-22解：以源平面任一点为原点建立一维直角坐标系，建立扩散方程：边界条件：i.;ii.;iii.;iv.;通解形式：，由条件i：由条件ii：（1）（2）由条件iii、iv：（3）（4）联系（1）可得：结合（2）可得：所以：3-23证明：以平板中线上任一点为原点建立一维直角坐标系，先考虑正半轴，建立扩散方程：即：，x>0边界条件：i.,ii.,iii.参考21题，可得通解形式：由条件ii可得：再由条件iii可得：所以：由于反曲余弦为偶函数，该解的形式对于整个坐标轴都是适用的。证毕。3-24设半径为R的均匀球体内，每秒每单位体积均匀产生S个中子，试求球体内的中子通量密度分布。解：以球心为原点建立球坐标系，建立扩散方程：即：边界条件：i.,ii..,iii.通解：由条件iii：再由条件ii：所以：（此时，）第四章4-1试求边长为a，b，c（包括外推距离）的长方体裸堆的几何曲率和中子通量密度分布。设有一边长a=b=c=0.5m，c=0.6m（包括外推距离）的长方体裸堆，L=0.0434m，τ=6cm2。（1）求达到临界时所必须的k∞；（2）如果功率为5000kW，Σf=4.01m-1，求中子通量密度分布。解：长方体的几何中心为原点建立坐标系，则单群稳态扩散方程为：边界条件：（以下解题过程中不再强调外推距离，可以认为所有外边界尺寸已包含了外推距离）因为三个方向的通量变化是相互独立的，利用分离变量法：将方程化为：设：先考虑x方向，利用通解：代入边界条件：同理可得：其中φ0是待定常数。其几何曲率：106.4(m-2)（1）应用修正单群理论，临界条件变为：其中：0.00248(m2)1.264（2）只须求出通量表达式中的常系数φ01.007×1018(m-2•s-1)4-2设一重水-铀反应堆堆芯的k∞=1.28，L2=1.8×10-2m2，τ=1.20×10-2m2。试按单群理论，修正单群理论的临界方程分别求出该芯部材料曲率和达到临界时总的中子不泄漏概率。解：对于单群理论：15.56(m-2)在临界条件下：0.7813(或用)对于单群修正理论：0.03(m2)9.33(m-2)在临界条件下：0.68\0.7813?（注意：这时仍能用，实际上在维持临界的前提条件下修正理论不会对不泄漏概率产生影响，但此时的几何曲率、几何尺寸已发生了变化，不再是之前的系统了）4-4解：=4.79×1024(m-3)，4.79×1028(m-3)堆总吸收截面：总裂变截面：=0.344(m-1)=0.280(m-1)=2.61×10-2(m2)=1.97则材料曲率：=37.3(m-2)在临界条件下：=0.514(m)考虑到外推距离：（如有同学用=0.018(m)也是正确的，但表达式相对复杂）再考虑到堆的平均密度：（或者由=957(kg/m3)）实际的临界质量：=156(kg)4-5证明：以球心为坐标原点建立球坐标系，单群稳态扩散方程：边界条件：i.；ii.；（如果不认为R2包括了外推距离的话，所得结果将与题意相悖）球域内方程通解：由条件i可得：由条件ii可得：由此可见，，证毕4-7一由纯235U金属（ρ=18.7×103kg/m3）组成的球形快中子堆，其周围包以无限厚的纯238U（ρ=19.0×103kg/m3），试用单群理论计算其临界质量，单群常数如下：235U：σf=1.5b,σa=1.78b,Σtr=35.4m-1,ν=2.51；238U：σf=0,σa=0.18b,Σtr=35.4m-1。解：以球心为坐标原点建立球坐标系，对于U-235和U-238分别列单群稳态扩散方程，设其分界面在半径为R处：U-235：U-238：方程1方程2边界条件：i.ii.iii.iv.令（在此临界条件下，既等于材料曲率，也等于几何曲率），球域内方程1通解：由条件i可知A5=0，所以：球域内方程2通解：由条件iv可知C8=0，所以：由条件ii可得：由条件iii可得：所以（由题目已知参数：）即：代入数据：4.79×10-28(m-3)4.81×10-28(m-3)2.1151.31×10-3(m2)29.17(m-1)0.1043(m)0.06474(m)21.3(kg)4-8证明：（1）如图4-8所示的柱坐标系下，单群稳态扩散方程可写为（临界条件下，几何曲率与材料曲率相等）：，（）边界条件（不考虑外推距离）：i.ii.iii.（注意，这里不能用线性微分方程解的存在唯一性定理：如果都是区间上的连续函数，则对于任一，及任意的，方程存在唯一解定义于区间上，且满足初值条件而此扩散方程并非线性微分方程。）对于表达式：，不难证明其满足上述全部三个边界条件。（）（2）将表达式代入方程，其中，已知如下关系：可推得：所以：所以：再有：所以方程化为：可知该表达式为方程的解。证毕。（也可如此推出解的形式：分离变量：方程变形：设：（n为任意实数），：变量替换：，此为n阶Bessel方程，通解为：由边界条件i可得，n须取使的值，在其中，我们只取基波，即n=1，相应的：相应的：由边界条件ii可得，对于z有：，由边界条件ii可得，所以：）4-10解：（1）对于均匀圆柱体裸堆，其几何曲率：可得，在临界条件下：临界体积：其取最小值时：，即：所以：0.5412（2）由上可得临界最小体积：由于临界条件下：所以：4-11设有一由纯239Pu（ρ=14.4×103kg/m3）组成的球形快中子临界裸堆，试用下列单群常数：ν=2.19,σf=1.85b,σr=0.26b,σtr=6.8b计算其临界半径与临界质量。解：4-11解：由已知条件可得：3.64×1028(m-3)1.921.77×10-3(m2)设临界半径为R，则由临界条件：，可得：0.138(m)对于这一实际问题，需考虑外推距离：0.0288(m)所以实际临界体积为：5.40×10-3(m3)临界质量：77.8(kg)4-12试求下列等效裸堆内热中子通量密度的最大值与平均值之比，即热中子通量密度的不均匀系数：（1）半径为R的球形堆，反射层节省为δT；（2）半径为R，高度为H的圆柱体堆，反射层节省分别为δr和δH；（3）边长为a，b，c的长方体堆，反射层节省分别为δx，δy，δz。解：可利用裸堆结论：球：圆柱：立方体：详细推导：根据97页表4-1裸堆的通解形式可得：球：圆柱：（与教材上数值的差异在于对立方体：所取的近似值的不同，在此取的是0.5191）4-16解：以平板厚度方向上的几何中心为原点建立坐标系，对两区分别建立单群稳态扩散方程（由于几何上的对称性，对于本题只需考虑一侧，如x为正一侧）：方程1方程2边界条件：i.;ii.由表3-1查得方程1的通解：其中第二项明显有悖于对称性条件，故CI=0，同理有：（由于本题是求解临界尺寸，默认的前提是几何曲率等于材料曲率，故以下不再对其进行区别，统一用B2表示）由条件ii可得：整个系统的临界条件为：即：（注意，此处的泄漏仅仅是II区外表面上的泄漏，I-II区之间的净流动是