粘弹性流动数值模拟方法及其在Poiseuille流中的验证.pdfVIP

16 DieandMouldTechnologyNo．42012 文章编号：1001—4934(2012)04—0016—04 粘弹性流动数值模拟方法及其在 Poiseuille流中的验证 张一丁 ，韩先洪 (上海交通大学 模具 CAD国家工程研究中心，上海 200030) 摘 要：设计和推 导了一个用 以模拟粘弹性流动 问题 的混合 型有限元方案 I—PS— DEVSS_CNBS。在基于Crank—Nicolson隐式差分 离散 的迭代型分步算法框架 内，引 入有限增量(FIC)压力稳定方法和离散粘弹应力分裂方法 (DEVSS)，分别用于克服空 间压力振荡和弹性失稳现象，本构方程采用非一致流线上风算法(SU)离散。该方案 可允许在粘弹性流动问题 中稳定有效地采用所期待的应力～速度一压力变量的等低 阶插值单元。方法的有效性将通过平面Poiseuille流进行验证。 关键词：粘弹性流体 ；有限元法；有限增量法；离散粘弹应力分裂法 中图分类号：0373 文献标识码 ：A A new numericalschemeforviscoelasticflow anditsverification in Poiseuilleproblem ZHANG Yi—ding，HAN Xian。hong Abstract：A mixedfiniteelementschemeabbreviated asI PS DEVSS CNBS scheme for — — — modelingviscoelasticflow problemsiSpresented．Thefinite incrementalcalculus (FIC) pressure stabilization process and the discrete elastic-—viscous stress--splitting method (DEVSS)areintroducedintothegeneralframeworkoftheiterativeversionofthefraction． alstepalgorithm with theuseoftheCrankNicolson-basedsplitting．Inconsistentstreamline upwindingmethod (SU)isemployedtospatiallydiscretizetheconstitutiveequationofvisco— elasticfluids．EquallOW．orderfiniteelementswhich violate theLBB compatibility condi． tionsaresuccessfullyusedintheproposedscheme．Inaddition，theproposedschemewillbe verifiedintheplanePoiseuilleflow problem quantitatiVely． Keywords：viscoelasticflow ；finiteelementmethod；FIC；DEVSS 0 引言 现象的基础 。速度 一压力形式的纯粘性牛 顿和非牛顿流在构造速度压力空间插值函数的 研究表明，在数值模型中引入粘弹性本构 时候需要满足 LBB条件 ，否则计算得到的压力 模型，是分析高聚物成型加工过程 中多种奇异 场会发生空间振荡。而描述粘弹性流动的控制 收稿 日期 ：2012-06-14 基