核主泵屏蔽電機(jī)內(nèi)部水路三維流場研究

李藏雪，趙博敏，路義萍，呂向平

(1.哈爾濱電氣動(dòng)力裝備有限公司，黑龍江哈爾濱150066 2.哈爾濱理工大學(xué)機(jī)械動(dòng)力工程學(xué)院，黑龍江 哈爾濱150080)

核主泵屏蔽電機(jī)內(nèi)部水路三維流場研究

李藏雪1，趙博敏2，路義萍2，呂向平1

(1.哈爾濱電氣動(dòng)力裝備有限公司，黑龍江哈爾濱150066 2.哈爾濱理工大學(xué)機(jī)械動(dòng)力工程學(xué)院，黑龍江 哈爾濱150080)

核主泵屏蔽電機(jī)是核電站核島一次回路系統(tǒng)的重要組成部分，而屏蔽電機(jī)內(nèi)部的冷卻水的流動(dòng)狀態(tài)直接影響電機(jī)內(nèi)的熱量傳遞的效果。對(duì)此建立了一種屏蔽電機(jī)整機(jī)的三維流場模型?；谟?jì)算流體力學(xué)理論中的有限體積法，利用Ansys Fluent軟件，以實(shí)驗(yàn)得到的額定工況時(shí)的一次水流速為邊界條件，反演計(jì)算出屏蔽電機(jī)內(nèi)部的流場速度分布。針對(duì)結(jié)果分別對(duì)屏蔽電機(jī)內(nèi)的關(guān)鍵部位的三維流場進(jìn)行了分析，得出屏蔽電機(jī)內(nèi)部的軸向及徑向速度分布規(guī)律。計(jì)算結(jié)果表明，在靠近旋轉(zhuǎn)壁面處，水速較大，且沿軸向以螺旋線形式向上流動(dòng)，在下飛輪、屏蔽套間隙流體沿徑向速度梯度較大。該研究可為后續(xù)的溫度場研究提供不同位置處速度分布和數(shù)值。

屏蔽電機(jī)；三維流場；有限體積法；計(jì)算流體力學(xué)

隨著我國環(huán)保意識(shí)的不斷增強(qiáng)和國內(nèi)面臨的日益嚴(yán)重的能源緊缺問題，對(duì)潔凈發(fā)電技術(shù)的發(fā)展需求不斷增加。核電作為一種新型能源在國民經(jīng)濟(jì)中的地位顯得越來越重要。核主泵屏蔽電機(jī)是核電站核島一次回路系統(tǒng)的重要組成部分，其安全穩(wěn)定運(yùn)行對(duì)冷卻堆芯以及防止發(fā)生核電站事故起著非常關(guān)鍵的作用[1]。而屏蔽電機(jī)內(nèi)部的冷卻水的流動(dòng)狀態(tài)直接影響電機(jī)內(nèi)的熱量傳遞的效果，因而對(duì)屏蔽電機(jī)內(nèi)部流場進(jìn)行研究非常必要。

