论文部分内容阅读
凝聚态炸药的爆炸流场广泛存在于国防和国民经济建设的许多领域,典型的爆炸流场中存在三种不同形式的运动界面:炸药/爆轰产物与空气的界面、空气与可动物体的界面、爆轰产物/空气与水的界面,本文采用数值方法对这三类运动界面问题进行研究。本文首先发展了一套动态混合网格的多组分流动数值模拟算法及软件,并以此作为研究工作的基础。该软件在动态混合网格上求解任意拉格朗日-欧拉形式的积分型控制方程,空间离散采用AUSM+-up和HLLC格式以及Venkatakrishnan限制器,时间离散采用Runge-Kutta格式,具有时空二阶精度,可以处理一维/二维/三维及球、柱、轴对称问题。计算方法不涉及特征矩阵的分解,适于计算各种形式的单流体/混合流体状态方程。在动网格方面,采用弹簧近似和网格重构相结合的方法模拟运动边界,采用虚拟网格通气技术解决运动物体从接触到分离过程中网格拓扑变化的难题。针对炸药/爆轰产物与空气的界面,采用流体混合型方法,提出了一种炸药/爆轰产物/空气的流体混合模型。该模型中炸药及爆轰产物采用JWL状态方程,空气采用完全气体模型,在固相(炸药)与气相(爆轰产物、空气)间采用等压假设,且体积可加;气相间满足等温假设及分压定理。该模型中混合流体的状态方程无需迭代求解,计算效率较高。将该模型与点火-生长反应模型结合,可以有效处理炸药/爆轰产物与空气间界面的复杂运动问题。数值模拟了空气中球形TNT装药的爆轰问题,计算得到了清晰的波系结构及流体界面,爆炸近区的峰值超压与实验符合较好。研究了驱动管中爆轰波的传播过程,得到了典型的有限直径药柱中的二维爆轰波结构,获得了管壁及管底的压力加载曲线,为驱动管设计提供了重要参考;分析了不同的驱动管排布方式以及延长段长度对激波管性能的影响。针对流体与可动物体间的动边界绕流问题,本文采用动网格跟踪运动边界,推导了ALE形式的流体控制方程,研究了动网格的几何守恒律,分析了不满足几何守恒律导致计算误差的原因,并提出了一种新的体积修正方法,可以消除网格运动引起的计算误差。将多体动力学求解方法与流场计算相耦合,解决了存在相互作用力的多个刚体的绕流的非定常求解方法。综合采用网格变形、网格重构、多体动力学耦合求解、虚拟网格通气等技术,仿真了内嵌式助推器的分离过程,发现了进气道整流罩打开后,进气道和燃烧室内存在衰减较慢的大幅压力振荡;助推器头部通过前体飞行器的喉部时,助推器所受的法向力剧烈振荡,给分离安全性带来一定影响。模拟了爆炸波驱动下卡车的运动,验证了本文方法对于含有运动物体的爆炸问题的模拟能力。针对水下爆炸中的水-气界面问题,本文采用ALE方法追踪流体界面,并将改进的虚拟流体方法应用于界面计算中。对气-气、气-水、爆轰产物-水等多个Riemann问题的模拟表明,相对于Euler型方法以及原始的ALE界面追踪方法,改进后的ALE方法具有更小的界面耗散以及更高的计算精度,界面压力无振荡。模拟了一维和二维水下爆炸问题,结果反映出了激波与水-气界面的相互作用,得到的气泡最大半径和首次脉动周期与实验值符合较好;二维水下爆炸中的气泡界面清晰,且不存在物理量的虚假振荡。最终,通过本文的工作,形成了模拟空中爆炸和水下爆炸的数值计算软件,可用于计算炸药爆轰、爆炸近区流场、爆炸波与刚性物体的相互作用,可以处理气-气、气-固、气-液等多种流体界面。软件具有较高的精度和较好的鲁棒性、通用性,并在工程实际问题中得到了应用。