多块结构网格上的Kershaw扩散格式.pdfVIP

计 算 物 理 V01．28．No．5 第 28卷 第 5期 CHINESE JOURNAL OFCOMPUTATIONAL PHYSICS Sep．，2011 2011年 9月 文章编号 ：1001—246X(2011)05-0641-08 多块结构网格上的Kershaw扩散格式 曾清红， 裴文兵， 成 娟， 勇 珩 (北京应用物理与计算数学研究所 ，北京 100094) 摘 要 ：Kershaw格式是在 四边形结构网格上求解扩散方程 的一种经典格式 ．基于对 Kershaw格式中 “流”的深 入理解 ，将其拓展到包含非结构点的多块结构网格 ，分别推导退化非结构点和强化非结构点情况下的Kershaw格 式，拓展的Kershaw格式满足流连续条件．三个数值算例的计算结果与精确解吻合得很好，表明将 Kershaw格式 拓展到多块结构网格的正确性和有效性 ． 关键词 ：惯性约束聚变 ；扩散方程 ；Kershaw格式；多块结构网格 ；非结构点 中图分类号 ：0241．8 文献标识码 ：A 0 引言 在惯性约束聚变 (ICF)的数值模拟 中，辐射能量方程的离散是辐射流体力学计算 的重要组成部分．一般 情况下，电子、离子热传导用非线性扩散 ，辐射输运使用多群扩散或多群输运来建模 ．扩散方程的精确求解对 于 ICF的数值模拟非常重要． 电子、离子或辐射 的扩散方程 ： 譬a： ．(、Kvr)， 、(1) 其中扩散系数 K是温度 的函数．1981年，Kershaw…使用变分原理和坐标变换 (物理空间变换到逻辑空 间)，在任意四边形结构网格上对方程 (1)建立差分格式，即Kershaw格式．在 Kershaw格式中，方程 (1)中的 扩散项用一个差分算子A来近似， · (，cT) AT=A ．肌，TK…， (2) 下标 K 表示对 网格 (K，L)及其周围的网格求和．Kershaw格式是求解扩散方程 的一种经典格式 ，被广泛 应用于各种扩散问题 ．一些大型程序 ，如 LARED—H程序 、美国罗彻斯特大学激光能量实验室 (LLE) 的DRACO程序 等 ，都采用 Kershaw格式求解非线性扩散方程． 计算表明，在大变形扭 曲网格上 ，Kershaw格式相比于其它格式具有较好 的精度和健壮性n。。]．但是随 着 ICF研究的不断深入和 ICF程序的不断发展 ，结构网格已经不能适应程序发展的需要 ，需要采用多块结构 网格来处理复杂的几何形状和介质分布．为充分利用多块结构 网格的优势 ，又保证原有程序的延续性 ，本文 将 Kershaw格式拓展到多块结构 网格 ，以满足 ICF大型程序可持续发展的需求．多块结构 网格的块内部是结 构网格，块与块之间可以非结构拼接 (图1)．多块网格 的非结构拼接将 出现非结构点，包括退化点和强化点． Kershaw格式拓展的主要工作就是如何在这两类非结构点处构造扩散格式． 1 结构网格上的Kershaw格式 首先简单回顾一下结构网格上的Kershaw格式．在结构 网格上 ，定义 网格和节点编号如 图2所示 ，网格 收稿 日期 ：2010—10—21；修 回日期 ：2011—02—18 基金项 目：国家 自然科学基金 10901021资助项 目 作者简介 ：曾清红 (1978一)，男，汉族 ，湖南长沙，副研究员，博士，主要从事辐射流体力学计算方法研究 E—mail：zeng — qinghong@ iapcm．ac．ca 通讯作者 第 5期 曾清红等 ：多块结构网格上的Kershaw扩散格式 2 Kershaw格式中的流 Kershaw格式是从纯数学角度推导出来的格式 ，物理意义上的 “流”没有清晰的定义．为 了推导 Kershaw 格式中流的准