陳 豪，石 磊，馬麗華，徐 景，盧朝梁，陳永安

（1.空軍工程大學信息與導航學院，陜西西安710077；2.95561部隊，西藏拉薩850002）

納秒激光等離子體減阻數值模擬

陳 豪1，石 磊1，馬麗華1，徐 景2，盧朝梁1，陳永安1

（1.空軍工程大學信息與導航學院，陜西西安710077；2.95561部隊，西藏拉薩850002）

針對激光等離子體減阻技術機理，采用納秒激光，數值模擬激光等離子體在流場中的演化過程，分析關鍵參數對納秒激光等離子體減阻性能的影響。結果表明：納秒激光能最大程度地提高激光等離子體減阻性能，阻力減小的百分比達99%，低阻力維持時間是入射激光持續(xù)時間的103倍；隨著來流馬赫數的增大，空氣來流強度增強，導致減小阻力的百分比減??；隨著激光能量的增加，激光引致的沖擊波強度增大，使得減小阻力的百分比增加；隨著激光聚焦擊穿位置的增大，減小阻力的百分比減小，但低阻力維持的時間明顯增加。

納秒激光；等離子體減阻；減小阻力的百分比；來流馬赫數；激光能量；激光聚焦位置

1 引 言

1978年，美國倫瑟勒工學院的Myrabo率先提出利用激光在飛行器前端聚焦擊穿空氣形成高溫高壓等離子體，減小超聲速飛行器氣動阻力的新概念［1］，其基本原理是在飛行器頭部通過激光等離子體形成高壓擾動，改變超聲速來流形成的弓形脫體激波波形，變?yōu)樾奔げ?，從而減小飛行器受到的阻力［2］。由于具有可控性好，減阻效率高等特點，受到國內外學者的廣泛關注［3］。近年來，基于微秒級CO2氣體激光，相繼開展了減阻機理、實驗和數值模擬研究［4－7］，討論分析了激光脈沖能量［8］、激光聚焦擊穿位置［9］及來流馬赫數［10］對減阻性能的影響，目前利用激光等離子體技術可以減小70%以上的波阻。

隨著激光器水平的發(fā)展，體積更小、脈沖更窄、峰值功率更高的納秒固體激光器在航空航天應用中的優(yōu)勢得以顯現。本文采用納秒激光，開展激光等離子體減阻技術數值模擬，分析納秒激光等離子體減阻性能。

2 計算模型

2.1 計算構型及計算區(qū)域劃分

數值計算中采用2維軸對稱計算模型，使用拋物形飛行器模型，滿足方程y2＋z2＝2px，p＝10為拋物參數，底部直徑為D＝30 mm，高為H＝11.25 mm，如圖1所示。為降低計算區(qū)域復雜度，基于計算分區(qū)思想，將計算區(qū)域劃分為3個區(qū)域，通過Gambit軟件生成計算網格。

圖1 計算模型Fig.1 Computation model

2.2 初邊值條件

數值模擬中邊界條件如圖1所示，其中inflow為強制高超聲速入口邊界條件，wall為滑移固壁邊界條件，outflow為出口邊界條件，axis為對稱軸。計算中假設環(huán)境氣體是距離地球表面30 km高空的大氣，具體參數如表1所示。

表1 環(huán)境氣體參數Tab.1 Parameters of ambient air

2.3 數值模擬方案

不考慮粘性和質量力，控制方程采用Navier-Stokes方程。數值模擬計算基于有限體積法，采用Roe-FDS通量格式，通過基于Van Albada通量限制器的單調逆風格式（MUSCL）方法獲得三階空間精度，時間積分采用二階精度的耦合隱式算法；用離散轉移方法（DTRM）分析等離子體輻射；使用標準κε雙方程模型進行湍流計算。

3 計算結果與分析

3.1 納秒激光等離子體減阻流場演變

首先以自由來流馬赫數Ma＝5.5為例，分析納秒激光等離子體減阻流場演變。計算中采用單脈沖激光能量E＝2J，脈沖寬度tp＝20ns的納秒激光，并假設激光聚焦位置L＝0.6D（L為激光聚焦點到飛行器頭部的距離，D為飛行器的直徑）。圖2顯示了計算時間t＝0～50μs內飛行器的阻力曲線。

