專利名稱：基于液態(tài)水分布的結(jié)冰表面粗糙度衡量方法
技術(shù)領(lǐng)域：
本發(fā)明涉及一種飛機(jī)結(jié)冰數(shù)值計(jì)算方法，尤其涉及一種基于液態(tài)水分布的結(jié)冰表面粗糙度衡量方法，屬于冰層表面粗糙度評(píng)估領(lǐng)域。
背景技術(shù)：
在飛機(jī)穿越含有過(guò)冷水滴的云層時(shí)，機(jī)翼前緣通常會(huì)發(fā)生結(jié)冰，造成飛機(jī)的飛行性能下降。為了保證飛機(jī)的飛行安全，并減少試驗(yàn)資源消耗，提高研究效率，必須大力發(fā)展飛機(jī)結(jié)冰數(shù)值模擬技術(shù)。而在結(jié)冰數(shù)值計(jì)算中，需要模擬冰層隨時(shí)間的增長(zhǎng)過(guò)程，因此就涉及到每個(gè)時(shí)間步長(zhǎng)上冰層表面的粗糙程度衡量問(wèn)題，即每個(gè)時(shí)間步長(zhǎng)后都必須根據(jù)當(dāng)前時(shí)刻表面積冰情 況更新粗糙度，合理衡量積冰對(duì)氣流壓力系數(shù)和水滴收集系數(shù)的影響，以及對(duì)冰層、過(guò)冷水及空氣三相的傳熱傳質(zhì)產(chǎn)生的影響。目前常規(guī)的結(jié)冰表面粗糙度衡量方法是等效沙粒粗糙高度模型，即把冰層表面的粗糙程度等效為等直徑排列的沙粒，計(jì)算得出常數(shù)粗糙度值，在計(jì)算過(guò)程中表面冰層外推后，不再進(jìn)行粗糙度更新，整個(gè)結(jié)冰計(jì)算中采用該常數(shù)作為冰層的粗糙度值。這種方法沒(méi)有把結(jié)冰粗糙度表達(dá)為空間和時(shí)間的函數(shù)，無(wú)法表達(dá)出冰層不斷形成及模型不同位置上粗糙度的差異，因此結(jié)冰計(jì)算的誤差較大。

發(fā)明內(nèi)容
本發(fā)明針對(duì)現(xiàn)有技術(shù)的缺陷，而提出一種基于液態(tài)水分布的結(jié)冰表面粗糙度衡量方法，以實(shí)現(xiàn)結(jié)冰過(guò)程中冰層表面粗糙度的實(shí)時(shí)更新。該方法的內(nèi)容如下根據(jù)控制體的質(zhì)量守恒原理獲得結(jié)冰表面的水膜高度hw，根據(jù)受力平衡模型，確定水膜發(fā)生破斷的臨界水膜高度hfi
「 I , r6·/^2 ·σ· (I-CosS)1^hf—r = [―^——--]3當(dāng)hw > hf_r時(shí),水膜連續(xù),粗糙度ks等于水膜表面波高,即:ks =0.75· —· I—；
Mw \ g當(dāng)hw < hf_r時(shí)，水膜破斷形成溪流，粗糙度ks等于當(dāng)?shù)叵鞲叨?，即ks =R(l_cos θ )；以上公式中yw為水動(dòng)力粘度，O為水膜表面張力，Θ為水在冰面上的接觸角，P w為水密度，τ a為水膜受到的氣流剪切應(yīng)力，g為重力加速度，R為溪流半徑。技術(shù)效果本方法可在每個(gè)時(shí)間步長(zhǎng)上對(duì)冰層表面的粗糙度進(jìn)行更新，很好地解決了結(jié)冰計(jì)算中冰層表面粗糙程度的衡量問(wèn)題，使得冰層表面粗糙度評(píng)估更接近實(shí)際結(jié)冰過(guò)程，為提高后續(xù)結(jié)冰傳熱傳質(zhì)和冰形的模擬精度奠定了基礎(chǔ)。


