CN105225710A - 循环回路系统的堆芯模拟体和蒸汽发生器模拟体设计方法 - Google Patents

Abstract

本发明公开了一种循环回路系统的堆芯模拟体和蒸汽发生器模拟体设计方法，包括堆芯模拟体的设计要求、结构设计、相似比例分析和设计，以及蒸汽发生器模拟体的设计要求、结构设计和设计结果，该方法根据自然循环实验回路的设计参数和功能需求，确定实验回路系统堆芯模拟体和蒸汽发生器模拟体的设计要求，再分别对堆芯模拟体和蒸汽发生器模拟体进行设计和分析计算。本发明提高了超设计基准事故设置的措施的可靠性和核电厂安全系统的多样化，具有很高的工程价值。

Description

循环回路系统的堆芯模拟体和蒸汽发生器模拟体设计方法

技术领域

本发明属于反应堆热工水利领域，尤其涉及一种循环回路系统的堆芯模拟体和蒸汽发生器模拟体设计方法。

背景技术

反应堆可以利用自然循环不依赖外部动力就将热量导出的能力，实现反应堆非能动安全设施在事故下的运行，从而提高反应堆的安全性。

自日本福岛事件后，国际和国内社会对核能安全提出了更高的要求，特别针对全厂断电和完全丧失冷却链等超设计基准事故缓解措施的可靠性给予了越来越多的关注。2012年6月国家核安全局对外发布的《福岛核事故后核电厂改进行动通用技术要求(试行)》中，多次提出对在运和在建核电机组在核电厂部分或全部安全系统功能丧失的情况下，如超设计基准洪水事件条件下，应该采取更多措施以带出余热。新近发布的《“十二五”期间新建核电厂安全要求和审评原则》中，明确要求十二五期间新建核电厂须增加反应堆堆芯的余热排出、应急冷却和最终热阱的考虑，应设置多样化的最终热阱。

ACPR1000项目目前的技术方案，尽管已经针对超设计基准事故设置了一些缓解措施，但在核电厂安全系统多样化设计方面尚有较大欠缺。根据目前确定论安全分析和PSA分析结果，蒸汽发生器二次侧相关事故具有重大贡献。因此，需要针对完全丧失给水、主蒸汽管道破裂及主给水管道破裂叠加辅助给水丧失、全厂断电、完全丧失冷却链等超设计基准事故设置可靠性更高的多样化缓解系统，而二次侧非能动余热排除系统正符合该要求。本发明正是针对上述问题而设计。

发明内容

本发明实施例的目的在于提供一种循环回路系统的堆芯模拟体和蒸汽发生器模拟体设计方法，旨在解决对超设计基准事故设置的措施可靠性不够高，核电厂安全系统多样化存在较大欠缺的问题。

本发明是这样实现的，一种循环回路系统的堆芯模拟体和蒸汽发生器模拟体设计方法，包括堆芯模拟体的设计要求、结构设计、相似比例分析和设计，以及蒸汽发生器模拟体的设计要求、结构设计和设计结果，该方法根据自然循环实验回路的设计参数和功能需求，确定实验回路系统堆芯模拟体和蒸汽发生器模拟体的设计要求，再分别对堆芯模拟体和蒸汽发生器模拟体进行设计和分析计算。

进一步，堆芯模拟体的设计以ACPR1000为原型，ACPR1000堆芯有157盒燃料组件，该堆芯模拟体模拟其中的一盒燃料组件，高度比取1:4，速度比为1:2，电加热元件的表面热流密度比为2:1；

堆芯模拟体包括铜棒、平盖、筒体端部、进出口接头、进出口法兰、筒体、吊篮、陶瓷、燃料组件模拟件、下封头、挂环、进出口反法兰；

燃料组件模拟件通过吊篮悬挂于挂环上，平盖上装有2根铜棒，铜棒中的两根铜棒一端分别和电源的正负极相连，另一端与燃料组件模拟件相连；

组件模拟件和金属吊篮之间用陶瓷填充。

进一步，堆芯模拟体堆芯模拟体管道用长度、水力学当量直径和流通截面积体现其流速、流动阻力。

进一步，堆芯模拟体设计中的无量纲动量方程为：

${\Π}_R{\mathrm{\Δ\ρ}}^{+}{l}_{hc}^{+}=\underset{i}{\Σ}\frac{l_i^{+}}{A_i^{+}}\frac{\∂u_{in}^{+}}{\∂\mathrm{\τ}}+\frac{1}{2}{\left(u_{in}^{+}\right)}^2{\Π}_F+{\left(u_{in}^{+}\right)}^2\underset{i}{\Σ}{\left(\frac{A_{in}}{A_i}\right)}^2\left(\frac{{\mathrm{\ρ}}_{in}^{+}}{{\mathrm{\ρ}}_i^{+}}\right)$

无量纲能量方程为：

$\left(\frac{{\mathrm{\ρ}}^{+}}{{\mathrm{\ρ}}_{in}^{+}}\right)\frac{\∂x}{\∂\mathrm{\τ}}+\frac{\∂x}{\∂s}={\Π}_{pch}{q}^{+}-{\Π}_C\frac{\∂\mathrm{\θ}}{\∂\mathrm{\τ}}$

动量和能量方程中无量纲参数的定义分别为：Richardson数表征浮升力与惯性力之间的关系；阻力系数表征回路的摩擦和形阻；热源数也称相变数表征燃料元件释热对回路焓升的影响；热容数表征结构的热容量与回路焓升的关系。

进一步，堆芯模拟体设计中的无量纲相变数为：

${\Π}_{pch}=\frac{q_0{l}_{hc0}{\mathrm{\ξ}}_i}{{\mathrm{\ρ}}_0{u}_0{A}_0{\mathrm{\Δh}}_{fg}}\left(\frac{{\mathrm{\Δ\ρ}}_{fg}}{{\mathrm{\ρ}}_g}\right)$

无量纲过冷度数为：

${\Π}_{sub}=\frac{{\mathrm{\Δh}}_{sub}}{h_{fg}}\left(\frac{{\mathrm{\Δ\ρ}}_{fg}}{{\mathrm{\ρ}}_g}\right).$

进一步，蒸汽发生器模拟体对蒸汽发生器进行模化，换热管高度比取1:4，采用原型的换热管径以及栅距，蒸汽发生器一次侧和二次侧的温度和压力都与原型一致，根据原型中蒸汽发生器的换热能力和换热面积，模化得到实验回路中蒸汽发生器的换热面积及换热管的结构尺寸；

蒸汽发生器模拟体包括上封头、端壳体、换热管、下端壳体、下封头，换热管通过固定管板固定于下端壳体上，换热管的不同位置装有足够数量的支撑板和挡板，下封头被隔板隔开隔开。

进一步，堆芯模拟体采用长1050mm，外径9.5mm的不锈钢管作为加热元件，电加热元件采用定位格架固定，定位格架拟采用原型中的定位格架，电加热元件的上下端和铜板相连，铜板分别和电源的正负极相连，堆芯模拟体的进出口管道用内径为105mm的不锈钢管，位置按原型布置，电加热元件和压力容器之间用陶瓷填充。

本发明提高了超设计基准事故设置的措施的可靠性和核电厂安全系统的多样化，具有很高的工程价值。

附图说明

图1是本发明实施例提供的堆芯模拟体主视图；

图2是本发明实施例提供的堆芯模拟体进出口管府视图；

图3是本发明实施例提供的堆芯模拟体管座俯视图；

图4是本发明实施例提供的换热管布管示意图；

图5是本发明实施例提供的换热管布管的B-B向剖视图；

图中：1、铜棒；2、平盖；3、筒体端部；4、进出口接头；5、进出口法兰；6、筒体；7、吊篮；8、陶瓷；9、燃料组件模拟件；10、下封头；11、挂环；12、进出口反法兰；13、蒸汽出口接管；14、上封头；15、上端壳体；16、换热管；17、下端壳体；18、固定管板；19、下封头；20、隔板；21、去离子水进出口接管；22、自来水进口接管；23、去离子水进口；24、去离子水出口；25、自来水进口；26、蒸汽出口。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本本发明，并不用于限定本发明。

1、堆芯模拟体的设计要求

堆芯模拟体是自然循环实验回路开展非能动余热排除实验研究的关键性设备之一。以ACPR1000为原型，ACPR1000堆芯有157盒燃料组件，本堆芯模拟体模拟其中的一盒燃料组件。高度比取1:4，速度比为1:2，电加热元件的表面热流密度比为2:1，并合理设计轴向和横向阻力特性。

2、堆芯模拟体结构设计

堆芯模拟体结构总图如图1-图3所示。堆芯模拟体包括铜棒1、平盖2、筒体端部3、进出口接头4、进出口法兰5、筒体6、吊篮7、陶瓷8、燃料组件模拟件9、下封头10、挂环11、进出口反法兰12。

所述燃料组件模拟件9采用长1050mm(燃料组件高度的1/4)，外径9.5mm(燃料组件的外径)的不锈钢管作为加热元件。为模拟燃料组件的定位格架，电加热元件也采用定位格架固定，其定位格架拟采用原型中的定位格架。燃料组件模拟件9通过吊篮7悬挂于挂环11上，实现固定。平盖2上装有2根铜棒1，铜棒1中的两根铜棒一端分别和电源的正负极相连，另一端与燃料组件模拟件9相连，实现电加热的目的。

为减小金属热容，组件模拟件9和金属吊篮7之间用陶瓷8填充。进出口接头4采用用内径为105mm的不锈钢管(约原型尺寸的1/7)，且位置按原型布置。

3、堆芯模拟体相似比例分析

反应堆一回路系统结构复杂，不论是堆芯还是蒸汽发生器，都需要通过一定的假设进行简化。反应堆一回路系统的设备可以简单的分为两种类型，即有热量传递和无热量传递；三种设备：管路、蒸汽发生器、堆芯。其中，此处的蒸汽发生器是指发生了进出口温度或密度变化的那部分，即蒸汽发生器管束。而反应堆堆芯则是发生了温度或密度变化的那部分堆芯，即加热部分，其余部分则视为管道。如果考虑到两相流型的影响，严格的说管道内流体的密度也可以发生变化。管道的特征是没有热量传递，只是流体的位置发生变化，并贡献阻力。反应堆管束中，如冷管和热管，在符合一维假设方面没有什么大的问题，而一些局部有突变的构件可以通过阻力作用来体现其在自然循环过程中的作用。因此，管道可以用长度、水力学当量直径和流通截面积体现其流速、流动阻力；通过在重力方向的投影体现其位置变化对重力作用的贡献。不论是几根管道，均可以一维化为上述特征的元件。

3.1基本方程

首先是一维假设，认为流动沿着设备或管道的轴向方向进行，而传热则在一维边界的径向进行；其次，可以通过采用特征参数，如热工水力直径等方法，对径向影响的作用进行评估并进行模拟；此外，还需要忽略热损失、轴向传热及方程中的高阶项，其控制方程为：

(1)连续性方程

$\frac{\∂\mathrm{\ρ}}{\∂t}+\frac{\∂\mathrm{\ρ}u}{\∂s}=0---\left(1\right)$

上式中ρ为流体密度，u为流体流速。

(2)动量方程

$-\frac{\∂P}{\∂s}=\frac{\∂\mathrm{\ρ}u}{\∂t}+\frac{\∂{\mathrm{\ρu}}^2}{\∂s}+\mathrm{\ρ}g-\frac{f}{De}\frac{1}{2}{\mathrm{\ρu}}^2---\left(2\right)$

上式中P为压力，f为阻力系数，D_e为水力学当量直径。

(3)能量方程

$\frac{\∂\mathrm{\ρ}h}{\∂t}+\frac{\∂\mathrm{\ρ}uh}{\∂s}=\frac{\mathrm{\ζ}}{A}h(T_w-T_f)---\left(3\right)$

上式中h为流体的焓，ζ为周长，A为截面积，T_w和T_f分别为壁面温度和流体温度。

(4)传热方程

$\frac{\∂\mathrm{\ρ}cT}{\∂t}+\frac{\mathrm{\ξ}}{A}h(T_w-T_f)=\left(\frac{\mathrm{\ξ}}{A}\right)q_s---\left(4\right)$

本发明对流换热系数不进行模拟，主要原因有：①金属构件结构复杂，本模型中的一维方法模拟不能准确描述；②两相传热温差很小，测量到的壁温不具有很大的价值；③燃料元件基本上是定热流密度传热，不受传热系数的影响，而燃料棒本身的壁温与材料、结构等有关，不必详细模拟，对其研究可以在其他情况下进行。这样，对基本方程中的能量方程和传热方程进行合并，其新的能量方程为：

$\frac{\∂\mathrm{\ρ}h}{\∂t}+\frac{\∂\mathrm{\ρ}uh}{\∂s}=\left(\frac{\mathrm{\ξ}}{A}\right)q_s-\frac{\∂\mathrm{\ρ}cT}{\∂t}---\left(5\right)$

将能量方程用质量含汽量表示

${\mathrm{\ρh}}_{fg}\frac{\∂x}{\∂t}+{\mathrm{\ρuh}}_{fg}\frac{\∂x}{\∂s}=\left(\frac{\mathrm{\ξ}}{A}\right)q_s-\mathrm{\ρ}c\frac{\∂T}{\∂t}---\left(6\right)$

3.2方程无量纲化

由能量方程(6)可知，回路流体的密度，主要受边界传热的控制。当堆芯传热有变化时，此时局部的密度发生变化，但该密度变化并不会立即对流速产生很大影响。只有该密度变化因流体流动传递到整个回路从而实质性的影响回路中的密度差，从而导致流速变化。因此，可以认为回路流体密度对于时间的偏微分项是一个可以忽略的项。因此，连续性方程变化为：

ρuA＝ρ_inu_inA_in＝const(7)

上式中ρ_in、u_in和A_in分别为堆芯入口处的密度、流速和流通截面积。

将守恒方程无量纲化，其无量纲参数分别为：速度长度密度 ${\mathrm{\ρ}}^{+}=\frac{\mathrm{\ρ}}{{\mathrm{\ρ}}_0};$ 密度差 ${\mathrm{\Δ\ρ}}^{+}=\frac{\mathrm{\Δ}\mathrm{\ρ}}{{\mathrm{\Δ\ρ}}_0};$ 流通面积 $A^{+}=\frac{A}{A_0};$ 时间 $\mathrm{\τ}=\frac{{\mathrm{tu}}_0}{l_{hc0}},$ 温度 $\mathrm{\θ}=\frac{T-T_0}{{\mathrm{\ΔT}}_0}$ 等，其角标0代表参数的初始值。

连续方程无量纲化：

$\frac{\mathrm{\ρ}uA}{{\mathrm{\ρ}}_0{u}_0{A}_0}=\frac{{\mathrm{\ρ}}_{in}{u}_{in}{A}_{in}}{{\mathrm{\ρ}}_0{u}_0{A}_0}---\left(8\right)$

即无量纲速度为 $u^{+}=\left(\frac{A_{in}^{+}}{A^{+}}\right)\left(\frac{{\mathrm{\ρ}}_{in}^{+}}{{\mathrm{\ρ}}^{+}}\right)u_{in}^{+}---\left(9\right)$

动量方程无量纲化：

(1)重力压降项Δρgl_hc

$\frac{{\mathrm{\Δ\ρgl}}_{hc}}{{\mathrm{\ρ}}_{in}{u}_0^2}=\frac{{\mathrm{\Δ\ρ}}_0{\mathrm{gl}}_{hc0}}{{\mathrm{\ρ}}_{in}{u}_0^2}\frac{\mathrm{\Δ}\mathrm{\ρ}}{{\mathrm{\Δ\ρ}}_0}\frac{l_{hc}}{l_{hc0}}={\Π}_R{\mathrm{\Δ\ρ}}^{+}{l}_{hc}^{+}---\left(10\right)$

(2)惯性项 $\underset{i}{\Σ}\frac{\∂\mathrm{\ρ}u}{\∂t}{l}_i=\underset{i}{\Σ}\frac{A_{in}}{A}\frac{\∂{\mathrm{\ρ}}_{in}{u}_{in}}{\∂t}{l}_i=\underset{i}{\Σ}\frac{A_{in}}{A_{in}}{\mathrm{\ρ}}_{in}\frac{\∂u_{in}}{\∂t}{l}_i$

由 $\mathrm{\τ}=\frac{{\mathrm{tu}}_0}{l_{hc0}}$ 可得 $t=\frac{{\mathrm{\τl}}_{hc0}}{u_0}$

$\begin{array}{c}\frac{\underset{i}{\Σ}\frac{\∂\mathrm{\ρ}u}{\∂t}{l}_i}{{\mathrm{\ρ}}_{in}{u}_0^2}=\frac{\underset{i}{\Σ}\frac{A_{in}}{A_i}{\mathrm{\ρ}}_{in}\frac{\∂u_{in}}{\∂t}{l}_i}{{\mathrm{\ρ}}_{in}{u}_0^2}=\frac{\underset{i}{\Σ}\frac{A_{in}}{A_i}{\mathrm{\ρ}}_{in}\frac{\∂u_{in}}{\∂\frac{{\mathrm{\τl}}_{hc0}}{u_0}}{l}_i}{{\mathrm{\ρ}}_{in}{u}_0^2}\\ =\underset{i}{\Σ}\frac{A_{in}}{A_i}\frac{\∂u_{in}^{+}}{\∂\mathrm{\τ}}{l}_i^{+}=\underset{i}{\Σ}\frac{l_i^{+}}{A_i^{+}}\frac{\∂u_{in}^{+}}{\∂\mathrm{\τ}}\end{array}---\left(11\right)$

(3)摩擦项 $\underset{i}{\Σ}(\frac{f}{De}{l}_i+K_i)\frac{1}{2}{\mathrm{\ρu}}^2$

$\begin{array}{c}\frac{\underset{i}{\Σ}(\frac{f}{De}{l}_i+K_i)\frac{1}{2}{\mathrm{\ρu}}^2}{{\mathrm{\ρ}}_{in}{u}_0^2}=\frac{1}{2}\frac{\underset{i}{\Σ}(\frac{1}{De}{l}_i+K_i)}{{\mathrm{\ρ}}_{in}{u}_0^2}\frac{A_{in}^2{\mathrm{\ρ}}_{in}^2{u}_{in}^2}{A^2\mathrm{\ρ}}\\ =\frac{1}{2}{\left(u_{in}^{+}\right)}^2\underset{i}{\Σ}{\left(\frac{A_{in}}{A_i}\right)}^2\left(\frac{{\mathrm{\ρ}}_{in}^{+}}{{\mathrm{\ρ}}_i^{+}}\right)(\frac{f}{De}{l}_i+K_i)=\frac{1}{2}{\left(u_{in}^{+}\right)}^2{\Π}_F\end{array}---\left(12\right)$

(4)迁移项 $\underset{i}{\Σ}\frac{\∂{\mathrm{\ρu}}^2}{\∂s}{l}_i$

$\frac{\underset{i}{\Σ}{\mathrm{\ρ}}_i{u}_i^2}{{\mathrm{\ρ}}_{in}{u}_0^2}=\frac{\underset{i}{\Σ}\frac{A_{in}^2{\mathrm{\ρ}}_{in}^2{u}_{in}^2}{A_i^2{\mathrm{\ρ}}_i}}{{\mathrm{\ρ}}_{in}{u}_0^2}={\left(u_{in}^{+}\right)}^2\underset{i}{\Σ}{\left(\frac{A_{in}}{A_i}\right)}^2\left(\frac{{\mathrm{\ρ}}_{in}^{+}}{{\mathrm{\ρ}}_i^{+}}\right)---\left(13\right)$

无量纲化后的动量方程为

${\Π}_R{\mathrm{\Δ\ρ}}^{+}{l}_{hc}^{+}=\underset{i}{\Σ}\frac{l_i^{+}}{A_i^{+}}\frac{\∂u_{in}^{+}}{\∂\mathrm{\τ}}+\frac{1}{2}{\left(u_{in}^{+}\right)}^2{\Π}_F+{\left(u_{in}^{+}\right)}^2\underset{i}{\Σ}{\left(\frac{A_{in}}{A_i}\right)}^2\left(\frac{{\mathrm{\ρ}}_{in}^{+}}{{\mathrm{\ρ}}_i^{+}}\right)---\left(14\right)$

无量纲化后的能量方程为(推导方法与动量方程类似)

$\left(\frac{{\mathrm{\ρ}}^{+}}{{\mathrm{\ρ}}_{in}^{+}}\right)\frac{\∂x}{\∂\mathrm{\τ}}+\frac{\∂x}{\∂s}={\Π}_{pch}{q}^{+}-{\Π}_C\frac{\∂\mathrm{\θ}}{\∂\mathrm{\τ}}---\left(15\right)$

动量和能量方程中无量纲参数的定义分别为：Richardson数表征浮升力与惯性力之间的关系；阻力系数表征回路的摩擦和形阻；热源数也称相变数表征燃料元件释热对回路焓升的影响；热容数表征结构的热容量与回路焓升的关系。

3.3相似准则

首先确定高度比相似

${\left(l_{hc}^{+}\right)}_R=1---\left(16\right)$

