三维多面体网格上扩散方程保正格式.pdfVIP

2015年 8月 计 算 数 学 第 37卷第 3期 Aug．，2015 MATHEMATICA NUMERICA SINICA Vo1．37，No．3 三维多面体网格上扩散方程的保正格式冰1) 王 帅 杭旭登 袁光伟 (IAPCM，北京应用物理与计算数学研究所计算物理实验室，北京100094) 摘 要 针对三维任意 (星形)多面体网格，本文构造了扩散方程的一种单元中心型非线性有限体积 格式，证明了该格式具有保正性．在该格式设计中，除引入网格中心量外，还引入网格节点量和网 格面中心量作为中间未知量，它们将用网格中心未知量线性组合表示，使得格式仅有网格中心未 知量作为基本未知量．在节点量计算中，利用网格面上的调和平均点，设计了一种适用于三维多面 体网格的局部显式加权方法．该格式适用于求解非平面的网格表面和间断扩散系数的问题．数值 例子验证了它对光滑解具有二阶精度和保正性． 关键词：多面体网格；保正格式；调和平均点；扩散方程 MR (2000)主题分类：65D17 1．引 言 辐射流体力学在惯性约束聚变 (ICF)、磁约束聚变和天体物理等诸多领域中都有着广泛 的应用．例如在惯性约束聚变 (ICF)中，腔壁由于受激而发射 x光，产生的光辐射在腔内传播， 到达聚变材料区后引起核聚变．其中腔内光辐射和物质能量的迁移现象由辐射流体力学方程 组中的能量守恒方程描述． 目前受实验条件和成本的制约，精密化数值模拟成为上述应用领域中主要的研究手段．由 于工程中对极高密度比的物质界面分辨精度的要求非常苛刻，所以，对于可压缩多介质辐射流 体力学问题，通常采用 Lagrange方法进行计算．随着流体的运动，计算网格会不断变形，能量 扩散方程的求解是在扭曲网格上进行的．因此，研究和构造扭曲网格上能量扩散方程的高精度 离散方法，对于提高Lagrange辐射流体力学数值模拟置信度至关重要．为了保证扩散计算具 有所需的置信度，不仅要求扩散格式具有所需要的精度和稳定性，还要求格式应尽可能多地保 持物理模型的基本性质，比如守恒性和保正性等．大量实际计算实践表明，在扭曲网格上，由 不满足保正性的扩散格式得到的数值解，会 出现热流从低温区流 向高温区的现象，这违背热力 学第二定律；在能量较低的区域，温度计算会出负；在温度变化剧烈的区域，会产生非物理振 荡 1【0j．对于能量扩散方程而言，方程的扩散系数为温度的非线性函数，它仅在温度非负时才 有定义，负的数值解将导致下一个非线性迭代步的扩散系数没有定义，程序运行被迫中止．实 际计算中常常采用”遇负置零”等方式进行修正，以保证计算顺利进行下去，但这样会破坏能 量的守恒性，难以保证长时间计算结果的置信度；并且我们的计算实践表明，当数值解出现非 物理振荡时，由自适应时间步长选取准则所得到的时间步长不断缩小，从而导致模拟效率的降 低． 对二维扭曲网格上扩散方程的求解，国内外已经进行了多年的研究，得到了很多各具特点 的有限体积格式．例如九点格式 1【】j支撑算子格式 [。]和多点流逼近格式 [。]等，这些格式都是 2014年 12月 13日收到． )基金项 目： 国家 自然科学基金 1140103411271054)资助项 目 计 算 数 学 2015芷 守恒的且具有一定的精度，但都不是保正格式，除非对网格和扩散系数添加较严格的限制 [46]． 在 [5]中给出了四边形网格上有限体积方法的保正性条件，其中指出当格式是守恒的并且对线 性问题精确时，不可能构造一个无条件保正且相容的九点格式．文 [7]针对三角形网格上的抛 物方程提出了一种非线．j生有限体积格式，它只在时间步长足够小时才是保正的．文 8『1进一步 发展和分析了 [7]中提出的保正格式，该格式的网格中心点的位置与扩散系数和网格形状有 关，因此它不能直接用于辐射流体力学计算．在辐射流体力学计算中，网格中心位置是由流体 力学计