鱼呈S形游动时的高推力产生机制研究.pdfVIP

鱼呈S形游动时的高推力产生机制* 王 亮 吴锤结一 (1．河海大学水利水电工程学院，南京，210098) (2．解放军理工大学理学院流体力学研究中心，南京，211101) 摘 要：本文研究了鱼呈S形路径游动时所产生的推力，并与相同条件下以 直线路径游动时所产生的推力进行了比较．结果表明：当鱼呈S形路径游动时， 鱼身体前部与来流有一个夹角，从而引起分离，形成前缘涡．通过调整鱼的摆动 频率和初始相位，可以令前缘涡与尾缘脱落涡发生相互作用，使反卡门涡街中涡 的强度增加了199．6％，推力提高了674．8％。 关键词：前缘涡 S形游动 反卡门涡街 动边界计算流体力学 1 引言 自然界在其几十亿年的演化过程中创造出千奇百怪的生物，它们以各自截然不同的、 优美的和最节约能量的运动方式往来于世界各地。但是，至今人们对生物运动过程中涉及 的流体力学机理还知之甚少。探索鱼类游动的奥秘，揭示鱼类运动的规律[4]，进而创造出 最优仿生运输工具，是人类不断追求的目标。 在对振荡翼型的研究[1]中发现：当在翼型的前缘附近产生分离时，就会形成前缘涡， 从而引起动态失速。在特定的组合参数条件下，前缘涡会和后缘涡发生合并。前缘涡(动 态失速涡)的形成不仅与高升力系数有关，而且还与高推进效率有关。这使人联想到自然 界中许多鱼在游动时，其运动轨迹不是直线而是呈S形的曲线。当鱼呈S形路径游动时， 鱼身体前部与来流有一个夹角，从而引起分离，形成前缘涡。本文的目的就是为了研究鱼 身体前部产生的涡对推力大小的影响。从广义上讲，是研究鱼如何利用来流中的涡提高 推力。 2数值算法 方程， 【j。声·矗ds—o 。 中国科学院力学研所非线力学国家重点实验室、国家自然科学基金资助项目(#lOl72095)和江苏省自然科 学基金重点课题(#BK2003208)资助项目 一 Email：cjwu@jtonline．corn或cjwu@mech．pku．edu．an 7 其中变量等符号与流体力学界常用习惯一致。对于时间离散，采用二阶精度的三层隐 式格式， ㈣ 似为： 。’ j。fdS一莩jScfds一莓五sc≈莓^sc，c—e,w,n,s 其中，-c为，在‘c’面上的平均值；^为厂在‘e’面中心的值。而体积分可以近似为： I gdl2≈gP△Q‘ (4) J0 其中，g。为被积函数在控制体中心P的值。因面积分或体积分的离散精度均为二阶，故 总精度为二阶。 鱼在摆动过程中，网格是在不断变化的。网格边界的运动会产生边界速度，如果它也 用来计算通量的话，就会导致质量等标量的不守恒。该误差可以通过应用体积守恒定律 Conservation (Space Law，简称SCL)来消除。SCL可以被看作流速为零的连续性方程， 其形式为： ’ 未．rndQ一-r。诜·矗ds—o (5) 该方程表明当控制体的形状和位置随时间变化时，体积具有守恒性。将上式离散得： …f—P，硼m。’’‘ 坠业兰塑一∑(磊．丢)…S t z一、。。 利用隐式欧拉格式，连续性方程可以离散为： 垒垒堕∑二∑二生!螋一∑，础。一0At ZJ。 ，c—P，叫，咒，s (7) 其中，非定常项必须用体积守恒定律来处理。对于不可压流体，网格运动对质量通量