轴对称坐标系下含曲率的水平集方程在非结构网格上的数值方法 Numerical scheme for level set equations containing curvature on axisymmetric unstructured meshes.pdfVIP

第32卷第2期 爆炸与冲击 VD】．32。No．2 ZOIZ年3月 EXPI．f)SI()NANDSHOCKWAVES Mal。，2012 文章编号：1001—1455(2012)02—0150—07 轴对称坐标系下含曲率的水平集方程 在非结构网格上的数值方法’ 程俊霞．胡晓棉 北京应用物理与计算数学研究所。北京100088) 摘要：为了在三角形和四边形网格上采用水平集方法模拟轴对称爆轰波阵面与曲率相关的运动．假设爆 轰渡阵面的法向速度是曲率的线性函数。通过坐标变换得到了轴对称坐标系下的水平集方程。水平集方程的 曲率无关项采用正格式离散，曲率项采用伽辽金等参有限元方法空间离散。时间离散采用半隐格式。算例表 明，在轴对称三角形网格和四边形网格上，台曲率的水平集方程的离散格式为强一阶精度，给出了三角形和 四边形混合网格上不光滑界面以曲率收缩的运动过程，收缩过程未出现不稳定现象。多个爆轰波阵面相互作 用的算例说明本文的格式可有效地模拟曲率相关的爆轰波的相互作用问题 关键词：爆炸力学；水平集方程；伽辽金等参有限元方法；爆轰波阵面；轴对称 中圈分类号：0382．4；0241．82国标学科代码：13035 文献标志码：A 与敏感炸药相比，钝感炸药的曲率效应较明显，因此在工程应用中必须考虑钝感炸药的曲率效应。 DSD方法…是一种考虑爆轰波曲率效应的工程爆轰流体力学数值方法。本文中基于DSD方法的基本 框架，致力于实现轴对称几何中曲率相关的爆轰波阵面的数值模拟方法。由于水平集方法容易处理界 面拓扑结构的变化．因此选择水平集方法模拟爆轰波阵面的运动。为了满足拉格朗日程序模拟爆轰波 阵面的需要，本文中在轴对称三角形和四边形网格上离散水平集方程。 到目前为止，大多数研究者心43在笛卡儿网格上采用各种精度的迎风型ENO差分格式离散水平集 方程。非笛卡儿网格上的水平集方程的离散方法并不多见，含曲率的高阶水平集方程在非笛卡儿网格 上的离散较困难。一种思路是把水平集方程转化为守恒率，从而可采用守恒率的方法。对于不可压流， E．Marchandise等¨1把与曲率无关的水平集演化方程写成守恒形式，采用间断有限元方法求解。F． 程，二阶导数项在笛卡儿网格上采用中心差分。针对含有高阶导数的方程，13．Cockburn等]提出了 LDG方法。把高阶方程写成一阶守恒型微分方程组，采用RKDG方法求解。此方法适用于对流占优的 扩散问题，如果扩散占优，时间步长将很受限制，效率降低。C．Q．Hu等旧1采用LDG方程求解了含曲 方法离散不含曲率的水平集方程，而含曲率的部分采用有限元方法二阶显式龙格一库塔法求解。与间断 有限元方法相比，此方法具有计算量小的优点，但显式求解含曲率的水平集方程要限制时问步长。K． Deckelnick等u“提出了含曲率的水平集方程的有限元离散格式，从格式的稳定性出发。提出了半隐离 散格式，求解了纽曼边界条件下三角形网格上的含曲率的水平集方程。水平集函数在演化过程中不能 保持距离函数的性质，梯度将变得过大或过小，因此需要重新初始化。F．Bornemann等¨L提出了水平 集函数在三角形网格上的重新初始化方法，适用于满足惠更斯原理的爆轰波传播问题。本文中拟基于 文献[9一11]．在三角形网格和四边形网格上求解轴对称几何中的含曲率的水平集方程。 -收稿日期：2010-07．09；修回日期：201卜03—25 基金项目：国家自然科学基金项目10931004) 作者简介：程俊霞(1977一)，女，博士研究生，副研究员。 万方数据 第2期 程俊霞等：轴对称坐标系下含曲率的水平集方程在非结构网格上的数值方法 151 1 轴对称几何中的含曲率的水平集方程 本文中求解的水平集方程在笛卡儿坐标系下可写为 “I+DJlVMI一口七lVMf=0 (1) 方程(1)高度非线性，是可能退化的对流扩散方程。