为满足原型和模型之间相似，(14)和(15)式的模型和原型的比值均等于1，所以有

${({\Π}_R=\frac{{\mathrm{\Δ\ρ}}_0{\mathrm{gl}}_{hc0}}{{\mathrm{\ρ}}_{in}{u}_0^2})}_R=1---\left(17\right)$

${({\Π}_F=\underset{i}{\Σ}{\left(\frac{A_{in}}{A_i}\right)}^2(\frac{{\mathrm{\ρ}}_{in}^{+}}{{\mathrm{\ρ}}_i^{+}}\left)\right(\frac{f}{De}{l}_i+K_i\left)\right)}_R=1---\left(18\right)$

${({\Π}_H=\frac{q_{0}lhc0\mathrm{\ξ}}{{\mathrm{\ρ}}_0{u}_0{A}_0{h}_{fg}}={\Π}_{pch})}_R=1---\left(19\right)$

${({\Π}_C=\frac{{\mathrm{\ρc}}_v{\mathrm{\ΔT}}_0}{{\mathrm{\ρ}}_0{h}_{fg}})}_R=1---\left(20\right)$

为满足连续方程相似，有

${\left(\frac{A_{in}^{+}}{A^{+}}\right)}_R{\left(\frac{{\mathrm{\ρ}}_{in}^{+}}{{\mathrm{\ρ}}^{+}}\right)}_R=1---\left(21\right)$

对于流动阻力，主要是因为气液两相流动模拟特性的需要，采用原型相同的棒束排列结构较为适宜。在相同介质(水-水模拟)且等压模拟中，首先，原型和模型的流通面积相似，即

$A_R=\frac{{(A/A_0)}_m}{{(A/A_0)}_p}=1---\left(22\right)$

由 ${\left(\frac{A_{in}^{+}}{A^{+}}\right)}_R{\left(\frac{{\mathrm{\ρ}}_{in}^{+}}{{\mathrm{\ρ}}^{+}}\right)}_R=1$ 及 $A_R=\frac{{(A/A_0)}_m}{{(A/A_0)}_p}=1$ 可得 ${\left(\frac{{\mathrm{\ρ}}_{in}^{+}}{{\mathrm{\ρ}}_i^{+}}\right)}_R=1---\left(23\right)$

根据式(23)有

${\left(\frac{{\mathrm{\ρ}}_{in}^{+}}{{\mathrm{\ρ}}_i^{+}}\right)}_R=1---\left(24\right)$

即这个准则数表征原型和模型的压力比为定值：

$\frac{P_m}{P_p}\≡Cons\tan t---\left(25\right)$

根据上式可以得到堆芯进出口密度存在以下关系

${\left(\frac{{\mathrm{\Δ\ρ}}_0}{{\mathrm{\ρ}}_0}\right)}_R=1---\left(26\right)$

根据式(17)、式(24)和式(26)可以得到

${\left({\Π}_R\right)}_R={\left(\frac{{\mathrm{\Δ\ρ}}_0{\mathrm{gl}}_{ch0}}{{\mathrm{\ρ}}_{in}{u}_0^2}\right)}_R={\left(\frac{{\mathrm{\Δ\ρ}}_{0}g}{{\mathrm{\ρ}}_{in}}\right)}_R{\left(\frac{l_{hc0}}{u_0^2}\right)}_R={\left(\frac{l_{hc0}}{u_0^2}\right)}_R=1---\left(27\right)$

上式表达的是相同物性下高度比和自然循环速度比的关系。

(Π_Fi)_R据式(21)可以简化为：

${\left({\Π}_{Fi}\right)}_R={\[\underset{i}{\Σ}\left(\frac{A_{in}}{A_i}\right)(\frac{fl}{D_e}+K)\]}_R=1---\left(28\right)$

由于模型和原型的压力比为定值，所以成立。

同样，可以将[(Π_H)_i]_R写成：

${\[{\left({\Π}_H\right)}_i\]}_R={\[\frac{q_0{l}_{hc0}{\mathrm{\ξ}}_i}{{\mathrm{\ρ}}_0{u}_0{A}_0{h}_{fg}}\]}_R={\[\frac{q_0{l}_{hc0}{\mathrm{\ξ}}_i}{{\mathrm{\ρ}}_0{u}_0{A}_0{h}_{fg}}\left(\frac{{\mathrm{\Δ\ρ}}_{fg}}{{\mathrm{\ρ}}_g}\right)\]}_R---\left(29\right)$

