不可压缩N-S方程的稳定分步算法及耦合离散.pdfVIP

第39卷第6期 力 学 学 报 V01．39．No．6 ChineseJournalofTheoreticMand Mechanics NOV．．2007 2007年11月 Applied 不可压缩N—S方程的稳定分步算法及耦合离散u 段庆林李锡夔2’ (大连理工大学工业装备结构分析国家重点实验室，大连116023) increment 摘要在有限增量微积分(finite calculus，FIC)的理论框架下，通过引入一个附加变量，发展了压 力稳定型分步算法，有效改善了经典分步算法的压力稳定性，同时还避免了标准FIC方法中存在的空间高阶 导数的计算．为保证数值方法同时具有较快的计算速度和较好的健壮性，发展了有限元与无网格的耦合空间离 散方法．该方案可在网格发生扭曲的区域采用无网格法空间离散以保证求解的精度和稳定性。而在网格质量较 好的区域以厦本质边界上保留使用有限元法空间离散以提高计算救率和便于施加本质边界蕞件．方腔流考题的 数值模拟结果突出地显示了所发展的压力稳定型分步算法比经典分步算法具有更好的压力稳定性，能够有效消 除速度一压力插值空间违反LBB条件而导致的压力场的虚假数值振荡．平面Poisseuille流动和一个典型型腔 充填过程的数值模拟结果，表明了发展的耦台离散方案相对于单一的有限元法和单一的无嘲格法在综合考虑计 算效率和算法健壮性方面的突出优点． 关键词分步算法，不可压缩N-S方程，LBB条件，压力稳定性，无网格法，FIC 中图分类号；0357．1文献标识码： A 文章编号：0459-1879(2007)06-0749-11 引 言 稳定化形式，然后对其采用标准Galerkin空间离散 即可消除压力场的虚假振荡．值得注意的是，广泛 应用于不可压缩流模拟的经典的压力非增量型分步 不可压缩N-S方程的准确求解历来是计算流体 算法能够允许速度、压力的同低阶插值，但其压力 力学的研究热点之一．其一大困难在于质量守恒方 稳定性依赖于时间步长的选择；另一方面，被认为具 程中不包含压力变量，因而整个方程组是混合形式 有更高精度的压力增量型分步算法却不允许速度、 的．这种混合特性要求速度，压力的插值空间满足所 压力的同低阶插值㈡Codk,m等【7J在仔细分析非 谓的LBB条件，否则，标准的Galerkin空同离散会导 增量型分步算法压力稳定机制的基础之上，通过引 致压力场的虚假振荡．另一方面，尽管少数几种构过 入一个附加变量(压力梯度投影)来改善经典分步算 较为复杂的插值模式(如有限元法中的T6P3，Q2PI 法的压力稳定性，发展了压力梯度投影法(pressure 单元等)能够满足这一相容条件，但是，无论有限元法 gradient projection，PGP)． 或无网格法中，从应用的角度看更为方便的速度、压 力的同低阶插值模式并不能满足该条件(k“．为解决 至于空间离散化方法的选择，目前主要有无网 这一问题，Hughes等川首先放弃标准Galerkin空间格法和有限元法，它们各有优缺点．有