近年來，國內(nèi)外專家學(xué)者對(duì)大型汽輪發(fā)電機(jī)和水輪發(fā)電機(jī)、風(fēng)力發(fā)電機(jī)內(nèi)溫度場及流體場進(jìn)行了大量的研究。姚若萍等基于有限元法計(jì)算了采用蒸發(fā)冷卻系統(tǒng)的大型水輪發(fā)電機(jī)在不同工況下的定子三維溫度場，得出了不同工況下定子鐵芯及股線的溫度分布[2]。路義萍等采用有限體積法對(duì)空冷汽輪發(fā)電機(jī)流場進(jìn)行了研究[3]。Kuosa等利用ABB公司的300MW空冷汽輪發(fā)電機(jī)研發(fā)測試實(shí)驗(yàn)臺(tái)，在充分發(fā)展的條件下，考慮轉(zhuǎn)子空氣射流，采用三孔“cobra”傳感器測量氣隙內(nèi)全壓與靜壓，并進(jìn)行數(shù)值模擬計(jì)算[4]。國內(nèi)外對(duì)大型空冷汽輪發(fā)電機(jī)、空冷電動(dòng)機(jī)等經(jīng)過數(shù)值化、精細(xì)化研究其三維流場與溫度場變化規(guī)律已經(jīng)相當(dāng)成熟，而在屏蔽電機(jī)方面，絕大多數(shù)是對(duì)其電磁參數(shù)，受力分析等方面研究較多，對(duì)整體的流場溫度場進(jìn)行的研究卻很少。如孫桂瑛等分析了屏蔽電機(jī)磁場的特點(diǎn)，用解析法求出了屏蔽電機(jī)的磁場分布，進(jìn)而計(jì)算出屏蔽層的渦流損耗，并通過實(shí)驗(yàn)驗(yàn)證了解析方法的可行性[5]。馮穎慧通過對(duì)屏蔽電機(jī)主泵電機(jī)內(nèi)屏蔽套的應(yīng)力分析，得出各載荷在電機(jī)內(nèi)屏蔽套上產(chǎn)生的應(yīng)力低于屏蔽套材料的屈服強(qiáng)度，可保證屏蔽套的完整性[6]。丁樹業(yè)等以一臺(tái)5500kW核主泵屏蔽電機(jī)為例，建立流體場與溫度場耦合的求解域物理數(shù)學(xué)模型，通過計(jì)算揭示了核主泵屏蔽電機(jī)內(nèi)溫度的分布規(guī)律[7]。成德等將核主泵電機(jī)內(nèi)部流場根據(jù)結(jié)構(gòu)特點(diǎn)分成四部分，針對(duì)每一部分分別建立物理模型，然后分別模擬計(jì)算[8]。本文不同于文獻(xiàn)[7-8]，以一種典型的外置換熱器的屏蔽電機(jī)為例，采用計(jì)算流體力學(xué)(CFD)方法，選用剪切應(yīng)力傳輸(SST)k-ω兩方程湍流模型等，研究額定運(yùn)行工況下電機(jī)內(nèi)的流場變化情況，此外，本文建立了屏蔽電機(jī)的整機(jī)模型，能夠更好地觀察其內(nèi)部的流場的變化規(guī)律，所得結(jié)論為屏蔽電機(jī)冷卻計(jì)算及設(shè)計(jì)提供參考。

1 物理模型與網(wǎng)格劃分

本文研究的電機(jī)為立式鼠籠三相四極異步屏蔽電動(dòng)機(jī)，額定轉(zhuǎn)速為1500r/min，頻率為50Hz，定子繞組形式為雙層短距。采用Solidworks軟件建立了整機(jī)的物理模型，如圖1所示。

圖1 整機(jī)模型Fig.1 The whole machine model

對(duì)建立的三維物理模型進(jìn)行布爾運(yùn)算得到整機(jī)的流場模型。外置換熱器中流出的冷卻水由機(jī)座底部軸芯孔處吸入后，向上流動(dòng)至輔葉輪處沿軸芯孔徑向依靠旋轉(zhuǎn)產(chǎn)生的離心力甩出，一部分向下流動(dòng)冷卻及潤滑轉(zhuǎn)子下部空腔中的下導(dǎo)軸承及下飛輪內(nèi)裝置后返回到軸心水冷通道中，與新進(jìn)入主循環(huán)冷水匯合返回軸芯孔下部向上流動(dòng)；另一部分直接向上流動(dòng)繼續(xù)冷卻定轉(zhuǎn)子屏蔽套及上部軸承，最終由集流腔經(jīng)封頭上的出口流出進(jìn)入外置冷卻器，如圖2所示。

圖2 水路示意圖Fig.2 Schematic diagram of water route

