PIV experimental study on vortex structures induced by free autorotation fall of a samaras
-
摘要: 本文利用粒子图像测速技术研究种子叶片自由旋转下落过程中不同涡系的相互作用机理。以典型枫树种子叶片为研究对象,通过对比不同叶片长度、叶面厚度、叶面宽度、下落锥角、自旋角速度、下落速度和下落初始角度等参数对过渡期和稳定期的影响,分析了诱发叶片自旋的外形特征和空间特性。结合特征参数分析,对叶片自旋稳定期开展了PIV流场测量实验,解析了无干扰下种子叶片下落过程中涡系的产生和演化机理。实验结果表明:稳定期气流在叶尖正面位置产生前缘涡(沿展向呈圆锥状结构),后缘位置产生反方向的后缘涡;两个涡发生相互耦合运动,前缘涡的强度大于后缘涡,从而导致叶面产生锥角。在前缘和叶尖前方观测到较高的速度向上的区域,而在后缘和叶根附近则出现较高的速度向下的区域,从而对种子产生向上的升力,使叶片实现自旋稳定下落。通过枫叶种子自由下落的无干扰PIV测量,初步获得了贴近叶片表面前缘涡的运动性状,验证了后缘涡的存在,结论对单翼型旋转叶片的设计有一定指导意义。Abstract: Taking the typical maple samara blade as the research object, the flow field of the free rotation and falling process is measured by the particle image velocimetry, to study the evolution and the law of the spin flow structure, which has certain guiding significance for the design of the single-wing aircraft. By comparing the effects of different blade lengths, thicknesses, widths, falling angles, spin angular velocities, falling velocities, and different falling attitudes on the transition period and stability period, the shape and spatial characteristics of the induced blade spin are obtained. Combined with the results of characteristic parameter analysis, PIV flow field measurement experiments are carried out for the period of blade spin stability, and the rules of vortex generation and evolution in the process of falling are obtained. The experimental results show that the leading-edge vortices (conical structure along the spanwise direction) are generated at the front tip during the stable period, the trailing edge vortices in the reverse direction are generated at the back-edge, and the two vortices are coupled with each other. The strength of the leading-edge vortex is greater than that of the trailing-edge vortex, which leads to the angle of attack. In front of the leading edge and the tip of the leaf, a higher upward velocity region was observed, and a higher downward velocity region was observed near the rear edge and the root of the leaf, which resulted in an upward lifting force on the seed and a stable falling of the blade.ade.
-
Keywords:
- maple samara /
- free autorotation fall /
- vortex /
- LEV /
- PIV
-
0 引 言
经过长期自然选择,鲨鱼体表形成了大面积错位排布的微沟槽型鳞片[1]。以往研究证实,该微结构降低了游动阻力,提高了鲨鱼游速。基于仿生学原理,采用现有CFD仿真技术[2]和微纳制造技术[3],可以制造出仿生鲨鱼皮结构,将其应用于飞行器减阻,能够有效提高飞行速度、降低能耗、延长续航时间。
目前,对仿生微纳减阻结构的研究主要集中于沟槽结构。Walsh[4]研究发现,三角形截面的沟槽结构减阻效果极佳,减阻率最高可达8%。SPEEDO公司生产的鲨鱼皮泳衣通过仿生鲨鱼皮微结构控制人体周围流场,使阻力降低7.5%左右[5]。刘志华等[6]对V型沟槽的研究表明,微沟槽尖峰具有扰流作用,其减阻率可达6.6%。田丽梅等[7]通过对仿生旋成体表面沟槽结构进行仿真计算发现,在44 m/s的速度下,棱纹结构可使旋成体总阻力降低11.12%。上述研究都证实了仿生微沟槽结构减阻的可行性。
本文在大量微纳减阻结构研究成果的基础上,采用CFD仿真与激光微纳制造相结合的方式,设计了仿生鲨鱼皮复合微纳减风阻结构,并在大气传感器半球头体模型表面完成了复合微纳结构的激光制备[8]。基于实验验证与减阻机理分析,对复合微纳结构进行了优化,使减风阻效果得到提升。
1 模型建立与理论分析
1.1 建立模型
1.1.1 建立鲨鱼皮鳞片结构模型
大量研究表明,鲨鱼皮鳞片结构具有很好的减阻效果,相关研究成果也证实了利用鲨鱼皮鳞片结构减阻的可行性。为了客观反映减阻效果及流体流动特性,根据鲨鱼皮鳞片的实际尺寸抽象出了单个鳞片的几何模型。单个鳞片及流场域参数如表1所示,其中,lc、ls分别为鳞片结构的横向长度和纵向长度,h1、h2分别为鳞片中间凸起的高度和两侧凸起的高度,Lx、Ly和Lz分别为计算域的流向长度、法向长度和展向长度。
表 1 实验模型参数Table 1 Experimental model parameterslc = ls h1 = h2 Lx Ly Lz 400 µm 150 µm 350 mm 350 mm 260 mm 采用作图软件UG构建了鲨鱼皮鳞片结构模型(图1),确定了鳞片排布组合方式,并在飞行器的大气传感器半球头体模型表面完成拼接覆盖(图2)。为尽可能在半球头体表面均匀铺满鳞片且不破坏鳞片间扰流的协同作用,基于计算结果将半球头体表面等分为6个区域,按照6个朝向呈放射状等距排列鳞片,每2个鳞片间距均为100 µm(图3)。
1.1.2 确定二级微沟槽尺寸
仅依靠鲨鱼皮鳞片结构进行减阻,减阻效果有限。在鲨鱼皮鳞片结构减阻基础上,采用二级微沟槽构造增强减阻效果,即在鲨鱼皮鳞片结构上覆盖尺度更加微小的沟槽结构。基于丛茜等[9]对多种仿生非光滑沟槽减阻效果的验证,同时考虑来流条件下结构的稳定性,选取V型沟槽作为二级微沟槽结构。引入微沟槽无量纲深度h+和无量纲宽度s+进行计算,获得具有最佳减阻效果的临界条件[10]:
$$ {s^ + } = \frac{{s{v_\infty }}}{\nu }\sqrt {\frac{{{c_f}}}{2}} $$ (1) $$ {h^ + } = \frac{{h{v_\infty }}}{\nu }\sqrt {\frac{{{c_f}}}{2}} $$ (2) 式中:s为槽宽,h为槽深;v∞为风速;ν为气流的运动黏度;摩擦系数cf = 0.072(Re)$ ^{ - \frac{1}{5}} $。基于Walsh等的研究以及制造的可行性,选取具有最佳减阻效果的临界条件s+ = h+,即s = h = 10 µm。二级微沟槽结构如图4(a)所示。
在气固两相绕流体动态模拟过程中[11-12],为减缓速度梯度过大导致的湍流,将微沟槽截面设计为具有周期结构的类正弦状图形[13-14]。二级微沟槽则是在鳞片结构上进一步加工的微结构。为有效叠加两级微沟槽的减阻效果,微元区域需选取为包含鳞片和两鳞片间距在内的正方形区域(边长500 µm),且该区域所有二级微沟槽沿鳞片结构纵向开槽并贯穿整个鳞片表面,其截面如图4(b)所示。采用微元法仿真小区域微结构的气流流动过程,得出该区域的减阻率,继而通过积分表征出带有二级微沟槽结构的整个半球头体的减阻水平。
1.1.3 选择求解模型
在风洞实验过程中,气流在半球头体表面产生摩擦阻力和压差阻力,并具有一定黏附性。由于湍流的存在,半球头体受力不断变化。本文选取k–ω SST模型进行求解[15],该模型对解的初始预测非常敏感,对边界层附近区域的网格划分也较为精细[16],计算精度高,对壁面气流流动黏性和分离界限的计算尤其精细,在工业应用中较为普适。
1.2 理论分析
在一定速度来流条件下,飞行器表面会出现湍流,需采用N–S方程进行求解与分析。式(3)即为该过程的N–S方程,式(4)~(6)为分解到x、y、z方向的动量方程[17]:
$$ {\text{div }}{ {σ}} {\text{ = 0}} $$ (3) $$ \frac{{\partial {\boldsymbol{u}}}}{{\partial t}} + {\text{div}}( {{\boldsymbol{u}}{ {σ}} } ) = - \frac{1}{\rho }\frac{{\partial {\boldsymbol{p}}}}{{\partial x}} + {\boldsymbol{v}}{\text{div}}\left[ {{\text{grad}}( {\boldsymbol{u}} )} \right] $$ (4) $$ \frac{{\partial {\boldsymbol{v}}}}{{\partial t}} + {\text{div}}( {{\boldsymbol{v}}{ {σ}} } ) = - \frac{1}{\rho }\frac{{\partial {\boldsymbol{p}}}}{{\partial y}} + {\boldsymbol{v}}{\text{div}}\left[ {{\text{grad}}( {\boldsymbol{v}} )} \right] $$ (5) $$ \frac{{\partial {\boldsymbol{w}}}}{{\partial t}} + {\text{div}}( {{{{\boldsymbol{w}}{ {σ}}}} } ) = - \frac{1}{\rho }\frac{{\partial {\boldsymbol{p}}}}{{\partial z}} + {\boldsymbol{v}}{\text{div}}\left[ {{\text{grad}}( {\boldsymbol{w}} )} \right] $$ (6) 将x方向的动量方程进行速度与压力分解,并代入N–S方程,整理后可得:
$$ \frac{{\partial {\boldsymbol{U}}}}{{\partial t}} + {\rm{div}}( {{\boldsymbol{U}}{ {τ}} } ) + {\rm{div}}( {\overline {{{\boldsymbol{u}}'}{{ {σ}} '}} } ) = - \frac{1}{\rho }\frac{{\partial {\boldsymbol{P}}}}{{\partial x}} + {\boldsymbol{v}}{\rm{div}}\left[ {{\rm{grad}}( {\boldsymbol{U}})} \right] $$ (7) 式中:${ {σ}} $为气流的速度矢量;${\boldsymbol{u}}、{\boldsymbol{v}}、{\boldsymbol{w}} $分别为${ {σ}}$分解到x、y、z方向的速度矢量;${ {τ}}$和${\boldsymbol{U}} $分别为气流时均速度和x方向的时均速度;${\boldsymbol{p}} $和${\boldsymbol{P}} $分别为气流对目标体的压力和时均压力;${ {σ}}'$和$ {\boldsymbol{u}}' $分别为气流脉动速度和x方向气流脉动速度;ρ为流体密度。
可以看到,式(7)增加了一个附加项$ \operatorname{div}(\overline{{\boldsymbol{u}} '{ {σ}}'}) $,即雷诺应力项。同理,其他两个方向也都增加了一个雷诺应力项。雷诺应力项的出现引起了湍流现象,说明该项即为湍流产生的条件[18]。
半球头体表面的速度梯度很大,因而涡量很大,这也是流场中旋涡形成的主要原因,边界层分离会导致旋涡的周期性脱落。图5(a)和(b)为仿生鲨鱼皮横、纵截面的气流流动视图,沟槽内回转的涡流受到鲨鱼皮鳞片凸起的阻碍,产生“二次涡”,减弱了微结构表面与低速条带相关联的涡流强度,抑制了低速条带的形成,并削弱了低速气流运动的不稳定性,即低速条带缓慢向上提升的过程被减弱,流体微团间的动量交换降低,从而达到减阻效果[19-20]。
2 实验与讨论
为更好地进行对比分析,将光滑半球头体与覆盖仿生鲨鱼皮鳞片结构的半球头体置于同一流场域中,在相同条件下进行仿真,风洞检测模型如图6所示。严格设定边界条件,将左右壁面定义为对称边界,上下表面无滑移,流体介质设置为空气,出入口设置为速度入口、压力出口。
网格划分是有限元处理的关键[21],网格质量决定了仿真结构正确与否。在CFD仿真过程中,覆盖于半球头体表面的每个仿生鲨鱼皮鳞片都需要极高精度的网格划分进行表征。由于计算机CPU、GPU运行速度有限,为提高计算速度及效率,将半球头体对称均分为24份,取其中之一代替整体进行风洞实验的仿真。图7为仿生鲨鱼皮的网格划分,总网格数为34300011,实体离散化节点7262815个,网格最小尺寸为10 µm,平均网格质量达0.88,网格质量良好,有助于计算收敛。
网格划分完成后进行仿真计算,根据空气阻力系数公式(8)进行求解:
$$ C = \frac{{2F}}{{\rho v_\infty ^2A}} $$ (8) 式中,C为空气阻力系数,F为阻力,ρ为空气密度,A为迎风面积,v∞为风速。
针对覆盖仿生鲨鱼皮鳞片的半球头体和对照组光滑半球头体,分别求解150 km/h(41.67 m/s)和160 km/h(44.44 m/s)风速下的阻力及阻力系数。当迭代曲线收敛后,根据积分原理还原出整个模型的受力情况,最后采用式(9)求得减阻率:
$$ {\text{Δ}} C = \frac{{\left| {{C_1} - {C_0}} \right|}}{{{C_0}}} $$ (9) 由仿真计算(计算方法详见文献[22])结果得到阻力值及阻力系数,如表2所示,2个风速下的减阻率分别为16.9%和18.4%。
表 2 CFD仿真参数Table 2 CFD simulation parameter样品 v∞/(km·h–1) F/N C ΔC 光滑半球头体 150 1.499 0.413 - 覆盖鳞片的半球头体 150 1.245 0.343 16.9% 光滑半球头体 160 1.731 0.419 - 覆盖鳞片的半球头体 160 1.407 0.342 18.4% 由式(8)可知,当空气密度、迎风面积和风速均相同时,阻力系数与阻力线性相关,两者趋势相同。迭代曲线稳定并收敛后,同一风速下覆盖仿生鲨鱼皮鳞片的半球头体的迭代次数比光滑半球头体明显减少。微结构的存在使迭代曲线更快收敛。为更直观地表征减阻效果,通过后处理得到光滑半球头体与覆盖仿生鲨鱼皮鳞片的半球头体同一位置的压力云图和速度矢量云图,如图8~10所示。
由图8和9可知,在半球头体前端,具有一定速度的气流在垂直方向(y方向)受到半球头体阻滞,瞬时速度降为0,全部动压转变为静压,在迎角微小的前端,存在一个较大正压区。图8中光滑球头体表面的边界层很薄,与之相比,由于鳞片的存在,图9中边界层厚度明显增大,混合层厚度减小,结合等高线对应数值跨度可以看到来流速度明显降低。
对比图10(a)与(b)可以看到:1/24光滑半球头体表面压力云图颜色明显分层,表明压力损失较大;而覆盖仿生鲨鱼鳞片的半球头体表面压降变化不明显,无明显压力损失,受力更加均匀,受力最大的正压区阻力也明显降低,可见仿生鲨鱼皮鳞片实现了很好的减阻效果。
在相同边界条件下对带有二级微沟槽的微元区域仿真计算,得到稳定的阻力系数为1.896 × 10−9,而光滑微元区域的阻力系数为1.986 × 10−9,减阻率为4.5%。经计算,两模型在风速为150和160 km/h时的雷诺数分别为2.7 × 106和2.9 × 106,且壁面第一层网格单元的y+均集中于0.8~2.0之间,进一步证明了边界层网格划分的合理性。
经后处理,得到微元区域二级微沟槽气流速度矢量云图,如图11所示。从x = 0、z = 0截面的云图可以看出:二级微沟槽近壁区边界层厚度明显大于光滑表面对应区域,且有增厚趋势;混合层厚度变薄,湍流被削弱,壁面速度梯度降低,同时压力也明显降低,从而导致阻力降低,实现了减阻效果。
基于仿真结果,按照设定尺寸对半球头体表面进行激光微纳加工,制造了一批成型钛合金(TC4)半球头样件并作预处理[23]。实验中使用了3个半球头体样件:表面光滑的半球头体样件A作为对照组;以相同工序在样件B和C表面激光扫描仿生鲨鱼皮微结构;进一步加工样件C,通过激光干涉加工二级微沟槽,形成复合结构。图12和13为制备完成后各样件表面微结构的SEM电镜图。
最后,通过风洞实验验证减阻效果。以压力传感器实时记录阻力值,并实时计算出减阻率,结果如表3所示。当风速为150 km/h时,与对照组样件A相比,样件B的阻力降低了8.7%,样件C的阻力降低了10.3%。半球头体表面的仿生鲨鱼皮结构及复合减阻结构均表现出了优秀的减阻性能。
表 3 风洞实验参数表Table 3 Parameters of wind tunnel experiment样品 v∞/(km·h–1) F/N C ΔC 光滑表面样件A 150 1.386 0.3823 - 非光滑表面样件B 150 1.266 0.3491 8.7% 非光滑表面样件C 150 1.243 0.3430 10.3% 光滑表面样件A 160 1.565 0.3804 - 非光滑表面样件B 160 1.442 0.3504 7.9% 非光滑表面样件C 160 1.407 0.3419 10.1% 对比表2和3可以发现:在残差精度范围内,覆盖微结构的样件的阻力系数风洞实验与仿真结果一致,验证了仿真方法的可行性。由此可见,减阻率出现偏差,原因在于光滑半球头体。对比表2和3的数据发现,光滑半球头体受力的仿真结果与实际检测结果存在偏差。从光滑半球头体风洞实验结果可知,微观角度的真实物体表面并非绝对光滑,其无规则的微观结构会导致阻力的变化。理想化模型无需考虑机械能转化和内能消耗,简化了流动过程,并能客观精确地反映流动特性。因此,可以认为仿真计算与风洞实验得到的规律一致。
3 结 论
1)仿生鲨鱼皮与二级微沟槽的复合结构具有很好的减阻效果。在仿生鲨鱼皮鳞片结构基础上均匀覆盖更小尺度的二级微沟槽结构,能够延长“二次涡”的回旋路径,降低沟槽底部气流速度梯度,减小气流对壁面的黏附力,削弱来流强度,增大边界层厚度,进一步提升减阻效果。
2)仿真结果与风洞实验结果一致。复合微纳减阻结构的减阻率可达10.3%,验证了仿真方法的可行性。仿真方法可以替代一部分风洞实验预测及评估减阻效果,从而降低成本、提高效率。
本文采用仿真与风洞实验手段进行研究,以验证结果的可靠性。受风洞实验条件所限,仅对当前所能达到的最高风速(即非常临近的2个风速条件150和160 km/h)进行了仿真与实验。在后续工作中,亟待在更高风速和不同来流方向条件下开展复合微纳减阻结构的减阻效果研究。
-
表 1 枫树种子平均外形特征参数统计
Table 1 Average shape characteristic parameters of maple samaras
类型 叶弦长
L/mm叶面最大宽度
b/mm叶面平均厚度
k/mm果实厚度
t/mm种子质量
M/mg小型 17.3 5.7 0.2 3.8 47.0 中型 20.0 6.7 0.2 4.0 49.3 大型 22.7 8.7 0.2 4.0 52.0 表 2 典型中型枫树种子表面平均特征参数
Table 2 Average parameters of typical surface characteristics
叶弦长
L/mm平均叶脉
数量叶面平均
厚度k/mm叶凹凸幅度/
mm叶凹凸高度/
mm前缘厚度/
mm20.0 70 0.2 0.15 0.18 0.53 表 3 枫树种子在稳定期中自旋运动的平均特性参数
Table 3 Average parameters of maple samaras in falling
类型 雷诺数
Re下落锥角
β/(°)下落速度
V/(m·s–1)自旋角速度
ω/(rad·s–1)过渡期距离/
cm小型 19 300 16.20 1.12 45.00 32.00 中型 22 510 22.61 1.13 34.14 33.75 大型 21 480 26.45 0.95 19.00 33.00 -
[1] NORBERG R Å. Autorotation, self-stability, and structure of single-winged fruits and seeds (samaras) with comparative remarks on animal flight[J]. Biological Reviews,1973,48(4):561-596. doi: 10.1111/j.1469-185x.1973.tb01569.x
[2] AZUMA A,YASUDA K. Flight performance of rotary seeds[J]. Journal of Theoretical Biology,1989,138(1):23-53. doi: 10.1016/S0022-5193(89)80176-6
[3] GREENE D F,JOHNSON E A. Seed mass and dispersal capacity in wind-dispersed diaspores[J]. Oikos,1993,67(1):69. doi: 10.2307/3545096
[4] BULLOCK J M, KENWARD R, et al. Dispersal Ecology[C]//Proc of 42nd Symposium of the British Ecological Society. 2008.
[5] MINAMI S,AZUMA A. Various flying modes of wind-dispersal seeds[J]. Journal of Theoretical Biology,2003,225(1):1-14. doi: 10.1016/S0022-5193(03)00216-9
[6] LEE I,CHOI H. Scaling law for the lift force of autorotating falling seeds at terminal velocity[J]. Journal of Fluid Mechanics,2018,835:406-420. doi: 10.1017/jfm.2017.746
[7] LEISHMAN J G. Principles of helicopter aerodynamics with CD extra[M]. UK: Cambridge University press, 2006.
[8] ULRICH E R,PINES D J,HUMBERT J S. From falling to flying: the path to powered flight of a robotic samara nano air vehicle[J]. Bioinspiration & Biomimetics,2010,5(4):045009. doi: 10.1088/1748-3182/5/4/045009
[9] LUGT H J. Autorotation[J]. Annual Review of Fluid Mechanics,1983,15(1):123-147. doi: 10.1146/annurev.fl.15.010183.001011
[10] SKEWS B W. Autorotation of many-sided bodies in an airstream[J]. Nature,1991,352(6335):512-513. doi: 10.1038/352512a0
[11] USHERWOOD J R,ELLINGTON C P. The aerodynamics of revolving wings II. Propeller force coefficients from mayfly to quail[J]. Journal of Experimental Biology,2002,205(11):1565-1576. doi: 10.1242/jeb.205.11.1565
[12] VARSHNEY K,CHANG S,WANG Z J. The kinematics of falling maple seeds and the initial transition to a helical motion[J]. Nonlinearity,2012,25(1):C1-C8. doi: 10.1088/0951-7715/25/1/c1
[13] SMITH E H. Autorotating wings: an experimental investigation[J]. Journal of Fluid Mechanics,1971,50(3):513-534. doi: 10.1017/s0022112071002738
[14] LAU E M,HUANG W X,XU C X. Progression of heavy plates from stable falling to tumbling flight[J]. Journal of Fluid Mechanics,2018,850:1009-1031. doi: 10.1017/jfm.2018.486
[15] LEE E J,LEE S J. Effect of initial attitude on autorotation flight of maple samaras (Acer palmatum)[J]. Journal of Mechanical Science and Technology,2016,30(2):741-747. doi: 10.1007/s12206-016-0129-2
[16] MYONG H S. Effects of the configuration characteristics on the motion parameters of autorotating flight of plant seeds[C]//Proc of the 5th International Conference on Experimental Fluid Mechanics. 2018.
[17] LENTINK D,DICKSON W B,VAN LEEUWEN J L,et al. Leading-edge vortices elevate lift of autorotating plant seeds[J]. Science,2009,324(5933):1438-1440. doi: 10.1126/science.1174196
[18] ENGELS T,KOLOMENSKIY D,SCHNEIDER K,et al. Bumblebee flight in heavy turbulence[J]. Physical Review Letters,2016,116(2):028103. doi: 10.1103/physrevlett.116.028103
[19] BIRCH J M,DICKSON W B,DICKINSON M H. Force production and flow structure of the leading edge vortex on flapping wings at high and low Reynolds numbers[J]. Journal of Experimental Biology,2004,207(7):1063-1072. doi: 10.1242/jeb.00848
[20] SHYY W,LIU H. Flapping wings and aerodynamic lift: the role of leading-edge vortices[J]. AIAA Journal,2007,45(12):2817-2819. doi: 10.2514/1.33205
[21] SALCEDO E,TREVIÑO C,VARGAS R O,et al. Stereoscopic particle image velocimetry measurements of the three-dimensional flow field of a descending autorotating mahogany seed (Swietenia macrophylla)[J]. The Journal of Experimental Biology,2013,216(11):2017-2030. doi: 10.1242/jeb.085407
[22] LEE S J,LEE E J,SOHN M H. Mechanism of autorotation flight of maple samaras (Acer palmatum)[J]. Experiments in Fluids,2014,55(4):1-9. doi: 10.1007/s00348-014-1718-4
[23] RAO M,HOYSALL D C,GOPALAN J. Mahogany seed - a step forward in deciphering autorotation[J]. Current Science(00113891),2014,106(8):1101-1109.
-
期刊类型引用(6)
1. 郑高杰,何小明,李东坡,谭慧俊,汪昆,吴祯龙,王德鹏. 尾部推进无人机双90°偏折进气道/蜗壳耦合流动特性. 航空学报. 2024(04): 178-191 . 百度学术
2. 肖传民,田立敏,曹远,丁冬. 航空发动机进气畸变发生器畸变特性研究. 航空精密制造技术. 2024(04): 25-27+55 . 百度学术
3. 王忠杰,黄丛磊,陈榴,戴韧. 涡流发生器抑制喷水推进器进水管总压畸变的实验研究. 水动力学研究与进展A辑. 2023(02): 270-277 . 百度学术
4. 戴韧,王忠杰,陈榴,樊丰凯,王宗龙. 基于VG/VGJ的喷水推进进口流道流动控制方法. 船舶. 2023(06): 43-49 . 百度学术
5. 郭君德,马宏伟,范聪聪,贠迪. 一种拓宽五孔压力探针测量范围的方法. 航空动力学报. 2022(11): 2659-2667 . 百度学术
6. 曹广州,李博,梁世波,谭红明. 机动飞行时长空靶机进气道气动性能数值研究(英文). Transactions of Nanjing University of Aeronautics and Astronautics. 2020(05): 739-749 . 百度学术
其他类型引用(2)