求解对流－扩散－反应问题的改进有限积分法-计算物理.pdf

第３５卷 第３期 计　 　 算　 　 物　 　 理 Ｖｏｌ．３５，Ｎｏ．３ ２０１８年５月 　 ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＣＯＭＰＵＴＡＴＩＯＮＡＬ ＰＨＹＳＩＣＳ　 Ｍａｙ，２０１８ 文章编号：１００１⁃２４６Ｘ（２０１８）０３⁃０２６９⁃０６ 求解对流－扩散－反应问题的改进有限积分法 孙婷婷，　 韩赛赛，　 许明田 （山东大学土建与水利学院工程力学系，山东 济南　 ２５００６１） 摘　 要：将求解偏微分方程的有限积分法应用于对流－扩散－反应问题，发现对于非对流占优的对流扩散问题，有限 积分法的精度比ＱＵＩＣＫ法高一个数量级，比传统的有限体积法高两个数量级．处理对流占优的对流－扩散－反应问题 时，对流项的离散时引进加权参数，通过调节该参数反映输运的方向性．结果表明这种改进的有限积分法的精度比传 统的有限体积法至少高四个数量级，同时明显改进了原来的有限积分法的精度和稳定性．对于对流占优的对流－扩散－ 反应问题，即使采用粗网格，计算结果也未出现非物理振荡现象，表明改进的有限积分法具有很好的稳定性． 关键词：对流－扩散－反应方程；有限积分法；有限体积法；数值模拟 中图分类号：Ｏ３８ 文献标志码：Ａ ＤＯＩ：１０．１９５９６／ ｊ．ｃｎｋｉ．１００１⁃２４６ｘ．７６３２ ０　 引言 对流－扩散－反应方程在化学、物理及工程领域有着广泛应用．流体力学、传热学、化学反应动力学以及 ［１］ 环境科学中的众多现象都可以归结为求解对流－扩散－反应方程问题 ． 因此研究对流－扩散－反应方程的 高效数值解法具有重要的应用价值．求解对流－扩散－反应方程的数值方法主要有有限差分法（ＦＤＭ）、有限 元法（ＦＥＭ）和有限体积法（ＤＶＭ）［２－６］． 其中有限体积法广泛用于求解流体流动和传热问题，与其它数值计算方法相比较，有限体积法具有 显著的特点．其得到的离散方程能更好地保持原微分方程的守恒性，且离散方程中的各项物理意义明确， 离散方程形式规范．但是传统的高阶有限体积法没有体现对流输运的方向性，在处理对流占优的问题时， 数值结果往往出现非物理的振荡．虽然迎风格式弥补了这一缺陷，但传统的一阶迎风格式精度不高，容易 导致假扩散的情况．ＱＵＩＣＫ格式虽然提高了一阶迎风格式的精度，但其对边界条件的近似处理往往会影 响其精度． ［７］ 为了构造数值模拟激波等具有间断或速度（或其它量）的梯度剧烈变化的流体力学问题，Ｈａｒｔｅｎ 通过 ［８－９］ 通量修正法构造了二阶ＴＶＤ（Ｔｏｔａｌ Ｖａｒｉａｔｉｏｎ Ｄｉｍｉｎｉｓｈｉｎｇ）格式；Ｖａｎ Ｌｅｅｒ 提出在传统的ＭＵＳＣＬ格式中添 加限制器的方法用于数值模拟复杂的流体力学问题．另外基于有限差分法的高阶紧致格式也成为计算流体 力学的热点［１０－１１］． ［１２］ ［１３］ 近年来，研究者们基于积分方程理论提出了一些新的数值方法 ．ＮＦＤＳ格式 在求解对流扩散问题时提 ［１４］ 高了计算精度，且节约了计算时间．Ｘｕ等 建立了数值模拟磁流体自激发电现象的积分方程法，展现出很高的 ［１５］ 精度和稳定性．Ｗｅｎ等 提出了一种有限积分方法求解偏微分方程，该方法首先通过积分将偏微分方程化为 积分方程，然后应用数值方法离散积分方程．本文拟改进该方法，并应用于求解对流－扩散－反应问题． １　 对流－扩散问题 首先考虑如下一维对流扩散问题 ｄ ｄ ｄＴ