徐路程，肖凱濤

?

基于CFD方法的紅外煙幕干擾性能研究

徐路程，肖凱濤

（防化研究院，北京 102205）

采用計(jì)算流體力學(xué)方法對(duì)一平坦開闊地域的風(fēng)場(chǎng)進(jìn)行數(shù)值模擬，運(yùn)用離散相模型對(duì)紅外煙幕在該區(qū)域的擴(kuò)散進(jìn)行計(jì)算，得到了紅外煙幕質(zhì)量濃度的三維空間分布，并由Lambert-Beer定律計(jì)算出紅外煙幕在跨風(fēng)方向上的有效遮蔽區(qū)域，研究了風(fēng)速對(duì)于紅外煙幕遮蔽區(qū)域的影響規(guī)律。

紅外煙幕；計(jì)算流體力學(xué)；離散相模型；Lambert-Beer定律

0 引言

紅外偵察和紅外制導(dǎo)由于其分辨率高和抗干擾能力強(qiáng)而倍受各國(guó)重視，工作在紅外波段的偵察和制導(dǎo)裝備迅速發(fā)展[1]。煙幕作為對(duì)抗精確制導(dǎo)武器和觀瞄器材的無(wú)源干擾手段在國(guó)內(nèi)外已有了很大發(fā)展[2]。按照遮蔽機(jī)制，煙幕可以分為輻射遮蔽型和衰減遮蔽型，目前的研究主要集中在衰減遮蔽型煙幕上[3]。這種煙幕主要靠煙幕對(duì)紅外的吸收和散射作用來(lái)實(shí)現(xiàn)遮蔽。煙幕的仿真對(duì)于戰(zhàn)術(shù)使用有較強(qiáng)的參考價(jià)值，文獻(xiàn)[3-4]分別用無(wú)風(fēng)的高斯模型和拉赫特曼擴(kuò)散理論計(jì)算了三維空間的濃度場(chǎng)，并以此為基礎(chǔ)利用Lambert-Beer定律計(jì)算了空間透光率。高斯模型中也可以包含風(fēng)速的影響，推導(dǎo)過(guò)程中用一點(diǎn)的風(fēng)速代表整個(gè)擴(kuò)散空間中的風(fēng)速，即認(rèn)為三維空間中的風(fēng)速分布是均勻的，但由于大氣邊界層中下墊面的作用，風(fēng)速在垂直方向上存在梯度；拉赫特曼擴(kuò)散理論中考慮了風(fēng)速梯度的影響，但適合于模擬瞬時(shí)體源的擴(kuò)散過(guò)程，不適用于連續(xù)源的模擬。本文將首先運(yùn)用計(jì)算流體力學(xué)方法對(duì)一平坦開闊地域的風(fēng)場(chǎng)進(jìn)行數(shù)值模擬，之后在計(jì)算得到的風(fēng)場(chǎng)的基礎(chǔ)上運(yùn)用離散相模型模擬某型紅外煙幕在該區(qū)域的擴(kuò)散，利用Lambert-Beer定律可以計(jì)算通過(guò)擴(kuò)散實(shí)現(xiàn)的三維空間分布的質(zhì)量濃度在跨風(fēng)方向上的透光率，進(jìn)而確定煙幕遮蔽的有效區(qū)域。利用以上方法，通過(guò)模擬不同風(fēng)速條件下的風(fēng)場(chǎng)，從而研究風(fēng)速對(duì)于紅外煙幕干擾性能的影響規(guī)律。

1 風(fēng)場(chǎng)模擬

1.1 幾何模型

風(fēng)場(chǎng)數(shù)值模擬的幾何模型如圖1所示，計(jì)算域?yàn)?000m×400m的長(zhǎng)方形區(qū)域，模擬的風(fēng)向?yàn)槲黠L(fēng)。由于氣溶膠擴(kuò)散在空間中是一個(gè)三維過(guò)程，并且擴(kuò)散過(guò)程一般都發(fā)生在大氣邊界層（一般指距離地面200m以內(nèi)的范圍[5]）中。因此，本模型中考慮垂直方向的計(jì)算高度為50m。于是，三維的計(jì)算區(qū)域?yàn)?000m×400m×40m，如圖2所示。

