陳修高，張 鵬，石 洋，胡服全，王 健

（1.國(guó)家電投集團(tuán)科學(xué)技術(shù)研究院有限公司,北京 102209；2.國(guó)核華清(北京)核電技術(shù)研發(fā)中心有限公司,北京 102209；3.北京合工仿真技術(shù)有限公司,北京 100192）

核電調(diào)節(jié)閥廣泛應(yīng)用于核電站核島及常規(guī)島的各個(gè)系統(tǒng)中。調(diào)節(jié)閥性能成為系統(tǒng)性能的組成部分，其中流量特性是核電調(diào)節(jié)閥最為關(guān)注的問(wèn)題之一。

目前研究核電調(diào)節(jié)閥性能主要以試驗(yàn)與數(shù)值模擬為主。近年來(lái)在核電閥門試驗(yàn)研究方面：于海峰等[1]采用試驗(yàn)驗(yàn)證核一級(jí)電動(dòng)截止閥樣機(jī)的性能，其試驗(yàn)數(shù)據(jù)表明，閥門滿足設(shè)計(jì)要求；符明海等[2]對(duì)ACP1000 核級(jí)主蒸汽安全閥的冷態(tài)和熱態(tài)性能進(jìn)行試驗(yàn)，其結(jié)果表明，熱態(tài)試驗(yàn)比冷態(tài)試驗(yàn)更能客觀反映閥門在實(shí)際工作狀態(tài)下的性能；文獻(xiàn)[3]也表明在條件允許的情況下應(yīng)進(jìn)行熱態(tài)試驗(yàn)，以保證核電站的安全性能；齊曉光等[4]開(kāi)展了核級(jí)電動(dòng)截止閥的流體阻斷試驗(yàn)，驗(yàn)證了開(kāi)展閥門流體阻斷試驗(yàn)的試驗(yàn)方法；穆冠宇[5]對(duì)核電廠多級(jí)籠式調(diào)節(jié)閥的布置位置進(jìn)行試驗(yàn)研究，其結(jié)果表明，為不使管道產(chǎn)生破壞性水錘，閥芯需完全布置于浸沒(méi)水的管道中，布置管段中不能有空氣存在；徐斌福等[6]介紹了核級(jí)閥門試驗(yàn)系統(tǒng)，該系統(tǒng)可滿足不同類型核級(jí)閥門產(chǎn)品的功能性試驗(yàn)。在數(shù)值模擬研究方面：劉雄偉等[7]利用數(shù)值模擬對(duì)AP1000 核電中壓調(diào)節(jié)閥所受氣動(dòng)力矩進(jìn)行數(shù)值模擬研究，其模擬結(jié)果與工程計(jì)算得到的結(jié)果基本吻合；Wang 等[8]利用數(shù)值模擬研究了迷宮調(diào)節(jié)閥的迷宮通道，其結(jié)果表明，迷宮通道可有效地使壓力下降，避免了空化現(xiàn)象的產(chǎn)生；何子昂等[9]通過(guò)數(shù)值模擬分析閥門內(nèi)部組件的沖刷損傷；文獻(xiàn)[10?11]采用數(shù)值模擬對(duì)大氣釋放閥前置隔離閥進(jìn)行開(kāi)啟特性的研究，其結(jié)果可為閥門的設(shè)計(jì)和改進(jìn)提供參考；劉立志等[12]對(duì)核能裝置管路系統(tǒng)截止閥流體激振特性進(jìn)行了研究，其結(jié)果表明直流式截止閥的流體激振特性優(yōu)于角式截止閥；陳志杰等[13]利用動(dòng)網(wǎng)格技術(shù)研究了核級(jí)定壓差止回閥的閥瓣運(yùn)動(dòng)，其結(jié)果與試驗(yàn)結(jié)果吻合程度良好；余航等[14]對(duì)船用核動(dòng)力裝置止回閥的泄漏問(wèn)題進(jìn)行了熱流固耦合研究，計(jì)算結(jié)果可對(duì)閥門優(yōu)化提供幫助；沈國(guó)強(qiáng)等[15]對(duì)DN100 的核級(jí)套筒式調(diào)節(jié)閥流量特性進(jìn)行了數(shù)值模擬，模擬值與試驗(yàn)值有較好的吻合；文獻(xiàn)[16?17]分析閥芯形狀對(duì)調(diào)節(jié)閥流量特性的影響，對(duì)進(jìn)一步研究核電站控制閥的流量控制和閥芯的優(yōu)化設(shè)計(jì)具有一定的參考；文獻(xiàn)[18?19]利用數(shù)值模擬和試驗(yàn)相結(jié)合的方式研究了DN100 調(diào)節(jié)閥局部流動(dòng)參數(shù)，其結(jié)果可用來(lái)開(kāi)發(fā)改進(jìn)調(diào)節(jié)閥閥內(nèi)件設(shè)計(jì)；祁崇可等[20]采用數(shù)值模擬方法研究核級(jí)主給水控制閥流道結(jié)構(gòu)，并通過(guò)流量試驗(yàn)驗(yàn)證了數(shù)值模擬方法的準(zhǔn)確性。

