《基于多参考系和滑移网格模型的吊舱推进器水动力性能研究》由会员分享，可在线阅读，更多相关《基于多参考系和滑移网格模型的吊舱推进器水动力性能研究（6页珍藏版）》请在金锄头文库上搜索。
1、Author:http:/www.5dw.org 青州 seoAuthor:http:/www.5dw.org 寿光 seo基于多参考系和滑移网格模型的吊舱推进器水动力性能研究闯振菊 1, 黄 胜 1, 胡 健 1 ，解学参 1（1,哈尔滨工程大学船舶工程学院，黑龙江省哈尔滨市 150001）摘要：利用 CFD 软件，结合多参考系方法和滑移网格技术对粘性流场中某型拖式吊舱推进器的定常及非定常水动力性能进行了研究，得到了吊舱推进器的敞水性征曲线，以及随着进速系数的变化吊舱推进器在 X,Y,Z 三个方向上推力和转矩的变化规律。在数学的建模过程中，利用 FORTRAN 语言编制了计算吊舱推进器型值点
2、的程序，然后把计算值导入 FLUENT 的前处理器 GAMBIT 建立了光滑的三维计算模型。文中给出了吊舱推进器表面的压力分布图，以及桨盘面处的轴向，径向和切向的速度分布图。计算结果表明，螺旋桨的推力及转矩呈周期性振荡变化，并且振荡的频率以一倍叶频为主，在支架的正前端存在一个轴向和切向速度的高峰区。关键词:吊舱推进器，CFD，多参考系模型，滑移网格模型，水动力性能Abstract: Use CFD method and combined with MRF and Moving Mesh technology to do some research on the podded propeller
3、s hydrodynamic performance in the viscous flow field. Obtained its characteristic curve and the regularity of the force and moment in X,Y,Z directions. The nodes of podded propeller were calculated by FORTRAN program, which were fitted through NUMBS method in GAMBIT so as to establish the calculatin
4、g model. This paper showed the pressure distribution on the blades, and the axial, radial and tangential velocity on the surface of the blade. The results showed that the force and moment of the podded propeller vibrates periodically, and there is a high axial and tangential velocity zone before the
5、 pod.Key Words: podded propeller, CFD, MRF, Moving Mesh, hydrodynamic performance1 引言吊舱式推进器是近年来发展起来的一种新型的船舶推进系统，是目前船舶推进系统研究开发领域引人瞩目的焦点。POD 推进器主要由支架、吊舱和螺旋桨等部件构成。其中，吊舱通过支架悬挂在船体下面，舱体内置电机直接驱动舱体前端和（或）后端的螺旋桨。其设计思想的革命性在于，它把螺旋桨驱动电机置于一个能 360o 回转的吊舱内，悬挂在船下，集推进装置和操舵装置于一体，省去了通常所使用的推进器轴系和舵。POD 推进器将推进系统置于船外，可以节省船
6、体内大量的空间，从而极大地增加了船舶设计、建造和使用的灵活性。目前国内在吊舱推进器方面主要关注其水动力性能，目前研究水动力性能主要有两种方法：势流方法和粘性流方法。势流方法的基本假设是把水看成一种无旋、无粘的理想流体，不考虑水的粘性作用，同时不考虑流体分离。粘性流方法相对于势流理论在上述方面有一定的优越性，可以近似的模拟流场的真实流动。鉴于此文章用选择 CFD 方法来计算吊舱推进器的水动力性能。2 CFD方法的基本理论2.1多参考系模型（MRF）的基本理论MRF 模型是旋转单元体的稳态近似。它求解出来的流场是一个充分发展的流场，这个流场再以一定的速度运动就可以得到实际的流场，多参考系模型方法是
7、近似的，当螺旋桨和舱体之间相互作用相对较弱时可以使用 MRF 模型求解吊舱推进器的定常水动力性能。2.1.1连续性方程旋转坐标系下的连续性方程表示为: Author:http:/www.5dw.org 青州 seoAuthor:http:/www.5dw.org 寿光 seo(1）()rmtVS写成张量形式为：（2）0jjuxt式中， 是密度， 是时间， 是速度矢量。此式是瞬态三维可压流体的质量守恒方程。若流体不可压tu，或者密度 为常数。则其散度形式为0dt div2.1.2 相对速度公式在 MRF 方法中计算区域分为不同的子域，每个子域的控制方程是关于子域参考系而写的。在两子域间的边界，子
8、域的控制方程的扩散项和其他项需要邻近子域的速度值，使用相对速度公式，每个子域的速度相对于子域的运动计算。速度和速度梯度从移动参考系如下描述的转换到绝对惯性系。计算区域旋转轴的初始位置的位置向量定义为：（3）0rx图 1 相对速度坐标系这里的 是笛卡儿坐标的位置向量， 是计算区域旋转轴的初始位置。移动参考系的相对速度可以x0x通过以下方程转换为绝对（静止）参考系的值：（4）()rtvv这里的速度 是绝对惯性参考系的速度， 是相对非惯性参考系的速度值， 是非惯性参考系的平移速vr tv度.根据定义的相对速度，绝对速度向量梯度以下式子给予;（5）)(vr2.2滑移网格模型（Moving Mesh）的
9、基本理论滑移网格模型可使在交界面两侧的网格相互滑动，而不要求交界面两侧的网格结点相互重合。但要计算交界面两侧的通量，并使其相等。为了计算交界面的通量，首先在每一新的时间步确定出交界面两Author:http:/www.5dw.org 青州 seoAuthor:http:/www.5dw.org 寿光 seo边交界区的重合面。基本上，通过网格重合面的通量，是由交界面两边交界区的重合面计算，而不是用整个交界面计算。非定常流动计算中，在流体机械的旋转部件出口与固定部件进口间形成网格滑移的交界面。当转子和定子的交互作用应用实时解法(而不是时均解法)时，必须用滑动网格模型计算非稳态流场。3数值计算模型吊
10、舱推进器螺旋桨的直径为 0.22m，4 叶，盘面比为 0.59，毂径比 0.29，侧斜角为 35 度，变螺距。椭球舱体的长度为 1.667D,最大直径 0.417D，支架高度为 0.79D，弦长 0.75D，厚度为 0.15D，剖面形状为椭圆，D 为螺旋桨直径.文中采用 FORTRAN 语言编制程序，计算出吊舱推进器表面的型值点，将原始的型值点数据转换为特定的格式数据，然后输入到与 Fluent 软件配套的前处理器 Gambit 软件，进行实体几何建模。在建模过程中使用的是直角坐标系 O-XYZ，X 轴方向代表来流方向，它沿着螺旋桨的旋转轴指向下游，Y 轴与螺旋桨的某一桨叶的叶面参考线一致，Z
11、 轴服从右手定则。如图 2，图 3 所示图 2 podded propulsor 三维模型 图 3 吊舱及控制域域的体网格生成建立好吊舱推进器的三维模型后，还要在螺旋桨的四周建立一个小的控制域以及整个吊舱的外部建立一个大的控制域，这样便于在划分网格时进行局部加密，提高计算结果的准确度。大域的直径约为螺旋桨直径的 5 倍，长度约为整个吊舱推进器的 5 倍。在舱体表面和控制域内布置网格，并且设置边界条件，边界条件的设置包括：速度入口，自由流出口，固壁边界等，具体的条件还可以在 FLUENT 中计算时根据需要进行详细的设置。本为划分的网格数目为 2337538。4 数值计算结果4.1 定常水动力性能
12、计算结果在 FLUENT 中利用多参考系方法，结合 RNG 湍流模型，进速系数 分别取为kJ0.2，0.3，0.4，0.5，0.6，0.7，0.8，0.9，0.95，1.0，螺旋桨转速为一定值， =600rpm，计算了敞n水拖式吊舱推进器的性征曲线，并且于实验值做比较，如图 4 所示：Author:http:/www.5dw.org 青州 seoAuthor:http:/www.5dw.org 寿光 seo0.20.30.40.50.60.70.80.91.01.0.5.10.5.20.5.30.5.40.5.0.5.60.5.70.5Kt,1q,J实 验 值 ： Kt 10Kq 计 算 值
13、： t 图 4 拖式吊舱推进器性能理论值与试验值的比较计算值与试验值的比较如图 4 所示，转矩系数 的计算结果与试验结果基本上一致，只是在斜率上TK稍有偏差；而推力系数 在 J=0.8,0.95 的情况下两者基本重合，实验得到的效率曲线要率高于计算值。QK总的来说吻合良好，验证了该计算方法的可靠性.吊舱表面的压力分布如图 5 所示图 5 J=0.4 时吊舱推进器表面的的压力分布从图 5 中可以看出螺旋桨桨叶边缘处的压力较小，总的来说桨叶表面的压力分布比较均匀，在吊舱支架的前端与舱体的交界处存在一个高压区，所以支架与舱体的连接处应该使其尽量光滑避免尖点。4.2 非定常水动力性能计算结果在 FLU
14、ENT 中使用 Moving Mesh 的方法假定流动是不稳定的，模拟出来的流场就是实际的流场，因此可以求解吊舱推进器的非定常水动力性能。如图 610 所示。Author:http:/www.5dw.org 青州 seoAuthor:http:/www.5dw.org 寿光 seo-5005010150202503035040-73.9-73.8-73.-73.6-73.5-73.4proelr foce X(N)(度 ) -50050101502025030350407.7.87.97.807.817.827.837.84pod frce X(N)度 图 6 螺旋桨在 X 方向的推力随旋转角
15、度变化 图 7 吊舱在 X 方向的推力随旋转角度变化-5005010150202503035040-6.05-6.0-65.9-65.90-65.8-65.80-65.7-65.70-65.-65.0pode proulsr foce X(N) (度 ) -50050101502025030350402.32.5342.5362.5382.5402.542.542.5462.548proelr moent X(N-M) (度 )图 8 吊舱推进器在 X 方向的推力随旋转角度变化 图 9 螺旋桨在 X 方向的转矩随旋转角度变化 -5005010150202503035040-0.-0.5-0.4
16、5-0.4-0.35-0.3-0.25-0.2-0.15-0.1proelr moent Y(NM)(度 )图10 螺旋桨在Y方向的转矩随旋转角度变化 图11 桨盘处的轴向速度分布图以上图示均为进速系数J=0.4时，吊舱推进器的的推力，转矩随旋转角度的变化规律。由图7，图8可以看出在螺旋桨工作的过程中其推力在73N附近呈四周期振荡，每当螺旋桨旋转90度就会出现一个推力的高峰值，本文计算的四叶桨最高值可达73.85N，最小值约为73.45N；吊舱在X方向产生的推力约为螺旋桨的十分之一，最高值可达7.83 N，最小值为7.77N,也呈四周期分布，但二者方向相反。所以整个吊舱推进器的推力范围为65.65N66N。由图9，图10可以看出螺旋桨在X和Y方向上的转矩分布，并且X方向上的转矩约为Y方向上转矩的100倍，并且二者方向相反。
