胡雪飛，羅 英，王留兵，李 浩，王仲輝，吳冰潔，曹奇鋒，張潤豪

(中國核動力研究設計院，四川 成都 610041)

在以“華龍一號”“玲龍一號”、AP1000等為代表的先進壓水堆型中，常采用一體化/集成式堆頂結(jié)構(gòu)方案。該方案能夠同時兼?zhèn)涠研咎綔y電纜敷設，控制棒驅(qū)動機構(gòu)(CRDM)定位支承、抗震保護、通風散熱等重要功能，還實現(xiàn)了對反應堆壓力容器頂蓋、頂蓋保溫層、CRDM、堆頂電纜等設備的集成式吊裝，簡化了換料工序，節(jié)約了換料時間[1]。在“玲龍一號”即ACP100堆型中，集成式堆頂?shù)耐L冷卻系統(tǒng)為CRDM磁軛線圈散發(fā)出的熱量提供了冷卻流道，如圖1所示：冷卻空氣通過堆頂圍筒上部的窗口進入CRDM線圈區(qū)域，以風冷的形式帶走熱量，再從冷卻圍板上的四個進口進入冷卻風管。ACP100的冷卻風管采用流線型、對稱式布置結(jié)構(gòu)，通過風道四并二、二并一的方式將冷卻空氣引導至出口。冷卻風管的壓降和流場分配能力，是其設計合理性的具體指標；同時，冷卻風管在整個堆頂結(jié)構(gòu)中的重量、重心，也關系到堆頂?shù)钠鸬酢?/p>

圖1 ACP100集成式堆頂結(jié)構(gòu)示意圖Fig.1 The structure of ACP100 integrated head

傳統(tǒng)的堆頂結(jié)構(gòu)通風性能分析多采用試驗的手段。而由于一體化/集成式堆頂結(jié)構(gòu)流場復雜、試驗成本高昂，更多的研究者們近年來逐漸采用CFD方法進行模擬分析。2013年，于浩等人使用CFX對CAP1000堆型一體化堆頂組件的風冷系統(tǒng)進行了模擬計算，驗證了其所選風機的合理性和一體化堆頂設計的科學性[2]。余志偉、何培峰、李燕等人進而耦合溫度場，對M310堆頂冷卻結(jié)構(gòu)流場和溫度場進行了模擬仿真，為冷卻結(jié)構(gòu)優(yōu)化設計提供了科學依據(jù)[3]。2016年，何培峰、許斌、羅英等人又對ACP100堆型中集成式堆頂結(jié)構(gòu)可壓縮氣體與不可壓縮氣體分別進行了模擬計算，研究了氣體可壓縮性對對流傳熱計算結(jié)果的影響[4]。

然而，盡管CFD數(shù)值模擬方法能夠帶來極大的便利，在實際中仍然需要花費大量的精力進行試驗設計和結(jié)構(gòu)調(diào)整。如何簡化選型過程、更高效地對堆頂結(jié)構(gòu)進行功能設計與優(yōu)化，是當前研究中的一個難點。Taguchi方法，又稱為田口法，是由日本田口玄一于20世紀發(fā)明出的關于質(zhì)量工程的數(shù)理方法。該方法以正交表設計試驗，以信噪比作為指標，通過在試驗中調(diào)整可控因素的水平，來改變噪音對質(zhì)量特性的影響[5]，進而提高試驗方案的抗干擾能力。進一步地，Taguchi方法還可以總結(jié)出各因素對質(zhì)量特性影響的顯著性程度，找出各因素的最優(yōu)水平方案，從而得出質(zhì)量特性的最優(yōu)解。

