非等温非牛顿黏性流体流动问题的修正光滑粒子动力学方法模拟精选.pdfVIP

物 理 学报 ActaPhys．Sin．Vo1．63，No．21(2014) 210203 非等温非牛顿黏性流体流动问题的修正光滑粒子 动力学方法模拟冰 蒋涛十 任金莲 徐磊 陆林广 (扬州大学数学科学学院，扬州大学水利与能源动力工程学院，扬州 225002) (2014年5月7日收到；2014年6月l1日收到修改稿) 为准确、有效地模拟非等温非牛顿黏性流体的流动 问题，本文基于一种不含核导数计算的核梯度修正格 式和不可压缩条件给出了一种改进光滑粒子动力学fSPH)离散格式，它较传统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 引 言 应用到聚合物流动过程中大变形问题的模拟． 近年来，无网格方法 [9-12】以其独特优势在非 非等温非牛顿黏性流体f主要指黏性聚合物熔 等温传热 问题的研究中得到了广泛应用，但大多 体1流动过程 【，】普遍存在于工业上聚合物材料的 研究都集 中在稳态传热 问题，对瞬态传热问题的 加工过程和 日常生活中，此过程包括复杂的热量传 研究还处于探索阶段，特别是非等温聚合物流体 递与转化，对熔体黏度及流动有重要影响 3[，41．因 流动问题．于是，本文采用光滑粒子动力学 (SPH) 此，非等温非牛顿黏性流体流动的研究对化工工 方法 [12—15]研究非等温热幂率模型和热Cross模 程、注塑成型及环境工程等产业具有重要意义． 型 [16]为主的聚合物流体流动 问题．SPH方法是一 非等温黏性聚合物流动的研究涉及瞬态温度 种纯无网格粒子方法，它完全摈弃了网格，避免了 控制方程 的求解、流体控制方程 的求解和 自由界面 网格类方法中存在的网格扭曲、网格重新划分等问 的捕捉．目前，已有很多学者采用网格类方法对 题．然而，目前SPH方法用于非等温流动的情况 比 非等温聚合物流动进行研究 2[-5]．然而，网格类方 较少，仅有Jiang等 [7J研究了非等温Couette流 问 法在模拟聚合物 自由面流动问题时存在一些不足， 题、Cleary等 [18]研究了熔融铝液的模具铸造 问题， 如有限差分法、有限体积法等网格类方法在实现过 本文作者L19J运用SPH方法研究了瞬态热传导问 程中需要网格重构，且需要借助界面追踪技术f如 题．然而上述非等温离散模型不适合模拟非牛顿黏 level—set方法、流体体积方法等 0【，】1来捕捉流体前 性流体的非等温流动过程．此外，传统SPH方法数 沿界面．此外，网格类方法难 以模拟复杂的大变形 值精度低、稳定性差，尤其是在粒子分布不均匀或 中国博士后科学基金 (批准号：2014M550310)、国家自然科学基金青年科学基金(批准号 513