使用标准的有限元方法求解非定常Navier-Stokes方程所得速度误差常受压力误差影响，且误差随粘性系数的减少而增大。为了增强压力的鲁棒性，本文引入grad-div稳定项，以提高近似解的精度，提出数值求解非定常Navier-Stokes方程的并行两水平grad-div稳定有限元算法，其时间和空间离散分别采用隐式Euler格式和Galerkin有限元方法。首先在全局粗网格上求解非线性grad-div稳定问题，然后在相互重叠的细网格子区域上并行求解grad-div稳定问题，以校正粗网格解。最后给出数值实验验证理论分析的正确性和算法的有效性。
