计算物理 ›› 2011, Vol. 28 ›› Issue (6): 853-860.

• 研究论文 • 上一篇    下一篇

基于对数构型压力稳定分步算法的黏弹性流动模拟

王宣平, 李锡夔   

  1. 大连理工大学工业装备结构分析国家重点实验室, 辽宁 大连 116023
  • 收稿日期:2010-11-24 修回日期:2011-03-31 出版日期:2011-11-25 发布日期:2011-11-25
  • 通讯作者: 李锡夔,E-mail:xikuifi@dlut.edu.cn
  • 作者简介:王宣平(1979-),男,辽宁大连,博士生。从事粘弹性流动数值模拟研究.E-mail:wangxp@mail.dlut.edu.cn
  • 基金资助:
    国家自然科学基金(11072046,90715011,10590354,10672033);国家973(2010CB731502)资助项目;大连理工大学高性能计算平台提供支持

Log-conformation-based Pressure-stabilized Fractional Step Algorithm for Viscoelastic Flows

WANG Xuanping, LI Xikui   

  1. State Key Laboratory of Structural Analysis of Industrial Equipment, Dalian University of Technology, Dalian 116023, China
  • Received:2010-11-24 Revised:2011-03-31 Online:2011-11-25 Published:2011-11-25

摘要: 利用近年来发展的对数构型方法,通过有限增量过程引入稳定化机制,提出对数构型压力稳定分步算法,显著降低每个迭代步的计算量,避免对构型张量梯度的差分近似.引入压力稳定化机制,不仅避开了LBB条件的限制、允许待求变量采用等低阶插值,且显著地增强了算法稳定性.黏弹性流体圆柱绕流问题的数值模拟表明,对数构型压力稳定分步算法具有优良的稳定性和精度.

关键词: 粘弹性, 对数构型张量, 分步算法, 压力稳定, 有限元

Abstract: Based on log-conformation method developed recently,a log-conformation-based fractional step algorithm with stabilization mechanism is proposed by virtue of the finite increment calculus(FIC)process.The algorithm remarkably reduces computational cost within each iteration step compared with existed coupled solution schemes.Meanwhile it bypasses finite difference approximation for gradient of conformation tensor.As FIC process is introduced as pressure stabilization mechanism,restrictions on u-p interpolations imposed by LBB compatibility conditions are circumvented and equal-lower-order interpdlations for u-p are allowed to apply.Numerical simulations for cylinder flow show excellent stabilization and accuracy of the proposed algorithm.

Key words: Viscoelasticity, log-conformation tensor, fractional step algorithm, pressure stabilization, finite element method

中图分类号: