隐-显积分因子间断Galerkin方法求解二维辐射扩散方程

计算物理 ›› 2012, Vol. 29 ›› Issue (5): 647-653.

隐-显积分因子间断Galerkin方法求解二维辐射扩散方程

张荣培¹, 蔚喜军², 崔霞², 冯涛²

1. 辽宁石油化工大学理学院, 抚顺 113001;
2. 北京应用物理与计算数学研究所计算物理实验室, 北京 100088

收稿日期:2011-12-06 修回日期:2012-04-04 出版日期:2012-09-25 发布日期:2012-09-25
作者简介:张荣培(1978-),男,博士,讲师,主要从事计算数学研究,E-mail:rongpeizhang@163.com
基金资助:
国家自然科学基金(11171038,11261035)资助项目

Implicit-explicit Integration Factor Discontinuous Galerkin Method for 2D Radiation Diffusion Equations、

ZHANG Rongpei¹, YU Xijun², CUI Xia², FENG Tao²

1. School oySciences, Liaoning ShiHua University, Fushtm 113001, China;
2. National Key Laboratory of Science and Technology on Computational Physics, Institute ofApplied Physics and Computational Mathematics, Beijing 100088, China

Received:2011-12-06 Revised:2012-04-04 Online:2012-09-25 Published:2012-09-25

摘要/Abstract

摘要： 提出一种求解二维非平衡辐射扩散方程的数值方法.空间离散上采用加权间断Galerkin有限元方法,其中数值流量的构造采用一种新的加权平均;时间离散上采用隐-显积分因子方法,将扩散系数线性化,然后用积分因子方法求解间断Galerkin方法离散后的非线性常微分方程组.数值试验中在非结构网格上求解了多介质的辐射扩散方程.结果表明:对于强非线性和强耦合的非线性扩散方程组,该方法是-种非常有效的数值算法.

关键词: 二维辐射扩散方程, 间断有限元, 加权平均, 隐-显积分因子方法, 非结构网格

Abstract: A numerical method is developed for two-dimensional nonequilibrium radiation diffusion equations.Discontinuous Galerkin method is applied in spatial diseretization in which numerical flux is constructed with weighted flux averages.Implicit-explicit integration factor method for time discretization is applied to nonlinear ordinary differential equations which is obtained with discontinuous Galerkin method. Radiation diffusion equations with multiple materials are solved on unstructured grids in numerical tests.It demonstrates that the method is effective for high nonlinear and tightly coupled radiation diffusion equations.

Key words: two-dimensional radiation diffusion equation, discontinuous Galerkin finite element method, weighted averages, implicit-explicit integration factor method, unstructured grid

中图分类号:

O241.82

张荣培, 蔚喜军, 崔霞, 冯涛. 隐-显积分因子间断Galerkin方法求解二维辐射扩散方程[J]. 计算物理, 2012, 29(5): 647-653.

ZHANG Rongpei, YU Xijun, CUI Xia, FENG Tao. Implicit-explicit Integration Factor Discontinuous Galerkin Method for 2D Radiation Diffusion Equations、[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2012, 29(5): 647-653.

[1]	高普阳. 聚合物充填过程中裹气现象的数值模拟[J]. 计算物理, 2021, 38(6): 693-706.
[2]	艾邦成, 张亮, 陈智. 低耗散非结构网格多维限制器[J]. 计算物理, 2018, 35(5): 545-553.
[3]	段茂昌, 蔚喜军, 陈大伟, 黄朝宝, 安娜. DG方法求解可压缩气固两相流动[J]. 计算物理, 2017, 34(6): 631-640.
[4]	林文洲, 林忠, 刘全. 非结构多边形网格滑移线开穴算法[J]. 计算物理, 2017, 34(3): 273-282.
[5]	赵国忠, 蔚喜军, 郭怀民. 二维Lagrangian坐标系下可压气动方程组的间断Petrov-Galerkin方法[J]. 计算物理, 2017, 34(3): 294-308.
[6]	郭虹平, 欧阳洁. 气液两相流的间断有限元模拟[J]. 计算物理, 2015, 32(2): 160-168.
[7]	赵国忠, 蔚喜军, 李珍珍. 多介质流模拟的Runge-Kutta控制体积间断有限元方法[J]. 计算物理, 2014, 31(3): 271-284.
[8]	冯涛, 蔚喜军, 安恒斌, 崔霞, 吴迪, 李珍珍. Lax-Wendroff时间离散的自适应间断有限元方法求解三维可压缩欧拉方程[J]. 计算物理, 2013, 30(6): 791-798.
[9]	张荣培, 蔚喜军, 崔霞, 冯涛. 间断有限元方法求解-维非平衡辐射扩散方程[J]. 计算物理, 2012, 29(5): 641-646.
[10]	魏军侠, 袁光伟, 阳述林, 申卫东. 二维柱几何非定常中子输运方程基于格式的界面预估校正并行算法[J]. 计算物理, 2012, 29(2): 198-204.
[11]	陈大伟, 秦承森, 王裴, 孙海权, 蔚喜军. 凝聚介质中斜激波的反射[J]. 计算物理, 2011, 28(6): 791-796.
[12]	雷国东, 李万爱, 任玉新. 求解可压缩流的高精度非结构网格WENO有限体积法[J]. 计算物理, 2011, 28(5): 633-640.
[13]	刘帅强, 欧阳洁, 阮春蕾. 非结构网格上Level Set方程的间断有限元解法[J]. 计算物理, 2011, 28(5): 649-658.
[14]	明平剑, 段文洋. 液舱横荡非结构网格高阶格式数值模拟[J]. 计算物理, 2011, 28(4): 507-514.
[15]	张来平, 刘伟, 贺立新, 赫新, 邓小刚. 基于静动态重构的高阶DG/FV混合格式在二维非结构网格中的推广[J]. 计算物理, 2011, 28(2): 188-198.