Taguchi方法在試驗設計、結(jié)構(gòu)優(yōu)化、工程質(zhì)量等領域有著廣泛的應用。G.Venkatesan等人利用該方法對波紋板汽水分離器進行了結(jié)構(gòu)優(yōu)化，篩選出了分離性能最優(yōu)的波紋板幾何結(jié)構(gòu)[6]。張啟明等人利用Taguchi方法獲得了神州系列某型號5A06鋁鎂合金過濾器產(chǎn)品電子束焊端面優(yōu)良焊縫，極大地減少了試驗次數(shù)，降低了研發(fā)成本[7]。史晨陽等人運用Taguchi方法設計試驗，優(yōu)化了曲面底板LED陣列的各項結(jié)構(gòu)參數(shù)[8]。王艾萌利用該方法實現(xiàn)了對“U”型轉(zhuǎn)子磁路結(jié)構(gòu)的內(nèi)置式永磁同步電機的結(jié)構(gòu)參數(shù)優(yōu)化設計，使各性能指標相較優(yōu)化前有了顯著的提升[9]。Taguchi方法在ACP1000一體化堆頂結(jié)構(gòu)的流場優(yōu)化中有過試算應用，符合工程實踐結(jié)果，顯示出了良好的適用性和科學性。

本文以ACP100集成式堆頂中的冷卻風管結(jié)構(gòu)為研究對象，利用Taguchi方法設計試驗，使用Fluent 14.5對各試驗方案進行數(shù)值模擬研究，通過Minitab 17進行后續(xù)數(shù)理計算，探討各結(jié)構(gòu)因素和水平對冷卻風管性能的影響，求出最優(yōu)因素水平組合并加以驗證，實現(xiàn)對風管結(jié)構(gòu)的優(yōu)化，提高工作性能。

1 Taguchi方法建模

1.1 物理模型

圖2所示的是ACP100集成式堆頂?shù)睦鋮s風管模型圖。冷卻風管外形結(jié)構(gòu)是以出風管口中心線為軸對稱的殼體，分為底層、中層和頂層三層流道。4個進風口跨中均布在同一平面上。冷卻空氣從進風口沿豎直方向向上進入風管流道，并沿著風管形狀流動。風管中的三層流道逐漸合并，冷卻空氣最終在風管頂部匯聚，從出風口排出。為了便于后續(xù)CFD建模，將出風管路合理延長，確保數(shù)值模擬過程中流場充分發(fā)展。

圖2 ACP100冷卻風管模型圖Fig.2 The integrated head structure of the cooling vent of ACP100

冷卻空氣通過風管后的壓降是衡量其性能的重要指標，在保證通風量滿足CRDM線圈散熱要求的前提下，盡可能地減少冷卻風阻，可以有效地提高堆頂?shù)慕?jīng)濟性能。

1.2 Taguchi試驗設計

由于冷卻風管隸屬于集成式堆頂冷卻系統(tǒng)，結(jié)構(gòu)優(yōu)化應在不影響相關外部接口的前提下進行。故冷卻風管外形尺寸、進風口尺寸與位置、出風口尺寸與位置應保持不變，不能作為Taguchi方法的可控因素。同時，選擇的因素應相互獨立，具備非相關性。如圖3所示，建立空間坐標軸，選取中層流道中心線角度坐標α，中層流道寬度L，中層流道高度H，出風口接管圓角半徑r作為可控因素，各因素及水平的選擇見表1。其中，case1冷卻風管為ACP100集成式堆頂結(jié)構(gòu)中采用的風管結(jié)構(gòu)。

圖3 冷卻風管四因素選擇示意圖Fig.3 Four factors chosen of cooling vent

表1 冷卻風管因素、水平選擇表Table 1 Selection of factors and levels for cooling vent

根據(jù)四因素、三水平選擇L9(34)正交表，并設計試驗方案如表2所示。9種試驗風管的三維流道模型如圖4所示。

圖4 9種試驗風管的三維流道模型圖Fig.4 3D model for 9 cases of the cooling vent

表2 L9(34)Taguchi方法試驗設計方案Table 2 Test design scheme of L9(34) Taguchi method

2 CFD模擬設置

2.1 計算網(wǎng)格模型