圖1 計(jì)算區(qū)域示意圖

圖2 三維幾何模型及邊界條件

1.2 網(wǎng)格劃分

垂直方向上，下邊界初始網(wǎng)格設(shè)置為0.5m，增長(zhǎng)比率1.2，總網(wǎng)格數(shù)為30；水平方向上，網(wǎng)格均勻分布，網(wǎng)格尺寸4m。本文使用前處理軟件ANSYS ICEM CFD采用計(jì)算精度更高的六面體結(jié)構(gòu)化網(wǎng)格劃分，計(jì)算域內(nèi)共劃分75萬(wàn)網(wǎng)格，網(wǎng)格及坐標(biāo)設(shè)置如圖3所示。

圖3 網(wǎng)格劃分示意圖

1.3 邊界條件

本文中考慮模擬以下情況：參考高度取為ref＝2m，參考風(fēng)速為ref，具體設(shè)置分別為：1.2m/s，2.0m/s，3.0m/s，4.0m/s，5.0m/s，6.0m/s，7.0m/s。根據(jù)文獻(xiàn)[6]，入口的風(fēng)速廓線可以表示為：

式中：＝0.42為卡曼常數(shù)；0＝0.01m為平坦地面的一種地面粗糙長(zhǎng)度；*為摩擦速度。對(duì)于本文中考慮的近地面層，動(dòng)量垂直通量隨高度的變化很小，可以認(rèn)為通量為常數(shù)[7]，即：

式中：表示湍流切應(yīng)力，根據(jù)文獻(xiàn)[6]，可以如下估算摩擦速度：

將各參數(shù)帶入上式中，可以得到*，將此式帶入式(1)即可確定速度廓線。常溫常壓條件下，空氣密度為＝1.225kg/m3，由式(2)可以估算近地面層湍流切應(yīng)力。

對(duì)于近地面層中的流動(dòng)模擬，應(yīng)該能夠滿足流動(dòng)方向上的均勻性，即保證入口的風(fēng)速廓線式(1)能夠在流動(dòng)過(guò)程得到保持。對(duì)于使用標(biāo)準(zhǔn)-模型的RANS模擬，已經(jīng)有學(xué)者提出了多種方法來(lái)保證水平均勻的速度廓線。最為常用的方法是基于Richard和Hoxey提出的邊界條件的方法[6]，他們建議速度、湍動(dòng)能、湍流耗散率都使用對(duì)數(shù)廓線表示：

1）入口邊界條件：采用速度入口（velocity-inlet）邊界，速度、湍動(dòng)能、湍動(dòng)能耗散率廓線采用式(1)，式(4)，式(5)的定義；

2）出口邊界條件：采用自由流（outflow）邊界，這種邊界假定在這個(gè)邊界上沒(méi)有流向的梯度[5]；

3）下邊界邊界條件：采用切應(yīng)力邊界條件，切應(yīng)力的形式如式(2)所示，文獻(xiàn)[5]指出這種方法能夠獲得非常小的入-出口誤差；

4）上邊界邊界條件：采用速度邊界條件，速度、湍動(dòng)能、湍動(dòng)能耗散率的取值按照式(1)，式(4)，式(5)在上邊界高度處的函數(shù)值來(lái)確定，從而保持物理量的連續(xù)性，相比于常見于其他文獻(xiàn)中用于上邊界的對(duì)稱邊界和滑移邊界，這種邊界能夠最少的引入流向上的梯度[5]，從而能夠更好的保證流向上的均勻性；

5）側(cè)邊界邊界條件：采用速度邊界條件，速度、湍動(dòng)能、湍動(dòng)能耗散率廓線采用與入口邊界相同的形式，方向沿流動(dòng)方向。

1.4 流場(chǎng)求解

本文使用Ansys Fluent 12.0對(duì)控制方程進(jìn)行求解。

1）控制方程離散

Fluent對(duì)控制方程采用有限體積法進(jìn)行離散。動(dòng)量方程、能量方程、湍動(dòng)能方程和湍流耗散率方程使用二階迎風(fēng)格式進(jìn)行離散；壓力使用body-force-weighted方法進(jìn)行離散。

