预处理JFNK方法求解非平衡辐射扩散方程组

计算物理 ›› 2013, Vol. 30 ›› Issue (4): 483-490.

预处理JFNK方法求解非平衡辐射扩散方程组

冯涛^1,2, 蔚喜军³, 安恒斌³, 张荣培²

1. 中国科学技术大学数学科学学院, 合肥 230052;
2. 中国工程物理研究院研究生部, 北京 100088;
3. 北京应用物理与计算数学研究所计算物理实验室, 北京 100094

收稿日期:2012-10-10 修回日期:2013-01-24 出版日期:2013-07-25 发布日期:2013-07-25
作者简介:冯涛(1984-),male,doctor,major in computational mathematics,E-mail:fengta02@mail.ustc.edu.cn
基金资助:
Supported by National Natural Science Foundation of China(11171039,11171038,11271054);National Basic Research Program of China(2011CB309702,2011CB309703);Science Foundation of CAEP(2012B0202026,2010A0202010);Foundation of National Key Laborotory of Science and Technology on Computation Physics and Defence Industrial Technology Development Program(B1520110011)

Preconditioned Jacobian-free Newton-Krylov Methods for Nonequilibrium Radiation Diffusion Equations

FENG Tao^1,2, YU Xijun³, AN Hengbin³, ZHANG Rongpei²

1. School of Mathmatical Sciences, University of Science and Technology of China, Hefei 230052, China;
2. Graduate School of China Academy Engineering Physics, Beijing 100088, China;
3. National Key Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing 100094, China

Received:2012-10-10 Revised:2013-01-24 Online:2013-07-25 Published:2013-07-25
Supported by:
Supported by National Natural Science Foundation of China(11171039,11171038,11271054);National Basic Research Program of China(2011CB309702,2011CB309703);Science Foundation of CAEP(2012B0202026,2010A0202010);Foundation of National Key Laborotory of Science and Technology on Computation Physics and Defence Industrial Technology Development Program(B1520110011)

摘要/Abstract

摘要： 分别采用四种半隐式离散方法构造预处理.针对一维辐射扩散方程组,采用预处理的Jacobian-free Newton-Krylov(PJFNK)求解.数值结果表明预处理方法能够很好地改进JFNK方法的收敛行为.

关键词: 非平衡辐射扩散, Newton-Krylov方法, 线性化, 预处理

Abstract: Four semi-implicit discretization schemes are used to construct preconditioners.And preconditioned Jacobian-free Newton-Krylov (JFNK) are presented to solve one-dimensional problems.Numerical results show that the preconditioning methods improve the convergence behavior of JFNK method dramatically.

Key words: nonequilibrium radiation diffusion, Newton-Krylov method, linearization, preconditioning

中图分类号:

O241.7

冯涛, 蔚喜军, 安恒斌, 张荣培. 预处理JFNK方法求解非平衡辐射扩散方程组[J]. 计算物理, 2013, 30(4): 483-490.

FENG Tao, YU Xijun, AN Hengbin, ZHANG Rongpei. Preconditioned Jacobian-free Newton-Krylov Methods for Nonequilibrium Radiation Diffusion Equations[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2013, 30(4): 483-490.

[1]	李凌霄, 翟传磊, 谢辉, 施意. 一种求解三维热辐射输运方程的整体预处理迭代方法及并行计算[J]. 计算物理, 2021, 38(3): 269-279.
[2]	李凌霄. 不可压缩流基于块预处理的并行有限元计算[J]. 计算物理, 2018, 35(2): 151-160.
[3]	刘中玉, 张明锋, 郑冠男, 杨国伟. 基于预处理HLLEW格式的全速域数值算法[J]. 计算物理, 2016, 33(3): 273-282.
[4]	杨容, 杭旭登, 翟传磊, 李双贵, 齐进, 李敬宏. 激光光路追踪模拟中的高精度求交算法[J]. 计算物理, 2015, 32(2): 207-213.
[5]	甘文彪, 周洲. 基于层流动能湍流模型的数值模拟方法[J]. 计算物理, 2013, 30(2): 169-179.
[6]	张荣培, 蔚喜军, 崔霞, 冯涛. 间断有限元方法求解-维非平衡辐射扩散方程[J]. 计算物理, 2012, 29(5): 641-646.
[7]	安恒斌, 莫则尧. JFNK方法迭代过程与物理约束[J]. 计算物理, 2012, 29(5): 654-660.
[8]	杨帅帅, 孙玉发. 应用AWE技术和等效偶极子法快速计算目标宽带RCS[J]. 计算物理, 2012, 29(3): 406-410.
[9]	韩忠华, 宋文萍, 乔志德. 一种隐式预处理方法及其在定常和非定常流动数值模拟中的应用[J]. 计算物理, 2009, 26(5): 679-684.
[10]	范宣华, 吴瑞安, 郝志明, 何颖波. 基于Tahoe框架的某夹具并行计算[J]. 计算物理, 2009, 26(5): 699-702.
[11]	欧平, 马汉东, 汪翼云. 任意马赫数非定常流动数值模拟的统一算法[J]. 计算物理, 2007, 24(2): 166-170.
[12]	马东军, 孙德军, 尹协远. 高阶谱元区域分解算法及其在流动稳定性中的应用[J]. 计算物理, 2007, 24(1): 7-12.
[13]	李晨, 吴雄华. Navier-Stokes方程的线性化微分求积法[J]. 计算物理, 2006, 23(1): 10-18.
[14]	谷同祥, 戴自换, 杭旭登, 符尚武, 刘兴平. 二维三温能量方程组的高效代数解法[J]. 计算物理, 2005, 22(6): 1-8.
[15]	吴建平, 刘兴平, 王正华, 戴自换, 李晓梅. 二维三温能量方程组离散求解的两个新预处理技术[J]. 计算物理, 2005, 22(4): 283-291.