飛行時，大氣中過冷水會於飛機翼表面及機翼前緣結晶形成積冰現象，進而改變飛機翼表面粗糙度及流場氣動力特性，造成飛機翼升阻比值驟降，導致飛機失速的情況發生。
為克服積冰問題，現今的飛機大都附屬配載除冰系統，但無人飛行載具(UAV)因酬載的限制，無法再配附除冰系統，因此探討積冰現象對UAV的影響格外重要。本研究以國外設計UAV常使用的GLC305及NLF0414 翼型為實驗機翼原型，探討其產生犄角狀透明積冰(Glaze Ice)現象後對氣動力的影響。本研究以風洞實驗及三維CFD計算求解，分析積冰生成前後氣動力特性之改變及流場變化等現象，並將所得之結果相互驗證比對。
研究結果顯示，積冰的現象造成機翼表面升力減少約46％、阻力增加達25％，失速攻角則由16度提前至8度。

During flying through the atmosphere with the existence of supercooled water, there will be ice accretion on the leading edge and surface of the aircraft airfoil. This iced airfoil will cause the change of surface roughness, lift/drag ratio, and aerodynamic characteristics of the airfoil, and finally results in aircraft stall.
Actually, almost all the aircrafts have the on-board ice protection system in order to overcome the problem of ice accretion. However due to the limitation of payload, UAV does not have a complete set of ice protection system in general. Therefore, the study of the effects of ice accretion on UAV is particularly important. The airfoils of UAV selected for the test model in this study were clean and iced (Glaze Ice) GLC305 and NLF0414. Both wind-tunnel experiments and 3-D CFD numerical approaches were used to investigate the aerodynamic and flowfield characteristics before and after ice accretion.
The results of this study indicate that ice accretion on airfoil will cause lift decreasing about 46% and drag increasing about 25%, and that the stalling angle will be changed from 16 degree to 8 degree.

誌謝 ii
摘要 iii
ABSTRACT iv
目錄 v
表目錄 viii
圖目錄 x
符號說明 xix
1. 前言 1
1.1 研究背景 1
1.2 文獻回顧 3
1.3 研究動機 8
2. 積冰介紹 10
2.1 積冰生成之條件 10
2.2 積冰種類 11
2.2.1 透明冰 12
2.2.2 霜狀冰 13
2.2.3 混合冰 14
2.2.4 積霜 15
2.3 預防積冰的產生 15
3. 研究方法 18
3.1 數值模擬 20
3.1.1 幾何構形 20
3.1.2 計算網格系統製作 24
3.1.3 邊界條件設定 27
3.1.4 統御方程式 29
3.1.5 紊流模式驗證 31
3.1.6 收斂標準與研究矩陣 32
3.1.7 程式驗證 33
3.2 風洞實驗實驗設備架設與流程 39
3.2.1 開放式低速風洞 39
3.2.2 六分量平衡儀 41
3.2.3 攻角調整機構 50
3.2.4 風洞測試區流場均勻度量測 54
3.2.5 實驗流程規劃 57
4. 結果與討論 61
4.1 風洞量測實驗 61
4.1.1 積冰對升力係數之影響 61
4.1.2 積冰對阻力係數之影響 66
4.2 數值模擬 70
4.2.1 GLC305翼型 70
4.2.2 NLF0414翼型 81
4.3 風洞實驗與數值模擬相互驗證 92
4.4 雷諾數2×〖10〗^6之數值模擬 98
4.4.1 GLC305翼型 99
4.4.2 NLF0414翼型 107
5. 結論與未來研究方向 116
5.1 結論 116
5.2 未來研究方向 117
參考文獻 118
附錄A 六分量平衡儀各分量校正電壓值 120
附錄B 六分量平衡儀校正矩陣 123
自傳 124

[1]http://www.aviationweather.ws/page093-1.jpg
[2]Berkowitz, B., “Air Florida,Inc.,Boeing 737-222,N62af,” NTSB Report, Washington, U.S.A , NTSB- AAR-82-08, pp.1-136, 1982.
[3]Mingione, G., “In-flight Icing Encounter and Loss of Control Simmons Airlines, d.b.a. American Eagle Flight 4184 Avions de Transport Regional (ATR) Model 72-212, N401AM, Roselawn, Indiana,” NTSB Report, Washington, U.S.A, PB96-91040I, pp.1-322, 1996
[4]Harold, E. A., Broeren, P. A., Zoeckler, J. G. and Lee, S., “A Wind Tunnel Study of Icing Effects on a Business Jet Airfoil,” NASA/TM, U.S.A, 2003-212124, pp1-17, 2003.
[5]Bragg, M. B., Broeren, A. P. and Blumenthal, L. A., “Iced-airfoil aerodynamics,” Progress in Aerospace Sciences, Vol. 41, pp. 323-362, 2005.
[6]Bragg, M. B. and Broeren, A. P., “Airfoil Ice-Accretion Aerodynamics Simulation,” NASA/TM, U.S.A, 2008-214830, pp. 1-21, 2007.
[7]Mirzaei, M., Ardekani, M. A. and Doosttalab, M., “Numerical and experimental study of flow field characteristics of an iced airfoil,” Aerospace Science and Technology, Vol. 13, pp. 267-276, 2009.
[8]鄭仁杰，林鴻佳，苗君易， “無人飛機長滯空高升力翼型設計及積冰對於翼型氣動力特性影響” ，航空、太空及民航學刊，第42卷，第5期，第1-10 頁，2010。
[9]朱家輝， “積冰現象之空氣動力特性分析-以ATR-72為例” ，碩士論文，淡江大學，台北，2008。
[10]張國治， “赴美國聖地牙哥參加飛航危險天氣與飛機積冰訓練課程” ， 行政院飛航安全委員會，台北，C09305531，pp1-22，2004。
[11]http://extice.cira.it/project/proposal-summary/icing-aerodynamics/figure/aero aerofig2.png
[12]Vargas, M. and Papadakis, M., “Ice Accretions on a Swept GLC-305 Airfoil,” NASA/TM, U.S.A, SAE 02GAA–43, pp3-14, 2002.
[13]http:// www.virtualskies.arc.nasa.gov/weather/4.html
[14]易賢， “飛機積冰的數值計算與積冰試驗相似準則研究” ，博士論文，中國空氣動力研究與發展中心，北京，2007。
[15]http://www.flightops.sita.aero/images/De-icing.jpg
[16]http://www.flyingway.com/airlecture/ice_rain_systems.html
[17]蘇晉興， “微飛行器之氣動力構型設計與分析” ，碩士論文，國防大學中正理工學院，桃園，2006。
[18]韓點忠，FLUENT 流體工程仿真計算案例與應用，北京理工大學出版社,第210-228頁，2004。
[19]徐渝昌， "微飛行器用低雷諾數機翼氣動力性能之風洞試驗研究" ，國防大學中正理工學院，桃園，2000。