切向流場影響下激光壁面加載特性

嚴向峰， 杜英杰

（成都航空職業技術學院 航空維修工程系，成都 610100）

0 引言

激光加載具有極高的過熱度，屬高度不平衡過程。在亞聲速飛行物目標抗激光防護技術的研究過程中發現，物體近壁區速度附面層和溫度附面層結構復雜，附面層中輸運現象旺盛，與激光輻照作用中的傳熱、傳質過程以及多相流強烈耦合，影響激光對壁面材料和結構的加熱、熔化、氣化、力學效應以及等離子現象。在前人的研究中［1-3］，鮮有考慮切向流場對于激光壁面加載機理的影響。本文將在激光壁面加載模型中引入亞音速切向流場，并將仿真計算持續至靶材表面出現熔池階段。通過對熔池發展的數據分析切向流場對激光壁面加載特性的影響，揭示存在切向氣流條件下的激光輻照效應機制。

1 模型建立

本文中設置的激光輻照區域為半徑5mm的圓面，且在計算進行過程中激光對應于靶材表面的加載位置固定不變。三維模型如圖1所示。其流場進口處氣流參數（壓力、靜溫與總溫）預先在更大的模型中計算得到，采用profile文件導入圖1所示的壓力進口。

模型中，速度為0.3倍音速的氣流通過壓力進口進入30 mm×30 mm×10 mm的氣相區域，流經輻照表面后通過壓力出口離開計算域。由于激光能量加載采用高斯分布，其光斑最大直徑處加載能量不及光斑中心處的十分之一，可假設激光輻照作用對加載面以外的部分影響極小。因此在計算模型中僅保留了激光光斑輻照區域所在的16 mm×16 mm×10 mm固相區域，其區域邊界距離激光光斑邊緣最近距離為3 mm。

圖1 計算模型

2 計算參數

計算模型中選擇的靶面材料為鋼鐵，在選擇其物性參數時部分參數根據鋼鐵中含碳量的不同略有差異。該材料常溫下的密度為6 957.5 kg/m3，熔點設置為1 630 K，熔化潛熱2.69×105J/kg，熔化后生成熔融金屬密度為6×10-3kg/（m·s）。模型中進口風速為0.3倍音速，可認為空氣是不可壓縮的。根據計算中的溫度變化范圍，空氣及靶材的部分熱物理參數變化如表1。

表1 材料熱物理性質隨溫度的變化

仿真模擬中使用UDF將激光輻照能量轉化為面熱源形式加載于靶材輻照區域表面，參考李勇等［4］的研究，結合實際試驗所采用的激光中心點能量強度，計算得輻照區域中心點處面熱源約為1.98×107W/m2。

3 結果分析

當模型中未引入切向氣流前，激光能量固定加載于靶材輻照表面，經過總時長為0.58 s的計算后得到此時固相材料熔池的發展狀況。如圖2所示，模型對稱剖面的溫度場及熔化邊界（每一個小網格均為邊長0.5 mm的正六面體），形成中心對稱，半徑4 mm的熔池。計算截止時，輻照中心處熔融金屬的最高溫度為2 538 K，緊靠靶材上表面輻照區域的空氣受熱明顯，氣液兩相之間換熱旺盛。而金屬物質間導熱效率明顯不及氣流的換熱效率，其熔池邊界正下方固相材料網格溫度尚不足750 K。

考察輻照表面溫度場分布后可知（如圖3所示），激光輻照區域內存在較大溫度梯度，但一旦離開此區域后，溫度梯度將迅速下降為零，可認定激光熱作用區域僅限于激光光斑輻照區域，因此在模型建立階段，采用縮減固相區域體積，降低模型網格總數，提升計算效率的方案是可行的。

當引入切向流場后，熔融金屬在氣動力作用下發生流動，出現了隨時間變化的兩相界面，激光能量輻照效應將受到兩相流動的影響。

圖2 熔池及溫度場參數

圖3 壁面溫度場分布

在氣流速度的3個分量中，Y軸分量占總速度的98%以上，直接影響流動，因此可用Vy在法線方向的速度梯度來考察氣流邊界層厚度的變化。如圖4所示，氣流在輻照區域及其后的流動過程中速度梯度保持一個相對小值，形成穩定發展的附面邊界層流場。隨著流動邊界層的形成與擴展，降低了氣流與壁面之間由于對流換熱形成的熱損失，有利于激光加載。圖5顯示氣流流經輻照區域時，近壁流場的溫度梯度出現劇烈變化。在接近輻照中心處，壁溫遠大于近壁氣流溫度。當氣流通過輻照中心后，壁面溫度與氣流溫度間的差距逐漸減小。在輻照區域末端，氣流溫度已高于壁面溫度，此時熱量將從流場反向傳遞至靶材表面。

