目录一、设计任务 (1)二、课程设计要求 (2)三、计算过程 (2)四、程序设计框图 (8)五、代码说明书 (9)六、热工设计准则和出错矫正 (10)七、重要的核心程序代码 (11)八、计算结果及分析 (17)一、设计任务某压水反应堆的冷却剂及慢化剂都是水，用二氧化铀作燃料，用Zr-4作包壳材料。

燃料组件无盒壁，燃料元件为棒状，正方形排列。

已知下列参数：系统压力 15.8MPa堆芯输出功率 1820MW冷却剂总流量 32100t/h反应堆进口温度287℃堆芯高度 3.66m燃料组件数 121燃料组件形式17×17每个组件燃料棒数 265燃料包壳直径 9.5mm燃料包壳内径 8.36mm燃料包壳厚度 0.57mm燃料芯块直径 8.19mm燃料棒间距（栅距） 12.6mm芯块密度 95%理论密度旁流系数 5%燃料元件发热占总发热的份额 97.4%径向核热管因子 1.35轴向核热管因子 1.528局部峰核热管因子 1.11交混因子 0.95热流量工程热点因子 1.03焓升工程热管因子 1.085堆芯入口局部阻力系数 0.75堆芯出口局部阻力系数 1.0堆芯定位隔架局部阻力系数 1.05若将堆芯自上而下划分为5个控制体，则其轴向归一化功率分布如下表：堆芯轴向归一化功率分布（轴向等分5个控制体）通过计算，得出1. 堆芯出口温度；2. 燃料棒表面平均热流及最大热流密度，平均线功率，最大线功率；3. 热管的焓，包壳表面温度，芯块中心温度随轴向的分布；4. 包壳表面最高温度，芯块中心最高温度；5. DNBR在轴向上的变化；6. 计算堆芯压降；二、课程设计要求1．设计时间为两周；2．独立编制程序计算；3．迭代误差为0.1%；4．计算机绘图；5．设计报告写作认真，条理清楚，页面整洁；6．设计报告中要附源程序。

三、计算过程目前，压水核反应堆的稳态热工设计准则有：（1）燃料元件芯块内最高温度应低于其相应燃耗下的熔化温度。

目前，压水堆大多采用UO2作为燃料。

二氧化铀的熔点约为2805 ±15℃，经辐照后，其熔点会有所降低。

燃耗每增加104兆瓦·日/吨铀，其熔点下降32℃。

在通常所达到的燃耗深度下，熔点将降至2650℃左右。

在稳态热工设计中，一般将燃料元件中心最高温度限制在2200～2450℃之间。

（2）燃料元件外表面不允许发生沸腾临界。

通常用临界热流密度比DNBR 来定量地表示这个限制条件。

DNBR 是根据堆内某处燃料元件周围的冷却剂状态使用专门的计算公式而得到的临界热流密度与该处燃料元件表面的实际热流密度的比值。

DNBR 随堆芯通道的长度是变化的，在整个堆芯内，DNBR 的最小值称为最小DNBR ，用MDNBR 或DNBRmin 表示。

为了确保燃料元件不烧毁，当计算的最大热功率下，MDNBR 不应低于某一规定值。

如果计算热流密度的公式没有误差，则当MDNBR=1 时，表示燃料元件表面要发生沸腾临界。

若该公式存在误差，则MDNBR 就要大于1。

例如，W-3 公式的误差为23%，所以当使用W-3 公式计算DNBR 时，就要求MDNBR≥ 1.3。

（3） 在稳态额定工况下，要求在计算的最大热功率下，不发生流动不稳定性。

对于压水堆，只要在堆芯最热通道出口附近冷却剂中的含气量不大于某一数值，就不会发生流动不稳定性。

（4） 必须保证正常运行工况下燃料元件和堆内构件能得到充分冷却；在事故工况下能提供足够的冷却剂以排出堆芯余热。

本计算根据《核反应堆热工分析》课程设计指导书中的计算提示，采用简单的C 语言编程计算。

将堆芯沿轴向划分为五个等分控制体进行计算以下是计算过程:3.1堆芯流体出口温度（平均管）℃按流体平均温度以及压力由表中查得。

3.2燃料表面平均热流密度W/m 2式中为堆芯燃料棒的总传热面积m 2,,..(1).f out f inFa Ntt t W Cpς=+-Cp ,,1()2f f in f out t t t =+q ./q Fa Nt F =总F 总=....cs F m n d Lπ总燃料棒表面最大热流密度q maxw/m 2燃料棒平均线功率W/m燃料棒最大线功率w/m3.3平均管的情况 平均管的流速Vm/s式中，堆芯内总流通面积n 0为燃料组件内正方形排列时的每一排（列）的燃料元件数由压力以及流体的平均温度查表得到：3.4为简化计算起见，假定热管内的流体流速V h 和平均管的V 相同。

（实际上，应该按照压降相等来求。

热管内的流体流速要小一些）。

则V h =V同样，热管四根燃料元件组成的单元通道内的流量3.5热管中的计算（按一个单元通道计算）（1）热管中的流体温度max q q ..N Eq q F F =l q .....cs l csq d Lq q d L ππ==,maxl q ,max ..N El l q q q q F F =(1)..f feff f W W V A A ςρρ-==22000.(.)4(.)42f cs A m n n p d m n s πδ=⎡⎤⎡⎤-+⨯⎢⎥⎢⎥⎣⎦⎣⎦fρft 1f fV ρ=(1)h hf W W A A ς-=224b csA S d π=-（2）第一个控制体出口处的包壳外壁温度式中：h(z)可以用来求。

