爆炸力学高精度数值模拟研究进展

爆炸问题的高精度数值模拟在武器弹药设计及毁伤评估、工业重大爆炸灾害等国防和民用问题中具有重要的理论价值和应用前景.本文对近十几年来爆炸力学高精度计算的研究作了回顾与评述,结合作者近年来的研究工作,着重总结了高精度数值模拟保正性、高精度边界条件、自适应网格技术和多物质界面处理等方面的研究进展.推广并发展了WENO格式和RKDG格式,构造了守恒型高精度保正计算格式,在复杂爆轰波问题的数值模拟中很好地解决了出现负密度和负压强的问题,清晰地捕捉了爆轰波阵面结构和流动特征;提出了基于笛卡尔网格的高精度Inverse Lax-Wendroff复杂边界处理方法,该方法通过Inverse Lax-Wendroff获得一阶的法向偏导数值,所有其他的高阶法向偏导数值由五阶WENO类型的外插方法获得,实现了含复杂边界爆轰问题的高精度数值模拟,并具有较好的稳定性和鲁棒性;将后验误差估计方法应用于求解Euler方程中,采用网格的后验误差作为细分的判据,网格细分时在空间上利用高阶WENO插值,时间上采用Hermite插值.对于不需要细分的网格,本文采用点对点的直接赋值方法进行合并;发展了多介质界面处理的Local Level Set方法,对计算域内速度场进行修正,有效避免了速度场间断在求解界面位置时引起的误差;将GFM(ghost fluid method)和RGFM(real ghost fluid method)界面处理方法相结合,避免了GFM方法处理物质界面强间断时引入的虚假物理解,克服了RGFM方法处理弱间断时,重复求解Riemann问题引起的计算资源浪费.最后,对爆炸力学高精度计算的发展前沿进行了展望.