微重力下航天器艙內(nèi)混合對(duì)流傳熱高性能數(shù)值研究

夏京京，張 超，王 嫻

微重力下航天器艙內(nèi)混合對(duì)流傳熱高性能數(shù)值研究

夏京京，張 超，王 嫻

（西安交通大學(xué)航天航空學(xué)院機(jī)械強(qiáng)度與振動(dòng)國(guó)家重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710049）

基于CPUs/GPUs架構(gòu)，采用混合熱格子-Boltzmann方法，對(duì)微重力環(huán)境下三維航天器艙簡(jiǎn)化模型內(nèi)的混合對(duì)流問題進(jìn)行了詳細(xì)的非穩(wěn)態(tài)數(shù)值模擬。結(jié)果表明：具有恒熱流邊界的元器件位于進(jìn)風(fēng)掠過(guò)的模型壁面中心時(shí)，熱壁面平均努塞爾數(shù)最大，換熱最強(qiáng)；保持雷諾數(shù)Re不變而縮小模型尺寸到1:5及以下時(shí)，模型在全重力情況下與原型在無(wú)重力時(shí)的熱壁面相同，流場(chǎng)和溫度場(chǎng)瞬態(tài)變化過(guò)程幾乎一致，可用來(lái)指導(dǎo)在地面進(jìn)行微重力下的實(shí)驗(yàn)研究；多GPUs結(jié)合格子Boltzmann方法計(jì)算的超高效率，可實(shí)現(xiàn)航天器電力系統(tǒng)和設(shè)備艙內(nèi)流動(dòng)換熱的實(shí)時(shí)預(yù)測(cè)。

航天器艙；微重力；混合對(duì)流；混合熱格子-Boltzmann法；GPUs

空間飛行器的發(fā)展對(duì)推進(jìn)技術(shù)提出了更高要求。相比傳統(tǒng)的化學(xué)推進(jìn)等技術(shù)，核推進(jìn)、太陽(yáng)能電推進(jìn)等新型推進(jìn)技術(shù)具備高能量密度，更適合為長(zhǎng)期駐留太空任務(wù)供電。艙內(nèi)電力系統(tǒng)是航天器的核心組成部分，因此組成電力系統(tǒng)的電子元器件的安全可靠性對(duì)航天器的安全運(yùn)行具有決定性作用。由于長(zhǎng)期服役于微重力環(huán)境，航天器艙內(nèi)自然對(duì)流效應(yīng)被削弱，熱均勻性下降，易造成局部熱流密度過(guò)大，從而損壞元器件，進(jìn)一步威脅航天器和推進(jìn)系統(tǒng)的安全。因此，研發(fā)艙內(nèi)熱管理技術(shù)意義重大。

然而，航天器在軌實(shí)驗(yàn)可能性小，利用傳統(tǒng)的落井等裝置創(chuàng)造的微重力時(shí)間短，無(wú)法滿足熱穩(wěn)定要求，實(shí)驗(yàn)研究存在困難。因此，利用數(shù)值模擬技術(shù)深入研究微重力環(huán)境下航天器艙內(nèi)對(duì)流換熱，對(duì)航天器一體化設(shè)計(jì)和指導(dǎo)在地面進(jìn)行微重力下的實(shí)驗(yàn)研究具有重要意義。

但以上研究涉及的通風(fēng)方式單一，且艙體內(nèi)部結(jié)構(gòu)簡(jiǎn)單，未考慮艙體內(nèi)置元器件，不夠接近真實(shí)物理問題。而針對(duì)腔體含熱源的研究頗多：Jami等人[15]對(duì)含導(dǎo)熱熱源方腔的自然對(duì)流問題進(jìn)行了研究；Gangawane等人[16-17]對(duì)內(nèi)置固體為等溫和等熱流2種邊界展開了討論。但基于微重力下的此類問題研究較少。此外，以往研究由于方法及硬件條件限制，采用的網(wǎng)格分辨率較低，且二維穩(wěn)態(tài)計(jì)算居多，對(duì)于三維流動(dòng)傳熱的非穩(wěn)態(tài)特性缺乏足夠的研究。

