韓 勇,龍新平,蔣治海,黃毅民,何 碧,吳 雄

（1.中國工程物理研究院化工材料研究所,四川綿陽 621900;

2.中國工程物理研究院,四川 綿陽 621900;

3.西安近代化學(xué)研究所,陜西 西安 710065）

1 引 言

炸藥爆轟性能參數(shù)的理論計算對炸藥配方設(shè)計十分重要,炸藥爆轟研究者對此開展了大量的工作。C.L.Mader[1]根據(jù)BKW爆轟產(chǎn)物狀態(tài)方程,編寫了BKW程序。其他研究者也開發(fā)了基于不同狀態(tài)方程的爆轟性能參數(shù)計算程序,如基于JCZ3狀態(tài)方程的Tiger[2]、CHEETAH[3]程序,基于WCA狀態(tài)方程的CHEQ[4]程序等。20世紀(jì)80年代,吳雄[5]應(yīng)用張光鑒[6]的相似理論,提出了一個以L-J勢函數(shù)為基礎(chǔ)的簡化維里模型的爆轟產(chǎn)物狀態(tài)方程,即VLW爆轟產(chǎn)物狀態(tài)方程,并編制了FORTRAN VLW程序。不同形式的狀態(tài)方程意味著對理論狀態(tài)方程不同方式的簡化,從而所采用勢參數(shù)也應(yīng)有所不同。JCZ3、WCA等狀態(tài)方程中所選用勢參數(shù)值均通過與實驗壓縮數(shù)據(jù)或分子動力學(xué)計算數(shù)據(jù)相匹配來調(diào)節(jié)獲得[7-9],而我國在VLW狀態(tài)方程方面尚未進行該方面的工作。在VLW程序中,部分爆轟產(chǎn)物勢參數(shù)引用文獻值,部分爆轟產(chǎn)物根據(jù)炸藥爆轟性能實驗值調(diào)整獲得[10],并未與該產(chǎn)物的靜壓或沖擊Hugoniot實驗數(shù)據(jù)對比。針對這一不足,本文中結(jié)合VLW狀態(tài)方程,初步開展不同勢參數(shù)值對計算炸藥爆轟產(chǎn)物組分之一的水的沖擊Hugoniot曲線的影響研究,并將計算結(jié)果與文獻實驗結(jié)果及采用BKW狀態(tài)方程的計算結(jié)果進行比較分析。

2 VLW狀態(tài)方程及沖擊Hugoniot計算原理

VLW狀態(tài)方程由維里狀態(tài)方程簡化而來,其假設(shè)高階維里系數(shù)與低階維里系數(shù)相似,高階維里系數(shù)可通過二階維里系數(shù)求得,進而將維里物態(tài)方程以簡便形式寫出

VLW狀態(tài)方程的第二階維里系數(shù)通過理論計算獲得

根據(jù)熱力學(xué)理論,若選取獨立變量溫度T、摩爾體積V,則熱力學(xué)函數(shù)p（T,V）、E（T,V）可從自由能F=F（T,V）導(dǎo)出

結(jié)合式（1）、（3）推導(dǎo)可得物質(zhì)自由能的表達式

式中 :q=B*/T*1/4,w=b0/V,φ0=[（G0-）0+（）0]/（R T）,其中（G0-）0為自由能,為焓。根據(jù)式（4）、（5）可得內(nèi)能E的表達式

物質(zhì)的沖擊Hugoniot關(guān)系為

式中:EH、pH、VH分別為物質(zhì)受沖擊后的內(nèi)能、壓力及體積,E0、p0、V0分別為物質(zhì)的初始內(nèi)能、壓力及體積。

由式（1）、（5）、（6）、（7）即可計算獲得物質(zhì)的沖擊Hugoniot曲線。

3 計算結(jié)果與討論

水的初始態(tài)參數(shù)引用文獻值[12]:T0=294.7 K,V0=18.014 cm3/mol,p0=100 kPa,E0=19.757 9 kJ/mol。ε、b0分別取不同參數(shù)值,計算獲得水的沖擊Hugoniot曲線如圖1所示。圖1（a）中b0=30.42,ε/K=80～ 400;圖 1（b）中 ε/K=120,b0=28.0 ～32.0。

圖1 水的沖擊Hugoniot曲線的計算值與文獻值比較Fig.1 Comparison between calculation data and experiment data about Hugoniot curve of water

由圖1可得,適當(dāng)調(diào)整勢參數(shù)的值,可改變計算所得水的Hugoniot曲線位置。當(dāng)勢參數(shù)ε/K=120、b0=30.42時,采用VLW狀態(tài)方程能夠與A.C.Mitchell等[13]的實驗結(jié)果和楊向東等[14]的理論計算結(jié)果符合較好。而由文獻[15]可得,采用VLW狀態(tài)方程計算炸藥爆轟性能參數(shù)時所采用的水的勢參數(shù)為:ε/K=180、b0=30.42,與我們計算所得的勢參數(shù)有所差異。然而,VLW程序仍然能夠成功計算炸藥的爆轟性能參數(shù),這一方面是由于所選用參數(shù)與優(yōu)選參數(shù)所得水的沖擊Hugoniot曲線相差不大,另一方面也與該程序的開發(fā)者在混合法則及部分爆轟產(chǎn)物勢參數(shù)的調(diào)整方面所作的努力分不開。

圖2 采用VLW和BKW狀態(tài)方程計算的水的沖擊Hugoniot曲線與文獻值的比較Fig.2 Comparison between calculation data by VLW and BKW EOS and experiment data about Hugoniot curve of water