在試驗(yàn)研究方面，國(guó)內(nèi)研究者已經(jīng)取得了一系列成果，但大多都是針對(duì)口徑較小的閥門或冷態(tài)進(jìn)行試驗(yàn)，對(duì)于口徑較大的核電閥門（公稱尺寸>DN500）一般采用數(shù)值模擬的方法研究其內(nèi)流的流動(dòng)特征。目前對(duì)于核電閥門的流場(chǎng)數(shù)值模擬研究基本是采用通用CFD 商業(yè)軟件如Fluent 等，雖然其基本能達(dá)到比較滿意的結(jié)果，但在閥門的動(dòng)態(tài)研究方面，對(duì)于復(fù)雜結(jié)構(gòu)的大口徑調(diào)節(jié)閥，動(dòng)網(wǎng)格技術(shù)實(shí)現(xiàn)起來(lái)較為困難，而將調(diào)節(jié)閥的行程等分開(kāi)度進(jìn)行穩(wěn)態(tài)研究，往往會(huì)忽略流量特性的細(xì)微特征。為此，本文采用泵閥類專用數(shù)值模擬軟件Simerics MP+對(duì)DN550 的核電主給水調(diào)節(jié)閥進(jìn)行流量特性研究，并與目前傳統(tǒng)CFD 軟件（Fluent 軟件）結(jié)果進(jìn)行對(duì)比分析。其結(jié)果表明，Simerics MP+在閥門動(dòng)態(tài)研究方面具有獨(dú)特的優(yōu)勢(shì)，可為核電調(diào)節(jié)閥的動(dòng)態(tài)研究提供參考。

1 數(shù)值計(jì)算理論與方法

為研究調(diào)節(jié)閥的流量特性，首先建立連續(xù)性方程和動(dòng)量守恒方程，在研究流量特性的過(guò)程中不考慮流體介質(zhì)的熱交換，所以不需要求解能量方程。流場(chǎng)數(shù)值計(jì)算的基本控制方程為：

式中：ρ為流體密度；i,j=1,2,3；ui表示與xi相關(guān)聯(lián)的速度分量；Sm為源項(xiàng)，是從分散的二級(jí)相中加入到連續(xù)相的質(zhì)量；p為靜壓；gi和Fi分別是i方向上的重力和外部體積力；τij為應(yīng)力張量，由式（3）給出。

在主給水調(diào)節(jié)閥流量特性的研究中，可以使用許多湍流模型來(lái)模擬流場(chǎng)。其中，標(biāo)準(zhǔn)k-ε模型適用范圍廣，計(jì)算成本低，精度合理，是基于湍流動(dòng)能及其耗散率的輸運(yùn)方程的模型。該模型是從實(shí)驗(yàn)現(xiàn)象中總結(jié)出來(lái)的半經(jīng)驗(yàn)公式，目前已成為計(jì)算流場(chǎng)的主要工具。標(biāo)準(zhǔn)k-ε模型能夠在保證結(jié)果準(zhǔn)確性的基礎(chǔ)上，提高計(jì)算效率[16]。因此，本文采用標(biāo)準(zhǔn)k-ε湍流模型進(jìn)行流場(chǎng)數(shù)值模擬研究。其控制方程為：

