火藥燃燒生成等離子體規(guī)律數(shù)值仿真研究

宋 鵬，毛保全，鐘孟春，李曉剛，蘭 圖

?

火藥燃燒生成等離子體規(guī)律數(shù)值仿真研究

宋 鵬1，毛保全2，鐘孟春2，李曉剛2，蘭 圖2

（1.中國(guó)華陰兵器試驗(yàn)中心，陜西華陰, 714200；2.裝甲兵工程學(xué)院兵器工程系，北京，100072）

為了研究火藥燃燒生成等離子體規(guī)律，分析了火藥燃燒時(shí)的物理化學(xué)過(guò)程，建立了火藥燃燒時(shí)生成等離子密度規(guī)律數(shù)學(xué)模型，并對(duì)火藥燃燒時(shí)等離子體生成規(guī)律進(jìn)行數(shù)值模擬仿真。將數(shù)值模擬的結(jié)果與試驗(yàn)結(jié)果進(jìn)行對(duì)比，確定模型的正確性與可行性。利用所建立的數(shù)學(xué)模型，分析不同的裝藥結(jié)構(gòu)參數(shù)對(duì)等離子體的影響，得到了等離子體密度變化的規(guī)律，為下一步提高等離子體生成濃度奠定基礎(chǔ)。

火藥；等離子體；燃燒；數(shù)值仿真

等離子體生成的方式主要包括氣體放電法、射線輻照法、光電離法、激光等離子體法、熱電離法、激波等離子體等。氣體放電包括電暈放電、輝光放電、電弧放電、介質(zhì)阻擋放電等，在這一領(lǐng)域的研究比較多。熱電離是指原子在熱運(yùn)動(dòng)中與其它粒子發(fā)生非彈性碰撞獲得足夠大的能量而產(chǎn)生的電離[1-2]，在這方面研究較少，只有關(guān)于化石燃料燃燒生成等離子體的研究。等離子體在火炮中的應(yīng)用現(xiàn)在只限于電熱化學(xué)炮。電熱化學(xué)炮中的等離子體是由等離子體發(fā)生器產(chǎn)生，主要用于火藥點(diǎn)火[3-5]。然而，關(guān)于火藥燃燒時(shí)生成等離子體的性能研究，以及利用火藥燃?xì)猱a(chǎn)生的等離子體改善火炮性能的研究卻未見(jiàn)報(bào)道。

本文以火藥燃燒產(chǎn)生等離子體為研究對(duì)象，火藥燃?xì)饽苓_(dá)到2 000~3 000K，這樣的溫度下會(huì)發(fā)生熱電離，產(chǎn)生等離子體[6]。通過(guò)對(duì)火藥燃燒時(shí)的物理化學(xué)過(guò)程進(jìn)行分析，建立火藥燃燒時(shí)生成等離子密度規(guī)律數(shù)學(xué)模型，然后編寫(xiě)程序進(jìn)行數(shù)值仿真，獲得了火藥燃燒時(shí)等離子體生成情況；并且通過(guò)改變裝藥量、藥室容積和火藥力等裝藥結(jié)構(gòu)參數(shù)，分析了其對(duì)生成等離子體規(guī)律的影響。

1 物理模型

以某型火炮為例，火炮發(fā)射時(shí)，通過(guò)擊針撞擊膛底，點(diǎn)燃底火，通過(guò)底火引燃主裝藥，主裝藥燃燒生成高溫高壓的燃?xì)?，通過(guò)燃?xì)庾龉?，推?dòng)彈丸運(yùn)動(dòng)，最終把彈丸發(fā)射出去[7-8]。根據(jù)火炮發(fā)射時(shí)火藥燃燒情況的特點(diǎn)，做出以下假設(shè)：（1）假設(shè)裝藥符合幾何燃燒規(guī)律；（2）火藥燃燒采用燃速指數(shù)公式，即；（3）假定火藥的燃燒和彈丸發(fā)射都是在平均壓力下進(jìn)行的；（4）所有的次要功用系數(shù)來(lái)計(jì)算，由于射擊過(guò)程中次要功主要是與動(dòng)能有關(guān)，所有的次要功都可以用系數(shù)來(lái)加以考慮；（5）火藥燃燒生成的燃?xì)夥现Z貝爾-阿貝爾狀態(tài)方程；（6）火藥燃燒以后生成的產(chǎn)物成分不變，彈丸由于燃?xì)獾耐屏Χ\(yùn)動(dòng)，燃?xì)庖驗(yàn)橥苿?dòng)彈丸做功而溫度下降，由于溫度的下降引起火藥其他參數(shù)（比如火藥力、余容等）也有略小的下降，但是對(duì)模擬結(jié)果影響微弱，所以把其他參數(shù)看作是常數(shù)。

