船舶操纵相关粘性流及水动力计算

【摘要】： 采用计算流体动力学(CFD)技术实现在船舶设计阶段对船舶水动力性能的精确预报，是船舶水动力学学科中一个具有重大理论和实用意义的、富有挑战性的课题。近十年来，随着船舶CFD计算方法和计算能力的发展，数值预报精度不断提高，其对船舶优化设计的指导作用也日益突出，并有与船模试验并驾齐驱、取代部分船模试验的趋势。三维船舶粘性流计算方法，具有准确捕捉船体周围粘性流动细节包括船模试验难以测量的流动形态的能力，已成功地应用于船舶快速性方面的阻力预报；在船舶操纵性方面，这类方法虽处于初始发展阶段，但也已获得重大进展，具有精确预报船舶操纵水动力的潜力。本文即在这种背景下，瞄准船舶操纵水动力预报方面的国际前沿和热点课题，通过对现代船舶粘性流计算方法的研究，自主开发了一个船舶操纵粘性流求解器，并将所开发的求解器成功地应用于一系列和船舶操纵问题相关的粘性流动与水动力计算，得到了令人满意的结果。 本文选取RANS方程作为控制方程，并用标准κ-ε湍流模式结合壁函数封闭方程；分析了常用的边界条件，特别是入流边界条件、壁面边界条件和轴边界条件，提出了处理这些边界条件的方法；在分区结构网格上，运用有限体积法(FVM)对RANS方程进行离散，其中，对流项采用混合迎风、中心差分格式，扩散项采用中心差分格式，源项采用中心差分格式并进行部分隐式化以增强离散方程系数矩阵主对角占优；压力、速度等采用非交错配置，并用SIMPLE法耦合求解；对离散后得到的代数方程组选用稳定和收敛性能佳的强隐式法(SIP)迭代求解，并在每步迭代之前对线性方程的矩阵系数进行亚松弛预处理以增强求解的稳定性；为提高数值解的稳定性和收敛性，采用了多重网格法。基于上述方法，采用C++程序设计语言，研究和开发了一个船舶操纵粘性流求解器(VSMAN)。 应用所开发的求解器，以NACA0015翼型舵为算例计算了船舵在不同雷诺 武汉理三「;大学博士学位论文 数下大舵角范围内的三维粘性流场及水动力，成功地预报了舵的失速角和最大升 力，并初步探讨了雷诺数对舵水动力的影响;计算结果与现有试验和计算数据比 较，吻合程度相当好，初步检验和验证了该求解器精确模拟粘性流动和计算水动 力的能力。 应用所开发的求解器，以6:1长椭球体为算例计算了回转体在大攻角下定常 斜航运动时的三维粘性流场及水动力，对层流流动和湍流流动分别进行了计算， 分析了层流与湍流分离流动和涡旋产生的特点及其对水动力的影响;计算结果与 他人的试验和计算数据比较，吻合程度良好，表明应用本求解器能够正确模拟这 种以层流/湍流分离流为主的复杂粘性流动，得到相当精确的水动力，检验和验 证了该求解器精确模拟回转体在大攻角下的粘性分离流动和计算水动力的能力。 应用所开发的求解器，以wigley船型为算例计算了大角度斜航船体粘性流 场和水动力，分析了漂角的变化对船体所受到的粘性水动力的影响，相当精确地 预报了以横流分离和般涡生成与泄出为特征的操纵运动船体特有流动形态及横 向水动力和转脂力矩，经与现有试验和计算数据比较，检验和验证了该求解器精 确模拟绕斜航运动船体的大尺度分离流动和计算非线性水动力的能力。 利用两个开放源代码的软件库FOX和VTK，设计开发了一个高效、友好的 船舶CFD可视化系统(SCFDVS)，实现了以多种方式绘制ZD/3D几何体和计算 网格、ZD/3D标量场和矢量场的功能，并在本文计算的前、后处理过程中得到了 成功的应用。 本文完成的工作缩短了我国在船舶操纵水动力预报方面与国际先进水平 的差距。由于操纵运动粘性流动的复杂性，粘性流方法目前在船舶操纵水动 力计算领域的开发与应用尚处于初始发展阶段，其应用能力和计算精度还不 能满足工程实用的要求。本文在研究基础薄弱、研究条件差的情况下对船舶 操纵粘性流动及水动力计算进行了有益的探索，取得了一定的成功，但仍有 许多问题值得进一步探讨和研究。