式中：Gk表示由于平均速度梯度而產(chǎn)生的湍流動(dòng)能；Gb是由于浮力產(chǎn)生的湍流動(dòng)能；YM表示可壓縮湍流中的脈動(dòng)膨脹對(duì)總耗散率的貢獻(xiàn)，本文流量系數(shù)研究的介質(zhì)為水，作為不可壓縮流體處理，所以可忽略此項(xiàng)的影響。

2 數(shù)值模擬幾何模型

2.1 三維結(jié)構(gòu)模型

圖1 是主給水調(diào)節(jié)閥的三維結(jié)構(gòu)示意圖。圖2是調(diào)節(jié)閥三維流道示意圖。為了使閥內(nèi)流體流動(dòng)充分均勻以及數(shù)值計(jì)算結(jié)果接近實(shí)際情況，根據(jù)標(biāo)準(zhǔn)IEC60534-2-3[21]的要求對(duì)調(diào)節(jié)閥進(jìn)出口兩端加長(zhǎng)管道。

圖1 主給水調(diào)節(jié)閥結(jié)構(gòu)示意圖

圖2 主給水調(diào)節(jié)閥流道示意圖

2.2 網(wǎng)格劃分

閥芯與套筒之間存在微小間隙0.2mm，在Simerics MP+中可以建立結(jié)構(gòu)化間隙網(wǎng)格，并進(jìn)行局部加密處理，如圖3 所示。在通用CFD 計(jì)算中為改善網(wǎng)格質(zhì)量，并提高計(jì)算效率，需要對(duì)其間隙和部分結(jié)構(gòu)進(jìn)行簡(jiǎn)化處理，采用四面體和六面體的混合網(wǎng)格對(duì)計(jì)算區(qū)域離散，如圖4 所示。

圖3 Simerics MP+網(wǎng)格劃分

圖4 通用CFD 網(wǎng)格劃分示意圖

2.3 網(wǎng)格無(wú)關(guān)性檢驗(yàn)

為了獲得較為經(jīng)濟(jì)的網(wǎng)格，數(shù)值計(jì)算前需要進(jìn)行網(wǎng)格無(wú)關(guān)性研究，通過(guò)比較計(jì)算模型的相關(guān)物理量隨網(wǎng)格密度變化的情況來(lái)選擇比較精確并且花費(fèi)時(shí)間較少的網(wǎng)格單元數(shù)進(jìn)行數(shù)值計(jì)算。以閥門全開(kāi)時(shí)為例，邊界條件設(shè)置為進(jìn)口壓力8MPa（絕壓8.1MPa），出口流量4867.1m3/h，計(jì)算出口壓力。表1 是Simerics MP+網(wǎng)格無(wú)關(guān)性檢驗(yàn)結(jié)果。計(jì)算結(jié)果表明，增加網(wǎng)格數(shù)量對(duì)計(jì)算結(jié)果影響已經(jīng)很小，所以選取方案2。通用CFD 方法的網(wǎng)格劃分?jǐn)?shù)量已達(dá)到3×106左右，計(jì)算結(jié)果變化已經(jīng)很小。

表1 網(wǎng)格無(wú)關(guān)性檢驗(yàn)

2.4 計(jì)算方法

在Simerics MP+中設(shè)置閥門剛性勻速位移，從全開(kāi)狀態(tài)到關(guān)閉，用時(shí)15s，位移隨時(shí)間變化如圖5所示。Simerics MP+集成有專用的閥門運(yùn)動(dòng)模板，無(wú)論是閥芯剛性勻速位移，還是閥芯與流體壓力等物理量成函數(shù)關(guān)系的動(dòng)態(tài)運(yùn)動(dòng)，都容易實(shí)現(xiàn)。通用CFD 軟件采用將閥門行程均分為10 個(gè)相對(duì)開(kāi)度的方法，每10%開(kāi)度進(jìn)行一組流量特性數(shù)值試驗(yàn)，分別建立模型劃分網(wǎng)格，并進(jìn)行穩(wěn)態(tài)計(jì)算。