本文利用Ansys icem cfd軟件進(jìn)行網(wǎng)格劃分。為了使計(jì)算結(jié)果更準(zhǔn)確，收斂速度更快，本文對(duì)整機(jī)采用了分塊結(jié)構(gòu)網(wǎng)格劃分，整機(jī)全部為六面體結(jié)構(gòu)網(wǎng)格，質(zhì)量較高。

2 數(shù)學(xué)模型與邊界條件

2.1 數(shù)學(xué)模型

本文采用有限體積法進(jìn)行數(shù)值求解，其基本思想是將計(jì)算區(qū)域劃分為網(wǎng)格，并使每個(gè)網(wǎng)格節(jié)點(diǎn)周圍有一個(gè)互不重復(fù)的控制體積，將待求解的微分方程對(duì)每一個(gè)控制體積積分，從而得出離散方程。

由于入口雷諾數(shù)大于2300，處于湍流狀態(tài)。計(jì)算域內(nèi)流體為水，視為不可壓縮流體，利用Fluent 軟件，選擇多重參考系計(jì)算方法。在旋轉(zhuǎn)坐標(biāo)系和固定坐標(biāo)系下，聯(lián)合求解描述電機(jī)內(nèi)部冷卻水速度、壓力等物理場的湍流時(shí)均質(zhì)量、動(dòng)量守恒方程。其中，旋轉(zhuǎn)坐標(biāo)系下湍流時(shí)均質(zhì)量、動(dòng)量守恒方程和絕對(duì)速度矢量u與相對(duì)速度矢量ur的關(guān)系如式(1)～式(3)所示。

(ρur)=0

(1)

(ρurur)+ρ(2Ω×ur+Ω×Ω×r)=
-p+τ+F

(2)

u=ur+Ω×r

(3)

固定坐標(biāo)系下湍流時(shí)均質(zhì)量、動(dòng)量守恒方程分別如式(4)、式(5)所示。

(ρu)=0

(4)

(5)

ρ——密度；

Ω——為旋轉(zhuǎn)角速度矢量；

r——轉(zhuǎn)動(dòng)坐標(biāo)系中微元體的位置矢量；

ρ(2Ω×ur+Ω×Ω×r)——科里奧里力；

F——微元體上的體積力；

τ——因分子粘性作用而產(chǎn)生的作用于微元體表面的粘性應(yīng)力；

p——壓力；

u、v、w——絕對(duì)速度矢量u在三個(gè)坐標(biāo)軸方向的速度分量。

關(guān)于湍流數(shù)學(xué)模型[9、10]，由于沿徑向流動(dòng)空間尺寸在毫米數(shù)量級(jí)，屬于邊界層粘性剪切力支配的繞圓柱的強(qiáng)制對(duì)流流動(dòng)區(qū)域內(nèi)，沿徑向壓力速度梯度非常大，因而本文采用了Shear-Stress Transport(SST)k-ω模型，表達(dá)式如下：

(6)

式中：Gk——層流速度梯度而產(chǎn)生的湍流動(dòng)能；

Gω——ω方程產(chǎn)生的湍流動(dòng)能；

Tk——k的擴(kuò)散率；

Tω——ω的擴(kuò)散率；

Yk——擴(kuò)散產(chǎn)生的湍流；

Yω——擴(kuò)散產(chǎn)生的湍流；

SK——用戶自定義；

Sω——用戶自定義。

2.2 邊界條件

電機(jī)內(nèi)部一次水入口采用速度入口邊界條件，大小根據(jù)實(shí)驗(yàn)測量得到的輔葉輪額定負(fù)荷下的工作點(diǎn)流量折算出具體數(shù)值，出口為自由出流；轉(zhuǎn)子轉(zhuǎn)速為1500r/min；針對(duì)本文屏蔽電機(jī)，在多重參考系下求解，為了符合計(jì)算模型要求，將流體區(qū)域分成了旋轉(zhuǎn)流體區(qū)與靜止(相對(duì)轉(zhuǎn)速為零)流體區(qū)兩部分。旋轉(zhuǎn)流體區(qū)：該區(qū)域位于轉(zhuǎn)軸體區(qū)域內(nèi)部孔隙中的水，隨主軸一起轉(zhuǎn)動(dòng)，其他均屬于靜止流體區(qū)；另外，固體壁面分為兩類，一類是靜止壁面，另一類是隨轉(zhuǎn)軸一起旋轉(zhuǎn)的壁面，即具有切向速度分量的旋轉(zhuǎn)壁面，分別設(shè)置轉(zhuǎn)速。

