王凱朋,孟祥堯,滕月慧,魏照宇
(1. 上海交通大學 船舶海洋與建筑工程學院,上海 200240;2. 上海交通大學 海洋學院,上海 200240;3. 海軍研究院,北京 102442;4. 水聲對抗技術重點實驗室 上海 201108)
??湛缬蚝叫衅魇强茖W家與工程師綜合考慮水下航行器和空中飛行器的優(yōu)點之后,提出的一種可以實現(xiàn)水空兩域作業(yè)并且可以多次自由穿越水氣界面的新型??諆蓷w行器[1]。
??湛缬蚝叫衅鬟@一概念最早來源于1934年蘇聯(lián)的鮑里斯·烏沙可夫提出的“會飛行的潛艇”項目[2],但受限于當時的技術水平,該項目并未實現(xiàn)。直到21世紀之后,隨著科學技術的迅速發(fā)展,??湛缬蚝叫衅鞑叛杆侔l(fā)展起來?,F(xiàn)有的海空跨域航行器主要可以分為兩大類:固定翼形式和多旋翼形式。如2014年美國海軍研究實驗室提出的Flimmer[3]、2016年倫敦帝國理工學院制作的AquaMAV[4]、2017年加拿大舍布魯克大學研制的SUWAVE[5]、2017年美國北卡羅萊納州立大學研制的EagleRay[6]等都屬于固定翼形式的??湛缬蚝叫衅鳎欢?015年美國羅格斯大學研制的Naviator[7]、2016年美國奧克蘭大學設計的Loon Copter[8]、2016年空軍工程大學提出的HUAUV[9]等都屬于采用多旋翼形式的??湛缬蚝叫衅鳌?/p>
與固定翼和多旋翼飛行器類似,固定翼和多旋翼的??湛缬蚝叫衅饕哺饔衅鋬?yōu)缺點。多旋翼的??湛缬蚝叫衅麟m然飛行速度較慢、飛行距離相對較短,但體積較小、重量輕、隱蔽性好的優(yōu)勢使其具有更大的應用場合。此外結構較簡單、易于操縱的特性使其成本也更低廉。固定翼的??湛缬蚝叫衅麟m然成本較多旋翼偏高,且運動過程中不能懸停在某處,但優(yōu)點也比較明顯,其具有更快的飛行速度、更高的飛行高度以及更長的飛行距離。
為了具有更快的飛行速度以及更長的飛行距離,本文提出一款采用固定翼形式的??湛缬蚝叫衅?。當通過試驗方法來獲取航行器流體動力特性時,需要較長的試驗準備期,同時也需要耗費較大的人力物力財力,因此在設計初期主要利用數(shù)值模擬的方法進來研究。廖保全等[10-11]提出了一種可變形的跨介質航行器,并通過Fluent軟件對該航行器的氣動/水動特性進行了模擬,驗證了其可以同時滿足水下航行和空中飛行的要求。馮歡等[12]基于LBM-LES方法對兩棲水翼航行器的水動力特性進行分析,驗證了可以通過改變航行器外形使得該航行器滿足水下航行和空中飛行的要求。
本文提出一種基于升力原理的固定翼形式??湛缬蚝叫衅鳎瑧肧TAR-CCM+軟件,基于重疊網(wǎng)格方法對該航行器的氣動及水動力特性進行模擬。通過對該航行器數(shù)值模擬結果的分析,驗證其可以滿足在空中飛行和水下航行的要求,從而驗證該設計的合理性,為??湛缬蚝叫衅鞯脑O計提供參考。
空氣和水這2種流體的物理性質差別巨大。海水密度約為空氣的800多倍,動力粘性系數(shù)約為空氣的60倍。這就導致空中飛行時需要保持較高的航速,以產生足夠的升力來克服航行器的重力;而在水下航行時由于航行器受到的阻力較大,一般只能保持較低的航速。
海空跨域航行器的外形設計應從水下航行器和常規(guī)飛行器的外形著手。水下航行器和常規(guī)飛行器的外形有著明顯區(qū)別。水下航行器為了在水下具有足夠的耐壓能力,通常采用圓柱、橢圓等回轉體外形,而常規(guī)飛行器一般采用流線型設計以減小空氣阻力。目前,海空跨域航行器外形主要參考常規(guī)飛行器的氣動布局設計[13]。
本文提出如圖1所示的純三角翼??湛缬蚝叫衅?。該航行器與現(xiàn)有的固定翼飛行器外形相似,翼展為1 200 mm,機身長為975 mm,機翼的翼型為NASA SC(2)0402翼型,機翼與機身間有2°的安裝攻角。此外,為了提高航行器的橫向穩(wěn)定性,還設計了垂直尾翼,其采用了NACA0010翼型。經過計算發(fā)現(xiàn),加上動力、控制、推進器等裝置后該航行器質量為3.2 kg,模型具有較大的凈浮力。
圖1 ??湛缬蚝叫衅鲙缀瓮庑蜦ig. 1 Geometry of the aerial underwater vehicle
處理湍流數(shù)值計算問題的方法主要有3種:直接數(shù)值模擬方法(DNS方法)、尺度解析模擬方法、雷諾時均模擬方法(RANS方法)。受限于計算方法和計算機內存及性能,實際主要應用RANS方法來求解。
RANS方法的基本原理為:在時間域上對流場物理量進行雷諾平均化處理,然后求解所得到的時均化控制方程。
對于流體的控制方程主要有連續(xù)性方程、動量方程、能量守恒方程等。其通用形式可以表示為:
式中:φ為通用變量;Γ為廣義擴散系數(shù);S為廣義源項。
當φ=1,Γ=0,S=0時,方程為連續(xù)性方程;當φ=ui,Γ=μ,S=??p/?xi+Si時,方程為動量方程[14]。
RANS方法在時間域上對流場進行平均化的過程中會失去很多流場的細節(jié)信息,為了找回這些信息并使湍流基本方程封閉,就必須引入湍流模型。最常見的湍流模型有零方程模型、一方程模型、二方程模型等。二方程模型中主要有k?ε模型和k?ω模型。
k?ε模型是目前應用最為廣泛的湍流模型,其包括3種形式:標準k?ε模型、RNGk?ε模型和可實現(xiàn)的k?ε模型。標準k?ε模型適用范圍廣,計算量適中,而且有較多數(shù)據(jù)積累和較好的精度,但只適合完全湍流的流動過程模擬。RNGk?ε模型和可實現(xiàn)的k?ε模型在強流線彎曲、漩渦和旋轉等復雜流動模擬中顯現(xiàn)出比標準k?ε模型有更好的表現(xiàn)[15]。
由于航行器在大攻角下航行時可能存在比較強的漩渦,故選用可實現(xiàn)的k?ε湍流模型。
在三維CFD仿真過程中,目前比較常見的是四面體、六面體以及棱柱等網(wǎng)格類型。六面體網(wǎng)格收斂性和精度比四面體網(wǎng)格好,但在劃分網(wǎng)格時需要耗費大量的時間和精力;四面體網(wǎng)格前處理時間較短,但計算時間和計算結果的精度上又不盡如人意。
多面體網(wǎng)格的出現(xiàn)使上述問題得以解決,其相比四面體網(wǎng)格來說,最大優(yōu)勢就是能大幅度減少網(wǎng)格的數(shù)量,從而能大大縮短計算時間、降低對電腦性能的要求;同時經過與實驗數(shù)據(jù)對比發(fā)現(xiàn),網(wǎng)格數(shù)量減少并不會對計算結果的精度造成太大影響。
目前多面體網(wǎng)格主要的問題是未經過大量的試驗數(shù)據(jù)驗證。因此本文在使用該方法進行模擬時先進行仿真驗證。
由于本文只研究海空跨域航行器的縱向特性,不分析側向特性,因此在計算時只需計算一半模型即可。同時對該模型進行等比例縮放,得到翼展為550 mm的航行器模型,這樣便能減少分析時的計算量,提高了數(shù)值模擬的運算效率。
本文在劃分網(wǎng)格時采用重疊網(wǎng)格方法,在背景網(wǎng)格區(qū)域與重疊網(wǎng)格區(qū)域分別劃分網(wǎng)格。所設置的背景網(wǎng)格區(qū)域為4.5 m×1 m ×1.7 m的長方體,其中背景網(wǎng)格區(qū)域的右邊距離航行器首部1 m、左邊距離航行器尾部3 m、前部的對稱面與航行器的對稱面重合、上部距離航行器垂直尾翼翼梢0.8 m、下部距離航行器底部0.7 m、后部距離航行器機翼翼梢0.725 m。所設置的背景網(wǎng)格為0.8 m× 0.4×0.7 m的長方體,其中重疊網(wǎng)格區(qū)域前部的對稱面與航行器的對稱面重合、右邊距離航行器前端0.15 m,左邊距離航行器尾部0.15 m,上部壁面距離航行器垂直尾翼翼梢0.2 m、下部壁面距離航行器底部0.3 m、后部壁面距離航行器機翼翼梢0.125 m。同時重疊網(wǎng)格區(qū)域附近還有1.6 m×0.6 m×1.2 m的長方體網(wǎng)格加密區(qū)域,這樣能使得重疊網(wǎng)格區(qū)域和背景網(wǎng)格區(qū)域間的網(wǎng)格單元長度不至于變化過大,從而影響計算結果。
本文在網(wǎng)格劃分時采用多面體網(wǎng)格類型。該網(wǎng)格類型能夠離散復雜外形區(qū)域,保證了網(wǎng)格的質量,同時也使得網(wǎng)格數(shù)量不至于過大,計算量不至于過高。這樣便能保證生成的網(wǎng)格能準確模擬航行器的運動狀態(tài),同時計算過程也能更迅速。劃分完成后的網(wǎng)格示意圖如圖2所示。
圖2 網(wǎng)格示意圖Fig. 2 The diagram of meshes
將背景網(wǎng)格區(qū)域中的右邊進流面設置為速度入口邊界條件(velocity inlet),并將空中航行的速度分別設置為15 m/s,30 m/s和45 m/s,將水下航行的速度分別設置為1 m/s,3 m/s和5 m/s;將背景網(wǎng)格區(qū)域中的左邊出流面設置為壓力出口邊界條件(pressure-outlet);將背景網(wǎng)格區(qū)域中的對稱面設置為對稱平面邊界條件(symmetry);將背景網(wǎng)格區(qū)域的其余3個面設置為無滑移壁面邊界條件(wall)。
將重疊網(wǎng)格區(qū)域中的對稱面設置為對稱平面邊界條件(symmetry);將重疊網(wǎng)格區(qū)域的其余5個面設置為重疊網(wǎng)格邊界條件(overset mesh);將??湛缬蚝叫衅靼肽P偷谋砻嬖O置為無滑移壁面邊界條件(wall)。設置好的計算域及邊界條件示意圖如圖3所示。
圖3 計算域及邊界條件示意圖Fig. 3 The diagram of computational domain and boundary conditions
在利用數(shù)值模擬方法對??湛缬蚝叫衅鬟M行氣動和水動力模擬之前,先利用該方法對MK46魚雷在攻角為0°,速度為20 m/s時的水動力特性進行了模擬,得到此時的阻力系數(shù)為0.093 52,這與文獻[16]中試驗值的誤差在可接受范圍。這樣便驗證了本文所述數(shù)值模擬方法的準確性。
在對??湛缬蚝叫衅鞯臍鈩犹匦赃M行研究時,速度分別選取15 m/s,30 m/s和45 m/s,在-8°~24°范圍內每隔2°選取一個攻角,得到共計51種工況時的流場。經過對模擬的結果進行整理,得到了海空跨域航行器半模型在不同工況下的氣動特性曲線如圖4所示,不同工況下的升阻力隨著攻角變化時的曲線如圖5所示。
由圖4(a)可知,??湛缬蚝叫衅髟诳罩酗w行時,隨著攻角的增大,升力系數(shù)也隨之增大,并在12°以下時基本為線性增長;之后升力系數(shù)曲線的斜率逐漸減小,并在20°左右時時達到最大值。隨著攻角繼續(xù)增大,升力系數(shù)又會減小,這主要是由于20°左右時航行器開始失速導致的。同時也可以看出,低速時雷諾數(shù)的變化不會使升力系數(shù)產生明顯變化,3個不同速度下的升力系數(shù)曲線基本重合。
由圖4(b)可知,阻力系數(shù)隨著攻角的增大先增大后減小,在-2°左右時達到最小。同時可以看出,低速時阻力系數(shù)曲線隨雷諾數(shù)變化也不明顯。
圖4 航行器氣動特性曲線Fig. 4 Aerodynamic characteristic curves of the vehicle
圖5 海空跨域航行器的升阻力曲線Fig. 5 Lift and drag curves of the vehicle
由圖4(c)可知,隨著攻角的增大,升阻比會迅速增大,在2°左右時達到最大值,之后又會減小。隨著速度變化,低速時升阻比曲線的發(fā)展趨勢一致,只是不同雷諾數(shù)時會對應不同的最值。
對于翼展為550 mm的航行器模型來說,與翼展為1 200 mm的航行器實物相比,其縮尺比為α=2.182,由動力學相似定律可知模型的重力應該是實物重力的1/α3倍。由于實物重量為3.2 kg,則可知對應模型的重量應該為0.31 kg。
通過以上分析,同時結合圖5(a)發(fā)現(xiàn),該航行器只需以較小的速度和攻角航行,產生的升力即可克服重力。同時,結合圖5(b)可知此時阻力非常小,可以很容易滿足航行器的阻力要求。
如圖6所示為航行器在45 m/s,攻角分別為0°和20°時航行器表面壓力云圖及對稱面的速度云圖。由圖6可知,航行器首部、機翼前端處壓力均較大,在設計時應特別注意。
與氣動分析類似,本文又對航行器進行了水動力分析。經過對模擬的結果進行整理,得到了海空跨域航行器在不同工況下的水動力特性曲線。由于與航行器的氣動特性曲線相似,故本文僅列出不同工況下的升阻力隨著攻角變化時的曲線,如圖7所示。
通過分析可知,模型重量約為0.32 kg。同時計算可知,模型的浮力應該為9.82 N,模型具有較大的凈浮力。
結合圖7(a)發(fā)現(xiàn),該航行器只需以較小的負攻角航行,其產生的向下的升力即可克服凈浮力,從而使得航行器可以實現(xiàn)水下航行。同時,結合圖7(b)可知此時阻力非常小,可以很容易滿足航行器的阻力要求。
圖6 ??湛缬蚝叫衅鞯牧鲌鰣DFig. 6 Flow field diagram of the vehicle
圖7 海空跨域航行器的升阻力曲線Fig. 7 Lift and drag curves of the vehicle
本文提出一種可實現(xiàn)水空兩棲作業(yè)且能多次自由穿越水氣界面的航行器—??湛缬蚝叫衅?,通過STAR-CCM+軟件對其進行數(shù)值模擬,分別得到了該航行器在水下和空中航行時各工況下的升阻力及升阻力系數(shù)。通過計算分析得出如下結論:
1)該航行器在空中/水下航行時,其產生的升阻力可以分別滿足空中/水下的航行要求,也就是參考空中飛行器設計的??湛缬蚝叫衅骺梢詫崿F(xiàn)水空作業(yè),可為海空跨域航行器的后續(xù)研究提供參考。
2)多面體網(wǎng)格的運用可以提高計算效率,同時計算精度也可以得到保證。