论文部分内容阅读
自由表面流动问题一直是水动力学研究中的热点问题。其中,物体入水问题由于有着广泛的工程应用背景而始终为许多研究者所关注。在对物体入水问题的模拟过程中,自由面的大变形问题始终是研究的难点所在。而易于处理自由面的大变形过程则是SPH方法最大的优点,因此本文采用SPH方法研究了物体入水等自由表面流动问题。本文推导了SPH方法的流体动力学控制方程,自主研发了基于SPH方法的流场计算软件,通过大量算例验证了自主开发的计算软件,并用以模拟物体的低速入水问题,给出了圆盘入水问题及楔形体入水问题的数值模拟结果。本文的主要工作如下:1.自行开发了基于SPH方法的并行化、水-气两相流流场计算软件,对各种不同的人工粘性项形式、密度重新初始化公式和边界处理方法等各种数值处理技术进行了测试,并测试了单相流SPH模型和两相流SPH模型。通过大量的算例对计算程序进行了验证。2.总结了SPH方法的各种常用固壁边界处理技术,并分析了各种处理方法的优缺点。本文对现有的处理方法进行了改进,提出了一种改进的固壁边界处理方法。在改进方法中,固壁边界点上的压强通过其附近的流体质点的压强插值得到,能够得到更准确的压强,进而得到更为准确的流场情况。通过具体的计算例子验证了本文提出的改进方法的有效性。3.提出了无反射边界处理方法,在固壁边界附近设置与数值波浪水槽中的消波层相类似的海绵层吸收到达固壁上的压力波,消除压力波在边界上的反射,避免反射波对流场造成干扰,从而可以采用较小的计算域进行计算,大大节省了计算量。4.实现了SPH模型的并行化方法,并对并行化SPH程序的并行效率进行了详细的分析讨论。5.应用本文自行开发的SPH流场计算程序研究了物体入水问题。首先采用柱坐标下轴对称形式的SPH方法研究了圆盘的入水过程,比较了不同Froude数下液面的封闭形式,研究了圆盘相对潜深与Froude数的关系,并与文献的实验结果进行了比较。采用并行化单相SPH程序模拟了楔形体入水初期过程。对比了计算及实验的边界压力分布及楔形体受力情况,验证了边界处理方法,研究了不同初始入水速度对楔形体入水初期最大阻力的影响。采用并行化单相及两相SPH程序分别模拟了楔形体入水的全过程,并进行了比较分析。最后通过自行设计的楔形体入水实验对计算结果进行了验证。