隨著計(jì)算流體動(dòng)力學(xué)（CFD）規(guī)模的增大，計(jì)算量呈非線性增加，刺激網(wǎng)格計(jì)算技術(shù)[18]和并行計(jì)算大力發(fā)展。硬件方面，GPU逐漸被用于非顯示的通用計(jì)算GPGPU（general purpose graphic process unit）。LBM具有天然并行性，CPUs并行效率近100%，其算法亦與具有單指令多線程并行模式的GPU完美匹配，但GPU在性能和功耗上均比CPU構(gòu)建的系統(tǒng)更具優(yōu)勢(shì)，故采用CPU+GPU異構(gòu)方式構(gòu)建的超級(jí)計(jì)算機(jī)已經(jīng)成為主流。研究表明，基于LBM的單GPU加速比可超過(guò)100倍[19-22]，在單GPU上運(yùn)行單精度LBM代碼，每秒可達(dá)10億網(wǎng)格更新[23]。此外，采用類似于CPU/GPU異構(gòu)架構(gòu)+CUDA化+MPI庫(kù)+LBM方法實(shí)現(xiàn)高性能計(jì)算（HPC）已成為開發(fā)高性能計(jì)算程序的流行趨勢(shì)。其中，Wang和Aoki[24]進(jìn)行了超大規(guī)模的三維流動(dòng)計(jì)算，上官燕琴等[25-26]已實(shí)現(xiàn)上億網(wǎng)格量的大規(guī)模湍流計(jì)算。

基于此，本文在CPU/GPU異構(gòu)集群上，采用混合熱格子Boltzmann方法（HTLBM，即LBM-FDM），基于CUDA C-MPI自行開發(fā)代碼，對(duì)航天器艙簡(jiǎn)化模型內(nèi)混合對(duì)流問題進(jìn)行實(shí)時(shí)動(dòng)態(tài)預(yù)測(cè)，旨在得到三維精細(xì)的流動(dòng)傳熱結(jié)果，為航天器一體化設(shè)計(jì)和在地面進(jìn)行微重力環(huán)境下的實(shí)驗(yàn)研究提供一定參考，并推動(dòng)GPU高性能計(jì)算在工程領(lǐng)域的應(yīng)用。

1 計(jì)算模型

圖1為航天器艙簡(jiǎn)化后的計(jì)算模型，艙體模型長(zhǎng)寬高比為::=1:1:1，艙體上壁面進(jìn)風(fēng)速度（c）均勻，進(jìn)口溫度恒定（c=300 K）。出風(fēng)口在腔體 兩側(cè)壁下端，且為高度的1/12，兩側(cè)壁等溫（h= 310 K），其余壁面絕熱。由于出口采用充分發(fā)展邊界，故所有模型在實(shí)際計(jì)算中出口均做延長(zhǎng)處理。

圖1 航天器艙幾何模型示意

2 數(shù)值方法

為了節(jié)省顯存消耗，采用HTLBM進(jìn)行模擬，即用等溫LBM模型求解速度場(chǎng)，用有限差分方法求解對(duì)流-擴(kuò)散方程從而得到溫度場(chǎng)結(jié)果。浮升力項(xiàng)采用Boussinesq假設(shè)[27]。硬件上采用CPU/GPU架構(gòu)集群（4節(jié)點(diǎn)，Tesla K20M GPU共12顆）。

2.1 流場(chǎng)計(jì)算

流場(chǎng)的模擬采用D3Q19模型（圖2），在LBE演化方程右端增加作用力項(xiàng)，以體現(xiàn)浮升力：

式中：fi和分別為粒子分布函數(shù)和平衡態(tài)分布函數(shù)，ei為粒子離散速度，d t為時(shí)間步長(zhǎng)，t 為無(wú)量綱馳豫時(shí)間，F(xiàn)i為離散速度空間的外力項(xiàng)。

外力項(xiàng)F采用Buick-Greated模型[27]

式中：浮力項(xiàng)=g(–c)，為流體密度，為重力加速度，為熱膨脹系數(shù)，為流體溫度，c為進(jìn)口流體溫度，為權(quán)系數(shù)，s為格子聲速。

通過(guò)宏觀N–S方程，可得到黏性系數(shù)與弛豫時(shí)間之間存在如下關(guān)系：

