传统光滑粒子动力学方法的适用性分析.pdfVIP

2013 年 第 58 卷 第 15 期： 1414 ~ 1421 《中国科学》杂志社 论 文 SCIENCE CHINA PRESS 传统光滑粒子动力学方法的适用性分析 * 周光正 , 葛蔚, 李静海 中国科学院过程工程研究所, 多相复杂系统国家重点实验室, 北京 100190 * 联系人, E-mail: gzzhou@ 2013-01-25 收稿, 2013-04-09 接受 国家自然科学基金和中国科学院战略性先导科技专项(XD资助 摘要 光滑粒子动力学 (smoothed particle hydrodynamics, SPH)是一种典型的无网格拉格朗 关键词 日型数值方法, 但它存在计算精度相对较低的缺陷. 为调查精度问题的物理本质, 将 SPH 粒 光滑粒子动力学 子类比于实体物理颗粒, 分析了SPH 粒子的弹性效应, 并提出了一个表征黏弹效应比的无量 流体力学 纲数. 通过二维Couette 流模拟, 考察了该无量纲数对粒子速度剖面与排布方式的影响. 研究 精度 Couette 流 结果表明, 系统的有效黏度随该无量纲数的增大而减小, 且粒子数的增加只能显著降低较小 粒子方法 无量纲数下的有效黏度. 名义黏性耗散、总耗散以及理论耗散之间的差异, 进一步验证了弹 性效应所导致的非物理耗散. 总之, 由于马赫数与该无量纲数的限制, SPH 模拟存在一个临 界雷诺数, 只有低于这个临界数, 才能基于合适的模型参数实现牛顿流体的精确模拟. 复杂流动的机理研究对很多科学和工程问题至 相对较低的计算精度、张力不稳定性以及不易处理边 关重要, 数值模拟在这方面扮演越来越重要的角色. 界条件[12]. 在数值精度的理论分析方面, SPH 的一致 光滑粒子动力学(smoothed particle hydrodynamics, 性问题已经被广泛研究[13,14]. 本质上, 由于SPH 粒子 SPH) 是一种典型的无网格粒子方法, 它在处理具有 的自适应特征, 不规则的粒子分布状态导致 SPH 不 复杂边界或变形界面的流动问题上具有独特优势. 满足变量及其导数近似的一致性条件[13,15]. 目前, 已 类似于其他粒子方法, SPH 将连续介质离散为大量具 存在多种改进传统 SPH 方法的技术方案. 再生核粒 有物理属性的数值粒子. 流体动力学守恒方程被 子法(reproducing kernel particle method, PKPM)[16]通 SPH 离散算法转化为类似于分子动力学(molecular 过构建校正函数改善 SPH 的一致性问题. 基于泰勒 [17] [18] dynamics, MD)描述的常微分方程. SPH 粒子基于这 级数展开, Chen 等人 与 Liu 等人 分别开发了针 些方程与邻域粒子相互作用, 并不断更新每一时刻 对固体力学与流体力学的修正光滑