黏性液滴变形过程的核梯度修正光滑粒子动力学模拟鄢-物理学报.pdfVIP

物理学报摇 粤糟贼葬孕澡赠泽援杂蚤灶援摇 灾燥造援远园袁晕燥援缘 渊圆园员员冤摇 园缘源苑园员 黏性液滴变形过程的核梯度修正 光滑粒子动力学模拟鄢 覮 蒋摇 涛摇 欧阳洁 摇 赵晓凯摇 任金莲 渊西北工业大学理学院应用数学系袁西安摇 苑员园员圆怨冤 渊圆园员园年苑月怨 日收到曰圆园员园年愿月圆圆 日收到修改稿冤 摇 摇 本文提出了一种核梯度改进光滑粒子动力学渊运郧悦鄄杂孕匀冤方法袁模拟了黏性液滴形变自由表面问题援首先袁通过 模拟等温黏性液滴拉伸和旋转变形袁验证了运郧悦鄄杂孕匀法较杂孕匀法具有较高精度和更好稳定性袁且能很好地保持总 角动量守恒援其次袁基于非等温增葬灶凿藻则宰葬葬造泽模型对平衡态圆形液滴的形成过程进行数值研究袁观察到小幅度振荡 现象袁并给出了一种新的克服张力不稳定性的方法和一种适合 运郧悦鄄杂孕匀方法的新的表面张力处理技术援最后袁研 究了增葬灶凿藻则宰葬葬造泽液滴的周期性振荡现象袁讨论了初始椭圆形液滴长短半轴比例尧雷诺数对液滴振幅及振荡周期 的影响援 关键词院光滑粒子动力学袁黏性液滴袁增葬灶凿藻则宰葬葬造泽模型袁表面张力 孕粤悦杂院源苑援员员援原躁袁源源援员园援垣蚤袁源苑援愿缘援原早 总角动量守恒援近年来袁针对 杂孕匀 的缺点已提出了 员郾 引 言 很多改进方法咱员苑要圆园暂 袁这些改进格式大多着眼于恢复 杂孕匀方法中核近似的相容性袁以改善其计算精度袁 咱员袁圆暂 咱猿袁源暂 自由表面流动 或大变形 现象广泛存在 同时忽略了粒子之间的相互作用援事实上袁利用传 于日常生活或工业生产之中袁目前人们已经发展了 统杂孕匀模拟流动问题时袁通常将流体控制方程中速 多种数值方法来追踪自由表面援 如质点粒子单元 度尧应力等物理量的导数形式近似成以核梯度形式 咱缘暂 咱远袁苑暂 咱愿要员园暂 法 和蕴葬贼贼蚤糟藻鄄月燥造贼扎皂葬灶灶 等无网格类方法 袁 的粒子离散形式援 因此袁控制方程离散格式中核梯 有限差分尧有限体积等网格类方法援光滑粒子动力 度的精度袁以及在此格式中是否考虑了粒子间的相 学渊泽皂燥燥贼澡藻凿责葬则贼蚤糟造藻澡赠凿则燥凿赠灶葬皂蚤糟泽袁杂孕匀冤方法是一 互作用将直接影响到整个数值模拟的精度尧稳定性 种纯粹的粒子方法袁该方法最初提出是为求解三维 等问题援 无边界空间的行星运动问题咱员员暂援在 员怨怨源年开始被 本文在黏性变形液滴的数值模拟过程中袁为提 应用于自由表面流问题咱员圆暂 袁并被广泛应用于黏性流 高传统杂孕匀离散格式的精度与稳定性袁基于梯度离 体咱员猿袁员源暂 袁多相流咱员缘袁员远暂等问题援与上述方法相比袁杂孕匀 散和栽葬赠造燥则级数展开思想袁对核梯度形式进行修正袁 方法计算空间导数时不需要使用任何网格袁自由面 并考虑粒子间相互作用袁提出了核梯度修正 杂孕匀 模拟过程中不需要额外追踪技术渊如灾韵云渊增燥造怎皂藻鄄 渊噪藻则灶藻造早则葬凿蚤藻灶贼糟燥则则藻糟贼藻凿杂孕匀袁运郧悦鄄杂孕匀冤方法援与 燥枣鄄枣造怎蚤凿冤方法冤曰能避免对流占优带来的数值困难问 传统杂孕匀法相比袁此方法具有较高的精度和更好的 题曰程序实现容易袁特别是三维问题援 稳定性袁并且能很好保持总角动量守恒援文中先是 传统杂孕匀方法虽然有很多优点袁但它也存在一 模拟了等温情况下黏性液滴拉伸和旋转变形过程袁 些缺点袁其缺点主要表现在精度不高和数值稳定性