圖5 Simerics MP+調(diào)節(jié)閥閥芯位移

流量系數(shù)Kv反映了調(diào)節(jié)閥的流通能力，是調(diào)節(jié)閥的主要評(píng)價(jià)指標(biāo)之一，無(wú)論是在調(diào)節(jié)閥選型還是實(shí)際應(yīng)用中都具有非常重要的作用。本文利用模擬結(jié)果來(lái)計(jì)算調(diào)節(jié)閥流量系數(shù)，從而進(jìn)行對(duì)比分析。根據(jù)標(biāo)準(zhǔn)IEC60534-2-3，調(diào)節(jié)閥的流量系數(shù)計(jì)算為

式中：利用國(guó)際單位制計(jì)算時(shí)，N1值為0.1；qv是通過(guò)調(diào)節(jié)閥的流量，單位為m3/h；Δp為被測(cè)調(diào)節(jié)閥前后壓差，單位為kPa；ρ1為實(shí)際工況下流體密度；ρ0為常溫條件下的流體密度。

2.5 仿真計(jì)算條件

在Simerics MP+和通用CFD 的計(jì)算中，邊界條件相同，具體如表2 所示。

表2 計(jì)算邊界條件

計(jì)算與試驗(yàn)所用的流體介質(zhì)一致，為去離子常溫水，其物性參數(shù)如表3 所示。

表3 流體介質(zhì)物性參數(shù)

3 結(jié)果與分析

3.1 定量結(jié)果對(duì)比

圖6 是2 種方法計(jì)算結(jié)果與試驗(yàn)值的對(duì)比曲線圖。流量系數(shù)曲線是根據(jù)相應(yīng)開(kāi)度的流量和壓差，由式（6）計(jì)算得到。由圖可知：主給水調(diào)節(jié)閥相對(duì)開(kāi)度小于70%時(shí)，為等百分比流量特性；相對(duì)開(kāi)度大于70%時(shí)，接近線性流量特性。由此可見(jiàn)，此主給水調(diào)節(jié)閥為復(fù)合式流量特性。相比通用CFD每10%開(kāi)度進(jìn)行一組流量特性數(shù)值試驗(yàn)得到的結(jié)果，Simerics MP+計(jì)算的是調(diào)節(jié)閥從開(kāi)啟到關(guān)閉的全過(guò)程，所得到的流量系數(shù)特性曲線為連續(xù)曲線，更能捕捉到調(diào)節(jié)閥流量特性細(xì)節(jié)。

圖6 流量特性計(jì)算結(jié)果對(duì)比圖

調(diào)節(jié)閥的試驗(yàn)在國(guó)家電投中央研究院牽頭研制的壓水堆核電廠主給水調(diào)節(jié)閥試驗(yàn)臺(tái)架上實(shí)施，如圖7 所示。由圖6 可知，Simerics MP+計(jì)算結(jié)果與試驗(yàn)值的相對(duì)誤差更小，更接近真實(shí)值。表4 是試驗(yàn)值與Simerics MP+、通用CFD 計(jì)算結(jié)果具體數(shù)值的對(duì)比。Simerics MP+的計(jì)算結(jié)果與試驗(yàn)值的最大相對(duì)差值為5.71%，最小相對(duì)差值為0.14%，遠(yuǎn)小于通用CFD 的計(jì)算結(jié)果，更接近于試驗(yàn)結(jié)果。

圖7 試驗(yàn)段樣機(jī)現(xiàn)場(chǎng)圖

表4 試驗(yàn)值與Simerics MP+、通用CFD 結(jié)果對(duì)比

通用CFD 計(jì)算，一方面可再對(duì)網(wǎng)格進(jìn)行細(xì)化處理，以減小誤差，但計(jì)算效率會(huì)大大降低；另一方面對(duì)模型進(jìn)行的簡(jiǎn)化處理會(huì)使相對(duì)誤差增大。

