三维爆轰波传播的LS方法研究

【摘要】： 本文首先回顾了爆轰传播理论发展历史，介绍了几种基本的理论模型及与爆 轰波传播有关的其它内容。然后对等位线方法（Level Set，LS）作了一个简单 的介绍，并根据该方法的原理，建立三维爆轰波传播的LS方程，研究该LS方 程的数值求解方法，编制了三维爆轰波传播计算程序LS3D，提出了处理边界条 件和爆轰波相互作用的方法，进行爆轰波相互作用多点起爆实验设计，完成各种 起爆条件下爆轰波发展过程的数值模拟，并与实验结果作了比较。 本文提出了在三维直角坐标系中爆轰波传播LS方程离散化方法，建立了相 应的差分方程，研究了对不同形状炸药在多种起爆形式（点、线、面、散心、聚 心）下初值条件的设定，以及各种边界条件的处理方法。根据本文建立的计算编 码LS3D，可以确定不同时刻爆轰波阵面的位置，再通过计算编码LS2DDRAW 和LS3DDRAW画出某一截面的波阵面曲线图，或在整个计算空间内的三维波阵 面曲面图。 本文的实验工作设计了两种多点起爆装置，采用了转镜式高速分幅照相和扫 描照相技术及SVR相机高速摄影技术，观测了几种多点起爆装置中爆轰波在炸 药柱（DG01A、B炸药和JB9014）底端面的发展过程及爆轰波相互作用后波形 的演变情况。 本文最后根据所建立的三维爆轰波传播的LS方程的基本理论模型和该方程 数值解法，编写了计算爆轰波传播的三维计算编码LS3D，并使用该编码对不同 形状炸药在不同起爆方式（点起爆、线起爆、面起爆、发散爆轰、聚心爆轰等） 下理想爆轰波的传播过程进行了计算，并且对钝感炸药JB9014几种多点起爆实 验及圆柱形和圆弧片JB9014的平面起爆条件下爆轰波在炸药中的传播进行了数 值模拟。并与实验结果进行比较，计算结果与实验结果符合得还比较好。另外， 本文还给出了其它几个算例，主要是聚心爆轰、爆轰波在多介质炸药的平面波透 镜中的传播及爆轰波的相互作用等。 本文计算结果表明，计算编码LS3D可以对不同起爆方式下不同形状炸药中 三维爆轰波的传播进行较准确的计算，适于处理不同形状、多介质的炸药装置中 三维爆更波传播的匕方法研究 爆轰传播问题，是工程计算的一种有力工具。但应指出的是，这类计算结果的精 度在很大程度上取决于所采用的爆速曲率关系是否准确，以及所采用的爆轰波在 炸药与空气或其它介质边界处的夹角是否合适等。这些关键参数的测量是近代爆 轰传播实验研究的主要课题。 三维爆轰波传播的数值模拟，采用流体力学编码计算是十分困难且很不准确 的。目前LS方法是一种较为现实和可靠的方法。本文的研究工作在这个重要方 向上取得了较大进展，提出了该方法关键问题的一些初步的解决方法，明确了为 发展和完善该方法应进一步探讨的难点。