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

2015年 8月 计 算 数 学 第 37卷第 3期 Aug．，2015 MATHEMATICA NUMERICA SINICA Vo1．37，No．3 二维柱坐标系辐射输运方程保球对称的 离散纵标格式木1) 宋 鹏 (北京应用物理与计算数学研究所，计算物理重点实验室，北京100094) 成 娟 (北京应用物理与计算数学研究所，北京100094) 摘 要 离散纵标格式是计算辐射输运方程的常用格式之一．但是，传统的离散纵标格式求解二 维柱坐标系辐射输运方程模拟一维球对称 问题时，会 出现明显的非对称现象，球对称性被破 坏．针对该问题，本文分析了传统离散纵标格式不能够保持球对称性的原因，提出了空间基于 柱坐标系、方 向基于球坐标系的辐射输运方程，并对该方程设计 了新的离散纵标格式，从理 论上证明了当空间网格取球对称剖分时该离散格式能够保持一维球对称性的充分必要条件． 通过对真空球区域辐射输运、与物质耦合辐射输运等球对称算例的数值模拟，验证了该格式 的保球对称性，球对称误差能够达到机器精度．非对称辐射驱动模型 以及非对称网格剖分条 件下的数值模拟等算例也取得了较好的结果． 关键词：辐射输运；离散纵标格式；柱坐标系；球对称 MR (2000)主题分类：65M99 1．引 言 辐射输运方程描述 了光子在真空或介质 中的传输、散射 以及与介质的相互作用，在 天体物理、高能量密度物理、核物理、生物医学等领域中有重要应用．辐射输运方程的 计算存在维数高、计算量大等难点，其数值计算方法一直都是偏微分方程数值解研究领 域的热 门课题．求解输运方程 的数值计算方法可 以分为确定论方法和随机模拟方法两类． 确定论方法主要有离散纵标方法、球谐函数方法等；随机模拟方法即蒙特卡罗方法．离 散纵标方法由Calson在1953年计算一维平板问题时首次提出[1]’并于1965年和Lathrop将 该方法扩展到二维 2【]．该方法的主要特点是采用选定的离散方向代替连续的角度变量， 其数值算法简洁、程序实现简单，是求解输运方程最常用的方法之一．很多研究人员对离 散纵标方法进行了改进和完善，包括离散纵标方程的空间离散格式 3【_6】，离散方向求积组 的选取方法 [，引，求解离散纵标方程 的源迭代加速方法 [。1，等等． 虽然离散纵标方法具有很多优点，且应用广泛，但是该方法 同样存在射线效应、假散 射、源迭代收敛速度慢等缺点．其中，射线效应是离散纵标方法的一个特有缺陷，一种明 显的表现是利用离散纵标方法求解二维柱坐标系下辐射输运方程模拟一维球对称模型 时，会 出现无法保持球对称性的问题．最直接 的减弱射线效应的方法就是增加离散方 向 ’2015年1月27日收到． )基金项 目： 国家 自然科 学基金、 国家高技 术研 究发展计划 (2012AA01A303)、国家基础研究项 目(2011CB309702)资助项 目． 300 计 算 数 学 2015庄 数，使离散方 向分布更接近连续方 向分布．但是 由于扰动的衰减与离散方 向数的开方成 正比，简单的依靠增加离散方 向数减弱射线效应会带来 巨大的计算量，在实际应用 中并 不可取． 针对二维柱坐标系辐射输运方程的离散纵标格式计算球对称模型时无法保持球对 称性的问题，本文给出了格式无法保持球对称性的原因分析，并且 由此提出了空间基于 柱坐标系、方 向基于球坐标系的辐射输运方程，证 明了针对该方程 的离散纵标格式在球 矾矾 ㈣ ㈣ 对称网格剖分条件下能够保持一维球对称性的充