飛行防冰與結(jié)冰中的流動(dòng)傳熱計(jì)算方法研究.pdf

1、結(jié)冰現(xiàn)象對(duì)于飛行安全來(lái)說(shuō)是一個(gè)重大的威脅。當(dāng)飛機(jī)在結(jié)冰氣象條件下飛行時(shí)，如不采取適當(dāng)?shù)姆辣胧?，飛機(jī)迎風(fēng)部件表面就會(huì)發(fā)生結(jié)冰。近年來(lái)，由于結(jié)冰引發(fā)的飛行安全事故仍在國(guó)內(nèi)外多次發(fā)生，因此飛行結(jié)冰和防冰依然是國(guó)內(nèi)外學(xué)者廣泛關(guān)注的一個(gè)熱點(diǎn)問(wèn)題。
　　作為飛行結(jié)冰和防冰問(wèn)題的重要研究手段，數(shù)值模擬技術(shù)在近年來(lái)的結(jié)冰和防冰研究中得到了十分廣泛的應(yīng)用和發(fā)展。本文針對(duì)目前國(guó)內(nèi)外在防冰計(jì)算領(lǐng)域的研究現(xiàn)狀及其發(fā)展特點(diǎn)和趨勢(shì)，以商業(yè)CFD軟件ANSY

2、S Fluent作為開(kāi)發(fā)平臺(tái)，圍繞防冰和結(jié)冰模擬的通用計(jì)算方法開(kāi)展了如下研究工作：
　　首先，本文基于ANSYS Fluent軟件的二次開(kāi)發(fā)功能，對(duì)水滴撞擊特性的計(jì)算方法進(jìn)行了詳細(xì)研究，成功地在ANSYS Fluent軟件中實(shí)現(xiàn)了水滴撞擊特性的拉格朗日法和歐拉法計(jì)算。對(duì)于前者，本文基于Fluent軟件的離散相模型獲得了水滴的運(yùn)動(dòng)軌跡，并采用自己發(fā)展的粒子統(tǒng)計(jì)方法實(shí)現(xiàn)了復(fù)雜三維表面局部水收集系數(shù)的拉格朗日法計(jì)算。對(duì)于后者，本文采用F

3、luent軟件的用戶(hù)自定義標(biāo)量方程對(duì)水滴相運(yùn)動(dòng)的控制方程進(jìn)行了求解，從而可以在Fluent軟件支持的任何坐標(biāo)系及邊界條件下獲得部件表面的局部水收集系數(shù)?；谏鲜龇椒?，本文對(duì)二維和三維翼型的局部水收集系數(shù)進(jìn)行了算例驗(yàn)證計(jì)算，得到了與文獻(xiàn)數(shù)據(jù)相符的計(jì)算結(jié)果。同時(shí)，本文對(duì)水滴撞擊特性的計(jì)算方法進(jìn)行了拓展研究，在Fluent軟件中實(shí)現(xiàn)了過(guò)冷大水滴條件和二維軸對(duì)稱(chēng)坐標(biāo)系下的局部水收集系數(shù)計(jì)算。此外，本文發(fā)展了求解水滴撞擊特性的拉格朗日-歐拉復(fù)合解

4、法，通過(guò)對(duì)算例計(jì)算的結(jié)果表明，該方法能夠有效提高歐拉法在求解水滴撞擊特性時(shí)的穩(wěn)定性和效率。
　　其次，本文開(kāi)展了防冰條件下部件表面流動(dòng)換熱的計(jì)算方法研究。本文從防冰計(jì)算中經(jīng)典的Messinger模型入手，對(duì)部件表面水膜流動(dòng)的質(zhì)量守恒和能量守恒關(guān)系進(jìn)行了分析，并發(fā)展了基于非結(jié)構(gòu)化網(wǎng)格體系的表面水膜質(zhì)量平衡計(jì)算方法，以及表面溢流水膜質(zhì)量的分配方法。在防冰部件表面溫度的計(jì)算上，本文通過(guò) Fluent用戶(hù)自定義函數(shù)將考慮溢流水膜的防冰部件

5、表面各項(xiàng)熱流添加到了 Fluent的能量方程中，并基于Fluent的耦合換熱壁面邊界實(shí)現(xiàn)了防冰熱氣、固壁、以及外流場(chǎng)在防冰條件下的耦合換熱計(jì)算?；谏鲜龇椒?，本文針對(duì)航空發(fā)動(dòng)機(jī)靜止部件在防冰條件下的表面溫度進(jìn)行了計(jì)算，并將得到的結(jié)果與本文開(kāi)展的防冰試驗(yàn)結(jié)果進(jìn)行了對(duì)比，結(jié)果令人滿(mǎn)意。此外，本文還對(duì)不同算例計(jì)算中得到的表面水膜流動(dòng)質(zhì)量和表面熱流進(jìn)行了分析。
　　再次，本文開(kāi)展了水膜破裂后表面溪流形態(tài)的計(jì)算研究。為了獲得破裂處水膜的厚度

6、，本文根據(jù)水膜很薄的客觀事實(shí)對(duì)水膜采用了速度線(xiàn)性分布假設(shè)，并對(duì)水膜流動(dòng)中的質(zhì)量守恒和動(dòng)量守恒關(guān)系進(jìn)行了建模計(jì)算。為了得到水膜破裂后的溪流厚度和寬度，本文假設(shè)溪流截面為圓的一部分，根據(jù)破裂前后水膜和溪流的質(zhì)量守恒和能量守恒關(guān)系，發(fā)展了破裂后溪流幾何參數(shù)的求解方法。通過(guò)與文獻(xiàn)實(shí)驗(yàn)結(jié)果的對(duì)比，本文方法能夠在一定誤差范圍內(nèi)對(duì)破裂后的溪流形狀進(jìn)行預(yù)測(cè)?；谏鲜龇椒ǎ疚膶?duì)不同壁面溫度和來(lái)流水滴直徑下的航空發(fā)動(dòng)機(jī)支板表面溪流寬度和厚度進(jìn)行了計(jì)算，研

7、究了這些因素對(duì)破裂后溪流幾何參數(shù)的影響。
　　最后，本文開(kāi)展了結(jié)冰冰形的數(shù)值預(yù)測(cè)研究。本文根據(jù)結(jié)冰表面的質(zhì)量守恒和能量守恒關(guān)系，對(duì)結(jié)冰表面的熱力學(xué)模型進(jìn)行了研究，并發(fā)展了基于非結(jié)構(gòu)化網(wǎng)格的表面冰增長(zhǎng)迭代計(jì)算方法。為了獲得結(jié)冰冰形，本文發(fā)展了結(jié)冰后壁面節(jié)點(diǎn)坐標(biāo)的計(jì)算方法，并采用ANSYS Fluent動(dòng)網(wǎng)格技術(shù)對(duì)結(jié)冰后流體的計(jì)算域進(jìn)行了自動(dòng)更新，從而實(shí)現(xiàn)了結(jié)冰過(guò)程的多步預(yù)測(cè)。通過(guò)對(duì)經(jīng)典算例的驗(yàn)算，本文的冰形預(yù)測(cè)方法得到了驗(yàn)證?；谏?/p>