论文部分内容阅读
不可压自由表面流问题主要是指与大气存在交界面的水流问题,在该交界面上,物理量会发生间断或不连续的变化。在船舶与海洋工程、水利工程、机械工程、石油化工及土木工程中,均存在大量的此类问题,因此对不可压自由表面流问题进行研究具有极为重要的工程意义。SPH法是一种完全无网格、纯拉格朗日的粒子法,其优点在于对流项直接通过粒子运动来模拟,完全消除自由界面上的数值发散,保证了自由液面追踪的准确性,且完全不需要网格,免去了网格生成的麻烦,更避免了网格扭曲与网格重构问题,能够十分方便地模拟具有自由液面大变形的流动问题。现有的研究工作已经使SPH法在自由表面流中的应用取得了较大的进展,但是仍然存在许多不足之处,本文即做进一步的研究,以推进该方法在不可压自由表面流中的应用。在基本理论方面,本文详细论述了SPH法基本思想、函数及其一、二阶导数的SPH积分表示方法,以及函数的SPH法粒子近似,建立了SPH法基本方程;对核函数、支持域与影响域等基本概念进行了分析,并阐述了MSPH法的基本原理和计算公式。然后在质量、动量、能量守恒的基础上得到笛卡儿坐标系、柱坐标系及球坐标系下的Lagrange型流体运动N-S方程,通过对该方程进行SPH空间离散,全面推导了各种坐标系下适用于广义流体动力学的SPH控制方程。在二维流体动力SPH法控制方程的基础上,研究了非冲击性不可压自由表面流SPH法数值模拟相关技术,其中首次详细讨论了初始布点方法与粒子间相互作用对称性问题,并提出了将SPH边界处理方法与一般流体流动数值模拟边界处理方法相结合的思想。通过液体受迫晃荡算例验证了所研究的数值技术的有效性。针对弱冲击性自由表面流问题的特点,根据MSPH法建立了修正后的二维SPH法流体控制方程,并对相关SPH法数值技术做了研究。通过溃坝算例对其合理性做了验证。对具有强冲击特性的自由表面流问题的SPH法数值模拟的实现做了深入研究。首先建立了引入黎曼解法的SPH流体控制方程以考虑强冲击性流动数值解的间断性,然后对SPH法相关数值技术进行了研究,其中提出了新的自由表面粒子确定方法,并首次在自由表面流SPH法模拟中考虑了张力不稳定性问题及人工应力项。最后通过液滴冲击液面算例,验证了所建方程及数值技术在强冲击性流动问题数值模拟中的有效性,并首次实现了液滴冲击液面问题的SPH法数值模拟。