非等温非牛顿黏性流体流动问题的修正光滑粒子动力学方法-物理学报.pdfVIP

物理学报 Acta Phys. Sin. Vol. 63, No. 21 (2014) 210203 非等温非牛顿黏性流体流动问题的修正光滑粒子 动力学方法模拟∗ 蒋涛 任金莲 徐磊 陆林广 (扬州大学数学科学学院, 扬州大学水利与能源动力工程学院, 扬州 225002) ( 2014 年5 月7 日收到; 2014 年6 月11 日收到修改稿) 为准确、有效地模拟非等温非牛顿黏性流体的流动问题, 本文基于一种不含核导数计算的核梯度修正格 式和不可压缩条件给出了一种改进光滑粒子动力学(SPH) 离散格式, 它较传统SPH 离散格式具有较高精度 和较好稳定性. 同时, 为准确地描述温度场的演化过程, 建立了非牛顿黏性的SPH 温度离散模型. 通过对等 温Poiseuille 流、喷射流和非等温Couette 流、4:1 收缩流进行模拟, 并与其他数值结果作对比, 分别验证了改进 SPH 方法模拟非牛顿黏性流动问题的可靠性和提出的SPH 温度离散模型求解非等温流动问题的有效性和准 确性. 随后, 运用改进SPH 方法结合SPH 温度离散模型对环形腔和C 形腔内非等温非牛顿黏性流体的充模 过程进行了试探性模拟研究, 分析了数值模拟的收敛性, 讨论了不同位置处热流参数对温度和流动的影响. 关键词: 非等温, 黏性流体, 光滑粒子动力学, 热黏度模型 PACS: 02.70.–c, 44.10.+i, 47.11.–j, 47.85.md DOI: 10.7498/aps.63.210203 8 现象, 如喷射弯曲、飞溅等现象 , 使其不易被推广 1 引 言 应用到聚合物流动过程中大变形问题的模拟. 近年来, 无网格方法912 以其独特优势在非 非等温非牛顿黏性流体(主要指黏性聚合物熔 等温传热问题的研究中得到了广泛应用, 但大多 体) 流动过程1,2 普遍存在于工业上聚合物材料的 研究都集中在稳态传热问题, 对瞬态传热问题的 加工过程和日常生活中, 此过程包括复杂的热量传 研究还处于探索阶段, 特别是非等温聚合物流体 递与转化, 对熔体黏度及流动有重要影响3,4 . 因 流动问题. 于是, 本文采用光滑粒子动力学(SPH) 此, 非等温非牛顿黏性流体流动的研究对化工工 方法1215 研究非等温热幂率模型和热Cross 模 程、注塑成型及环境工程等产业具有重要意义. 16 型 为主的聚合物流体流动问题. SPH 方法是一 非等温黏性聚合物流动的研究涉及瞬态温度 种纯无网格粒子方法, 它完全摈弃了网格, 避免了 控制方程的求解、流体控制方程的求解和自由界面 网格类方法中存在的网格扭曲、网格重新划分等问 的捕捉. 目前, 已有很多学者采用网格类方法对 题. 然而, 目前SPH 方法用于非等温流动的情况比 非等温聚合物流动进行研究25 . 然而, 网格类方 17 较少, 仅有Jiang 等 研究了非等温Couette 流问 法在模拟聚合物自由面流动问题时存在一些不足, 18 题、Cle