3.2 定性結(jié)果對(duì)比

圖8、圖9 分別是以主給水調(diào)節(jié)閥100%開(kāi)度為例，采用2 種計(jì)算方法得到的壓力場(chǎng)云圖和速度場(chǎng)云圖。由圖8 可知，在壓力場(chǎng)分布方面，兩者的計(jì)算結(jié)果基本吻合，閥門上游處高壓，經(jīng)調(diào)節(jié)閥套筒窗口處節(jié)流降壓，壓降約為0.518MPa，進(jìn)口到出口的壓力變化較為均勻。

圖8 壓力云圖對(duì)比

圖9 速度云圖對(duì)比

由圖9 可知，在速度場(chǎng)分布方面，兩者計(jì)算結(jié)果略有不同，表現(xiàn)在閥門下游區(qū)域流速分布不一致。閥門上游區(qū)域速度變化比較均勻，閥門位置速度變化明顯變大，在閥門套筒窗口處區(qū)域出現(xiàn)了高流速區(qū)，高流速區(qū)對(duì)閥體產(chǎn)生流體沖擊。由于閥門窗口節(jié)流處流通面積較小，結(jié)構(gòu)復(fù)雜，流體流過(guò)窗口時(shí)所受阻滯作用較大，導(dǎo)致閥門下腔室壓力大，閥門上腔室壓力較低，流體流速增大，因而在閥門下腔室節(jié)流處和閥門上腔室均形成渦流，且一直影響到閥門出口。

3.3 對(duì)比結(jié)果分析

通過(guò)對(duì)比，Simerics MP+計(jì)算得到的流量特性曲線能反映更多的流量特性細(xì)節(jié)。由于其計(jì)算的是調(diào)節(jié)閥開(kāi)啟到關(guān)閉的連續(xù)過(guò)程，Simerics MP+憑借其集成的閥門運(yùn)動(dòng)模板對(duì)調(diào)節(jié)閥動(dòng)態(tài)開(kāi)啟或關(guān)閉過(guò)程的實(shí)現(xiàn)比通用CFD 計(jì)算的動(dòng)網(wǎng)格更容易。通用CFD 每10%開(kāi)度進(jìn)行一組流量特性數(shù)值試驗(yàn)的計(jì)算結(jié)果會(huì)忽略一些細(xì)節(jié)。

與本文所采用的通用CFD 計(jì)算結(jié)果相比，Simerics MP+計(jì)算結(jié)果與試驗(yàn)值的相對(duì)誤差更小，更接近真實(shí)值。通用CFD 可以通過(guò)加密網(wǎng)格減小誤差，但會(huì)降低相應(yīng)的計(jì)算效率。相比來(lái)說(shuō)，Simerics MP+能以相對(duì)較少的網(wǎng)格數(shù)量獲得更高的精確度。

4 結(jié)論

本文利用Simerics MP+和通用CFD 仿真數(shù)值分析方法對(duì)核電廠主給水調(diào)節(jié)閥的流量特性進(jìn)行了研究，得到如下結(jié)論。

1）主給水調(diào)節(jié)閥相對(duì)開(kāi)度小于70%時(shí)，為等百分比流量特性；相對(duì)開(kāi)度大于70%時(shí)，接近線性流量特性。此主給水調(diào)節(jié)閥為復(fù)合式流量特性。

2）從分析結(jié)果來(lái)看，與本文所采用的通用CFD計(jì)算相比，Simerics MP+計(jì)算結(jié)果與試驗(yàn)值的相對(duì)誤差更小，最大相對(duì)誤差值為5.71%，表現(xiàn)更優(yōu)。

3）通用CFD 可以通過(guò)加密網(wǎng)格減小誤差，但也會(huì)降低其計(jì)算效率，相比而言，Simerics MP+能以相對(duì)較少的網(wǎng)格數(shù)量獲得更高的精確度，且在閥門動(dòng)態(tài)研究方面具有易實(shí)現(xiàn)的優(yōu)勢(shì)。