短艙防冰系統(tǒng)三維內(nèi)外流耦合計算方法

鄭梅，馮麗娟，秦娜，尹金鴿

中國航發(fā)商用航空發(fā)動機有限責(zé)任公司，上海 200241

當(dāng)飛機穿越由過冷水滴組成的云層時，迎風(fēng)部件表面會發(fā)生結(jié)冰現(xiàn)象［1-3］。發(fā)動機短艙進氣道表面產(chǎn)生的冰聚集會改變發(fā)動機進氣流場和進氣流量，可能導(dǎo)致發(fā)動機推力下降、喘振或失速［4］；當(dāng)表面積冰發(fā)生破碎時，冰片可能被吸入發(fā)動機內(nèi)部，撞擊風(fēng)扇葉片，造成發(fā)動機機械損傷，嚴(yán)重危害飛行安全［5-7］。目前，在役的民用航空發(fā)動機普遍采用熱氣防冰結(jié)構(gòu)，其中笛形管熱氣防冰系統(tǒng)由于其結(jié)構(gòu)簡單、可靠性強、技術(shù)成熟度高等優(yōu)點，被廣泛應(yīng)用在發(fā)動機短艙進氣道上［8-9］。在短艙笛形管熱氣防冰系統(tǒng)設(shè)計過程中，防冰性能仿真分析是設(shè)計流程的一個重要環(huán)節(jié)，構(gòu)建有效的防冰計算方法可以大幅提升設(shè)計效率，為防冰性能評估及防冰結(jié)構(gòu)優(yōu)化設(shè)計提供支撐。

笛形管熱氣防冰系統(tǒng)的性能仿真涉及復(fù)雜的內(nèi)外流耦合傳熱［10-13］?；诜辣到y(tǒng)的運行環(huán)境與工作原理，防冰腔內(nèi)、外流動換熱特性差異顯著［14-18］，為仿真計算帶來較大難度。文獻中關(guān)于笛形管防冰系統(tǒng)的耦合仿真大多針對機翼防冰腔構(gòu)型［19-23］，所采用的耦合計算方法可分為強固耦合方法［24-25］和松散耦合方法［26-27］。強固耦合計算通常將防冰腔外流域、固體域及內(nèi)流域進行一體化網(wǎng)格劃分，在計算過程中同時求解內(nèi)外流場及固體域?qū)幔?8-29］；而松散耦合方法則通常將計算域進行拆分并選取合適的迭代策略以實現(xiàn)內(nèi)外流域的耦合計算。由于松散耦合并不要求內(nèi)外表面網(wǎng)格節(jié)點完全一致，因此在減少計算網(wǎng)格和提升計算效率上具有較大優(yōu)勢［30］。對于短艙進氣道笛形管防冰系統(tǒng)，由于全環(huán)模型尺寸大、笛形管上射流孔數(shù)多，且進氣道型面具有周向非對稱特性，無法參考機翼模型截取一段直翼防冰腔進行模擬。綜合計算網(wǎng)格與計算效率等因素，松散耦合方法更為適用。松散耦合迭代過程中涉及內(nèi)外流計算域耦合面的數(shù)據(jù)交互。相較于干空氣條件，濕空氣條件下防冰表面熱質(zhì)傳遞過程更為復(fù)雜，且進氣道表面還需考慮周向非對稱流場導(dǎo)致的三維溢流效應(yīng)對防冰表面熱載荷的影響。

基于此，針對短艙全尺寸模型的笛形管熱氣防冰系統(tǒng)，采用松散耦合方法開展了三維防冰內(nèi)外流仿真計算方法研究。將三維短艙進氣道計算模型中的外流場獨立為一個計算域，將蒙皮固體域與防冰腔內(nèi)流場組合為另一個計算域，以防冰腔蒙皮外表面作為兩個計算域的耦合面，通過耦合面溫度和換熱系數(shù)（Heat Transfer Coeffi?cient， HTC）相互賦值為邊界條件來實現(xiàn)內(nèi)外流計算域的數(shù)據(jù)交互。根據(jù)防冰計算流程，依次構(gòu)建干、濕空氣條件下的防冰內(nèi)外流松散耦合迭代策略，并在濕空氣條件下引入短艙周向非均勻性對防冰表面溢流流動換熱特性的影響。最終獲得迭代收斂后的短艙進氣道防冰腔內(nèi)外表面溫度及換熱系數(shù)分布。

1 三維防冰內(nèi)外流耦合計算方法