由于碳?xì)溲醯仍氐碾婋x電位比較高，火藥燃燒的溫度只有3 000K左右，難以使燃燒產(chǎn)物電離。對(duì)于空氣而言，在低于1atm的條件下，需要達(dá)到6 000K以上的溫度，才能具有可觀的等離子體濃度。為了增加電離度，通常采用添加電離種子的方法，既在燃燒物中添加一部分電離電位比較低的物質(zhì)，使燃燒產(chǎn)物在較低的溫度下能獲得較高密度的等離子體。與碳?xì)溲醯啾龋瑝A金屬具有更低的電離電位，所以一般選用堿金屬鹽（比如鉀鹽或者銫鹽等），在燃燒物中添加少量的碳酸鉀，增加燃燒產(chǎn)物的熱電離，獲得一定濃度的等離子體。根據(jù)熱電離理論，現(xiàn)做出如下假設(shè)：（1）碳酸鉀在高溫環(huán)境中完全分解；（2）假定燃燒產(chǎn)物是均勻的；（3）不考慮熱電離過(guò)程中出現(xiàn)概率極小的二次電離或者多重電離的情況。

2 數(shù)學(xué)模型

2.1 燃?xì)鉁囟裙?/p>

溫度是與壓力、彈丸行程、火藥已燃百分比有關(guān)的函數(shù)。這些數(shù)據(jù)在內(nèi)彈道方程組中都可以求得。

2.2 電子密度公式

通過(guò)燃?xì)鉁囟裙?，可以獲得火藥燃?xì)獾娜紲?。由于火藥燃?xì)庵猩傻入x子體的方式是熱電離，所以可以應(yīng)用薩哈方程計(jì)算燃?xì)庵须娮用芏龋?/p>

式（2）中，為熱力學(xué)溫度；E為離子的電離電位；g為原子基態(tài)的統(tǒng)計(jì)權(quán)重；g為離子基態(tài)的統(tǒng)計(jì)權(quán)重；m為電子質(zhì)量；n為電子密度；n為離子密度；n為原子密度。對(duì)于堿金屬來(lái)說(shuō)，的值約為1，其他氣體一般為2。

2.3 火藥燃燒時(shí)生成等離子體密度數(shù)學(xué)模型

火藥燃燒時(shí)生成等離子體密度方程組由3部分組成：內(nèi)彈道方程組、燃?xì)鉁囟裙胶碗娮用芏裙絒9-10]。

3 數(shù)值求解及結(jié)果分析

采用四階龍格庫(kù)塔法求解上述偏微分方程組，以MATLAB語(yǔ)言編寫(xiě)內(nèi)彈道計(jì)算程序。以某型坦克炮發(fā)射裝藥為例進(jìn)行數(shù)值模擬。

3.1 仿真結(jié)果

圖1為火炮內(nèi)彈道仿真結(jié)果與實(shí)驗(yàn)數(shù)據(jù)對(duì)比圖。由圖1可知數(shù)值仿真獲得的數(shù)據(jù)與實(shí)驗(yàn)結(jié)果基本吻合，圖像基本相近，驗(yàn)證了所建模型的正確性。圖2是等離子體密度與溫度曲線，由圖2可知在溫度達(dá)到2 000K以上時(shí)，電離種子碳酸鉀就會(huì)分解，生成鉀原子，鉀原子又發(fā)生熱電離，生成電子和離子。從圖2中可以看出隨著溫度的上升，電子密度增大，而且溫度越高，電子密度增加的速率越快，也就是熱電離的程度更加劇烈。

圖1 實(shí)驗(yàn)曲線與仿真曲線對(duì)比

圖2 等離子體密度—溫度曲線

3.2 裝藥量變化對(duì)生成等離子體密度的影響

在火炮設(shè)計(jì)的過(guò)程中，基于對(duì)膛壓和初速的要求，通常采用改變裝藥量的方法，來(lái)達(dá)到初速和膛壓的指標(biāo)。裝藥量的增加，實(shí)際上就是燃燒生成火藥燃?xì)饪偭康脑黾?，其?huì)導(dǎo)致膛內(nèi)壓力的增加，而壓力的增大會(huì)使彈丸的推力變大，從而使彈丸的初速增加。裝藥量的變化對(duì)膛壓和初速的影響見(jiàn)表1。

