二维三温热传导方程求解中的非线性迭代初值选取 - 北京应用物理与 .pdf

第 24 卷 第 2 期 计 　　算 　　物 　　理 Vol. 24 ,No. 2 　 2007 年 3 月 　CHINESE JOURNAL OF COMPUTATIONAL PHYSICS 　 Mar. , 2007 ( ) [文章编号] 　1001246X 2007 二维三温热传导方程求解中的非线性迭代初值选取 1 2 3 安恒斌 , 　莫则尧 , 　徐小文 ( 1 北京应用物理与计算数学研究所 计算物理实验室 , 北京　100088 ; 2 北京应用物理与计算数字研究所 高性能计算中心 , 北京　100088 ; 　3 中国工程物理研究院研究生部, 北京　100088) [摘 　要] 　由于二维三温热传导方程具有很强的非线性特性, 因此采用全隐格式对该方程离散后 , 所得非线性 代数方程组的求解将变得非常困难. 针对二维三温热传导方程离散所得非线性代数方程组的迭代求解 , 提出了 一种有效的选取初值的方法. 对两种不同性质的介质进行数值实验 , 结果表明 , 所设计的初值选取方法不仅大 大提高了计算效率 , 而且能够降低非线性解法器对时间步长的影响. [ 关键词] 　三温热传导方程 ; 非线性方程组 ; 非精确Newton 法 ; 迭代初值 [ 中图分类号] 　O2417 [文献标识码] 　A 0 　引言 在惯性约束聚变( ) ICF 数值模拟中 ,对二维辐射流体力学方程组的求解具有非常重要的意义. 而在二维 [9 ] 辐射流体力学方程组的数值求解中 ,二维三温能量方程的求解占整个计算量的 80 %以上 , 因此 ,设计有效 的求解二维三温能量方程的数值方法 ,对于提高整个数值模拟有着非常重要的作用. 二维三温能量方程反映 电子、离子和光子在介质中的能量传播规律. 随着时间的推进 ,既有能量扩散 ,又有光子和电子 , 电子和离子 之间的能量交换 ,三者的温度耦合在一起 ,是需要求解的未知量. 对二维三温能量方程进行适当简化 ,就得到 二维三温热传导方程. 这一方程在本质上与二维三温能量方程是相同的[10 ] ,其定义为 Te 1 Δ ( Δ ) ω ( ) ω ( ) c - ·K T = T - T + T - T , Ve e e ei i e er r e t ρ