1.1 防冰計算流程

通常，防冰計算流程包括以下3個步驟：流場計算、水滴撞擊特性計算和防冰熱分析計算。

1.1.1 流場計算

采用ANSYS Fluent商用軟件進行流場計算，通過求解三維、定常、黏性流動的Navier-Stokes （N-S）方程，獲得防冰腔外部冷空氣流場、防冰腔內(nèi)部熱空氣流場以及蒙皮固壁溫度場。此時，為干空氣條件下的流場計算，防冰腔內(nèi)、外表面對流換熱系數(shù)均可通過CFD方法計算得到。其中，防冰腔外表面對流換熱系數(shù)以外流場環(huán)境溫度為參考溫度換算得到；而防冰腔內(nèi)表面對流換熱系數(shù)則以熱氣總溫作為參考溫度進行提取。

1.1.2 水滴撞擊特性計算

水滴撞擊特性計算采用歐拉法［31］，引入水滴相體積分?jǐn)?shù)，通過求解水滴相控制方程獲得水滴運動軌跡及防冰腔外表面局部水收集系數(shù)分布。

水滴相的連續(xù)性方程和動量方程可描述為

式中：α為水滴相體積分?jǐn)?shù)；ρd為水滴密度，由于水滴是不可壓縮的，因此式（1）中水滴密度項可約去；ud為水滴速度；CD為阻力系數(shù)；ρa為空氣密度；dd為水滴直徑；ua為空氣速度。

根據(jù)定義，表面局部水收集系數(shù)β的計算公式為

式中：n為表面撞擊位置的單位法向量；α∞與u∞分別為遠(yuǎn)場區(qū)域位置水滴相體積分?jǐn)?shù)和氣流速度。

將水滴相控制方程展開為Fluent用戶自定義標(biāo)量方程形式，采用Fluent自帶的求解器，實現(xiàn)水滴相x、y、z方向速度與水滴相體積分?jǐn)?shù)的迭代求解。

1.1.3 防冰熱分析計算

防冰熱分析計算采用改進的Messinger模型，考慮過冷水滴撞擊防冰腔外表面后未完全蒸發(fā)而形成的溢流流動影響，但暫不考慮防冰不足導(dǎo)致溢流冰形成的情況。此時防冰腔外表面的熱質(zhì)傳遞如圖1［32］所示。

對于防冰腔外表面任意控制體積，在穩(wěn)態(tài)情況下均滿足質(zhì)量守恒與能量守恒，即

圖1 考慮溢流水流動的防冰腔外表面熱質(zhì)傳遞［32］Fig.1 Heat and mass transfer on external surface of anti-icing chamber considering runback flow［32］

式中：為溢流流入質(zhì)量流量；為過冷水滴撞擊質(zhì)量流量；為溢流流出質(zhì)量流量；為蒸發(fā)質(zhì)量 流量；為防冰系 統(tǒng)外表 面加熱量；為水滴撞擊帶入的總焓值，包括水滴撞擊動能及水滴自身所攜帶的焓；為溢流流入帶入的滯止 焓；為表面 對 流 換 熱量；為 表 面蒸發(fā)帶走的總焓值；為溢流流出帶走的滯止焓。以上各項具體表達式詳見文獻［32］。

由于防冰表面溢流流動主要受氣動剪切力與表面壓力梯度驅(qū)動，在短艙進氣道周向非對稱流場的作用下，會呈現(xiàn)較強的三維特性。如圖2所示，圖中CVi為控制體積，Ai為面法向量，為流入當(dāng)前控制體積的質(zhì)量流量，為流出當(dāng)前控制體積的質(zhì)量流量，i=1,2,…,4，j=1,2，防冰表面任意控制體積與其相鄰的上下游控制體積之間可能存在多個溢流流入流出面，每個面上的流入流出質(zhì)量流量可根據(jù)質(zhì)量守恒并結(jié)合氣流流動特性計算得到。

圖2 三維防冰表面溢流分配示意圖Fig.2 Schematic of runback distribution on 3D antiicing surface

通過Fluent用戶自定義函數(shù)（User-Defined Functions ，UDFs），將水滴撞擊、溢流流動及表面蒸發(fā)等能量源項引入Fluent計算的能量方程中，根據(jù)能量守恒在給定的防冰腔外表面溫度邊界條件下求解表面加熱量Q?anti?icing，wet。