表 1 裝藥量的變化與膛壓和初速的關(guān)系

Tab.1 Relation between the changes of charge mass and muzzle velocity, as well as chamber pressure

從表1中可以看出，隨著裝藥量的變化，膛壓和炮口初速都會(huì)增大，但是膛壓的增速更快。根據(jù)燃速方程，壓力增加會(huì)導(dǎo)致火藥燃速增大，結(jié)果導(dǎo)致火藥更早地結(jié)束燃燒，而當(dāng)火藥燃燒結(jié)束以后，膛內(nèi)燃?xì)鉁囟燃眲∠陆?。圖3為不同裝藥量下，火藥燃燒百分比和時(shí)間的關(guān)系曲線。

圖3 不同裝藥量下火藥燃燒百分比——時(shí)間的曲線

圖4是在不同裝藥量情況下，等離子體生成密度的變化曲線。為了便于觀察，將等離子體密度取對(duì)數(shù)。從圖3~4中可以看出，隨著裝藥量的增加，等離子體的密度略有下降，等離子體存在的時(shí)間也逐漸減少。

圖4 不同裝藥量下等離子體密度——時(shí)間曲線

3.3 藥室容積變化對(duì)生成等離子體密度的影響

在裝藥量不變的情況下，藥室容積的變化既裝填密度的變化。藥室容積分別取132cm3、152cm3、172 cm3，進(jìn)行仿真試驗(yàn)。圖5為火藥燃燒百分比——時(shí)間的曲線圖。

圖5 不同藥室容積下火藥燃燒百分比——時(shí)間的曲線

從圖5中可以看出隨著藥室容積的增加火藥燃燒結(jié)束時(shí)間越來(lái)越晚，這表明火藥的燃速變得越來(lái)越慢。圖6是改變藥室容積，得到的等離子體密度——時(shí)間曲線。

圖6 不同藥室容積下等離子體密度——時(shí)間曲線

圖6表明隨著藥室容積的增大，火藥燃燒速度降低，使得火藥燃?xì)獗３指L(zhǎng)時(shí)間的高溫。相應(yīng)地，由于火藥燃?xì)鉁囟认陆禍p緩，等離子體密度的下降也減緩，既在藥室容積大的情況下，等離子體的密度也會(huì)大。

3.4 火藥力的變化對(duì)生成等離子體密度的影響

不同的火藥成分對(duì)應(yīng)不同的火藥力，分別取火藥力為900kJ/kg、950 kJ/kg、1 000 kJ/kg的火藥，進(jìn)行仿真試驗(yàn)。火藥力的增加會(huì)導(dǎo)致膛壓和初速的增加，而膛壓的增大會(huì)使火藥燃速更快。圖7~8是改變火藥力得到的溫度——時(shí)間曲線和等離子體密度——時(shí)間曲線。從圖7中可以看出隨著火藥力的增加，燃?xì)獾臏囟纫苍黾?。和裝藥量與藥室容積變化不同的是，火藥力的變化會(huì)導(dǎo)致燃?xì)獬跏紲囟鹊母淖?，這是因?yàn)榛鹚幜Φ母淖円簿褪腔鹚幊煞值母淖?，而燃?xì)獾谋瑴厥桥c火藥成分相關(guān)的。而火藥力的增加也就是增加燃?xì)饽芰康脑黾?，必然?huì)導(dǎo)致燃?xì)鉁囟鹊脑黾印８鶕?jù)等離子體密度方程，燃?xì)鉁囟鹊纳仙謺?huì)導(dǎo)致等離子體密度的增加，如圖8所示。

圖7 不同火藥力下溫度——時(shí)間曲線

圖8 不同火藥力下等離子體密度——時(shí)間曲線

4 結(jié)束語(yǔ)

根據(jù)火藥燃燒時(shí)的物理化學(xué)過(guò)程，利用內(nèi)彈道理論、燃燒學(xué)理論和等離子體物理理論建立了火藥燃燒時(shí)生成等離子體密度的數(shù)學(xué)模型。通過(guò)數(shù)值模擬，獲得的模擬計(jì)算結(jié)果和實(shí)驗(yàn)結(jié)果吻合較好，并計(jì)算了火藥燃燒時(shí)等離子體生成情況，最后通過(guò)分析裝藥量、藥室容積和火藥力等影響等離子體密度的因素，總結(jié)了等離子體密度變化的規(guī)律。本研究對(duì)提高火藥燃燒時(shí)生成的等離子體密度和進(jìn)一步研究等離子體在火炮方面的應(yīng)用具有重大意義。

