多面体网格和多边形网格在计算流体力学的应用近十年才兴起,其能够克服传统网格的一些缺点,包括更多的临近单元使得梯度计算和局部流动状况更准确,并且同样的计算域对网格数目要求更少,从而减少计算量。相对于四面体网格,多面体网格减少了网格变形对于计算精度的影响。由于目前能够直接生成和优化处理多面体网格的工具较少,相关文档细节资料也比较缺乏,造成多面体网格应用的不便。本文主要根据近些年来计算图形学领域的成果,综合前人的对多面体网格的处理方案,针对多面体网格的空间多边形面单元的翘曲问题提出一套抑制算法。该算法基于构建在空间多边形上的离散拉普拉斯算子,组建体现网格全局翘曲程度的二次泛函能量,并根据拟牛顿法L-BFGS(Limited-memory Broyden-Fletcher-Goldfarb-Shanno)优化算法求解泛函能量的极值。该算法在Code_Saturne开源通用计算流体软件平台上实现,并通过U形管算例验证了该算法的性能,以及对于实际流体计算结果的影响。结果表明本文提出的模型针对面单元含有任意边数的非结构网格均具有良好的网格面抑制翘曲的效果。