計(jì)算過程中，所有網(wǎng)格節(jié)點(diǎn)的離散方程組采用分離、隱式求解，其中，壓力與速度耦合方程采用SIMPLEC算法，最終，方程組獲得穩(wěn)定的網(wǎng)格獨(dú)立收斂解。計(jì)算過程中，殘差取1×10-3。

3 計(jì)算結(jié)果與分析

3.1 速度跡線分析

圖3所示為計(jì)算域內(nèi)用速度數(shù)值表示的內(nèi)部一次水的流動(dòng)跡線圖，圖中箭頭所指方向?yàn)榱黧w流動(dòng)方向。從圖中可以看出流體由機(jī)座底部軸芯孔處吸入后，向上流動(dòng)至輔葉輪處沿軸芯孔徑向甩出，一部分向下流動(dòng)，返回到軸芯水冷通道中，與新進(jìn)入主循環(huán)冷卻水匯合；另一部分直接向上呈螺旋線流動(dòng)形式流經(jīng)定轉(zhuǎn)子屏蔽套的空腔及上部導(dǎo)軸承，最終由集流腔經(jīng)封頭上的出口流出進(jìn)入外置冷卻器。流動(dòng)跡線與實(shí)際工況下流體流動(dòng)相一致。圖中最大速度82.44m/s位于下飛輪間隙中，為了分析方便，以下按由下到上順序分別輸出重要局部部件中的流動(dòng)特征并加以分析。

圖3 計(jì)算域內(nèi)用速度數(shù)值表示的跡線圖Fig.3 Trace map stated as speed value of computational domain

3.2 輔葉輪及軸芯孔隙內(nèi)流體分析

圖4所示為輔葉輪及軸芯孔隙內(nèi)流體運(yùn)動(dòng)速度矢量圖。輔葉輪處于轉(zhuǎn)子轉(zhuǎn)軸內(nèi)，隨著轉(zhuǎn)子的轉(zhuǎn)動(dòng)，輔葉輪處的流體也跟著旋轉(zhuǎn)，從矢量圖中可以明顯看出流體的轉(zhuǎn)動(dòng)方向，速度最大處為輔葉輪出口處，達(dá)到39.76m/s，即出口處半徑r最大，線速度達(dá)到最大，與理論線速度計(jì)算一致。

圖4 輔葉輪及軸芯孔隙內(nèi)流體運(yùn)動(dòng)速度矢量圖Fig.4 Velocity vector diagram of fluid in auxiliary impeller and shaft poles

圖5所示為輔葉輪及軸芯孔隙內(nèi)流體靜壓分布云圖，從圖中可以看出最大壓力為輔葉輪出口，即旋轉(zhuǎn)半徑最大處，相對(duì)靜壓出口的靜壓數(shù)值為378809.75Pa，最低壓力在輔葉輪旋轉(zhuǎn)中心處。由于輔葉輪的旋轉(zhuǎn)使旋轉(zhuǎn)中心處形成了負(fù)壓，將下部底蓋內(nèi)部的水吸入軸芯孔中，并向上運(yùn)動(dòng)，到達(dá)輔葉輪后，在離心力的作用下，克服流動(dòng)阻力，被沿徑向甩出?？梢?，輔葉輪起到了生成負(fù)壓將流體吸入，同時(shí)，使甩出的水升壓，能夠克服后續(xù)流動(dòng)阻力的作用。