使用ANSYS ICEM為圖4 9種風管模型進行計算網(wǎng)格劃分。按照模型的幾何形狀，合理劃分block區(qū)域，采用六面體結(jié)構(gòu)化網(wǎng)格。以案例(1)為例，分別劃分出308 616、556 301和716 103數(shù)量的網(wǎng)格，結(jié)合2.2節(jié)中的求解器設置，對計算收斂后進出口的壓降采用面積加權(quán)平均法進行比較，不同數(shù)量的網(wǎng)格壓降值及相對誤差如表3所示?？芍?，案例(1)取556 301網(wǎng)格數(shù)時已經(jīng)能夠較好地滿足精度要求，因此選擇該網(wǎng)格作為后續(xù)試驗中案例(1)風管的網(wǎng)格。其他8種風管結(jié)構(gòu)采用相近的網(wǎng)格劃分方法，整體網(wǎng)格數(shù)量保證在52～57萬，網(wǎng)格雅各比質(zhì)量均在0.8以上，局部網(wǎng)格結(jié)構(gòu)細節(jié)如圖5所示。

表3 案例(1)網(wǎng)格獨立性檢驗及模擬精度Table 3 Grid independence test andsimulation precision of case(1)

圖5 案例(1)局部網(wǎng)格細節(jié)圖Fig.5 The details of part grids of case(1)

2.2 求解器設置

使用ANSYS Fluent14.5對冷卻風管進行數(shù)值模擬。流體選擇對應工況溫度70 ℃下的空氣，認為流體不可壓縮。選取能對管內(nèi)流動、邊界層流動等多類型湍流進行高精度求解的Realizablek-ε方程。

在設置邊界條件時，模擬實際工況中抽風機的工作方式，將出風口邊界條件設置為速度流量出口，速度值為20 m/s，對應實際抽風機41 500 m3/h的風量。設置4個入口為壓力入口，壓力值為0 Pa。

選取SIMPLE算法進行流體的壓力和速度耦合，選擇二階迎風格式作為空間離散格式。在監(jiān)測各迭代參數(shù)殘差時，控制殘差值≤10-4。

3 結(jié)果與討論

3.1 冷卻風管流場特性

通過對9種風管流場進行模擬，可以得到各風管內(nèi)空氣的峰值速度對比，如圖6所示。從圖6中可以看到，案例(4)風管的峰值速度最高，達到了60.1 m/s。按照式(1)計算當前工況下風管內(nèi)空氣的密度變化情況[10]。

圖6 9種風管模型峰值速度對比Fig.6 Peak velocities comparison for 9 cases of the cooling vent

(1)

式中:ρ——峰值速度下空氣的密度；

ρs——當前溫度工況下空氣靜止時的密度；

u——峰值速度;

空氣比熱γ=1.399;

空氣常數(shù)R=287 m2/s2·K;

工況溫度T=(273+70) K。

計算得ρ/ρs=0.987，即空氣以峰值速度60.1 m/s運動時，與靜止狀態(tài)下的空氣密度相對變化約1.3%，因此在2.2節(jié)中認為流體不可壓縮的條件設置是合理的。

圖7 9種風管模型內(nèi)部流場速度云圖Fig.7 Velocity magnitudes for 9 cases of the cooling vent

圖7所示的，是9種風管模型的流場速度云圖。冷卻空氣從同一側(cè)的兩個進風口進入底層流道后，向中層流道流動并發(fā)生匯聚；在頂層流道與出風口的交點處，來自不同側(cè)的兩股流體再次匯聚，向出口管口流出。受總體接口尺寸的影響，中層流道高度減小，底層流道的腔室空間增大。當流體從進風口進入并向上方流動時，會在中層流道交匯區(qū)下方形成漩渦，底層流道的腔室空間越大，漩渦的尺度也益發(fā)明顯，如案例(1)、案例(3)、案例(6)、案例(7)所示。

空氣在中層流道匯聚的過程中，隨著風管壁面形狀的幾何彎曲，流體運動方向也發(fā)生急速偏轉(zhuǎn)，并在發(fā)生偏轉(zhuǎn)且靠近風管壁面區(qū)域出現(xiàn)較高的速度，這意味著較大的能量損失。一般情況下，由于進風口2和4在布置方向上更靠近出風管路，從這兩個入口進入風管的流體在整個運動過程中的速度方向偏轉(zhuǎn)是最為激烈的；而進風口1和3遠離出風管路，從這兩個入口進入的空氣在偏轉(zhuǎn)過程中，不會出現(xiàn)180°急速變向，運動較為舒緩。因此，在中層流道彎曲近壁面靠近進風口2和4的位置附近常會出現(xiàn)流場的峰值速度。以案例(1)為例，達到了53.8 m/s。

