Grüneisen状态方程的数值解及其应用

计算物理 ›› 2007, Vol. 24 ›› Issue (3): 361-366.

Grüneisen状态方程的数值解及其应用

王庭辉¹, 宋顺成^1,2

1. 西南交通大学应用力学与工程系, 四川成都 610031;
2. 北京理工大学爆炸科学与技术国家重点实验室, 北京 100081

收稿日期:2006-02-23 修回日期:2006-06-29 出版日期:2007-05-25 发布日期:2007-05-25
作者简介:王庭辉(1982-),男,河南,博士生,从事高速冲击动力学方面的研究.

Numerical Study of Grüneisen State Equation and Applications

WANG Tinghui¹, SONG Shuncheng^1,2

1. Department of Applied Mechanics and Engineering, Southwest Jiaotong University, Chengdu 610031, China;
2. State Key Lab of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China

Received:2006-02-23 Revised:2006-06-29 Online:2007-05-25 Published:2007-05-25

摘要/Abstract

摘要： 从Grüneisen状态方程、冲击绝热关系出发,给出求解相关微分方程组的Runge-Kutta法格式.研究钨材料的内部压力、Grüneisen系数、冷压、冷能与比体积之间的数值关系,得出Grüneisen系数、冷压、冷能及常温下压力随比体积的变化曲线.计算给出的压力~比体积关系同实验数据符合较好.在此基础上给出两个压缩状态(比体积分别为V₁,V₂)下压力随温度变化的曲线,说明温度的变化对压力的影响小于比体积变化的影响.

关键词: 高压状态方程, Runge-Kutta法, P~V关系

Abstract: A computational scheme of Runge-Kutta method is presented for Grüneisen state equation with adiabatic impact relation.Thermodynamic state variables,including pressure,Grüneisen coefficient,cold pressure,cold energy and pressure,as a function of specific volume are calculated.Numerical states parameters of tungsten at high-pressure are given.The curves of γ~V,E_K~V,P_K~V,P~V of tungsten are obtained and the P~V result agrees well with the experiment.Numerical solutions of f(P,T,V₁)=0 and f(P,T,V₂)=0 in two compressive states show that the effect of temperature on pressure is less than that of specific volume.

Key words: EOS with high-pressure, Runge-Kutta method, relation of P-V

中图分类号:

O521.22

王庭辉, 宋顺成. Grüneisen状态方程的数值解及其应用[J]. 计算物理, 2007, 24(3): 361-366.

WANG Tinghui, SONG Shuncheng. Numerical Study of Grüneisen State Equation and Applications[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2007, 24(3): 361-366.

[1]	王凯歌, 牛憨笨, 周燕. 电子束在微束斑X射线源中运动轨迹的计算[J]. 计算物理, 2003, 20(1): 21-24.
[2]	肖家鑫, 林家逖, 郑熙铭. 修正Burridge-Knopoff模型的地震前兆场数值模拟[J]. 计算物理, 1997, 14(S1): 600-602.