迎浪航行条件下船舶耐波性的计算分析.pdfVIP

迎浪航行条件下船舶耐波性的计算分析 查若思，叶海轩，万德成 (上海交通大学船舶海洋与建筑工程学院，海洋工程国家重点实验室，上海200240) 摘要 Fluid DynaIllic)方法对某型模型尺度下的海洋巡逻船在规则波中迎浪航 本文采用CFD(Computational 行的耐波性问题进行了计算分析。采用基于开源CFD软件OpenFOAM自主开发的求解器 naoe．FOAM．SJTU进行计算，包括完整的造波和消波模块，求解船体六自由度运动方程，运用有限体积法 (FVM)离散不可压RANS方程，并使用VOF方法捕捉自由面，最终得到不同工况下的船体垂荡和纵摇 运动响应。计算结果可以为船型优化和设计提供参考。 关键词：运动响应；粘流；耐波性；naoe．FOAM．SJTU求解器 1引言 船舶在波浪中航行时，受到波浪和流的影响，船体会发生一定的摇荡运动，会对船舶的耐波性 和快速性造成不利影响。因此，对于船舶在波浪中的运动响应预报是船舶耐、波Jf生的一个重要课题。 本文利用CFD数值计算的方法研究了某型海洋巡逻船在不同规则波条件下迎浪航行的运动响应情 况。传统的研究方法主要有基于势流理论的计算和模型试验方法，相对于势流方法简单方便的特点， CFD方法充分考虑了流场的粘性，并且能求解传统方法难以解决的强非线性问题。随着计算机硬件 和数值方法的发展，CFD粘流计算的方法越来越成为一个有力的研究工具。相对于直接使用商业软 件，基于开源软件进行自主开发的方式具有更高的灵活性和拓展性。本文的数值模拟使用了基于开 块，该数值水池能够模拟不同海况和不同波浪条件，与理论值符合较好。另外，该求解器基于 OpenFOAM进一步开发了六自由度运动模块，能够处理船舶在波浪作用下的六自由度运动问题。该 求解器适用于求解多种复杂的船舶与海洋工程水动力学问题，文献[3，4，8，9，10，1l，12，13】利用 行了数值模拟，并与试验结果进行了对比；同时在迎浪条件下考虑船体纵倾和升沉运动，计算了船 舶波浪增阻，验证了该求解器求解船舶在波浪中强非线性大幅度运动问题的有效性。 2数学模型和数值方法 本文采用两相不可压RANS方程作为控制方程，其形式可写为： V·U=0 (1) (2) 其中(1)为连续性方程，(2)为动量守恒方程：U代表速度场，％表示网格移动速度。pd为动压力， 其数值等于总压力值减去静水压力，P为液体或者气体的密度，窖为重力加速度向量。∥。ff表示有效 动力粘性，1，表示运动粘度，1，。表示涡粘度。厶为表面张力项，五是用于消波的源项。 本文采用SSTk—gO湍流模型璐1进行方程的闭合。k表示流体质点的湍动能，@表示特征耗散 率。该湍流模型综合了标准k一国模型和七一占模型的优点，既不受自由面的影响，又能保证在壁面 处求解的精确性和可靠性。 of 本文使用带有人工可压缩项的VOF(VolumeFluid)方法[61来处理自由面。VOF输运方程定 义为： 警+v．[(u—U)cr】+口Iu。(1一a)aJ：0 (3) 其中，珥为用于压缩界面的速度场，仅为两相流体的体积分数，代表液体部分所占体积的百分比， 口=0表示空气，口=l表示水，0口l表示界面。 可以通过设置入口边界条件来实现模拟造波板运动，从而实现数值造波门恸能，本文采用的是 考虑遭遇频率的一阶stokes深水波用来产生不同规则波的工况。消波功能是通过在数值水池的末端 设置海绵层来实现的，用来避免波浪反射，造成数值计算结果的居《烈震荡。 本文中采用的naoe—FOAM．SJTU求解器包含了完整的六自由度运动【8’9】求解模块，并采用了动网 格技术，用来进行运动响应预报。坐标系包括船体坐标系和大地坐标系，其中船体坐标系原点位于 船体重心，并随船一起运动。在两种坐标系下的速度可以变换矩阵进行相互转化。船体