[1] 徐學(xué)基,諸定昌.氣體放電物理[M].上海:復(fù)旦大學(xué)山版社,1996.

[2] 袁忠才,時(shí)家明,黃勇,馬柳. 低溫等離子體數(shù)值模擬方法的分析比較[J]. 核聚變與等離子體物理,2008(3):278-284.

[3] 黃慧,胡向南,馬鐵華,文治國(guó). 電熱化學(xué)炮膛壓測(cè)試系統(tǒng)的電磁兼容設(shè)計(jì)[J]. 傳感技術(shù)學(xué)報(bào),2014(9):1 226-1 231.

[4] 金涌.電熱等離子體對(duì)固體火藥的輻射點(diǎn)火及燃燒特性研究[D].南京:南京理工大學(xué),2014.

[5] 朱艷明.電熱化學(xué)炮內(nèi)彈道過(guò)程無(wú)網(wǎng)格法數(shù)值模擬[D].南京:南京理工大學(xué),2012.

[6] 胡嵐,張皋,王婧娜,高朗華.火藥燃燒氣體產(chǎn)物檢測(cè)方法研究[J].含能材料,2008(5):527-530.

[7] 孫艷馥,鮑雪. 高膛壓火炮內(nèi)彈道仿真研究[J]. 裝備制造技術(shù),2013(11):167-168.

[8] 李克婧.某艦炮彈道過(guò)程模擬仿真及優(yōu)化[D].南京:南京理工大學(xué),2012.

[9] 劉林,樊成軍,胡敬坤,撒彥成. 混合裝藥內(nèi)彈道仿真及多參數(shù)符合計(jì)算[J].火炮發(fā)射與控制學(xué)報(bào),2014(1):48-51.

[10] 李鋒,拜云山,馮曉偉. 一級(jí)氣體炮內(nèi)彈道方程修正及驗(yàn)證[J].彈道學(xué)報(bào),2016(1):14-18.

[11] 翟曉超,陳亞軍,姜毅.動(dòng)網(wǎng)格在仿真發(fā)射裝置內(nèi)彈道中的應(yīng)用[J].現(xiàn)代防御技術(shù),2006(2):24-28.

[12] 蘇菁.新型高能高燃速點(diǎn)傳火藥的研究[D].南京:南京理工大學(xué),2013.

[13] 劉東堯,趙振宇,余永剛,陸欣,張領(lǐng)科,周彥煌.阻燃劑包覆藥與常規(guī)火藥混合物的密閉爆發(fā)器燃燒試驗(yàn)及數(shù)值計(jì)算[J]. 火炸藥學(xué)報(bào),2011(4):83-86.

[14] 劉志濤,徐濱,南風(fēng)強(qiáng),廖昕.低溫感包覆火藥燃燒過(guò)程破孔規(guī)律研究[J].南京理工大學(xué)學(xué)報(bào),2011(5):709-713.

[15] 王靜.火藥燃燒催化劑的制備及其應(yīng)用研究[D].南京:南京理工大學(xué),2009.

[16] 胡嵐,張皋,王婧娜,高朗華.火藥燃燒氣體產(chǎn)物檢測(cè)方法研究[J].含能材料,2008(5):527-530.

[17] 謝玉樹(shù).等離子體與火藥相互作用的研究[D].南京:南京理工大學(xué),2005.

Numerical Simulation Study on the Law of Plasma Generated by Propellant Combustion

SONG Peng1，MAO Bao-quan2，ZHONG Meng-chun2，LI Xiao-gang2，LAN Tu2

(1. Huayin Weapon Test Center in China，Huayin，714200；2.Academy of Armored Forces Engineering，Beijing，100072)

To research the rules of plasma generated by propellant combustion, the physical and chemical process of the propellant combustion was analyzed. Mathematical model of plasma density was established, then the numerical simulation analysis of the gunpowder combustion generating plasma rules were carried out. The numerical simulation results were compared with the experimental results, and the correctness and feasibility of the model were determined. Finally, using the established mathematical model, the effect of different structural parameters of charge on the plasma was analyzed, and the plasma density change rule was obtained. The study laid the foundation for the next step to improve plasma concentration.

Propellant；Plasma；Combustion；Numerical simulation

1003-1480（2017）01-0042-04

TQ562

A

2016-10-17

宋鵬（1990-），男，在讀碩士研究生，主要從事武器系統(tǒng)設(shè)計(jì)研究。

武器裝備軍內(nèi)科研項(xiàng)目（項(xiàng)目號(hào)2014ZB03）。