以外流環(huán)境溫度為參考溫度，根據(jù)防冰腔外表面加熱量換算得到濕空氣條件下防冰腔外表面等效換熱系數(shù)：

式中：Tw，wet為濕空氣條件下防冰腔外表面溫度；T∞為外流場環(huán)境溫度；A為換熱面積。

考慮固體域?qū)?，穩(wěn)態(tài)情況下防冰腔內(nèi)外表面達到熱平衡：

式中：為防冰腔內(nèi)表面對流換熱量；為固體域?qū)崃?。濕空氣條件下防冰腔內(nèi)表面對流換熱系數(shù)計算方法與干空氣條件下保持一致，僅防冰腔外表面熱載荷發(fā)生變化，由此可迭代得到濕空氣條件下防冰腔固體域溫度場。

1.2 內(nèi)外流松散耦合迭代策略

1.2.1 干空氣流場計算

圖3 干空氣流場計算的內(nèi)外流松散耦合迭代策略Fig.3 Loose-coupling iterative strategy of internal and external flow for flow calculation under dry air condition

干空氣條件下內(nèi)外流松散耦合迭代策略如圖3所示，每一迭代輪次中均包含外流計算、內(nèi)流與固體域計算。除第1輪次（Run 1）為初場計算外，其余輪次（Run 2~N）內(nèi)均先開展內(nèi)流與固體域計算后再進行外流計算。同一迭代輪次內(nèi)，內(nèi)外流耦合的數(shù)據(jù)交互為：內(nèi)流與固體域計算獲得的防冰腔外表面溫度分布賦值給外流計算作為第2類邊界條件。相鄰輪次耦合迭代計算的數(shù)據(jù)傳遞為：前一輪次（RunN?1）外流計算獲得的防冰腔外表面對流換熱系數(shù)分布作為后一輪次（RunN）內(nèi)流與固體域計算的第3類邊界條件。計算過程中，不斷更新防冰腔外表面對流換熱系數(shù)分布與溫度分布，當(dāng)相鄰迭代輪次的表面平均溫度偏差不超過±0.5 ℃時，可判定迭代收斂。

1.2.2 濕空氣防冰熱分析計算

圖4 濕空氣流場計算的內(nèi)外流松散耦合迭代策略Fig.4 Loose-coupling iterative strategy of internal and external flow for flow calculation under wet air condition

濕空氣條件下，內(nèi)外流松散耦合迭代策略如圖4所示。以干空氣條件下內(nèi)外流耦合計算的收斂結(jié)果作為初場，引入水滴撞擊對防冰表面換熱特性的影響，外流流場計算的同時還需開展防冰表面熱分析計算。同一迭代輪次中，先開展外流防冰熱分析計算，提取防冰腔外表面等效換熱系數(shù)作為第3類邊界條件賦值給內(nèi)流與固體域計算。相鄰迭代輪次中，傳遞的數(shù)據(jù)則為內(nèi)流與固體域計算獲得的防冰腔外表面溫度分布。在反復(fù)多次耦合迭代后，不斷更新防冰腔外表面參數(shù)，迭代收斂依據(jù)仍為相鄰迭代輪次防冰腔外表面平均溫度偏差不超過±0.5 ℃。

2 計算模型與網(wǎng)格劃分

2.1 計算模型與計算域

短艙笛形管熱氣防冰系統(tǒng)三維全尺寸模型如圖5所示。將短艙所處外流域獨立為一個計算域，防冰腔固體域和內(nèi)流域合并為另一個計算域。以防冰腔蒙皮外表面作為兩個計算域的耦合面，通過耦合面數(shù)據(jù)交互實現(xiàn)內(nèi)外流域的耦合迭代。根據(jù)計算域拆分方式，依次完成外流計算域和內(nèi)流計算域的網(wǎng)格劃分。

圖5 計算模型與內(nèi)外流計算域Fig.5 Computational model and computational do?mains of internal and external flow

2.2 外流計算域網(wǎng)格與邊界條件

采用ANSYS ICEM軟件對計算域進行網(wǎng)格劃分。外流計算域網(wǎng)格如圖6所示，總網(wǎng)格數(shù)約444萬。為了準(zhǔn)確模擬附面層內(nèi)的復(fù)雜流動特性，短艙及其延長段的近壁面網(wǎng)格均進行了局部加密處理，附面層第1層網(wǎng)格高度約為0.01 mm，保證壁面y+≈1，網(wǎng)格尺度與選用的剪切應(yīng)力輸運（Shear Stress Transfer， SST）κ-ω湍流模型相匹配；另外，在進氣道唇口附近等流場參數(shù)變化較為劇烈的區(qū)域也適當(dāng)進行了網(wǎng)格加密處理。

