悬停昆虫的空气动力学性能在很大程度上是由于其翅膀顶部存在稳定附着的前缘涡流(LEV)。虽然LEV已经在真实的、物理建模的和模拟的昆虫上可视化,但对其稳定性的物理机制还知之甚少。为了从根本上了解扑翼飞翼的LEV稳定性,我们在附于机翼表面的旋转参照系中表达了Navier–Stokes方程。利用这些方程,我们表明扑翼上的LEV动力学受三个项控制:角加速度、向心加速度和科里奥利加速度。我们对悬停条件的分析表明,角加速度与无量纲冲程振幅的倒数成正比,而科里奥利加速度和向心加速度与罗斯比数的倒数呈正比。使用一个动态缩放的果蝇翅膀扑动机器人模型,系统地改变这些无量纲数,我们确定了三个加速度中哪一个影响LEV稳定性。我们的力测量和流动可视化表明,LEV通过“准静态”向心加速度和科里奥利加速度来稳定,这些加速度在低罗斯比数下出现,并由机翼的螺旋桨式扫掠产生。相反,扑翼前后运动产生的非定常角加速度似乎在LEV的稳定附着中不起作用。然而,角加速度对LEV完整性至关重要,因为我们发现它可以调节LEV螺旋破裂,即高雷诺数效应。我们的分析和实验进一步表明,LEV稳定性的机制不依赖于雷诺数,至少在与昆虫飞行最相关的范围内(100<重新<14,000). LEV是稳定的,即使在爆炸时也会继续增强力量。这些以及雷诺数更高的螺旋桨和风力涡轮机的类似发现表明,即使是大型飞行动物也可能在慢速悬停飞行、起飞或着陆期间利用基于LEV的力量增强。我们根据300多种昆虫、鸟类、蝙蝠、自转种子和鱼类胸鳍的单翼长宽比计算出了罗斯比数。我们发现,平均而言,翅膀和鳍的罗斯比数与苍蝇接近(Ro公司=3). 从理论上讲,这些动物中的许多应该能够产生稳定的LEV,这一预测得到了一些昆虫、一只蝙蝠、一只鸟和一条鱼的最新发现的支持。这表明,通过稳定附着的(前缘)涡流增强力可以代表在相当大的尺寸范围内产生高流体力的收敛解决方案。

稳定前缘涡流(LEV)的存在是昆虫翅膀在悬停飞行中出乎意料的高性能的一个关键特征(麦克斯沃西,1979年Ellington等人,1996年Dickinson等人,1999年Srygley和Thomas,2002年). 而在昆虫翅膀的平移2D模型上,LEV在经过几个弦长的移动后脱落(迪金森和戈茨,1993年迪金森,1994Miller和Peskin,2004年Lentink等人,2008年),它仍然稳定地附着在围绕其基座旋转的3D模型机翼上(Dickinson等人,1999年Usherwood和Ellington,2002年Birch等人,2004年). 范登伯格和艾灵顿(Van Den Berg和Ellington,1997年)注意,由天蛾翅膀的机械模型产生的螺旋LEV与三角翼和后掠翼产生的螺旋LEV非常相似(Ellington等人,1996年Van Den Berg和Ellington,1997年). 这种后掠翼上的螺旋LEV通过机翼后掠翼诱导的展向流动来稳定,这表明展向流动对昆虫翅膀上LEV的稳定性同样至关重要(Ellington等人,1996年Van Den Berg和Ellington,1997年). 具体地说,昆虫翅膀上的LEV的生长可以通过LEV核心的展向流动、与扑翼速度梯度相关联的动态压力梯度、边界层中的“离心”加速度或螺旋涡线的诱导速度场来稳定(Ellington等人,1996年). 另一个假设是,由低展弦比昆虫翅膀的强大叶尖涡诱导的气流通过大大降低有效迎角来稳定LEV(Birch和Dickinson,2001年). 使用各种挡板阻止翼展方向流动的尝试对LEV强度或稳定性几乎没有影响(Birch和Dickinson,2001年)但这些实验并没有明确地确定LEV稳定性的独特解释。

在本研究中,我们首先通过实验表明,无论是后掠机翼的类比还是由叶尖涡引起的诱导流动都不能完全解释飞行机翼的LEV稳定性。基于旋转昆虫翅膀也能稳定LEV的概念(Dickinson等人,1999年Usherwood和Ellington,2002年Birch等人,2004年)然后,我们使用坐标变换将参考系附加到扑翼的表面,将纳维–斯托克斯方程应用于扑翼[这些方程已在随附的论文中导出(Lentink和Dickinson,2009年)]. 该分析显示了机翼运动学如何可能稳定旋转机翼上的LEV。基于两个控制无量纲数,Rossby数(Ro公司),用于测量科里奥利加速度和无量纲冲程振幅(A类*),这是流动不稳定性的度量,然后我们进行了一组实验,以确定这些无量纲数中是否有任何一个通过相应的加速度调节LEV稳定性。使用流动可视化和力测量,我们显示Ro公司,而不是A类*,似乎可以解释LEV的稳定性。然后,我们将我们的理论和实验结果与自然界和技术中其他翅膀和鳍的文献进行比较。

前面已经描述过动态缩放昆虫翅膀的基本方法(Dickson和Dickinson,2004年). 我们构建了一个模型黑腹果蝇2.0 mm厚透明丙烯酸板制成的机翼,带(单)翼展b条0.187 m和表面积S公司0.0167米2.平均弦宽c(c)定义为S公司/b条这里和其他地方我们的机翼参数(b条S公司c(c))指的是单翼,而不是双翼。力传感器将机器人手臂与机翼底座连接,从而形成翼尖半径R(右)该机翼连接到一个与三自由度执行器串联的力传感器上,执行器连接到浸没在1 m×1 m×2 m油箱中的移动机器人臂上,油箱中装满了油或水。雷诺数计算为:重新=()/ν、 其中c(c)是平均弦长,U型回转半径处的平均速度R(右)(艾灵顿,1984年)和ν运动粘度。油箱中装满了厚厚的矿物油(密度ρ=840 kg m–3; ν=140×10–62–1)以获得重新=110,稀矿物油(ρ=830 kg m–3; ν=11.0×10–62–1)以获得重新=1400和水(ρ=998 kg m–3; ν=1.004×10–62–1)以获得重新=14,000.

在我们的实验中,我们使用了以下扑翼运动模式:冲程位置的正弦运动和攻角运动的平滑梯形(Dickinson等人,1999年). 我们基于六只缓慢悬停果蝇的自由飞行运动学实验中的70°冲程振幅(Fry等人,2003年). 昆虫翅膀运动学具有先进的翅膀旋转、U形冲程平面和攻角运动中的小驼峰等特征,我们在这里忽略了这些特征,尽管已知其中一些特征可以提高昆虫的飞行性能(例如。Dickinson等人,1999年). 注意,我们使用数学意义上的振幅,它等于Ellington定义的总机翼振幅的一半(艾灵顿,1984年). 在每次扑动试验中,机器人扑动了六个完整的周期[平均扑动频率为0.22Hz重新=110,0.23赫兹重新=1400和0.20 Hz重新=14,000]. 几何攻角振幅α0从0°到90°变化,步长为4.5°[有关扑翼运动学的定义,请参见Sane和Dickinson,以及Lentink和Dicknson(Sane和Dickinson,2001年Lentink和Dickinson,2008年)]. 单向旋转和平移机翼运动学由恒定速度冲程组成,冲程开始和结束时具有恒定加速度和减速度。旋转和平移机翼的加速持续时间均为冲程的10%。与扑翼试验一样,α0以4.5°的步长从0°变为90°。旋转机翼扫过320°的弧线;计算平移翼的移动距离,使其在回转半径处移动到与旋转翼相似的距离。

我们生成了一系列Rossby数字(Ro公司)对于特定的雷诺数(重新=1400),将机器人手臂伸长1.27和1.53倍,增加了Ro公司回转半径分别为3.6和4.4。单向和双向平移机翼运动学(Ro公司=∞)。The stroke amplitudes for theRo公司在执行机构盘面积为(斯蒂芬妮夫斯基和凯斯,1984年)在划水过程中被机翼扫过,精度小于0.1%。平均雷诺数R(右)单向运动变化<0.5%,往复运动变化<5%。

使用气泡的流动可视化

我们在机翼前缘(相距约25 mm)和后缘(相隔约30 mm)向油中释放了小气泡(在重新=110和1400)。使用一根2 mm的细管将空气输送至前缘和后缘,该细管与2 mm厚机翼的边缘平齐,从而使气流受到管道的干扰最小。我们用昆虫别针在管子上打了个洞,上面画了一个白点。通过这些孔释放的气泡向上上升,因为它们不是中性浮力,因此不会形成完美的条纹。虽然我们能够通过用最小的昆虫针(针对果蝇)刺穿管子来最小化气泡大小,从而最小化向上速度,但我们无法准确控制气泡大小,导致气泡大小不同,上升速度也不同。最后,我们注意到气泡最初执行轨迹跟踪,与单个气泡相比,一组气泡的上升速度更高。然而,通常情况下,流动可视化与之前执行的粒子图像测速测量很好地匹配。因此,该方法特别适用于简单的3D LEV可视化,因为气泡将被吸入具有低压核心的强旋涡中,从而形成紧密的螺旋,而较弱且较宽的旋涡将形成较宽的气泡螺旋。此外,由于气泡密度低,在离心载荷作用下,气泡将优先向内(从翼尖到基座)驱动。因此,向外流动的气泡可靠地表示向外流动。

我们用单向或往复冲程运动学对α的平移和旋转机翼周围的流动进行了可视化0=0,18、36、45、54、72和90°对于图像记录,我们使用数字单色Basler相机:656×491,100帧采样–1对于扑翼,我们获得了六个扑翼周期的可视化结果。我们排除了第一个周期,并在重新=1400.破裂时刻被定义为我们首次注意到漩涡中出现白色多云螺旋状气泡堆积的时刻(这是一个定性定义)。使用Photoshop(8.0,Adobe)的自动对比功能增强了图像,并使用VirtualDub(1.5.10,网址:www.cole2k.net).

升力和阻力测量

如前所述,使用定制的力传感器测量作用在机翼上的升力和阻力(Dickinson等人,1999年Birch和Dickinson,2001年Birch等人,2004年Dickson和Dickinson,2004年). 对于后处理,我们以300 Hz的频率对测量值进行下采样,这仍然大约是摆动频率的1400倍。使用零相位延迟低通4极数字巴特沃斯滤波器对力信号进行离线滤波。确定截止频率,使其与0.3弦长的平均移动距离相对应R(右)该距离比已知LEV脱落的距离低至少10倍(迪金森,1994迪金森和戈茨,1993年). 具有往复运动的机翼(六个襟翼)的力在四个周期(第二到第五个周期)内平均,而对于具有单向运动的机翼,其平均值在冲程周期的70%到90%之间(排除启动和停止瞬态)。最终值是三次试验的平均值,但单向平移情况除外重新=110和后掠机翼极N个=1.

