间断有限元方法求解-维非平衡辐射扩散方程

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

间断有限元方法求解-维非平衡辐射扩散方程

张荣培^1,3, 蔚喜军², 崔霞², 冯涛³

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

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

Discontinuous Finite Element Method for 1D Non-equilibrium Radiation Diffusion Equations

ZHANG Rongpei^1,3, YU Xijun², CUI Xia², FENG Tao³

1. School of Sciences, Liaoning ShiHua University, Fushun 113001, China;
2. Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing 100088, China;
3. Graduate School, China Academy of Engineering Physics, Beijing 100088, China

Received:2011-11-07 Revised:2012-02-03 Online:2012-09-25 Published:2012-09-25

摘要/Abstract

摘要： 研究一维非平衡辐射扩散方程的数值方法.通过求解间断系数热传导方程的广义黎曼问题,得到种一带加权数值流量,基于该数值流量构造了一类新型的间断有限元方法.在时间离散上采用向后Euler方法,形成的非线性方程组采用Picard迭代求解.数值试验表明该方法具有捕捉大梯度的能力,而且能适应扩散系数间断的情形.

关键词: 非平衡辐射扩散方程, 间断有限元, 间断系数

Abstract: We discuss numerical simulation of one-dimensional non-equilibrium radiation diffusion equations.A weighted numerical flux between adjacent grid cells is obtained by solving heat conduction equation with discontinuous coefficient.With this numerical flux of diffusive generalized Riemann problem(dGRP),a discontinuous finite element method is proposed for radiation diffusion equations. A backward Euler time diseretization is applied for semi-discrete form and a Picard iteration is used to solve nonlinear system of equations.Numerical results demonstrate that the method has a capability of capturing strong gradients and can be accommodated to discontinuous diffusion coefficient.

Key words: non-equilibrium radiation diffusion equation, discontinuous finite element method, discontinuous coefficient

中图分类号:

O241.82

张荣培, 蔚喜军, 崔霞, 冯涛. 间断有限元方法求解-维非平衡辐射扩散方程[J]. 计算物理, 2012, 29(5): 641-646.

ZHANG Rongpei, YU Xijun, CUI Xia, FENG Tao. Discontinuous Finite Element Method for 1D Non-equilibrium Radiation Diffusion Equations[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2012, 29(5): 641-646.

[1]	高普阳. 聚合物充填过程中裹气现象的数值模拟[J]. 计算物理, 2021, 38(6): 693-706.
[2]	段茂昌, 蔚喜军, 陈大伟, 黄朝宝, 安娜. DG方法求解可压缩气固两相流动[J]. 计算物理, 2017, 34(6): 631-640.
[3]	赵国忠, 蔚喜军, 郭怀民. 二维Lagrangian坐标系下可压气动方程组的间断Petrov-Galerkin方法[J]. 计算物理, 2017, 34(3): 294-308.
[4]	郭虹平, 欧阳洁. 气液两相流的间断有限元模拟[J]. 计算物理, 2015, 32(2): 160-168.
[5]	赵国忠, 蔚喜军, 李珍珍. 多介质流模拟的Runge-Kutta控制体积间断有限元方法[J]. 计算物理, 2014, 31(3): 271-284.
[6]	冯涛, 蔚喜军, 安恒斌, 崔霞, 吴迪, 李珍珍. Lax-Wendroff时间离散的自适应间断有限元方法求解三维可压缩欧拉方程[J]. 计算物理, 2013, 30(6): 791-798.
[7]	张荣培, 蔚喜军, 崔霞, 冯涛. 隐-显积分因子间断Galerkin方法求解二维辐射扩散方程[J]. 计算物理, 2012, 29(5): 647-653.
[8]	魏军侠, 袁光伟, 阳述林, 申卫东. 二维柱几何非定常中子输运方程基于格式的界面预估校正并行算法[J]. 计算物理, 2012, 29(2): 198-204.
[9]	陈大伟, 秦承森, 王裴, 孙海权, 蔚喜军. 凝聚介质中斜激波的反射[J]. 计算物理, 2011, 28(6): 791-796.
[10]	刘帅强, 欧阳洁, 阮春蕾. 非结构网格上Level Set方程的间断有限元解法[J]. 计算物理, 2011, 28(5): 649-658.
[11]	吴迪, 蔚喜军, 徐云. 局部时间步长间断有限元方法求解三维欧拉方程[J]. 计算物理, 2011, 28(1): 1-9.
[12]	吴迪, 蔚喜军. 自适应间断有限元方法求解三维欧拉方程[J]. 计算物理, 2010, 27(4): 492-500.
[13]	张磊, 袁礼. 龙格库塔间断有限元方法在计算爆轰问题中的应用[J]. 计算物理, 2010, 27(4): 509-517.
[14]	魏军侠, 阳述林, 傅连祥. 二维柱几何中子输运方程的并行区域分解方法[J]. 计算物理, 2010, 27(1): 1-7.
[15]	洪振英, 袁光伟. 粒子输运方程的线性间断有限元方法[J]. 计算物理, 2009, 26(3): 325-334.