圖2 阻力曲線Fig.2 Air drag curve in the case of Ma＝5.5

從圖中可知，來流馬赫數Ma＝5.5條件下，飛行器受到的穩(wěn)態(tài)阻力為57.86N，激光作用后的最小阻力僅有0.53N，減小阻力的百分比ηD達99%，明顯高于微秒激光70%的減阻效果［3］，另一方面低阻力（小于穩(wěn)態(tài)阻力）維持時間達23μs，是入射激光持續(xù)時間20 ns的1150倍。

分析曲線變化趨勢可知，納秒激光經聚焦擊穿后產生的等離子體沖擊波在4.3μs時傳播到飛行器外壁面，并開始產生作用力，導致阻力迅速升高，在5.8μs時達到峰值，最大值為215 N。此后阻力迅速減小，在t＝14μs時下降趨勢變緩，在21μs時達到最小值0.53 N。在t＝22μs時，阻力再次迅速升高，然而在t＝23μs后，納秒激光等離子體的影響幾乎消失，在自由空氣來流作用下，阻力緩慢升高，而后下降，收斂至穩(wěn)態(tài)阻力。

圖3給出了t＝0～50μs之間流場內的壓強等值線。

在t＝4.3μs時，激光等離子體沖擊波剛剛傳播到飛行器外壁面，波形呈非對稱性傳播，前端波陣面（逆空氣來流方向傳播）受空氣來流影響傳播較慢，后端波陣面?zhèn)鞑ポ^快。

在t＝5.8μs時，阻力迅速升高達到最大值。原因是納秒激光引致的沖擊波強波陣面與飛行器相互作用面積增加，使得飛行器外壁面壓力增大，導致出現阻力峰值。

在t＝10μs時，飛行器受到的阻力迅速下降。原因是沖擊波與飛行器頭部弓形激波相互作用，改變激波結構，形成較弱的斜激波，同時伴隨沖擊波波陣面強度的衰減，飛行器受到的阻力迅速減小。

在t＝14μs時，沖擊波傳播到飛行器外側，同時沖擊波經外壁面發(fā)射，在飛行器頭部形成低壓區(qū)，使得阻力下降趨勢變緩。在t＝21μs時達到最小值，僅為0.53 N。

圖3 不同時間壓強等值線Fig.3 Pressure contours at the different time

在t＝22μs時，阻力迅速升高，原因是受空氣來流影響，沖擊波前端波陣面及反射波陣面被吹向下游，并傳播到飛行器外壁面，產生作用力，導致阻力瞬時升高。

在t＝23μs瞬時達到31 N后，伴隨沖擊波強度的衰減，激光引致的沖擊波對飛行器的影響幾乎消失，在自由空氣來流的作用下，阻力緩慢升高，在t＝36μs時達到第二阻力峰值，隨后逐漸減小，趨于穩(wěn)態(tài)。

3.2 來流馬赫數的影響

在相同激光參數及聚焦點位置條件下，計算比較分析來流馬赫數Ma＝5～9范圍內高超音速來流對激光等離子體減阻性能的影響。計算結果如表2所示。

表2 不同來流條件下阻力減小情況Tab.2 Air drag reduction with different inflow

從表中可以看出，隨著來流馬赫數的增大，飛行器受到的穩(wěn)態(tài)阻力快速增大，激光作用后的最小阻力增大，但減小的阻力也快速增大。

圖4給出了不同來流馬赫數條件下飛行器受到的阻力隨時間變化曲線。

圖4 不同來流條件下阻力曲線Fig.4 Air drag curveswith different Mach number

從圖中可知，穩(wěn)態(tài)條件下，飛行器受到的阻力隨著來流馬赫數的增大明顯增大。在納秒激光作用下，隨著來流馬赫數的增大，空氣來流加速等離子體沖擊波吹向下游，使得沖擊波到達飛行器外壁面的時間提前，阻力峰值顯著增加，最小阻力增大。從低阻力維持時間看，隨來流馬赫數的增大，等離子體沖擊波的作用時間減小，導致低阻力維持時間縮短，阻力收斂時間提前。