上述演化方程的求解需通過(guò)碰撞和遷移2個(gè)過(guò)程來(lái)完成（假設(shè)t=1）：

碰撞步：

遷移步：

為恢復(fù)相應(yīng)的宏觀方程，平衡態(tài)分布函數(shù)的選取至關(guān)重要。本文采取如下平衡態(tài)分布函數(shù)：

D3Q19的離散速度分別為0(0, 0, 0)、1(1, 0, 0)、2(–1, 0, 0)、3(0, 1, 0)、4(0, –1, 0)、5(0, 0, 1)、6(0, 0, –1)、7(1, 1, 0)、8(–1, 1, 0)、9(1, –1, 0)、10(–1,–1, 0)、11(1, 0, 0)、12(–1, 0, 1)、13(1, 0, –1)、14(–1, 0, –1)、15(0, 1, 1)、16(0, –1, 1)、17(0, 1, –1)、18(0, –1, –1)。權(quán)系數(shù)具體取值為0=1/3，1~6=1/18，7~8=1/36。

Buick-Greated模型[28]中宏觀速度?。?/p>

此外，流場(chǎng)計(jì)算中，邊界均采用非平衡外推格式[28]進(jìn)行處理。

2.2 溫度場(chǎng)計(jì)算

由于流動(dòng)的不可壓縮性，且忽略黏性效應(yīng)的影響，溫度場(chǎng)可進(jìn)行獨(dú)立求解。將其采用有限差分格式離散可得

式中：為流體溫度；ΔFD為有限差分的時(shí)間步長(zhǎng)；為熱擴(kuò)散系數(shù)，=/c，其中和p分別為流體的導(dǎo)熱系數(shù)和比定壓熱容；腳標(biāo)f代表流體。對(duì)流擴(kuò)散方程的對(duì)流項(xiàng)和擴(kuò)散項(xiàng)分別采用QUICK和二階中心差分格式。

為了實(shí)現(xiàn)與LB方程的協(xié)同，采用無(wú)量綱化對(duì)流擴(kuò)散方程

其中無(wú)量綱變量定義如下：

為保證數(shù)值計(jì)算的穩(wěn)定性，時(shí)間步長(zhǎng)ΔFD應(yīng)滿足：

CFL條件限制：

擴(kuò)散項(xiàng)限制條件：

由此確定最終時(shí)間步長(zhǎng)應(yīng)滿足

取進(jìn)口速度為特征速度：c=in，腔體高度為特征長(zhǎng)度，相關(guān)無(wú)量綱參數(shù)定義如下：

此外，本文將以以下參數(shù)作為流動(dòng)換熱性能的評(píng)價(jià)指標(biāo)：

2）局部努塞爾數(shù)1：

3）速度不均勻系數(shù)v：

其中v越大，表征速度越不均勻。

4）溫度不均勻系數(shù)T：

其中T越大，表征溫度分布越不均勻。

3 CUDA-MPI并行程序設(shè)計(jì)與開發(fā)

采用多GPUs開展并行計(jì)算，網(wǎng)格分辨率為192×192×192，采用CUDA C針對(duì)GPU編寫程序，采用MPI實(shí)現(xiàn)跨節(jié)點(diǎn)GPU間數(shù)據(jù)傳輸。

3.1 GPU一維分割設(shè)置

對(duì)三維區(qū)域采用方向上的一維分割。根據(jù)以往研究[24]，當(dāng)GPU數(shù)量較少時(shí)，計(jì)算區(qū)域的一維分割比多維分割編程更容易，且并行效率并未降低。本文采用12顆GPU并行，計(jì)算區(qū)域在方向進(jìn)行分塊。關(guān)于數(shù)據(jù)傳遞的具體細(xì)節(jié)詳見文獻(xiàn)[24]。此外，計(jì)算采用1個(gè)時(shí)間步長(zhǎng)內(nèi)，先遷移后碰撞的順序，使GPU計(jì)算時(shí)減少1次讀寫顯存的過(guò)程，計(jì)算效率比先碰撞后遷移的順序提高30%~40%。計(jì)算采用708萬(wàn)網(wǎng)格，512 000個(gè)時(shí)間步長(zhǎng)的非穩(wěn)態(tài)計(jì)算時(shí)間約為2 970 s，其中每個(gè)步長(zhǎng)中包含流動(dòng)和傳熱的計(jì)算，計(jì)算性能約為1 220.7 MLUPS。

