基于矩近似的辐射流体力学计算方法研究.ppt

最近的研究成果 自由流问题（free streaming limit） 最近的研究成果 Su-Olson problem 谢谢！ 辐射流体力学数值计算方法研究 闫 凯 2013.11 内容 辐射流体方程 国内外研究进展 最近的研究成果 辐射流体方程 辐射流体方程的一般表达式： 其中： 辐射流体方程 辐射强度是时间、空间、角度和频率的函数。 频率问题:多群、灰体。 多群和灰体方程之间的求解方法差异 不大，下面的讨论仅限于灰体近似下的辐射流体方程。 其中： 辐射流体方程 求解难点： 辐射流体方程内含三种重要的动力学尺度：物质波波速、辐射波波速以及源项相互作用的速度； 不论做何种近似，这三种尺度始终是存在的。该问题是限制辐射流体力学求解效率的关键。 相对论效应问题：关于该问题的讨论可以见Caster J I和Mihalas D等人的著作。 辐射流体方程不论统一在何种坐标系下，客观上都增加了方程组的求解难度。 其中： 辐射流体方程 求解难点： 不透明度问题。 由于介质不透明度和介质的温度和密度相关，不透明度的变化会极大影响辐射输运过程，辐射输运带来的介质密度和温度的变化又会进一步影响介质不透明度。目前关于该问题的讨论不多。 局部热力学平衡问题。 局部热力学平衡条件的适用范围？在温度急剧变化的区域是否适用？ 其中： 辐射流体方程 辐射输运模型： 光子波尔兹曼方程 直接求解该方程的难点在于角度问题的处理，目前多采用SN方法或者Mc方法。 SN方法中空间步长的选取是值得商榷的问题。按照JI Caetor的说法，空间步长不能远大于光子平均自由程，否则讨论角度问题毫无意义。 其中： 辐射流体方程 辐射输运模型： 矩近似方程 目前最常用的是一阶矩近似方法。该方法又称为Eddington近似、M1近似。目前的火球计算即采用该方法。如果对一阶矩近似进一步简化，可以得到P1近似方程，限流扩散方程（FLD). 矩方法消除了角度的影响，对空间步长不敏感 。但是需要添加一个封闭条件。封闭条件的选取往往要带来一定的误差。 其中： 辐射流体方程 辐射输运模型： 几个概念 辐射流体方程 数值方法： 根据辐射输运模型的不同，具体的数值方法可以分为SN方法、PN方法和MC方法。在这里只讲一阶矩近似方程（下面简称为M1方程）的数值方法。 由于辐射流体方程组的多尺度特性，直接求解是非常困难的，基本思想是分裂、隐式、迭代。具体的数值格式在第二部分详细说明。 其中： 国内外相关领域工作 国内的应用： 辐射流体力学在国内多应用于惯性约束聚变或者激光辐照问题； 其采用的基本方程： 三温方程及其演变形式；光子波尔兹曼方程与流体方程耦合求解； 目前形成的计算程序： 基于JASMIN框架下的LARED-H,RH-2D； 国内的天体物理领域很少有关于RHD的讨论。 国内外相关领域工作 国外的应用： 辐射流体力学在国外更多应用于天体物理的数值模拟； 目前形成的计算程序非常多，很多程序已经实现了三维、AMR、并行化等要求； 天体物理中辐射输运模型一般采用FLD近似或者M1近似； 目前关于RHD方程讨论比较热门的话题有：辐射输运模型问题、刚性源项问题、相对论效应引起的坐标系问题、数值格式问题； 关于核爆火球数值模拟的文章依旧很少，仅有的几篇全部来自LA报告，且绝大部分发表时间在2000年以前。 国内外相关领域工作 计算程序 辐射输运模型 方程组维数 是否AMR 坐标系 CASTRO FLD 3D Yes Mixed-frame Orion FLD 3D Yes Mixed-frame RAMSES FLD 3D Yes Moved-frame CRASH FLD 3D Yes Moved-frame ZEUS FLD 2D No Moved-frame HERACLES M1 3D No Moved-frame Nike M1 1D No Mixed-frame 国内外相关领域工作 M1近似的数值格式： HERACLES 国内外相关领域工作 M1近似的数值格式： Nike Step 1:求解带源项的辐射输运方程(implicit HLLE scheme) Step2:求解带源项的流体方程（semi-impli