论文部分内容阅读
随着信息化社会科技的飞速发展,人们对可视化的需求也越来越高:影片中洪水、海浪、漩涡等精彩的视觉特效;宇宙飞船在太空中的飞行过程;大气运动形成的雨、雪、风等自然场景的仿真;自然灾害一旦发生时对建筑物和人类的冲击力和伤害程度的真实效果模拟;现实物理场景雨滴、喷泉、雪飘、烟雾等的计算机可视化效果的模拟;都是计算流体力学的重要研究课题。所谓计算流体力学就是用计算机和离散化的数值方法对一些物理现象包括流体流动、热传导、物质扩散等现象进行的数值分析和模拟的一门新兴学科。这门学科在工程及科学领域中已经得到了广泛的应用。这种数值模拟方法相对于流体实验方法而言有许多优势:数值模拟能够研究流体实验难以控制的系统;能够超出人类行为的极限,研究危险和人类不可及条件下的系统;可以大幅度的减少设计所需的时间和成本;而且操作自由灵活,可以无限量的提供研究结果的细节。要了解计算流体力学的研究方法,首先要了解流体力学的基本方程的性质和分类。流体力学的基本方程是由C.-L.-M.-H纳维和G.G.斯托克斯等人在19世纪上半叶建立的Navier-Stokes方程,简称N-S方程。求解该方程的数值方法可分为网格方法和无网格方法两类。本文主要介绍无网格方法中的一种——移动粒子半隐式(MPS)法。首先,本文详细介绍了移动粒子半隐式(MPS)法的理论模型和数值求解过程,该方法是一种基于Lagrange观点的数值方法。流体是由有限数量的移动粒子表示,这些粒子同时还携带着材料性质,可以在外力和内力的相互作用下运动。控制方程是在粒子间相互作用模型所表示的梯度模型、拉普拉斯模型和自由表面模型的基础上离散而得的。该方法不需要建立网格,因此不受网格扭曲的影响,适合模拟带有大变形的水力学问题。MPS法不同于其他一些粒子法,如SPH法是使用弱可压缩模型来近似模拟不可压缩流体流动性问题。MPS法最主要的特点就是引入了粒子数密度的概念,为流体的不卡压缩性提供依据,使MPS法成为模拟真正不可压流体问题的方法。同时MPS法引入了半隐式的算法,有效地保持了流体的不可压缩性。而且该方法属于完全Lagrange描述,不存在因对流而形成的数值扩散。之后本文还详细介绍了邻域粒子的搜索方法,并与传统的搜索方法进行了比较。其次,本文在能量守恒方面对MPS法做了详细的分析,如线动量、角动量的守恒性问题。针对这个问题,本文还简要介绍了一种MPS法的改进方法(CMPS法),改进后的CMPS法使整个计算过程能够更好的满足线动量和角动量的守恒。此外,本文通过使用MPS法做了一些数值模拟实验。模拟了单水滴的三维可视化效果和二维溃坝现象,模拟效果真实感较强,有效地证明了方法的高仿真性。