水利机械计算水动力学------提交报告.docx

一、求解N-S方程的数值方法求解 N-S 方程的数值方法可分为有网格和无网格两类。有网格方法的求解思路是将连续空间划分为离散网格 (或单元)，并在网格点 (或单元) 上离散控制方程，从而寻求控制方程的数值解；无网格方法则将着眼点回归到流体运动的基本单元质点上，通过追踪质点运动性质的变化过程来获得流体解答。网格数值模拟方法发展较早，求解方法较成熟，应用也较广；而无网格方法则刚刚起步，具有很大的发展空间。下面分别就两类中较有代表性的方法予以介绍。1.1 网格数值模拟方法数值计算网格可分为不随流体运动而改变形状的欧拉网格和随流体运动改变形状的拉格朗日网格两类。欧拉网格划分较简单，但在用于具有自由表面流体运动数值计算中需要对自由表面位置进行适时追踪；拉格朗日网格虽然无需进行自由表面追踪，但需要在计算过程中随时调整计算网格，以适应流体自由表面形状，额外增加了计算量，并且在流体变形剧烈处还会出现网格失效现象。有限差分法 (Finite Difference Method, FDM)和边界元法 (BoundaryElement Method, BEM) 都是欧拉网格方法中常用的方法。在拉格朗日网格法中，有限元法和边界元法可直接应用，而有限差分法一般需对坐标做适体变换。用欧拉网格法求解含有自由表面的水波动力问题难点在于追踪自由表面。目前比较有代表性的自由表面追踪方法有PIC (Particle-I