所以，式中：流体的k(z)、μ（z ）和Pr 数根据流体的压力好温度由表查得。

（k=λ 传热系数）如果流体已经达到过冷沸腾，用Jens-Lottes 公式：.,(),0..(.)()N E EzR H H m csf h z f in zh q F F F d t t z d W Cpπϕ∆∆=+⎰1,(),(),()cs h z f h z f h z t t θ=+∆,()..()N Eq q f h z q F F t h z =+,()..().()N ER q f h z q F z F t h z ϕ=+0.80.4()0.023Re Pr ()eh z D Nu k z ==0.80.4()()0.023Re Pr .e k z h z D =.Re .eh e b G D W D A μμ==224()44cs be csS d A D U d ππ-==0.256.26..().25.10p NE Rqw s q F z Ft t eϕ⎛⎫- ⎪⎝⎭⎛⎫-= ⎪ ⎪⎝⎭()()()()()2,(),(),f h z w f h z w s s f h z t z t t z t t t θ∆=-=-+-当时，用前面的式子当时，用（3）第一个控制体出口处的包壳内壁温度式中：Zr-4的W/m.℃(4)第一个控制体出口处的UO 2芯块外表面温度（5）第一个控制体出口处的UO 2芯块中心温度用积分热导求解的方法，即其他2个控制体的计算方法相同，重复上述过程即可。

3.6热管中的用w-3公式计算，同样对3个控制体都算()()()0.256.2,6..25.10p NE Rqs f h z q F Z F et t ϕ⎛⎫- ⎪⎝⎭⎛⎫=+- ⎪ ⎪⎝⎭()()22,,f h z f h z θθ∆∆p ()()22,,f h z f h z θθ∆∆f ()()21,,f h z f h z θθ∆∆代替()()(),,,t ci h z cs h z c h z t θ=+∆()()(),...ln 2N E R q cs cs h z C ci q F F z d t k z d ϕπ⎛⎫=+⎪⎝⎭()0.00547 1.83213.8c c k t =⨯⨯++()()()()12c cs ci z t z t z t =+()()(),,..2N E R q u h z ci h z ci u gq F F Z t t d d h ϕπ=++()()()()()0,,0..4h z u h z N E t t R q u u q F F Z k t dt k t dt ϕπ=+⎰⎰(),DNB h z q3.7DNBR 的计算3.8计算热管中的压降 3.9单相流体的摩擦压降式中：单相流体加速压降：单相流体提升压降：局部压降，出口：进口：定位格架出口压降：其中，比容v 按相应的流体压力和温度，由表查得。

()()()(),,...DNB h z DNB h z N E R q q q DNBR q z q F F z ϕ==22..2..2f e e L V p f D L G v f D ρ∆==0.640.250.316..Re nw w iso f f f f μμμμ⎛⎫⎛⎫⎛⎫== ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭ (2)in outel p g L g Lρρρ+∆==2222out out outout outoutV G v p k k ρ∆==2222in in inin ininV G v p k k ρ∆==()221222out in gr gr gr grG v v V p k k ρ+∆==四、程序设计框图五、代码说明书本代码主要由五个小部分组成。

堆芯出口温度计算、堆热流量计算、堆平均参数计算、第一至第六控制体各量计算、热管的压降计算。

（1）堆芯出口温度计算：此段根据任务书给出的基本参数和热量与流量之间关系，运用迭代的算法，求出堆芯的出口温度。

（2）堆热流量计算：先根据堆芯的输出功率和释热率以及总的传热面积，求出燃料元件表面平均热流量，再根据热管因子求出最大热流量。

再求出平均线功率和最大线功率。

（3）堆平均参数计算：根据基本的尺寸，求出堆体的流通截面积和一个栅元的流通截面积。

然后再求出流经栅元的流量。

依据上面的温度结果，查出热物性参数，再求出冷却剂的流速。

（4）第一至第五控制体的各量计算：因为五个控制体的计算过程类似，这里只说明第一个控制体的计算过程。

在现有的参数下，根据热流量与流量的关系和迭代算法，求出该控制体的出口温度。

通过流通截面积与湿周的关系求出栅元的当量直径。

再根据上面的温度，查出对应的热物性参数由雷诺数与努尔数的关系，解出控制体出口处的对流换热系数。

因为不知该处的流体状态，分别用单相强迫对流放热公式和詹斯-洛特斯传热方程算出各自的膜温压，取较小的值加上出口处的流体温度即是包壳的外表面温度。

由包壳的外表面的温度再根据圆管的传热方程运用迭代算法解出包壳内表面的温度。

芯块与包壳内表面之间的导热问题，根据间隙导热模型，即可解出芯块表面的温度，根据内热源的导热模型，依据积分热导率与温度的对应关系和插值方法，解出芯块中心的温度。

接下来依据冷却剂的温度，得出的控制体出口处的含汽量。

进而依据W-3公式求出该出的临界热流量，最后得出该出的烧毁比DNBR。