圖5 輔葉輪及軸芯孔隙內(nèi)流體靜壓分布云圖Fig.5 Contours of static pressure of fluid in auxiliary impeller and shaft pores

3.3 下飛輪間隙流體分析

圖6所示為下飛輪間隙處的流體運(yùn)動(dòng)速度矢量圖。該區(qū)域的固體部件本身沒有熱源，流體主要用來冷卻下飛輪在運(yùn)行過程中水摩擦產(chǎn)生的熱量，流體速度主要由高速旋轉(zhuǎn)的下飛輪帶動(dòng)，因而該部分流體在流動(dòng)過程中也在不斷轉(zhuǎn)動(dòng)。從矢量圖中可以明顯看出流體的旋轉(zhuǎn)，內(nèi)側(cè)即靠近旋轉(zhuǎn)壁面處的水流速為屏蔽電機(jī)中速度峰值，其數(shù)值為82.44m/s；而外側(cè)即靠近靜止壁面流體速度較小，壁面位置處水的速度為0，與理論分析一致。因此，下飛輪中水的剪切摩擦損耗將較大，合理設(shè)計(jì)結(jié)構(gòu)，保證下飛輪中的再循環(huán)冷卻潤滑水的流量很重要。

圖6 飛輪間隙流體運(yùn)動(dòng)速度矢量圖Fig.6 Velocity vector diagram of fluid in the flywheel clearance

圖7 屏蔽套間隙流體運(yùn)動(dòng)速度云圖Fig.7 Contours of velocity of fluid in shielding clearance(a) 軸向速度變化;(b) 徑向速度變化

3.4 屏蔽套間隙流體分析

圖7所示為定轉(zhuǎn)子屏蔽套間隙內(nèi)的流體運(yùn)動(dòng)速度云圖。圖7(a)為流體軸向速度分布云圖，圖中流體從下向上流動(dòng)，與定轉(zhuǎn)子屏蔽套間進(jìn)行對(duì)流換熱，帶走熱量。圖7(b)為流體徑向速度分布云圖?？拷D(zhuǎn)子屏蔽套側(cè)的流體速度最大，達(dá)到48.31m/s。由計(jì)算流體力學(xué)理論可知轉(zhuǎn)子屏蔽套外側(cè)的線速度應(yīng)與屏蔽套間隙內(nèi)流體運(yùn)動(dòng)的最大速度相等，從圖7(b)中可以看出流體最大速度與理論計(jì)算值相等，最小速度0m/s發(fā)生在靠近定子屏蔽套壁面位置。

3.5 上導(dǎo)軸承間隙流體分析

圖8為上導(dǎo)軸承間隙流體運(yùn)動(dòng)速度分布云圖。圖8(a)圖8(b)分別為軸承內(nèi)外側(cè)間隙流體速度分布云圖。該部分流體主要是用來潤滑并冷卻軸承在工作過程中產(chǎn)生的摩擦損耗。從圖中可以看出流體在進(jìn)入上導(dǎo)軸承前的匯流腔內(nèi)速度較大，進(jìn)入導(dǎo)軸承后速度明顯減小，這是由于在匯流腔內(nèi)靠近轉(zhuǎn)軸處，存在無滑移邊界，在黏性力作用下，存在一流體薄層，該層內(nèi)速度梯度較大，流體最大速度在靠近轉(zhuǎn)子軸旋處，而在上導(dǎo)軸承間隙內(nèi)壁面都是靜止壁面，因而速度較小。

圖8 上導(dǎo)軸承間隙流體運(yùn)動(dòng)速度云圖Fig.8 Contours of velocity of fluid in the guide bearing(a) 軸承內(nèi)側(cè);(b) 軸承外側(cè)

3.6 封頭內(nèi)流體分析

圖9(a)所示為封頭內(nèi)的流體速度云圖，圖9(b)為該處出口徑向截面處的流體速度分布云圖。

圖9 封頭內(nèi)流體速度分布云圖Fig.9 Contours of velocity of fluid in the head(a) 封頭內(nèi)整體;(b) 徑向截面