该式就是无量纲相变数同理，由入口含气量的定义可以推导到无量纲过冷度数

从能量方程得知，堆芯内的质量含气量变化相同，如果保证堆芯出口的含气量相等(X_e)_R＝1，则(X_in)_R＝1，有

(Π_sub)_R＝1(30)

对于热容数有：

${\[{\left({\Π}_c\right)}_i\]}_R={\[{\left(\frac{{\mathrm{\ρc}}_v{\mathrm{\ΔT}}_{0}A}{{\mathrm{\ρ}}_0{A}_i{h}_{fg}}\right)}_i\]}_R={\[{\left(\frac{{\mathrm{\ρc}}_v{\mathrm{\ΔT}}_0}{{\mathrm{\ρ}}_0{\mathrm{\Δh}}_{fg}}\right)}_i\]}_R{\[\left(\frac{A_0}{A_i}\right)\]}_R---\left(31\right)$

上式在流体和金属材料的热物性相同时，有：

${\[{\left({\Π}_c\right)}_i\]}_R={\[{\left(\frac{A_s}{A_i}\right)}_i\]}_R=1---\left(32\right)$

该式表示在每一局部区域，金属构件与回路流体之间的质量关系。

综上所述，等压水-水模拟准则关系如表1所示：

表1等压水-水模拟准则数

3.4堆芯模拟体

本实验拟模拟一盒燃料组件，且原型中有157盒燃料组件，原型中堆芯入口流通截面积为由(22)式可得

${\left(\frac{1}{A_0=\frac{1}{4}{\mathrm{\τD}}^2}\right)}_m={\left(\frac{157}{A_0=\frac{1}{4}\mathrm{\π}{760}^2\×3}\right)}_p---\left(32\right)$

即模型中堆芯入口管径约为105mm。本实验中，直径比取1:7.23，面积比1:52，高度比取1:4，即速度比为1:2，电加热元件的表面热流密度比为2:1。堆芯模拟体的设计结果为：

a)堆芯高度1:4设计；

b)流通面积按照1:52设计，堆芯流速比为1:2；

c)采用直接加热方式，棒的表面热流密度比为2:1；

d)采用原型的棒直径，以及栅距；

e)合理设计轴向和横向阻力特性(模拟定位格架)；

f)适当模拟横向流动特性。

4、蒸汽发生器设计要求

与堆芯模拟体类似，对蒸汽发生器进行模化，换热管高度比取1:4，采用原型的换热管径以及栅距。由模化可知，蒸汽发生器的热工参数尽应和原型保持一致，即蒸汽发生器一次侧和二次侧的温度和压力都应与原型一致。根据原型中蒸汽发生器的换热能力和换热面积，模化就可得到实验回路中蒸汽发生器的换热面积及换热管的结构尺寸。

5、蒸汽发生器结构设计

如图4和5所示为蒸汽发生器结构设计图。蒸汽发生器包括上封头14、端壳体15、换热管16、下端壳体17、下封头19。去离子水由去离子水进出口23进入换热管16后与冷凝体换热，并由去离子水出口24流出，所述换热管16通过固定管板18固定于下端壳体17上，换热管16的不同位置应装有足够数量的支撑板和挡板，下封头19被隔板20隔开隔开。冷流体自来水由自来水进口25流入，流入后沉积在下端壳体17的空腔内与换热管16内的热流体发生热交换后气化为蒸汽，并由蒸汽出口26排出进行循环。

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (7)

1.一种循环回路系统的堆芯模拟体和蒸汽发生器模拟体设计方法，其特征在于，所述的循环回路系统的堆芯模拟体和蒸汽发生器模拟体设计方法包括堆芯模拟体的设计要求、结构设计、相似比例分析和设计，以及蒸汽发生器模拟体的设计要求、结构设计和设计结果，该方法根据自然循环实验回路的设计参数和功能需求，确定实验回路系统堆芯模拟体和蒸汽发生器模拟体的设计要求，再分别对堆芯模拟体和蒸汽发生器模拟体进行设计和分析计算。

2.如权利要求1所述的循环回路系统的堆芯模拟体和蒸汽发生器模拟体设计方法，其特征在于，堆芯模拟体的设计以ACPR1000为原型，ACPR1000堆芯有157盒燃料组件，该堆芯模拟体模拟其中的一盒燃料组件，高度比取1:4，速度比为1:2，电加热元件的表面热流密度比为2:1；