圖6 外流計算域網(wǎng)格Fig.6 Mesh for external computational domain

外流計算域邊界條件設(shè)置如下：外流邊界為壓力遠(yuǎn)場；風(fēng)扇進口為壓力出口；進氣道防冰腔外表面為內(nèi)外流松散耦合迭代時的數(shù)據(jù)交互面，盡管該表面邊界類型始終為壁面邊界條件，但賦值的溫度分布會隨迭代輪次不斷更新；除上述邊界以外，其余均為絕熱壁面。

外流計算選用密度基求解器，流域中空氣物性參數(shù)選用理想氣體模型進行計算。松散耦合迭代過程中每一輪次的外流計算均同時監(jiān)控流場殘差、防冰腔外表面平均換熱系數(shù)及風(fēng)扇進口流量，計算收斂后流場結(jié)果才可用于后續(xù)迭代。計算的收斂標(biāo)準(zhǔn)為各殘差指標(biāo)下降3個數(shù)量級或殘差曲線不再下降，并且外表面平均換熱系數(shù)和風(fēng)扇進口流量保持穩(wěn)定。

2.3 內(nèi)流計算域網(wǎng)格與邊界條件

內(nèi)流與固體域組成的計算網(wǎng)格如圖7所示。為精確模擬熱氣由射流孔噴射沖擊至防冰腔內(nèi)表面的過程，需對每個小孔的射流方向及其沿程周向進行網(wǎng)格局部加密。同時，由于防冰腔內(nèi)表面射流沖擊流動換熱特性較為復(fù)雜，為提高計算精度，需對內(nèi)表面近壁面網(wǎng)格也進行局部加密，第1層網(wǎng)格高度約為0.01 mm，確保壁面y+與SSTκ-ω湍流模型相匹配。三維全環(huán)模型的內(nèi)流與固體域計算網(wǎng)格總數(shù)約為1 822萬。

圖7 內(nèi)流計算域網(wǎng)格Fig.7 Mesh for internal computational domain

內(nèi)流與固體域計算邊界條件如下：防冰腔和前隔框的內(nèi)表面均為流固耦合面；前隔框外表面為絕熱壁面；笛形管管壁為定壁溫邊界條件；管壁上均布的射流孔為壓力入口邊界條件；防冰腔外表面仍為內(nèi)外流松散耦合計算的數(shù)據(jù)交互面，壁面邊界條件賦值的換熱系數(shù)也會隨耦合迭代過程不斷更新。

內(nèi)流與固體域計算時，同樣選用密度基求解器，流域中熱空氣物性參數(shù)選用理想氣體模型進行計算。松散耦合迭代過程中每一輪次內(nèi)流與固體域計算均同時監(jiān)控流場殘差和防冰腔外表面平均溫度。計算的收斂標(biāo)準(zhǔn)同樣要求流場殘差與監(jiān)控的表面參數(shù)收斂曲線保持不變。

2.4 計算工況

計算工況如表1所示，包括飛行狀態(tài)參數(shù)、結(jié)冰氣象條件以及射流孔出流參數(shù)。

表1 計算工況Table 1 Computational conditions

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

3.1 算例驗證

截取短艙進氣道12點鐘方向局部模型，依次開展強固耦合計算與松散耦合計算，通過對比分析以驗證本文松散耦合計算方法及耦合迭代策略的有效性。圖8為驗證模型的幾何及計算網(wǎng)格。

圖8 計算方法驗證模型及網(wǎng)格Fig.8 Model and mesh for computational method validation

圖9 干空氣條件下強固耦合與松散耦合計算結(jié)果對比 （z=0）Fig.9 Comparison between tight-coupling and loosecoupling methods under dry air condition （z=0）

圖10 濕空氣條件下強固耦合與松散耦合計算結(jié)果對比（z=0）Fig.10 Comparison between tight-coupling and loosecoupling methods under wet air condition （z=0）

