基于SPH方法的瞬态粘弹性流体的数值模拟.pdfVIP

第 27卷 第 5期 计 算 物 理 CHINESE JOURNALOFCOMPUTATIONALPHYSICS Vo1．27．NO．5 2010年 9月 Sep．，2010 文章编号 ：1001—246X(2010)05 679—06 基于 SPH方法的瞬态粘弹性流体的数值模拟 杨 波 ， 欧阳洁 (西北工业大学应用数学系 ，陕西 西安 710072) 摘 要 ：运用 SPH (SmoothedParticleHydrodynamics)方法模拟基于 Oldroyd—B模型的平面突然起动 Couette流 。 通过数值解与解析解 的比较 ，验证 SPH方法模拟瞬态粘弹性流动的准确性 ；且对基于 Oldroyd—B模型的方腔驱动 流进行 SPH模拟．采用一种新 的固壁边界处理方法 ，有效地防止 了粒子穿透 ，提高数值计算 的准确性 ．用数值 算例验证 SPH方法对粘弹性流体模拟的有效性和稳定性 ． 关键词 ：光滑粒子动力学 ；粘弹性 ；Oldroyd—B；瞬态 ；Couette流 ；方腔驱动流 中图分类号 ：0373 文献标识码 ：A U 引 舌 近年来 ，随着工业生产中成型技术 的 日益发展 ，粘弹性流体流变特性的研究已成为 国内外学者关注的热 点．目前 ，模拟粘弹性流体 的数值格 式大多基于 网格 ，如有 限差分法 、有限元法等．光滑粒子流体动力学 (SmoothedParticleHydrodynamics，简称 为 SPH)方 法 是一种 基 于粒子 的纯 Lagrange无 网格方 法，最早 由 Lucy”和 Gingold 等于 1977年提出．与其它数值方法相比，SPH方法完全不需要 网格，不仅免去了生成网 格的麻烦 ，而且避免 了通常 Lagrange方法 中的网格缠结、网格扭 曲以及网格重构等问题 。 ． SPH方法用于粘弹性流体 的模拟是近年来探索研究 ．目前对粘性和粘弹性模型的 Couette流 ，均已 有了很多基于网格的数值模拟结果．然而 ，对于方腔驱动流 ，大部分研究集 中在粘性流动问题 ，粘弹性流动方 面 的研究较少．至今 尚未见到基于 Oldroyd．B模型对粘弹性 Couette流和方腔驱动流 问题进行 SPH方法 研究的报道． 本文将 SPH方法首次应用到基于 Oldroyd—B本构模型的突然起动平面 Couette流和方腔驱动流的数值 模拟中，验证了SPH方法模拟粘弹性流体 的准确性和可行性 ，并对问题中流变参数对流动的影响进行了讨 论 ． 1 控制方程 在 Lagrange坐标系下 ，等温粘弹性流体的控制方程为 + v．1，：0， (1) P Dt ’ P =一 V叩 + V·T+F， (2) 其 中P， ，P和 分别表示流体的密度 、速度、压力和偏应力张量 ，F表示流体受到的外力 ，D／Dt表示物质 导数． 1．1 状态方程 本文采用的状态方程为 P(P)=c2p， (3) 收稿 日期：2009—06—09；修回 日期 ：2010一O1—15 基金项 目：国家 自然科学基金重大 项 目(批 准号 ，国家 自然 科 学基金 (批 准号 和 国家重点基础 研究发展计 划 (973)项 目(批准号 ：2005CB321704)资助项 目 作者简介 ：杨波 (1985一)，女 ，江苏泰州 ，硕士生 ，从事复杂流体数值模拟及 自由面流动 问题 的无网格粒子方法研究． 680