水动力学数值模拟MPS方法的改进及应用.pdfVIP

中文摘要 阐述了目前水动力学研究中数值模拟方法的内容，按照拉格朗日观点和欧拉 观点将数值方法划分为网格方法和无网格方法，选取各自具有代表性的数值方法 进行介绍。其中MPS(movingparticlesemi．implicit)方法是一种新兴的无网格方 法，在处理边界复杂、自由表面剧烈变形的流体问题上具有很大优势，但其运行 效率低，限制了其实际应用，急需消除瓶颈。 在MPS方法中，引入预定候选粒子集的概念，邻接粒子在预定候选粒子集 中搜索生成。以数值算例说明了候选粒子集的引入可以显著节省生成邻接粒子集 的时间。从理论上讨论分析了更新候选粒子集的时间步数k与生成邻接粒子集平 均耗费时间的关系，指出了存在k使得平均耗费时间最小。 在求解压力泊松方程的程序代码中，采用一维数组只记录方程系数矩阵中的 非零元素，降低了程序运行占用的内存量，简化了矩阵运算。采用稳定性好，收 敛速度快的Bi．CGSTAB算法求解方程，提高了求解速度。将改进前后的方程求 解过程耗时进行对比，表明该部分占用的机时显著减少。 根据实验资料设计数值孤立波水槽，以改进的MPS模型模拟孤立波在水槽 中的传播过程，并首次将MPS模拟的自由水面随时间变化过程与实测资料对比， 结果吻合较好。 发现MPS模型模拟所得的压力随时间变化存在不规则的非物理振荡，针对 该问题，建议深入研究粒子间的相互作用方式，选取更合适的核子函数，或者引 入ALE模型作为辅助运算，以使压力变化更光滑。 关键词：MPS方法，预定候选粒子集，压力泊松方程，Bi．CGS仉心算法，孤j 波 ABSTRACT Thenumericalsimulationmethodson are in are reviewed hydrodynamics detail．They intotwo methodand the classifiedcategories-gridgridless to method，corresponding andEulerian of Lagarangian the methods approach,respectively,andrepresentative each areintroduced．MPSmethodisanew of in category approachgreatadvantage witll of andsituationswhereviolent dealingproblemscomplexboundaryshape deformationoffreesurfaceoccurs．However，itismoreorlessdeficientincalculation duetothe amountof the ofMPSmethod speed large operations．Thusefficiency must be for improvedpracticalapplication． The ofcandidatelistis into to concept particle introducedMFSme也odgenerate