于洋 劉東 宋小明 王雅峰 董豎彪

摘 要本文基于節(jié)點(diǎn)壓力法，建立了核動(dòng)力系統(tǒng)流體網(wǎng)絡(luò)數(shù)學(xué)模型，開發(fā)了適用于可壓縮流體（以水或蒸汽為工質(zhì)）的單相熱工流體網(wǎng)絡(luò)計(jì)算軟件。仿真結(jié)果表明，該軟件可應(yīng)用于復(fù)雜管網(wǎng)系統(tǒng)的模擬，具有良好的穩(wěn)定性。

關(guān)鍵詞節(jié)點(diǎn)壓力法;單相;熱工流體網(wǎng)絡(luò)

0 引言

在核動(dòng)力系統(tǒng)中，流體網(wǎng)絡(luò)（簡(jiǎn)稱“流網(wǎng)”）仿真技術(shù)由于具有計(jì)算速度快、穩(wěn)定性高的特點(diǎn)，可用于輔助系統(tǒng)和二回路系統(tǒng)的研發(fā)設(shè)計(jì)、運(yùn)行研究。流網(wǎng)仿真計(jì)算一般采用節(jié)點(diǎn)壓力法，將流網(wǎng)系統(tǒng)模型簡(jiǎn)化為常微分方程組和代數(shù)方程組進(jìn)行求解。

許多熱工流網(wǎng)系統(tǒng)的建模與仿真可簡(jiǎn)化為單相流網(wǎng)系統(tǒng)的求解問題。相對(duì)于多相流網(wǎng)系統(tǒng)、單相流網(wǎng)系統(tǒng)數(shù)學(xué)模型較簡(jiǎn)單，求解速度更快，更容易滿足實(shí)時(shí)性要求。本文針對(duì)核動(dòng)力系統(tǒng)開發(fā)了適用于可壓縮流體（以水或蒸汽為工質(zhì)）的單相流網(wǎng)計(jì)算軟件。

1 單相流網(wǎng)計(jì)算模型

實(shí)際的核動(dòng)力系統(tǒng)由閥門、泵、汽輪機(jī)、換熱器等多種設(shè)備及連接管道所組成。為便于計(jì)算方便，基于集總參數(shù)法將熱工流體系統(tǒng)劃分為節(jié)點(diǎn)和流線，然后針對(duì)節(jié)點(diǎn)和流線分別建立數(shù)學(xué)模型，最終轉(zhuǎn)化為求解節(jié)點(diǎn)壓力和流線流量的問題，如圖1所示。在建模時(shí)做如下假設(shè)[1]：

（1）節(jié)點(diǎn)具有一定容積，且處于均勻狀態(tài)，與外界進(jìn)行質(zhì)量和能量交換;

（2）流線具有一定長(zhǎng)度和流通面積，且有連接方向，但與外界不發(fā)生質(zhì)量和能量交換。

1.1 節(jié)點(diǎn)模型

1.1.1 質(zhì)量守恒方程

流體在流動(dòng)過程中，其流動(dòng)方向可能發(fā)生改變，本文采用帶方向的關(guān)聯(lián)矩陣記錄流網(wǎng)的流動(dòng)方向。矩陣D、DE、DF分別是節(jié)點(diǎn)之間、節(jié)點(diǎn)與壓力邊界之間以及節(jié)點(diǎn)與流量邊界之間的關(guān)聯(lián)矩陣，其中D為Y階方陣，DE為Y行N列的矩陣，而DF則為Y行M列的矩陣，矩陣元素用于存儲(chǔ)流網(wǎng)的流動(dòng)方向，含義如下：

此外，流量值具有方向，其絕對(duì)值代表流量大小，而符號(hào)則表示流動(dòng)方向，當(dāng)實(shí)際流動(dòng)方向與關(guān)聯(lián)矩陣中的流動(dòng)方向一致時(shí)，流量值為正，否則為負(fù)。

1.1.2 能量守恒方程

1.2 流線模型

根據(jù)連接關(guān)系，可將流線分為內(nèi)部流線、壓力邊界流線和流量邊界流線。

1.3 數(shù)值算法

本文采用壓力-焓值分離算法進(jìn)行求解，忽略方程（6）、（7）的非穩(wěn)態(tài)項(xiàng)，并將其按流量半隱式差分，得到壓差和流量的線性關(guān)系式：

2 軟件總體設(shè)計(jì)與計(jì)算流程

2.1 軟件設(shè)計(jì)

