论文部分内容阅读
水下爆炸的计算研究在军事领域、工程领域有广泛的应用前景。在军事领域,水雷和深水炸弹的研制需要对相应的水下爆炸波传播和超压峰值进行计算研究。在工程领域,港口航道的建设需要对水下爆炸进行计算研究。水下爆炸的计算研究可以为实验研究提供参考,特别是在实验条件不具备的情况下,更是可以为水下兵器的研制和工程爆破的设计提供理论指导并具备重要的参考意义。而对水下爆炸计算方法的研究和改进也具有重要的学术价值。目前,关于水下爆炸的计算研究,方法有有限元方法、SPH方法(Smoothed Particle Hydrodynamics)、边界元方法、GFM(Ghost Fluid Method)以及一些欧拉-拉格朗日耦合类的方法。其中,GFM又包括OGFM(Original GFM),NGFM(New Version of GFM),MGFM (Modified GFM),RGFM(Real GFM)和SMGFM(Solid MGFM)等等。目前,国内用GFM对水下爆炸进行的计算研究仍然很少。本文选用GFM对水下爆炸进行计算研究。本文根据不同的问题和计算条件选用不同的GFM子方法进行数值模拟。本文主要对水下爆炸的爆炸波传播进行数值模拟、对水下爆轰流场的超压峰值进行数值计算以及针对水下爆炸的数值计算改进GFM中的计算步骤。本文的主要研究内容和创新性工作如下:1)对柱形装药水下爆炸问题利用Level set方法完成了精确的三维数值模拟,提高了计算精度。柱形装药在水下爆炸中较为常见。对于三维圆柱问题的数值模拟一般是在计算域可以切割成二维轴对称区域的情况下,将三维问题做为二维轴对称问题处理。但是一些较复杂的计算域并不能切割成二维轴对称问题处理,这时就需要求解三维情况下的Level set函数初值。以往对初始时刻Level set函数的定义使用的是重新初始化技术,该技术会造成界面位置的误差。本文通过计算域分块定义的方法,对三维柱形装药水下爆炸数值模拟中的Level set函数初值进行精确求解,避免了使用重新初始化技术,从而减小了数值模拟结果的误差。2)本文将插值法耦合进ARPS(Approximate Riemann Problem Solver),从而减小计算误差。用ARPS预估界面状态是GFM中的一个重要步骤。原有的ARPS在用于多维问题时,使用最小夹角法则来选择界面两侧的网格节点,而最小夹角法则会造成较大的数值误差。本文用插值法替代最小夹角法则,减小了计算误差,并将耦合方法的计算结果同原有ARPS的计算结果进行对比分析。对比表明,耦合插值法的ARPS比原有的ARPS计算精度更高。3)本文用两种不同的方法数值模拟复杂计算域的爆炸波传播,并通过对比研究分析和总结了各自的优缺点。一是在贴体坐标网格下使用任意坐标系方法。二是在规则笛卡尔网格使用NGFM。研究结果表明,在规则笛卡尔网格中使用NGFM可以得到较稳定的计算结果;但是无法刚好得到壁面上网格节点的压力值,在实际计算中只能用紧邻壁面的网格节点的压力值做为壁面压力值的参考,不过计算结果表明在加密网格的情况下,壁面内侧网格节点的压力值和壁面压力值非常接近。相比较而言,在贴体坐标网格中使用任意坐标系方法,可以得到壁面网格节点的压力值;但是利用Tait方程作为水的状态方程时,流场常常会会出现较大范围的数值振荡,另一方面,网格导数的引入会影响计算结果的精度。4)本文推导了一个基于Euler-弹性动力学基本方程的流固ARPS,从而扩大了原有SMGFM的适用范围。用流固ARPS预估流固界面状态,是SMGFM中一个重要步骤。原有的流固ARPS的缺点在于是用Euler方程和Naviers方程推导的,而Naviers方程对固体的数值计算而言应用起来效果并不太好。本文用弹性动力学基本方程描述固体,推导一个基于Euler-弹性动力学基本方程的流固ARPS,并将计算结果同文献中刚壁计算结果、可压缩固体计算结果进行了对比分析。数值计算的结果表明本文推导的ARPS是正确的。同时,数值计算的结果也表明线弹性假设得到的流体一侧爆炸波传播结果和可压缩固体假设得到的结果相近,壁面的压力动态响应也较好地吻合。