3.2 程序驗(yàn)證

為驗(yàn)證本文所開發(fā)的HTLBM-GPUs程序的準(zhǔn)確性，利用該程序?qū)θS頂蓋驅(qū)動(dòng)混合熱對(duì)流問題進(jìn)行模擬，物理模型如圖3所示。上壁面為高溫壁面h，下壁面為低溫壁面c，其他壁面絕熱，所有壁面均是無(wú)滑移固體邊界。對(duì)=1.0（表示自然對(duì)流與強(qiáng)迫對(duì)流的相對(duì)強(qiáng)度），=100、=400、=1 000共3種工況進(jìn)行計(jì)算。

圖3 三維頂蓋驅(qū)動(dòng)混合熱對(duì)流幾何模型

圖4為本文數(shù)值計(jì)算得到不同工況下的流線圖與文獻(xiàn)[29]的對(duì)比，其中，文獻(xiàn)[29]采用MAC（Marker and Cell）方法求解N-S方程，網(wǎng)格數(shù)為81×81×81。

圖4 本文與文獻(xiàn)[29]計(jì)算得到的不同工況下的流線圖

由圖4可見：隨著的增大，主渦的尺寸逐漸減小，而二級(jí)渦的渦核高度逐漸增大；本文結(jié)果與文獻(xiàn)[29]的計(jì)算結(jié)果基本一致。此外，本文的數(shù)值計(jì)算結(jié)果表明，當(dāng)增大到400與1 000時(shí)，在下壁面附近還生成了尺寸較小的三級(jí)渦，而文 獻(xiàn)[29]的結(jié)果并未給出三級(jí)渦，可能與其所使用的網(wǎng)格分辨率較低有關(guān)。

Tab.1 The average Nusselt number of the hot wall obtained by the present work and reference [29-30]

4 結(jié)果與討論

4.1 發(fā)熱元器件位置對(duì)流動(dòng)換熱影響

圖6為=100、=0下發(fā)熱元器件位于艙體不同位置時(shí)的等溫線圖、溫度等值面圖和流線圖。其中，除元器件位于前（后）壁面中心時(shí)為/=1/16處的截面，其余均為展向中截面/=1/2。由圖6a)和圖6b)可看到：發(fā)熱元器件周圍流體的溫度最高，進(jìn)口處溫度最低，且除元器件位于艙體右（左）壁面中心外，其他幾種情況下的溫度場(chǎng)分布均關(guān)于/=1/2對(duì)稱。此外，由圖6c)可觀察到，當(dāng)元器件位于艙體正中心和上壁面中心時(shí)，在上壁面附近存在一對(duì)反對(duì)稱渦，但位于上壁面中心的元器件阻礙了這對(duì)反對(duì)稱渦的生長(zhǎng)，故此時(shí)的反對(duì)稱渦的尺寸較小。當(dāng)元器件位于艙體下壁面中心時(shí)，有一對(duì)尺寸較小且渦強(qiáng)度較弱的反對(duì)稱渦也形成于上壁面附近。這主要是因?yàn)?，元器件距離進(jìn)風(fēng)口較遠(yuǎn)，對(duì)冷氣射流的阻塞效應(yīng)較弱，使得這對(duì)反向旋轉(zhuǎn)渦的尺寸及強(qiáng)度都較小。而當(dāng)元器件位于右（左）壁面中心時(shí)，與其位于上壁面中心時(shí)在相同的位置產(chǎn)生一對(duì)反向旋轉(zhuǎn)渦，但由于障礙物的存在，渦尺寸較小，且由于元器件的位置改變了流場(chǎng)的左右對(duì)稱性，使上壁面附近反向旋轉(zhuǎn)渦分布不對(duì)稱。當(dāng)元器件位于艙體前（后）壁面中心時(shí)，流線分布特性與元器件位于艙體正中心時(shí)很相似。而且，由于位于艙體前（后）壁面中心的元器件對(duì)冷氣射流的阻塞效應(yīng)較強(qiáng)，所以反對(duì)稱渦的渦強(qiáng)及尺寸比元器件位于艙體正中心的大。