2）壓力速度耦合方式

本文中選擇SIMPLEC算法對(duì)壓力速度進(jìn)行耦合。

3）殘差設(shè)定

為了保證計(jì)算結(jié)果的收斂性，將能量方程的收斂準(zhǔn)則設(shè)定為10－6（標(biāo)準(zhǔn)化殘差），其他方程的收斂準(zhǔn)則設(shè)定為10－3（標(biāo)準(zhǔn)化殘差）。

2 擴(kuò)散相模擬

紅外煙幕的施放源的尺寸相比整個(gè)計(jì)算域的尺寸可以忽略不計(jì)，因此研究中考慮使用離散相模型中的點(diǎn)源來(lái)模擬施放源，點(diǎn)源的屬性如下：釋放點(diǎn)的水平位置如圖1所示，施放點(diǎn)的高度高于地面2m，空間坐標(biāo)為(100, 200, 2)；施放方向?yàn)檠刂L(fēng)速方向；施放時(shí)間為0s到90s；粒子粒徑為10mm；錐形點(diǎn)源的半角設(shè)定為15°；源強(qiáng)為0.15kg/s。

擴(kuò)散模擬的總時(shí)間為90s，時(shí)間步長(zhǎng)取為1s，共90個(gè)時(shí)間步。

3 透光率計(jì)算

當(dāng)紅外輻射通過(guò)煙幕氣溶膠時(shí)，被煙幕氣溶膠所散射和吸收而消弱，其規(guī)律符合Lambert-Beer定律，即：

式中：為通過(guò)煙幕后的輻射強(qiáng)度，W/cm2；0為通過(guò)煙幕前的輻射強(qiáng)度，W/cm2；為煙幕濃度，g/m3；為紅外輻射通過(guò)濃度的光程，m；e為消光系數(shù)，m2/g，它與發(fā)煙劑的性質(zhì)、煙幕粒子的大小、輻射波等有關(guān)，經(jīng)煙箱測(cè)定，本文中使用的紅外煙幕對(duì)于中紅外和遠(yuǎn)紅外的消光系數(shù)分別為1.4m2/g和1.3m2/g，本文中對(duì)中紅外進(jìn)行模擬。

通過(guò)一段大氣路徑的透光率定義為前后輻射通量密度之比[7]，即：

對(duì)于煙幕氣溶膠非均勻分布的情況，透光率可以表示為：

計(jì)算流體力學(xué)計(jì)算得到的濃度場(chǎng)不是顯式表示的函數(shù)，而是離散的空間分布，積分可采用數(shù)值積分的方法進(jìn)行：

本文對(duì)跨風(fēng)方向的透光率進(jìn)行計(jì)算，于是積分方向?yàn)檠刈鴺?biāo)系正軸方向，積分范圍為[100,300]，積分步長(zhǎng)取為D＝0.5m。

4 結(jié)果統(tǒng)計(jì)與分析

在垂直于軸方向的平面上，選取方向[50, 1000]×方向[0, 20]取間距為1m的共951×21個(gè)點(diǎn)，利用FLUENT中編寫journal程序的功能對(duì)∈[100, 300]內(nèi)每0.5m一個(gè)步長(zhǎng)的數(shù)據(jù)點(diǎn)的濃度進(jìn)行統(tǒng)計(jì)并輸出，之后按照式(9)進(jìn)行線積分，得到各離散點(diǎn)處的透光率。運(yùn)用Surfer軟件繪制透光率為0.15的等值線區(qū)域，并測(cè)量有效遮蔽區(qū)域（透光率小于0.15）的面積和長(zhǎng)、寬尺寸，各種風(fēng)速條件下的有效遮蔽區(qū)域如圖4～圖10所示，有效遮蔽區(qū)域的長(zhǎng)度、寬度、面積與風(fēng)速的關(guān)系如圖11～圖13所示。

