摆动水翼绕流的数值研究1).pdfVIP

摆动水翼绕流的数值研究1) 何国毅* 张星† 2) 张曙光* * （北京航空航天大学航空科学与工程学院，北京 100083 ） † （中国科学院力学研究所LNM室，北京 100080 ） 摘要 鱼类经常采用垂直流向的摆动进行游动，这种摆动可以用行进波来表示。本文用浸入边界方法 （Immersed Boundary Method ）模拟了低雷诺数条件下摆动水翼NACA65-010 的绕流问题并研究了雷诺数 对水生动物推进效率的影响。结果表明：随着雷诺数的增大，推力系数和推进效率增大，而功率系数减小； 在 Re 小于 20 时，推力系数，推进效率和功率系数的变化尤为剧烈。随雷诺数增加，由于水翼摆动诱导的 流场变化也更加复杂，水翼后缘处的涡量场强度逐渐增强。摆动诱导反卡门涡街产生推力。 关键词 水翼，摆动，推进，浸入边界方法，反卡门涡街 中图分类号：O257, O351 引言1 要动网格技术，大大减少了数值模拟的成本。因 此我们将采用浸入边界方法研究生物的游动。本 在一定的雷诺数范围内，海洋和水生动物经 文将利用浸入边界法模拟NACA65-010 水翼的摆 常采用身体的摆动进行游动[1] 。这些动物因身体 动绕流，计算水翼摆动产生的推力、升力和效率， 做垂直于流向的摆动，导致动物沿水平方向运 并分析水翼绕流的流动特征，研究的推进机理。 动。在Lighthill 的先驱性工作中[2]，他把垂向的摆 动表示为行进波，发现改变行进波的方式，可以 1 控制方程和数值方法 减少游动的阻力，增加推力，提高推进效率。进 一步的试验[3][4]表明推进效率与垂向摆动的频率 1.1 水翼的数学模型 有关。然而，现有的实验方法很难测量海洋和水 我们以NACA65-010 水翼为例研究摆动推 生游动时的阻力和升力，也很难捕捉摆动物体绕 进的性能与机理。NACA65-010 是一种典型的低 流的流动特征。因此，数值模拟成为研究水生动 速水翼，该水翼的形状和初始位置见图 1。我们 物游动的有效手段。 取坐标系原点为水翼的头部端点，x1为水平方向， 水生动物的游动是一个典型的动边界问题。 x2 为垂直方向。我们把水翼放在长宽各为 7.5L和 通常采用的贴体坐标和动网格方法，这种方法需 2.8L 的矩形区域，攻角为零度，其前缘端点的位 要跟踪移动边界然后生成贴体网格，大大增加了 置到c1 的距离为 1.5L，到c3 的距离为 1.4L，这 数值模拟的复杂性。文[5][6]分别用有限元方法数 里L表示水翼的弦长。水平来流的速度为U0 。 值模拟了水翼摆动的绕流，比较了不同Re数时水 翼的推进效率。本文采用浸入边界方法[7][8][9][10] 模拟水翼的摆动绕流。在浸入边界方法中，不需 1）国家自然科学基金资助项目和 和 973 项目非线性科学（2007CB814800 ）。 2 ） Email ：zhangx@lnm.imech.ac.cn 1.2 流场的控制方程 二 维 不 可 压 缩 流 场 的 控 制 方 程 为 Navier-Stokes 方程。在直角坐标系下，它的形式 为