矩形前缘剪切流对圆柱气动和绕流特性影响规律的研究.pdf

第 20 卷 第 10 期2023 年 10 月铁道科学与工程学报Journal of Railway Science and EngineeringVolume 20 Number 10October 2023矩形前缘剪切流对圆柱气动和绕流特性影响规律的研究李文舒1，王汉封1,2，李欢1,2，何旭辉1,2(1.中南大学 土木工程学院，湖南 长沙 410075；2.中南大学 高速铁路建造技术国家工程研究中心，湖南 长沙 410075)摘要：横风中高速列车桥梁系统存在复杂气动干扰效应，主要表现为桥面分离流与列车绕流间的相互作用。该效应可能显著改变列车风荷载，并最终影响桥上行车稳定性。为揭示列车桥梁系统气动干扰效应机理，采用错列布置的二维圆柱和矩形钝体结构为抽象模型，通过风洞实验系统研究两者相对位置对圆柱气动力的影响规律。圆柱直径D=60 mm，矩形迎风面高度H=48 mm，基于来流风速U和D的雷诺数Re=61041105。实验发现，圆柱距矩形前沿水平距离L与两者垂直方向的间隙G对圆柱气动力均有显著影响。当G/D较小时，圆柱上侧产生较大负压区；随G/D逐渐增大，圆柱下侧产生较大负压区；当G/D进一步增大到某一临界值时，圆柱风压基本呈对称分布，圆柱尾流中开始产生周期性旋涡脱落，且涡脱频率低于孤立圆柱。随着L/D的增大，该G/D临界值逐渐增大，圆柱涡脱频率逐渐下降，涡脱强度逐渐减弱。当圆柱和矩形处在特定的相对位置时，圆柱气动力出现明显的雷诺数效应，圆柱的瞬时风压会产生2种模态。当G/D=0.4时，圆柱绕流会产生双稳态现象，且随着L/D的增大双稳态现象逐渐减弱。该现象是由于矩形前缘分离流的加速冲击在圆柱上，促使局部绕流提前转捩，从而使圆柱在亚临界区产生临界区的气动特性。研究结果可为实际工程中桥上行车的安全性和稳定性分析提供参考。关键词：圆柱；矩形；异构钝体；雷诺数效应；双稳态中图分类号：TU311 文献标志码：A 开放科学(资源服务)标识码(OSID)文章编号：1672-7029（2023）10-3638-10Effects of the leading-edge shear layer of a rectangular on the aerodynamics and flow of a circular cylinderLI Wenshu1,WANG Hanfeng1,2,LI Huan1,2,HE Xuhui1,2(1.School of Civil Engineering,Central South University,Changsha 410075,China;2.National Engineering Research Center for High-speed Railway Construction,Central South University,Changsha 410075,China)Abstract:The high-speed train and bridge system in cross wind involves complex aerodynamic interference,which is characterized by the interaction between the separated flow from the bridge and the flow around the train.This phenomenon has significant effects on the train aerodynamic load,which may threaten the running stability of the train.To reveal the mechanism of the aerodynamic interference of the train-bridge system,wind 收稿日期：2022-11-04基金项目：国家自然科学基金资助项目(52078505，52208514，U1934209，51925808)；湖南省自然科学基金资助项目(2021JJ40744)通信作者：李欢(1988)，男，河南信阳人，副教授，博士，从事桥梁风工程及风致振动控制研究；Email：DOI:10.19713/ki.43-1423/u.T20222113第 10 期李文舒，等：矩形前缘剪切流对圆柱气动和绕流特性影响规律的研究tunnel experiments were conducted with a mismatched arrangement of a 2D circular cylinder and a blunt rectangular structure,which can be considered as the simplified model of train and bridge.The effects of relative position on the wind load of the circular cylinder were systematically revealed.The cylinder diameter was D=60 mm,and the windward height of the rectangle was H=48 mm.The tested Reynolds number based on D and oncoming flow velocity U ranged from 6104 to 1105.The experimental results indicate both the horizontal space L between the cylinder and the rectangle leading edge and the gap G between them have remarkable effects on the wind load of the cylinder.When G/D is small,strong negative pressure appears on the upper side of the cylinder,while the strong negative pressure shifts to the lower side of the cylinder with increasing G/D.The pressure around the cylinder becomes symmetric when G/D is larger than a critical value,accompanied by periodic vortex shedding in the cylinder wake,with a shedding frequency lower than that of an isolated cylinder.The critical G/D enlarges gradually with increasing L/D,along with the reduction in vortex shedding frequency and strength.For particular relative positions,the instantaneous pressure of the cylinder presents two typical modes,which is similar to Reynolds number effects.For G/D=0.4,flow around the cylinder is bistable,which becomes less obvious with increasing L/D.This observation is ascribed to the accelerated shear flow from the leading edge of rectangle,which impacts on the cylinder and promotes local transition of the flow around the cylinder,even when the Reynolds number still belongs to subcritical regime.The present investigation provides references for analyzing the safety and stability of high-speed train running on a bridge in real case.Key words:circular cylinder;rectangle;isomeric bluff body system;Reynolds number effect;bistable phenomenon 桥梁电缆、电力线路、冷却塔、高层建筑及海底管道等多种圆柱形结构在土木工程、机械工程和海洋工程中得到了广泛应用。圆柱形结构周围的流场包含剪切层、尾涡和卡门涡街12之间的相互作用，这些相互作用在圆柱体周围产生了各种流动模式。在过去的几十年里，对双圆柱绕流36的气动特性和流动特性已经有了相当深入的研究，与串列和并列相比，交错排列的双圆柱在现实生活中的应用更为广泛。在不同的中心间距比 P/D(P为圆柱圆心之间的距离，D为圆柱直径)和入射角(来流与圆心连线的夹角)条件下，上游圆柱的剪切层打在下游圆柱的不同位置上，对下游圆柱的气动力和流场形态79产生明显的影响。PRICE 等10研究了在雷诺数 Re=3.21047.4104，P/D=1.1254.0，=090的条件下，2个直径相等的交错圆柱体的平均空气动力系数和涡脱频率。文章根据中心间距比对双圆柱流型进行分类。当P/D=1.1251.25时，属于紧密间距，其双圆柱的流型类似于单个钝体；当P/D=1.52.5时，属于中等间距，当入射角较小时会产生负升力以及相应的负阻力和较高的 St数(涡脱频率)；当 P/D=34时，属于大间距，2个圆柱的流场互不干扰，气动特性类似于孤立圆柱。WU等11采用大涡模拟的方法，研究了高亚临界雷诺数(Re=1.4105)条件下，2个交错圆柱在P/D=1.54，=090条件下的绕流问题。研究发现，当 P/D=1.52，=010时，高速间隙流的出现导致下游圆柱迎风侧局部区域产生强大的负压，从而产生负阻力；当P/D=1.53，=10时，下游圆柱会产生负升力；当P/D=34，=20时，下游圆柱产生正升力。GU等12研究了在高亚临界雷诺数条件(Re=2.2105)下，2个相同圆柱交错排列时的压力分布干扰。根据不同的压力分布模式，按照入射角将下游圆柱的流型分为尾迹型、剪切层型和领域干涉型。在=920时，会产生剪切层流型，此时下游圆柱的间隙面出现较大负压区，这是由于上游圆柱间隙侧分离的剪切层高速侧重新附着在下游圆柱间隙侧，以及下游圆柱间隙侧高频率有规则的旋涡脱3639铁 道 科 学 与 工 程 学 报2023 年 10月落。由于圆柱气动特性和流动特性易受雷诺数的影响，已有研究选取的雷诺数大小不一，故很难得出一致的结论。相比圆柱，矩形板的剪切层较为稳定，基本不受雷诺数影响，因此如果将上游圆柱换成无限长的矩形板，可以使下游圆柱受到稳定的流动干扰，有利于对其进行气动特性分析。此外，高速列车、桥梁缆索和风屏障檩条等结构常常浸没在桥梁前缘剪切流中，矩形板和圆柱组成的异构抽象模型类似于实际工程中的车桥模型，探究矩形前缘剪切流对下游圆柱的影响对研究桥上行车安全性和稳定性有一定的参考价值。鉴于此，本文通过测压试验，研究圆柱的气动特性与圆柱和矩形前缘相对位置以及雷诺数之间的关系，进一步揭示和探讨矩形前缘剪切流引发的圆柱绕流的雷诺数效应。1 风洞试验介绍试验在中南大学风洞实验室CSU-3风洞进行，该风洞为闭口环流式，试验段尺寸为 800 mm120 mm，风速控制范围为025 m/s，湍流度小于0.5%。本试验为刚性模型测压试验。模型如图 1 所示。圆柱直径为 D=60 mm，矩形板厚度为 H=48 mm。测压试验中，圆柱压力测点沿模型跨中位置周向布置，1号测点正对来流。在圆柱发生流动分离附近(813号测点、2732号测点)每8设置一个测点，其余位置每隔10设置一个测点。圆柱模型材质为铝合金，矩形板模型材质为 ABS 板。试验模型安装如图2所示。试验探究3组参数G/D，L/D和Re对圆柱气动特性的影响。其中，圆柱直径 D=60 mm，G为圆柱下表面距矩形板上表面的垂直流向距离，L为矩形板前缘距圆柱圆心的水平流向距离。具体含义图1模型尺寸及测点布置Fig.1Model size and measurement point arrangement图2试验照片Fig.2Test photo3640第 10 期李文舒，等：矩形前缘剪切流对圆柱气动和绕流特性影响规律的研究如图1所示。具体试验工况如表1所示。实验中自由来流风速 U为 1025 m/s，对应Re为41041105。模型表面的压力测压点平均风压系数Cp与脉动风压系数Cp的定义如下1314：Cp=(P-P)/0.5U2(1)Cp=Prms/0.5U2(2)式中：P为各测点风压的平均值；P为风洞静压；Prms为脉动风压的均方根值。由此亦可通过面积积分的形式计算模型总体的平均阻力系数Cd，平均升 力 系 数 Cl，脉 动 阻 力 系 数Cd和 脉 动 升 力系数Cl。2 结果与讨论2.1气动力系数2.1.1不同 G/D 和不同 L/D 条件下圆柱气动力的变化图3所示为圆柱和矩形板的相对位置对圆柱气动力特性的影响规律，图中的圆点为试验工况的位置。Cd的分布规律与矩形前缘形成的剪切层的形状相似。越靠近矩形板的上壁面，圆柱的Cd越小，即Cd随着G/D的减小而减小，这一现象与近壁面圆柱绕流中圆柱的阻力特性类似。当G/D较小时，Cd较小而Cl较大，且Cl随着G/D的增大而减小，当G/D增大到一定值时，圆柱会产生负升力。Cd总体变化不大，圆柱越靠近矩形前缘，其Cd越小。圆柱的Cl的变化规律与近壁面圆柱的类似，即Cl随着G/D的减小而减小，这一现象也与参考文献1517中G/D较小时圆柱的涡脱会受到抑制这一结论吻合。2.1.2Re对圆柱气动力的影响图4所示为当G/D=0.1时，雷诺数对圆柱气动力的影响。很明显，在本文研究的雷诺数范围内，当圆柱距离矩形板前缘很近(L/D=0.60)时，圆柱气动力的雷诺数效应不明显。由于近壁面的影响，圆柱在距离矩形前缘不同位置处的阻力系数都低于孤立圆柱的阻力系数(Cd=1.2)，当L/D2.02时，Re对Cd的影响很小，当L/D3.00时，Cd随着Re的增大而减小。当L/D=1.60和L/D=2.02时，升力系数的绝对值随着 Re的增大而增大，随着L/D的增大，Cl随着Re的增大而增大，并且会逐渐由负值恢复为正值，且 Re越大，零升力对应的 L/D表1模型测试工况Table 1Model test conditions参数G/DL/DRe工况0,0.1,0.2,0.3,0.40.60,1.60,2.02,3.00,3.56,4.40,5.804104,6104,8104,1105(a)Cd；(b)Cl；(c)Cd；(d)Cl图3气动力系数随圆柱和矩形板相对位置的变化Fig.3Aerodynamic coefficient varies with the relative position of the cylinder and rectangular plates3641铁 道 科 学 与 工 程 学 报2023 年 10月越小。2.2平均风压系数2.2.1G/D对Cp的影响图 5 以 L/D=1.60 和 L/D=4.40 为例，表明了在Re=1105的工况下，圆柱的平均风压系数随 G/D的变化情况。随着 G/D的增大，圆柱产生的较大负压区由上侧逐渐转移至下侧，最后平均风压实现基本对称分布，这与参考文献1819中近壁面圆柱绕流的规律相似。因此可以认为，圆柱发生湍流转捩的位置与矩形前缘剪切流打在圆柱上的位置有关。在L/D相同的条件下，圆柱的风压驻点随着G/D的减小而向下游移动，背压随着G/D的增大而增大。相较于圆柱下部的分离点区域，圆柱上部分离点区域的平均风压受 G/D 的影响较小，且随着G/D的增大，圆柱上部的负压极值点逐渐向上游移动。从图5的变化可以看出，在G/D较大的情况下，当 L/D 较小时，圆柱平均风压基本上呈对称分布，而随着 L/D 的增大，其对称性逐渐减 弱，因 此，下 一 步 探 究 L/D 对 圆 柱 流 动 的影响。2.2.2L/D对Cp的影响图6以G/D=0和G/D=0.2为例，表明了在Re=1105的工况下，圆柱的平均风压系数随L/D的变化情况。当 G/D=0时，在圆柱与矩形板上表面接触的位置(=270)，平均风压的突变程度在 L/D=0.60时最大，之后随着L/D的增大，突变程度先增大后减小。随着L/D的增大，风压驻点从圆柱下侧偏移至圆柱上侧，且驻点处的风压值逐渐减小，(a)Cd与Re和L/D的关系；(b)Cl与Re和L/D的关系图4当G/D=0.1时，雷诺数对圆柱气动力的影响Fig.4Effect of Reynolds numbers on cylinder aerodynamics when G/D=0.1(a)L/D=1.60；(b)L/D=4.40图5平均风压系数随G/D的变化Fig.5G/D dependence of Cp-mean3642第 10 期李文舒，等：矩形前缘剪切流对圆柱气动和绕流特性影响规律的研究并伴随有圆柱上部流动分离点逐渐向下游移动的现象。由图5和图6可知，圆柱在相对矩形板前缘不同位置产生不同气动特性的机理与剪切流打在圆柱上的位置有关。2.2.3Re对Cp的影响图 7 以 L/D=0.60 和 L/D=3.56 为例，表明了当G/D=0.1时，圆柱平均风压随 Re的变化情况。可以看出，在本文研究的雷诺数范围内，当圆柱距离矩形板前缘很近(L/D=0.60)时，圆柱风压的雷诺数效应不明显。随着L/D的增大，雷诺数主要影响圆柱的背压以及流动分离区的负压，随着Re增大，负压极值点向下游移动，负压极值增大，背压减小，从而导致分离泡的体积随着 Re增大而增大。有趣的是，在圆柱与矩形前缘特殊的相对位置处会出现雷诺数影响圆柱湍流转变的情况，即当G/D=0.10 和 L/D=3.00，G/D=0.10 和 L/D=3.56，G/D=0.10 和 L/D=4.40 时，Re 会影响湍流转变的位置，在这些工况下，随着雷诺数的增大，圆柱发生湍流转变的位置会由下侧转移至上侧，且不会出现由上侧转移至下侧的相反情况。2.3涡脱频率图8所示为圆柱涡脱频率随L/D变化的规律，其中横坐标代表圆柱涡脱频率，f为周期性流动频率，D为圆柱直径，U为来流速度，纵坐标为升力功率谱。从图8(a)8(c)可以看出，当G/D较小时，圆柱的涡脱会被抑制，升力功率谱没有明显的峰(a)G/D=0；(b)G/D=0.2图6平均风压系数随L/D的变化Fig.6L/D dependence of Cp-mean(a)L/D=0.60；(b)L/D=3.56图7平均风压系数随Re的变化Fig.7Reynolds number dependence of Cp-mean3643铁 道 科 学 与 工 程 学 报2023 年 10月值，这一现象与近壁面圆柱类似，这是由于当G/D很小时，壁面的存在阻碍了圆柱尾流区下部剪切层与外侧流动间的相互作用，抑制了尾流区下部大尺度涡的卷起，减少了剪切层内涡量，使下部涡旋得不到充分发展，从而涡脱受到抑制。当L/D=0.60 时，涡脱抑制的临界间隙比(G/D)cr=0.30；当L/D=3.56时，(G/D)cr=0.40；当L/D=5.80时，(G/D)cr0.40。因此，涡脱抑制的临界间隙比(G/D)cr会随着L/D的增大而增大。2.4瞬时风压2.4.1圆柱瞬时风压的双稳态现象图 9 所示为圆柱双稳态现象随 L/D 的变化规律，其中横坐标表示对时间的无量化，纵坐标为圆柱平均风压系数。双稳态现象只在G/D=0.4时发生，且随着 L/D 的增大，双稳态现象逐渐减弱。ALAM等13曾在交错双圆柱绕流试验中观察到，当P/D=1.12.3，=10(为双圆柱圆心连线和来流方向的夹角)时，下游圆柱会出现双稳态现象，并指出该现象与下游圆柱间隙侧表面分离泡的形成和破裂有关。对于本文出现的双稳态现象的机理将在后续的试验中进一步研究。(a)L/D=0.60；(b)L/D=3.56；(c)L/D=5.80图8圆柱涡脱频率随L/D和G/D的变化规律Fig.8Regularity of cylindrical vortex frequency with L/D and G/D(a)L/D=0.60；(b)L/D=1.60；(c)L/D=3.00；(d)L/D=4.40图9圆柱发生双稳态现象对应的风压时程Fig.9Wind pressure time range corresponding to the phenomenon of bistable state occurring in the cylinder3644第 10 期李文舒，等：矩形前缘剪切流对圆柱气动和绕流特性影响规律的研究2.4.2圆柱瞬时风压的转捩现象图 10 所示为在 L/D=4.40，G/D=0.20 的工况下，圆柱的瞬时风压随雷诺数的变化，左侧是风压时程图，右侧是对应的概率密度函数图。其中，实线代表高斯分布，虚线代表测点风压的概率密度分布，分别取各工况下平均风压负压极值和脉动风压最大的点进行分析，纵坐标与左侧图相对应，表示测点的瞬时风压，曲线与纵坐标围成的面积为1。通过分析，这些工况下发生转捩的位置都在9号测点附近(如图1所示)。在图10(c)中(Re=8104)，右侧对应的概率密度函数出现了双峰特征，表明圆柱出现在临界雷诺数区才会出现的时程转捩现象2021。若降低或增大来流雷诺数，如图10(a)，10(c)和10(d)，双峰现象会消失，非高斯分布特性也逐渐减弱，因此圆柱的转捩现象只在特定雷诺数且只在其相对于矩形板的特定位置时才(a)Re=4104；(b)Re=6104；(c)Re=8104；(d)Re=1105图10当L/D=4.40,G/D=0.20时，圆柱的风压时程和概率密度函数随雷诺数的变化Fig.10Wind pressure time history and probability density function of a cylinder change with Reynolds number when L/D=4.40,G/D=0.203645铁 道 科 学 与 工 程 学 报2023 年 10月会发生。由图10可知，随着雷诺数的增大，圆柱经历了由转捩开始到正在转捩再到转捩完成的一个过程。除了图 10所提到的L/D=4.40，G/D=0.20这一工况，在L/D=3.00，G/D=0.10；L/D=3.56，G/D=0.20；L/D=5.80，G/D=0.30时，同样会出现类似图10(c)的现象。可以看出，这些发生转捩现象的工况，其圆柱与矩形板相对位置的变化特点是L/D随着G/D的增大而增大，这与矩形板前缘产生的剪切流发展的形式类似。因此，从风压时程分布曲线可以进一步确定，由于矩形板前缘产生的剪切流在圆柱流动分离的位置产生了局部流动加速，从而使圆柱局部雷诺数增大，使得圆柱在亚临界雷诺数区间呈现出临界雷诺数区间的流动特性。在杜晓庆等22的双圆柱绕流相关文献中提到，在特定的双圆柱相对位置的情况下，由于上游圆柱形成的加速剪切层打在下游圆柱上，导致下游圆柱的局部位置产生较大负压，阻力系数突降。类比得到，矩形板与圆柱在特定的相对位置也会由于矩形板前缘剪切层加速，从而使圆柱在亚临界雷诺数区间产生临界雷诺数区间才会有的气动现象。3 结论1)圆柱在相对矩形板前缘不同位置时，会产生不同气动特性的机理，这与剪切流打在圆柱上的位置有关，圆柱相对矩形板处于不同位置时，会产生不同的风压分布、气动力分布以及涡脱特性。2)圆柱与矩形板处于特定的相对位置时，矩形前缘剪切流的加速侧打在圆柱上，会使圆柱局部雷诺数增大，从而在亚临界区产生临界区的气动特性，且该相对位置的特点是G/D随着L/D的增大而增大。3)当G/D=0.4时，圆柱会产生双稳态现象，且随着L/D的增大，双稳态现象逐渐减弱。参考文献：1张明超,王汉封,彭思.近平板圆柱体气动力雷诺数效应 J.水 动 力 学 研 究 与 进 展(A 辑),2018,33(5):666673.ZHANG Mingchao,WANG Hanfeng,PENG Si.Reynolds number effects on the aerodynamic forces of a circular cylinder near a flat plateJ.Chinese Journal of Hydrodynamics,2018,33(5):666673.2王汉封,邹超,王启文,等.均匀流中悬臂圆柱体气动力雷诺数效应J.湖南大学学报(自然科学版),2015,42(5):6571.WANG Hanfeng,ZOU Chao,WANG Qiwen,et al.Reynolds number effects on the aerodynamic forces of a cantilevered circular cylinder in uniform flowJ.Journal of Hunan University(Natural Sciences),2015,42(5):6571.3BEARMAN P W,ZDRAVKOVICH M M.Flow around a circular cylinder near a plane boundaryJ.Journal of Fluid Mechanics,1978,89(1):3347.4LIN W J,LIN Chang,HSIEH S C,et al.Flow characteristics around a circular cylinder placed horizontally above a plane boundaryJ.Journal of Engineering Mechanics,2009,135(7):697716.5NIEMANN H J,HLSCHER N.A review of recent experiments on the flow past circular cylindersJ.Journal of Wind Engineering and Industrial Aerodynamics,1990,33(1/2):197209.6MICHELIS T,KOTSONIS M.Interaction of an off-surface cylinder with separated flow from a bluff body leading edgeJ.Experimental Thermal and Fluid Science,2015,63:91105.7周磊,许慧,喜冠南.间隙比对近壁单圆柱绕流影响的PIV实验研究J.机械设计与制造,2020(11):2125.ZHOU Lei,XU Hui,XI Guannan.Effect of gap ratio on flow passed a cylinder near wall via PIV experimentJ.Machinery Design&Manufacture,2020(11):2125.8WANG X K,TAN S K.Near-wake flow characteristics of a circular cylinder close to a wallJ.Journal of Fluids&Structures,2008,24(5):605627.9TANIGUCHI S,MIYAKOSHI K.Fluctuating fluid forces acting on a circular cylinder and interference with a plane wallJ.Experiments in Fluids,1990,9(4):197204.10 PRICE S J,SUMNER D,SMITH J G,et al.Flow visualization around a circular cylinder near to a plane wallJ.Journal of Fluids and Structures,2002,16(2):175191.11 WU Gefei,LIN Weiqun,DU Xiaoqing,et al.On the flip-flopping phenomenon of two side-by-side circular cylinders at a high subcritical Reynolds number of 3646第 10 期李文舒，等：矩形前缘剪切流对圆柱气动和绕流特性影响规律的研究1.4105J.Physics of Fluids,2020,32(9):181192.12 GU Zhifu,SUN Tianfeng.On interference between two circular cylinders in staggered arrangement at high subcritical Reynolds numbersJ.Journal of Wind Engineering and Industrial Aerodynamics,1999,80(3):287309.13 ALAM M M,SAKAMOTO H,ZHOU Y.Determination of flow configurations and fluid forces acting on two staggered circular cylinders of equal diameter in cross-flowJ.Journal of Fluids and Structures,2005,21(4):363394.14 黄伯城.临界雷诺数区光滑圆柱气动力特性研究D.石家庄:石家庄铁道大学,2018.HUANG Bocheng.Study on aerodynamic characteristics of smooth cylinder in critical Reynolds number regionD.Shijiazhuang:Shijiazhuang Tiedao University,2018.15 邹云峰,何旭辉,邹思敏,等.一种测试移动列车桥梁系统气动特性的U形滑道和方法J.铁道科学与工程学报,2018,15(7):17471753.ZOU Yunfeng,HE Xuhui,ZOU Simin,et al.A U-shape launching ramp and method for aerodynamic characteristic testing of running train-bridge systemJ.Journal of Railway Science and Engineering,2018,15(7):17471753.16 HE Guosheng,WANG Jinjun,PAN Chong,et al.Vortex dynamics for flow over a circular cylinder in proximity to a wallJ.Journal of Fluid Mechanics,2017,812:698720.17 BURESTI G,LANCIOTTI A.Mean and fluctuating forces on a circular cylinder in cross-flow near a plane surfaceJ.Journal of Wind Engineering and Industrial Aerodynamics,1992,41(1/2/3):639650.18 LEI C,CHENG L,KAVANAGH K.Re-examination of the effect of a plane boundary on force and vortex shedding of a circular cylinderJ.Journal of Wind Engineering&Industrial Aerodynamics,1999,80(3):263286.19 MARYAMI R,ALIHAN SHOWKAT ALI S A,AZARPEYVAND M,et al.Turbulent flow interaction with a circular cylinderJ.Physics of Fluids,2020,32(1):015105.20 DESAI A,MITTAL S,MITTAL S.Experimental investigation of vortex shedding past a circular cylinder in the high subcritical regimeJ.Physics of Fluids,2020,32(1):014105.21 WU Gefei,DU Xiaoqing,WANG Yuliang.LES of flow around two staggered circular cylinders at a high subcritical Reynolds number of 1.4105J.Journal of Wind Engineering and Industrial Aerodynamics,2020,196:104044.22 杜晓庆,林伟群,施春林,等.高雷诺数下并列双圆柱绕流的大涡模拟J.哈尔滨工业大学学报,2019,51(6):193200.DU Xiaoqing,LIN Weiqun,SHI Chunlin,et al.Large eddy simulation of flow around two side-by-side circular cylinders at a high Reynolds numberJ.Journal of Harbin Institute of Technology,2019,51(6):193200.(编辑 杨碧峰)3647