- 相關推薦
高速飛行器熱流固耦合光傳輸分析論文
引言
高速飛行器在大氣層內(nèi)高速飛行時，其光學頭罩與來流相互作用產(chǎn)生大量的熱，形成嚴重的氣動光學效應。一方面，光學頭罩外空氣受到突然壓縮和摩擦，其溫度、壓力、密度和化學成分發(fā)生變化，形成高速復雜流場; 另一方面，光學頭罩的溫度和應變量會隨著飛行時間發(fā)生變化，進而引起晶體材料折射率的改變。來自目標的光線到達探測器將受到流場與光學頭罩的雙重影響。出射波面發(fā)生畸變，對成像探測系統(tǒng)造成光傳輸干擾，使目標圖像出現(xiàn)模糊、抖動等現(xiàn)象，嚴重影響探測結(jié)果。因此，除了需事先對紅外成像系統(tǒng)的性能進行測試，還需對高速飛行器氣動光學效應進行全面的研究，為系統(tǒng)的矯正提供數(shù)據(jù)基礎。
目前，對氣動環(huán)境下的光傳輸影響分析主要是針對高速流場或光學頭罩其中一個方面進行的，對其耦合光傳輸?shù)挠绊懛治鲚^少，本文針對這一情況，通過建立氣動環(huán)境下高超聲速飛行器在飛行過程中流場與頭罩的耦合參數(shù)變化模型，計算CFD 網(wǎng)格數(shù)據(jù)，采用光線追跡計算熱流固耦合光傳輸影響。
1 飛行過程的熱流固耦合分析
1. 1 繞流流場的數(shù)值模擬
湍流的形成是由于流體中的慣性力對流體的影響占到主導地位。流體的流動受到物理守恒定律的支配。由于流場結(jié)構(gòu)復雜，無法用統(tǒng)一的方程來描述，為方便研究，建立CFD 網(wǎng)格，在流場中任取正六面體微元作為研究對象，該微元和周圍流體之間存在相互作用力以及質(zhì)量和能量的交換。基于所取微元的任意性，對該微元建立可代表流體的運動方程，稱為Navier - Stokes 方程組:
1.2 窗口的熱流固耦合分析
高速飛行器在飛行時與周圍的空氣劇烈摩擦，動能轉(zhuǎn)化為熱能，導致流場的溫度急劇升高，并加熱頭罩。在熱流固耦合數(shù)值模擬中，繞流流場作用于頭罩，使后者產(chǎn)生溫度和應力應變的動態(tài)變化。因此，光學頭罩在氣動環(huán)境下的加熱為持續(xù)的非穩(wěn)態(tài)過程，在計算過程中應將流場與頭罩進行直接耦合計算。
1. 2. 1 表面壓力場的計算根據(jù)流體力學模型可以假設，當氣體質(zhì)點與頭罩表面碰撞后，氣體沿頭罩表面的法向動量全部損失掉，形成施加于物體上的力，通過計算動量的變化，可以計算出作用在頭罩表面的壓力:
1.2.2 溫度場分析
為保證流場與頭罩之間的氣動對流換熱，頭罩內(nèi)部傳熱主要以傳導和輻射方式進行。在直角坐標系中，溫度分布的基本方程為
2 耦合光傳輸過程
2.1折射率場的計算
光線通過非均勻的折射率場會發(fā)生偏折，并產(chǎn)生附加相位，影響光線傳輸。精確計算折射率場對光線追跡計算光傳輸效應具有重要的意義。根據(jù)之前對繞流流場數(shù)值模擬可以計算出流場的壓力、密度等信息，通過物質(zhì)原子理論和Gladston - Dale 定律計算流場的折射率場分布:
氣動熱環(huán)境下的光學窗口由于受到非均勻的溫度場和應力場影響，材料的密度和極化率發(fā)生改變，進而引起折射率的變化。主要表現(xiàn)為熱光效應和彈光效應。式( 10) 中的第二項表明光學材料的折射率是隨著溫度的變化而產(chǎn)生的變化量，此過程稱為熱光效應。第三項表明由應力應變的作用引起窗口材料的光學性質(zhì)改變，主要是指折射率發(fā)生變化，稱為彈光效應。改變后的折射率表示為:
在熱流固耦合條件下，不同飛行時間下z = 0處的對稱截面上流場與頭罩的折射率場的分布如圖所示。
2. 2 光線追跡法
由于流場與頭罩的折射率場情況復雜并且分布無規(guī)律，本文通過數(shù)值方法求解光線方程實現(xiàn)光線追跡:
為保證計算精度和運算速度滿足工程需要，本文采用四階Runge - Kutta 法計算此光線方程。通過引入外推參量，將方程轉(zhuǎn)化為一階方程組，并以三維表示如下:
2. 3 交界面數(shù)據(jù)處理
基于應力場與溫度場的存在，頭罩在飛行過程中發(fā)生形變，對頭罩表面進行光線追跡時采用最小二乘法對發(fā)生形變后的表面網(wǎng)格進行曲面擬合。由于流體的可壓縮性大，沒有固定的形狀，不能用計算固體材料應變理論計算流體的形狀分布，導致近頭罩處的流場與頭罩膨脹部分的坐標數(shù)據(jù)出現(xiàn)交疊，這在計算耦合光傳輸時是不合理的，因此，需要對該部分的CFD 網(wǎng)格數(shù)據(jù)做進一步計算處理。
現(xiàn)在從理論上解決該問題的方法分為兩種:一是動網(wǎng)格技術(shù)，二是根據(jù)仿真結(jié)果對數(shù)據(jù)做插值擬合處理。動網(wǎng)格技術(shù)可用來模擬流體的形狀隨時間和邊界運動而發(fā)生改變的情況，更加真實地模擬出物體運動過程的變化，目前被用于魚雷等低速武器上。對于高速和高超聲速飛行器的仿真，由于缺乏有力的運動模型支持，迭代仿真參數(shù)難以收斂。本文采用第二種方法，通過數(shù)值手段，對交疊的網(wǎng)格數(shù)據(jù)進行插值，使計算后的流場能夠緊密包裹住光學頭罩。
2. 4 交界面光傳輸處理
當光線在兩種折射率不同的透明介質(zhì)中傳播時，根據(jù)入射情況不同，會發(fā)生折射、反射、全反射等現(xiàn)象。除此之外，在氣動熱環(huán)境下，光學頭罩會發(fā)生不可忽略的形變，因此，精確計算光線在流場與頭罩交界面處的傳播情況是分析耦合光傳輸效應的重要內(nèi)容。
2. 5 像質(zhì)評價指標
通過合理的光學質(zhì)量評價指標可以直觀的描述非均勻折射率場對光束質(zhì)量的影響，為圖像畸變校正提供依據(jù)，本文通過選擇光學傳遞函數(shù)、能量集中度以及斯特列爾比等三項指標來進行像質(zhì)評價，具體計算方法如下:
3 仿真結(jié)果分析
設來流的靜壓為101 325 Pa，靜溫為300 K。飛行器采用半球形硫化鋅頭罩，頭罩厚度為5 mm，飛行高度為20 km，飛行速度分別為3Ma 和5Ma，攻角都為0°。其余的光學特性參數(shù)如 所示。利用光線追跡計算得到不同時刻的耦合光學傳遞函數(shù)和能量集中度情況如 所示。由 可以看出，飛行時間較短時，系統(tǒng)的傳遞函數(shù)較好，此時光學系統(tǒng)與理想情況比較接近，失真程度較小。當飛行速度達到5Ma，飛行時間為120 s 時，4Hz 的目標傳遞函數(shù)已經(jīng)衰減到接近10%，系統(tǒng)幾乎失效。
而從能量集中度可以看出，隨著飛行速度的增大以及飛行時間變長，達到90% 能量時的彌散半徑越來越大，系統(tǒng)失真程度增大。
不同飛行條件下的Strehl 比同樣反映圖像失真情況，如 所示，隨著飛行條件越苛刻，其比值越小，能量分布越分散，探測器接受到的圖像成像質(zhì)量越低。
4 結(jié)論
本文分析了高速飛行器在飛行過程中流場與頭罩的熱流固耦合情況，計算相應條件下的折射率場分布情況。通過該耦合折射率場光程差情況得出相應的光學傳遞函數(shù)、能量集中度以及斯特列爾比，根據(jù)計算結(jié)果可知:
隨著飛行速度增大，耦合折射率場對光傳輸?shù)挠绊戇M一步加大，成像質(zhì)量迅速下降并且隨著飛行條件的惡劣，能量衰減嚴重，彌散半徑增大。
本文所建立的氣動環(huán)境下熱流固耦合的數(shù)值仿真法，更加真實地模擬飛行過程中流場與頭罩產(chǎn)生的非均勻折射率場的變化。通過計算光線在該折射率場中的傳輸情況，直觀地顯示了該動態(tài)過程對成像系統(tǒng)造成的影響。
【高速飛行器熱流固耦合光傳輸分析論文】相關文章：
移動網(wǎng)絡中OTN傳輸技術(shù)的應用論文05-27
利用光纜解決音頻信號傳輸問題的論文04-25
高速公路水泥探討的論文05-22
論文寫作特點分析05-08
復雜網(wǎng)絡上的傳播和耦合動力學過程研究論文提綱06-11
有線電視信號傳輸中的干擾問題論文08-01
對于OMNEST 的以太網(wǎng)傳輸丟包問題研究論文04-27
GPRS問題的分析與解決論文04-27