通過改變中層流道幾何特性(調(diào)整α，L，H)，可以明顯變更峰值速度區(qū)的出現(xiàn)位置。以案例(4)為例，中層流道中心線位置向出口管路方向偏轉(zhuǎn)15°后，從進風口2和4進入的空氣不再出現(xiàn)激烈的偏折，速度梯度變化也相對緩和。而由于底層空間高度的壓縮，從進風口1和3進入的空氣流道更加狹窄，峰值速度的出現(xiàn)位置改在此側(cè)的中層流道彎曲近壁面附近出現(xiàn)。

在頂層流道與出口接管相貫處增加圓角，也可以明顯改善該區(qū)域的壓降損失。案例(1)、案例(5)和案例(9)均采用了ACP100風管的0圓角結(jié)構(gòu)，從圖7中可以發(fā)現(xiàn)，由于相貫壁面幾何尺寸的尖銳變化，空氣在經(jīng)過尖銳點時會受到較大的擾動，出現(xiàn)高速點，進而在流場后方出現(xiàn)狹長的低速區(qū)，使出口管路內(nèi)的流場遲遲不能發(fā)展均勻，這意味著明顯的壓降損失。而設置圓角過度后，流體圓滑經(jīng)過管路相交處，后續(xù)流場可以很快發(fā)展均勻，如案例(2)、案例(4)和案例(6)所示。

通過監(jiān)測9種風管的進出口靜壓力，采用面積加權(quán)平均法處理，得到9種風管的進出口壓降如表4所示。從表中可以看到，9種風管的進出口壓降差別顯著，其中，案例(5)風管的壓降值最高，為1 396.07 Pa；案例(7)風管的壓降值最低，為594.58 Pa。案例(1)風管，即當前ACP100堆頂冷卻風管結(jié)構(gòu)并不能顯示出明顯的優(yōu)勢，壓降高達1 181.83 Pa。

表4 9種風管方案的進出口壓降對比Table 4 Pressure drops comparison for 9cases of the cooling vent

3.2 Taguchi試驗結(jié)果

為了進一步分析各因素、水平對風管性能的影響，使用Minitab 17軟件中的Taguchi方法進行下一步處理。由于壓降越低，風管性能越佳，因此本案例屬于Taguchi方法中的望小特性，其信噪比公式如式(2)所示：

(2)

式中，S/N——信噪比；

n——試驗次數(shù)；

yi——輸出性能特性

以9種方案的壓降作為輸出性能特性，代入式(2)進行計算，得到了Taguchi信噪比主效應圖，如圖8所示。信噪比均值為-59.39，減小中層流道中心角度、減小中層流道寬度、增大中層流道高度都會使信噪比明顯降低，即增大風管的壓降損失；相反，則可增大信噪比，改善風管內(nèi)的壓降情況。究其原因，這三個因素基本決定了風管內(nèi)的流道曲折形狀，顯著影響空氣的流程阻力。利用Taguchi方法，通過比較不同因素的信噪比極差，可以得到本試驗中各因素對壓降的影響程度排秩：出口圓角r>中層流道高度H>中層流道寬度L>中層流道中心線角度坐標α。

圖8 Taguchi信噪比主效應圖Fig.8 Signal-to-noise ratio in Taguchi

從圖8中還可看到，α3L3H1r2為最優(yōu)的因素水平組合方案。該組合與此前9種方案完全不同。使用Taguchi方法對該組合方案進行預測，信噪比高達-54.19。為了進一步驗證該方案的合理性，對α3L3H1r2組合進行建模，并采用同樣的數(shù)值模擬方法進行試驗。優(yōu)化后的風管三維流道模型如圖9所示。

圖9 使用Taguchi方法優(yōu)化后的風管三維流道模型Fig.9 3D model forthe cooling vent by Taguchi optimization

3.3 優(yōu)化后風管的性能驗證

