二维柱坐标系辐射输运方程保球对称的离散纵标格式

来源：计算数学杂志2015年第3期北京时间：

作者：宋鹏, 成娟

单位：1. 北京应用物理与计算数学研究所, 计算物理重点实验室, 北京 100094;\r\n2. 北京应用物理与计算数学研究所, 北京 100094

摘要：离散纵标格式是计算辐射输运方程的常用格式之一. 但是, 传统的离散纵标格式求解二维柱坐标系辐射输运方程模拟一维球对称问题时, 会出现明显的非对称现象, 球对称性被破坏. 针对该问题, 本文分析了传统离散纵标格式不能够保持球对称性的原因, 提出了空间基于柱坐标系、方向基于球坐标系的辐射输运方程, 并对该方程设计了新的离散纵标格式, 从理论上证明了当空间网格取球对称剖分时该离散格式能够保持一维球对称性的充分必要条件. 通过对真空球区域辐射输运、与物质耦合辐射输运等球对称算例的数值模拟, 验证了该格式的保球对称性, 球对称误差能够达到机器精度. 非对称辐射驱动模型以及非对称网格剖分条件下的数值模拟等算例也取得了较好的结果.

关键词：辐射输运,离散纵标格式,柱坐标系,球对称

基金资助：<p>国家自然科学基金(11401033, 11471049, 91330205)、国家高技术研究发展计划(2012AA01A303)、国家基础研究项目(2011CB309702)资助项目</p>