如圖4～圖10所示，在較短的施放時(shí)間（30s）后，紅外煙幕的遮蔽區(qū)域在空間中連續(xù)、完整，并且煙幕的有效長(zhǎng)度能夠隨著風(fēng)速的增加而逐漸增長(zhǎng)；在稍長(zhǎng)的施放時(shí)間（60s, 90s）后，適中的風(fēng)速（1.2m/s～4m/s）條件下有效長(zhǎng)度仍能夠較穩(wěn)定增長(zhǎng)，但較大風(fēng)速（5m/s～7m/s）條件下，在煙幕施放的遠(yuǎn)端達(dá)到閾值透光率的等值線已經(jīng)支離破碎，難以達(dá)到遮蔽的效果，有效長(zhǎng)度明顯下降。風(fēng)場(chǎng)引起的輸運(yùn)和擴(kuò)散是煙幕形成有效遮蔽區(qū)域的主要?jiǎng)恿?lái)源，風(fēng)速的加強(qiáng)一方面使煙幕氣溶膠在相同的時(shí)間內(nèi)運(yùn)動(dòng)到更遠(yuǎn)的區(qū)域，但同時(shí)也使其空間濃度分布更為分散，而湍流的存在會(huì)加強(qiáng)這種分散作用，因此在風(fēng)速增大的過(guò)程中會(huì)出現(xiàn)有效遮蔽區(qū)域長(zhǎng)度先增后降的規(guī)律。

煙幕的有效遮蔽區(qū)域高度在各個(gè)時(shí)間段內(nèi)都表現(xiàn)出隨風(fēng)速增加先遞增再遞減的規(guī)律，高度的峰值出現(xiàn)在風(fēng)速在2m/s～3m/s的條件，這種現(xiàn)象產(chǎn)生的主要原因是擴(kuò)散作用：風(fēng)速較低時(shí)垂直方向上的速度脈動(dòng)量較小，煙幕聚集在相對(duì)較低的區(qū)域，導(dǎo)致有效遮蔽高度較低，而隨著風(fēng)速的增加，垂直方向上的脈動(dòng)速度也有所增加，這導(dǎo)致了一方面煙幕由于更大的風(fēng)速快速地向下風(fēng)向運(yùn)動(dòng)，另一方面垂直向上的輸送也有所增加，使有效遮蔽高度增加，但過(guò)大的風(fēng)速使煙幕較為分散難以達(dá)到有效線積分濃度，于是有效遮蔽高度下降。

圖4 風(fēng)速為1.2m/s時(shí)3個(gè)時(shí)刻煙幕遮蔽的有效區(qū)域

圖5 風(fēng)速為2.0m/s時(shí)3個(gè)時(shí)刻煙幕遮蔽的有效區(qū)域

Fig.5 Effective area of smoke shielding at 3 moments (wind velocity＝2.0m/s)

圖6 風(fēng)速為3.0m/s時(shí)3個(gè)時(shí)刻煙幕遮蔽的有效區(qū)域

圖7 風(fēng)速為4.0m/s時(shí)3個(gè)時(shí)刻煙幕遮蔽的有效區(qū)域

Fig.7 Effective area of smoke shielding at 3 moments (wind velocity＝4.0m/s)

圖8 風(fēng)速為5.0m/s時(shí)3個(gè)時(shí)刻煙幕遮蔽的有效區(qū)域

圖9 風(fēng)速為6.0m/s時(shí)3個(gè)時(shí)刻煙幕遮蔽的有效區(qū)域

圖10 風(fēng)速為7.0m/s時(shí)3個(gè)時(shí)刻煙幕遮蔽的有效區(qū)域

圖11 有效遮蔽區(qū)域長(zhǎng)度與風(fēng)速關(guān)系

and wind velocity

圖12 有效遮蔽區(qū)域高度與風(fēng)速關(guān)系

圖13 有效遮蔽區(qū)域面積與風(fēng)速關(guān)系

Fig.13 Relationship between acreage of effective shielding area and wind velocity

風(fēng)速對(duì)有效遮蔽區(qū)域長(zhǎng)度和高度的影響綜合體現(xiàn)在對(duì)于面積的影響上：風(fēng)速對(duì)于面積的影響同樣表現(xiàn)出先遞增再遞減的規(guī)律，但由于長(zhǎng)度和高度峰值對(duì)應(yīng)的風(fēng)速不同，有效面積的峰值出現(xiàn)在風(fēng)速為3m/s～4m/s的條件下。

