三维爆轰波传播的LS方法研究
|
摘要
本文首先回顾了爆轰传播理论发展历史，介绍了几种基本的理论模型及与爆
轰波传播有关的其它内容。然后对等位线方法（Level Set，LS）作了一个简单
的介绍，并根据该方法的原理，建立三维爆轰波传播的LS方程，研究该LS方
程的数值求解方法，编制了三维爆轰波传播计算程序LS3D，提出了处理边界条
件和爆轰波相互作用的方法，进行爆轰波相互作用多点起爆实验设计，完成各种
起爆条件下爆轰波发展过程的数值模拟，并与实验结果作了比较。
本文提出了在三维直角坐标系中爆轰波传播LS方程离散化方法，建立了相
应的差分方程，研究了对不同形状炸药在多种起爆形式（点、线、面、散心、聚
心）下初值条件的设定，以及各种边界条件的处理方法。根据本文建立的计算编
码LS3D，可以确定不同时刻爆轰波阵面的位置，再通过计算编码LS2DDRAW
和LS3DDRAW画出某一截面的波阵面曲线图，或在整个计算空间内的三维波阵
面曲面图。
本文的实验工作设计了两种多点起爆装置，采用了转镜式高速分幅照相和扫
描照相技术及SVR相机高速摄影技术，观测了几种多点起爆装置中爆轰波在炸
药柱（DG01A、B炸药和JB9014）底端面的发展过程及爆轰波相互作用后波形
的演变情况。
本文最后根据所建立的三维爆轰波传播的LS方程的基本理论模型和该方程
数值解法，编写了计算爆轰波传播的三维计算编码LS3D，并使用该编码对不同
形状炸药在不同起爆方式（点起爆、线起爆、面起爆、发散爆轰、聚心爆轰等）
下理想爆轰波的传播过程进行了计算，并且对钝感炸药JB9014几种多点起爆实
验及圆柱形和圆弧片JB9014的平面起爆条件下爆轰波在炸药中的传播进行了数
值模拟。并与实验结果进行比较，计算结果与实验结果符合得还比较好。另外，
本文还给出了其它几个算例，主要是聚心爆轰、爆轰波在多介质炸药的平面波透
镜中的传播及爆轰波的相互作用等。
本文计算结果表明，计算编码LS3D可以对不同起爆方式下不同形状炸药中
三维爆轰波的传播进行较准确的计算，适于处理不同形状、多介质的炸药装置中
三维爆更波传播的匕方法研究
爆轰传播问题，是工程计算的一种有力工具。但应指出的是，这类计算结果的精
度在很大程度上取决于所采用的爆速曲率关系是否准确，以及所采用的爆轰波在
炸药与空气或其它介质边界处的夹角是否合适等。这些关键参数的测量是近代爆
轰传播实验研究的主要课题。
三维爆轰波传播的数值模拟，采用流体力学编码计算是十分困难且很不准确
的。目前LS方法是一种较为现实和可靠的方法。本文的研究工作在这个重要方
向上取得了较大进展，提出了该方法关键问题的一些初步的解决方法，明确了为
发展和完善该方法应进一步探讨的难点。
The development history of detonation wave propagation is reviewed first in this
thesis. Some basic theory model and other contents about detonation wave
propagation are introduced. Then the level set method is introduced. By means of this
method, the properties of three dimensional detonation waves are studied. The main
contents of this thesis includes: the LS equations of three dimensional detonation
wave propagation; the numerical solution method of LS equations; the experimental
design of multi-point initiation; the numerical simulation results of detonation wave
propagating on different initiation conditions.
In this paper, the numerical methods for the LS equation of three-dimensional
detonation wave are introduced. Detail are introduced as following: transfering the LS
equation to difference equation, setting the initial conditions for different shape
explosive at different initiation way (point, line, plane, divergent or convergent); the
dealing method of different edge conditions; how to draw the curve graph in two
dimensional space or surface graph in three dimensional space with the calculated
data.
Two kinds of initiation way are designed in this paper. By means of high-speed
frame camera, streak camera and SVR electrical camera, the propagation process and
the interaction of the detonation wave are recorded for several kinds of explosive
(DGO IA, B and JB90 14) at different imitation way.
Finally, based on the basic theory model of LS method applied to three
dimensional detonation wave propagation, we compile the calculation code LS3D. By
employing the code LS3D, the propagation processes of ideal detonation wave at
different shape explosives initiated by different way (point, line, plane) are calculated.
The detonation propagation problems of insensitive high explosive JB9014 with
different shape (cylinder, arc, pie) initiated from one plane are numerical simulated.
The calculated results agree with the experimental results at some degree. In addition,
the problems of convergent detonation and detonation wave propagating in planewave
lens with multi-medium and the interaction of detonation waves are calculated.
All the calculated results show that: the code LS3D can be used for calculating
the detonation wave propagation in different explosive initiated by different way. It
can deal with different shape, multi-medium explosive. The veracity of calculated
results is depended on the relation of detonation speed with the curvature of
detonation front and the edge angle of explosive with air or other medium.
Measurement for these key parameters is an important study aspect of detonation
wave propagation.
It is very difficult and not veracious to simulate the three dimensional detonation
wave propagation by means of fluid mechanics code. But the Level Set method is a
practicable technology. In this paper, some primary methods of resolution for some
key problems of LS method are given, and some propose are given for improving this
method.
轰波传播有关的其它内容。然后对等位线方法（Level Set，LS）作了一个简单
的介绍，并根据该方法的原理，建立三维爆轰波传播的LS方程，研究该LS方
程的数值求解方法，编制了三维爆轰波传播计算程序LS3D，提出了处理边界条
件和爆轰波相互作用的方法，进行爆轰波相互作用多点起爆实验设计，完成各种
起爆条件下爆轰波发展过程的数值模拟，并与实验结果作了比较。
本文提出了在三维直角坐标系中爆轰波传播LS方程离散化方法，建立了相
应的差分方程，研究了对不同形状炸药在多种起爆形式（点、线、面、散心、聚
心）下初值条件的设定，以及各种边界条件的处理方法。根据本文建立的计算编
码LS3D，可以确定不同时刻爆轰波阵面的位置，再通过计算编码LS2DDRAW
和LS3DDRAW画出某一截面的波阵面曲线图，或在整个计算空间内的三维波阵
面曲面图。
本文的实验工作设计了两种多点起爆装置，采用了转镜式高速分幅照相和扫
描照相技术及SVR相机高速摄影技术，观测了几种多点起爆装置中爆轰波在炸
药柱（DG01A、B炸药和JB9014）底端面的发展过程及爆轰波相互作用后波形
的演变情况。
本文最后根据所建立的三维爆轰波传播的LS方程的基本理论模型和该方程
数值解法，编写了计算爆轰波传播的三维计算编码LS3D，并使用该编码对不同
形状炸药在不同起爆方式（点起爆、线起爆、面起爆、发散爆轰、聚心爆轰等）
下理想爆轰波的传播过程进行了计算，并且对钝感炸药JB9014几种多点起爆实
验及圆柱形和圆弧片JB9014的平面起爆条件下爆轰波在炸药中的传播进行了数
值模拟。并与实验结果进行比较，计算结果与实验结果符合得还比较好。另外，
本文还给出了其它几个算例，主要是聚心爆轰、爆轰波在多介质炸药的平面波透
镜中的传播及爆轰波的相互作用等。
本文计算结果表明，计算编码LS3D可以对不同起爆方式下不同形状炸药中
三维爆轰波的传播进行较准确的计算，适于处理不同形状、多介质的炸药装置中
三维爆更波传播的匕方法研究
爆轰传播问题，是工程计算的一种有力工具。但应指出的是，这类计算结果的精
度在很大程度上取决于所采用的爆速曲率关系是否准确，以及所采用的爆轰波在
炸药与空气或其它介质边界处的夹角是否合适等。这些关键参数的测量是近代爆
轰传播实验研究的主要课题。
三维爆轰波传播的数值模拟，采用流体力学编码计算是十分困难且很不准确
的。目前LS方法是一种较为现实和可靠的方法。本文的研究工作在这个重要方
向上取得了较大进展，提出了该方法关键问题的一些初步的解决方法，明确了为
发展和完善该方法应进一步探讨的难点。
The development history of detonation wave propagation is reviewed first in this
thesis. Some basic theory model and other contents about detonation wave
propagation are introduced. Then the level set method is introduced. By means of this
method, the properties of three dimensional detonation waves are studied. The main
contents of this thesis includes: the LS equations of three dimensional detonation
wave propagation; the numerical solution method of LS equations; the experimental
design of multi-point initiation; the numerical simulation results of detonation wave
propagating on different initiation conditions.
In this paper, the numerical methods for the LS equation of three-dimensional
detonation wave are introduced. Detail are introduced as following: transfering the LS
equation to difference equation, setting the initial conditions for different shape
explosive at different initiation way (point, line, plane, divergent or convergent); the
dealing method of different edge conditions; how to draw the curve graph in two
dimensional space or surface graph in three dimensional space with the calculated
data.
Two kinds of initiation way are designed in this paper. By means of high-speed
frame camera, streak camera and SVR electrical camera, the propagation process and
the interaction of the detonation wave are recorded for several kinds of explosive
(DGO IA, B and JB90 14) at different imitation way.
Finally, based on the basic theory model of LS method applied to three
dimensional detonation wave propagation, we compile the calculation code LS3D. By
employing the code LS3D, the propagation processes of ideal detonation wave at
different shape explosives initiated by different way (point, line, plane) are calculated.
The detonation propagation problems of insensitive high explosive JB9014 with
different shape (cylinder, arc, pie) initiated from one plane are numerical simulated.
The calculated results agree with the experimental results at some degree. In addition,
the problems of convergent detonation and detonation wave propagating in planewave
lens with multi-medium and the interaction of detonation waves are calculated.
All the calculated results show that: the code LS3D can be used for calculating
the detonation wave propagation in different explosive initiated by different way. It
can deal with different shape, multi-medium explosive. The veracity of calculated
results is depended on the relation of detonation speed with the curvature of
detonation front and the edge angle of explosive with air or other medium.
Measurement for these key parameters is an important study aspect of detonation
wave propagation.
It is very difficult and not veracious to simulate the three dimensional detonation
wave propagation by means of fluid mechanics code. But the Level Set method is a
practicable technology. In this paper, some primary methods of resolution for some
key problems of LS method are given, and some propose are given for improving this
method.
引文
[1] Eyring H，Powell R E，Duffy G H et al．The Stability of Detonation[J]．Chem．Rev．， 1949，45：69-181
[2] Campbell A W et al．Proc．of the 6th Symp．(Int．)on Detonation，1976，642-652．
[3] Leiper G A et al．Proc．of the 9th Symp．(Int．)on Detonation，1989，197-208
[4] W Fickett and W C Davis．Detonation．Univ．of California Press，1979
[5] Jaimin Lee．Detonation Shock Dynamics of Composite Energetic Materials[Ph．D Thesis]．Dept．of Materials and Metallurgical Engineering，New Mexico Institute of Mining and Technology，1990
[6] Bdzil J B．[J]．Fluid Mech．，v．108，1981，195-226．
[7] Bdzil J B and Stewart D S．Phys．Fluids A，v．1，n．7，1989，1261-1267．
[8] Stewart D S and Bdzil J B．T&AM ReportNo．481，Univ．of Illinois，Urbana，1986
[9] Bdzil J B，Fickett W and Stewart D S．the 9th Syrup．(Int．)on Detonation，1989， 730-742．
[10] Stewart D S．Investigations on Detonation Shock Dynamics and Related Topics．Final Report for DOE-LANL-9-XG8-3931-P-1，1992．
[11] Stewart D S，Bdzil D S．The Shock Dynamics of Stable Multidimensional Detonation[J]． Combust．and Flame，1988，72：311-323．
[12] Stewart D S，Bdzil J B．Examples of Detonation Shock Dynamics for Detonation Wave Spread Application．Proc．of 9th Symp．(Int．)on Detonation，1989，773-783．
[13] 张星，孙锦山．计算发散爆轰波传播的DSD方法．爆炸与冲击，1997，17(2) ，97-103．
[14] 孙承纬，赵锋，高文．研究直径效应的爆轰冲击波动力学方法．爆炸与冲击，1996，16 (3)
[15] 孙承纬．二维爆轰波传播的DSD方法．第三次全国爆轰学术会议论文集，1992．
[16] Mulder W，Osher S，and Sethian J A．Computing interface motion in compressible gas dynamics[J]．Comput．Phys．，100，1992，209-228．
[17] Osher S，Sethian J A．Front propagation with curvature-dependent speed：Algorithm based on hamilton-Jacobi formulations[J]．Comput．Phys．，79(1) ，1988，12-49．
[18] Adalsteinsson D and Sethian J A．A Fast Level Set Method for Propagating Interfaces[J]． Comp．Phys．，Vol．118(2) ，269-277．
[19] Aslam T D，Bdzil J B and Stewart D S．Level set methods applied to modeling detonation shock dynamics[J]．Comput．Phys．，126，1996，390-409．
[20] Michel Debruyne and Francois Enot．The level set method applied to the relaxed Jouguet model：numerical results vs．experiments．Proc．of 11th Symp．(Iat．)on Detonation，1998，924-932．
[21] 张旭．曲线坐标系中钝感炸药非理想爆轰波传播的Level Set方法：[硕士学位论文]． 绵阳：中国工程物理研究院流体物理研究所，1999．
[22] 陈森华，张旭，柏劲松．贴体坐标系中多维非理想爆轰波传播计算的位标函数法．高 压物理学报，2000，14(4) ，280-284．
[23] Hill L G，Bdzil J B and Aslam T D．Front curvature rate stick measurements and detonation shock dynamics calibration for PBX9502 over a wide temperature range． Proc．of 11th Symp．(Int．)on Detonation，1998，1029-1037．
[2] Campbell A W et al．Proc．of the 6th Symp．(Int．)on Detonation，1976，642-652．
[3] Leiper G A et al．Proc．of the 9th Symp．(Int．)on Detonation，1989，197-208
[4] W Fickett and W C Davis．Detonation．Univ．of California Press，1979
[5] Jaimin Lee．Detonation Shock Dynamics of Composite Energetic Materials[Ph．D Thesis]．Dept．of Materials and Metallurgical Engineering，New Mexico Institute of Mining and Technology，1990
[6] Bdzil J B．[J]．Fluid Mech．，v．108，1981，195-226．
[7] Bdzil J B and Stewart D S．Phys．Fluids A，v．1，n．7，1989，1261-1267．
[8] Stewart D S and Bdzil J B．T&AM ReportNo．481，Univ．of Illinois，Urbana，1986
[9] Bdzil J B，Fickett W and Stewart D S．the 9th Syrup．(Int．)on Detonation，1989， 730-742．
[10] Stewart D S．Investigations on Detonation Shock Dynamics and Related Topics．Final Report for DOE-LANL-9-XG8-3931-P-1，1992．
[11] Stewart D S，Bdzil D S．The Shock Dynamics of Stable Multidimensional Detonation[J]． Combust．and Flame，1988，72：311-323．
[12] Stewart D S，Bdzil J B．Examples of Detonation Shock Dynamics for Detonation Wave Spread Application．Proc．of 9th Symp．(Int．)on Detonation，1989，773-783．
[13] 张星，孙锦山．计算发散爆轰波传播的DSD方法．爆炸与冲击，1997，17(2) ，97-103．
[14] 孙承纬，赵锋，高文．研究直径效应的爆轰冲击波动力学方法．爆炸与冲击，1996，16 (3)
[15] 孙承纬．二维爆轰波传播的DSD方法．第三次全国爆轰学术会议论文集，1992．
[16] Mulder W，Osher S，and Sethian J A．Computing interface motion in compressible gas dynamics[J]．Comput．Phys．，100，1992，209-228．
[17] Osher S，Sethian J A．Front propagation with curvature-dependent speed：Algorithm based on hamilton-Jacobi formulations[J]．Comput．Phys．，79(1) ，1988，12-49．
[18] Adalsteinsson D and Sethian J A．A Fast Level Set Method for Propagating Interfaces[J]． Comp．Phys．，Vol．118(2) ，269-277．
[19] Aslam T D，Bdzil J B and Stewart D S．Level set methods applied to modeling detonation shock dynamics[J]．Comput．Phys．，126，1996，390-409．
[20] Michel Debruyne and Francois Enot．The level set method applied to the relaxed Jouguet model：numerical results vs．experiments．Proc．of 11th Symp．(Iat．)on Detonation，1998，924-932．
[21] 张旭．曲线坐标系中钝感炸药非理想爆轰波传播的Level Set方法：[硕士学位论文]． 绵阳：中国工程物理研究院流体物理研究所，1999．
[22] 陈森华，张旭，柏劲松．贴体坐标系中多维非理想爆轰波传播计算的位标函数法．高 压物理学报，2000，14(4) ，280-284．
[23] Hill L G，Bdzil J B and Aslam T D．Front curvature rate stick measurements and detonation shock dynamics calibration for PBX9502 over a wide temperature range． Proc．of 11th Symp．(Int．)on Detonation，1998，1029-1037．
© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号 地址：北京市海淀区学院路29号 邮编：100083 电话：办公室：(+86 10)66554848；文献借阅、咨询服务、科技查新：66554700 |