轴对称坐标系下含曲率的水平集方程在非结构网格上的.pdf

第 卷 第 期 爆 炸 与 冲 击 , 32 2 Vol.32 No.2 年 月 , 2012 3 EXPLOSIONANDSHOCKWAVES Mar. 2012 文章编号: ( ) 1001-1455201202-0150-07 轴对称坐标系下含曲率的水平集方程 在非结构网格上的数值方法* , 程俊霞 胡晓棉 ( , ) 北京应用物理与计算数学研究所 北京 100088 : , 摘要 为了在三角形和四边形网格上采用水平集方法模拟轴对称爆轰波阵面与曲率相关的运动 假设爆 , 。 轰波阵面的法向速度是曲率的线性函数 通过坐标变换得到了轴对称坐标系下的水平集方程 水平集方程的 , , 。 曲率无关项采用正格式离散 曲率项采用伽辽金等参有限元方法空间离散 时间离散采用半隐格式 算例表 , , 。 明 在轴对称三角形网格和四边形网格上 含曲率的水平集方程的离散格式为强一阶精度 给出了三角形和 , 。 四边形混合网格上不光滑界面以曲率收缩的运动过程 收缩过程未出现不稳定现象 多个爆轰波阵面相互作 用的算例说明本文的格式可有效地模拟曲率相关的爆轰波的相互作用问题 : ; ; ; ; 关键词 爆炸力学 水平集方程 伽辽金等参有限元方法 爆轰波阵面 轴对称 中图分类号: ; 国标学科代码: 文献标志码: O382.4O241.82 13035 A , , 。 与敏感炸药相比 钝感炸药的曲率效应较明显 因此在工程应用中必须考虑钝感炸药的曲率效应 [] 1 。 DSD方法 是一种考虑爆轰波曲率效应的工程爆轰流体力学数值方法 本文中基于DSD方法的基本 , 。 框架 致力于实现轴对称几何中曲率相关的爆轰波阵面的数值模拟方法 由于水平集方法容易处理界 , 。 面拓扑结构的变化 因此选择水平集方法模拟爆轰波阵面的运动 为了满足拉格朗日程序模拟爆轰波 , 。 阵面的需要 本文中