用于不可压缩流高阶有限元离散的静凝聚迭代罚方法

我们引入并分析了一种静凝聚迭代罚（SCIP）方法，用于求解以$p$阶Scott-Vogelius元离散的不可压缩流问题。虽然标准迭代罚方法常作为计算离散解的首选算法，但它在每次迭代中需要求解具有$\mathcal{O}(p^{d})$个未知量的线性系统。SCIP方法将该系统的规模缩减至$\mathcal{O}(p^{d-1})$个未知量，同时保持迭代罚方法的几何收敛速率。将SCIP应用于Kovasznay流与Moffatt涡表明其与理论吻合良好。