带自由面船体周围粘性流场数值模拟.pdfVIP

第 29卷 第2期 江苏船舶 Vol_29 No．2 2012年 04月 JIANGSU SHIP Apr．2012 带 自由面船体周 围粘性流场数值模拟 于 刚，刘建闽，蔡振雄 (集美大学轮机工程学院，福建 厦门361021) 摘 要：利用FLUENT流体软件，通过求解三维粘性不可压RANS方程和使用标准k．s湍流模型对带自由液面船 体进行了周围粘性流场的数值模拟。自由液面采用 VOF方法处理 ，得出了不同弗劳德数下船体的粘性阻力系 数 ，并与 ITrC1957给出的平板摩擦阻力系数计算公式以及巴普米尔粘压阻力系数公式丽者所求的粘性阻力系 数相 比较 ，验证了本文方法的可行性 ，为船舶阻力预报提供了更为精确的流体分析方案。同时 ，利用 VOF方法 可以很好的捕捉自由液面波形，反应了船体的兴波现象，为船体优化设计提供重要的参考。 关键词：FLUENT；数值模拟 ；粘性流场；自由液面 中图分类号：U661．1 文献标识码：A 0 引言 1 计算模型 船舶设计中，船体的阻力预报精度具有极其重 1．1 物理模型 要的作用，对设计性能 良好的船舶具有重要意义。 本文采用的三维船体模型是根据文献 [2]的 船舶是一个复杂结构体，在海上航行时受到船体周 WIGELY船模，通过三维建模软件 UG生成的船体 围粘性流场作用和波浪对船体的兴波作用。在这种 曲面模型，转换文件格式导人 ICEM．CFD进行网格 流固耦合系统中可看作是非线性不可压缩流体的作 划分。船体建模的主要参数见表 1。其中湿水面积 用，并符合连续性方程和 Navier．Stokes(N．S)方程 为FLUENT求得 。船体外壳示意图如图1所示。 式。由于这种复杂作用，考虑带 自由液面船体周围 表 1 模型船体主尺度 粘性流场情况下的预报 ，在水动力研究中仍是一个 船长／m 船宽／m 吃水／m 湿水面积／m 难题 。很长一段时间内，粘性流场绕流和 自由液面 10 1 O．5 12．43 兴波阻力的计算是分开考虑的，即用势流理论方法 处理 自由面，另外通过求解 RANS方程来计算船体 的粘性边界层…，这种分离方法忽略了自由面对粘 性的影响。近年来 由于高速计算机的迅速发展和普 及，加上数值方法的进步，基于计算流体力学理论的 数值模拟成为研究复杂湍流的必要手段。 本文采用VOF方法，对 WIGELY船在不同弗劳 德数情况下周围粘性流场绕流问题进行数值模拟。 通过与 ITTC1957给出的平板摩擦阻力系数计算公 式 以及 巴普米尔粘压阻力系数公式两者所求的粘性 流体控制区域为一长方体，区域大小根据实际 阻力系数相比较，验证了此方法的可靠性，为解决带 计算决定。若区域太大，则计算网格数量庞大，增加 自由液面船体粘性绕流预报提供了有效的方法。 了不必要的计算量。若区域太小，则不能满足计算 精度。因此综合以上考虑，计算区域为纵向从船艏 收稿 日期：2011—11—18 向前部 1倍船长，船艉向后3倍船长。横向上向两 基金项 目：福建省自然科学基金项 目(2009J01010) 边分别延伸5倍半宽。垂向上 自自由液面向上延伸 作者简介：于刚(1982一)，男，硕士研究生 ，主要从事流固耦合、船体