基于变分原理的三维扩散方程离散化格式

计算物理 ›› 2003, Vol. 20 ›› Issue (4): 291-297.

基于变分原理的三维扩散方程离散化格式

陈光南¹, 李德元¹, 万正苏²

1. 北京虚用物理与计算数学研究所计算物理实验室, 100088;
2. 上海大学数学系, 上海 200436

收稿日期:2002-03-27 修回日期:2002-08-22 出版日期:2003-07-25 发布日期:2003-07-25
作者简介:陈光南(1937-),男,江苏常州,研究员,博士生导师,从事计算数学方面的研究,北京8009信箱12分箱100088.
基金资助:
国家高技术863计划惯性约束聚变委员会和计算物理国家重点实验室基金资助项目

Discrete Scheme Based on Variation Principle for Three Dimensional Diffusion Equation

CHEN Guang-nan¹, LI De-yuan¹, WAN Zheng-su²

1. Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing 100088, China;
2. Department of Mathematics, Shanghai University, Shanghai 200436, China

Received:2002-03-27 Revised:2002-08-22 Online:2003-07-25 Published:2003-07-25

摘要/Abstract

摘要： 研究三维扩散方程的数值模拟.在非正规六面体网格上,用变分原理建立扩散流形式的差分格式.将扩散流通量作为未知函数,求其泛函极值,并与原函数温度联立求解.详细推导了计算公式,并给出了数值试验结果.

关键词: 变分原理, 非正规六面体网格, 三维扩散方程, 差分格式

Abstract: Numerical simulation for three dimensional diffusion equation is studied.Using variation principle,the difference discrete scheme with diffusion flux form on the irregular hexahedral grid is constructed.The flux of diffusion as an unknown function is obtained by minimizing the functional and solved with the equation of temperature function simultaneously.Computational formulas are deduced and numerical test is given.

Key words: variation principle, irregular hexahedral grid, three-dimensional diffusion equation, difference scheme

中图分类号:

O241.82

陈光南, 李德元, 万正苏. 基于变分原理的三维扩散方程离散化格式[J]. 计算物理, 2003, 20(4): 291-297.

CHEN Guang-nan, LI De-yuan, WAN Zheng-su. Discrete Scheme Based on Variation Principle for Three Dimensional Diffusion Equation[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2003, 20(4): 291-297.

[1]	曹启伟, 肖德龙, 杨显俊, 王建国. 磁瑞利-泰勒不稳定性非线性演化数值模拟[J]. 计算物理, 2021, 38(1): 5-15.
[2]	叶珍宝, 周海京. 基于高阶叠层矢量基函数的E-H时域有限元方法分析谐振腔和波导结构[J]. 计算物理, 2016, 33(3): 333-340.
[3]	王林雪, 宗谨, 王雪玲, 石玉仁. mBBM方程的双扭结孤立波及其动力学稳定性[J]. 计算物理, 2016, 33(2): 212-220.
[4]	况晓静, 王道平, 张量, 吴先良, 沈晶, 沈勐. 基于紧致差分格式的高效时域有限差分算法[J]. 计算物理, 2014, 31(1): 91-95.
[5]	左风丽, 崔霞, 袁光伟. 热传导方程三层并行差分格式初始条件的计算[J]. 计算物理, 2011, 28(4): 488-492.
[6]	孙洪广, 陈文, 蔡行. 空间分数阶导数“反常”扩散方程数值算法的比较[J]. 计算物理, 2009, 26(5): 719-724.
[7]	刘国昭, 张树道. 气相爆轰高阶中心差分-WENO组合格式自适应网格方法[J]. 计算物理, 2008, 25(4): 387-395.
[8]	黄志祥, 沙威, 吴先良, 陈明生. 基于多步高阶差分格式求解时域Maxwell方程[J]. 计算物理, 2008, 25(3): 263-268.
[9]	胡彦梅, 陈建忠, 封建湖. 双曲型守恒律的一种五阶半离散中心迎风格式[J]. 计算物理, 2008, 25(1): 29-35.
[10]	杭旭登, 李敬宏, 袁光伟. 大长宽比扭曲网格上辅助网格差分目的[J]. 计算物理, 2007, 24(3): 268-276.
[11]	欧平, 马汉东, 汪翼云. 任意马赫数非定常流动数值模拟的统一算法[J]. 计算物理, 2007, 24(2): 166-170.
[12]	刘伟, 袁益让. 二维半导体问题在复合网格上的有限差分格式[J]. 计算物理, 2006, 23(6): 721-730.
[13]	郭晓虎, 张林波. 求解可压Navier-Stokes方程的对称超紧致差分格式及其并行算法[J]. 计算物理, 2006, 23(3): 281-289.
[14]	段毅, 杨永. 混合通量差分格式的比较[J]. 计算物理, 2006, 23(3): 355-360.
[15]	涂国华, 罗俊荣. 利用通量限制思想改进紧致格式[J]. 计算物理, 2005, 22(4): 329-336.