粘弹性流体模拟的DCQ-QUICK格式-.pdfVIP

第卷 第期 工 程 数 学 学 报 年月 文章编号 粘弹性流体模拟的 格式∗ 1 1 1 2 周 文 , 解 岩 , 欧阳洁 , 李 强 西北工业大学应用数学系，西安 河南理工大学数学与信息科学学院，焦作 摘 要: 流场中对流项的离散是其数值求解的一大难点．本文基于非结构同位网格格心有 限体积法，针对流场守恒方程与Oldroyd-B 本构方程的对流项，提出了一种耦合高 阶Q-QUICK 格式的延迟修正格式．通过平面Poiseuille 流在不同W e 数下数值解与解 析解的比较，验证了该方法具有较高的精度和较好的数值稳定性．通过4:1 粘弹性收 缩流的数值模拟，揭示了不同Re、W e 数下流场中压力、应力变化及角涡生长情况， 同时也表明了该方法可有效扩大W e 数的计算范围． 关键词: 粘弹性；Q-QUICK 格式；W e 数；Oldroyd-B 分类号: (2000) 74S10; 76A10 中图分类号: O24 文献标识码: A 引引引言言言 粘弹性流体具有复杂的流变特性，实验模拟与解析求解非常困难，故数值模拟成为研 究其流动行为的一种有效手段． 粘弹性流体数值模拟方面，应用最广泛的当属有限元方法．传统有限元方法的显著 优势在于精度高，但是稳定性方面略显不足．近年来，出现了兼顾精度和稳定性的有 限元法，如SUPG 方法、EVSS 方法、DEVSS 方法等，但是计算效率普遍较低．相对而 [1] 言，有限体积法在占用存储空间、计算速度及数值稳定性方面 优于有限元方法，因而 在计算流体力学领域逐渐受到学者们的青睐．1991 年，Yoo[2] 等首次采用基于交错网格 的SIMPLER 算法对Oldroyd-B 流体在平面4:1 收缩流道中的流动进行了模拟，该工作充 分说明有限体积法模拟粘弹流动问题的可行性．2004 年，Edussuriya[3] 首次采用非结构 网格格心有限体积法推导了Oldroyd-B 流体本构方程的离散格式，并对其在平面管道和收 缩流管道中的流动做了模拟，给出可计算的最大临界W e 数为3.0，但实际上，他只给出 了矩形正交网格上的模拟结果．2011 年，Mehmet[4] 基于非结构四边形网格，采用有限体 积法结合交错网格存储技术对Oldroyd-B 流体的流动行为进行了数值模拟，但是该方法存 储麻烦，并且可计算的最大临界W e 数较小．事实上，制约粘弹性流体数值模拟的因素主 要来自两个方面：一方面，工程问题中求解区域十分复杂；另一方面，当Re 数(表征粘性 影响的相似准数) 和W e 数(表征流体弹性力和粘性力之比) 达到一定数值后，流动控制方 收稿日期 作者简介 周文年月生，男，博士 研究方向：复杂流体的高精度有限体积法 基金项目 国家重点基础研究发展计划 第期 周文，等：粘弹性流体模拟的 格式 51 程兼有椭圆型与双曲型的特点，且本构方程刚性过大，数值求解易出现数值不稳定现象． 因此针对粘弹流动模拟，寻找区域适应性强、精度高、稳定性好的数值方法非常必要． 数值模拟粘弹性流体流动问题时，对流项的离散精度在很大程度上决定了最终数 值结果的精度．本文针对对流项，发展了一种耦合高阶Q-QUICK[5] (Quasi-Quadratic Upwind Interpolation of Convective Kinematics) 格式和一阶迎风格式的高精度延迟修正 [6] 格式 ，简记为DCQ-QUICK (Deferred Co