2011-郑州-年会-修改稿.pdf

辐射扩散方程非线性迭代方法 袁光伟 岳晶岩 北京应用物理与计算数学研究所 内容摘要 研究背景 迭代方法  数值算例  总结讨论 2 研究背景 辐射流体力学中能量输运由如下非线性热传 导方程描述: T ( , , ) ( , , ) C  x t T T S x t T v   t 热传导系数的特点： 非线性，间断，由数据库形式插值给出 3 两温模型 E 4 T 4 K T T E (DE )  (T E ) (  )  (  ) a a t t 其中E 辐射能量，T 物质温度， 1 z 3 ( ) D T  (T ) (光子吸收截面) （不限流） a 3 3a (T ) T 1 D T E E 2 5 2 ( , , ) 3a (T )  E E  （限流）K (T ) 110 T z 为原子质量数，与物质有关 z 的取值会有间断，从而扩散系数会产生间断。 4 研究背景：计算效率 高度非线性，多时间尺度  稳定性要求采用全隐格式，如何设计？ 大变形网格，间断系数  非线性迭代方法，如何线性化？  线性迭代 T n1 T n n1 n1 x T T  (( , ) ) 0 h h t 5 研究背景：常用的迭代方法  算