2. 上海交通大学 核科学与工程学院,上海 200240
2. Shanghai Jiao Tong University, Shanghai 200240, China
螺旋管式直流蒸汽发生器由于体积小、换热效率高等特性[1-2],是一体化反应堆重要设备之一。为研究一体化反应堆的仿真系统及控制方法,需要建立螺旋管式直流蒸汽发生器的动态模型[1-8]。另外,由于螺旋管式直流蒸汽发生器二回路水容积小,蒸汽流量变化过程中蒸汽压力极易发生变化,且若给水跟不上,很容易对二回路设备造成影响[9],因此需要对螺旋管式直流蒸汽发生器开展控制方案研究,以保证二次侧蒸汽出口压力恒定。
本文将螺旋管式直流蒸汽发生器分为一次侧、金属管壁、二次侧3个部分,每个部分基于移动边界理论,将整个轴向分为过冷段、两相段、过热段3个部分,结合每段的质量、能量、动量三大守恒方程,给出了螺旋管式直流蒸汽发生器的状态方程。以国际革新与安全反应堆(IRIS)为参考对象,进行稳态计算和动态仿真实验。当蒸汽流量变化时,蒸汽发生器蒸汽出口压力发生变化,此时蒸汽发生器蒸汽出口压力的偏差作为PID控制器的输入,进而改变给水流量,从而使蒸汽出口压力回到给定值范围内。
1 物理模型根据螺旋管内流体的流动换热过程,对螺旋管式直流蒸汽发生器模型做如下假设:1)由于螺旋管直流蒸汽发生器结构复杂,在建立热工水力数学模型时可忽略对热工特性影响较小的结构,主要分析一次侧冷却剂与螺旋管内二回路工质在螺旋管区域热交换过程,未考虑蒸汽发生器进出口直管段换热以及其他热损失;2)在建立模型时,假设一次侧工质不发生相变,为单相液流动;3)将螺旋管做单管模型处理,等效方法如下:总质量流量为单根螺旋管的质量流量与螺旋管数目的积,总流通面积为单根螺旋管的流通面积与螺旋管数目的积,总换热面积为单根螺旋管的换热面积与螺旋管数目的积;4)使用一维模型,忽略轴向导热;5)一次侧密度沿着轴向维度不发生变化;6)一次侧不考虑轴向压降,且压力沿着时间维度保持不变;7)忽略金属管壁密度的变化;8)一次侧不考虑控制体外部压力做功。
基于以上假设,根据二次侧工质的状态,将蒸汽发生器分为一次侧、金属管壁、二次侧3部分。基于可移动边界理论,将二次侧沿轴向管长分为3个区域,即过冷段ssc、两相段ssb和过热段sss;同理将一次侧分为
Download:
|
|
在图1中,二次侧过冷段为1~3区域、两相段为3~5区域、过热段为5~7区域,其中1、3、5、7为各段边界,2、4、6为各段控制体中心代表点,3为饱和液相点,5为饱和汽相点。
2 数学模型 2.1 基本方程螺旋管式直流蒸汽发生器是由很多根螺旋管并联组成的表面换热设备。二次侧管内工质在流动时,不断地吸热,导致热力状态发生变化,基本的质量、能量和动量三大守恒方程为
$\frac{{\partial \rho }}{{\partial t}} + \frac{{\partial D}}{{\partial l}} = 0$ | (1) |
$\frac{{\partial \left( {\rho h} \right)}}{{\partial t}} = \frac{{\partial \left( {Dh} \right)}}{{\partial l}} + \frac{Q}{A} + \frac{{\partial P}}{{\partial t}}$ | (2) |
$\frac{{\partial D}}{{\partial t}} + \frac{\partial }{{\partial l}}\left( {\frac{{{D^2}}}{\rho }} \right) + \frac{{\partial P}}{{\partial t}} + \frac{{f{D^2}S}}{{2\rho A}} + \rho g = 0$ | (3) |
式中:ρ为流体密度;D为单位面积质量流速;h为流体比焓;Q为单位长度热流率;A为流通面积;P为流体压力;f为摩擦系数;g为重力加速度。
2.2 可移动边界理论设控制体流通面为单位面积A,j、j+1分别为上、下边界的节点,上界面为Zj+1,下界面为Zj,
Download:
|
|
由图2可知,控制体质量变化为
$\frac{{{\rm{d}}\left( {\overline \rho L} \right)}}{{{\rm{d}}t}} = {W_j} - {W_{j + 1}} + {\rho _{j + 1}}\frac{{{\rm{d}}{Z_{j + 1}}}}{{{\rm{d}}t}} - {\rho _j}\frac{{{\rm{d}}{Z_j}}}{{{\rm{d}}t}} + S$ | (4) |
又由于控制体对外膨胀做功为
$\begin{split}& \dfrac{{{\rm{d}}\left( {\overline \rho L\overline h } \right)}}{{{\rm{d}}t}} - \dfrac{{{\rm{d}}\left( {\overline P L} \right)}}{{{\rm{d}}t}} = {W_j}{h_j} - {W_{j + 1}}{h_{j + 1}} + {P_j}\dfrac{{{\rm{d}}{Z_j}}}{{{\rm{d}}t}} - \\ & \quad\quad {P_{j + 1}}\dfrac{{{\rm{d}}{Z_{j + 1}}}}{{{\rm{d}}t}} + {\rho _{j + 1}}{h_{j + 1}}\dfrac{{{\rm{d}}{Z_{j + 1}}}}{{{\rm{d}}t}} - {\rho _j}{h_j}\dfrac{{{\rm{d}}{Z_j}}}{{{\rm{d}}t}} + Q \end{split} $ | (5) |
一次侧通过使用二次侧定义的移动边界,将一次侧分成3个区域,忽略每段进出口流量变化,假设在瞬态过程中,一次侧工质水为单相不可压缩流体,根据移动边界差分理论,得到能量守恒方程。设下标i、o、p和m分别代表进口、出口、一次侧和金属管壁;M为流体质量;
由于
$ \begin{split}&{C}_{{p}_{\text{Ⅰ}}}\dfrac{{\rm{d}}\left({M}_{{p}_{\text{Ⅰ}}}{T}_{{p}_{\text{Ⅰ}}}\right)}{{\rm{d}}t}-{A}_{\rm{p}}\dfrac{{\rm{d}}\left({P}_{\rm{p}}{L}_{{s}_{\rm{ss}}}\right)}{{\rm{d}}t}={W}_{\rm{p}}{C}_{{p}_{\text{Ⅰ}}}\left({T}_{{p}_{\text{Ⅰ}},{\rm{i}}}-{T}_{{p}_{\text{Ⅰ}},{\rm{o}}}\right)-\\ &\quad\quad\quad\quad\quad {K}_{{p}_{\text{Ⅰ}}{m}_{\text{Ⅰ}}}\left({T}_{{p}_{\text{Ⅰ}}}-{T}_{{m}_{\text{Ⅰ}}}\right){\rm{{\text{π}} }}{D}_{\rm{o}}{L}_{{s}_{\rm{ss}}}N-\\ &\quad\quad\quad\quad\quad {\rho }_{{p}_{\text{Ⅰ}}}{A}_{\rm{p}}{C}_{{p}_{\text{Ⅰ}}}{T}_{{p}_{\text{Ⅰ}},{\rm{o}}}\left(\dfrac{{\rm{d}}{L}_{{s}_{\rm{sb}}}}{{\rm{d}}t}+\dfrac{{\rm{d}}{L}_{{s}_{\rm{sc}}}}{{\rm{d}}t}\right)\\[-15pt] \end{split}$ | (6) |
同理可得
$\begin{split}& {C}_{{p}_{\text{Ⅱ}}}\dfrac{{\rm{d}}\left({M}_{{p}_{\text{Ⅱ}}}{T}_{{p}_{\text{Ⅱ}}}\right)}{{\rm{d}}t}-{A}_{\rm{p}}\dfrac{{\rm{d}}\left({P}_{\rm{p}}{L}_{{s}_{\rm{sb}}}\right)}{{\rm{d}}t}={W}_{\rm{p}}{C}_{{p}_{\text{Ⅱ}}}({T}_{{p}_{\text{Ⅰ}}}-\\ &\quad\quad {T}_{{p}_{\text{Ⅱ}}})-{K}_{{p}_{\text{Ⅱ}}{m}_{\text{Ⅱ}}}\left({T}_{{p}_{\text{Ⅱ}}}-{T}_{{m}_{\text{Ⅱ}}}\right){\rm{{\text{π}} }}{D}_{\rm{o}}{L}_{{s}_{\rm{sb}}}N+\\ &\quad\quad\quad {\rho }_{{p}_{\text{Ⅰ}}}{A}_{\rm{p}}{C}_{{p}_{\text{Ⅰ}}}{T}_{{p}_{\text{Ⅰ}},{\rm{o}}}\left(\dfrac{{\rm{d}}{L}_{{s}_{\rm{sb}}}}{{\rm{d}}t}+\dfrac{{\rm{d}}{L}_{{s}_{\rm{sc}}}}{{\rm{d}}t}\right)-\\ &\quad\quad\quad\quad\quad\quad {\rho }_{{p}_{\text{Ⅱ}}}{A}_{\rm{p}}{C}_{{p}_{\text{Ⅱ}}}{T}_{{p}_{\text{Ⅱ}},{\rm{o}}}\dfrac{{\rm{d}}{L}_{{s}_{\rm{sc}}}}{{\rm{d}}t}\end{split}$ | (7) |
$ \begin{split}&{C}_{{p}_{\text{Ⅲ}}}\dfrac{{\rm{d}}\left({M}_{{p}_{\text{Ⅲ}}}{T}_{{p}_{\text{Ⅲ}}}\right)}{{\rm{d}}t}-{A}_{\rm{p}}\dfrac{{\rm{d}}\left({P}_{\rm{p}}{L}_{{s}_{\rm{sc}}}\right)}{{\rm{d}}t}={W}_{\rm{p}}{C}_{{p}_{\text{Ⅲ}}}\left({T}_{{p}_{\text{Ⅱ}},{\rm{o}}}-{T}_{{p}_{\text{Ⅲ}},{\rm{o}}}\right)-\\ &{K}_{{p}_{\text{Ⅲ}}{m}_{\text{Ⅲ}}}\left({T}_{{p}_{\text{Ⅲ}}}-{T}_{{m}_{\text{Ⅲ}}}\right)\left({\rm{{\text{π}} }}{D}_{\rm{{\rm{o}}}}{L}_{{s}_{\rm{sc}}}N\right)+{\rho }_{{p}_{\text{Ⅱ}}}{A}_{\rm{p}}{C}_{{p}_{\text{Ⅱ}}}{T}_{{p}_{\text{Ⅱ}},{\rm{o}}}\dfrac{{\rm{d}}{L}_{{s}_{\rm{sc}}}}{{\rm{d}}t}\\[-15pt]\end{split}$ | (8) |
金属管壁侧分为3个节点,利用可移动边界法,在推导出每个金属管壁节点能量方程时,必须考虑由于边界移动引起的能量传递。图3为金属管壁的节点能量传递图。
Download:
|
|
设Am为金属管壁传热面积;Cm为金属管壁比热容;
$\begin{split}& {C}_{\rm{m}}\dfrac{{\rm{d}}\left({M}_{{m}_{\text{Ⅰ}}}{T}_{{m}_{\text{Ⅰ}}}\right)}{{\rm{d}}t}={K}_{{p}_{\text{Ⅰ}}{m}_{\text{Ⅰ}}}\left({T}_{{p}_{\text{Ⅰ}}}-{T}_{{m}_{\text{Ⅰ}}}\right){\rm{{\text{π}} }}{D}_{\rm{o}}{L}_{{s}_{\rm{ss}}}N-\\ & \quad\quad\quad\quad {K}_{{m}_{\text{Ⅰ}}{s}_{\rm{{ss}}}}\left({T}_{{m}_{\text{Ⅰ}}}-{T}_{{s}_{\rm{ss}}}\right){\rm{{\text{π}} }}{D}_{\rm{i}}{L}_{{s}_{\rm{ss}}}N-\\ &\quad {\rho }_{\rm{m}}{A}_{\rm{m}}{C}_{\rm{m}}\left(\dfrac{{T}_{{m}_{\text{Ⅰ}}}+{T}_{{m}_{\text{Ⅱ}}}}{\rm{2}}\right)\left(\dfrac{{\rm{d}}{L}_{{s}_{\rm{sb}}}}{{\rm{d}}t}+\dfrac{{\rm{d}}{L}_{{s}_{\rm{sc}}}}{{\rm{d}}t}\right)\end{split}$ | (9) |
同理,金属管壁
$\begin{split} &{C}_{\rm{m}}\dfrac{{\rm{d}}\left({M}_{{m}_{\text{Ⅱ}}}{T}_{{m}_{\text{Ⅱ}}}\right)}{{\rm{d}}t}={K}_{{p}_{\text{Ⅱ}}{m}_{\text{Ⅱ}}}\left({T}_{{p}_{\text{Ⅱ}}}-{T}_{{m}_{\text{Ⅱ}}}\right){\rm{{\text{π}} }}{D}_{\rm{o}}{L}_{{s}_{\rm{sb}}}N-\\ &\quad\quad\quad\quad {K}_{{m}_{\text{Ⅱ}}{s}_{{\rm{sb}}}}\left({T}_{{m}_{\text{Ⅱ}}}-{T}_{{s}_{\rm{sb}}}\right){\rm{{\text{π}} }}{D}_{\rm{i}}{L}_{{s}_{\rm{sb}}}N-\\ &\quad\quad\quad\quad {\rho }_{\rm{m}}{A}_{\rm{m}}{C}_{\rm{m}}\left(\dfrac{{T}_{{m}_{{\rm I}}}+{T}_{{m}_{\text{Ⅱ}}}}{2}\right)\dfrac{{\rm{d}}{L}_{{s}_{\rm{sc}}}}{{\rm{d}}t}+\\ &\quad +{\rho }_{\rm{m}}{A}_{\rm{m}}{C}_{\rm{m}}\left(\dfrac{{T}_{{m}_{{\rm I}}}+{T}_{{m}_{\text{Ⅱ}}}}{2}\right)\left(\dfrac{{\rm{d}}{L}_{{s}_{\rm{sb}}}}{{\rm{d}}t}+\dfrac{{\rm{d}}{L}_{{s}_{\rm{sc}}}}{{\rm{d}}t}\right)\end{split}$ | (10) |
$\begin{split}& {C}_{\rm{m}}\dfrac{{\rm{d}}\left({M}_{{m}_{\text{Ⅲ}}}{T}_{{m}_{\text{Ⅲ}}}\right)}{{\rm{d}}t}={K}_{{p}_{\text{Ⅲ}}{m}_{\text{Ⅲ}}}\left({T}_{{p}_{\text{Ⅲ}}}-{T}_{{m}_{\text{Ⅲ}}}\right){\rm{{\text{π}} }}{D}_{\rm{o}}{L}_{{s}_{\rm{sc}}}N-\\ &\quad\quad\quad\quad\quad {K}_{{m}_{\text{Ⅲ}}{s}_{{\rm{sc}}}}\left({T}_{{m}_{\text{Ⅲ}}}-{T}_{{s}_{\rm{sc}}}\right){\rm{{\text{π}} }}{D}_{\rm{i}}{L}_{{s}_{\rm{sc}}}N+\\ &\quad\quad\quad\quad\quad {\rho }_{\rm{m}}{A}_{\rm{m}}{C}_{\rm{m}}\left(\dfrac{{T}_{{m}_{\text{Ⅱ}}}+{T}_{{m}_{\text{Ⅲ}}}}{2}\right)\dfrac{{\rm{d}}{L}_{{s}_{\rm{sc}}}}{{\rm{d}}t}\end{split}$ | (11) |
根据螺旋管式直流蒸汽发生器的传热分区,将式(1)~(3)在各传热段内进行积分,得质量方程和能量方程。其中,控制体的温度使用控制体平均温度进行计算,控制体的压力和流量使用控制体平均压力和平均流量进行计算。过冷段进口密度和饱和温度都是压力函数。设下标f、g、s及fw分别代表液相、气相、饱和及给水;As为二次侧流通面积;Zsv为二次侧蒸汽阀门开度;Csv为主蒸汽阀门流量系数;F1, F2,…,F6为二次侧各节点摩擦系数[11];P1, P2,…,P7为二次侧各节点压力;W1, W2,…,W7为二次侧各节点质量流量,且
1)过冷段
由于过冷段质量的变化为
${A_{\rm{s}}}\frac{{{\rm{d}}\left( {{\rho _{{s_{{\rm{sc}}}}}}{L_{{s_{{\rm{sc}}}}}}} \right)}}{{{\rm{d}}t}} = {W_{{\rm{fw}}}} - {W_{\rm{3}}} + {\rho _{\rm{f}}}{A_{\rm{s}}}\frac{{{\rm{d}}{L_{{s_{{\rm{sc}}}}}}}}{{{\rm{d}}t}}$ | (12) |
由于过冷段能量变化为
$\begin{split} &{C}_{{s}_{\rm{sc}}}{A}_{\rm{s}}\dfrac{{\rm{d}}\left(\dfrac{{T}_{\rm{fw}}+{T}_{\rm{f}}}{\rm{2}}{\rho }_{{s}_{\rm{sc}}}{L}_{{s}_{\rm{sc}}}\right)}{{\rm{d}}t}-{A}_{\rm{s}}\dfrac{{\rm{d}}\left({P}_{\rm{2}}{L}_{{s}_{\rm{sc}}}\right)}{{\rm{d}}t}={W}_{\rm{fw}}{h}_{\rm{fw}}-\\ &\quad\quad\quad {W}_{\rm{3}}{h}_{\rm{sf}}+{K}_{{m}_{\text{Ⅲ}}{s}_{\rm{sc}}}\left({T}_{{m}_{\text{Ⅲ}}}-{T}_{{s}_{\rm{sc}}}\right){\rm{{\text{π}} }}{D}_{\rm{i}}{L}_{{s}_{\rm{sc}}}N+\\ &\quad\quad\quad\quad\quad {\rho }_{\rm{f}}{C}_{\rm{f}}{A}_{\rm{s}}{T}_{\rm{f}}\dfrac{{\rm{d}}{L}_{{s}_{\rm{sc}}}}{{\rm{d}}t}-{P}_{\rm{3}}{A}_{\rm{s}}\dfrac{{\rm{d}}{L}_{{s}_{\rm{sc}}}}{{\rm{d}}t}\end{split}$ | (13) |
由于在动量方程中,惯性项比摩擦项小很多[7-8],故简化的动量守恒方程为
${P_{\rm{3}}} = {P_{\rm{1}}} - {\rho _{\rm{2}}}g{L_{{s_{{\rm{sc}}}}}}\sin {\rm{ }}\theta - {F_{\rm{2}}}{L_{{s_{{\rm{sc}}}}}}D_{\rm{2}}^{\rm{2}}$ | (14) |
两边同时对时间t求导数,写成微分格式有
$ \dfrac{{{\rm{d}}{P_{\rm{3}}}}}{{{\rm{d}}t}} = \dfrac{{{\rm{d}}{P_{\rm{1}}}}}{{{\rm{d}}t}} - \left( {{\rho _{\rm{2}}}g\sin {\rm{ }}\theta + {F_{\rm{2}}}D_{\rm{2}}^2} \right)\dfrac{{{\rm{d}}{L_{{s_{{\rm{sc}}}}}}}}{{{\rm{d}}t}} - {\rm{2}}{F_{\rm{2}}}{L_{{s_{{\rm{sc}}}}}}{D_{\rm{2}}}\dfrac{{{\rm{d}}{D_{\rm{2}}}}}{{{\rm{d}}t}} $ | (15) |
2)两相段
同理可得两相段的质量、能量和动量方程为
$ \dfrac{{{\rm{d}}{M_{{s_{{\rm{sb}}}}}}}}{{{\rm{d}}t}} = {W_{{s_{{\rm{sc}}}},o}} - {W_{{s_{{\rm{sb}}}},o}} - {\rho _{\rm{f}}}{A_{\rm{s}}}\dfrac{{{\rm{d}}{L_{{s_{{\rm{sc}}}}}}}}{{{\rm{d}}t}}+ {\rho _{\rm{g}}}{A_{\rm{s}}}\dfrac{{{\rm{d}}\left( {{L_{{s_{{\rm{sc}}}}}} + {L_{{s_{{\rm{sb}}}}}}} \right)}}{{{\rm{d}}t}} \\ $ | (16) |
$\begin{split} &\dfrac{{\rm{d}}\left({M}_{{s}_{\rm{sb}}}{h}_{\rm{4}}-{P}_{\rm{4}}{A}_{\rm{s}}{L}_{{s}_{\rm{sb}}}\right)}{{\rm{d}}t}={W}_{{s}_{\rm{sc}},o}{h}_{\rm{f}}-{W}_{{s}_{\rm{sb}},o}{h}_{\rm{g}}+\\ &{K}_{{m}_{\text{Ⅱ}}{s}_{\rm{sb}}}\left({T}_{{m}_{\text{Ⅱ}}}-{T}_{{s}_{\rm{sb}}}\right){\text{π}} {D}_{\rm{i}}{L}_{{s}_{\rm{sb}}}N-{\rho }_{\rm{f}}{A}_{\rm{s}}{h}_{\rm{f}}\dfrac{{\rm{d}}{L}_{{s}_{\rm{sc}}}}{{\rm{d}}t}+\\ &\quad{\rho }_{\rm{g}}{A}_{\rm{s}}{h}_{\rm{g}}\dfrac{{\rm{d}}\left({L}_{{s}_{\rm{sc}}}+{L}_{{s}_{\rm{sb}}}\right)}{{\rm{d}}t}+{P}_{\rm{3}}{A}_{\rm{s}}\dfrac{{\rm{d}}{L}_{{s}_{\rm{sc}}}}{{\rm{d}}t}-\\ &\quad\quad\quad\quad\quad{P}_{\rm{5}}{A}_{\rm{s}}\dfrac{{\rm{d}}\left({L}_{{s}_{\rm{sc}}}+{L}_{{s}_{\rm{sb}}}\right)}{{\rm{d}}t}\end{split}$ | (17) |
$ \dfrac{{{\rm{d}}{P_{\rm{5}}}}}{{{\rm{d}}t}} = \dfrac{{{\rm{d}}{P_{\rm{3}}}}}{{{\rm{d}}t}} - \left( {{\rho _{\rm{4}}}g\sin {\rm{ }}\theta + {F_{\rm{4}}}D_{\rm{4}}^{\rm{2}}} \right)\dfrac{{{\rm{d}}{L_{{s_{{\rm{sb}}}}}}}}{{{\rm{d}}t}} - {\rm{2}}{F_{\rm{4}}}{L_{{s_{{\rm{sb}}}}}}{D_{\rm{4}}}\dfrac{{{\rm{d}}{D_{\rm{4}}}}}{{{\rm{d}}t}} $ | (18) |
对密度、压力和焓做如下处理:
$\begin{aligned} & {\rho _{\rm{4}}} = {X_{{s_{{\rm{sb}}}}}}{\rho _{\rm{5}}} + \left( {{\rm{1}} - {X_{{s_{{\rm{sb}}}}}}} \right){\rho _{\rm{3}}} \\ &{h_{{s_{{\rm{sb}}}}}} = {X_{{s_{{\rm{sb}}}}}}{h_{\rm{g}}} + \left( {{\rm{1}} - {X_{{s_{{\rm{sb}}}}}}} \right){h_{\rm{f}}} \end{aligned} $ | (19) |
饱和水和饱和蒸汽的密度、饱和水和饱和蒸汽的焓都可以写成压力的函数。这里h4对应
3)过热段
同理可得过热段的质量、能量和动量方程为
$\dfrac{{{\rm{d}}{M_{{s_{{\rm{ss}}}}}}}}{{{\rm{d}}t}} = {W_{{s_{{\rm{sb}}}},{\rm{o}}}} - {W_{{s_{{\rm{ss}}}},{\rm{o}}}} - {\rho _{\rm{5}}}{A_{\rm{s}}}\dfrac{{{\rm{d}}\left( {{L_{{s_{{\rm{sc}}}}}} + {L_{{s_{{\rm{sb}}}}}}} \right)}}{{{\rm{d}}t}}$ | (20) |
$\begin{split}& \dfrac{{\rm{d}}\left({M}_{{s}_{\rm{ss}}}{h}_{\rm{so}}-{P}_{\rm{6}}{A}_{\rm{s}}{L}_{{s}_{\rm{ss}}}\right)}{{\rm{d}}t}={W}_{{s}_{\rm{sb}},\rm{o}}{h}_{\rm{sg}}-{W}_{{s}_{\rm{ss}},o}{h}_{\rm{so}}+\\ &\quad\quad\quad\quad {K}_{{m}_{\text{Ⅰ}}{s}_{\rm{sc}}}\left({T}_{{m}_{\text{Ⅰ}}}-{T}_{{s}_{\rm{ss}}}\right){\rm{{\text{π}} }}{D}_{\rm{i}}{L}_{{s}_{\rm{ss}}}N-\\ &\quad\quad {\rho }_{\rm{5}}{A}_{\rm{s}}{h}_{\rm{5}}\dfrac{{\rm{d}}\left({L}_{{s}_{\rm{sc}}}+{L}_{{s}_{\rm{sb}}}\right)}{{\rm{d}}t}+{P}_{\rm{5}}{A}_{\rm{s}}\dfrac{{\rm{d}}\left({L}_{{s}_{\rm{sc}}}+{L}_{{s}_{\rm{sb}}}\right)}{{\rm{d}}t}\end{split}$ | (21) |
$\; \dfrac{{{\rm{d}}{P_{\rm{7}}}}}{{{\rm{d}}t}} = \dfrac{{{\rm{d}}{P_{\rm{5}}}}}{{{\rm{d}}t}} - \left( {{\rho _{\rm{6}}}g\sin \left( \theta \right) + {F_{\rm{6}}}D_{\rm{6}}^2} \right)\dfrac{{{\rm{d}}{L_{{s_{{\rm{ss}}}}}}}}{{{\rm{d}}t}} - {\rm{2}}{F_{\rm{6}}}{L_{{s_{{\rm{ss}}}}}}{D_6}\dfrac{{{\rm{d}}{D_{\rm{6}}}}}{{{\rm{d}}t}} $ | (22) |
过热段蒸汽出口密度和焓、压力均相关。蒸汽出口流量与蒸汽压力关系为
${W_{{s_{{\rm{ss}}}},o}} = {C_{{\rm{sv}}}}{Z_{{\rm{sv}}}}{P_{\rm{7}}}$ | (23) |
由式(1)~(23),螺旋管式直流蒸汽发生器的输入为一次侧流量Wp、一次侧压力Pp、一次侧入口温度
$M\left( {{x}} \right)\dot {{x}} = F\left( {{{x}},{{u}}} \right)$ | (24) |
式中:
$ \begin{array}{c} {\dot{{{x}}}}^{\rm{T}}=\Bigg[\dfrac{{\rm{d}}{T}_{{p}_{\text{Ⅰ}},{\rm{o}}}}{{\rm{d}}t}\quad \dfrac{{\rm{d}}{T}_{{p}_{\text{Ⅱ}},{\rm{o}}}}{{\rm{d}}t}\quad \dfrac{{\rm{d}}{T}_{{p}_{\text{Ⅲ}},{\rm{o}}}}{{\rm{d}}t}\quad \dfrac{{\rm{d}}{T}_{{m}_{\text{Ⅰ}}}}{{\rm{d}}t}\quad \dfrac{{\rm{d}}{T}_{{m}_{\text{Ⅱ}}}}{{\rm{d}}t}\quad \\ \quad \quad \dfrac{{\rm{d}}{T}_{{m}_{\text{Ⅲ}}}}{{\rm{d}}t} \quad \dfrac{{\rm{d}}{P}_{3}}{{\rm{d}}t}\quad \dfrac{{\rm{d}}{P}_{5}}{{\rm{d}}t} \quad \dfrac{{\rm{d}}{P}_{7}}{{\rm{d}}t}\quad \dfrac{{\rm{d}}{L}_{{s}_{\rm{sc}}}}{{\rm{d}}t}\quad \dfrac{{\rm{d}}{L}_{{s}_{\rm{sb}}}}{{\rm{d}}t}\quad \\ \quad \quad \dfrac{{\rm{d}}{W}_{3}}{{\rm{d}}t}\quad\dfrac{{\rm{d}}{W}_{5}}{{\rm{d}}t}\quad \dfrac{{\rm{d}}{h}_{\rm{so}}}{{\rm{d}}t}\quad \dfrac{{\rm{d}}{P}_{1}}{{\rm{d}}t}\Bigg]\\ {{u}}=\left[{W}_{\rm{p}}\quad {P}_{\rm{p}}\quad {T}_{{p}_{\text{Ⅰ}},{\rm{i}}}\quad {W}_{\rm{fw}}\quad {h}_{\rm{fw}}\quad {Z}_{\rm{sv}}\right] \end{array} $ |
为了检验模型的正确性和适用程度,根据IRIS的设计和运行参数,使用Matlab程序语言SIMULINK仿真系统,构建仿真计算平台,对模型进行稳态结果和动态结果的检验。稳态结果可以和已有的稳态设计数据进行定量比较,文献[1,12]给出了IRIS在100%负荷工况下的稳态设计值,将此作为模型稳态检验的基础,以文献中的结果为参考值,进行模型稳态结果的定量检验。对于形式如式(24)的动态特性方程,因为在稳定状态时
表1列出了稳态程序计算结果与IRIS设计值[12]的对比情况。从表1可知,通过本文模型计算得到的稳态温度、压降等数据与设计值符合较好,相对偏差最大项均出现在二次侧压降上,且计算结果小于1%。仿真结果表明本文建立的模型具有代表性,物理机理关系式选用合理,且求解方法正确有效。
为验证模型的动态特性,本文进行了阶跃仿真实验。
考察第1个典型的动态工况:在100%额定负荷下,保持一次侧流量、一次侧压力、一次侧入口温度、二次侧给水温度和焓值以及蒸汽调节阀开度等输入参数不变,二次侧给水流量阶跃减小5%,系统状态量变化如图4所示。
Download:
|
|
瞬态开始时,二次侧流量阶跃降低,此时一次侧与二次侧的换热并不会瞬间降低,即二次侧的吸热能力并不会瞬间降低。根据能量守恒定理Q=cm∆t可知,m在减小,Q保持不变,故二次侧进出口温差∆t增加,所以二次侧出口温度增加。随着时间的推移,流量降低到一定程度并且稳定下来,流量的降低,使得二次侧换热系数减少,进而使得一、二次侧的换热量减少,一次侧出口温度增加;二次侧给水加热到饱和水和饱和蒸汽所需的热量减少,使得过冷段和两相段的长度减小,过热段长度增加。由于过热段长度的增加使得二次侧出口蒸汽温度升高。由于进出口质量守恒,所以二次侧给水流量和出口蒸汽流量相等,进而由式(23)可知,蒸汽出口压力与蒸汽流量成正比,出口压力会降低。
当系统各参数经过100 s的瞬态后,一次侧和二次侧的能量再次达到平衡,各出口参数重新稳定在一个新的水平上。
考察第2个典型的动态工况:在100%额定负荷下,保持一次侧入口温度、一次侧压力、二次侧给水流量、二次侧给水温度和焓值以及蒸汽调节阀开度等输入参数不变,在10 s时,一次侧入口流量阶跃增加10%,观察系统状态参数的如图5所示变化。
Download:
|
|
瞬态开始时,一次侧流量阶跃增加,但是此时一次侧与二次侧的换热并不会瞬间增加,即二次侧的吸热能力并不会瞬间增加。根据能量守恒定理可知,一次侧流量在增加,热量保持不变,故一次侧进出口温差降低,所以一次侧温差减小,一次侧出口温度增加。随着时间的推移,流量增加到一定程度,并且稳定下来,流量的增加使得一次侧传热系数增加,进而使得一、二次侧的换热量增加,过热段长度增加,过冷段和两相段的长度减小。由于过热段长度的增加使得二次侧出口蒸汽温度升高;由于两相段长度的减少使得二次侧压降降低。
当系统各参数经过100 s的瞬态后,一次侧和二次侧的能量再次达到平衡,各出口参数重新稳定在一个新的水平上。
4 控制系统仿真 4.1 控制器设计本文采用PID控制器建立螺旋管式直流蒸汽发生器蒸汽出口压力控制系统。当改变蒸汽出口流量时,通过调节二次侧给水流量来保证蒸汽出口压力恒定,给水流量调节的输入信号是蒸汽出口压力的偏差。
$W = {K_{\rm{P}}}\left[ {\Delta {P_{\rm{s}}} + \frac{1}{{{T_{\rm{I}}}}}\int_0^t {\Delta {P_{\rm{s}}}{\rm{d}}t + {T_{\rm{D}}}\frac{{{\rm{d}}\Delta {P_{\rm{s}}}}}{{{\rm{d}}t}}} } \right]$ |
式中:KP为比例系数;TI为积分时间常数;TD为微分时间常数[13]。
由于螺旋管式直流蒸汽发生器二回路水容积小,因此蒸汽流量变化过程中,蒸汽压力极易发生变化,且若给水跟不上,很容易对二回路设备造成影响。因此,为快速响应蒸汽流量的变化,本文提出用蒸汽压力偏差信号和蒸汽流量信号对二次侧给水流量进行调节。其中,蒸汽压力偏差信号为主控信号,蒸汽流量作为辅控信号,共同调节给水流量。控制器结构如图6所示。需求给水流量表达式为
${W_{{\rm{fw}}}} = {K_{\rm{1}}}{W_{\rm{s}}} + {K_{\rm{2}}}\Delta {P_{\rm{s}}} + \frac{{{K_{\rm{2}}}}}{{{\tau _{\rm{I}}}}}\int {\Delta {P_{\rm{s}}}{\rm{d}}t} + {K_{\rm{2}}}{\tau _{\rm{D}}}\frac{{{\rm{d}}\Delta {P_{\rm{s}}}}}{{{\rm{d}}t}}$ |
式中:K1为蒸汽流量与需求给水流量的转换系数;K2为比例常数;
Download:
|
|
在100%额定负荷下,保持一次侧入口温度、一次侧压力、一次侧流量、二次侧给水流量以及二次侧给水温度和焓值等输入参数不变,前10 s为满功率稳态运行;在10 s时,给蒸汽调节阀开度斜坡输入,令变化速率为1%/s,使得在90 s时完成把蒸汽调节阀开度线性降低为满功率时的20%。由式(23)可知,蒸汽调节阀开度与蒸汽流量成正比,故有蒸汽流量线性降低为满功率的80%,此时控制仿真结果如图7所示,最终蒸汽出口压力可以稳定到所要的值5.8 MPa。
Download:
|
|
当蒸汽出口流量降低时,蒸汽温度和蒸汽压力会升高。此时,二次侧从一次侧带走的热量减少,一次侧出口温度增加。随着流量的降低,把二次侧给水加热到饱和水和饱和蒸汽所需的热量减小,使得过冷段和两相段的长度减小,过热段长度增加。各出口参数均在250 s内达到稳定值。
5 结论对已建立螺旋管式直流蒸汽发生器设备的数学模型进行了满功率稳态计算,结果表明本文建立的数学模型具有代表性,物理机理关系式选用合理,且求解方法正确有效。
1)对一次侧流量和二次侧给水流量进行阶跃扰动,系统的动态响应过程符合热工水力学基本规律。
2)针对螺旋管式直流蒸汽发生器的结构特点和动态特性,本文提出了基于PID调节给水流量,使得蒸汽出口压力保持恒定的控制方案。大幅度甩负荷的实验结果表明了控制方法的有效性,也验证模型的合理性。
[1] | 袁媛, 彭敏俊, 夏庚磊, 等. 螺旋管式直流蒸汽发生器热工水力分析模型[J]. 原子能科学技术, 2014, 48(S1): 251-256. (0) |
[2] | 杨自强. 小型模块化反应堆螺旋管式直流蒸汽发生器热工水力研究[D]. 重庆: 重庆大学, 2018. (0) |
[3] | ABDALLA M A. A four-region, moving-boundary model of a once-through, helical-coil steam generator[J]. Annals of nuclear energy, 1994, 21(9): 541-562. DOI:10.1016/0306-4549(94)90078-7 (0) |
[4] | DE SÁ MAGALHÃES M A, DE OLIVEIRA LIRA C A B, DA SILVA M A B, et al. Boron transient effects on the behavior of IRIS reactor using dynamic modeling[J]. Progress in nuclear energy, 2013, 69: 44-52. DOI:10.1016/j.pnucene.2013.05.006 (0) |
[5] | ARDA S E, HOLBERT K E. A dynamic model of a passively cooled small modular reactor for controller design purposes[J]. Nuclear engineering and design, 2015, 289: 218-230. DOI:10.1016/j.nucengdes.2015.04.026 (0) |
[6] | 付明玉, 夏国清, 常宗虎, 等. 直流蒸汽发生器静态和动态特性分析[J]. 哈尔滨工程大学学报, 2002, 23(3): 28-32. DOI:10.3969/j.issn.1006-7043.2002.03.008 (0) |
[7] | 李海鹏, 黄晓津, 张良驹. 螺旋管式直流蒸汽发生器的集总参数动态模型[J]. 原子能科学技术, 2008, 42(8): 729-733. (0) |
[8] | LI Haipeng, HUANG Xiaojin, ZHANG Liangju. A lumped parameter dynamic model of the helical coiled once-through steam generator with movable boundaries[J]. Nuclear engineering and design, 2008, 238(7): 1657-1663. DOI:10.1016/j.nucengdes.2008.01.009 (0) |
[9] | 许志斌, 吴婕, 全正庭, 等. 直流式蒸汽发生器压力控制仿真[J]. 热力发电, 2015, 44(12): 98-102. DOI:10.3969/j.issn.1002-3364.2015.12.098 (0) |
[10] | LEE M H. Nonlinear dynamic modeling of a once-through steam generator[D]. Knoxville: University of Tennessee, 1978. (0) |
[11] | 黄晓津. 10 MW高温气冷堆的动态建模及动态特性的仿真研究[D]. 北京: 清华大学, 1998. (0) |
[12] | IRIS Consortium. International Reactor Innovative and Secure[R]. 2002. (0) |
[13] | 刘妍, 夏虹. 基于模糊PID控制器一体化反应堆协调控制技术[J]. 哈尔滨商业大学学报(自然科学版), 2013, 29(5): 619-622, 628. (0) |