Wigley船波浪增阻的研究.docVIP

Wigley船波浪增阻的研究 　　摘 要：以黏性理论为基础，基于N- S方程和标准k-ε湍流模型， 采用VOF方法追踪自由面， 利用CFD方法建立了具有造波和消波功能的数值波浪水池，模拟了波浪的生成与传播。同时，计算了Wigley在规则波浪中的水动力性能。计算结果表明，本文的三维数值波浪水池性能良好，使用的方法可以准确地模拟波浪的运动，计算的波浪增阻为船舶波浪阻力的预报提供了有益的参考。 　　关键词：CFD；数值波浪水池；波浪增阻 　　中图分类号：U661.31+1 文献标识码：A 　　Abstract：Based on the viscosity theory， the N-S equations and the standard k-ε turbulence model are used to calculate the added resistance， the VOF method is used to track the free surface， the numerical wave tank is established to simulate the wave generation and propagation by CFD method. At the same time， the hydrodynamic performance of Wigley in regular waves is calculated. The numerical results show that the performance of the 3D numerical wave tank is good and the method can be used to simulate the wave motion accurately. 　　Keywords：CFD； Numerical wave tank； Added resistance in waves 　　1 引言 　　当海水受海风的作用和气压变化等影响会形成海上波浪，波浪是一种周期性的起伏运动。船舶在波浪中的水动力性能直接影响其能耗的利用，可见波浪对船舶的水动力性能影响至关重要。为了在提高船舶能耗量利用率，减小CO2的排放，船舶耐波性研究成为了21世纪重要的内容之一。到目前为止，船舶在波浪中的阻力性能研究成果已经得到了广泛的展开，如：试验方法[1]、基于势流理论[2]和黏性流理论[3]方法等。随着计算机技术的飞速发展，基于CFD方法建立的数值波浪水池不仅费用低廉而且易于实现。李宏伟[4]采用摇板式造波建立了二 维数值波浪水池，模拟了波浪的生成和传播；封星[5]采用推板法成功建立了二维数值波浪水池；方昭昭[6]采用速度造波法建立数值波浪水池，模拟了顶浪航行船舶的水动力性能；查若思[7]采用Openform自主开发程序，计算了DTMB5512在波浪中的运动响应，得到了很好的计算结果。本文基于前人的基础，以黏性流理论为计算，以FLUENT为基本平台，并对其进行二次开发建立三维数值波浪水池。造波方法采用速度造波方法，消波方法采用人工粘性消波。计算了Wigley在波浪中的水动力性能，采用动网格技术计算了该船在波浪中的纵倾和升沉运动，并与试验值[8]相比，得到了很好的计算结果。本文所使用的波浪阻力计算方法可以准确的预报Wigley的波浪增阻，为船舶在波浪中运动的研究提供了参考。 　　2 基本原理 　　2.1控制方程 　　利用GAMBIT软件建立船体模型，建立长方体作为外流场流域，计算域尺度为5L×1.2L×1.3L。其中：计算域的入口距船首一倍船长；出口距船尾三倍船长；水深为一倍波长。网格划分的好坏直接影响数值计算的结果，因此整个计算模型采用结构网格对计算域进行划分。在水面处网格进行细化，以捕捉自由液面；远离船体部分适当加大网格尺寸，以减小网格数量；船体附近划分边界层网格，以处理船体附近的湍流流动。最后模型网格总数为564 562。 　　5.2 边界条件设置 　　模型入口采用速度入口，并按照造波公式定义波浪垂向和横向的速度，同时在入口处加上船速。出口处设置成压力出口。船体侧面设置成对称边界，以减小计算网格数量。船体表面、上下边界和右边界设置壁面边界。 　　5.3 数值模拟与分析 　　本文利用CFD方法计算了Wigley船在波浪中的水动力性能，其中波浪参数工况如表2所示。造波区波面如图1所示。从图中可以看出造波区的波形稳定，波浪规则，衰减较慢，因此该水流域适合波浪增阻的研究。 　　从图中可以看出：短波长纵摇和升沉传递函数计算的结果较准确，随着波长增加，纵摇和升沉传递函数与试验值差距较大；波浪增阻在短波长计算结果误差偏大，但是总体与实验数据吻合较好，可以