粘性流体和刚性多联体结构相互作用的数值模拟.pdfVIP

数学、力学、物理学、高新技术研究进展——2008(12)卷 粘性流体与刚性多联体结构相互作用的数值模拟 白 洁‘张 星 ， 中国科学院力学研究所非线性力学重点实验室，北京，100190 摘要Childress最近在实验中发现低雷诺数下刚性体能够在振荡来流中悬停。本文主要通 过数值模拟的方法去研究这种刚性多联体结构拍动产生悬停的机制。这种多联体结 构由两个刚性平板铰支连接，在静止流场中作上下周期运动。流场的计算采用浸入 边界算法。多联体的运动采用求解欧拉．拉格朗日方程的方法。从数值结果中发现 多联体结构拍动可以产生升力，并且产生的升力的大小与振幅和频率有一定的关 系。 关键宇刚性多联体，浸入边界方法，弹性系数，欧拉啦格朗日方程 1．引言 自20世纪90年代以来，人们就致力于研究拍翼非定常空气动力学，探索昆虫在低 雷诺数下的空气动力学机制。近10年来，前人在活体观测，模型翼拍动试验和数值模拟 三个方面做了很多工作，力图揭开高升力之谜。Sun[1l等按照果蝇悬停状态拍翼模型的外 形及运动学数据，做了三维、非定常数值模拟研究，给出了拍动周期中变化的气动力以 及尾涡结构的流动细节，并提出了产生非定常升力的三个作用机制。吴江浩12J通过数值 计算比较了第一周期(不存在尾迹影响)和第四周期(存在尾迹影响)的气动力大小，得出了 “尾迹捕获”这一高升力机制一般并不存在的结论。继Sun对果蝇的悬停状态的研究之后， 近来研究昆虫和鸟类通过扑动翅膀而产生悬停机制的人越来越多。Stephenl3，4J等人通过 实验的研究，测量悬停状态与振幅和频率之间的关系。 2．数值方法 本文中模拟的是如图l所示的多联体结构，由第一段和第二段两个单元组成，其中 由一个节点将两个单元连接起来，在两段单元之间存在一个角弹簧，弹簧系数如图所示 为K，不考虑阻尼作用，将该结构放人不可压缩的黏性流中，给定固体的铰接点在竖直 方向为y=A·sin(2·兀·．厂幸f)的正弦谐波运动。采用浸入边界算法实现流体与固体的耦 合，通过计算流体网格中固体所占的体积分数来计算固体所受的力和力矩，再将固体所 受的力和力矩代人固体本构关系。固体计算采用与【5】类似的‘‘N段摆”模型，根据欧拉．拉 格朗日方程，运用四阶定步长龙格库塔算法16J，对该问题进行求解，具体方法请参考文 献【7·引。求解区域为矩形，采用均匀网格。流场初始速度为零，计算区域的四个边界都是 无滑移边界条件。 角度鼠，第二段的角度幺，以及连接点在】，方向的运动。其中岔，岛分别表示第一段， 粘性流体与刚性多联体结构相互作用的数值模拟 281 第二段与水平x轴负方向的夹角。 l F 图l多联体结构拍动的结构图 3．计算结果与分析 为了研究多联体拍动给升力带来的影响，采用了一个刚性模型进行对比，在两段多 联体连接处，采用弹性系数很大的结构，抑制两段翼的拍动。这样两段翼上的所有的点 的运动都与连接点指定的运动相同，即与刚体运动相似。 图2 从图2可以发现对于振荡的刚体运动，流体作用在每一段上的向上的力与向下的力 的大小是相等的。因而在一个时间周期内进行积分，可以发现不产生向上的升力。然而 当多联体的弹性系数KS=50的时候，它能够随着连接点的运动而产生拍动时，将能产生 升力，即随着时间的变化，流体对多联体产生向上方向的力大于向下方向的力。时间平 均的升力为20(特征长度为一段刚性板长，特征速度为振动速度幅值)。 282 数学、力学、物理学、高新技术研究进展——2008(12)卷 lV 图3多联体向上运动与向下运动的角度变化关系图 图3描述的是多联体向上和向下的运动过程，当多联体向上运动时，受到流体给的 向下的力，从而使两翼向下拍动，产生升力，拍动角度(岛一岛一岛)减小。当多联体向 下运动时，会受到流体给的向上的推力，从而使两翼张开，从而产生负升力，