軟件模塊關(guān)系如圖1所示。在數(shù)據(jù)輸入輸出過程中，存儲(chǔ)模塊作為中轉(zhuǎn)站，為讀寫模塊和計(jì)算模塊提供高效的數(shù)據(jù)寄存平臺(tái);在軟件運(yùn)行控制過程中，控制模塊統(tǒng)籌各模塊的運(yùn)行;在計(jì)算過程中，流網(wǎng)計(jì)算模塊通過調(diào)用水與蒸汽物性庫、矩陣解法庫、設(shè)備模塊仿真計(jì)算。

2.2 軟件計(jì)算流程

單相流網(wǎng)計(jì)算軟件的計(jì)算流程見圖2，其計(jì)算過程為：

（1）從指定路徑讀取輸入卡，給定邊界條件、初始條件和結(jié)構(gòu)參數(shù);

（2）對(duì)流網(wǎng)進(jìn)行拓?fù)渥R(shí)別，并將流網(wǎng)中不存在質(zhì)量交換的小流網(wǎng)劃分為若干子流網(wǎng);

（3）采用壓力焓值分離式算法計(jì)算流網(wǎng)的熱工水力參數(shù)，并根據(jù)熱板蓄熱量計(jì)算熱板溫度;

（4）計(jì)算收斂判定;完成后，將計(jì)算結(jié)果以csv文件形式保存至指定路徑。

3 仿真實(shí)例驗(yàn)證

以華龍一號(hào)中壓安注和低壓安注系統(tǒng)為對(duì)象，針對(duì)中壓安注系統(tǒng)運(yùn)行、低壓安注系統(tǒng)處于備用狀態(tài)瞬態(tài)工況進(jìn)行了仿真用于驗(yàn)證軟件模擬復(fù)雜管網(wǎng)（止回閥+并聯(lián)管網(wǎng)+回流管網(wǎng)）的能力，流程圖如圖3所示。

假設(shè)在事故工況下，當(dāng)一回路系統(tǒng)壓力為3.5MPa時(shí)，RCS（主冷卻劑系統(tǒng)）壓力高于低壓安注泵注入壓頭而低于中壓安注泵注入壓頭時(shí)，中壓安注系統(tǒng)的小流量管線自動(dòng)隔離，熱段注入管線閥門關(guān)閉，中壓安注泵從IRWST取水后從冷段注入，而低壓安注系統(tǒng)則通過小流量管線維持泵的運(yùn)行，此時(shí)低壓安注系統(tǒng)處于備用狀態(tài)。在軟件計(jì)算穩(wěn)定后，使RCS壓力躍遷變化為3MPa，理論上此時(shí)低壓安注系統(tǒng)仍處于備用狀態(tài)，而中壓安注系統(tǒng)的安注流量增加。設(shè)計(jì)資料顯示，系統(tǒng)設(shè)計(jì)流通形式如圖4所示。

仿真結(jié)果顯示，RCS系統(tǒng)壓力降低后，低壓泵從IRWST取水，然后全部經(jīng)旁通管線回流至IRWST，即低壓安注系統(tǒng)處于備用狀態(tài)。同時(shí)，中壓安注泵的取水流量增加，中壓安注泵出口壓力下降，且中壓安注并聯(lián)注入管線的壓力和流量的瞬態(tài)變化趨勢(shì)一致且穩(wěn)態(tài)值相同，這與設(shè)計(jì)流通形式是一致的。此外，曲線變化趨勢(shì)平穩(wěn)，表明軟件在仿真復(fù)雜流網(wǎng)系統(tǒng)具有良好的穩(wěn)定性，結(jié)果如圖5所示。

4 結(jié)論

本文開發(fā)了適用于可壓縮流體的單相熱工流體網(wǎng)絡(luò)計(jì)算軟件。針對(duì)華龍一號(hào)中壓安注系統(tǒng)運(yùn)行、低壓安注系統(tǒng)處于備用狀態(tài)瞬態(tài)工況進(jìn)行了仿真驗(yàn)證，軟件具備模擬復(fù)雜管網(wǎng)系統(tǒng)的能力。

參考文獻(xiàn)

[1]蔡瑞忠，謝茂清.熱工流體網(wǎng)絡(luò)的實(shí)時(shí)仿真模型及其算法[J].系統(tǒng)仿真學(xué)報(bào)，1992，4（4）：13-18，26.

[2]Sylvie Lorente， Wishsanuruk Wechsatol， Adrian Bejan. Fundamentals of tree-shaped networks of insulated pipes for hot water and exergy[J]. Exergy， an International Journal， 2002， 2： 227-236.

[3]周盛煌.兩相流體網(wǎng)絡(luò)水力熱力耦合模型及計(jì)算方法研究[D]. 哈爾濱工程大學(xué)，2014.