节块展开法求解对流扩散方程的稳定性和数值耗散特性分析

Stability and Numerical Diffusion Analysis of Nodal Expansion Method for Convection-diffusion Equation

  • 摘要: 本文研究了节块展开法求解对流扩散方程的稳定性和数值耗散特性。通过离散方程精确解和数值实验方法分析不同阶节块展开法的稳定性和数值耗散特性,并将其与有限体积法中的中心差分和一阶迎风格式进行对比。结果表明:偶数阶节块展开法的稳定是有条件的,即Peclet数(Pe)小于限值,且Pe限值会随展开阶数的增大而增大,其稳定性范围和精度均优于中心差分格式;奇数阶节块展开法是无条件稳定的,但随Pe的增大,数值耗散增大、计算误差增大,且当Pe大于一定值后,产生的数值耗散大于一阶迎风格式。

     

    Abstract: The stability and numerical diffusion analysis of nodal expansion method (NEM) for convection-diffusion equation was studied. The stabilities and numerical diffusion analyses for NEM with different order basis functions by exact solution of the discretization equation and numerical experiment results were done, and the results from NEM were compared with the results from both the center difference scheme and the first order upwind difference scheme of the finite volume method. The results show that the even order NEM is conditionally stable with Peclet number (Pe) being less than limit value,Pe limit value also increases with the order of expansion functions, and stability range and accuracy are better than those of center difference scheme. The odd order of NEM is unconditionally stable, however, with the increase of Pe, the numerical diffusion becomes lager and the calculation error becomes bigger. In addition, when Pe reaches certain value, the numerical diffusion even exceeds that of the first order upwind difference scheme.

     

/

返回文章
返回