1、第四章1.试求边长为 （包括外推距离）的长方体裸堆的几何曲率和中子通量密度的分布。,abc设有一边长 （包括外推距离）的长方体裸堆， 0.5,.6m 0.43,Lm。 （1）求达到临界时所必须的 ；（2 ）如果功率为426k，求中子通量密度分布。50,.fkW解：长方体的几何中心为原点建立坐标系，则单群稳态扩散方程为： 22()0aaDkxyz边界条件： /,)(,/(,/2)abxyc（以下解题过程都不再强调外推距离，可认为所有外边界尺寸已包含了外推距离）因为三个方向的通量拜年话是相互独立的，利用分离变量法： (,)()(xyzXYZz将方程化为：2221kXL设：222,xyzYBB想考虑
2、 X 方向，利用通解： ()cosinxxAC代入边界条件： 1cos0,3.5,xn xaABaa同理可得： (,)()s()(yzyz其中 是待定常数。0其几何曲率： 222()()106.4gBmabc（1 ）应用修正单群理论，临界条件变为： 2gkBM其中： 220.48ML1.64k（2 ）只须求出通量表达式中的常系数 032220 02cos()cos()cos()()abfff fVPEdxdydzdEabc 318202().7f mabc2.设一重水铀反应堆的堆芯 。试按单群理222.,.0,1.0kLm论，修正单群理论的临界方程分别求出该芯部的材料曲率和达到临界时候的总的中
3、子不泄露几率。解：对于单群理论：在临界条件下： 2210.7813gmBL（或用 ）k对于单群修正理论： 22.M193在临界条件下： 220.7813gmB（注意：这时能用 ，实际上在维持临界的前提条件下修正理论不1k会对不泄露几率产生影响，但此时的几何曲率、几何尺寸已发生了变化，不再是之前的系统了。 ）4. 设有圆柱形铀-水栅装置，R=0.50 米，水位高度 H=1.0 米，设栅格参数为：k =1.19，L 2=6.610-4米 2，=0.5010 -2米 2。 （a）试求该装置的有效增殖系数k；（b）当该装置恰好达临界时，水位高度 H 等于多少？（c）设某压水堆以该铀-水栅格作为芯部，堆
4、芯的尺寸为 R=1.66 米，H=3.50 米，若反射层节省估算为 r=0.07 米， H=0.1米。试求反应堆的初始反应性 以及快中子不泄漏几率和热中子不泄漏几率。5.一个球壳形反应堆，内半径为 ,外半径为 ，如果球的内、外均为真空，求证单群理论1R2的临界条件为： 12tantBR解答：以球心为坐标原点建立球坐标系，单群稳态扩散方程：22r边界条件：i. 1lim0;xRJii. 2()（如果不 包括了外推距离的话，所得结果将与题意相悖）球域内方程通解： cosin()BrrrAC由条件 i 可得：11 11112 2isicoslm 0rRrR RBRJDC111cosintaniBCA
5、B由条件 ii 可得： 由此可见， ，证毕。12tantBR7.一由纯 金属 组成的球形快中子堆，其周围包以无限厚的纯235U3(18.70/)kgm，试用单群理论计算其临界质量，单群常数如下：2383(9.0/)k。5 1238 1:,.,5.4,.;:0,.18,35.4fatr fatrbvUbm 解：以球心为左边原点建立球左边系，对于 U-235 和 U-238 分别列单群稳态扩散方程，设其分界面在半径为 R 处：方程 125523:kUL方程 281边界条件：i. ii. 50limr58()Riii. iv. 8rRrRDlim0r令 （.在此临界条件下，既等于材料曲率，也等于几何
6、曲率） ，球域内251kBL方程 1 通解： 55cosin()BrrrAC由条件 i 可知 ，所以：50s()球域内方程 2 通解： 8888exp/exp(/)rLrLr由条件 iv 可知,所以： ()()A由条件 ii 可得： 88e/e(/)sinsinRRBCCAB由条件 iii 可得：885 8285(1)exp()cosin1()()exp()sincosRDLRDBDALRLB所以（由题目已知参数 ）,5,58,5,13trtrtrtr8885 8(1)exp()exp(/)sincos(1)sinsincossiRLRLRAABBBDL即： 8inB8cot(1/)1coss
7、arRL代入数据：328358104.7910ANmM8.,5,52 325,5,21258,5,8835.11.039.70.433cot(1/)/2arctn(1/)0.647.ffaafafvkLmkBLrLBLR mBmVRkg8.试证明有限高半圆形反应堆中子通量密度分布和几何曲率 1(,)()sinco()xrzrzAJH22gBR其中： 是 的第一个零点，即。13.89x1()Jx证明：（1）书上图 4-8 所示的柱坐标系下，单群稳态扩散方程可写为（临界条件下，几何曲率与材料曲率相等）：222,(0,/2/)gBrHZrrz边界条件（不考虑外推距离）：i. 0rRrII. III.
8、 /2/2zHz(注意,这里不能用线性微分方程解的存在唯一性定理：如果 都是区间 上的连续函数，则对于任)(1,2),(iatnft ,ab一 及任意的 方程：0b01)(2(1)0,nxx()()1nnxaft存在唯一解 ()t定义于区间 上，且满足初值条件,b()()0,1),kkxtn而此扩散方程并非线性微分方程。 ）对于表达式： 11(,)()sinco(),3.89rzrzAJRH不难证明其满足上述全部三个边界条件。 )0J（2 ）将表达式代入方程，其中，已知如下条件：101,nnxJxJ可推得： 10 0010011101 12 222()JxJJxJJxx x 1110()()(
9、)()rr rRR12 22111 111021 ()()()xJx xrxrJ JxrrRR 所以：22111021()()xJJrRRxr101()()xJJrrR212()xJrr所以： 1122211 1002221()()1()xrxrJJxrxRRJrRr 再有：22 2cos()(zFHz所以方程为：221gxBR可知该表达式为方程的解。证毕。（也可如此推出解的形式：分离变量： (,)()(rzrQZz方程变形：2221gddQZrBr设： （ 为任意实数） ， ；22dQn22zdZB22 22 21 ()0gzr rddrBnr变量替换： 2,(),rxxx此为 阶方程，通解
10、为nesl()()nnrJYxB由边界条件 i 可得,n 须取使 的值，在其中，我们只去基波，即 ，0n1n相应的 ：1rBR1()/)rJxR相应的： ()siiQAC由边界条件 ii 可得： 0,sinA对于 z 有： in()co()zzzzZB由边界条件 ii 可得， ,/,cos(/)zHZCH所以: 1/siJxrR10.设有均匀圆柱形裸堆，其材料曲率等于，试求：（1）使临界体积为最小的 的值；/（2）最小临界体积 V 与 的关系。2mB解：（1）对于均匀圆柱体裸堆，其几何曲率： 222.405()gBHR可得，在临界条件下： 22.405()gR临界体积：3222.gVHB其取最
11、小值时: ,即：0d22232222.4053.45 33()() ggggggHHBHB 22 2.0.4053/ggggRRB所以： 2.45.1（2）由上可得临界最小体积：222 3.4053.405gggVRHBB由于临界条件下： ，所以：2gmB18./m11.设有意纯 组成的球形快中子临界裸堆，试用下列单群常数：239Pu3(14.0/)k计算其临界半径与临界质量。.1,85,6,.frtrvbb解：由已知条件可得：32836410ANM1.92ffarvk2 321.703()trtrfrDL m设临界半径为 ，则临界条件： ，可得：R2gmB210.381kLHk对于这一实际问
12、题，需要考虑外推距离： 0.74.1.8trtrdN所以实际临界体积为： 334()5.401VRdm临界质量： 7.mkg12.试求下列等效裸堆内热中子通量密度的最大值与平均值，即热中子通量密度的不均匀系数:（1）半径为 的球形堆，反射层节省为 ；RT（2）半径为 ，高度为 的圆柱形堆，反射层节省分别为 和 ；HrH（3）边长为 的长方形堆，反射层节省分别为 。,abc,xyz解：可利用裸堆的结论，球： 3, 204/.7sin()HbareRKrd3T圆柱：2,2/ 03.62.45cos()()Hbare RHzdJrd23.6rHK立方体：3,/2/2/2.8cos()cos()cos
13、()Hbarebaxdydzd 3()()()822HxyzaaK详细推导：据 97 页 4-1 裸堆的通解形式可得：球： 1()sin()TrArRmax0 0cos()lisi()lim1Tr rTT TrRAA 34/VR200sinsin()TRTdAdrrd(co)icos()T TrR |004s()()cos()TTRT TRArrdrR |22 20()()cos(4()TTTxA 323max4()1()TH TVRAKd圆柱： 02.5(,)cos()2TzrzJrRHmax0.4lir zzAA2VH2/00.5()cos()rRHTzdJrdd 2/12.4() in.
14、5rR HT zzA |20.920.8637()(). z TzAR 2 2max .4()()1.8637()()HTzTzVRHKAHd立方体： (,)coscos)cos()2xyzyzaaamax0licos()cos()cos()222xyzyzAAaaa Vb/2/2/2cos()cos()cos()x y zabVx y ydddzd 3() 2yzxyzAAabmax3 ()()()128()(2)H xyzxyzVabc cKd 16.设有如图 4-9 所示的 一维无限平板反应堆。中间区域（ ）的 ，厚度 为已知，I1Ik2b两侧区域（ ）的 ，试用单群理论导出确定临界尺寸
15、 的公式及临界时中子通I1Ika量密度的分布。说明尺寸 对临界尺寸有无影响及其理由。b解：以平板厚度方向上的几何中心为原点建立坐标系，对两区分别建立单群稳态扩散方程（由于几何上的对称性，对于本体只需考虑一侧，如 X 为正一侧）：方程 122,0IIIkxxL方程 221,IIIba边界条件：i. ii. ();II()0Ib由表 3-1 查得方程 1 的通解： )cosinIIIIxABxCx其中第二项明显有悖于对称性条件，故 ，同理有：()cosIIIxABx（由于本体是求解临界尺寸，默认的前提是几何曲率等于材料曲率，故以下不再对其进行区别，统一用 表示）2有条件 ii 可得： cos()02()II IbaBab整个系统的临界条件为：=中子率/（中子泄漏率 +中子吸收率）=1efk即： 020 1bbaIIIIII affIef I IaaIIaxabIb bkdxkdxvRdVvJSRdV 20 0bb aIIIIIIIaaab bxxkxkx20(1)(1)(1)abb aIIIIIIII aab bDBdkddd ()/Ik
