论文部分内容阅读
在船舶与海洋工程领域中,流体力学无论在理论研究还是在工程应用方面均具有重要的意义,而随着数值计算方法的和计算机硬件的发展,计算流体力学已经成为船舶与海洋结构物水动力性能计算和预报的重要工具。对于具有复杂形状问题域的问题,采用结构化网格需要花费大量的时间,相比之下,采用非结构网格的数值计算方法更具优势,而如何计算非结构网格下计算流体力学中的偏微分方程成为开发基于非结构网格数值方法的关键。此外,海洋工程的研究对象往往具有跨越多个数量级的几何尺度,在单一尺度下的数值方法很难同时满足不同尺度下计算精度和成本的需要,而建立宏观和介观数值方法的耦合体系,能够很好地解决这一难题。近年来发展的梯度光滑法(Gradient Smoothing Method,GSM)基于适用于复杂问题域剖分的非结构化网格,采用梯度光滑技术,具有灵活、准确且对网格畸变不敏感等优点。因此,本文开展了非结构网格下梯度光滑法在计算流体力学方面的研究。论文的主要工作如下:(1)在非结构网格下,采用梯度光滑法对对流方程进行数值计算。本文回顾了现存的主要对流格式,并进行了详尽地分析,特别是对于TVD(Total Variation Diminishing)和NVD(Normalized Variable Diagram)算法,对比研究了二者之间的联系。为了能够将基于结构网格上提出的TVD和NVD算法扩展至非结构网格下的梯度光滑法,本文提出了一种基于梯度光滑技术来计算迎风变量的插值方法,并在梯度光滑法的框架下进行计算验证。通过定义迎风点的位置来判断其所在单元,然后根据不同梯度光滑域(节点光滑域、中点光滑域和中心点光滑域)提出了三种插值计算迎风变量的方法,即nGSM(node-based gradient smoothing method),mGSM(midpoint-based gradient smoothing method)和cGSM(centroid-based gradient smoothing method)。在数值实验中,既包括间断问题和连续问题,也包括稳态问题和瞬态问题,并通过与之前方法对比验证了本文方法的准确性。(2)为实现非结构网格下对自由液面的模拟,利用梯度光滑法对VOF(Volume of Fluid)模型进行数值计算。VOF模型的控制方程为对流方程,在结构网格下,通常采用几何重构的方法,但这种方法难以应用于非结构网格。为了克服这一问题,本文采用了基于NVD(Normalized Variable Diagram)概念构造的高精度离散格式,如CICSAM(Compressive Interface Capturing Scheme for Arbitrary Meshes),FBICS(Flux-Blending Interface-Capturing Scheme)以及 CUIBS(Cubic Upwind Interpolation based Blending Scheme),并利用cGSM计算这些高精度格式在非结构网格下所需的迎风变量。数值结果表明非结构网格下采用高精度格式的GSM能够对自由液面进行准确的数值模拟,能够准确预测液面形状并保持界面的锐利性。(3)不可压缩流的数值模拟一直是CFD(Computational Fluid Dynamics,CFD)研究的核心问题,通过求解Navier-Stokes控制方程能够对结构物的水动力性能进行准确预报。在本文中,利用GSM开展对非结构网格下不可压缩流数值计算的研究。应用非结构网格,一方面降低网格划分的时间成本,另一方面通过合理的网格布置提高计算效率。为了解决不可压缩流中的速度和压力耦合问题,在控制方程中引入了人工压缩性项,并通过构造相应的光滑域,利用梯度光滑技术对对流项与粘流项进行离散。在数值算例中,GSM能够灵活地进行网格划分并得到准确的数值结果。此外,还将GSM应用于经典的钝体绕流分析,数值结果显示了不同形状的钝体对尾流的影响,并对比讨论了在定常流动与非定常流动情况下圆柱和三角柱在阻力系数、升力系数以及斯特劳哈尔数随雷诺数的变化趋势。计算结果证明了非结构网格下GSM能够准确、有效地解决基础水动力问题。(4)为了解决多尺度问题,本文在GSM对不可压缩流数值模拟的基础上,引入格子玻尔兹曼方法(Lattice Boltzmann Method,LBM),提出了 GSM-LBM耦合算法。在耦合算法中,将计算域划分为利用GSM计算的宏观子计算域和采用LBM计算的介观子计算域,两种方法通过耦合区域进行流动信息传递。本文提出的GSM-LBM耦合算法在宏观计算域采用了非结构网格,并改进了空间耦合方式。通过数值算例验证了 GSM-LBM耦合算法的准确性和有效性。计算结果表明GSM-LBM耦合算法在计算效率上要高于LBM,而且该方法不仅能够给出整个流场的流动信息,而且还能够描述介观尺度的流动特性。由于GSM采用了非结构网格,可以通过优化网格布置如局部网格加密,进一步提高GSM-LBM耦合算法的计算效率,此外,也有利于模拟计算域形状复杂的多尺度问题。