堆芯模拟体包括铜棒、平盖、筒体端部、进出口接头、进出口法兰、筒体、吊篮、陶瓷、燃料组件模拟件、下封头、挂环、进出口反法兰；

燃料组件模拟件通过吊篮悬挂于挂环上，平盖上装有2根铜棒，铜棒中的两根铜棒一端分别和电源的正负极相连，另一端与燃料组件模拟件相连；

组件模拟件和金属吊篮之间用陶瓷填充。

3.如权利要求1所述的循环回路系统的堆芯模拟体和蒸汽发生器模拟体设计方法，其特征在于，堆芯模拟体堆芯模拟体管道用长度、水力学当量直径和流通截面积体现其流速、流动阻力。

4.如权利要求1所述的循环回路系统的堆芯模拟体和蒸汽发生器模拟体设计方法，其特征在于，堆芯模拟体设计中的无量纲动量方程为：

${\Π}_R{\mathrm{\Δ\ρ}}^{+}{l}_{hc}^{+}=\underset{i}{\Σ}\frac{l_i^{+}}{A_i^{+}}\frac{\∂u_{in}^{+}}{\∂\mathrm{\τ}}+\frac{1}{2}{\left(u_{in}^{+}\right)}^2{\Π}_F+{\left(u_{in}^{+}\right)}^2\underset{i}{\Σ}{\left(\frac{A_{in}}{A_i}\right)}^2\left(\frac{{\mathrm{\ρ}}_{in}^{+}}{{\mathrm{\ρ}}_i^{+}}\right)$

无量纲能量方程为：

$\left(\frac{{\mathrm{\ρ}}^{+}}{{\mathrm{\ρ}}_{in}^{+}}\right)\frac{\∂x}{\∂\mathrm{\τ}}+\frac{\∂x}{\∂s}={\Π}_{pch}{q}^{+}-{\Π}_C\frac{\∂\mathrm{\θ}}{\∂\mathrm{\τ}}$

动量和能量方程中无量纲参数的定义分别为：Richardson数表征浮升力与惯性力之间的关系；阻力系数表征回路的摩擦和形阻；热源数也称相变数表征燃料元件释热对回路焓升的影响；热容数表征结构的热容量与回路焓升的关系。

5.如权利要求1所述的循环回路系统的堆芯模拟体和蒸汽发生器模拟体设计方法，其特征在于，堆芯模拟体设计中的无量纲相变数为：

${\Π}_{pch}=\frac{q_0{l}_{hc0}{\mathrm{\ξ}}_i}{{\mathrm{\ρ}}_0{u}_0{A}_0{\mathrm{\Δh}}_{fg}}\left(\frac{{\mathrm{\Δ\ρ}}_{fg}}{{\mathrm{\ρ}}_g}\right)$

无量纲过冷度数为：

${\Π}_{sub}=\frac{{\mathrm{\Δh}}_{sub}}{h_{fg}}\left(\frac{{\mathrm{\Δ\ρ}}_{fg}}{{\mathrm{\ρ}}_g}\right).$

6.如权利要求1所述的循环回路系统的堆芯模拟体和蒸汽发生器模拟体设计方法，其特征在于，蒸汽发生器模拟体对蒸汽发生器进行模化，换热管高度比取1:4，采用原型的换热管径以及栅距，蒸汽发生器一次侧和二次侧的温度和压力都与原型一致，根据原型中蒸汽发生器的换热能力和换热面积，模化得到实验回路中蒸汽发生器的换热面积及换热管的结构尺寸；

蒸汽发生器模拟体包括上封头、端壳体、换热管、下端壳体、下封头，换热管通过固定管板固定于下端壳体上，换热管的不同位置装有足够数量的支撑板和挡板，下封头被隔板隔开隔开。

7.如权利要求1所述的循环回路系统的堆芯模拟体和蒸汽发生器模拟体设计方法，其特征在于，堆芯模拟体采用长1050mm，外径9.5mm的不锈钢管作为加热元件，电加热元件采用定位格架固定，定位格架拟采用原型中的定位格架，电加热元件的上下端和铜板相连，铜板分别和电源的正负极相连，堆芯模拟体的进出口管道用内径为105mm的不锈钢管，位置按原型布置，电加热元件和压力容器之间用陶瓷填充。