圖5 發(fā)熱元器件位于航天器艙內(nèi)不同位置的模型

Tab.2 The , and of the hot wall with the heating component located at different positions

4.2 模型比例對(duì)流動(dòng)換熱影響

航天器服役于微重力環(huán)境，在地面實(shí)驗(yàn)條件下，由于重力的作用會(huì)產(chǎn)生自然對(duì)流，故需予以消除。本小節(jié)采用通過(guò)縮小模型尺寸、增大速度來(lái)保持不變的方法，指導(dǎo)圖1計(jì)算模型在地面實(shí)現(xiàn)微重力環(huán)境下的實(shí)驗(yàn)研究。

Tab.3 The average Nusselt number of hot wall under different model scales

圖7展示了不同R下的側(cè)壁面局部1分布，其中/為沿側(cè)壁高度方向的無(wú)量綱距離。當(dāng)R=1:1時(shí)，自然對(duì)流作用很強(qiáng)，換熱能力和無(wú)重力時(shí)相差很大；當(dāng)R變小時(shí)，ln1分布逐漸趨于無(wú)重力狀態(tài)，而R≤1:5之后，無(wú)重力和有重力的側(cè)壁面ln1分布幾乎重合。

圖7 不同模型比例下熱壁面局部分布

圖8和圖9分別為=100時(shí)，不同R下展向中截面（/=1/2）上的等溫線圖及流線圖。

由圖8可以看出：當(dāng)R=1:1時(shí)，等溫線出現(xiàn)“S”型反轉(zhuǎn)。此時(shí)熱傳遞方式為對(duì)流傳熱，且自然對(duì)流作用占主導(dǎo)。對(duì)應(yīng)地，該模型比例下艙體內(nèi)出現(xiàn)一對(duì)尺寸大且渦強(qiáng)大的反向旋轉(zhuǎn)渦。隨著R減小，等溫線逐漸呈現(xiàn)近似平行于熱壁面的一簇曲線。此時(shí)混合對(duì)流中強(qiáng)迫對(duì)流占主導(dǎo)地位，自然對(duì)流作用減弱。

相應(yīng)地，由圖9可看到艙體內(nèi)的反向旋轉(zhuǎn)渦尺寸隨著R的減小而減小，并最終消失；而當(dāng)R=1:3時(shí)，在艙體上壁面附近開始形成一對(duì)新的反向旋轉(zhuǎn)渦，這是由于進(jìn)口射流兩側(cè)壓強(qiáng)減小，在頂部形成負(fù)壓區(qū)，進(jìn)而卷吸周圍流體形成旋渦；此外，隨著R減小，這對(duì)小旋渦逐漸向中部靠攏，表明自然對(duì)流作用進(jìn)一步減小，上升氣流對(duì)射流的壓迫作用逐漸消失。對(duì)比原模型在無(wú)重力時(shí)的等溫線及流線，發(fā)現(xiàn)在R≤1:5時(shí)，二者幾乎完全相同。

圖10為模型艙內(nèi)流動(dòng)換熱瞬態(tài)結(jié)果。從換熱演變過(guò)程可以看到，有重力、R=1:5的流動(dòng)換熱瞬態(tài)結(jié)果與無(wú)重力、原模型的結(jié)果幾乎完全一致。熱量從左、右兩側(cè)的熱壁面向流體傳遞，使得中間的流體溫度升高。

圖8 不同模型比下展向中截面(y/L=1/2)的等溫線圖

圖9 不同模型比例下展向中截面(y/L=1/2)的流線圖

圖10 模型艙內(nèi)流動(dòng)換熱瞬態(tài)結(jié)果

5 結(jié) 論

