螺旋桨诱导流场的数值分析.docxVIP

螺旋桨诱导流场的数值分析 海洋螺母的尾流场一直是人们关注的研究对象之一。在研究螺旋尾流场驱动性能以及螺旋桨与船尾、船体和其他固定附件之间的干扰时，应详细研究螺钉的尾流场。船舶螺钉后尾流旋转水流带的总能量约为总能量的15%30%。由于后缸的重量比后缸的加速运动减少，该功率因反向激励螺钉、电机底板、桨后固定导叶、反应舵、动力船等装置的功率，我们必须对螺钉的尾流场进行预处理，以便对这些装置进行改进。同时，计算螺钉尾流场是分析螺旋桨振动和相关干扰的关键之一。 关于螺旋桨尾流场的计算,从发表的研究成果来看,多为势流理论.文献用面元法计算了P4119桨的尾流场.文献用准非线性升力线理论计算MAU桨后流场速度分布;文献用升力面方法计算螺旋桨的尾流场;该文用计算流体力学(CFD)对螺旋桨尾流场进行数值预报,并对计算结果作了相应的分析. 1 基于简单算法的流场价值计算 1.1 压力场的求解过程 求解压力耦合方程组的半隐式方法 (semi-implicit method for pressure-linked equations, SIMPLE)算法的基本思想可描述如下:对于给定的压力场(它可以是假定的值,或是上一次迭代计算所得到的结果),求解离散形式的动量方程,得出速度场.因为压力场是假定的或不精确的,这样,由此得到的速度场一般不满足连续方程,因此,必须对给定的压力场加以修正.修正的原则是:与修正后的压力场相对应的速度场能满足这一迭代层次上的连续方程.据此原则,把由动量方程的离散形式所规定的压力与速度的关系代入连续方程的离散形式,从而得到压力修正方程,由压力修正方程得出压力修正值.接着,根据修正后的压力场,求得新的速度场.然后检查速度场是否收敛.若不收敛,用修正后的压力值作为给定的压力场,开始下一层次的计算.如此反复,直到获得收敛的解. 在上述求解过程中,如何获得压力修正值(即如何构造压力修正方程),以及如何根据压力修正值确定正确的速度(即如何构造速度修正方程),是SIMPLE算法的2个关键问题. 1.2 压力场与速度场 设有初始的猜测压力场p*,动量方程的离散方程可借助该压力场得以求解,从而求出相应的速度分量v*和v*,根据动量方程的离散方程,有 假设压力修正值p′为正确的压力场p与猜测的压力场p*之差,有 同样,定义速度修正值u′和v′,以联系正确的速度场(u,v)与猜测的速度场(u′,v′),有 将正确的压力场p代入动量离散方程,得到正确的速度场(u,v).假定源项b不变,并引入压力修正与速度修正的表达式.可得 可以看出,由压力修正p′可求出速度修正(u′,v′),便可对猜测的速度场作出相应的速度修正,得到正确的速度场(u,v). 1.3 瞬态误差的压力修正方程 在速度修正方程中,只考虑了动量方程,实际上,速度场还受连续方程的约束.连续方程的离散方程为 将正确的速度值代入(5),整理可得压力修正值p′的离散方程并简记为 求解方程(9),可得到空间所有位置的压力修正值p′.在瞬态问题中,动量方程变化不大,压力修正方程需要重新建立.将连续方程在二维空间的一个标量控制体积上进行积分有 从上式导出的压力修正方程,必须要包含代表瞬态特性的项. 1.4 压力场及速度场 基于交错网格的SIMPLE算法的求解步骤: 1) 假设一个速度分布u*,v*,?*(u,v,?)用于计算首次迭代时的动量离散方程中的系数和常数项. 2) 假设一个压力场,即给定压力猜测值p*. 3) 根据当前速度场及压力场,计算动量离散方程中的系数和常数项. 4) 求解动量离散方程(1)、(2). 5) 根据速度u*、v*,求解压力修正方程(9). 6) 修正压力与速度. 7) 求解所有其他的离散化输运方程(视需要进行). 8) 判断收敛性,如收敛,结束;否则,令p*=p,u*=u,v*=v,?*=?,转入第3)步,重新迭代计算,直到收敛为止. 对于瞬态问题的流动计算,借助隐式时间积分方案,在每个时间步内进行迭代,就好像在调用SIMPLE算法进行普通稳态问题的迭代计算一样,直到取得本时间步的收敛解,然后转入下个时间步继续重复上述过程. 2 标准k-模型 高速旋转叶轮的雷诺数超过相应的临界雷诺数时会诱发湍流.湍流是一种复杂的非稳态流动,要使方程组封闭,必须做出假设,即建立模型.这种模型把未知的更高阶的时间平均值表示成较低阶的在计算中可以确定的量的函数.采用Reynolds(雷诺)时均方程方法(RANS)研究湍流模型. 根据时间平均法,用时均值与脉动值之和代替流动变量,即 并将式(11)代入瞬时状态下的连续方程和动量方程,可得用张量的形式表示的时均连续方程和Reynolds方程如下: 脉动值乘积项?(ρu′iu′jˉˉˉˉˉˉˉˉ)-(ρu′iu′jˉ)称作雷诺应力,是一个二阶张量,正是