计算物理 ›› 2018, Vol. 35 ›› Issue (3): 269-274.DOI: 10.19596/j.cnki.1001-246x.7632

• • 上一篇    下一篇

求解对流-扩散-反应问题的改进有限积分法

孙婷婷, 韩赛赛, 许明田   

  1. 山东大学土建与水利学院 工程力学系, 山东 济南 250061
  • 收稿日期:2017-01-22 修回日期:2017-04-18 出版日期:2018-05-25 发布日期:2018-05-25
  • 作者简介:孙婷婷(1990-),女,硕士,主要从事流体力学对流扩散方程数值模拟方法研究,E-mail:stt2014@126.com
  • 基金资助:
    国家自然科学基金(11272187)资助项目

A Modified Finite Integration Method for Convection-Diffusion-Reaction Problems

SUN Tingting, HAN Saisai, XU Mingtian   

  1. Department of Engineering Mechanics, Shandong University, Jinan 250061, China
  • Received:2017-01-22 Revised:2017-04-18 Online:2018-05-25 Published:2018-05-25

摘要: 将求解偏微分方程的有限积分法应用于对流-扩散-反应问题,发现对于非对流占优的对流扩散问题,有限积分法的精度比QUICK法高一个数量级,比传统的有限体积法高两个数量级.处理对流占优的对流-扩散-反应问题时,对流项的离散时引进加权参数,通过调节该参数反映输运的方向性.结果表明这种改进的有限积分法的精度比传统的有限体积法至少高四个数量级,同时明显改进了原来的有限积分法的精度和稳定性.对于对流占优的对流-扩散-反应问题,即使采用粗网格,计算结果也未出现非物理振荡现象,表明改进的有限积分法具有很好的稳定性.

关键词: 对流-扩散-反应方程, 有限积分法, 有限体积法, 数值模拟

Abstract: Finite integration method is employed to solve convection-diffusion equations. For diffusion-dominated convection-diffusion problems, numerical results show that finite integration method is more than one and two orders of magnitude better in accuracy than QUICK scheme and central differencing scheme, respectively. We find that finite integration method can be improved by introducing a weighting parameter in numerical quadrature. By appropriately adjusting the parameter, modified finite integration method is at least more than four orders of magnitude better in accuracy than traditional finite volume method for convection-dominated flows. Furthermore, even coarse grids are utilized to divide flow domain modified finite integration method has not lead to unphysical oscillations in convection-dominated convection-diffusion problem, which demonstrates a good stability of the method.

Key words: convection-diffusion-reaction equation, finite integration method, finite volume method, numerical simulation

中图分类号: