基于弱可压与不可压光滑粒子动力学方法的封闭方腔自然对流数值.PDF

物理学报 Acta Phys. Sin. Vol. 63, No. 22 (2014) 224701 基于弱可压与不可压光滑粒子动力学方法的封闭 方腔自然对流数值模拟及算法对比 雷娟棉 杨浩 黄灿 (北京理工大学宇航学院, 北京 100081) ( 2014 年5 月9 日收到; 2014 年6 月2 日收到修改稿) 为了对比研究弱可压光滑粒子动力学(WCSPH) 方法和不可压光滑粒子动力学(ISPH) 方法在模拟封闭 方腔自然对流问题时的特性, 采用粒子位移技术有效地解决了高瑞利数条件下, 拉格朗日型SPH 方法模拟封 闭方腔自然对流时流体域内的粒子聚集和空穴问题, 将拉格朗日型SPH 方法求解封闭方腔自然对流问题的 最高瑞利数提高到了10 ; 进而通过对比瑞利数分别为10 , 10 , 10 的条件下, 采用拉格朗日型WCSPH、拉 格朗日型ISPH、欧拉型ISPH 三种SPH 方法模拟得到的封闭方腔速度分布云图、温度分布云图、壁面努赛尔 特数分布曲线和平均努塞尔特数, 分析了三种SPH 方法在模拟封闭方腔自然对流时的精度、稳定性和计算效 率. 结果表明: 在低瑞利数条件下, 以上三种SPH 方法都可以较好地模拟此问题, 在高瑞利数条件下, 欧拉型 ISPH 方法的模拟结果最为精确; 拉格朗日型WCSPH 方法模拟所得结果比拉格朗日型ISPH 方法模拟所得结 果稍好些. 关键词: 光滑粒子动力学, 不可压光滑粒子动力学, 粒子位移技术, 自然对流 PACS: 47.11.–j, 44.25.+f, 47.85.–g DOI: 10.7498/aps.63.224701 SPH, WCSPH) 方法与不可压光滑粒子动力学(in- 1 引 言 compressible SPH, ISPH) 方法为SPH 方法的两个 重要分支. 两者的区别在于, WCSPH 方法通过 光滑粒子流体动力学(smoothed particle hy- 使用人工压缩方程求解压力, ISPH 方法通过构造 drodynamics, SPH) 方法是一种基于粒子的无网 压力泊松方程求解压力. 近几年, 国内外已有多 格数值模拟方法. 它最初于1977 年由Gingold 和 名学者开展了WCSPH 方法和ISPH 方法的对比研 Monaghan1 以及Lucy2 提出, 用于求解天体物理 16 究. 2008 年, Lee 等 使用以上两种方法对溃坝和 学问题. 该方法在提出后, 引起了流体力学、固体 块体绕流问题进行了模拟, 在多个较高雷诺数下考 力学领域的广泛关注, 并被应用于解决各种复杂问 3,4 5 察了传统WCSPH 方法与ISPH 方法模拟所得的速 题和高度非线性问题 , 包括流固耦合问题 、 凝固问题6 7 度场、压力场以及流体自由表面形状等. 研究表明: 、射流问题、液滴溅落问题 、多相流问 8,9