曲率相关的爆轰波在非结构四边形网格上的数值方法

计算物理 ›› 2011, Vol. 28 ›› Issue (2): 199-206.

曲率相关的爆轰波在非结构四边形网格上的数值方法

程俊霞

北京应用物理与计算数学研究所, 北京 100088

收稿日期:2010-04-01 修回日期:2010-11-09 出版日期:2011-03-25 发布日期:2011-03-25
作者简介:程俊霞(1977-),女,河南,副研究员,博士生,从事计算流体力学工作,北京8009信箱16分箱100088.
基金资助:
中物院科技发展基金(2007A09006);实验室基金(9140C690101070C69,9140C6901030803);国家重大基础研究基金(2005CB32170);国家自然科学基金(10901022)资助项目

Numerical Schemes for Detonation Front on Curvature on Unstructured Quadrilateral Meshes

CHENG Junxia

Institute of Applied Physics and Computational Mathernatics, Beijing 100088, China

Received:2010-04-01 Revised:2010-11-09 Online:2011-03-25 Published:2011-03-25

摘要/Abstract

摘要： 假设爆轰波阵面的法向速度是曲率的线性函数,在非结构四边形网格上采用水平集方法模拟爆轰波阵面的运动过程.水平集方程的曲率无关项采用正格式离散,曲率项采用伽辽金等参有限元方法空间离散,时间离散采用半隐格式.在笛卡儿网格和随机网格上,含曲率的水平集方程的离散格式为强一阶精度,重新初始化方程的离散格式精度为近似一阶精度.曲率收缩的不光滑界面和多个爆轰波阵面相互作用的算例说明格式可有效地模拟爆轰波与曲率相关的运动.

关键词: 爆轰波阵面, 水平集方程, 非结构四边形网格, 伽辽金等参有限元方法

Abstract: Assuming that detonation normal velocity is a linear function of curvature,we study propagation of detonation front on unstructured quadrilaterals with level set method.Convection term in leve set equation is solved by positive scheme,and curvature term is solved by Galerkin isoparametric finite element method and semi-implicit time stepping.On Cartesian meshes and random meshes,the scheme of level set equations containing curvature is more than first order accuracy in L₂ and L_∝ norms.Examples of nonsmooth level sets shortening and three detonation fronts colliding show that the schemes can be used to simulate propagation of detonation front on curvature.

Key words: detonation front, level set equations, unstructured quadrilateral meshes, Galerkin isoparametric finite element method

中图分类号:

程俊霞. 曲率相关的爆轰波在非结构四边形网格上的数值方法[J]. 计算物理, 2011, 28(2): 199-206.

CHENG Junxia. Numerical Schemes for Detonation Front on Curvature on Unstructured Quadrilateral Meshes[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2011, 28(2): 199-206.

[1]	赵丰祥, 潘亮, 王双虎. 基于非结构四边形网格的WENO格式[J]. 计算物理, 2018, 35(5): 525-534.
[2]	聂存云, 舒适, 盛志强. 非结构四边形网格下的一类保对称有限体元格式[J]. 计算物理, 2009, 26(2): 175-183.
[3]	王瑞利, 姚彦忠, 林忠, 葛全文. 一种新的非结构四边形网格生成方法及软件[J]. 计算物理, 2007, 24(1): 13-18.