基于弱可压及不可压光滑粒子动力学方法的封闭方腔自然对流数值模拟及算法对比.pdfVIP

物 理 学报 ActaPhys．Sin．Vo1．63，No．22(~ola) 224701 基于弱可压与不可压光滑粒子动力学方法的封闭 方腔 自然对流数值模拟及算法对比 雷娟棉 杨浩十 黄灿 (北京理工大学宇航学院，北京 100081) (2014年5月9日收到；2014年6月2日收到修改稿) 为了对 比研究弱可压光滑粒子动力学 (WCSPH)方法和不可压光滑粒子动力学 (ISPH)方法在模拟封 闭 方腔 自然对流 问题时的特性，采用粒子位移技术有效地解决了高瑞利数条件下，拉格朗日型SPH方法模拟封 闭方腔 自然对流时流体域内的粒子聚集和空穴问题，将拉格朗日型SPH方法求解封闭方腔 自然对流 问题的 最高瑞利数提高到了100；进而通过对 比瑞利数分别为 10，10，100的条件下，采用拉格 朗日型WCSPH、拉 格朗日型ISPH、欧拉型ISPH三种SPH方法模拟得到的封 闭方腔速度分布云图、温度分布云图、壁面努赛尔 特数分布曲线和平均努塞尔特数，分析了三种SPH方法在模拟封闭方腔 自然对流时的精度、稳定性和计算效 率．结果表明：在低瑞利数条件下，以上三种SPH方法都可以较好地模拟此问题，在高瑞利数条件下，欧拉型 ISPH方法的模拟结果最为精确；拉格朗日型WCSPH方法模拟所得结果比拉格朗 日型ISPH方法模拟所得结 果稍好些． 关键词：光滑粒子动力学，不可压光滑粒子动力学，粒子位移技术，自然对流 PACS：47．11．-j44．25．+f47．85．一g DOI：10．7498／aps．63．224701 SPH，WCSPH)方法与不可压光滑粒子动力学fin— 1 引 言 compressibleSPH，ISPH)方法为SPH方法的两个 重要分支．两者 的区别在于，WCSPH 方法通过 光滑粒子流体动力学fsmoothedparticlehy— 使用人工压缩方程求解压力，ISPH方法通过构造 drodynamics，SPH1方法是一种基于粒子 的无 网 压力泊松方程求解压力．近几年，国内外已有多 格数值模拟方法．它最初于 1977年 由Gingold和 名学者开展了WCSPH方法和 ISPH方法 的对 比研 Monaghan[】以及Lucy[。]提出，用于求解天体物理 究．2008年，Lee等 [16]使用 以上两种方法对溃坝和 学 问题．该方法在提出后，引起 了流体力学、固体 块体绕流问题进行了模拟，在多个较高雷诺数下考 力学领域的广泛关注、并被应用于解决各种复杂问 察了传统WCSPH方法与ISPH方法模拟所得的速 题和高度非线性 问题 [0，，包括流固耦合 问题 【]、 凝 固问题 0【]、射流 问题、液滴溅落 问题 [】、多相流问 度场、压力场 以及流体 自由表面形状等．研究表明： 题 Is，9]、热传导 问题 [1ol等．近些年，光滑粒子动力 ISPH方法在 以上几方面的精度都优于WCSPH方 学方法作为最经典的无网格方法被不断完善，多种 法．2010年，Hughes和Graham[17]比较 了以上两 修正技术被提出并应用，并且发展 了相应 的前处 种方法模拟溃坝与水波冲击竖直壁面这两个经 理 [11]与后处理技术，其 已经成为较为完善的粒子 典 自由表面流 问题时的数值结果，其在传统WC— 型无网格