强迫运动柱体附加质量与阻尼系数CFD 计算.docVIP

. . 　　 本文档下载自HYPERLINK 360文档中心，更多免费文档请访问网址  HYPERLINK /doc/info-07e9f0e9e009581b6bd9ebb7.html /doc/info-07e9f0e9e009581b6bd9ebb7.html 强迫运动柱体附加质量与阻尼系数的CFD计算 DOI：10.3969/j.issn1000-4874.2011.04.016 510 水 动 力 学 研 究 与 进 展 A辑2011年第4期 1 引言 对船舶运动方程中的水动力系数进行求解，是准确预报船舶在波浪中的运动状态的前提。目前船舶耐波性的预报都主要基于势流理论。对于纵向运动，它们受到的粘性影响较小，基于势流理论的纵向运动预报尚能满足工程精度要求，而横摇运动受粘性影响大，其运动预报精度取决于横摇阻尼的准确确定。横摇阻尼系数的确定一直是人们研究的重点，它是影响其摇荡幅值的重要参数。工程应用上粘性横摇阻尼通常采用模型试验或半经验公式得到[1,2]。 近几十年来，随着计算机的硬件设施和计算技术不断更新，计算流体力学发展迅猛，它广泛应用于船舶性能的各个领域，现已能计算粘性绕流场和船舶的总阻力[4]，并对螺旋桨的性能预报也有十分重要的应用[5,6]。在操纵性和耐波性范畴内，往往伴随着航行状态和运动速度的变化，多为非定常运动。已有不少学者运用CFD方法对船舶操纵性水动力导数进行了求解，得到了较好的结果[7]。在耐波性方面，国内外也有许多学者对粘性流中船舶的运动受力进行了研究[8,9]。在船舶六自由度摇荡中，横摇是受粘性影响最大的运动，它的产生有众多非线性因素，暂时难以得到其精确解，而运用CFD方法对其进行研究已成一种趋势[10-12]。在此基础上，现已逐渐发展成了一种新的方法，它通过求解RANS方程数值模拟传统物理试验[13,14]。该方法还比较新，有较大的发展空间。基于CFD方法的数值模拟试验费用低，无尺度比效应，而且能够方便地对整个试验中的条件和参数进行修改，在一定程度上能够给最优化设计提供更好的参考依据。 本文基于商业软件Fluent对船舶剖面强迫摇荡运动进行数值模拟。通过监测其所受流体作用力，对运动方程中的附加质量和阻尼系数进行预报和分析，并研究了舭龙骨对横摇运动的影响，为进一步细致描述船舶周围流场和预报船舶在波浪中的运动打下了坚实的基础。 ?ρ? (ρui)=0, ?t?xi ???p(ρui) (ρuiuj)=? ?t?xj?xi ??xj ???u??i ?ρuiuj?? ρfi (1) ?μ??????????xj 式中：fi为质量力；ui为流体质点速度在i方向的分量；对于两相流来说，流体平均密度可定义为 ρ=∑aqρq，其中体积分数aq表示单元内第q相流 q=1 2 体占的体积与单元体积的比例，故有∑aq=1；μ为 q=1 2 相体积分数平均的动力粘性系数，其定义与密度的定义类似。 本文通过VOF方法对自由面的位置进行捕捉。整个水槽宽度深度都取得尽量大，而且距船体较远的区域网格尺寸较大，它起到了一定的消波作用，可以避免波遇到固壁面后产生二次反射干扰。 2.2 数值模拟结果分析表达 船体剖面在第k个模态上的运动可表示为 xk=xk0sinωt (2) 其中xk0表示运动的幅值，ω为摇荡频率。 剖面所受流体作用力表示为： xkμjk xkλjk Fjk=0 (3) 其中：k,j=2,3,4分别为横荡、升沉和横摇；xk为第k个模态的位移；μjk与λjk分别为各模态的附加质量与阻尼系数；Fjk表示流体动力或力矩。在后面，力矩用M来表示。 流体对物体的作用力，可以由压力沿着物体的湿表面进行积分得到，该部分数据处理由fluent软件UDF模块编程实现。物面单元的压力值p由动压力pd和静压力ps组成，p=pd ps，其中当y≥0时有ps=0，当ylt;0时有ps=?ρgy，而y为每个物面单元中心的垂向坐标。将pd沿物面积分，可以得到船体在各个方向所受的流体动力和力矩。例如x方向的作用力为 ??? 2 数值计算方法 2.1 控制方程 基于粘性流理论，本文模拟船舶在二维数值水槽内作强迫运动，其控制方程为雷诺平均的N-S方程（RANS）， 罗敏莉，等：强迫运动柱体附加质量与阻尼系数的CFD计算 511 Fx= 所受力矩为 C0 ∫pndl d x 以横摇运动为例，船舶所受流体动力满足 ??? ? ?μ44x4 λ44x4 M44=0 (4) ???? ?