长江科学院院报 ›› 2020, Vol. 37 ›› Issue (7): 175-182.DOI: 10.11988/ckyyb.20190907

• 独立覆盖形行法专栏 • 上一篇    

一维对流扩散方程的独立覆盖分析方法

刘亚军1, 龚亚琦1,2, 苏海东1,2   

  1. 1.长江科学院 材料与结构研究所,武汉 430010;
    2.长江科学院 水利部水工程安全与病害防治工程技术研究中心,武汉 430010
  • 收稿日期:2019-07-29 出版日期:2020-07-01 发布日期:2020-08-06
  • 通讯作者: 苏海东(1968-),男,湖北武汉人,教授级高级工程师,博士,从事水工结构数值分析工作和计算方法研究。E-mail:suhd@mail.crsri.cn
  • 作者简介:刘亚军(1994-),男,江西赣州人,硕士研究生,从事水工结构数值分析工作和计算方法研究。E-mail:liuyj963@163.com
  • 基金资助:
    国家自然科学基金项目(51409012)

Numerical Method for Solving One-dimensional Convection-Diffusion Equation Based on Independent Covers

LIU Ya-jun1, GONG Ya-qi1,2, SU Hai-dong1,2   

  1. 1. Material and Engineering Structure Department, Yangtze River Scientific Research Institute, Wuhan 430010, China;
    2. Research Center on Water Engineering Safety and Disaster Prevention of Ministry of Water Resources, Yangtze River Scientific Research Institute, Wuhan 430010, China
  • Received:2019-07-29 Online:2020-07-01 Published:2020-08-06

摘要: 针对现有的各种数值方法在求解一维对流扩散方程时容易出现的数值振荡、假扩散等计算稳定性和计算精度不足问题,提出应用独立覆盖流形法进行数值求解的新思路,即分区的多项式级数逼近。基于标准的伽辽金法推导一维对流扩散方程的独立覆盖流形法求解公式。采用场变量的一阶导数在独立覆盖之间的窄条形覆盖重叠区域是否连续的后验误差估计方法,通过覆盖加密和级数升阶的h-p型混合自适应进行自动求解。给出的稳态和非稳态分析算例结果表明:分区级数的数值解稳定地逼近于精确解,最终两者很好地吻合;对于对流占优问题,自适应求解可以有效避免数值振荡。另外还尝试了将数值解代回微分方程计算残差作为误差指标,如果能使微分方程逐点满足,那么将是对数值解最严格的误差判断。

关键词: 对流扩散方程, 数值振荡, 自适应分析, 数值流形法, 独立覆盖

Abstract: In solving one-dimensional convection-diffusion equations, present numerical methods are prone to suffer from stability and accuracy problems caused by numerical oscillation and pseudo-diffusion. In view of this, an idea of applying Numerical Manifold Method (NMM) based on independent covers (the approximation using polynomial series piecewise-defined) to the numerical solution is proposed. The solution formula of the one-dimensional convection-diffusion equation is derived based on the standard Galerkin method. The posterior error estimation method about the continuity of the first-order derivative of field variable in the narrow overlapping area between independent covers is used for the automatic solving by h-p hybrid self-adaptive analysis with mesh refinement and ascending series order. The results of the steady-state and unsteady-state analysis examples show that the numerical solution of the piecewise-defined series steadily approximates and finally well fits the exact solution. For the convection-dominated problem, the adaptive solution effectively avoids numerical oscillation. In addition, the error index of the residual by substituting the numerical result back to the differential equation is successfully attempted. If the differential equation is solved point by point, the method would the most stringent error judgment for the numerical solution so far.

Key words: convection-diffusion equation, numerical oscillation, self-adaptive analysis, Numerical Manifold Method (NMM), independent covers

中图分类号: