一维对流.docVIP

求解一维稳态对流扩散方程的一种 基于交错网格的有限体积法 袁晓彬1 李仕华2 ,袁富贵1，郭保生1 (1.广东白云学院 艺术与建筑学院，广州 610065; 2. 西京学院，西安学院 710123) 摘 要：对于一维稳态对流扩散方程，本文提出了一种新的有限体积格式。通过三个相邻节点组成有限体积，采用交错网格，对场变量用二阶拉格朗日插值离散，将一维稳态对流扩散方程离散成线性方程组求解。该法具有守恒性、无输运性、在一定条件下有界、二次代数精度的特点。通过与数值算例的解析解对比分析，表明该法是准确可靠的。 关键词：一维稳态对流扩散方程 有限体积 交错网格 二阶拉格朗日方程插值 中图分类号：O351.文献标识码：A 0 引言 有限体积法是解决流体力学问题的重要计算方法。有限体积法自提出后得到了快速发展。Godunov[1]提出了一维对流扩散方程的迎风差分格式，Leonard[2]则给出了QUICK格式。有限体积法也被运用于固体力Zhang Xiaofeng等[5]运用待定系数法,将一维纯对流下的HAUC2格式推广应用到一维对流扩散方程，给出了一种高精度算法。吴腾等[6]对传统迎风格式进行了改进，运用待定系数法,引入量纲-参数来抑制格式的频散和耗散,该格式具有较高的精度。陈翠霞等[7]针对常系数对流扩散方程,基于微分算子分裂算法思想,构造了一种新的一维对流扩散方程的数值求解格式。 本文基于交错网格，通过三个相邻节点组成有限体积，对场变量用二阶拉格朗日插值, 收稿日期: 2013 -03-01．修改稿收到日期: 作者简介: 袁晓彬( 1970—) ，男，博士，讲师 （Email:cocoyxbyxb@）（１） 式中为流体密度，为场变量，为在方向的流动速度，为扩散系数，为源项。满足连续性方程（2） 采用图1的控制容积，其中为控制容积主节点，、为控制容积边节点。 图一 有限体积网格示意图 Fig.1 Finite volume mesh diagram 对控制容积积分，有（3） 同理（4） 由奥高公式，控制容积内积分的对流扩散方程可写成（5） 即 （6） 类似地，连续性方程可写成（7） 即（8）对控制容积内三点、、场变量采用二阶拉格朗日插值有[8]（9） 式中，为的场变量，为的场变量，为的场变量。 其一阶导数为：（10） 取，，代入并化简得（11） 、处的一阶导数分别为：（12）（13） 将式（12）、（13）代入式（6）、（8）并令、、，，式（6）、（8）分别化为（14）（15） 整理得（16） 记（17）（18）（19）（20） 式（16）简化为（21） 式（21）即为最终的求解公式。 2 基本性质的讨论 （1）守恒性 在有限体积内采用相邻3节点拉格朗日插值，虽然采用交错二阶插值，但场变量在界面处在各节点是是协调的，从整个求解域来看扩散流动应是守恒的。 （2）有界性 由连续性方程知，，因此 ，在节点内部满足收敛性条件。若则是收敛的。得到。若满足该条件，算法是有界的。 （3）输运性 由二阶拉格朗日插值知处场变量对所有相邻节点的影响一样，没有反应处扩散与对流的区别，因此该法不具有输运特征。 （4）计算精度 由二阶拉格朗日插值场变量有三次截断精度，其导数有二次截断精度。参照Lamber[9]关于计算精度阶数定义，本法应具有二次截断精度。 3 算例 某场变量经过对流扩散从一维区域的点输运到点，输运控制方程为 流体密度，，扩散系数。已知此问题的解析解为 求：（1）当流速为，中间离散为5个节点网格时，在区域内的分布；（2）当流速为，中间离散为5个节点网格时，在区域内的分布； （3）当流速为，中间离散为20个节点网格时，在区域内的分布。 对问题1，中间取5个节点，由于，由图2可见数值解和精确吻合得很好；当取，中间取5个节点，数值解与精确解误差较大，数值解出现振荡，而通过节点加密，当中间取20个节点时，数值解与精确解吻合得相当好。由此可见，本法的精确性很好，计算结果令人满意。图二 不同问题下数值解与解析解的对比分析 Fig.2 The example numerical solutions comparison with analytic one under different conditions 4 结论 基于上述分析和计算，得出如下结论： （1）本文在原有有限体积基础上将控制容积扩大到三相邻节点上，对场变量利用二阶拉格朗日插值构造出了求解一维稳态对流扩散方程的一类新的算法。 （2）本法具有守恒性、当是有界的、无输运性、具有二次代数精度的特点。 （3）利用本法求解了一个一维稳态对流扩散方程问题。数值算例表明：本法有很好的计算精度，在小步长比较接近解析解，因而本法是正确可靠的。 参考文献