該處的流體經(jīng)由上導(dǎo)軸承的軸瓦與軸承縫隙的流體匯集到封頭內(nèi)的匯流腔內(nèi)，最后通過封頭出口流出電機(jī)。在匯流腔內(nèi)的流動(dòng)的水能夠帶走上部飛輪傳遞下來的熱量，同時(shí)也在冷卻封頭。該處流體即將到達(dá)整個(gè)電機(jī)的出口，從圖中可以看出流體貼壁處速度較低，只有在中心處靠近轉(zhuǎn)子軸處，旋轉(zhuǎn)壁面帶動(dòng)與之相鄰的水運(yùn)動(dòng)，最大速度達(dá)到32.67m/s。

根據(jù)以上各部件的速度云圖及矢量圖，可以看出速度軸向及徑向分布特點(diǎn)及數(shù)值范圍，在采用有限元法進(jìn)行溫度場計(jì)算時(shí)，可以采用對(duì)應(yīng)的對(duì)流傳熱系數(shù)實(shí)驗(yàn)關(guān)聯(lián)式算出所要計(jì)算的壁面上的平均對(duì)流傳熱系數(shù)，進(jìn)行后續(xù)計(jì)算，或者在采用有限體積法進(jìn)行溫度場計(jì)算時(shí)，輸出某些壁面上的對(duì)流傳熱系數(shù)，與實(shí)驗(yàn)關(guān)聯(lián)式算出的數(shù)值相互對(duì)比，分析并檢驗(yàn)計(jì)算結(jié)果的準(zhǔn)確性。

3.7 相鄰部件銜接處的流體分析

此外，本物理模型的另一優(yōu)點(diǎn)是可以得到相鄰兩部件銜接處的流體通道中，流體過渡流動(dòng)的特征，見圖10，該特征是文獻(xiàn)[8]所不能提供的。圖10所示為飛輪軸承輔葉輪銜接處流體的速度矢量圖。在該處流動(dòng)復(fù)雜，水由軸芯孔進(jìn)入輔葉輪、從輔葉輪甩出到進(jìn)入?yún)R流腔、再分成兩部分，整個(gè)過程非常清晰，可以得到水的速度大小、方向等詳細(xì)信息，避免了人為的中間截面上的假設(shè)。這與文獻(xiàn)[8]分開建模，分別設(shè)置入口邊界條件進(jìn)行分析有著很大的不同。

圖10 飛輪輔葉輪軸承銜接處速度矢量圖Fig.10 Velocity vector diagram of fluid in the flywheel、auxiliary impeller and guide bearing

4 結(jié)論

本文利用有限體積法，采用了Shear-Stress Transport(SST)k-ω模型，針對(duì)核主泵屏蔽電機(jī)在正常運(yùn)行工況下內(nèi)部復(fù)雜的流場變化進(jìn)行了分析，在轉(zhuǎn)子轉(zhuǎn)速和入口速度已知的條件下，利用Ansys Fluent軟件計(jì)算屏蔽電機(jī)內(nèi)部流場的速度壓力分布。針對(duì)結(jié)果分別對(duì)屏蔽電機(jī)內(nèi)的關(guān)鍵部位的流場進(jìn)行了分析，所得結(jié)果與理論分析相一致。得出的結(jié)論如下：

最大壓力在輔葉輪出口，即旋轉(zhuǎn)半徑最大處，最低壓力在輔葉輪旋轉(zhuǎn)中心處，最大速度在飛輪貼近旋轉(zhuǎn)壁面處，水的流動(dòng)跡線表明，在靠近旋轉(zhuǎn)壁面處，水速較大，且沿軸向以螺旋線形式向上流動(dòng)，在下飛輪、屏蔽套流體間隙中沿徑向速度梯度較大。得出的速度分布及數(shù)值可為屏蔽電機(jī)冷卻計(jì)算及設(shè)計(jì)時(shí)平均表面對(duì)流傳熱系數(shù)計(jì)算提供參考數(shù)據(jù)。此外，可為流場、溫度場耦合計(jì)算提供檢驗(yàn)或比較基準(zhǔn)。

