提出了两种基于完全重叠区域分解的并行两步有限元算法,用于求解二维/三维时变自然对流问题。时间离散化分别采用一阶隐式欧拉公式和二阶克朗克-尼科尔森公式。这些算法基于低阶元对计算速度、压力和温度的非线性,基于同一网格上的高阶元对求解线性近似,并对算法的稳定性进行了分析。通过数值实验验证了算法的有效性和准确性。