鱼类波状摆动推进的数值模拟.pdfVIP

鱼类波状摆动推进的数值模拟 夏全新 鲁传敬 吴磊 (上海交通大学 船舶海洋与建筑工程学院，上海，200030) 摘要：本文采用有限体积法并利用动网格技术，数值模拟了鱼类的波状摆动推进。首先分析了二维 零厚度波状板的波状摆动推进，研究了在高雷诺数下匀速巡游的流场细节，发现了尾迹中涡的形态是呈之 字形的涡环链。并同刚性平板的绕流结果比较，验证了波状摆动可以充分抑制湍流，减小表面所受的粘性 阻力。其次使用NACA 标准翼型模拟了有厚度的鱼体波状摆动，得到了尾部推力型的反卡门涡街。还分析 了不同运动参数以及鱼体厚度对推力和阻力的影响。最后计算了边界条件对鱼类运动的影响。 关键词：外部生物流体力学，动网格，漩涡脱落，反卡门涡街 1 引言 海洋中的鱼类经历了长时间的自然选择，已经充分发展和完善了在较低能耗条件下，快 速前进的能力。关于鱼类运动的流体力学原理的研究，有助于对水生动物形态进化的认识， 并有利于实现在高效推进下如何保持高机动性、稳定性及低噪声等问题。 在上世纪 50年代之前，人们对于水生动物运动的研究限于实验和观测，Taylor[1]率先对 微生物和蠕虫作细长体运动分析，发展了抗力理论。而从二十世纪六十年代开始，Lighthill[2] 和 Wu, T.Y. [3,4] （吴耀祖）在鱼类运动的理论研究方面作了开创性的研究。Lighthill 根据鱼类 身体细长的体形特点，利用细长体理论对鱼类运动的动力学问题作了系统的研究，并利用二 维非定常机翼理论研究了月牙尾的运动。Wu, T.Y.系统地发展了二维波动板理论，并讨论了 [5] [6] 鱼类的最佳运动方式。国内的童秉纲研究组（Chen，Zhang和 Tong ，童秉纲和庄礼贤 ） 利用半解析和半数值的方法提出了三维波动板理论，求解了鱼类运动的正问题和反问题。 在实验方面，国内外很多学者利用了 3D PIV 仪器对鱼类运动的整个流场进行了运动显 示。实验结果都显示，鱼类在匀速巡游状态时脱落的尾涡形态为反卡门涡街，此时鱼类运动 的效率极高，同时还有低噪声，高推进效率等很好的水力特性。并且MIT 的M.Triantafyliou 研究群成功制造了可以完全模拟波状摆动的机器鱼。此外，Liu ，Wassersug 和 Kawachi[7]使 用三维模型数值分析蝌蚪的推进问题。Chen 和 Doi[8,9]利用 MAC 方法数值模拟了鱼类的波 状摆动问题。以上各种结果都表明，鱼类的波状摆动能取得很好的机动性，并降低体表的湍 流度，在运动中取得很好的运动效率。 本文利用 Fluent 6.1 中的动网格技术，应用有限体积法对二维鱼类的波状摆动问题进行 了数值模拟。研究了二维波动板在高雷诺数下所受的推力和阻力与板的运动方程中各个参数 之间的关系。计算得到了鱼类在高雷诺数下匀速巡游的流场。此后又用标准 NACA 翼型模 拟了鱼体，研究了厚度对鱼游动的影响。还计算了各种边界条件对鱼类运动的影响。 2 数学模型及数值算法 2.1 数学模型 鱼类运动所处的流体——水一般可以假定为不可压缩流体，控制方程为连续性方程及与 时间相关的N-S 方程。 ∂W ∂F ∂G ∂H ∂F ∂G ∂H l l l V V v + + + + + ∂t ∂x ∂y ∂z ∂x ∂y ∂z 其中 ρ ⎡0 ⎤ ⎡ 0 ⎤ ⎡0 ⎤ ⎡ ⎤ ⎢ ⎥ ⎢τ ⎥ ⎢τ ⎥ ⎢τ ⎥ ρu ⎢ xx ⎥ ⎢ yx ⎥ ⎢ zx ⎥ W ⎢ ⎥， F ，G ，H ⎢ ⎥ V ⎢τ ⎥ V ⎢τ ⎥ V