分別采用VLW狀態(tài)方程及BKW狀態(tài)方程計算水沖擊Hugoniot曲線的結(jié)果與文獻值的比較如圖2所示。VLW狀態(tài)方程選取2套勢參數(shù)（ε/K,b0）,分別為（120,30.42）、（180,30.42）;BKW 狀態(tài)方程選取2種余容值k=250、420。其中,k=250為C.L.Mader[1]擬合水的沖擊Hugoniot曲線所選用值,k=420為根據(jù)分子結(jié)構(gòu)計算的幾何余容。由圖2可得,即使當(dāng)勢參數(shù)ε/K=180、b0=30.42時,在爆轟條件（壓力p一般小于50 GPa）下,VLW狀態(tài)方程仍然能夠比BKW狀態(tài)方程更準(zhǔn)確描述水的沖擊 Hugoniot關(guān)系,這也進一步證明了采用該套參數(shù)預(yù)測炸藥爆轟性能的有效性。根據(jù)2種狀態(tài)方程的形式分析可得,采用VLW狀態(tài)方程之所以能夠比BKW狀態(tài)方程更準(zhǔn)確地描述水的沖擊Hugoniot關(guān)系,是因為BKW狀態(tài)方程是基于維里方程的一級近似建立的,而VLW狀態(tài)方程的前3項是維里方程的理論式,第4項以后才是經(jīng)驗式,VLW狀態(tài)方程更趨近于理論狀態(tài)方程[16]。

4 結(jié) 論

（1）通過對勢函數(shù)參數(shù)的調(diào)整,采用VLW狀態(tài)方程可以較好地擬合炸藥爆轟產(chǎn)物組分 水的沖擊Hugoniot關(guān)系曲線;（2）與BKW狀態(tài)方程相比,VLW狀態(tài)方程能夠更準(zhǔn)確地描述水的沖擊壓縮狀態(tài);（3）通過對炸藥其他爆轟產(chǎn)物沖擊Hugoniot實驗數(shù)據(jù)的擬合,可建立一套基于基礎(chǔ)實驗、并適用于VLW狀態(tài)方程的勢函數(shù)參數(shù)。

[1] Mader C L.Numerical modeling of explosives and propellants[M].New York:CRC Press,1998.

[2] Cowperthwaite M,Zwisler W H.The JCZ equations of state for detonation products and their incorporation into the Tiger code[C]//Proceedings of the 6th International Symposium on Detonation.Coronado,California,1976:162-172.

[3] Howard W M,Fried L E,Souers P K.Kinetic modeling of non-ideal explosives with CHEETAH[R].DE200114716,2001.

[4] Ree F H.Statistical mechanical theory of chemically reacting multiphase mixtures:Application to detonation properties of PETN[J].Journal of Chemical Physics,1984,81（3）:1251-1263.

[5] WU Xiong.Detonation performance of condensed explosive computed with the VLW equation of state[C]//Proceedings of the 8th International Symposium on Detonation.Albuquerque,New Mexico,1985:796-804.

[6] 張光鑒.相似論[M].南京:江蘇科技出版社,1992.

[7] Hobbs M L,Baer M R,Mcgee B C.Extension of the JCZ product species database[C]//Proceedings of the 11th International Symposium on Detonation.Snowmass,Colorado,1998:958-968.

[8] Fried L E,Howard W M,Souers P C.Exp6:A new equation of state library for high pressure thermochemistry[C]//Proceedings of the 12th International Symposium on Detonation.San Diego,California,2002.

[9] Fried L E,Nir G,I-Feng W K,et al.Non-molecular phases of H2O and HF under detonation-like conditions[C]//Proceedings of the 13th International Symposium on Detonation.Virginia,Norfolk,2006.

[10] 龍新平,何碧,蔣小華,等.論VLW 狀態(tài)方程[J].高壓物理學(xué)報,2003,17（4）:247-254.

LONG Xin-ping,HE Bi,JIANG Xiao-hua,et al.Discussions on the VLW equation of state[J].Chinese Journal of High Pressure Physics,2003,17（4）:247-254.

[11] Hirschfelder J O,Curtiss C F,Bird R B.Molecular theory of gases and liquids[M].New York:Wiley,1964.[12] Ree F H.Molecular interaction of dense water at high temperature[J].Journal of Chemical Physics,1982,76（12）:6287-6302.

[13] Mitchell A C,Nellis W J.Equation of state and electrical conductivity of water and ammonia shocked to the 100 GPa（1 M bar）pressure range[J].Journal of Chemical Physics,1982,76（12）:6273-6281.

[14] 楊向東,胡棟,經(jīng)福謙.炸藥爆轟產(chǎn)物液氮液氦和水狀態(tài)方程研究[J].高壓物理學(xué)報,1999,13（2）:93-102.

YANG Xiang-dong,HU Dong,JING Fu-qian.Studies of EOS for detonation products:Liquid nitrogen,liquid helium and water[J].Chinese Journal of High Pressure Physics,1999,13（2）:93-102.

[15] 吳雄.新型爆轟產(chǎn)物物態(tài)方程[J].高壓物理學(xué)報,1991,5（2）:98-103.

WU Xiong.A new equation of state for detonation products[J].Chinese Journal of High Pressure Physics,1991,5（2）:98-103.

[16] 吳雄,龍新平,何碧,等.VLW 狀態(tài)方程的回顧與展望[J].高壓物理學(xué)報,1999,13（1）:55-58.

WU Xiong,LONG Xin-ping,HE Bi,et al.Review and look forward to the progress of VLW equation of state[J].Chinese Journal of High Pressure Physics,1999,13（1）:55-58.