圖9和圖10分別為干、濕空氣條件下強固耦合計算方法與松散耦合計算方法的對比結(jié)果，其中h0為對流換熱系數(shù)基準(zhǔn)值，T0為溫度基準(zhǔn)值，圖9所示干空氣條件下松散耦合計算結(jié)果為迭代3個輪次后的收斂結(jié)果，圖10所示濕空氣條件下計算結(jié)果為迭代4個輪次后的收斂結(jié)果?？梢钥闯鰞煞N計算方法所得到的防冰腔內(nèi)外表面各參數(shù)基本保持一致。由此可知，所采用的松散耦合計算方法及迭代策略可以達到與強固耦合計算方法基本一致的仿真精度。

3.2 干空氣流場

圖11 干空氣條件下短艙進氣道防冰腔外表面對流換熱系數(shù)分布Fig.11 Convective heat transfer coefficient distributions on external surface of nacelle inlet anti-icing chamber under dry air condition

圖12 干空氣條件下短艙進氣道防冰腔外表面溫度分布Fig.12 Temperature distributions on external surface of nacelle inlet anti-icing chamber under dry air condition

圖13 干空氣條件下短艙進氣道防冰腔內(nèi)表面對流換熱系數(shù)分布Fig.13 Convective heat transfer coefficient distributions on internal surface of nacelle inlet anti-icing chamber under dry air condition

圖11~圖13為干空氣條件下短艙進氣道笛形管防冰腔外表面對流換熱系數(shù)、溫度以及內(nèi)表面對流換熱系數(shù)在內(nèi)外流松散耦合不同迭代輪次下的計算結(jié)果。在經(jīng)過3個輪次的內(nèi)外流耦合迭代計算后，防冰腔內(nèi)外表面各參數(shù)分布逐漸趨于收斂。其中，由于干空氣條件下第1輪次為初場計算，防冰腔外表面溫度設(shè)置為定壁溫，因此圖12中并未給出該輪次下的溫度分布，但增加了第4輪次內(nèi)流與固體域計算結(jié)果。對比耦合迭代的第3輪次和第4輪次防冰腔外表面溫度分布可以看出，兩者基本達到一致。

圖14 干空氣條件下短艙進氣道防冰腔內(nèi)外表面參數(shù)分布曲線（12點鐘方向）Fig.14 Parameter distribution curves on internal and ex?ternal surfaces of nacelle inlet anti-icing chamber under dry air condition （12 o’clock direction）

如圖14所示，取短艙進氣道防冰腔在12點鐘方向上的參數(shù)分布曲線，可以更加準(zhǔn)確地分析干空氣條件下內(nèi)外流松散耦合迭代計算的收斂過程。可以看出，第1輪次防冰腔內(nèi)外表面對流換熱系數(shù)分布均與第2輪次計算結(jié)果有較為明顯的差距，但第2輪次和第3輪次各參數(shù)分布曲線基本重合；第3、4輪次的表面溫度曲線吻合良好，可達到迭代收斂要求；防冰腔內(nèi)表面對流換熱系數(shù)分布曲線在迭代過程中變化不大，僅射流駐點區(qū)域隨迭代輪次略有下降。綜上，干空氣條件下防冰腔內(nèi)、外表面平均對流換熱系數(shù)在耦合迭代過程中均呈現(xiàn)先減后增再逼近收斂結(jié)果的變化規(guī)律，而防冰腔外表面平均溫度則先增后減至迭代收斂。

3.3 水滴撞擊特性

圖15 短艙進氣道防冰腔外表面局部水收集系數(shù)分布Fig.15 Local collection efficiency distributions on ex?ternal surface of nacelle inlet anti-icing chamber

圖15為短艙進氣道防冰腔外表面局部水收集系數(shù)分布，圖中β0為局部水收集系數(shù)基準(zhǔn)值?？梢钥闯?，水滴撞擊區(qū)域主要集中在唇口前緣駐點位置附近，且唇口內(nèi)表面撞擊極限要明顯大于唇口外表面，沿唇口周向各截面上局部水收集系數(shù)分布并不均勻，12點鐘方向具有最大的局部水收集系數(shù)峰值和撞擊范圍，而6點鐘方向撞擊范圍明顯減小，3點鐘和9點鐘方向左右對稱，具有相同的局部水收集系數(shù)分布。由于水滴撞擊特性計算結(jié)果是后續(xù)防冰熱分析計算的輸入，唇口表面局部水收集系數(shù)的周向分布不均勻性很大程度上會影響濕空氣條件下防冰腔表面的換熱特性及溫度分布。

3.4 濕空氣防冰熱分析

