SPH法理论基础详解.doc

摘要 1 第一章 绪论 2 1.1 基于网格的方法 2 1.1.1 拉格朗日网格 3 1.1.2 欧拉网格 3 1.1.3 基于网格的数值方法的局限性 3 1.2 无网格方法 4 1.3光滑粒子动力学方法 4 第二章 SPH方法在流体动力学问题中的应用 5 2.1 光滑粒子动力学原理 5 2.2 光滑粒子动力学基本方程 5 2.2.1 函数的积分表达 5 2.2.2 函数的粒子表达 7 2.2.3光滑函数 7 2.3 拉格朗日型的Navier-Stokes 方程 7 2.4 Navier-Stokes方程的SPH表达式 8 第三章 SPH的计算实施方法 9 3.1 粒子的密度近似 9 3.2 核函数 10 3.3 状态方程 11 3.4 人工粘度 12 3.5 边界处理 13 3.6 时间积分 13 摘要 计算机数值仿真逐渐成为解决现代工程和科学问题的一条重要途径。数值仿真能为理论提供测试和检验，有助于对复杂的物理问题加深认识，甚至还能帮助解释和发现新现象。 基于网格的数值方法虽然已经有广泛的应用，但是在很多方面仍存在不足之处，比如在计算流体动力学的大变形、运动物质交界面、自由表面等问题时，由于网格产生畸变导致计算误差过大或无法进行，从而使其在许多问题的应用上受到限制。 近年来，无网格法倍受关注，这种方法在许多应用中都优于传统的基于网格的有限元法、有限差分法以及有限体积法等数值方法。本文主要研究新一代无网格方法---光滑粒子流体动力学方法（SPH）。 第一章 绪论 1.1 基于网格的方法 数值计算方法通常可以分为两种：基于网格的方法和无网格方法。通常对于物理控制方程的描述有两种基本方法：欧拉描述法和拉格朗日描述法。欧拉描述法是对空间的描述方法，其典型代表是有限差分法；拉格朗日描述法是对物质点的描述方法，其典型代表是有限元法。欧拉描述和拉格朗日描述对应着两种不同的区域离散化网格：欧拉网格和拉格朗日网格。针对不同类型的问题，这两种网格在数值方法中都得到广泛应用。 1.1.1 拉格朗日网格 拉格朗日网格，基于拉格朗日网格的数值方法在整个计算过程中网格是固定附着于物质上的，网格会随着物质的运动而运动，所以在物质点上的所有场变量的整个时间历程都可以很容易地追踪，常见的方法如有限元法等。 基于拉格朗日网格方法的优点是：由于在相关的偏微分方程里不存在迁移项，所以程序在方案设计上会变得相对简单而且运行较快；由于只需在问题域内布置网格，问题域外不需要布置，所以计算效率很高；不规则或者复杂的几何形状可以用不规则的网格来处理。由于具有以上这些优点，拉格朗日方法得到广泛的应用，并且能成功地求解计算固体力学问题。 然后，基于拉格朗日网格的方法难以应用于具有极大网格变形的情况，因为其公式的形式是以网格为基础的，当网格变形太大的时候，公式的精度和求解都会受到很大的影响。另外，由于时间步长是由最小单元尺寸所控制，若网格太小就会影响计算的效率，甚至会导致计算失败。 1.1.2 欧拉网格 欧拉网格，相对于拉格朗日网格，欧拉网格刚好相反，它是固定在模拟对象所处的空间上，模拟对象在固定网格单元上运动。因此，在物质流过网格时，所有网格节点和网格单元依然固定在空间上而且不会随时间的改变而改变。通过模拟质量、动量和能量经过网格单元边界的通量，可以计算质量、速度和能量等等的分布。在整个计算过程中网格单元的形状和体积都保持不变。 由于欧拉网格在时间和空间上都是固定的，物体的大变形不会引起网格本身的任何变化，所以再以物质流动为主体的计算流体力学领域里，较为广泛的应用欧拉法。 但是，欧拉法依然有很多缺点。由于欧拉法不能用固定网格来追踪物质的运动，所以很难分析物质上固定点的场变量的变化情况而只能得到空间上固定的欧拉网格的场变量的变化情况；由于欧拉法追踪的是流过网格单元边界的质量、动量和能量的通量，所以自由表面、变形边界和运动物质交界面的位置就很难精确确定；由于欧拉法需要在整个计算区域上都覆盖网格，所以有时为了提高计算效率而使用较粗糙的网格，这样就会降低离散化区域的求解精度。 1.1.3 基于网格的数值方法的局限性 传统的基于网格的数值方法如有限差分法和有限元法已经广泛地应用在计算流体力学和计算固体力学的各个领域中，是现在进行区域离散化和数值离散化模拟的主要方法。但是基于网格的数值方法在很多方面仍存在不足之处，比如在计算流体动力学中的大变形、运动物质交界面、自由表面等问题时，由于网格产生畸变导致计算误差过大或无法进行，使其在许多问题的应用上受到限制。 在基于网格的数值方法中，数值模拟的先决条件就是在问题域生成网格。对于基于欧拉网格的方法，在固定的欧拉网格上要精确确定自由表面、变形边界、运动交界面和不均匀物质之间的位置是一项非常困难的工作。并且欧拉方法也不适用于研究如粒子流动这类问题。