圖I為結(jié)冰表面水膜流動(dòng)示意圖，圖中字母含義g表示重力加速度，gs表示重力加速度沿表面切線方向的分量，gy表示重力加速度沿表面法線方向的分量，Pa表示氣流作用于水膜表面的壓強(qiáng)，Ta(S)表示氣流作用于水膜表面的剪切應(yīng)力，u(s，y)表示水膜內(nèi)速度分布，hj表冰層厚度，hw表水膜厚度。圖2為水膜破斷過(guò)程示意圖，圖中dz表示展向微元長(zhǎng)度。圖3為水膜破斷位置受力分析示意圖，圖中字母含義Θ為水在冰面上的接觸角，F(xiàn)p為上游動(dòng)能轉(zhuǎn)化成的水膜破斷位置的壓力，F(xiàn)"、F02為水膜破斷位置弧形輪廓線上所受的表面張力。圖4為溪流模型示意圖，圖中字母含義Θ為水在冰面上的接觸角，R為溪流半徑，X為溪流間距。圖5為結(jié)冰模擬計(jì)算結(jié)果對(duì)比圖。
具體實(shí)施例方式下面對(duì)本發(fā)明作進(jìn)一步說(shuō)明。結(jié)冰表面粗糙度是由表面的結(jié)冰過(guò)程所誘導(dǎo)產(chǎn)生的，根據(jù)試驗(yàn)現(xiàn)象，在結(jié)冰氣象條件下機(jī)翼結(jié)冰表面存在一層薄水膜，如圖I所示。由于水膜在向下游流動(dòng)的過(guò)程中不斷與氣流及冰層之間發(fā)生熱質(zhì)交換，其高度會(huì)不斷減小，當(dāng)其高度達(dá)到臨界高度時(shí)，水膜破斷形成溪流，如圖2所示。因此，本方法針對(duì)結(jié)冰表面不同的液態(tài)水分布狀態(tài)，建立分區(qū)的粗糙度模型，具體內(nèi)容如下根據(jù)控制體的質(zhì)量守恒原理獲得機(jī)翼結(jié)冰表面任意處的水膜高度hw，根據(jù)受力平
衡模型，確定水膜發(fā)生破斷的臨界水膜高度hf_p水膜破斷位置的受力分析如圖3所示。
權(quán)利要求
1.ー種基于液態(tài)水分布的結(jié)冰表面粗糙度衡量方法，其特征在于 該方法的內(nèi)容如下 根據(jù)控制體的質(zhì)量守恒原理獲得結(jié)冰表面的水膜高度hw，根據(jù)受カ平衡模型，確定水膜發(fā)生破斷的臨界水膜高度hfi
全文摘要
本發(fā)明公開(kāi)了一種基于液態(tài)水分布的結(jié)冰表面粗糙度衡量方法，屬于冰層表面粗糙度評(píng)估領(lǐng)域。該方法是針對(duì)結(jié)冰表面不同的液態(tài)水分布狀態(tài)，建立分區(qū)的粗糙度模型，即把結(jié)冰表面的水膜高度與水膜發(fā)生破斷的臨界水膜高度進(jìn)行比較，當(dāng)前者較大時(shí)，粗糙度等于水膜表面波高，當(dāng)后者較大時(shí)，粗糙度等于水膜破斷形成的溪流高度。本方法在每個(gè)時(shí)間步長(zhǎng)上都根據(jù)實(shí)際冰層表面液態(tài)水特征更新粗糙度值，使粗糙度評(píng)估更接近實(shí)際的結(jié)冰過(guò)程。
文檔編號(hào)G06F19/00GK102663238SQ20121007588
公開(kāi)日2012年9月12日 申請(qǐng)日期2012年3月21日 優(yōu)先權(quán)日2012年3月21日
發(fā)明者孫志國(guó), 朱春玲, 檀妹靜 申請(qǐng)人:南京航空航天大學(xué)