圖16~圖18為濕空氣條件下短艙進氣道笛形管防冰腔外表面等效換熱系數(shù)、溫度以及內(nèi)表面對流換熱系數(shù)在內(nèi)外流松散耦合不同迭代輪次下的計算結(jié)果。相比于干空氣條件，濕空氣條件下耦合迭代輪次明顯增加，這是由于該工況條件下表面撞擊水量不能完全蒸發(fā)并在防冰腔表面形成周向不均勻分布的溢流水流動，如圖19所示，其中m0為溢流水質(zhì)量流量基準(zhǔn)值。溢流水量分布與表面換熱特性相互影響，導(dǎo)致表面溫度分布收斂速度減緩。防冰腔外表面由于過冷水滴撞擊引入的熱質(zhì)傳遞源項，濕空氣條件下外表面等效換熱系數(shù)要明顯大于干空氣條件；而對于防冰腔內(nèi)表面對流換熱系數(shù)，干、濕空氣條件下數(shù)值變化并不大。此外，相較于干空氣條件，濕空氣條件下進氣道防冰腔外表面溫度明顯下降，且表面溫度均勻性有所提升。

圖17 濕空氣條件下短艙進氣道防冰腔外表面溫度分布Fig.17 Temperature distributions on external surface of nacelle inlet anti-icing chamber under wet air condition

圖18 濕空氣條件下短艙進氣道防冰腔內(nèi)表面對流換熱系數(shù)分布Fig.18 Convective heat transfer coefficient distributions on internal surface of nacelle inlet anti-icing chamber under wet air condition

圖19 濕空氣條件下短艙進氣道防冰腔外表面溢流水質(zhì)量流量分布（第8輪次迭代）Fig.19 Runback mass flux on external surface of na?celle inlet anti-icing chamber under wet air con?dition （Run 8）

圖20 濕空氣條件下短艙進氣道防冰腔內(nèi)外表面參數(shù)分布曲線（12點鐘方向）Fig.20 Parameter distributions curves on internal and external surfaces of nacelle inlet anti-icing chamber under wet air condition （12 o’clock direction）

結(jié)合圖20所示的短艙進氣道防冰腔內(nèi)外表面各參數(shù)在12點鐘方向上的分布曲線可以看出，在撞擊區(qū)域內(nèi)松散耦合迭代的第1輪次到第4輪次防冰腔外表面溫度下降明顯，且外表面等效換熱系數(shù)分布也發(fā)生較大的變化，但從第5輪次開始到第8輪次計算結(jié)果都較為接近，尤其在唇口內(nèi)表面各參數(shù)曲線吻合良好，但在唇口外表面撞擊極限之后的蒙皮外表面等效換熱系數(shù)及溫度分布還略有差異?？偟膩碚f，濕空氣條件下防冰腔內(nèi)外表面參數(shù)在耦合迭代過程中所呈現(xiàn)的變化規(guī)律與干空氣條件恰好相反。防冰腔外表面平均等效換熱系數(shù)與內(nèi)表面平均對流換熱系數(shù)在迭代過程中均先增加后減小再逼近收斂，而防冰腔外表面平均溫度則先降低后升高直至迭代收斂。

4 結(jié) 論

針對三維全尺寸短艙進氣道笛形管熱氣防冰系統(tǒng)，開展了內(nèi)外流耦合計算方法研究，結(jié)合進氣道防冰計算流程，確定了干、濕空氣條件下短艙防冰系統(tǒng)內(nèi)外流松散耦合迭代策略及計算收斂所需的迭代輪次，得到的主要結(jié)論如下：

1）干、濕空氣條件下內(nèi)外流松散耦合迭代均以防冰腔外表面作為內(nèi)、外計算域的耦合面，并通過耦合面溫度和換熱系數(shù)相互賦值為邊界條件來實現(xiàn)內(nèi)外計算域的數(shù)據(jù)交互，其中濕空氣條件下考慮過冷水滴撞擊引入的熱質(zhì)傳遞影響及三維溢流效應(yīng)，采用防冰腔外表面等效換熱系數(shù)作為耦合面邊界條件。

2）干空氣條件下短艙進氣道防冰腔內(nèi)外流松散耦合僅需3個輪次的迭代計算即可收斂。

3）濕空氣條件下進氣道防冰腔內(nèi)表面在松散耦合迭代4個輪次后趨于收斂，但受表面溢流流動換熱影響，防冰腔外表面溫度分布在耦合迭代過程中仍有一定波動。后續(xù)將針對該問題，對內(nèi)外流耦合迭代策略作進一步優(yōu)化。