船舶水动力的粘势流耦合方法初探-船舶与海洋工程专业论文.docx

万方数据 万方数据 Dissertation for the Master’s Degree in Engineering Study on Viscous/Potential Method for ship hydrodynamics Author: Ji Zhao Specialty: Naval Architecture and Ocean Engineering Supervisor: Prof. Zhu Renchuan School of Naval Architecture, Ocean and Civil Engineering Shanghai Jiao Tong Univerisity January.2015 上海交通大学硕士学位论文 上海交通大学硕士学位论文 摘 要 船舶水动力性能预报是船舶生产设计过程中的重要环节，船舶 水动力数值计算模拟也逐渐成为辅助船舶设计的重要手段，应用高 效准确的数值计算理论方法预报船舶水动力性能将会提升船舶设计 能力，提高设计质量、缩短设计周期。目前随着计算机的发展，基 于势流理论的船舶航行兴波、切片理论方法等，在生产设计中已被 广泛用于预报船舶阻力与耐波性能，尽管计算流体动力学（CFD） 方法更为全面真实的考虑了流体粘性属性和流动非线性，但应用 CFD 方法计算船体水动力性能时占有大量的计算机资源，计算时间 长，甚至导致生产设计周期过长，使其在船舶设计中应用受到较大 的限制。本文采用速度分解方法，充分利用了势流理论快速计算与 CFD 方法全面真实的特性，结合两者优势给出了粘势流耦合求解船 舶水动力问题的方法，基于开源 CFD 软件 OpenFOAM 平台，开发 了相应的粘势流耦合求解软件，并验证了理论方法和所开发的软件 的正确性。 首先，论文采用 Helmholtz 速度分解，将流场速度分解为势流 部分和非势流部分，剥离势流部分后得到粘势流耦合的变形 NS 方程， 并进一步给出了耦合势流理论求解变形 NS 方程的数值计算方法。以 二维圆柱绕流为例，在 OpenFOAM 平台上自主编程，实现了耦合二 维圆柱绕流解析解的变形 NS 方程的求解。对流场中的速度、压力， 与文献中的试验及传统 CFD 等的结果进行仔细比较和分析，验证了 该方法理论的正确性以及计算结果的准确性。 在此基础上，论文将该理论的应用推广至三维，并就三维潜体 绕流进行了粘势流耦合求解。三维潜体绕流的势流问题直接采用了 -I - 三维简单格林函数法进行了数值计算求解，并提出了源偶偏移法以 确定变形 NS 方程相应的定解边界条件。自主编程在 OpenFOAM 平 台上实现了耦合三维潜体绕流势流解的变形 NS 方程的求解。通过流 场细节分析和最终计算的阻力与实验数据的对比，验证了粘势流耦 合方法用于计算三维绕流问题的准确性。 最后论文以粘势流耦合方法计算 Wigley 船体阻力为例，初步探 索了其在船舶水动力性能计算上的应用。船体阻力为兴波阻力和粘 性阻力之和。Wigley 船航行兴波阻力采用势流方法求解，粘性阻力 采用了粘势流耦合方法求解叠模绕流获得。船体叠模绕流的势流解 先由简单格林函数法计算获得，以此为基础进一步求解了变形 NS 方程的定解问题。获得的船体阻力，通过与实验及文献上相关数据 的进行了比较和分析，结果吻合良好。 论文研究表明该方法正确，适应性强，为改善粘流计算效率提 供了新的途径。 关键字：粘势流耦合，OpenFOAM，变形的 NS 方程，Helmholtz 分 解 -II Study on Viscous/Potential Method for ship hydrodynamics ABSTRACT The prediction of ship hydrodynamic performance is an important part in the process of ship design and production while numerical simulation of ship hydrodynamic has becoming a common method to assistant ship design, so it is vital to explore new numerical algorithms with better efficiency. The potential theory, which includes the using in wavemaking and Strip theory, has been widely used in prediction of ship resistance and seakeeping performance al