圖5給出了減小阻力的百分比ηD與來流馬赫數的關系。

圖5 減小阻力的百分比隨來流馬赫數的變化關系Fig.5 The percentage of air drag reduction vs.Mach number

從圖中擬合曲線可以看出，隨著來流馬赫數的增大，阻力減小的百分比呈下降趨勢。

3.3 納秒激光能量的影響

以來流馬赫數Ma＝5.5為例，假設激光聚焦位置L＝0.6D，分析納秒激光能量為0～2J范圍內對減阻性能的影響。圖6給出了不同激光能量條件下，飛行器受到的阻力隨時間變化曲線。

圖6 不同激光能量條件下阻力曲線Fig.6 Air drag curves with different laser energy

從圖中可以看出，不同激光能量條件下，阻力曲線的變化規(guī)律是一致的。區(qū)別在于，隨著激光能量的增加，激光引致的沖擊波強度增大，沖擊波到達飛行器外壁面的時間縮短，阻力峰值增大，阻力最小值減小，低阻力維持的時間延長，阻力收斂時間推遲。

圖7 減小阻力的百分比隨激光能量的變化關系Fig.7 The percentage of air drag reduction vs.laser enengy

圖7 給出了減小阻力的百分比ηD與納秒激光能量之間的變化關系。

從圖中可以看出，在納秒激光能量E＜0.6J時，減小阻力的百分比ηD隨著激光能量的增加而快速增加；而當激光能量E＞0.6J時，減小阻力的百分比ηD達到97%以上，變化緩慢。

3.4 聚焦位置的影響

以來流馬赫數Ma＝5.5為例，采用激光能量為E＝2 J，脈沖寬度為tp＝20 ns的納秒激光，分析激光聚焦擊穿位置對減阻性能的影響。圖8給出了不同點火位置的飛行器受到的阻力隨時間變化曲線。

圖8 不同激光聚焦位置條件下阻力曲線Fig.8 Air drag curveswith different focusing location

從圖中可知，相同來流環(huán)境及激光參數條件下，聚焦位置離飛行器頭部越近，等離子體沖擊

圖9 減小阻力的百分比隨激光聚焦位置的變化Fig.9 The percentage of air drag reduction vs.laser focusing location

波到達飛行器外壁面的時間越早，達到的阻力峰值越大，最小值越小。但從低阻力持續(xù)時間上看，聚焦位置越遠，低阻力維持時間越長。

圖9給出了減小阻力的百分比隨點火位置的變化規(guī)律。

從圖中擬合曲線可知，隨著激光聚焦位置的增大，減小阻力的百分比ηD減小，當L／D＜0.8時，ηD變化較為緩慢；當L／D＞0.8時，ηD快速下降。

4 結 論

以實驗室已有的飛行器為例，采用納秒激光，開展激光等離子體減阻數值模擬仿真，簡要分析了來流馬赫數、脈沖激光能量和激光聚焦點對減小阻力的百分比的影響。主要結論如下：

（1）納秒激光等離子體能夠很大程度提高減小阻力的百分比，在來流馬赫數為Ma＝5.5條件下，達99.08%，并能夠增加低阻力維持的時間。

（2）隨著來流馬赫數的增大，空氣來流對飛行器的阻力增大，對沖擊波的作用增大，導致阻力最大值升高，將沖擊波加速吹向下游，使得低阻力維持時間縮短，減小阻力的百分比減小。

（3）隨著納秒激光能量的增加，激光引致的沖擊波強度增大，阻力最小值減小，減小阻力的百分比ηD增加，當E＜0.6J時，ηD迅速增加，而當E＞0.6J時，ηD達97%以上，變化緩慢。

（4）隨著聚焦位置的增大，減小阻力的百分比ηD減小，但低阻力維持的時間明顯增加。

［1］ Myrabo.Solar-Powered global air transportation［R］. AIAA，1978：78－0689.