圖10是優(yōu)化后的風管流場速度云圖。由于中層流道中心線角度坐標向遠離出風口方向偏移，且拓寬了中層流道寬度，從進風口1和3流進的空氣不需要經(jīng)過明顯的偏折就可以相繼進入中層和頂層流道，減小了流動過程的能量損耗。同時，從進風口2和4處流進的流體在底層和中層流道中有足夠的長行程實現(xiàn)偏折，流程變化緩和，故峰值速度低至41.5 m/s。此外，由于中層流道布置位置的變化，流體在頂層流道中的行程也相應增長，截面速度能有足夠的行程趨于穩(wěn)定分布。最后，伴隨著出口圓角的設置，兩側(cè)的流體在頂層匯聚進入出口管路后，未出現(xiàn)案例(1)中的尖銳速度點與狹長低速區(qū)。

圖10 經(jīng)Taguchi方法優(yōu)化后的風管流場速度云圖Fig.10 Velocity magnitude in the cooling vent by Taguchi optimization

將優(yōu)化后風管的進出口壓降同前述的9種方案進行對比，可知優(yōu)化后的壓降低至566.74 Pa，低于9種方案的壓降，符合Taguchi方法的預測結(jié)果，如圖11所示。除了使壓降盡可能小之外，優(yōu)良的冷卻風管結(jié)構(gòu)還需要考慮對堆頂圍筒內(nèi)部及CRDM線圈附近的流場均勻分配的功能。使用面積加權(quán)平均法統(tǒng)計各方案中四個進風口的截面速度，計算其標準差作為風管對流場分配效果的衡量指標。仍由圖11所示，優(yōu)化后風管四進風口截面氣速標準差為0.87，僅稍高于案例(1)(標準差0.41)和案例(2)(標準差0.31)，明顯低于案例(3)至案例(9)的對應結(jié)果。

圖11 各方案風管的壓降及速度分配特性Fig.11 Pressure drops and velocity distributions for various cases of the cooling vent

圖12所示的，是各方案風管的水平重心與豎直重心坐標位置對比。因各風管結(jié)構(gòu)以X坐標軸為軸對稱，故X軸重心即為水平方向重心。由圖可知，案例(1)的水平重心坐標為289.9 mm，豎直重心坐標為972.5 mm，受到頂層流道結(jié)構(gòu)產(chǎn)生的力臂影響，整個風管重心處于偏心狀態(tài)。案例(2)至案例(9)風管受到中層流道結(jié)構(gòu)的影響，水平重心坐標處于214.8 mm至381.2 mm之間。而經(jīng)Taguchi方法優(yōu)化后的風管豎直重心坐標為1 018.5 mm，水平重心僅為207.8 mm，水平重心小于全部9種方案。對比案例(1)風管，優(yōu)化后的風管雖然豎直方向重心增大了4.73%，但偏心量減少了28.3%，這對提高整個堆頂?shù)目拐鹉芰?、減小起吊工況下堆頂?shù)蹙叩钠妮d荷有著積極的作用。

圖12 各方案風管的重心坐標對比Fig.12 Gravity center comparison for various cases of the cooling vent

因此，經(jīng)過Taguchi方法優(yōu)化后的風管結(jié)構(gòu)具有優(yōu)異的性能，滿足優(yōu)化目的。

4 結(jié)論

本文應用Taguchi方法及CFD數(shù)值模擬方法，以ACP100集成式堆頂結(jié)構(gòu)冷卻風管為研究對象，得到了主要結(jié)論如下：

1)驗證了Taguchi方法在ACP100集成式堆頂冷卻風管結(jié)構(gòu)研究和優(yōu)化過程中的合理性與可行性。

2)得到了冷卻風管結(jié)構(gòu)因素對壓降的影響程度排秩為：出口圓角r>中層流道高度H>中層流道寬度L>中層流道中心線角度坐標α。

3)設計出了信噪比最優(yōu)的風管結(jié)構(gòu)優(yōu)化方案，并通過CFD數(shù)值模擬的方法驗證其合理性。最終認為，全新風管結(jié)構(gòu)具備壓降低、堆頂流場速度分配均勻度高以及重心偏心量低的優(yōu)點，可作為ACP100集成式堆頂結(jié)構(gòu)中風管結(jié)構(gòu)優(yōu)化的較好方案。