在上述作用的影響下，靶材表面熔池邊緣的發展也呈現出明顯的非對稱性。在圖6～圖8中，記錄了沿氣流方向上熔池半徑隨加載過程進行的發展及熔池深度的變化過程。

在靶材熔化前及熔化前期，切向氣流所造成的靶面輻照區域前后緣熱環境不對稱對熔池發展起主導作用。輻照區域后緣部分由于流動邊界層增大，對流換熱效率降低，且流經此處的氣流溫度在輻照初期甚至高于壁面溫度的多重影響之下，靶材表面溫度達到熔點以及后續熔化邊界的推移速度都要快于前緣。

但在圖8中可以發現，隨著能量加載的持續進行，熔化深入靶材金屬材料內部，熔池發展開始呈現一種相反的特點。隨著出現少量熔融金屬滯留于已形成的熔池后半端，其對于激光輻照的隔絕作用將遠遠強于氣流的加熱效果，此時熔化邊界向前緣的推移速度將高于向后緣的推移速度。

在與無切流場影響下所得的模擬結果做對比的過程中發現，熔池向靶材內部的發展在輻照初期（0.4 s內）呈現明顯的加速趨勢。當附面層氣流不存在時，經0.58 s的輻照時間，熔池深度為1 mm。當加入附面層氣流的影響后，熔池深度由0發展至0.5 mm時用時0.248 s，由0.5 mm發展至1 mm用時0.208 25 s，由1 mm發展至1.5 mm時僅用時0.065 s。由于激光熱效應的淺表特性，熔融金屬覆蓋于輻照表面所產生的保護作用對熔池形狀的發展有較大影響。當存在切向氣流作用時，淺表熔池中的大多數熔融金屬很快被氣流帶離輻照區域，熔化向靶材內部發展的速度相對于無風模型有所增加，輻照區域中出現的凹陷也有利于激光能量的加載。在激光加載過程中，若能夠通過一個短暫預脈沖致使靶材表面出現淺表熔化，那么對于隨后而來的主脈沖加載無疑十分有利，這一特點能夠應用于激光武器防護性能研究。

4 結論

本文采用計算流體力學方法分析了切向流場對激光壁面加載特性的影響：

1）激光加載效應具備明顯的淺表特性，雖然金屬具有較大的熱導率，但仍然可以認為激光加載的熱作用僅限于輻照光斑所在的位置。

2）切向流場受激光輻照加熱及兩相流流動的多重影響，其邊界層厚度在輻照區域穩定增加，并伴隨著巨大的溫度梯度變化，從對流換熱的角度考慮，這種變化降低了對流換熱效率，有利于激光加載。

3）當加入切向流場的影響后，熔池的發展趨勢與無風模型有著顯著區別。在對流換熱效率沿流向不斷下降的前提條件下，輻照區域后緣壁面材料在初期更早達到熔點。但隨著加載過程的持續進行，熔化邊界逐漸深入至靶材內部，駐留于熔池后部的熔融金屬對激光輻照能量的遮蔽效果將導致熔池向輻照前緣的發展速度增加。

4）熔化深入靶材內部的速度在整個能量加載過程中呈現不斷加速的趨勢，即當輻照能量直接深入材料內部時，其熱效應將大大增強。

圖4 附面層速度梯度

圖5 附面層溫度梯度

圖6 熔池發展趨勢（熔池深度0.5 mm）

圖7 熔池發展趨勢（熔池深度1 mm）

圖8 熔池發展趨勢（熔池深度1.5 mm）

［1］ 王偉平，劉常齡，王春彥，等.切向氣流對激光加熱材料的影響［J］.強激光與粒子束，1996，8（3）：373-377.

［2］ 張健，黃晨光.外部流場對激光加熱運動目標影響的數值模擬［J］.強激光與粒子束，2007，19（11）：1817-1821.

［3］ 楊立新，彭曉峰，王補宣.激光加工溶池熔化和凝固過程數值分析［J］.工程熱物理學報，2001，22（增刊 1）：133-136.

［4］ 李勇，汪民樂，張均.高能激光武器對彈道導彈毀傷能力研究［J］.紅外與激光工程，2006，10，35（5）：588-592.