本文研究了二维高雷诺数情形下,非定常不可压Navier-Stokes方程的Euler隐/显格式子格涡旋粘性非协调有限元方法。隐式处理线性项,避免了时间步长的苛刻限制,显式处理非线性项,使得对所有时间层求解时,系统矩阵为同一常数矩阵;时间项做向后Euler差分离散,空间用C-R非协调有限元逼近,构造子格涡旋粘性有限元方法,克服了在情形下Galerkin有限元方法的不稳定现象。本文改善了稳定性下对时间步长的限制,并给出了不依赖粘性系数的速度和压力误差估计。
参考文献