綜上所述，本文中研究的發(fā)煙劑對(duì)于中紅外波段，在3m/s～4m/s的風(fēng)速條件下能夠形成最為穩(wěn)定、遮蔽效果最為理想的紅外煙幕。

5 結(jié)論

本文提出了一種將計(jì)算流體力學(xué)方法、離散相模型、Lambert-Beer定律相結(jié)合的紅外煙幕干擾性能的數(shù)值模擬方法；運(yùn)用以上方法，對(duì)多種風(fēng)速條件下紅外煙幕擴(kuò)散進(jìn)行了數(shù)值模擬，得到了風(fēng)速對(duì)于該種煙幕對(duì)中紅外透光率的影響規(guī)律。該種方法為煙幕干擾性能的數(shù)值模擬提供了一種新的思路，對(duì)煙幕技術(shù)的實(shí)際使用具有指導(dǎo)意義。主要結(jié)論如下：

1）煙幕的有效遮蔽長(zhǎng)度在適中風(fēng)速（2m/s～4m/s）條件下能夠較穩(wěn)定增長(zhǎng)，較大風(fēng)速（5m/s～7m/s）條件下，有效長(zhǎng)度明顯下降；

2）煙幕的有效遮蔽高度隨風(fēng)速增加先遞增、再遞減，高度峰值對(duì)應(yīng)風(fēng)速為2m/s～3m/s；

3）煙幕的有效遮蔽面積隨風(fēng)速增加先遞增、再遞減，面積峰值對(duì)應(yīng)風(fēng)速為3m/s～4m/s；

4）綜合考慮煙幕有效遮蔽區(qū)域的長(zhǎng)度、高度、面積，本文中研究的發(fā)煙劑對(duì)于中紅外波段在3m/s～4m/s的風(fēng)速條件下使用能夠形成最為穩(wěn)定、最為理想的遮蔽效果。

[1] 肖凱濤, 宋偉偉, 何友金. 基于海陸交界地區(qū)特征的紅外煙幕干擾性能研究[J]. 紅外技術(shù), 2013, 35(2): 113-118.

[2] 姚祿玖, 高鈞麟, 肖凱濤. 煙幕理論與測(cè)試技術(shù)[M]. 北京: 國(guó)防工業(yè)出版社, 2004.

[3] 張楠, 張科, 王紅梅, 等.一種紅外煙幕三維建模仿真方法[J]. 電光與控制, 2014, 21(6): 48-51.

[4] 邱繼進(jìn), 趙曉哲. 紅外煙幕干擾效果的計(jì)算與模擬[J]. 電光與控制, 2006, 13(4): 31-33.

[5] B. Blocken, T. Stathopoulos, J. Carmeliet. CFD simulation of the atmospheric boundary layer:wall function problems[J]., 2007, 41: 238-252.

[6] P. J. Richards, R. P. Hoxey. Appropriate boundary conditions for computational wind engineering models using the k-eturbulence models[J]., 1993, 46-47: 145-153.

[7] 盛裴軒, 毛節(jié)泰, 李建國(guó), 等. 大氣物理學(xué)[M]. 北京: 北京大學(xué)出版社, 2003.

CFD-based Study on Countermeasure Performance of Anti-infrared Smoke Screen

XU Lu-cheng，XIAO Kai-tao

(,102205,)

In this paper, computational fluid dynamics method was applied for the numerical simulation of a flat and open wind field. The dispersion of anti-infrared smoke screen in this area was evaluated by discrete phase model, and the three-dimensional distribution of the concentration of anti-infrared smoke screen was obtained. The efficient shielding region of anti-infrared smoke screen across the wind direction was calculated by using Lambert-Beer law. The relationship between wind velocity and shielding region of anti-infrared smoke screen was also studied in this paper.

anti-infrared smoke screen，computational fluid dynamics，discrete phase model，Lambert-Beer law

O434.3

A

1001-8891(2015)04-0337-05

2014-12-22；

2015-03-13.

徐路程（1990-），男，吉林通化人，碩士研究生，主要從事大氣擴(kuò)散及煙幕干擾特性研究。

總裝備部“十二五”預(yù)先研究課題項(xiàng)目。