［2］ Sun Zhongxiang.Progress in plasma assisted drag reduction technology［J］.Advances in Mechanics，2003，33（1）：87－94.（in Chinese）孫宗祥.等離子體減阻技術的研究進展［J］.力學進展，2003，33（1）：87－94.

［3］ Hong Yanji，Li Qian，Fang Juan，et al.Advances in study of laser plasma drag reduction technology［J］.Acta aeronautica et astronautica sinica，2010，31（1）：93－101.（in Chinese）洪延姬，李倩，方娟，等.激光等離子體減阻技術研究進展［J］.航空學報，2010，31（1）：93－101.

［4］ Fang Juan，Hong Yanji，LiQian.Numerical analysis of interaction between single-pulse laser-induced plasma and bow shock in a supersonic flow［J］.Plasma Science and Technology，2012，14（8）：741－746.

［5］ David Sperber，Fabian Schmid，Hans-Albert Eckel.Wave drag reduction of blunt bodies using laser sustained energy deposition in argon atmosphere［C］.6thAIAA Flow Control Conference，2012：2012－2815.

［6］ HyoungJin Lee，InSeuck Jeung，SangHun Lee，et al.Experiment of flow control using laser energy deposition around high speed propulsion system［C］.Beamed Energy Propulsion AIPConf.Proc.，2011，1402：416－423.

［7］ Fang Juan，Hong Yanji，Li Qian，et al.The influence of flight altitude on supersonic drag reduction with laser energy depositions［C］.Beamed Energy Propulsion AIP Conf.Proc.，2011，1402：430－436.

［8］ Li Qian，Jin Xing，Cao Zhengrui，et al.Effects on aerodynamic drag of incident laser energy in technology of laser plasma point source drag reduction［J］.Journal of propulsion technology，2010，31（3）：377－380.（in Chinese）李倩，金星，曹正蕊等.激光等離子體點源減阻技術中入射能量對氣動阻力的影響［J］.推進技術，2010，31（3）：377－380.

［9］ Fang Juan，Hong Yanji，Huang Hui，et al.Influence of ignition location on drag reduction by laser plasma［J］. High power laser and particle beams，2010，22（9）：2059－2062.（in Chinese）方娟，洪延姬，黃輝等.點火位置對激光等離子體減阻效能的影響［J］.強激光與粒子束，2010，22（9）：2059－2062.

［10］Liu Zhun，Dou Zhiguo，Huang Hui，et al.The influence of the inflow mach number on laser plasma drag reduction［J］.Development＆Innovation ofMachinery＆Electrical Products，2010，23（4）：18－20.（in Chinese）劉準，竇志國，黃輝，等.來流馬赫數對激光等離子體減阻性能影響［J］.機電產品開發(fā)與創(chuàng)新，2010，23（4）：18－20.

Numerical simulation of nanosecond laser plasma drag reduction

CHEN Hao1，SHILei1，MA Li-hua1，XU Jing2，LU Chao-liang1，CHEN Yong-an1
（1.School of Information and Navigation，Air Force Engineering University，Xi′an 710077，China；2.Army 95561 of PLA，Lhasa 850002，China）

Aiming at the mechanism of laser plasma drag reduction technology，the propagation of the shock wave caused by the nanosecond laser is numerically simulated.The influence of the key parameters on the air drag reduction performance is analyzed.It can be concluded from the simulation results that the performance of the plasma drag reduction is improved greatly by the nanosecond laser，where the percentage of air drag reduction is close to 99%and the time of low air drag is103times of the incident laser time.The percentage of air drag reduction decreaseswith the Mach number increasing because the intensity of the air inflow increased.The percentage increases with the laser energy increasing because the shock wave caused by the nanosecond laser increases.The percentage decreaseswith the focusing location being far，but the time of low air drag obviously increases.

nanosecond laser；plasma drag reduction；percentage of drag reduction；Mach number；laser energy；focusing location

TN249；V439

A

10.3969／j.issn.1001-5078.2014.02.0

1001-5078（2014）02-0131-05

國家自然科學基金資助項目（No.61108068）。

陳 豪（1989－），男，碩士研究生，主要從事激光等離子體減阻技術研究。E-mail：chenhao_8889＠163.com

2013-07-08；

2013-07-27