[1] 李良.AP1000核反應(yīng)堆用冷卻劑泵水力模型的設(shè)計(jì)與研究[D].大連：大連理工大學(xué)，2010.

[2] 姚若萍，饒芳權(quán)．蒸發(fā)冷卻水輪發(fā)電機(jī)定子溫度場研究[J].中國電機(jī)工程學(xué)報(bào)，2003，23(6)：87-90.

[3] 路義萍，洪光宇，湯璐，韓家德.多風(fēng)路大型空冷汽輪發(fā)電機(jī)三維流場計(jì)算[J].中國電機(jī)工程學(xué)報(bào)，2013，33(3)：133-139.

[4] KuosaM，Sallinen P，Larjola J.Numerical and experimental modelling of gas flow and heat transfer in the air gap of an electric machine[J].Journal of Thermal Science，2004，13(3)：264-278.

[5] 孫桂瑛，鄒繼斌，李連鋒，李夢啟，索文旭，李春廷.屏蔽電機(jī)的磁場與屏蔽層損耗的解析計(jì)算[J].哈爾濱工業(yè)大學(xué)學(xué)報(bào)，1999，31(4)：30-33.

[6] 馮穎慧.AP100020屏蔽電機(jī)主泵屏蔽套力學(xué)分析[D].上海交通大學(xué)，11.

[7] 丁樹業(yè)，孟繁東，葛云中.核主泵屏蔽電機(jī)溫度場研究[J].中國電機(jī)工程學(xué)報(bào)，2012，32(36)：149-155.

[8] 成德，薛亞波，沈洪，姚振強(qiáng)，呂向平，李藏雪.屏蔽式電動(dòng)泵內(nèi)部流場的分析與可視化[J].核科學(xué)與工程，2013，33(4)：398-403.

[9] 陶文銓.數(shù)值傳熱學(xué)[M].西安：西安交通大學(xué)出版社，2001：1-6.

[10] 王福軍.計(jì)算流體動(dòng)力學(xué)分析[M]．北京：清華大學(xué)出版社，2004：1-23，69-131.

Research of 3D Flow Field of the Water route of Canned Primary Pump Motors in Nuclear Power Stations

LI Cang-xue1，ZHAO Bo-min2，LU Yi-ping2，LV Xiang-ping1

(1.Harbin Electric Machinery Company Limited，Harbin of Heilongjiang Prov.150066，China；2.School of Mechanical & Power Engineering，Harbin University of Science and Technology，Harbin of Heilongjiang Prov. 150080，China)

The canned primary pump motor is an important part of a loop system of the nuclear island in nuclear power stations，while the flow states of the cooling water inside the canned motor directly affect the heat transfer within the motor. Consequently，a three dimensional model of flow field of the whole motor was established. Based on the finite volume method in the computational fluid dynamics theory，the velocity distribution in canned motor was calculated under boundary condition of the velocity inlet of primary water at the rated working point using the Ansys Fluent software . The axial and radial velocity distribution inside the canned motor was obtained by respectively analyzing the flow field of key parts of the canned motor. The results show that，near the rotating wall surface，the water with high velocity flow upward along the axial direction in the form of the spiral line. In the flywheel and the shielding clearance，the fluid velocity along the radial velocity gradient is larger. The research can provide the velocity magnitude and distribution in different locations for the followed temperature field investigation.

The canned motor；Three dimensional flow field；Finite volume method；Computational Fluid Dynamics

2016-05-22

國家科技重大專項(xiàng)(2013ZX06002002)

李藏雪(1966—)，女，學(xué)士，高級(jí)工程師，從事多種電機(jī)優(yōu)化設(shè)計(jì)及多場特性研究

TK121

A

0258-0918(2016)04-0533-06