1）微重力場(chǎng)中，當(dāng)發(fā)熱元器件位于進(jìn)風(fēng)掠過(guò)的壁面中心時(shí)，最利于散熱；2）對(duì)于所研究模型，當(dāng)模型尺寸縮小至1:5及以下時(shí)，自然對(duì)流作用被抑制，可利用該方法在地面模擬微重力環(huán)境；3）多GPUs結(jié)合LBM計(jì)算的超高效率可實(shí)現(xiàn)航天器電力系統(tǒng)和設(shè)備艙內(nèi)流動(dòng)換熱的實(shí)時(shí)預(yù)測(cè)。

［1］ 李勁東, 邵東兵. 航天器密封艙微重力對(duì)流換熱的數(shù)值分析[J]. 工程熱物理學(xué)報(bào), 2002, 23(3): 333-335. LI Jingdong, SHAO Dongbing. Numeric analysis on a micro-gravity heat convection for a spacecraft pressured cabin[J]. Journal of Engineering Thermophysics, 2002, 23(3): 333-335.

［2］ 劉云龍. 空間站艙內(nèi)流動(dòng)和傳熱的地面模擬及數(shù)值分析[D]. 北京: 清華大學(xué), 1997: 20-29. LIU Yunlong. Ground simulation and numerical analysis of flow and heat transfer in a space station cabin[D]. Beijing: Tsinghua University, 1997: 20-29.

［3］ 劉云龍, 任建勛. 抑制自然對(duì)流的數(shù)調(diào)整技術(shù)[J]. 清華大學(xué)學(xué)報(bào)(自然科學(xué)版), 1998(10): 116-118. LIU Yunlong, REN Jianxun. Grashof number adjustment technique to suppress natural convection[J]. Journal of Tsinghua University (Science and Technology), 1998(10): 116-118.

［4］ 任建勛, 梁新剛, 陳澤敬, 等. 抑制自然對(duì)流的可視化實(shí)驗(yàn)研究[J]. 清華大學(xué)學(xué)報(bào)(自然科學(xué)版), 1998(8): 57. REN Jianxun, LIANG Xingang, CHEN Zejing, et al. Vsualization experimentresearch on gravity elimination technique[J]. Journal of Tsinghua University (Science and Technology), 1998(8): 57.

［5］ 馬蘭新, 譚建宇. 航天器密封艙降壓法抑制自然對(duì)流的數(shù)值研究[J]. 節(jié)能技術(shù), 2013, 31(6): 487-490. MA Lanxin, TAN Jianyu. Numerical simulation of a spacecraft pressurized cabin using lowering pressure method[J]. Energy Conservation Technology, 2013, 31(6): 487-490.

［6］ 張學(xué)學(xué), 劉靜, 蒙毅. 載人飛船座艙內(nèi)空氣對(duì)流換熱地面模擬[J]. 清華大學(xué)學(xué)報(bào)(自然科學(xué)版), 1999, 39(6): 67-70. ZHANG Xuexue, LIU Jing, MENG Yi. Ground-test simulating of air heat convection in manned space capsule[J]. Journal of Tsinghua University (Science and Technology), 1999, 39(6): 67-70.

［7］ 姬朝玥, 任建勛, 梁新剛, 等. 微重力下對(duì)流換熱的地面縮比—減壓模擬技術(shù)[J]. 清華大學(xué)學(xué)報(bào)(自然科學(xué)版), 2001, 41(2): 92-94. JI Chaoyue, REN Jianxun, LIANG Xingang, et al. Ground-based scale and pressure reduction technology for simulation of convection in micro-gravity[J]. Journal of Tsinghua University (Science and Technology), 2001, 41(2): 92-94.

［8］ 吳宏, 李娜, 梁新剛, 等. 空間站實(shí)驗(yàn)艙進(jìn)出風(fēng)口大小影響的數(shù)值研究[J]. 航天醫(yī)學(xué)與醫(yī)學(xué)工程, 2001, 14(4): 268-271. WU Hong, LI Na, LIANG Xingang, et al. Numerical study of the effects of size of inlets and outlets of ventilation system in experimental space station cabin[J]. Space Medicine & Medical Engineering, 2001, 14(4): 268-271.

［9］ 鄭忠海, 張吉禮, 梁珍. 空間站艙內(nèi)側(cè)上送側(cè)下回對(duì)稱通風(fēng)方式數(shù)值模擬[J]. 哈爾濱工業(yè)大學(xué)學(xué)報(bào), 2007, 39(2): 270-273. ZHENG Zhonghai, ZHANG Jili, LIANG Zhen. Numerical simulation of upside air-supply bottom-side air-return symmetric ventilation in space station cabin[J]. Journal of Harbin Institute of Technology, 2007, 39(2): 270-273.

