A numerical solution method with high accuracy is proposed herein for the instability problems of a two -dimensional steady convection-dominated convection-diffusion equation.The new scheme is a stabilized spectral el-ement method based on the consistent approximate upwind method and the Chebyshev spectral element method.The feasibility of the method is verified by numerical examples.Moreover,the relation among the calculation accuracy, convergence rate,boundary layer approximation, and the order of the element interpolation is discussed.The re-sults show that compared with the spectral element method,the stabilized spectral element method greatly enlarges the stable solution region for the convection-diffusion equation.The numerical oscillations in the computational do-main and near the boundary layers in the convection-diffusion problem are completely eliminated.In addition,the internal and external boundary layers are excellently approximated.The computational accuracy and the effect of approximation to the boundary layers are both rapidly enhanced with the increase of the element interpolation order.%为了求解对流项占优的对流扩散方程的非稳定性问题,本文提出了二维稳态对流扩散方程的一种高精度求解方法.该方法将一致逼近迎风方法和Chebyshev谱元方法相结合,形成了稳定化谱元方法.通过数值算例对该方法的可行性进行了验证,讨论了计算精度、收敛速度及边界层逼近和单元插值阶数的关系.研究表明:和谱元方法相比,稳定化谱元方法扩大了对流扩散方程的稳定求解区域;消除了计算区域内部及边界层附近的数值振荡,对内部和外部边界层都进行了很好的逼近;单元插值阶数的增加使计算精度及边界层的逼近效果都获得了迅速的提高.
展开▼