为了在实验之间进行公平的比较,我们计算了机翼在等动压力下产生力的效率,这是工程和动物飞行文献中的标准方法。这些力系数是根据平均动压力计算的

\(1{/}2{\rho}\上划线{V{\mathrm{g}}^{2}}\)
我们用叶片单元法计算(艾灵顿,1984年Dickinson等人,1999年Sane和Dickinson,2001年Dickson和Dickinson,2004年). 升力系数计算如下
\(C_{\mathrm{L}}=2\overline{L}{/}{\rho}S\overline{V{\mathr{g}}^{2}})
而阻力系数计算为
\(C_{\mathrm{D}}=2\上横线{数字视频}_{\mathrm{g}}{/}{\rho}S(上划线{V{\mathr{g}{^{2}})^{3{/}2}\)
在哪儿V(V)是回转半径处的速度,ρ是密度S公司机翼表面积。根据这个定义,C类D类将平动运动学的阻力系数降为经典阻力系数。对于旋转机翼运动学,C类D类是平均阻力系数,等于局部阻力系数沿机翼半径不变化的情况下的功率系数。

现有假设的测试

使用动态缩放的飞翼(黑腹果蝇)(Dickson等人,1999年;Dickson和Dickinson,2004年)为了测量力和可视化流动,我们首先测试了单靠机翼后掠角和翼尖效应是否能够稳定以固定速度平移(但不旋转)的飞翼上的LEV。我们在大范围的攻角(0至90°)上系统地改变了0至60°的扫描角。结果,在重新110和1400,表明机翼后掠角无法稳定LEV(图1A; 补充材料电影1)。此外,零扫掠角下的结果表明,叶尖涡的存在也不足以稳定LEV(图1B; 补充材料电影2)。重要的是要注意,当以恒定角速度围绕其基座旋转时,完全相同的机翼会产生稳定的LEV和升高的力(图1C补充材料电影3)其他人发现的(Dickinson等人,1999年Usherwood和Ellington,2002年Birch等人,2004年). 在雷诺数110和1400时,平移后掠翼不仅迅速降低了LEV,而且实际产生的升力比未后掠翼小(图2). 因此,在后掠翼飞机上稳定LEV的机制的严格类比似乎不适用于昆虫翅膀。此外,无后掠平动机翼的脱落LEV和低性能表明,单靠翼尖效应无法产生稳定的LEV,至少在我们的模型飞翼的展弦比下是如此。然而,翼尖效应似乎可以解释低展弦比机翼的LEV稳定性(1936年冬季林盖特,2007年).

这些初步实验促使我们明确研究旋转、螺旋桨状运动在LEV稳定性中的作用。当悬停时,大多数昆虫在大致水平的飞行平面上前后摆动翅膀。每次反转时,机翼都会迅速翻转并改变方向,在此期间,力和气流高度不稳定。在一个冲程开始时产生的LEV脱落,当机翼翻转并反转方向时,形成一个新的反向LEV(Poelma等人,2006年). 然而,在每半冲程(即上冲程和下冲程)中,机翼的运动是“螺旋桨式”的,因为机翼以大致恒定的攻角围绕其基座旋转(Usherwood和Ellington,2002年). 我们的工作假设是,这种旋转类螺旋桨运动固有的流体动力学特性是LEV稳定性的原因。我们首先确定了理论上可能导致LEV稳定性的基于旋转的流体加速度的完整列表,然后对其进行了实验测试,从而探索了这一假设。

图1。

前缘涡(LEV)从平移模型昆虫翅膀上脱落,无论其扫掠角度如何,而当翅膀旋转时,前缘涡仍保持稳定附着。LEV以雷诺数可视化重新=110和1400,机翼前缘和后缘释放小气泡。距离()机翼的行程以弦长表示c(c)在回转半径处R(右)(A)后掠机翼上的LEV不稳定,如图所示,在α=36°和重新=1400.(B)在没有后掠角的情况下,同一机翼上的LEV也不稳定。(C) 旋转同一个机翼会产生一个保持稳定连接的LEV,为了清晰起见,如图所示重新=110和=8c,其中s是以弦长表示的行驶距离。

图1。

前缘涡(LEV)从平移模型昆虫翅膀上脱落,无论其扫掠角度如何,而当翅膀旋转时,前缘涡保持稳定附着。LEV在雷诺数下可视化重新=110和1400,机翼前缘和后缘释放小气泡。距离()机翼的行程以弦长表示c(c)在回转半径处R(右)(A)后掠机翼上的LEV不稳定,如图所示,在α=36°和重新=1400.(B)在没有后掠角的情况下,同一机翼上的LEV也不稳定。(C) 旋转同一个机翼可使LEV保持稳定连接,如图所示重新=110和=8c,其中s是以弦长表示的行驶距离。

扑翼附近流动的Navier–Stokes方程

我们开发了一个简单的理论框架,以确定悬停状态下可能影响单向(类螺旋桨)和往复(类昆虫)运动机翼附近流动的无量纲数。我们分析的一个关键特征是,它能够适应从纯旋转运动到纯平移运动的连续冲程运动学范围。分析是连续的,因为平移代表了机翼绕无限半径以无穷小角度旋转的极限情况(图3A); 从这个意义上说,平移机翼围绕无限转弯半径执行悬停飞行。为了在实验之间进行一致的比较,三个关键条件得到了很好的近似值。首先,旋转翼扫过的面积保持不变,从而保持恒定的弗劳德效率(斯蒂芬妮夫斯基和凯斯,1984年). 第二,无量纲冲程振幅(A类*)在回转半径处(艾灵顿,1984年)保持恒定,以确保翼-尾流相互作用(Birch和Dickinson,2003年)都是类似的。最后,重新在回转半径处也保持不变(图3A). 这种分析最方便的理论框架是纳维-斯托克斯方程的无量纲形式,用固定在悬停飞行中旋转机翼上的非惯性参考系表示[推导见Lentink和Dickinson(Lentink和Dickinson,2008年)](另请参见Vanyo,1993年Greitzer等人,2004年). 惯性系中作用在流体“粒子”上的净粘性和压力引起的无量纲流体加速度惰性的与旋转框架中的相关本地由(巴鲁,1999年):
\[\\mathbf{a}_{\mathrm{惰性}}=\mathbf{a}_{\mathrm{loc}}+(\mathbf{a}_{\mathrm{ang}}+\mathbf{a}_{\mathrm{cen}}+\mathbf{a}_{\mathrm{Cor}}),\]
(1)
哪里:
\[\\mathbf{a}_{\mathrm{ang}}=1{/}A^{\ast}}{\cdot}{\dot{{\Omega}}}{\times}\mathbf{r},\]
(2)
\[\\mathbf{a}_{\mathrm{cen}}=1{/}Ro{\cdot}{\Omega}{\times}({\Omega}{\times}\mathbf{r})\]
(3)
\[\\mathbf{a}_{\mathrm{Cor}}=1{/}Ro{\cdot}2{\Omega}{\times}\mathbf{u}_{\mathrm{loc}}.\]
(4)
这里,Ω是角速度
\({\dot{{\Omega}}\)
是旋转框架的角加速度,以及第页u个本地分别是旋转框架中流体体积的位置和速度(图3B). 角加速度与A类*,这是无量纲冲程振幅的度量:
\[\A^{{\ast}}={\Phi}_{0}右{/}c,\]
(5)
其中Φ0是以弧度表示的冲程振幅,R(右)是机翼长度,并且c(c)是平均弦长。该术语将振幅表示为移动的弦长数。
图2。

扫掠不会增加平移果蝇翅膀产生的升力重新=110和1400。(A) 电梯(C类L(左))–拖动(C类D类)显示的极轴重新=110(B)和重新=1400(C)表示无量纲升力(L(左))和弦向阻力(D类和弦)以4.5°的步长将攻角(α)从0°变为90°获得的力。我们测试了从0°(用黑色表示)到60°(用红色表示)的机翼扫掠,扫掠增量为10°。

图2。

扫掠不会增加平移果蝇翅膀产生的升力重新=110和1400。(A) 电梯(C类L(左))–拖动(C类D类)显示的极轴重新=110(B)和重新=1400(C)表示无量纲升力(L(左))和弦向阻力(D类和弦)以4.5°的步长将攻角(α)从0°变为90°获得的力。我们测试了从0°(用黑色表示)到60°(用红色表示)的机翼扫掠,扫掠增量为10°。

图3。

(A) 用于分析平移和旋转机翼上流体加速度的框架。执行机构盘的面积在模型内是恒定的。R(右)是回转半径,在以下位置移动的弦长数R(右)在全冲程期间,以及V(V)(第页)是沿机翼半径的速度分布第页(B)在附翼框架中,由于机翼的冲程运动学,靠近机翼的流体经历三次加速度:角加速度英国,向心加速度中心和科里奥利加速度科尔。请注意u个本地是机翼固定框架中的局部速度,Ω(打、击等的)一下是由于行程引起的角速度,

\({\dot{{\Omega}}\)
冲程是由于冲程和第页是流体颗粒在旋转框架中的位置。

图3。

(A) 用于分析平移和旋转机翼上流体加速度的框架。执行机构盘的面积在模型内是恒定的。R(右)是回转半径,在以下位置移动的弦长数R(右)在全冲程期间,以及V(V)(第页)是沿机翼半径的速度分布第页(B)在附翼框架中,由于机翼的冲程运动学,靠近机翼的流体经历三次加速度:角加速度英国,向心加速度中心和科里奥利加速度科尔。请注意u个本地是机翼固定框架中的局部速度,Ω(打、击等的)一下是由冲程引起的角速度,

\({\dot{{\Omega}}\)
冲程是由于冲程和第页是流体颗粒在旋转框架中的位置。

中括号内的三个术语等式1是有角度的(英国),向心(中心)和科里奥利(科尔)加速。从物理上讲,这三个加速度是由机翼的运动学产生的,并在靠近机翼表面的空气中产生,而机翼表面既不能流过,也不能相对于机翼滑动(Vanyo,1993年Greitzer等人,2004年). 我们举例说明了由冲程引起的主导角速度引起的三个旋转加速度(Lentink和Dickinson,2008年)英寸图3B第一个组件,英国,是机翼围绕其基座的角加速度的表现形式,这导致局部弦向加速度(图3B). 在以恒定角速度单向旋转的机翼上没有这个术语,但如果角速度发生变化,如往复运动,这个术语就会出现(图3A). 第二任期,中心,表示向心加速度,其沿翼展方向朝向机翼底部(图3B). 第三学期,科尔,表示科里奥利加速度;它的方向取决于局部流体速度的方向u个本地(图3B). 向心加速度和科里奥利加速度(中心科尔)是“准静态”的,因为它们取决于机翼角速度Ω的瞬时值。这与角加速度相反(英国),这取决于角速度的变化

\({\dot{{\Omega}}\)
注意,我们考虑加速度(方程式24图3B)而不是指向相反方向的类似“虚拟力”(Vanyo,1993年Greitzer等人,2004年).

三个加速度项的大小英国中心科尔根据流体的对流加速度(在局部框架中)进行缩放,从而得出等式中的三个独立无量纲数24在悬停飞行的特殊情况下,向心加速度和科里奥利加速度与罗斯比数成反比Ro公司(罗斯比,1936年Lentink和Dickinson,2008年). 从现在开始,我们将使用无量纲长度刻度A类*Ro公司量化角加速度、向心加速度和科里奥利加速度。对于旋转机翼,Ro公司等于R(右)/c(c),回转半径除以平均弦长的比率。Ro公司对于平移的机翼来说是无限的,因为旋转半径是无限大的(图3A). 可以方便地计算与翼尖半径相关的罗斯比数,R(右),而不是回转半径:
\[\Ro=R{/}c,\]
(6)
因为这个值相当于单个机翼的展弦比应收账并且很容易从生物文献中提取(有关详细信息,请参阅Lentink和Dickinson,2008年). 的典型值Ro公司(基于翼尖半径)3附近昆虫翅膀簇(图12),这立即表明旋转加速度可能很大(注R(右)/c(c)≈1.5,因为R(右)≈0.5R(右)用于昆虫)。对于往复式机翼,A类*等于A类/c(c); 冲程振幅比A类平均弦长c(c)再次,为了简单起见,我们考虑机翼尖端的冲程振幅,而不是回转半径。注意,对于单向旋转机翼A类,因此A类*,是无限的。

什么是相对重要性英国中心科尔昆虫翅膀?在悬停飞行中A类*Ro公司(等式的商56)是Φ0,其中Φ0是定义机翼往复运动的调和函数的振幅(弧度)(参见A类*Ro公司在等式中56). Φ0昆虫的范围约为0.6至1.5(艾灵顿,1984年). 因此A类*Ro公司在昆虫中,排序为一,这表明英国中心科尔也有类似的震级。这不仅适用于昆虫,也适用于持续或短暂悬停条件下的大型动物。假设来回运动大致协调,可以粗略评估每个冲程开始、中间和结束时的旋转加速度,这对许多昆虫来说是一个合理的假设,从而获得进一步的见解(艾灵顿,1984年). 在冲程结束和开始时.a个中心科尔是最小的,因为Ω为零,而

\({\dot{{\Omega}}\)
因此英国最大值,按1缩放/A类*然而,这些条件对LEV稳定性的影响可能不大,因为在中风逆转期间,LEV脱落并重新形成(具有相反的意义)(例如,参见Poelma等人,2006年). 中游时,当LEV稳定性出现问题时,中心科尔最大值,按1缩放/Ro公司,而英国接近零。这个简单的分析表明,LEV的稳定性可能由旋转加速度介导中心科尔而不是不稳定的加速度英国。以下实验分析的主要目的是明确测试此理论预测。

LEV动力学对无量纲数的依赖性

我们对旋转体进行了一系列流动可视化和力测量(Ro公司=2.9)和翻译(Ro公司=∞)单向飞行的机翼(A类*=∞)和往复式(A类*=3.5)攻角振幅在0至90°之间的运动Ro公司A类*是自由飞行中缓慢盘旋的果蝇的代表(Fry等人,2003年)与许多昆虫的平均值接近。(注意:Ro公司=∞对应于1/Ro公司=0,即科里奥利加速度和向心加速度为零。同样,A类*=∞表示零角加速度。)

图4。

罗斯比数的影响Ro公司LEV稳定性。(A) 在旋转往复运动的机翼上显示的LEV(Ro公司=2.9,A类*=3.5)或平移往复运动(Ro公司=∞)。可视化是在=4c,在回转半径处,接近冲程末端(4.4c(c)). Re=1400时,LEV沿机翼中部爆发(b)(Ro公司=2.9),但保持稳定连接。(B) 顶部面板:LEV在重新=1400,行程中途,α=45°。底部面板:在中冲程α=18°冲程结束时,我们观察到爆裂LEV呈λ形;它在顶端附近分裂成两个“双涡”,其中红线以下的底部涡呈螺旋状破裂(b)。(C) 观察到稳定LEV的出现和往复旋转机翼冲程内螺旋破裂的开始(Ro公司=2.9)重新=1400.LEV在所有攻角下都是稳定的,并且在α>18°的情况下表现出螺旋式中冲程,其中起点用圆圈表示(直径大于五冲程计算的s.d.)。

图4。

罗斯比数的影响Ro公司LEV稳定性。(A) 在旋转往复运动的机翼上显示的LEV(Ro公司=2.9,A类*=3.5)或平移往复运动(Ro公司=∞)。可视化是在=4c,在回转半径处,接近冲程末端(4.4c(c)). Re=1400时,LEV沿机翼中部爆发(b)(Ro公司=2.9),但仍保持稳定附着。(B) 顶部面板:LEV在重新=1400,行程中途,α=45°。底部面板:在中冲程α=18°冲程结束时,我们观察到爆裂LEV呈λ形;它在顶端附近分裂成两个“双涡”,其中红线以下的底部涡呈螺旋状破裂(b)。(C) 观察到稳定LEV的出现和往复旋转机翼冲程内螺旋破裂的开始(Ro公司=2.9)重新=1400.LEV在所有攻角下都是稳定的,并且在α>18°的情况下表现出螺旋式中冲程,其中起点用圆圈表示(直径大于五冲程计算的s.d.)。

LEV在单向旋转机翼上是稳定的(Ro公司=2.9,A类*=∞;图1C补充材料电影3)如之前的研究所发现(Dickinson等人,1999年Usherwood和Ellington,2002年Birch等人,2004年)但不适用于单向平移机翼(Ro公司=∞,A类*=∞;图1B补充材料电影2). LEV在往复旋转机翼上也很稳定(Ro公司=2.9,A类*=3.5),但不适用于反向平移机翼(Ro公司=∞,A类*=3.5;图4A补充材料电影4和5). 这些结果在重新110(果蝇秤)和1400(家蝇或蜜蜂秤;补充材料电影6). 总之,往复运动(有限A类*)不足以稳定LEV。相反,LEV稳定性似乎只需要低Ro公司由旋转的螺旋桨状运动产生的。

虽然我们在两个位置都观察到稳定的LEV重新单向和双向旋转机翼上的110和1400,其流动结构并不相同重新数字,如前所述(Birch等人,2004年). 特别是,在重新1400在单向和往复运动下呈现螺旋破裂(图4B; 补充材料电影7)。螺旋涡的“爆发”是一种被认为是由核心流减速引起的现象(格林维尔,2002)并描述了三角翼在重新约1000人。对于单向旋转机翼,LEV在启动后以大于18°的攻角立即爆发,而对于往复式机翼,当机翼开始减速时,LEV几乎在机翼到达中点位置后爆发(图4C). 在我们的昆虫模型翅膀上,LEV在爆炸后保持一致,导致旋转流体的“湍流”体积,其相对于翅膀的位置保持稳定。螺旋爆破示例如所示图4A在中划和中划时,在迎角α=36°的情况下划臂结束图4B对于α=45°的中游[注:α=90°–α0用于拍打翅膀(Sane和Dickinson,2001年)]. 类似于Lu等人(Lu等人,2006年)在少数情况下,我们观察到双LEV结构,在较大的突发LEV前有一个较小的LEV(图4B).

图5。

流动画将我们的流可视化实验总结为重新Ro公司A类*.低Ro公司(2.9)产生稳定的LEV,A类*(3.5,∞)不修改此值和更高值重新(110至1400)会引起涡流爆发,但不会影响LEV相对于机翼的稳定连接。重新=110代表果蝇和重新=1400只家蝇。三角形表示单向冲程运动学,圆表示往复冲程运动学。Ro用颜色表示:黄色,Ro公司=2.9; 蓝色,Ro公司=∞,如中所用图6.

图5。

流动画将我们的流可视化实验总结为重新Ro公司A类*.低Ro公司(2.9)产生稳定的LEV,A类*(3.5,∞)不修改此值和更高值重新(110至1400)诱导涡流爆发,但不影响LEV相对于机翼的稳定附着。重新=110代表果蝇和重新=1400只家蝇。三角形表示单向冲程运动学,圆表示往复冲程运动学。Ro用颜色表示:黄色,Ro公司=2.9; 蓝色,Ro公司=∞,如中所用图6.

图6。

旋转机翼的升力和阻力系数(低Ro公司)比平移机翼产生的尺寸大(高Ro公司)在重新=110、1400和14000。电梯(C类L(左))–拖动(C类D类)所示的极坐标是通过将攻角从0°改变到90°来获得的。步长为4.5°。三角形表示单向冲程运动学,圆表示往复冲程运动学。Ro公司用颜色表示:黄色,Ro公司=2.9; 蓝色,Ro公司=∞; 白色,Ro公司=3.6和Ro公司=4.4. (A,B)在重新=1400,升力和阻力系数直接取决于Ro公司用于单向和往复冲程运动学。与单向平动机翼(a)的性能相比,B中平动往复机翼上不稳定LEV的力增加仍然很大。(C,D)低强度下的增力Ro公司从中找到重新=110至14000。

图6。

旋转机翼的升力和阻力系数(低Ro公司)比平移机翼产生的尺寸大(高Ro公司)在重新=110、1400和14000。电梯(C类L(左))–拖动(C类D类)所示极点是通过将攻角从0°变为90°来获得的。步长为4.5°。三角形表示单向冲程运动学,圆表示往复冲程运动学。Ro公司由颜色表示:黄色,Ro公司=2.9; 蓝色,Ro公司=∞; 白色,Ro公司=3.6和Ro公司=4.4. (A,B)在重新=1400,升力和阻力系数直接取决于Ro公司用于单向和往复冲程运动学。与单向平动机翼(a)的性能相比,B中平动往复机翼上不稳定LEV的力增加仍然很大。(C,D)低强度下的增力Ro公司从中找到重新=110至14000。

我们在中总结了基本的流可视化图5使用卡通来指示不同值下的基本流动结构A类*Ro公司重新它表明旋转的翅膀(Ro公司=2.9)中等紧凑和稳定的螺旋LEV,而LEV对于平动机翼来说是不稳定的(Ro公司=∞). 往复运动(A类*)不会改变LEV的稳定性,但在小冲程振幅下,它会保持LEV靠近机翼,因为它会及时脱落不稳定的LEV,并在每次冲程反转期间形成一个新的LEV(具有相反的意义)(A类*=3.5). 对于较高的冲程振幅,这不起作用,因为LEV在冲程反转之前脱落;单向平动是极限情况(A类*=∞). 雷诺数似乎不会影响所检查范围内的LEV稳定性。然而,雷诺数的增加确实会改变LEV的完整性,因为它会导致旋转机翼上的涡流爆发。在平移机翼时,我们没有观察到旋涡破裂,但在形成紧密旋涡并开始与机翼分离后,气流确实变得更加不稳定,这表明存在向湍流的过渡。

图7。

LEV的升力和阻力增加随雷诺数变化不大。(A) 单向平移机翼。(B) 单向旋转机翼。(C) 往复旋转机翼;星号代表果蝇的运动学。部队在重新攻角小于45°时,1400和14000重叠。低雷诺数(110)的主要影响是它会同时阻尼升力和阻力。

图7。

LEV的升力和阻力增加随雷诺数变化不大。(A) 单向平移机翼。(B) 单向旋转机翼。(C) 可相互旋转的机翼;星号代表果蝇的运动学。力的极值为重新攻角小于45°时,1400和14000重叠。低雷诺数(110)的主要影响是它会同时阻尼升力和阻力。

图8。

升阻极轴与性能极轴的关系图。减小最小阻力系数,这是为了增加雷诺数,可以提高性能(从黑色开始,到灰色结束,我们为零升力时的阻力系数绘制极点C类D0(数字0)=1.6, 0.8, 0.4, 0.2, 0.1). 圆圈表示最大功率因数和平方最大滑动数。

图8。

升阻极轴与性能极轴的关系图。减小最小阻力系数,这是为了增加雷诺数,可以提高性能(从黑色开始,到灰色结束,我们为零升力时的阻力系数绘制极点C类D0(数字0)=1.6, 0.8, 0.4, 0.2, 0.1). 圆圈表示最大功率因数和平方最大滑动数。

力系数对无量纲数的依赖性

评估LEV在不同条件下的疗效时的一个关键问题重新爆炸LEV是否仍会增加机翼产生的力。为了解决这个问题,我们测量了模型机翼在不同运动条件下产生的力。比较升力和阻力系数(C类D类C类L(左))生成于重新=110和1400,在上述四种运动条件下(Ro公司A类*在里面图5)表明爆裂不会导致力增加损失(图6A–C图7). 相反,力系数实际上在较高时升高重新,正如之前报道的那样(Birch等人,2004年). 低液位下稳定附着的LEV的存在Ro公司(旋转运动)在所有情况下都伴随着相对于Ro公司=∞情况(平动)。

测试LEV稳定性和力增强是否直接取决于Ro公司,我们改变了Ro公司对于单向和双向旋转机翼。这是通过实验实现的,方法是将机翼从其旋转轴向外延伸,以创建Ro公司值3.6和4.4(重新=1400)。力量增加随着增加而减少Ro公司(图6A、B)与LEV稳定在低水平的一般预测一致Ro公司然而,在这些条件下,我们在部队或视频记录中没有发现LEV脱落的证据。我们推测在低的允许范围内Ro公司旋转加速度的精确大小可能会影响平衡条件,并决定LEV的强度和功效。

观察到LEV即使在爆裂后也保持稳定,这鼓励我们测试LEV力的增加是否会扩大到更高的程度重新。在重新14000(蜂鸟等级)我们继续发现力量增强(图6D). 图7可以看出,扑翼、旋转和平移机翼的气动力极仅弱依赖于重新.重要重新然而,在零升力条件下,这些影响仍然作为最小阻力系数存在,C类D0(数字0),随着雷诺数的增加而减小。

扑翼、旋转和平移飞翼的功效

在过去,昆虫翅膀的性能分析主要集中在最大升力生产上,而LEV可以增加最大升力,但这种高升力机制的效率如何?为了评估空气动力学效能,我们使用两个相关指标构建了“性能极点”:滑翔数、,C类L(左)/C类D类和功率因数,C类3/2L(左)/C类D类(例如。Ruijgrok,1994年王,2008). 一定提升量所需的功率随着功率因数的增加而降低。图8显示了力极和性能极之间的关系。我们进一步说明了减少C类D0(数字0),从而提高性能;相应的性能最大值出现在低攻角时C类D0(数字0).功率因数最大(图8,圆形),攻角略高于滑翔次数最大的攻角(图8,方形)(Ruijgrok,1994年).

昆虫使用相互旋转的翅膀来产生升力,而直升机则使用旋转的叶片,飞机则简单地通过空气来转动翅膀。基于图7我们可以很容易地推断,扑动和旋转的翅膀产生的升力和阻力是平移的翅膀产生的升力和阻力的两倍<重新<14,000. 但哪种运动学产生升力最有效?对于重新=110–14000,我们发现旋转的飞翼比拍打的飞翼性能好100%,比平移的飞翼的性能好50%,以功率因数衡量(图9). 此外,1400<重新<14000我们发现平动翅膀的表现也优于扑翼,而扑翼飞行的翅膀在重新=110(果蝇秤)。我们的测量表明,果蝇实际上以接近最佳功率因数的运动学摆动翅膀;它们的机翼运动学导致悬停性能,与我们简化的机器人运动学的峰值性能非常匹配重新=110 (图9). 随着雷诺数的增加,与最小功率相对应的迎角减小,表明不太突出的LEV导致最大悬停效率。然而,该分析假设功率因数是衡量一般性能最合适的指标。如果最大升力受到限制,人们会得出截然不同的结论,因为最大升力大约发生在α=45°,与重新.

图9。

旋转飞翼的表现优于平移和拍打飞翼,达到110<重新<14000(A–C)。使用与图中相同的符号67最大功率因数表示悬停所需的最小空气动力。果蝇的运动学用星号表示;请注意,这大约对应于悬停所需的最小功率(最大功率因数)重新=110,在该值下运行。在更高的重新,果蝇运动学也很适合在最小功率下悬停,尽管功率因数略低于最佳值(B,C)。

图9。

110时,旋转的翅膀胜过平移和拍打的翅膀<重新<14000(A–C)。使用与图中相同的符号67最大功率因数表示悬停所需的最小空气动力。果蝇的运动学以星形表示;请注意,这大约对应于悬停所需的最小功率(最大功率因数)重新=110,在该值下运行。在更高的重新,果蝇运动学也很适合在最小功率下悬停,尽管功率因数略低于最佳值(B,C)。

使用动态缩放的机器人飞翼,我们可视化了气流,并测量了110度时机翼运动学范围产生的相应升力和阻力<重新<1400.我们测试了后掠机翼和旋转平移机翼的单向和往复运动。这使我们能够确定哪种运动学能够产生稳定的LEV、最大升力增强和最大的空气动力学性能。最终,该测试使我们能够确定哪个无量纲数和相应的旋转加速度最能预测LEV稳定性。

一阶Rossby数介导稳定LEV

我们的结果表明向心中心和科里奥利科尔加速度调节单向和双向旋转昆虫翅膀上LEV的稳定性(图5)这些旋转加速度与Rossby数成反比(Ro公司). 力量增加的减少Ro公司=2.9至3.6至4.4(图6A、B)表明LEV稳定性仅限于Ro公司一级或以下。然而,高展弦比旋转机翼在径向位置根部附近仍会经历显著的旋转加速度第页与弦长相比较小c(c)因此可以在本地支持稳定的LEV。这种影响可能是导致高展弦比风力涡轮机叶片轮毂附近的局部力高于预期的原因Ro公司=第页/c(c)小于3(Tangler,2004年)计算流体动力学(CFD)模拟证实了这一现象(Beom-Seok等人,2002年). 本地位于第页/c(c)<3科里奥利加速度和向心加速度较高,因为局部值与局部罗斯比数成反比第页/c(c),不考虑代表整个机翼的Rossby数(R(右)/c(c)),这对于高纵横比机翼来说很高。Lu等人的LEV可视化(Lu等人,2006年)关于大展弦比机翼的扑翼(Ro公司=1.3–10)表明,最显著的LEV(他们发现重新>640)实际上被限制在局部Ro公司大约低于3。这证实了这样一种观点,即细长的机翼可以在基地附近局部支撑稳定的LEV,而基地附近的局部Rossby数第页/c(c)<3,这会产生局部显著的旋转流加速度。

我们的理论预测和实验证实表明,旋转加速度介导了LEV的稳定性,但中心科尔实际参与?在风力涡轮机上,在重新订单106中心被归因于离心泵送(例如。林登堡,2004Vanyo,1993年Greitzer等人,2004年),从而导致轮毂附近的向外展向流动。在轮毂处,叶片经历所谓的“3D失速”,并产生升高的升力,导致局部升力系数远高于2(Tangler,2004年). 风力涡轮机轮毂区域的流动模式,其中Ro公司与昆虫翅膀的形状相似,与更远处的模式不同,在那里细长的叶片被称为“2D失速”(Beom-Seok等人,2002年Tangler,2004年林登堡,2004).

图10。

旋转圆盘上的边界层[在Vanyo之后(Vanyo,1993年)]. 旋转边界层通常被称为埃克曼层,而平衡径向流所需的朝向机翼的轴向流被称为埃克曼泵送。边界层速度剖面是自相似的,因为它以角速度乘以径向距离进行缩放。因此,所示速度剖面描述了整个速度场。压力场是旋转对称的,不随半径变化;它只随离圆盘的轴向距离而变化。因此,支撑科里奥利和离心加速度所需的力完全是由摩擦力引起的,摩擦力与速度梯度成正比。虽然在惯性(实验室)参考系中求解旋转机翼上的Navier–Stokes方程实际上是不可能的,因为需要进行表面跟踪,但在旋转圆盘的特殊情况下,最容易在惯性(实验)参考系下求解方程。原因是圆盘表面填满了一个无限平面,不需要跟踪,这大大简化了数学。注意,Ω是圆盘和粘附在圆盘上的流体颗粒的角速度,Ω′<Ω;第页,径向矢量;u个拉德,径向速度;中心,向心加速度;科尔,科里奥利加速度;(f)科尔,归一化科里奥利力;(f)中心,归一化向心力。

图10。

旋转圆盘上的边界层[在Vanyo之后(Vanyo,1993年)]. 旋转边界层通常称为埃克曼层,而朝向机翼的轴向流,需要平衡径向流,称为埃克曼抽吸。边界层速度剖面是自相似的,因为它以角速度乘以径向距离进行缩放。因此,所示速度剖面描述了整个速度场。压力场是旋转对称的,不随半径变化;它只随与阀盘的轴向距离而变化。因此,支撑科里奥利加速度和离心加速度所需的力仅由摩擦引起,摩擦与速度梯度成比例。虽然在惯性(实验室)参考系中求解旋转机翼上的Navier–Stokes方程实际上是不可能的,因为需要进行表面跟踪,但在旋转圆盘的特殊情况下,最容易在惯性(实验)参考系下求解方程。原因是圆盘表面填满了一个无限平面,不需要跟踪,这大大简化了数学。注意,Ω是圆盘和粘附在圆盘上的流体颗粒的角速度,Ω′<Ω;第页,径向矢量;u个拉德,径向速度;中心,向心加速度;科尔,科里奥利加速度;(f)科尔,归一化科里奥利力;(f)中心,归一化向心力。

对于简单的旋转圆盘,由圆盘向心加速度介导的向外径向流体流动是众所周知的(例如。Vanyo,1993年). 径向流被限制在称为Ekman层的边界区域(图10). 相应的埃克曼数,埃克=ν/Ωc(c)2,是该边界层中粘性力与科里奥利加速度之比的度量。对于旋转圆盘和机翼,埃克等于Ro公司/重新; 因此,人们可以独立地选择鄂、罗重新作为特征无量纲数的集合(Vanyo,1993年Greitzer等人,2004年). 第一盘,重新Ro公司,更适合于生物飞行分析,因为重新在文献中有更广泛的报道。在Ekman层中,旋转圆盘表面的流体具有与圆盘相同的角速度(由于表面无滑移条件),因此会经历向心加速度中心与旋转圆盘的相同,由径向摩擦力支撑(f)中心。在阀盘表面上方稍高的位置,由于摩擦,流体被拉向切向,但同时径向向外滑动。它之所以滑动,是因为没有足够大的摩擦力(来自径向速度梯度)来支持流体在圆盘表面获得的全部向心加速度。当流体粒子径向向外滑动时,会经历科里奥利加速度科尔在切向,因为它加速以匹配向外更高的切向速度并改变方向(它在被圆盘拉动时旋转)。科里奥利加速度由切向摩擦力支撑(f)科尔这是由轴向切向流动梯度引起的。在圆盘表面更高的地方,边界层就不复存在了。由于质量守恒,向外径向流动的流体必须供应“新鲜”流体,以轴向流向阀瓣。这个过程称为埃克曼抽水。

旋转圆盘顶部向外径向流动的区域也是昆虫翅膀顶部的一个显著特征(麦克斯沃西,1979年Ellington等人,1996年Birch和Dickinson,2001年Birch等人,2004年Poelma等人,2006年)靠近两个螺旋桨的轮毂(希梅尔斯坎普,1947年)和风力涡轮机(Tangler,2004年林登堡,2004)其中Ro公司局部偏低。事实上,已经计算了风力涡轮机叶片的类埃克曼边界层剖面(Dumitrescu和Cardos,2003年). 低速螺旋桨展向流动的重要历史证据Ro公司可以在Himmelskamp找到(希梅尔斯坎普,1947年),是风力涡轮机文献中的经典参考文献。图11,我们重现了他对低展弦比螺旋桨翼展方向流动的可视化。早期昆虫翅膀展向流动的可视化和测量表明,展向流动主要局限于LEV核心占据的流体区域(麦克斯沃西,1979年Ellington等人,1996年). 后来对单向和往复旋转模型果蝇翅膀进行的数字粒子图像测速(DPIV)测量表明,翼展方向的流动主要不限于LEV,而是一直延伸到后缘(Birch和Dickinson,2001年Birch等人,2004年Poelma等人,2006年). 与展向流动同时,存在平行于旋转轴的正交方向的显著流动,这表明存在埃克曼式泵送。甚至在离机翼表面很远的地方也能看到这种气流(Birch等人,2004年Poelma等人,2006年). 倾斜往复和单向旋转昆虫翅膀上朝向翅膀“轮毂”的径向流和正交流区域表明存在类埃克曼边界层,尽管由于“流动分离”,其空间范围要大得多(图12A)(Vanyo,1993年).

我们对旋转圆盘、风力涡轮机和螺旋桨的测量和观察表明,翼展方向的流动是由“离心泵送”造成的,因为靠近机翼表面的流体径向向外滑动。流体滑动是因为没有足够的摩擦力来支持完全的离心加速度中心机翼表面需要的流体(图12A、B). 在这些情况下,向心加速度也是由朝向机翼基座的净比力产生的((f)中心在里面图12B)它可能由压力分量和摩擦分量组成,摩擦分量会出现在旋转圆盘上(图10). 所有这些发现强烈表明,由于机翼在低高度旋转而产生的离心泵送Ro公司可以找到低和高重新–从昆虫翅膀到风力涡轮机叶片。

观察到分离的展向流动区域(包括LEV,图5)稳定地附在旋转机翼上意味着它不仅受到向心加速度的影响中心而且在弦向上持续加速。这种运动状态所需的加速度正是科里奥利项科尔(类似于旋转圆盘上的流动)。这意味着由中心相对于机翼稳定通过科尔,必须由指向行驶方向的净特定弦向力提供((f)科尔在里面图12B). 林登堡(林登堡,2004)表明如果科尔由作用于翼展方向流动区域的前后压力梯度支撑,这可能是风力涡轮机轮毂附近产生的提升力的原因。这种增力模型是否能定量地用于扑动昆虫翅膀尚不清楚,但可以在未来的研究中进行检验。

图11。

Ro公司<2和重新=280,000. (A) Himmelskamp制作的这些独特的基于簇的流可视化(希梅尔斯坎普,1947年)从未在期刊上发表,因此在此转载。预付款比率J型螺旋桨的速度计算为前进速度与翼尖速度之比;它的变化范围从0.124(几乎悬停状态)到0.459(向前飞行状态)。簇状物表明,在接近悬停状态的低推进比下,沿翼展方向的流动越来越强。此沿翼展方向的流动对应于升高的升力。(B) 为了完整起见,我们还根据Schlichting出版的径向站压力测量结果,复制了螺旋桨草图和测得的截面升力系数(Schlichting,1979年). 请注意R(右)是翼尖半径,第页局部半径,α迎角和C类截面升力系数。最大截面升力系数3远高于风洞中相同翼型产生的最大升力系数(虚线)。内侧部分,其中Ro公司是最低的,对应于最大升力,这也是由于螺旋桨叶片的扭曲。

图11。

Ro公司<2和重新=280,000. (A) Himmelskamp制作的这些独特的基于簇的流可视化(希梅尔斯坎普,1947年)从未在期刊上发表,因此在此转载。预付款比率J型螺旋桨的速度计算为前进速度与翼尖速度之比;它的变化范围从0.124(几乎悬停状态)到0.459(向前飞行状态)。簇状物表明,在接近悬停状态的低推进比下,沿翼展方向的流动越来越强。此沿翼展方向的流动对应于升高的升力。(B) 为了完整起见,我们还根据Schlichting出版的径向站压力测量结果,复制了螺旋桨草图和测得的截面升力系数(Schlichting,1979年). 请注意R(右)是翼尖半径,第页局部半径,α攻角和C类截面升力系数。最大截面升力系数3远高于风洞中相同翼型产生的最大升力系数(虚线)。内侧部分,其中Ro公司是最低的,对应于最大升力,这也是由于螺旋桨叶片的扭曲。

我们对向心加速度和科里奥利加速度的定性分析适用于整个展向流动区域,但可能不太详细。为了计算详细的加速度分布和相应的局部方向,需要速度场。孙和武(孙和武,2004)计算了单向旋转昆虫翅膀周围的速度场重新=480使用CFD。根据速度场,他们计算了“虚拟力”,该力指向与相应旋转加速度相反的方向(Vanyo,1993年Greitzer等人,2004年). 他们对攻角为40°的边界层流动中“科里奥利力”和“离心力”沿翼展方向分量的计算证实,科里奥利和离心加速度确实很重要。他们进一步发现,由于线性展向速度分布,径向压力梯度力甚至更大,并集中在LEV中,因此可以解释在重新=480。但压力梯度力在其他没有顺翼展流动的区域也很重要。鉴于此,我们注意到Birch等人的DPIV测量(Birch等人,2004年)在重新=110和1400表明,在无翼展方向流动的情况下,在飞翼下方、前方和上方的延伸区域,弦向速度显著。这一点很重要,因为在该区域,来自机翼表面的粘性效应可以忽略不计,因此粘性“离心泵送”机制无法工作。然而,该区域的压力梯度力是非零的,因为弦向流从机翼根部到叶尖径向变化,如果这种梯度是有效的,那么它应该再次驱动展向流动,但事实似乎并非如此。因此,翼展方向的流动和整个机翼周围流场中的压力梯度力在高低两种情况下似乎没有一一对应的关系重新。由于压力梯度,很可能沿翼展流动主要局限于LEV的核心。Aono和同事(Aono等人,2008年)进一步表明,果蝇的LEV存在较小的展向流动和压力梯度重新=134,但他们确实发现LEV后面有显著的顺翼展流动。对于天蛾重新=6300,他们确实在LEV中发现了强大的顺翼展流动和压力梯度。这些翼展方向流动的发现与我们的发现以及Birch等人的发现一致(Birch等人,2004年). Aono及其同事进一步指出,科里奥利和向心加速度可能是解释果蝇翅膀上LEV后面的展向流动的候选因素。

综合所有这些观察结果,包括我们的理论预测和实验验证,我们得出结论:(1)昆虫LEV核心的展向流动,当存在时,很可能是由展向压力梯度驱动的;(2) LEV后面扩展粘性流区域中的展向流动可以通过离心泵得到最好的解释,这与旋转圆盘、螺旋桨叶片和风力涡轮机叶片上的流动直接相似。

LEV完整性由以下因素调节重新A类*

与van den Berg和Ellington的旋涡破裂发现类似(范登伯格和艾灵顿,1997年)和Lu等人(Lu等人,2006年)我们发现飞翼上的LEV在重新1400.螺旋破裂是三角翼在临界点以上运行时出现的一种现象重新大约1000,被认为是由堆芯流量减速引起的(格林维尔,2002). 如果我们假设昆虫翅膀(包括LEV的核心)周围的气流在翅膀加速时加速,我们可以定性地理解为什么LEV在接近中游时开始破裂,而此时翅膀以及LEV的中心开始减速。这表明,尽管A类*不影响LEV的稳定性,相应的角加速度可以介导LEV的完整性。涡破裂可以解释在类似条件下先前定量流量测量中发现的LEV中不稳定的速度矢量(Birch等人,2004年). 最后,我们的测量表明,力系数不会因涡流破裂而降低(图56). 这表明,基于LEV的兵力增强对高重新数字效应。

新旧LEV稳定性假设的比较

我们的发现与之前导致后掠机翼类比和叶尖涡假设的发现有何关联?我们发现叶尖涡诱导的LEV稳定性似乎不适用于重新=110和1400(图12)尽管有证据表明,较高的雷诺数表明,当展弦比大约为一(或更少)时,平移短而粗的机翼确实会产生稳定的LEV(1936年冬季林盖特,2007年). 艾灵顿及其同事(Ellington等人,1996年)提出三种机制可以潜在地解释如何产生顺翼展流动。我们的理论分析、实验测试和文献调查表明,边界层中的“离心”加速度是低层大气的可能机制重新果蝇。处于高位重新,压力梯度力可以解释LEV堆芯的展向流动,而离心泵可以解释LEF堆芯后的展向流。在低雷诺数和高雷诺数下,展向流动区域显然不仅限于LEV核心,而且其在倾斜机翼上方的空间分布在很大程度上取决于雷诺数(Birch等人,2004年; Aona等人,2008年)。我们的分析进一步表明,科里奥利加速度在LEV的稳定附着中同样重要,因为该加速度是相对于旋转机翼稳定附着展向流的不可或缺的运动学条件。然而,所有这些实验都支持这样一种假设,即沿翼展流动平衡了前缘涡度的形成,并将其排入叶尖涡(麦克斯沃西,1979年Ellington等人,1996年). 的重要性Ro公司,而不是重新,在确定LEV稳定性时,进一步表明LEV后面大的“类埃克曼”分离流区域中的涡度输送在维持低平衡方面同样至关重要重新与高空LEV中发现的涡度输送相比重新我们从先前的一项研究中推断出这一点,该研究表明,地核内向外的螺旋流是重新从属的(Birch等人,2004年),而这些结果以及我们自己的结果表明,LEV的稳定性并非如此。此外,对于在同样低雷诺数下运行的形状相同的机翼,后掠和旋转机翼上的LEV之间的直接类比似乎不成立,因为我们没有观察到后掠果蝇翅膀的稳定LEV或力增强(图12). 在其他高展弦比扫掠鸟翅膀的实验中(雨燕)发现稳定的LEV(Videler等人,2004年Lentink等人,2007年),但12000无明显的力增加(升力系数低于1)<重新<77000(Lentink等人,2007年). 这表明,后掠机翼和旋转机翼上的LEV之间的直接类比在较高雷诺数下也不成立。

图12。

提出了稳定LEV的三维流动结构及其在飞行动物中通用的可行性。(A) 旋转果蝇翅膀局部流场示意图重新=110和1400,基于先前发布的粒子图像测速(PIV)数据(Birch等人,2004年)和图145翼展方向的气流存在于机翼顶部的整个区域,并将LEV中积累的涡度对流至叶尖涡。(B) 机翼的大展向流动区域由旋转加速度调节中心科尔这会形成一个类似埃克曼的边界层,类似于旋转圆盘和风力涡轮机叶片上的边界层。加速度由相应的归一化向心力平衡((f)中心)和科里奥利((f)科尔)流体中的力可以由压力和摩擦力组成。(C) 假定零推进比为体重函数的动物翅膀的罗斯比数(N个=总计319)。平均值接近果蝇(Ro公司=2.9):昆虫Ro公司=3.1(标准偏差=1.1,N个=98),蜂鸟Ro公司=3.7(标准偏差=0.3,N个=65),拍Ro公司=3.3(标准偏差=0.4,N个=39)和鸟类Ro公司=3.2(标准偏差=1.18,N个=117)。带有黑色轮廓的圆表示直接基于一个机翼纵横比的值。没有黑色轮廓的圆圈表示我们修正成对机翼和总机翼表面之间翼尖到翼尖距离的展弦比的值,以便我们获得可与翼尖Rossby数计算相比较的单翼展弦比。这相当于减去了机翼底座之间的机身宽度所产生的距离。关于这一点以及所有昆虫、鸟类、蜂鸟和蝙蝠翅膀形态参考的信息可以在附录1中找到。常数-雷诺数线的计算采用平均值C类L(左)1.5,代表我们的实验和Himmelskamp的雷诺数(希梅尔斯坎普,1947年); 详情见附录2。

图12。

提出了稳定LEV的三维流动结构及其在飞行动物中通用的可行性。(A) 旋转果蝇翅膀局部流场示意图重新=110和1400,基于先前发布的粒子图像测速(PIV)数据(Birch等人,2004年)和图145翼展方向的气流存在于机翼顶部的整个区域,并将LEV中积累的涡度对流至叶尖涡。(B) 机翼的大展向流动区域由旋转加速度调节中心科尔这会形成一个类似埃克曼的边界层,类似于旋转圆盘和风力涡轮机叶片上的边界层。加速度由相应的归一化向心力平衡((f)中心)和科里奥利((f)科尔)流体中的力可以由压力和摩擦力组成。(C) 假定零推进比为体重函数的动物翅膀的罗斯比数(N个=总计319)。平均值接近果蝇(Ro公司=2.9):昆虫Ro公司=3.1(标准偏差=1.1,N个=98),蜂鸟Ro公司=3.7(标准偏差=0.3,N个=65),拍Ro公司=3.3(标准偏差=0.4,N个=39)和鸟类Ro公司=3.2(标准偏差=1.18,N个=117). 带有黑色轮廓的圆表示直接基于一个机翼纵横比的值。没有黑色轮廓的圆圈表示我们修正成对机翼和总机翼表面之间翼尖到翼尖距离的展弦比的值,以便我们获得可与翼尖Rossby数计算相比较的单翼展弦比。这相当于减去了机翼底座之间的机身宽度所产生的距离。有关这方面的信息以及所有昆虫、鸟类、蜂鸟和蝙蝠翅膀形态学参考资料,请参阅附录1。常数-雷诺数线的计算采用平均值C类L(左)1.5,代表我们的实验和Himmelskamp的雷诺数(希梅尔斯坎普,1947年); 详情见附录2。

Rossby数和“准静态”升力理论之间的联系

这项研究支持早期的观点,即昆虫产生的空气动力可能被视为“准静态”,不包括中风逆转期间发生的并发症(Dickinson等人,1999年Sane和Dickinson,2001年Usherwood和Ellington,2002年). 准静态意味着,对于理解和预测空气动力,尤其是机翼升力,流速的瞬时值比其瞬时变化率更重要(Sane和Dickinson,2001年). 我们的理论框架支持这一观点,因为“准静态”旋转加速度中心科尔负责LEV稳定性,而非稳定角加速度英国不是。这并不是说不稳定性在昆虫飞行中并不重要。正如Wang等人发现的那样,平移往复情况导致LEV不稳定,仍然显著增加了力(Wang等人,2004年),但这个力小于旋转机翼的力(图6). 此外,与果蝇相比,LEV稳定性对中风幅度小得多的昆虫来说不那么重要。例如,空载悬停蜜蜂使用窄冲程振幅(Φ0=0.78 rad,而果蝇为1.2 rad),导致无量纲冲程振幅较低。最近的实验表明,蜜蜂在中风开始和结束时更加依赖于非定常力的增加(Altshuler等人,2005年),此时角加速度最大,旋转加速度最小。这意味着非定常升力机制,如附加质量效应和尾流捕获(Dickinson等人,1999年)与基于稳定LEV的“准静态”力相比,当无量纲冲程振幅减小时,变得越来越重要。

雷诺数高于14000时的升力增加

我们的实验(110<重新<14000)和其他人在螺旋桨上的实验(重新=280000)和风力涡轮机(重新订单106)表明对于一阶Rossby数,在惯性流区升力增加是连续的。艾灵顿和乌舍伍德(Ellington和Usherwood,2001年)然而,发现旋转模型鹰蛾翅膀的运行速度为10000<重新<50000未能产生高扬程。原因尚不清楚,他们的研究结果重新=10000与我们的重新=14000,鹌鹑翅膀在重新=26000,由Usherwood和Ellington提供(Usherwood和Ellington,2002年). 很可能,翼型在10000米的薄而锋利的鹰式机翼的明显升力危机中起着关键作用<重新<50,000. 施密茨(施密茨,1942年)在同一范围内通过风洞实验发现了翼型的这种现象。他发现薄而尖的翼型在低空时表现优于钝而厚的翼型重新反之亦然处于高位重新对于这些翼型,存在一个临界值重新其升力降低,阻力急剧增加。这是由于存在层流分离气泡和过渡到湍流边界层流动,而湍流边界层在中间层控制翼型性能重新这是相关的,因为Ellington和Usherwood测试的旋转机翼(Ellington和Usherwood,2001年)(中间位置重新)具有薄而尖锐的亚临界翼型,而螺旋桨(希梅尔斯坎普,1947年)和风力涡轮机叶片处于高位重新(Tangler,2004年)具有厚而钝的超临界翼型。艾灵顿和乌瑟伍德(Ellington和Usherwood,2001年)初步得出结论,LEV在旋转机翼上不稳定重新≥10000,因为缺乏稳定LEV的展向流量。Himmelskamp的展向流动显示(希梅尔斯坎普,1947年)在螺旋桨上重新=280,000 (图11)风力涡轮机叶片上稳定LEV的可视化(Beom-Seok等人,2002年),螺旋桨式风扇(例如。Simonich等人,1992年)和船用螺钉(例如。科尔文,1986年)在相似的攻角下运行Ro公司,但更高重新,反驳这个观点。Hubel及其同事的工作(Hubel,2006年)结果表明,模型鹅在中间雷诺数范围内扑动,确实可以在前飞过程中产生稳定的LEV。因此,我们得出的结论是,有明显证据表明,从旋转果蝇翅膀到风力涡轮机叶片,基于LEV的力增强可以持续存在Ro公司,但更多研究翼型对10000年中期LEV生成、稳定性和力增强的影响<重新需要<100000。

自然界和技术中旋转机翼的增力

总之,LEV稳定性和最大力增强的单一条件似乎足够低Ro公司因此,使用LEV增强力量在游泳和飞行动物中的分布可能比以前更广泛。图12C,我们显示了文献调查绘图的结果Ro公司用于悬停翅膀作为体重的函数。分布表明许多大型动物的翅膀高度足够低Ro公司创建稳定的LEV(注意Ro公司最重要的回转半径大约比Ro公司在翼梢处计算图12C). 这并不是说所有比昆虫大的动物都可以在巡航速度飞行时悬停或产生LEV,但没有空气动力学原因说明它们为什么不能在慢速悬停飞行或短距离起飞和降落时,在它们的推进比很小和Ro公司为一级(例如,参见图11,随着推进比的增加,展向流量减小)。LEV,低电平Ro公司蝙蝠在悬停飞行时翅膀上确实发现了这种情况(Muijres等人,2008年). 我们测量的提升力重新14000是一个适用于蜂鸟的值,与蜂鸟LEV的假定观测结果一致(Altshuler等人,2004年Warrick等人,2005年)和鹌鹑(Usherwood和Ellington,2002年)翅膀。最后,Hubel和他的同事在前飞中的一只扑翼鹅模型的底部附近发现了一个稳定的LEV(Hubel,2006年). 与风力涡轮机叶片一样,靠近机翼基座的局部Rossby数要低得多(另请参阅Lentink和Dickinson,2008年)因此,即使在相对较高的提前率下,也可以在本地支持稳定的LEV。

由于雷诺数和冲程幅度不是LEV稳定性的关键特征,我们认为LEV可能是一种高效的高升力机构,适用于大小缓慢悬停的动物。我们的实验表明,当雷诺数越来越高时,在较小攻角下产生的较小LEV的气动效率最高(图89). 然而,总的悬停效率不仅取决于空气动力学效率,还取决于驱动机翼的肌肉的效率。悬停时,动物的重量由垂直推力平衡,垂直推力与升力系数、扑翼频率平方和冲程振幅平方的乘积成正比。注意,为了提高肌肉效率,拍打频率被限制在一个狭窄的频带内(麦克马洪,1984年)冲程幅度限制在180°或更低,高最大升力系数明显有助于适应悬停期间平衡重量所需的高垂直推力和巡航期间所需的低得多的推力。

我们的理论框架同样很好地描述了空气和水。许多游泳动物的胸鳍类似于飞行动物的翅膀,不仅用于产生升力,还用于产生机动阻力。图6结果表明,旋转机翼(单向和往复式)不仅能为其表面积产生更多升力,而且在非常大的迎角(远远超过45°)下也能产生更多阻力。对于45°以上的迎角,附连的前缘和后缘涡(LEV和TEV)模型上的翅膀与最近在太阳鱼胸鳍上观察到的翅膀惊人地相似[兰黛和麦登在Bandyopadhyay等人(Bandyopadhyay等人,2008年)]. 基于鱼类胸鳍的单翼长宽比的罗斯比数通常很低。对于文献中描述的七个物种,我们发现Ro公司平均=2.5,s.d.=0.7(参考附录1),足够低的值Ro公司用于稳定的LEV。为了完整性,我们还估计了平均值Ro公司自动旋转种子的值为3.7(s.d.=1.14,N个=26,参考见附录1)。这表明稳定的LEV也可以解释这些植物结构的升力升高(Azuma和Yasuda,1989年). 我们最近使用3D DPIV对此进行了测试,并将在其他地方报告结果。在技术方面,我们设想微型飞行器可以更容易地模拟自然,通过简单地采用稳定LEV和力增强的唯一约束条件,产生稳定LEV,即以低罗斯比数连续旋转机翼,这比扑动同一机翼更有效。

从昆虫到鸟类的动物翅膀的罗斯比数

表A1 

表A1。

机翼数据来源

*

Puget Sound大学Slater自然历史博物馆的在线Wing收藏,http://www.ups.edu/x5662.xml

之间的关系Ro公司、质量和重新

悬停飞行中的雷诺数和罗斯比数(基于机翼半径;单翼翼展)如下所示:
\[\Re_{\mathrm{g}}=\frac{{\rho}\上划线{U}(U)_{\mathrm{g}}c}{{\mu}}(3.1)\]
(A1)
\[\Ro=AR_{\mathrm{s}}=\frac{b^{2}}{2S}=\frac{b}{2\上划线{c}},\]
(A2)
哪里重新是回转半径处的雷诺数。请注意b条是翼展,S公司机翼表面积和应收账单翼展弦比。悬停飞行中的力平衡要求以下关系成立:
\[\W=\overline{L}{\rightarrow}\ overline}{C_{\mathrm{L}}}1{/}2{\rho}\ overrine{U_{\mathrm{g}}^{2}}S=mg,\]
(A3)
哪里W公司是重量,时间平均升力,
\(上划线{C_{\mathrm{L}}}\)
时间平均升力系数,质量和引力常数。现在,我们用r.m.s.时间平均速度近似时间平均速度,从而得出:
\[\overline{U{\mathrm{g}}{\approx}\sqrt{\overline{U{\tathrm{C}}^{2}}=\sqrt{\frac{mg}{\overrine{C{\mathrm{L}}}1{/}2{\rho}S}}.\]
(A4)
组合公式A1类A2类A4(A4),我们获得了雷诺数作为悬停飞行总质量和翼尖Rossby数的函数:
\[\Re_{\mathrm{g}}=\frac{{\rho}\overline{U_{\methrm{c}}c}{\mu}}{\approx}\sqrt{\frac}mg}{\overline{c{\mathrm{L}}}}{\ times}\frac[{\rho}}{\\mu}^2}{\times}=\sqrt{\frac{mg}{\overline{c_{\mathrm{L}}}{\times}\frac}{\rho}}{{\mu}^{2}}{\times}\frac{1}{Ro}}.\]
(A5)
我们绘制了悬停动物翅膀的雷诺数等值线图12C假设从昆虫到鸟类的整个重量范围的时间平均升力系数为1.5。在这样做的过程中,我们假设动物正在使用LEV,并以高升力系数运行C类L(左)=1.5. 该近似值足以满足雷诺数的精确值相关性较小的情况。
图A1。

罗斯比的数字只需要对昆虫和鸟类进行校正。对于昆虫,我们使用了昆虫展开翅膀的照片,这导致了准确的修正(黄色圆圈没有轮廓图12C). 对于校正的鸟值(蓝色圆圈,无轮廓图12C)我们无法获得如此准确的照片,于是采取了以下行动。首先,我们获得了一组基于全翼展的展弦比数据(Tennekes,1997年). 随后,我们将这些值与我们也获得了准确的单翼展弦比值的值进行了比较(斯莱特博物馆,2005-2006年,美国华盛顿州塔科马市普吉特湾大学斯莱特自然历史博物馆的在线机翼收藏;http://www.ups.edu/x5662.xml)并确定了所需的修正系数C类基于差异。我们关联了校正因子C类整个机翼的展弦比除以2,应收账/2.该系数用于通过保守推断估算剩余鸟类翅膀的单翼展弦比第页信天翁的价值(应收账/2<7.5),接近1。

图A1。

Rossby数字只需针对昆虫和鸟类进行校正。对于昆虫,我们使用了昆虫展开翅膀的照片,这导致了准确的修正(黄色圆圈没有轮廓图12C). 对于校正的鸟值(蓝色圆圈,无轮廓图12C)我们无法获得如此准确的照片,于是采取了以下行动。首先,我们获得了一组基于全翼展的展弦比数据(Tennekes,1997年). 随后,我们将这些值与我们也获得单翼展弦比准确值的值进行了比较(斯莱特博物馆,2005-2006,斯莱特自然历史博物馆的在线机翼收藏,普吉特湾大学,华盛顿州塔科马,美国;http://www.ups.edu/x5662.xml)并确定了所需的修正系数C类基于差异。我们关联了校正系数C类整个机翼的展弦比除以2,应收账/2.该系数用于通过保守推断估算剩余鸟类翅膀的单翼展弦比第页信天翁的价值(应收账/2<7.5),接近1。

缩写词列表

     
  • α

    机翼迎角

  •  
  • α0

    机翼迎角振幅

  •  
  • Φ0

    机翼冲程振幅(总冲程振幅的一半Φ)

  •  
  • ν

    运动粘度

  •  
  • ρ

    流体密度

  •  
  • Ω

    旋转框架的角速度

  •  
  • Ω冲程

    机翼冲程引起的角速度

  •  
  • \({\dot{{\Omega}}\)

    旋转框架的角加速度

  •  
  • \({\dot{{\Omega}}}_{\mathrm{stroke}}\)

    机翼冲程引起的角加速度

  •  
  • Ω′

    从飞翼分离的流体的角速度

  •  
  • 英国

    角加速度

  •  
  • 中心

    向心加速度

  •  
  • 科尔

    科里奥利加速度

  •  
  • 惰性的

    相对于惯性坐标系的加速度

  •  
  • 本地

    相对于局部坐标系的加速度

  •  
  • A类*

    冲程振幅

  •  
  • 应收账

    翼尖-翼尖展弦比

  •  
  • 应收账

    单翼展弦比

  •  
  • b条

    爆裂涡(图4)

  •  
  • b条

    单翼翼展

  •  
  • c(c)

    平均翼弦长度

  •  
  • C类

    修正系数机翼展弦比数据

  •  
  • C类L(左)

    升力系数

  •  
  • C类D类

    阻力系数

  •  
  • C类D0(数字0)

    零升力阻力系数

  •  
  • 计算流体力学

    计算流体动力学

  •  
  • D类

    阻力

  •  
  • D类和弦

    弦向阻力

  •  
  • DPIV公司

    数字粒子图像测速

  •  
  • 埃克

    埃克曼数

  •  
  • (f)中心

    比向心力(单位体积)

  •  
  • (f)科尔

    科里奥利比力(单位体积)

  •  
  • 引力常数

  •  
  • J型

    预付款比率

  •  
  • 时间平均升力

  •  
  • L(左)

    升力

  •  
  • 列维

    前缘涡

  •  
  • 群众

  •  
  • N个

    实验次数

  •  
  • 第页

    半径矢量大小

  •  
  • 第页

    流体颗粒在旋转框架中的位置

  •  
  • R(右)

    机翼半径

  •  
  • R(右)

    机翼回转半径

  •  
  • 重新

    雷诺数

  •  
  • 重新

    回转半径处的雷诺数

  •  
  • Ro公司

    罗斯比数

  •  
  • 有效值。

    均方根

  •  
  • 以弦长表示的行驶距离

  •  
  • S公司

    单翼面积

  •  
  • 某人

    螺旋爆裂涡(图4)

  •  
  • 标准差。

    标准偏差

  •  
  • 在以下位置移动的弦长数R(右)在全冲程期间(图3)

  •  
  • U型

    回转半径处的平均速度

  •  
  • u个本地

    局部坐标系中的速度

  •  
  • V(V)

    沿机翼半径的速度

  •  
  • W公司

    重量

  •  
  • (x、 y,z)

    局部坐标系

我们感谢Will Dickson在实验设置、宝贵建议和手稿校对方面的帮助,以及Andrew Straw在视频设置方面的帮助。我们还感谢道格拉斯·阿尔赫舒勒(Douglas Althshuler)的宝贵意见,并借给了他的防水力传感器。我们感谢乌尔里克·穆勒、吉姆·乌舍伍德、米斯·穆勒、约翰·达比里和格特扬·范·海伊斯特的宝贵意见。我们感谢科特·林登堡校对手稿和数学推导。我们感谢约翰·范·吕文(Johan van Leeuwen)对各种版本手稿的衷心支持、鼓励和校对。最后,D.L.希望感谢彼得·巴克尔、海丝特·比吉尔和巴斯·范·乌德尤斯登帮助他获得本研究的旅行奖学金。这项研究得到了荷兰科学研究组织旅行奖学金的支持实验生物学杂志以及J.M.Burgerscentrum流体动力学研究和NWO-ALW向D.L.提供的817.02.012号拨款,以及国家科学基金会(IBN-0217229)和帕卡德基金会(2001-17741A)向M.H.D.提供的拨款。

Altshuler,D.、Dudley,R.和Ellington,C.P。(
2004
). 蜂鸟旋转翅膀和翅膀模型的空气动力。
J.佐尔。伦敦。
264
327
-332.
Altshuler,D.L。(
2001
).
蜂鸟沿海拔梯度飞行的生态生理学:一种综合方法
美国得克萨斯大学奥斯汀分校博士论文。
Altshuler,D.L.、Dickson,W.B.、Vance,J.T.、Roberts,S.P.和Dickinson,M.H。(
2005
). 短幅高频机翼划动决定了蜜蜂飞行的空气动力学。
程序。国家。阿卡德。科学。美国
102
18213
-18218.
Aono,H.,Liang,F.和Liu,H。(
2008
). 昆虫悬停飞行中的近场和远场空气动力学:综合计算研究。
实验生物学杂志。
211
239
-257.
A.阿祖玛。(
2006
).
飞行和游泳的生物动力学
弗吉尼亚州雷斯顿:AIAA Books。
Azuma,A.和Yasuda,K。(
1989
). 旋转种子的飞行性能。
J.西奥。生物。
138
23
-53.
Bandyopadhyay,P.R.、Beal,D.N.和Menozzi,A。(
2008
). 生物机器人对动物如何游泳的见解。
实验生物学杂志。
211
206
-214.
H·巴鲁。(
1999
).
分析动力学
纽约:McGraw-Hill。
Beom-Seok,K,Jeong-Hwan,K,Koji,K,van Rooij,R.P.和Young-Ho,L。(
2002
).
Delft大学T40/50模型水平轴风力机功率特性的三维数值预测
日本名古屋:第五届JSME-KSME流体工程会议。
Birch,J.M.,Dickinson,M.H。(
2001
). 翼展流动和前缘涡流在昆虫翅膀上的附着。
自然
412
729
-733.
Birch,J.M.,Dickinson,M.H。(
2003
). 扑翼飞行中翼面相互作用对产生空气动力的影响。
实验生物学杂志。
206
2257
-2272.
Birch,J.M.、Dickson,W.B.和Dickinson,M.H。(
2004
). 高、低雷诺数下扑翼前缘涡的产力和流动结构
实验生物学杂志。
207
1063
-1072.
Chai,P.和Millard,D。(
1997
). 飞行和尺寸限制:大型蜂鸟在最大载荷下的悬停性能。
实验生物学杂志。
200
2757
-2763.
Combes,S.A.和Daniel,T.L。(
2001
). 形状、拍打和弯曲:用于向前飞行的机翼和尾翼设计。
实验生物学杂志。
204
2073
-2085.
M·H·狄金森。(
1994
). 低雷诺数下机翼旋转对非定常气动性能的影响。
实验生物学杂志。
192
179
-206.
Dickinson,M.H.和Go¨tz,K.G。(
1993
). 低雷诺数下模型机翼的非定常气动性能。
实验生物学杂志。
174
45
-64.
Dickinson,M.H.,Lehmann,F.O.和Sane,S.P。(
1999
). 翅膀旋转和昆虫飞行的空气动力学基础。
科学类
284
1954
-1960.
Dickson,W.B.和Dickinson,M.H。(
2004
). 推进比对旋转机翼空气动力学的影响。
实验生物学杂志。
207
4269
-4281.
Dumitrescu,H.和Cardos,V。(
2003
). 风力涡轮机边界层流动的旋转效应。
美国汽车协会J。
42
408
-411.
C.P.艾灵顿。(
1984
). 昆虫飞行的空气动力学I-VI。
菲洛斯。事务处理。R.Soc.伦敦。生物科学B。
305
1
-181.
Ellington,C.P.和Usherwood,J.R。(
2001
). 旋转和扑翼的升力和阻力特性。
微型飞行器应用中的固定翼和扑翼空气动力学
(编辑:T.J.Mueller)。弗吉尼亚州雷斯顿:AIAA图书。
Ellington,C.P.、Van den Berg,C.、Willmott,A.P.和Thomas,A.L.R。(
1996
). 昆虫飞行中的前缘旋涡。
自然
384
626
-630.
Fry,S.N.、Sayaman,R.和Dickinson,M.H。(
2003
). 果蝇自由飞行机动的空气动力学。
科学类
300
495
-498.
D.I.格林维尔。(
2002
). 三角翼涡破裂的简单工程模型。
J.飞机
40
402
-405.
Greitzer,E.M.、Tan,C.S.和Graf,M.B。(
2004
).
内部流概念和应用
剑桥:剑桥大学出版社。
F.A.哈特曼。(
1963
). 蝙蝠的一些飞行力学。
俄亥俄州科学杂志。
53
59
-65之间。
Himmelskamp,H。(
1947
).
旋转螺旋桨的剖面研究
《1945年Goéttingen论文》,报告与翻译第832号。
Hove,J.R.、O’Bryan,L.M.、Gordon,M.S.、Webb,P.W.和Weihs,D。(
2001
). Boxfish(真骨纲:介形科)作为鱼类多鳍游泳的模型系统:运动学。
实验生物学杂志。
204
1459
-1471.
T·胡贝尔。(
2006
).
Untersuchungen zur instationa¨ren Aerodynamik an einem vogela¨hnlichen Flu¨gelschlagmodell公司
,博士论文,TU Darmstadt,Fachbereich Biologie。
Jones,G.、Webb,P.I.、Sedgeley,J.A.和O'Donnell,C.F.J。(
2003
). 神秘神秘蝙蝠:新西兰短尾蝙蝠(结核神秘蝙蝠)如何定位昆虫猎物。
实验生物学杂志。
206
4209
-4216中。
科尔文,J.E。(
1986
). 船用螺旋桨。
每年。流体力学版次。
18
367
-403.
F.O.莱曼和M.H.迪金森。(
1998
). 果蝇(Drosophila spp.)翅膀运动学和飞行力的控制。
实验生物学杂志。
201
385
-401.
Lentink,D.和Dickinson,M.H。(
2009
). 扑翼、旋转和平移鳍和翅膀的生物流体动力学缩放。
实验生物学杂志。
212
2691
-2704.
Lentink,D.,Muéller,UK,Stamhuis,E.J.,de Kat,R.,van Gestel,W.,Veldhuis,L.L.M.,Henningsson,P.,Hedenstro­M,A.,Videler,J.J.和van Leeuwen,J.L。(
2007
). 雨燕如何通过变形翅膀控制滑翔性能。
自然
446
1082
-1085.
Lentink,D.、Muijres,F.T.、Donker-Duyvis,F.J.和van Leeuwen,J.L。(
2008
). 模拟动物游泳和飞行的拍击箔的涡尾相互作用。
实验生物学杂志。
211
267
-273.
林登堡,C。(
2004
).
基于工程考虑和测量的旋转增强建模
伦敦:欧洲风能会议。
Lu,Y.,Shen,G.X.和Lai,G.J。(
2006
). 扑翼上的双前缘旋涡。
实验生物学杂志。
209
5005
-5016.
H.J.马登。(
1987
). 飞行动物起飞时产生的最大升力。
实验生物学杂志。
130
235
-258.
T·麦克斯沃西。(
1979
). 悬停飞行昆虫产生升力的Weis-Fogh机理实验。第1部分:。“投掷”的动态。
J.流体力学。
93
47
-63.
麦克马洪,T.A。(
1984
).
肌肉、反射和运动
新泽西州普林斯顿:普林斯顿大学出版社。
Miller,L.A.和Peskin,C.S。(
2004
). 当漩涡粘住时:微小昆虫飞行的空气动力学转变。
实验生物学杂志。
207
3073
-3088.
Muijres,F.T.、Johansson,L.C.、Barfield,R.、Wolf,M.、Spedding,G.R.和Hedenstro­M,A。(
2008
). 前缘涡流提高了慢速飞行蝙蝠的升力。
科学类
319
1250
-1253.
Norberg,U.M.、Brooke,A.P.和Trewhella,W.J。(
2000
). 翼手类飞翔和非飞翔蝙蝠(飞狐,翼手类):翅膀形态和飞行性能。
实验生物学杂志。
203
651
-664.
Poelma,C.、Dickson,W.B.和Dickinson,M.H。(
2006
). 动态缩放扑翼周围全速度场的时间分辨率重建。
实验流体
41
213
-225.
林盖特(M.Ringuette)、米兰(M.Milano)和加里布(M.Gharib)。(
2007
). 叶尖涡在低展弦比法向平板产生力中的作用。
J.流体力学。
581
453
-468.
C.G.罗斯比。(
1936
). 关于海面上的动量传递。第一部分论文
物理学。Oceanogr公司。美托洛尔。
3号。
Ruijgrok,G.J.J。(
1994
).
飞机性能要素
荷兰代尔夫特:代尔夫脱大学出版社。
Sane,S.P.和Dickinson,M.H。(
2001
). 扑翼对飞行力的控制:升力和阻力的产生。
实验生物学杂志。
204
2607
-2626.
H·施利钦。(
1979
).
边界层理论,第7版。
纽约:McGraw-Hill。
F.W.施密茨。(
1942
).
C.J.E.Volckmann Nachf,Aerodynamik des flugmodells公司。E类
柏林威特。
Simonich,J.、McCormick,D.和Haas,R。(
1992
). 起飞时螺旋桨前缘旋涡的流动可视化。
艾亚
0386
.
Srygley,R.B.和Thomas,A.L.R。(
2002
). 自由飞行蝴蝶的非传统升力产生机制。
自然
420
660
-664.
Stepniewski,W.Z.和Keys,C.N。(
1984
).
旋转摆动空气动力学
纽约州米诺拉:多佛出版公司。
Stiles,F.G.、Altshuler,D.L.和Dudley,R。(
2005
). 一些北美蜂鸟的翅膀形态和飞行行为。
奥克
122
872
-886.
Sun,M.和WU,J。(
2004
). 低雷诺数下扫翼上的大空气动力。
力学学报
20
24
-31.
汤格勒,J.L。(
2004
). 深入了解风力涡轮机失速和失速后空气动力学。
风能
7
247
-260.
Tennekes H。(
1997
).
飞行的简单科学
马萨诸塞州剑桥:麻省理工学院出版社。
Usherwood,J.R.和Ellington,C.P。(
2002
). 旋转机翼的空气动力学I-II。
实验生物学杂志。
205
1547
-1576.
Van Den Berg,C.和Ellington,C.P。(
1997
). “悬停”模型鹰蛾的三维前缘旋涡。
菲洛斯。事务处理。R.Soc.伦敦。生物科学B。
352
329
-340.
J.P.Vanyo。(
1993
).
工程与科学中的旋转流体
纽约州米诺拉:多佛出版公司。
Vaughan,N.、Parsons,S.、Barlow,K.E.和Gannon,M.R。(
2004
). 西印度群岛蝙蝠的回声定位叫声和翅膀形态。
《Chiropter学报》。
6
75
-90.
维德勒,J.J.、斯塔姆休斯,E.J.和波维尔,G.D.E。(
2004
). 前沿旋涡提升雨燕。
科学类
306
1960
-1962.
沃克,J.A。(
2004
). 极端划桨鱼类胸鳍划动的动力学:三棘棘棘棘鱼(Gasterosteus aculetus)
实验生物学杂志。
207
1925
-1939.
J.A.Walker和M.W.Westnight。(
2002
). 实验室推进的性能极限与鳍形状和运动相关。
实验生物学杂志。
205
177
-187.
王振杰。(
2008
). 扑翼飞行的气动效率:一个二冲程模型的分析
实验生物学杂志。
211
234
-238.
Wang,Z.J.,Birch,J.M.和Dickinson,M.H。(
2004
). 低雷诺数悬停飞行中的非定常力和流动:二维计算与机器人机翼实验
实验生物学杂志。
207
449
-460.
D.R.沃里克、B.W.托巴尔斯克和D.R.鲍尔斯。(
2005
). 悬停蜂鸟的空气动力学。
自然
435
1094
-1097年。
H·温特。(
1936
).
短跨距平板和翼型上的流动现象
NACA报告798。
山田,T.和铃木,E。(
1999
). 亚洲梧桐科带翅一水硬铝石的比较形态和异速生长。
J.特罗普。经济。
15
619
-635.