流速對(duì)魚(yú)骨式與板式調(diào)溫?zé)岢列阅艿挠绊?/h1>
2022-11-01 11:21:40劉智勇顏昌林祁雪琴
蘭州交通大學(xué)學(xué)報(bào) 2022年5期
馮 寧，劉智勇*，顏昌林，祁雪琴
（1.蘭州交通大學(xué)環(huán)境與市政工程學(xué)院，蘭州 730070；2.蘭州真空設(shè)備有限責(zé)任公司華宇分公司，蘭州 730013；3.中科院近代物理研究所，蘭州 730099）
為了保證航天器研制和發(fā)射的安全，有必要在地面開(kāi)展空間環(huán)模試驗(yàn).在熱真空試驗(yàn)箱中，能夠模擬高、低溫交替變換的空間環(huán)境設(shè)備稱(chēng)為熱沉.具有溫度調(diào)節(jié)功能的熱沉被稱(chēng)為調(diào)溫?zé)岢?，即操縱可控的載冷劑溫度或流量于熱沉中，使熱沉溫度在一定范圍內(nèi)變化，從而模擬出溫度變化的外界環(huán)境［1］.
目前常見(jiàn)的調(diào)溫?zé)岢两Y(jié)構(gòu)有魚(yú)骨式和板式兩類(lèi)［2-5］.魚(yú)骨式熱沉是采用不銹鋼管作為載冷劑介質(zhì)管路，焊接在不銹鋼管上的紫銅當(dāng)作散熱片，以充分發(fā)揮不銹鋼良好的低溫耐腐蝕性能、紫銅良好的導(dǎo)熱性能，所組成的熱沉［6-7］.而板式熱沉則是將兩個(gè)不銹鋼板疊加，用激光焊接成具有矩形樣式的整體板，再通過(guò)壓力膨脹處理而成［8-9］.圖1和圖2分別是常見(jiàn)的魚(yú)骨式熱沉結(jié)構(gòu)與板式熱沉結(jié)構(gòu)的實(shí)物圖.文獻(xiàn)［10-12］在計(jì)算機(jī)上利用成熟的有限元分析軟件（ANSYS）對(duì)影響?hù)~(yú)骨式熱沉結(jié)構(gòu)均勻性的主要結(jié)構(gòu)因素（速度、溫度、壓力）進(jìn)行了數(shù)值模擬分析.徐照武等［13］對(duì)魚(yú)骨式熱沉進(jìn)行了溫度調(diào)節(jié)系統(tǒng)的設(shè)計(jì)、試驗(yàn)研究以及關(guān)鍵參數(shù)分析計(jì)算，并對(duì)系統(tǒng)進(jìn)行了調(diào)試.文獻(xiàn)［14-16］對(duì)影響新型板式熱沉結(jié)構(gòu)氣氮熱沉均勻性的重要結(jié)構(gòu)因素（壁板間距、氮?dú)鈮毫Α⒉煌M(jìn)、出口布置形式、流速等）進(jìn)行速度、溫度、壓力數(shù)值模擬分析.魏龍濤等［17］研制了一種新型密閉循環(huán)的回?zé)崾綒獾{(diào)溫系統(tǒng)，并采用數(shù)值模擬與試驗(yàn)研究相結(jié)合的方法，通過(guò)改變循環(huán)氮?dú)鈮毫Γ瑢?duì)此調(diào)溫系統(tǒng)的溫度場(chǎng)、壓力場(chǎng)、升、降溫過(guò)程進(jìn)行研究.目前國(guó)內(nèi)外大多數(shù)的研究還只是針對(duì)環(huán)境模擬器魚(yú)骨式調(diào)溫?zé)岢恋臄?shù)值模擬，在設(shè)定熱沉?xí)r往往根據(jù)經(jīng)驗(yàn)和仿真軟件來(lái)確定，對(duì)其試驗(yàn)的驗(yàn)證相對(duì)較少，而對(duì)于新型調(diào)溫?zé)岢翜囟染鶆蛐缘膶?shí)驗(yàn)研究更少.因此對(duì)空間環(huán)境模擬設(shè)備熱沉均勻性的研究還有很大的發(fā)展空間.
圖1 魚(yú)骨式熱沉結(jié)構(gòu)Fig．1 Fishbone heat sink structure
圖2 板式熱沉結(jié)構(gòu)Fig．2 Plate heat sink structure
本文研究的目的就是通過(guò)數(shù)值模擬的方法從循環(huán)氮?dú)膺M(jìn)口速度及升、降溫速率來(lái)比較分析兩種熱沉結(jié)構(gòu)的性能優(yōu)劣，提出可行的建議，以提高熱沉溫度分布的均勻性.在熱沉溫度均勻性滿(mǎn)足指標(biāo)要求的前提下，對(duì)不同熱沉結(jié)構(gòu)的動(dòng)態(tài)響應(yīng)性能進(jìn)行了分析和評(píng)價(jià)，為今后的熱沉結(jié)構(gòu)設(shè)計(jì)提供理論基礎(chǔ).
1 熱沉的幾何模型與數(shù)值計(jì)算
1.1 熱沉的幾何模型
兩種熱沉的規(guī)格尺寸均以現(xiàn)場(chǎng)實(shí)驗(yàn)裝置為計(jì)算分析模型.魚(yú)骨式結(jié)構(gòu)如圖3（a）所示，框架主要由主管與支管焊接形成，材質(zhì)均為不銹鋼，框架的內(nèi)外兩側(cè)皆焊接有銅翅片，以吸收熱輻射，熱沉大小為Φ1.0 m×1.5 m，主管內(nèi)徑100 mm，支管內(nèi)徑16 mm，翅片寬度為80 mm，翅片厚度為2 mm.下方主管為進(jìn)液管，上方主管為出液管，為更容易使熱沉溫度達(dá)到均勻的狀態(tài)，采用載冷劑下進(jìn)上出的方式，這樣可以起到避免產(chǎn)生空氣阻力的作用.將兩塊鋼板模壓出凹坑再點(diǎn)焊連接，形成夾層結(jié)構(gòu).板式熱沉結(jié)構(gòu)如圖3（b）所示，空間直徑1 000 mm、長(zhǎng)1 500 mm，不銹鋼板厚度為2 mm.不銹鋼板的間距為10 mm，焊點(diǎn)直徑為12 mm，焊點(diǎn)正常分布間距為60 mm.熱沉整體分成兩部分，每部分用折流板隔成3部分并呈現(xiàn)“蛇形”形狀的長(zhǎng)流道，每個(gè)流道的進(jìn)、出口設(shè)有一個(gè)半徑為20 mm的進(jìn)液支管與出液支管，分別與進(jìn)液主管和出液主管相通，下側(cè)為進(jìn)液管，上側(cè)為出液管，且進(jìn)、出液口設(shè)置在同一側(cè).
圖3 仿真模型Fig．3 Simulation model
1.2 熱沉的數(shù)值計(jì)算
本文對(duì)建模時(shí)需要載冷劑換熱的數(shù)學(xué)模型進(jìn)行以下幾種假設(shè)：
1）若熱真空試驗(yàn)箱中的壓力在10-3Pa以下，可忽略對(duì)流換熱；2）載冷劑是一種不可壓縮流體；3）忽略載冷劑密度差引起的浮力；4）由于載冷劑在熱沉管網(wǎng)中的流速較低，因此可以忽略粘性耗散產(chǎn)生的熱效應(yīng)；5）載冷劑的物性參數(shù)保持不變，呈氣態(tài)；6）忽略熱沉各翅片間的輻射換熱.
翅片管中載冷劑地流動(dòng)和傳熱的守恒型控制方程描述如下：
1）連續(xù)性方程：
式中：u，v，w為 x，y，z方向上速度矢量的分量．
2）直角坐標(biāo)系中i方向上的動(dòng)量方程：
式中：ρ為流體密度；p為流體壓力；ν為流體的運(yùn)動(dòng)黏度；Ui為i方向的速度分量．
3）能量守恒方程：
式中：a為熱擴(kuò)散率；T為流體的溫度．
本文介紹的兩種熱沉中載冷劑的流態(tài)均為湍流，湍流模型采用RNG模型中的κ-ε方程，其計(jì)算公式如下：
式中：Gk為平均速度的湍流動(dòng)能；μeff為黏度系數(shù)；C1ε、C2ε為模型常量，分別為 1．42、1．68；αk、αε為 κ方程和ε方程的湍流Prandtl數(shù)；Sk、Sε為源項(xiàng)，根據(jù)具體條件定義．
本論文使用FLUENT計(jì)算模塊，以有限容積方法為基礎(chǔ)，通過(guò)數(shù)值模擬方法來(lái)求解問(wèn)題.而離散格式是一種離散形式，它是數(shù)值計(jì)算中微分和偏微分導(dǎo)數(shù)的離散化方法，本文使用的松弛因子如表1所列.
表1 離散方程中采用的松弛因子Tab．1 Relaxation factors used in discrete equations
1.3 邊界條件
用Ansys Fluent 19.0仿真計(jì)算時(shí)，采用速度入口邊界條件、自由出流邊界條件、壁面邊界條件，中心平面設(shè)置為symmetry.在低溫調(diào)節(jié)過(guò)程中，入口處氮?dú)馑俣仍O(shè)置為14.227 m／s，氮?dú)鉁囟仍O(shè)置為143 K，熱沉外壁表面僅接收來(lái)自真空容器內(nèi)壁面的輻射，溫度設(shè)置為298 K，外發(fā)射率設(shè)定為0.2，負(fù)載時(shí)熱沉內(nèi)表面也僅接受被試件的輻射，溫度設(shè)定為298 K，內(nèi)發(fā)射率設(shè)定為0.9，入口壓力為60 kPa.湍流模型使用κ-ε標(biāo)準(zhǔn)模型，壓力-速度耦合采用 Simple算法［18］.保證板式熱沉與魚(yú)骨式熱沉邊界條件一致.
1.4 網(wǎng)格劃分及模型驗(yàn)證
1.4.1 網(wǎng)格劃分
本次模擬利用ICEM CFD對(duì)魚(yú)骨式熱沉及板式熱沉整體結(jié)構(gòu)進(jìn)行了網(wǎng)格劃分，熱沉模型網(wǎng)格劃分結(jié)果如圖4所示.為了避免網(wǎng)格密度對(duì)仿真結(jié)果造成影響，減少對(duì)計(jì)算機(jī)精度的要求，縮短計(jì)算時(shí)間，在正式仿真前需要?jiǎng)澐侄嘟M不同密度的網(wǎng)格.本文采用了3組不同的網(wǎng)格（3 982 992、4 701 048、5 534 122）進(jìn)行傳熱模擬，模擬結(jié)果如圖5所示.對(duì)計(jì)算結(jié)果進(jìn)行比較分析，從3組網(wǎng)格獲得的最大相對(duì)誤差僅為0.015%，達(dá)到了工程設(shè)計(jì)的精度要求.另外考慮到配備的計(jì)算機(jī)性能，最終選擇了4 701 048的網(wǎng)格數(shù)，以便進(jìn)行后續(xù)的計(jì)算.
圖4 熱沉溫度場(chǎng)網(wǎng)格劃分示意圖Fig．4 Schematic diagram of heat sink temperature field meshing
圖5 網(wǎng)格獨(dú)立性驗(yàn)證Fig．5 Verification of grid independence
1.4.2 模型驗(yàn)證
為驗(yàn)證數(shù)值模擬方法的準(zhǔn)確性，本文參照文獻(xiàn)［18］對(duì)調(diào)溫?zé)岢恋脑囼?yàn)研究進(jìn)行模型驗(yàn)證.其中氣氮調(diào)溫系統(tǒng)的工作原理為：系統(tǒng)以氮?dú)鉃槔鋮s介質(zhì)，將一定壓力和溫度的氮?dú)庵苯虞斔偷綗岢林?，并通過(guò)不斷調(diào)節(jié)氮?dú)獾臏囟葋?lái)控制熱沉的溫度，以實(shí)現(xiàn)熱沉溫度在一定范圍內(nèi)調(diào)節(jié)，進(jìn)而實(shí)現(xiàn)熱沉溫度連續(xù)可調(diào)，并間接對(duì)試驗(yàn)件的溫度進(jìn)行控制.圖6給出了調(diào)溫循環(huán)系統(tǒng)原理圖.試驗(yàn)臺(tái)由熱沉、回?zé)崞?、加熱器、混合器、循環(huán)風(fēng)機(jī)、液氮儲(chǔ)槽等組成.循環(huán)過(guò)程如下：將氣氮裝入系統(tǒng)，進(jìn)行閉式循環(huán)，高溫循環(huán)時(shí)，接通加熱器，循環(huán)氣氮由加熱器進(jìn)行加熱處理，高溫氣氮進(jìn)入熱沉與熱沉外壁面交換熱量，使其溫度升高.低溫循環(huán)期間，高溫氣氮和低溫氣氮進(jìn)行充分的熱交換后降低氣氮的溫度，低溫氣氮進(jìn)入熱沉再與熱沉外壁面進(jìn)行熱交換降低熱沉溫度.最后將數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果進(jìn)行了比較.調(diào)溫?zé)岢翜囟葴y(cè)點(diǎn)設(shè)在熱沉筒體上，如圖7所示，一共安裝12個(gè)測(cè)溫點(diǎn)，均勻分布在3個(gè)回路上，每個(gè)回路中有4個(gè)測(cè)溫點(diǎn)，控制界面上以TR01～TR12進(jìn)行編號(hào)，試驗(yàn)環(huán)境溫度為298 K.
圖6 調(diào)溫?zé)岢料到y(tǒng)原理圖Fig．6 Schematic diagram of tem perature regulating heat sink system
圖7 溫度測(cè)點(diǎn)分布圖Fig．7 Distribution diagram of tem perature measurement points
將-130℃和+150℃平衡態(tài)時(shí)的數(shù)值模擬結(jié)果與試驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比，驗(yàn)證了該數(shù)值模擬過(guò)程的準(zhǔn)確性.圖8和圖9分別為-130℃、+150℃平衡態(tài)時(shí)試驗(yàn)結(jié)果與模擬結(jié)果的比較.通過(guò)比較分析發(fā)現(xiàn)，平衡態(tài)時(shí)熱沉外壁面試驗(yàn)與模擬溫度趨勢(shì)基本吻合，溫度相差最大為2.04 K、2.399 K.結(jié)果表明，可以用該方法模擬熱沉溫度場(chǎng).
圖8 －130℃平衡態(tài)試驗(yàn)與模擬結(jié)果對(duì)比Fig．8 Comparison between equilibrium test and simulation results at－130℃
圖9 ＋150℃平衡態(tài)試驗(yàn)與模擬結(jié)果對(duì)比Fig．9 Comparison between equilibrium test and simulation results at＋150℃
2 循環(huán)氮?dú)膺M(jìn)口速度對(duì)熱沉溫度均勻性的影響
模擬空間環(huán)境中的溫度變化，可以使用將溫度或流量可控制的冷劑傳遞到熱沉的方法將熱沉溫度控制在指定時(shí)間間隔內(nèi).本文保持熱沉載冷劑進(jìn)口溫度不變，分析氮?dú)膺M(jìn)口流量改變對(duì)不同熱沉結(jié)構(gòu)溫度均勻性的影響，其邊界條件同1.3.
圖10和圖11分別為不同流速時(shí)兩種熱沉結(jié)構(gòu)的外壁面溫度分布云圖.
圖10 不同流速時(shí)魚(yú)骨式熱沉外壁面溫度分布云圖Fig．10 Cloud diagram of tem perature distribution on outer wall of fishbone heat sink at different flow rates
圖11 不同流速時(shí)板式熱沉外壁面溫度分布云圖Fig．11 Cloud diagram of temperature distribution on outer wall of plate heat sink at different flow rates
從圖12可以看出，熱沉翅片外壁面溫度隨著進(jìn)口速度的增長(zhǎng)而有所下降，并且隨著流速的增加，兩種熱沉結(jié)構(gòu)溫度分布越均勻，熱沉壁面溫差也越小，溫度越接近控制點(diǎn)溫度.圖12（a）為在魚(yú)骨式熱沉中進(jìn)口速度由10 m／s增大到35 m／s時(shí)的熱沉外壁面溫度分布曲線(xiàn)，從圖中可以看出，熱沉翅片外壁面溫度隨進(jìn)口速度的增長(zhǎng)而有所下降，接近143 K，并且外壁面上側(cè)總是高于下側(cè).圖12（b）為在板式熱沉中不同流速下熱沉外壁面溫度分布曲線(xiàn)，不同的流速下，熱沉外壁溫度均勻性分別為±4.389 5℃、±3.144 5℃、±2.627 3℃，隨著流速的增加，熱沉分布越均勻，熱沉壁面溫差越小，且越接近控制點(diǎn)溫度.在相同熱沉尺寸、相同邊界條件、相同冷劑消耗且不考慮熱沉折流板尖角的情況下，板式熱沉與魚(yú)骨式熱沉相比，當(dāng)流速增加到15 m／s時(shí)，熱沉外壁面溫度已經(jīng)接近控制點(diǎn)溫度，不需要再提高流速，且溫度均勻度比魚(yú)骨式熱沉更好，流量消耗更少.板式熱沉具有更高的傳熱效率和更好的溫度均勻性，這是由于載冷劑與熱沉的內(nèi)表面直接接觸進(jìn)行熱交換，并且熱沉夾層中有許多小的圓柱體，可以起到擾流的作用.
圖12 不同進(jìn)口流速下熱沉外壁面溫度的分布曲線(xiàn)Fig．12 Tem perature distribution curve of outer wall of heat sink at different in let flow rates
3 熱沉溫度場(chǎng)降溫過(guò)程模擬計(jì)算
本次模擬的是兩種熱沉結(jié)構(gòu)在+150～-130℃范圍內(nèi)的降溫速率.進(jìn)行非穩(wěn)態(tài)計(jì)算時(shí)，設(shè)置氮?dú)馊肟跍囟葹?30 K，入口流速為14.227 m／s，熱沉外壁面僅進(jìn)行輻射換熱，溫度設(shè)置為423 K，外部發(fā)射率設(shè)置為0.2，熱沉的初始環(huán)境為423 K，計(jì)算步長(zhǎng)為6 s，計(jì)算步數(shù)為500，計(jì)算結(jié)果如圖13和圖14所示.
圖13 和圖14為熱沉外壁面溫度非穩(wěn)態(tài)降溫過(guò)程.從魚(yú)骨式熱沉結(jié)構(gòu)云圖看出，模擬開(kāi)始后，入口附近的翅片表面溫度開(kāi)始降低.當(dāng)仿真時(shí)間達(dá)到3 000 s時(shí)，熱沉翅片表面溫度基本不發(fā)生變化，翅片表面溫度分布在156 K到168 K之間，溫度分布較之前更加均勻.從板式熱沉結(jié)構(gòu)云圖看出，流體首先流入第1個(gè)進(jìn)口，第1段熱沉溫度開(kāi)始下降，熱沉第2段溫度依舊很高，隨著流體在進(jìn)口支管內(nèi)流動(dòng)，流體流入第2段熱沉，第2段熱沉開(kāi)始變化，第1段熱沉變化快于第2段熱沉.隨著時(shí)間的增加，熱沉溫度逐漸下降到設(shè)定溫度，趨向于均勻，且在折流板附近逐漸聚集，出現(xiàn)溫度過(guò)高的現(xiàn)象.
圖13 魚(yú)骨式熱沉外壁面溫度隨時(shí)間變化云圖Fig．13 Cloud diagram of temperature variation w ith time on outer wall of fishbone heat sink
圖14 板式熱沉外壁面溫度隨時(shí)間變化云圖Fig．14 Cloud diagram of tem perature variation w ith time on outer wall of plate heat sink
翅片外壁面溫度隨時(shí)間變化曲線(xiàn)如圖15所示.從圖15可看出，熱沉外壁面的溫度隨著載冷劑的流入開(kāi)始下降，而當(dāng)熱沉剛開(kāi)始工作時(shí)，初始時(shí)刻的載冷劑和熱沉之間存在較大的溫度差.在1～15 min時(shí)，魚(yú)骨式熱沉結(jié)構(gòu)的降溫速率為9.84 K／min；15～35 min時(shí)，熱沉的降溫速率達(dá)到3.41 K／min；當(dāng)時(shí)間增加至60 min時(shí)，熱沉的降溫速率僅為1.07 K／min，此時(shí)熱沉溫度場(chǎng)已經(jīng)達(dá)到了穩(wěn)定狀態(tài)，溫度穩(wěn)定在143 K左右.板式熱沉第1段表面的溫度下降很明顯，第1段熱沉溫度平均下降速率為25.328 K／min，而第2段熱沉溫度平均下降速率為16.276 K／min；隨著時(shí)間的增加，熱沉溫度下降速率也在減緩，在5～15 min時(shí)，熱沉溫度下降的平均速率為10.548 K／min；而30 min后，熱沉外壁面的溫度完全穩(wěn)定.相同進(jìn)口流量下，板式熱沉溫度下降速率遠(yuǎn)大于魚(yú)骨式熱沉，板式熱沉在第30 min就能降至設(shè)定溫度.而魚(yú)骨式熱沉雖然30 min后降溫速率滿(mǎn)足指標(biāo)要求，但溫度并未達(dá)到設(shè)定溫度，降溫還在繼續(xù).相對(duì)來(lái)講，板式熱沉消耗的載冷劑量更少，消耗的時(shí)間也更少.
圖15 熱沉外壁面溫度隨時(shí)間變化曲線(xiàn)Fig．15 Curve of tem perature w ith time on outer wall of heat sink
4 結(jié)論
通過(guò)以上的數(shù)值模擬分析及試驗(yàn)研究得出以下結(jié)論：
1）板式結(jié)構(gòu)熱沉氮?dú)獾牧魉購(gòu)?0 m／s增加到25 m／s的過(guò)程中，隨著流速的增加其溫度分布越均勻，均勻性從4.39℃減小至2.63℃，接近控制點(diǎn)溫度；魚(yú)骨式結(jié)構(gòu)熱沉氮?dú)饬魉購(gòu)?0 m／s增加到35 m／s的過(guò)程中，溫度均勻性從1.71℃減小至0.97℃，二者均滿(mǎn)足溫度均勻性的設(shè)計(jì)要求.
2）在板式熱沉中，流速增加到15 m／s時(shí)，熱沉外壁面溫度已經(jīng)在設(shè)定點(diǎn)溫度附近，而魚(yú)骨式熱沉流速要到35 m／s時(shí)才到達(dá)設(shè)定點(diǎn)溫度附近，板式熱沉結(jié)構(gòu)相比于魚(yú)骨式結(jié)構(gòu)，不需要再提高流速，系統(tǒng)所需氮?dú)獾牧髁枯^少.
3）板式熱沉溫度下降速率較大，在相同熱沉尺寸、相同邊界條件、相同冷劑消耗且不考慮熱沉折流板尖角的情況下，板式熱沉前30 min的降溫速率為9.33 K／min，30 min后熱沉外壁溫度完全穩(wěn)定；魚(yú)骨式熱沉結(jié)構(gòu)的降溫速率較慢，魚(yú)骨式熱沉前30 min的降溫速率為6.03 K／min，60 min后溫度場(chǎng)達(dá)到穩(wěn)定狀態(tài).
猜你喜歡
天然氣液化工廠冷劑儲(chǔ)存區(qū)的布置淺析云南化工(2021年7期)2021-12-21 07:27:50 滑蓋產(chǎn)品的三板式注塑模具設(shè)計(jì)裝備制造技術(shù)(2020年12期)2020-05-22 09:25:06 奶奶愛(ài)拼魚(yú)骨畫(huà)保健醫(yī)苑(2019年2期)2019-04-03 06:28:52 板式吊耳設(shè)計(jì)中的相關(guān)計(jì)算石油化工建設(shè)(2018年4期)2018-11-30 02:04:02 基于混料設(shè)計(jì)的MRC工藝中混合冷劑配比的優(yōu)化天然氣化工—C1化學(xué)與化工(2018年5期)2018-11-15 01:29:52 CRTSⅡ型板式無(wú)砟軌道板間接縫開(kāi)裂修補(bǔ)技術(shù)的研究進(jìn)展江西建材(2018年4期)2018-04-10 12:36:48 一星期沒(méi)換水的夢(mèng)境詩(shī)潮(2018年3期)2018-03-26 12:29:30 淺析吸收式制冷機(jī)組冷劑水故障問(wèn)題西部皮革(2018年2期)2018-02-17 06:51:01 魚(yú)骨千萬(wàn)別丟 它能幫你增壽工友(2016年4期)2016-09-18 05:57:49 某天然氣液化工藝?yán)湎浞e液分析及排除①石油與天然氣化工(2014年4期)2014-09-11 05:54:44
馮 寧，劉智勇*，顏昌林，祁雪琴
（1.蘭州交通大學(xué)環(huán)境與市政工程學(xué)院，蘭州 730070；2.蘭州真空設(shè)備有限責(zé)任公司華宇分公司，蘭州 730013；3.中科院近代物理研究所，蘭州 730099）
為了保證航天器研制和發(fā)射的安全，有必要在地面開(kāi)展空間環(huán)模試驗(yàn).在熱真空試驗(yàn)箱中，能夠模擬高、低溫交替變換的空間環(huán)境設(shè)備稱(chēng)為熱沉.具有溫度調(diào)節(jié)功能的熱沉被稱(chēng)為調(diào)溫?zé)岢?，即操縱可控的載冷劑溫度或流量于熱沉中，使熱沉溫度在一定范圍內(nèi)變化，從而模擬出溫度變化的外界環(huán)境［1］.
目前常見(jiàn)的調(diào)溫?zé)岢两Y(jié)構(gòu)有魚(yú)骨式和板式兩類(lèi)［2-5］.魚(yú)骨式熱沉是采用不銹鋼管作為載冷劑介質(zhì)管路，焊接在不銹鋼管上的紫銅當(dāng)作散熱片，以充分發(fā)揮不銹鋼良好的低溫耐腐蝕性能、紫銅良好的導(dǎo)熱性能，所組成的熱沉［6-7］.而板式熱沉則是將兩個(gè)不銹鋼板疊加，用激光焊接成具有矩形樣式的整體板，再通過(guò)壓力膨脹處理而成［8-9］.圖1和圖2分別是常見(jiàn)的魚(yú)骨式熱沉結(jié)構(gòu)與板式熱沉結(jié)構(gòu)的實(shí)物圖.文獻(xiàn)［10-12］在計(jì)算機(jī)上利用成熟的有限元分析軟件（ANSYS）對(duì)影響?hù)~(yú)骨式熱沉結(jié)構(gòu)均勻性的主要結(jié)構(gòu)因素（速度、溫度、壓力）進(jìn)行了數(shù)值模擬分析.徐照武等［13］對(duì)魚(yú)骨式熱沉進(jìn)行了溫度調(diào)節(jié)系統(tǒng)的設(shè)計(jì)、試驗(yàn)研究以及關(guān)鍵參數(shù)分析計(jì)算，并對(duì)系統(tǒng)進(jìn)行了調(diào)試.文獻(xiàn)［14-16］對(duì)影響新型板式熱沉結(jié)構(gòu)氣氮熱沉均勻性的重要結(jié)構(gòu)因素（壁板間距、氮?dú)鈮毫Α⒉煌M(jìn)、出口布置形式、流速等）進(jìn)行速度、溫度、壓力數(shù)值模擬分析.魏龍濤等［17］研制了一種新型密閉循環(huán)的回?zé)崾綒獾{(diào)溫系統(tǒng)，并采用數(shù)值模擬與試驗(yàn)研究相結(jié)合的方法，通過(guò)改變循環(huán)氮?dú)鈮毫Γ瑢?duì)此調(diào)溫系統(tǒng)的溫度場(chǎng)、壓力場(chǎng)、升、降溫過(guò)程進(jìn)行研究.目前國(guó)內(nèi)外大多數(shù)的研究還只是針對(duì)環(huán)境模擬器魚(yú)骨式調(diào)溫?zé)岢恋臄?shù)值模擬，在設(shè)定熱沉?xí)r往往根據(jù)經(jīng)驗(yàn)和仿真軟件來(lái)確定，對(duì)其試驗(yàn)的驗(yàn)證相對(duì)較少，而對(duì)于新型調(diào)溫?zé)岢翜囟染鶆蛐缘膶?shí)驗(yàn)研究更少.因此對(duì)空間環(huán)境模擬設(shè)備熱沉均勻性的研究還有很大的發(fā)展空間.
圖1 魚(yú)骨式熱沉結(jié)構(gòu)Fig．1 Fishbone heat sink structure
圖2 板式熱沉結(jié)構(gòu)Fig．2 Plate heat sink structure
本文研究的目的就是通過(guò)數(shù)值模擬的方法從循環(huán)氮?dú)膺M(jìn)口速度及升、降溫速率來(lái)比較分析兩種熱沉結(jié)構(gòu)的性能優(yōu)劣，提出可行的建議，以提高熱沉溫度分布的均勻性.在熱沉溫度均勻性滿(mǎn)足指標(biāo)要求的前提下，對(duì)不同熱沉結(jié)構(gòu)的動(dòng)態(tài)響應(yīng)性能進(jìn)行了分析和評(píng)價(jià)，為今后的熱沉結(jié)構(gòu)設(shè)計(jì)提供理論基礎(chǔ).
1 熱沉的幾何模型與數(shù)值計(jì)算
1.1 熱沉的幾何模型
兩種熱沉的規(guī)格尺寸均以現(xiàn)場(chǎng)實(shí)驗(yàn)裝置為計(jì)算分析模型.魚(yú)骨式結(jié)構(gòu)如圖3（a）所示，框架主要由主管與支管焊接形成，材質(zhì)均為不銹鋼，框架的內(nèi)外兩側(cè)皆焊接有銅翅片，以吸收熱輻射，熱沉大小為Φ1.0 m×1.5 m，主管內(nèi)徑100 mm，支管內(nèi)徑16 mm，翅片寬度為80 mm，翅片厚度為2 mm.下方主管為進(jìn)液管，上方主管為出液管，為更容易使熱沉溫度達(dá)到均勻的狀態(tài)，采用載冷劑下進(jìn)上出的方式，這樣可以起到避免產(chǎn)生空氣阻力的作用.將兩塊鋼板模壓出凹坑再點(diǎn)焊連接，形成夾層結(jié)構(gòu).板式熱沉結(jié)構(gòu)如圖3（b）所示，空間直徑1 000 mm、長(zhǎng)1 500 mm，不銹鋼板厚度為2 mm.不銹鋼板的間距為10 mm，焊點(diǎn)直徑為12 mm，焊點(diǎn)正常分布間距為60 mm.熱沉整體分成兩部分，每部分用折流板隔成3部分并呈現(xiàn)“蛇形”形狀的長(zhǎng)流道，每個(gè)流道的進(jìn)、出口設(shè)有一個(gè)半徑為20 mm的進(jìn)液支管與出液支管，分別與進(jìn)液主管和出液主管相通，下側(cè)為進(jìn)液管，上側(cè)為出液管，且進(jìn)、出液口設(shè)置在同一側(cè).
圖3 仿真模型Fig．3 Simulation model
1.2 熱沉的數(shù)值計(jì)算
本文對(duì)建模時(shí)需要載冷劑換熱的數(shù)學(xué)模型進(jìn)行以下幾種假設(shè)：
1）若熱真空試驗(yàn)箱中的壓力在10-3Pa以下，可忽略對(duì)流換熱；2）載冷劑是一種不可壓縮流體；3）忽略載冷劑密度差引起的浮力；4）由于載冷劑在熱沉管網(wǎng)中的流速較低，因此可以忽略粘性耗散產(chǎn)生的熱效應(yīng)；5）載冷劑的物性參數(shù)保持不變，呈氣態(tài)；6）忽略熱沉各翅片間的輻射換熱.
翅片管中載冷劑地流動(dòng)和傳熱的守恒型控制方程描述如下：
1）連續(xù)性方程：
式中：u，v，w為 x，y，z方向上速度矢量的分量．
2）直角坐標(biāo)系中i方向上的動(dòng)量方程：
式中：ρ為流體密度；p為流體壓力；ν為流體的運(yùn)動(dòng)黏度；Ui為i方向的速度分量．
3）能量守恒方程：
式中：a為熱擴(kuò)散率；T為流體的溫度．
本文介紹的兩種熱沉中載冷劑的流態(tài)均為湍流，湍流模型采用RNG模型中的κ-ε方程，其計(jì)算公式如下：
式中：Gk為平均速度的湍流動(dòng)能；μeff為黏度系數(shù)；C1ε、C2ε為模型常量，分別為 1．42、1．68；αk、αε為 κ方程和ε方程的湍流Prandtl數(shù)；Sk、Sε為源項(xiàng)，根據(jù)具體條件定義．
本論文使用FLUENT計(jì)算模塊，以有限容積方法為基礎(chǔ)，通過(guò)數(shù)值模擬方法來(lái)求解問(wèn)題.而離散格式是一種離散形式，它是數(shù)值計(jì)算中微分和偏微分導(dǎo)數(shù)的離散化方法，本文使用的松弛因子如表1所列.
表1 離散方程中采用的松弛因子Tab．1 Relaxation factors used in discrete equations
1.3 邊界條件
用Ansys Fluent 19.0仿真計(jì)算時(shí)，采用速度入口邊界條件、自由出流邊界條件、壁面邊界條件，中心平面設(shè)置為symmetry.在低溫調(diào)節(jié)過(guò)程中，入口處氮?dú)馑俣仍O(shè)置為14.227 m／s，氮?dú)鉁囟仍O(shè)置為143 K，熱沉外壁表面僅接收來(lái)自真空容器內(nèi)壁面的輻射，溫度設(shè)置為298 K，外發(fā)射率設(shè)定為0.2，負(fù)載時(shí)熱沉內(nèi)表面也僅接受被試件的輻射，溫度設(shè)定為298 K，內(nèi)發(fā)射率設(shè)定為0.9，入口壓力為60 kPa.湍流模型使用κ-ε標(biāo)準(zhǔn)模型，壓力-速度耦合采用 Simple算法［18］.保證板式熱沉與魚(yú)骨式熱沉邊界條件一致.
1.4 網(wǎng)格劃分及模型驗(yàn)證
1.4.1 網(wǎng)格劃分
本次模擬利用ICEM CFD對(duì)魚(yú)骨式熱沉及板式熱沉整體結(jié)構(gòu)進(jìn)行了網(wǎng)格劃分，熱沉模型網(wǎng)格劃分結(jié)果如圖4所示.為了避免網(wǎng)格密度對(duì)仿真結(jié)果造成影響，減少對(duì)計(jì)算機(jī)精度的要求，縮短計(jì)算時(shí)間，在正式仿真前需要?jiǎng)澐侄嘟M不同密度的網(wǎng)格.本文采用了3組不同的網(wǎng)格（3 982 992、4 701 048、5 534 122）進(jìn)行傳熱模擬，模擬結(jié)果如圖5所示.對(duì)計(jì)算結(jié)果進(jìn)行比較分析，從3組網(wǎng)格獲得的最大相對(duì)誤差僅為0.015%，達(dá)到了工程設(shè)計(jì)的精度要求.另外考慮到配備的計(jì)算機(jī)性能，最終選擇了4 701 048的網(wǎng)格數(shù)，以便進(jìn)行后續(xù)的計(jì)算.
圖4 熱沉溫度場(chǎng)網(wǎng)格劃分示意圖Fig．4 Schematic diagram of heat sink temperature field meshing
圖5 網(wǎng)格獨(dú)立性驗(yàn)證Fig．5 Verification of grid independence
1.4.2 模型驗(yàn)證
為驗(yàn)證數(shù)值模擬方法的準(zhǔn)確性，本文參照文獻(xiàn)［18］對(duì)調(diào)溫?zé)岢恋脑囼?yàn)研究進(jìn)行模型驗(yàn)證.其中氣氮調(diào)溫系統(tǒng)的工作原理為：系統(tǒng)以氮?dú)鉃槔鋮s介質(zhì)，將一定壓力和溫度的氮?dú)庵苯虞斔偷綗岢林?，并通過(guò)不斷調(diào)節(jié)氮?dú)獾臏囟葋?lái)控制熱沉的溫度，以實(shí)現(xiàn)熱沉溫度在一定范圍內(nèi)調(diào)節(jié)，進(jìn)而實(shí)現(xiàn)熱沉溫度連續(xù)可調(diào)，并間接對(duì)試驗(yàn)件的溫度進(jìn)行控制.圖6給出了調(diào)溫循環(huán)系統(tǒng)原理圖.試驗(yàn)臺(tái)由熱沉、回?zé)崞?、加熱器、混合器、循環(huán)風(fēng)機(jī)、液氮儲(chǔ)槽等組成.循環(huán)過(guò)程如下：將氣氮裝入系統(tǒng)，進(jìn)行閉式循環(huán)，高溫循環(huán)時(shí)，接通加熱器，循環(huán)氣氮由加熱器進(jìn)行加熱處理，高溫氣氮進(jìn)入熱沉與熱沉外壁面交換熱量，使其溫度升高.低溫循環(huán)期間，高溫氣氮和低溫氣氮進(jìn)行充分的熱交換后降低氣氮的溫度，低溫氣氮進(jìn)入熱沉再與熱沉外壁面進(jìn)行熱交換降低熱沉溫度.最后將數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果進(jìn)行了比較.調(diào)溫?zé)岢翜囟葴y(cè)點(diǎn)設(shè)在熱沉筒體上，如圖7所示，一共安裝12個(gè)測(cè)溫點(diǎn)，均勻分布在3個(gè)回路上，每個(gè)回路中有4個(gè)測(cè)溫點(diǎn)，控制界面上以TR01～TR12進(jìn)行編號(hào)，試驗(yàn)環(huán)境溫度為298 K.
圖6 調(diào)溫?zé)岢料到y(tǒng)原理圖Fig．6 Schematic diagram of tem perature regulating heat sink system
圖7 溫度測(cè)點(diǎn)分布圖Fig．7 Distribution diagram of tem perature measurement points
將-130℃和+150℃平衡態(tài)時(shí)的數(shù)值模擬結(jié)果與試驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比，驗(yàn)證了該數(shù)值模擬過(guò)程的準(zhǔn)確性.圖8和圖9分別為-130℃、+150℃平衡態(tài)時(shí)試驗(yàn)結(jié)果與模擬結(jié)果的比較.通過(guò)比較分析發(fā)現(xiàn)，平衡態(tài)時(shí)熱沉外壁面試驗(yàn)與模擬溫度趨勢(shì)基本吻合，溫度相差最大為2.04 K、2.399 K.結(jié)果表明，可以用該方法模擬熱沉溫度場(chǎng).
圖8 －130℃平衡態(tài)試驗(yàn)與模擬結(jié)果對(duì)比Fig．8 Comparison between equilibrium test and simulation results at－130℃
圖9 ＋150℃平衡態(tài)試驗(yàn)與模擬結(jié)果對(duì)比Fig．9 Comparison between equilibrium test and simulation results at＋150℃
2 循環(huán)氮?dú)膺M(jìn)口速度對(duì)熱沉溫度均勻性的影響
模擬空間環(huán)境中的溫度變化，可以使用將溫度或流量可控制的冷劑傳遞到熱沉的方法將熱沉溫度控制在指定時(shí)間間隔內(nèi).本文保持熱沉載冷劑進(jìn)口溫度不變，分析氮?dú)膺M(jìn)口流量改變對(duì)不同熱沉結(jié)構(gòu)溫度均勻性的影響，其邊界條件同1.3.
圖10和圖11分別為不同流速時(shí)兩種熱沉結(jié)構(gòu)的外壁面溫度分布云圖.
圖10 不同流速時(shí)魚(yú)骨式熱沉外壁面溫度分布云圖Fig．10 Cloud diagram of tem perature distribution on outer wall of fishbone heat sink at different flow rates
圖11 不同流速時(shí)板式熱沉外壁面溫度分布云圖Fig．11 Cloud diagram of temperature distribution on outer wall of plate heat sink at different flow rates
從圖12可以看出，熱沉翅片外壁面溫度隨著進(jìn)口速度的增長(zhǎng)而有所下降，并且隨著流速的增加，兩種熱沉結(jié)構(gòu)溫度分布越均勻，熱沉壁面溫差也越小，溫度越接近控制點(diǎn)溫度.圖12（a）為在魚(yú)骨式熱沉中進(jìn)口速度由10 m／s增大到35 m／s時(shí)的熱沉外壁面溫度分布曲線(xiàn)，從圖中可以看出，熱沉翅片外壁面溫度隨進(jìn)口速度的增長(zhǎng)而有所下降，接近143 K，并且外壁面上側(cè)總是高于下側(cè).圖12（b）為在板式熱沉中不同流速下熱沉外壁面溫度分布曲線(xiàn)，不同的流速下，熱沉外壁溫度均勻性分別為±4.389 5℃、±3.144 5℃、±2.627 3℃，隨著流速的增加，熱沉分布越均勻，熱沉壁面溫差越小，且越接近控制點(diǎn)溫度.在相同熱沉尺寸、相同邊界條件、相同冷劑消耗且不考慮熱沉折流板尖角的情況下，板式熱沉與魚(yú)骨式熱沉相比，當(dāng)流速增加到15 m／s時(shí)，熱沉外壁面溫度已經(jīng)接近控制點(diǎn)溫度，不需要再提高流速，且溫度均勻度比魚(yú)骨式熱沉更好，流量消耗更少.板式熱沉具有更高的傳熱效率和更好的溫度均勻性，這是由于載冷劑與熱沉的內(nèi)表面直接接觸進(jìn)行熱交換，并且熱沉夾層中有許多小的圓柱體，可以起到擾流的作用.
圖12 不同進(jìn)口流速下熱沉外壁面溫度的分布曲線(xiàn)Fig．12 Tem perature distribution curve of outer wall of heat sink at different in let flow rates
3 熱沉溫度場(chǎng)降溫過(guò)程模擬計(jì)算
本次模擬的是兩種熱沉結(jié)構(gòu)在+150～-130℃范圍內(nèi)的降溫速率.進(jìn)行非穩(wěn)態(tài)計(jì)算時(shí)，設(shè)置氮?dú)馊肟跍囟葹?30 K，入口流速為14.227 m／s，熱沉外壁面僅進(jìn)行輻射換熱，溫度設(shè)置為423 K，外部發(fā)射率設(shè)置為0.2，熱沉的初始環(huán)境為423 K，計(jì)算步長(zhǎng)為6 s，計(jì)算步數(shù)為500，計(jì)算結(jié)果如圖13和圖14所示.
圖13 和圖14為熱沉外壁面溫度非穩(wěn)態(tài)降溫過(guò)程.從魚(yú)骨式熱沉結(jié)構(gòu)云圖看出，模擬開(kāi)始后，入口附近的翅片表面溫度開(kāi)始降低.當(dāng)仿真時(shí)間達(dá)到3 000 s時(shí)，熱沉翅片表面溫度基本不發(fā)生變化，翅片表面溫度分布在156 K到168 K之間，溫度分布較之前更加均勻.從板式熱沉結(jié)構(gòu)云圖看出，流體首先流入第1個(gè)進(jìn)口，第1段熱沉溫度開(kāi)始下降，熱沉第2段溫度依舊很高，隨著流體在進(jìn)口支管內(nèi)流動(dòng)，流體流入第2段熱沉，第2段熱沉開(kāi)始變化，第1段熱沉變化快于第2段熱沉.隨著時(shí)間的增加，熱沉溫度逐漸下降到設(shè)定溫度，趨向于均勻，且在折流板附近逐漸聚集，出現(xiàn)溫度過(guò)高的現(xiàn)象.
圖13 魚(yú)骨式熱沉外壁面溫度隨時(shí)間變化云圖Fig．13 Cloud diagram of temperature variation w ith time on outer wall of fishbone heat sink
圖14 板式熱沉外壁面溫度隨時(shí)間變化云圖Fig．14 Cloud diagram of tem perature variation w ith time on outer wall of plate heat sink
翅片外壁面溫度隨時(shí)間變化曲線(xiàn)如圖15所示.從圖15可看出，熱沉外壁面的溫度隨著載冷劑的流入開(kāi)始下降，而當(dāng)熱沉剛開(kāi)始工作時(shí)，初始時(shí)刻的載冷劑和熱沉之間存在較大的溫度差.在1～15 min時(shí)，魚(yú)骨式熱沉結(jié)構(gòu)的降溫速率為9.84 K／min；15～35 min時(shí)，熱沉的降溫速率達(dá)到3.41 K／min；當(dāng)時(shí)間增加至60 min時(shí)，熱沉的降溫速率僅為1.07 K／min，此時(shí)熱沉溫度場(chǎng)已經(jīng)達(dá)到了穩(wěn)定狀態(tài)，溫度穩(wěn)定在143 K左右.板式熱沉第1段表面的溫度下降很明顯，第1段熱沉溫度平均下降速率為25.328 K／min，而第2段熱沉溫度平均下降速率為16.276 K／min；隨著時(shí)間的增加，熱沉溫度下降速率也在減緩，在5～15 min時(shí)，熱沉溫度下降的平均速率為10.548 K／min；而30 min后，熱沉外壁面的溫度完全穩(wěn)定.相同進(jìn)口流量下，板式熱沉溫度下降速率遠(yuǎn)大于魚(yú)骨式熱沉，板式熱沉在第30 min就能降至設(shè)定溫度.而魚(yú)骨式熱沉雖然30 min后降溫速率滿(mǎn)足指標(biāo)要求，但溫度并未達(dá)到設(shè)定溫度，降溫還在繼續(xù).相對(duì)來(lái)講，板式熱沉消耗的載冷劑量更少，消耗的時(shí)間也更少.
圖15 熱沉外壁面溫度隨時(shí)間變化曲線(xiàn)Fig．15 Curve of tem perature w ith time on outer wall of heat sink
4 結(jié)論
通過(guò)以上的數(shù)值模擬分析及試驗(yàn)研究得出以下結(jié)論：
1）板式結(jié)構(gòu)熱沉氮?dú)獾牧魉購(gòu)?0 m／s增加到25 m／s的過(guò)程中，隨著流速的增加其溫度分布越均勻，均勻性從4.39℃減小至2.63℃，接近控制點(diǎn)溫度；魚(yú)骨式結(jié)構(gòu)熱沉氮?dú)饬魉購(gòu)?0 m／s增加到35 m／s的過(guò)程中，溫度均勻性從1.71℃減小至0.97℃，二者均滿(mǎn)足溫度均勻性的設(shè)計(jì)要求.
2）在板式熱沉中，流速增加到15 m／s時(shí)，熱沉外壁面溫度已經(jīng)在設(shè)定點(diǎn)溫度附近，而魚(yú)骨式熱沉流速要到35 m／s時(shí)才到達(dá)設(shè)定點(diǎn)溫度附近，板式熱沉結(jié)構(gòu)相比于魚(yú)骨式結(jié)構(gòu)，不需要再提高流速，系統(tǒng)所需氮?dú)獾牧髁枯^少.
3）板式熱沉溫度下降速率較大，在相同熱沉尺寸、相同邊界條件、相同冷劑消耗且不考慮熱沉折流板尖角的情況下，板式熱沉前30 min的降溫速率為9.33 K／min，30 min后熱沉外壁溫度完全穩(wěn)定；魚(yú)骨式熱沉結(jié)構(gòu)的降溫速率較慢，魚(yú)骨式熱沉前30 min的降溫速率為6.03 K／min，60 min后溫度場(chǎng)達(dá)到穩(wěn)定狀態(tài).
猜你喜歡
云南化工(2021年7期)2021-12-21 07:27:50
裝備制造技術(shù)(2020年12期)2020-05-22 09:25:06
保健醫(yī)苑(2019年2期)2019-04-03 06:28:52
石油化工建設(shè)(2018年4期)2018-11-30 02:04:02
天然氣化工—C1化學(xué)與化工(2018年5期)2018-11-15 01:29:52
江西建材(2018年4期)2018-04-10 12:36:48
詩(shī)潮(2018年3期)2018-03-26 12:29:30
西部皮革(2018年2期)2018-02-17 06:51:01
工友(2016年4期)2016-09-18 05:57:49
石油與天然氣化工(2014年4期)2014-09-11 05:54:44
