论文部分内容阅读
涉及运动物体的多相流动是近年来船舶与海洋工程水动力学中研究的热点问题。对于任意复杂形状的多体运动,非贴体笛卡尔网格比贴体网格更加具有优势,其中切割网格有限体积法不仅保留了笛卡尔网格的优势,同时又具有守恒特性,非常适合于描述多相流中的运动界面。然而,由于典型的多体、多相流问题中的强非线性自由表面和高雷诺数等具有强对流特征,导致求解的压力、速度不连续或出现尖锐梯度解。传统的线性高阶对流离散格式在这种间断或大梯度解附近会产生非物理振荡现象,因此需要采用非线性的高精度格式,其中TVD通量限制器格式应用广泛且适用于规则的网格。因而本文对限制器格式进行了研究,并在切割网格法中引入了TVD通量限制器格式,抑制间断或大梯度解附近的非物理振荡。文章的主要内容如下:第1章阐述了笛卡尔切割网格有限体积法对于研究非线性多体两相流问题的优势,并突出在数值求解该类问题的过程中高精度TVD限制器格式的意义。总结了限制器格式的特性研究现状以及基于各类网格的限制器格式在CFD中的研究现状。并且对求解多相流问题的三大类湍流模型进行分析,表明了大涡模拟的优势。第2章分析了对流项离散及通量的求解思路,总结了切割网格中有限体积法的求解及目前在切割单元中的梯度重构方法。其次比较了TVD限制器格式相对其他高精度格式的优势,对有限体积TVD限制器格式的原理、表达式构成、约束条件进行了详细阐述。此外,还将斜率限制器转变为通量限制器,分析了通量限制器函数约束区域图的内涵及TVD方法与NVD方法的联系和差异,对两种框架中的限制器函数进行特征分类。最后简单描述了多维问题的维度分裂标量限制器算法,以及方程组求解中基于特征分裂作用于数值波的有限体积限制器格式。第3章对本实验室大涡模拟平台进行介绍,阐述了多体质量力模型的原理与特点、分析了隐式切割网格有限体积算法的具体求解过程,并重点提出了基于切割网格法的高精度TVD限制器格式。第4章基于本实验室的大涡模拟求解器,结合质量力模型和切割网格TVD限制器格式对物体入水的案例进行数值模拟,并将限制器格式的结果与实验进行了比较,验证了本文的限制器格式用于切割网格法的有效性。本文基于规则体的多体质量力模型切割网格有限体积法Navier-Stokes方程求解器,通过切割网格TVD通量限制器高精度格式,自适应调整速度梯度的变化,实现计算的高精度和强稳定性。与目前的切割网格中基于梯度重构的斜率限制器不同,本文的切割网格TVD通量限制器格式与切割网格单元的尺寸无关,其思想是对反扩散通量进行修正,具有概念简单、实施方便且确保离散方程对角占优的特性。