［10］ 吳群剛, 梁新剛, 任建勛, 等. 方形空間通風(fēng)換熱的氣流組織性能數(shù)值分析[J]. 航天醫(yī)學(xué)與醫(yī)學(xué)工程, 2000, 13(3): 174-178. WU Qungang, LIANG Xingang, REN Jianxun, et al. Numerical study of ventilation and heat transfer performance in a square space[J]. Space Medicine & Medical Engineering, 2000, 13(3): 174-178.

［11］ KAREEM A K, GAO S. Mixed convection heat transfer of turbulent flow in a three-dimensional lid-driven cavity with a rotating cylinder[J]. International Journal of Heat & Mass Transfer, 2017, 112: 185-200.

［12］ VERGARA J B, CUESTA I, SALUENA C. Numerical study of laminar mixed convection in a square open cavity[J]. International Journal of Heat & Mass Transfer, 2016, 99: 599-612.

［13］ KARIMIPOUR A, NEZHAD A H, D’ORAZIO A, et al. Investigation of the gravity effects on the mixed convection heat transfer in a microchannel using lattice Boltzmann method[J]. International Journal of Thermal Sciences, 2012, 54(19): 142-152.

［14］ 陽(yáng)海成, 楊帆, 郭雪巖. 方腔內(nèi)自然對(duì)流與混合對(duì)流的有限差分格子Boltzmann研究[J]. 力學(xué)季刊, 2016(1): 33-44. YANG Haicheng, YANG Fan, GUO Xueyan. Finite difference lattice Boltzmann study of natural convection and mixed convection in a square cavity[J]. Chinese Quarterly of Mechanics, 2016(1): 33-44.

［15］ JAMI M, MEZRHAB A, BOUZIDI M, et al. Lattice Boltzmann method applied to the laminar natural convection in an enclosure with a heat-generating cylinder conducting body[J]. International Journal of Thermal Sciences, 2007, 46(1): 38-47.

［16］ GANGAWANE K M, MANIKANDAN B. Mixed convection characteristics in lid-driven cavity containing heated triangular block[J]. Chinese Journal of Chemical Engineering, 2017, 25(10): 1381-1394.

［17］ GANGAWANE K M. Computational analysis of mixed convection heat transfer characteristics in lid-driven cavity containing triangular block with constant heat flux: effect of prandtl and grashof numbers[J]. International Journal of Heat & Mass Transfer, 2017, 105: 34-57.

［18］ 蔡巧言, 杜濤. 網(wǎng)格計(jì)算技術(shù)在大規(guī)模CFD計(jì)算上的應(yīng)用[J]. 空氣動(dòng)力學(xué)學(xué)報(bào), 2009, 27(6): 696-700. CAI Qiaoyan, DU Tao. Application of grid computing technology to large-scale CFD computing[J]. Acta Aerodynamica Sinica, 2009, 27(6): 696-700.

［19］ JACOBSEN D A, THIBAULT J C, SENOCAK I. An MPI-CUDA implementation for massively parallel incompressible flow computations on multi-GPU clusters[M]. 48th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition. 2010: 1-16.

［20］ WEI X, ZHAO Y, FAN Z, et al. Lattice-based flow field modeling[J]. IEEE Transactions on Visualization & Computer Graphics, 2004, 10(6): 719-729.

［21］ FAN Z, QIU F, KAUFMAN A, et al. GPU cluster for high performance computing[C]. Supercomputing, ACM/IEEE Science Conference. IEEE, 2004.

［22］ CHANG H W, HONG P Y, LIN L S, et al. Simulations of three-dimensional cavity flows with multi relaxation time lattice Boltzmann method and graphic processing units[J]. Procedia Engineering, 2013, 61: 94-99.

［23］ KUZNIK F, OBRECHT C, RUSAOUEN G, et al. LBM based flow simulation using GPU computing processor[J]. Computers & Mathematics with Applications, 2010, 59(7): 2380-2392.

［24］ WANG X, AOKI T. Multi-GPU performance of incompressible flow computation by lattice Boltzmann method on GPU cluster[J]. Parallel Computing, 2011, 37(9): 521.

［25］ SHANGGUAN Y, WANG X, LI Y. Investigation on the mixing mechanism of single-jet film cooling with various blowing ratios based on hybrid thermal lattice Boltzmann method[J]. International Journal of Heat & Mass Transfer, 2016, 97: 880-890.

［26］ SHANGGUAN Y, WANG X, ZHANG H, et al. Analysis on the mechanism of evolutionary process of counter-rotating vortex pair in film cooling based on hybrid thermal lattice Boltzmann method[J]. International Communications in Heat & Mass Transfer, 2017, 87: 72.

［27］ BUICK J M, GREATED C A. Gravity in a lattice Boltzmann model[J]. Physical Review E Statistical Physics Plasmas Fluids & Related Interdisciplinary Topics, 2000, 61(5A): 5307-5319.

［28］ GUO Z L, ZHENG C G, SHI B C. Non-equilibrium extrapolation method for velocity and pressure boundary conditions in the lattice Boltzmann method[J]. Chinese Physics B, 2002, 11(4): 366-374.

［29］ IWATSU R, HYUN J M. Three-dimensional driven-cavity flows with a vertical temperature gradient[J]. International Journal of Heat & Mass Transfer, 1995, 38(18): 3319-3328.

［30］ DU H Y, CHAI Z H, SHI B C. Lattice Boltzmann study of mixed convection in a cubic cavity[J]. Communications in Theoretical Physics, 2011, 56: 144-150.

High-performance numerical study on mixed convection and heat transfer in spacecraft cabin under microgravity

XIA Jingjing, ZHANG Chao, WANG Xian

(State Key Laboratory for Strength and Vibration of Mechanical Structures, School of Aerospace, Xi’an Jiaotong University, Xi’an 710049, China)

On the basis of architecture of CPUs/GPUs and hybrid thermal lattice Boltzmann method, unsteady state numerical simulations on mixed convection and heat transfer in three-dimensional simplified model of the spacecraft cabin under microgravity are performed. The results show that, when the heating component with constant heat flow boundary is located in the center of the model wall which is swept by the inlet air, the average Nusselt number of the hot wall reaches the largest and the heat transfer is the strongest. When the model scale is or below 1:5 and the Reynolds number Re is kept unchanged, the average Nusselt number of the hot wall, and the transient evolution process of the flow and temperature field of the scaled model under gravity environment is the same as that of the prototype under microgravity environment, which can be used to guide the experimental study under microgravity on the ground. With ultra high efficiency, the multiple GPUs combined with LBM can realize real-time prediction of flow and heat transfer in electric system and equipment cabin of the spacecraft.

spacecraft cabin, microgravity, mixed convection, hybrid thermal lattice Boltzmann method, GPUs

TK124

A

10.19666/j.rlfd.201904040

夏京京, 張超, 王嫻. 微重力下航天器艙內(nèi)混合對(duì)流傳熱高性能數(shù)值研究[J]. 熱力發(fā)電, 2019, 48(8): 27-35. XIA Jingjing, ZHANG Chao, WANG Xian. High-performance numerical study on mixed convection and heat transfer in spacecraft cabin under microgravity[J]. Thermal Power Generation, 2019, 48(8): 27-35.

2019-04-01

面向E級(jí)計(jì)算機(jī)的大型流體機(jī)械并行計(jì)算軟件系統(tǒng)及示范(2016YFB0200901)

Supported by：Parallel Computing Software System and Demonstration of Large Fluid Machinery for E-Class Computer (2016YFB0200901)

夏京京(1991—)，女，碩士研究生，主要研究方向?yàn)槲⒅亓ο屡搩?nèi)熱管理及其高性能程序開發(fā)，xjj20080808@stu.xjtu.edu.cn。

王嫻(1977—)，女，博士，副教授，主要研究方向?yàn)榱鲃?dòng)傳熱高性能計(jì)算，wangxian@mail.xjtu.edu.cn。

（責(zé)任編輯 李園）