类方柱的气动性能和流场特性

杜晓庆 刘延泰 董浩天 施定军

摘   要：以优化方柱气动外形为目的，采用大涡模拟方法，在雷诺数为22 000时，研究了采用不同角、边平面形态及其组合（包括尖角、圆角、直边、凸边和凹边等组合）的类方柱的气动性能和流场特性，重点分析了凸边圆角柱的气动性能随边部曲率半径的变化规律，并通过分析流场结构揭示了平面形态修正对类方柱气动性能影响的作用机理. 研究发现，角、边形状修正可显著改变绕流场结构，最终影响类方柱的气动性能：尖角柱体圆角化、直边柱体凸边化能显著降低气动力、负压区强度和涡脱强度，并伴随斯托罗哈数上升;而直边柱体凹边化后气动性能变化趋势相反;不同角边形状的组合中，凸边和圆角的组合可导致分离剪切层更紧贴柱体壁面，上、下侧回流区范围变小，尾流回流区长度增大，涡脱强度减弱，气动力下降. 对凸边圆角柱的进一步研究表明，柱体的气动性能对边部曲率半径非常敏感;并存在一最优曲率半径，此时凸边圆角柱的平均和脉动气动力最小，斯托罗哈数最大，尾流长度最长，涡脱强度最低.

关键词：方柱绕流;大涡模拟;气动优化;高雷诺数;截面修正

中图分类号：TU973                                  文献标志码：A

Abstract：To optimize the shape of square cylinder， Large-eddy Simulation is applied in the study of flow past square cylinders at a high Reynolds Number Re of 22 000. The cylinders have various combinations of sharp/round corners and flat/concave/bulging edges. Particularly，the aerodynamic characteristics of cylinders with a combination of round corners and bulge edges are studied with the changing curvature radius of the edge. Besides， the flow patterns are discussed in detail to reveal the mechanism of aerodynamic optimization. The results show that corner and side modifications can significantly alter the structure of flow around cylinders and influence the aerodynamic performance of cylinders. In comparison to the concave modification， the round modification and bulge modification can reduce the aerodynamic forces， overall surface pressures and vortex shedding intensity of cylinders considerably. Among cylinders with various combinations of shape modifications， the cylinder with rounded corner and bulge edge leads to a thinner shear layer and a smaller recirculation region at the upper and lower sides of cylinder. Besides，it also causes a longer wake length， a weaker vortex shedding and reduces the aerodynamic forces. A further research shows that the aerodynamics of the round-bulge cylinder is extremely sensitive to the curvature radius of the edge. There exists an optimal curvature radius， which results in the lowest aerodynamic forces，the highest Strouhal number， the longest length of the recirculation region， and the weakest vortex shedding intensity.

Key words：flow past square cylinders;Large-eddy Simulation（LES）;aerodynamic optimization;high Reynolds number;cross-section modification

大長细比方形截面柱体结构在土木工程中有广泛应用，如超高层建筑和大跨度缆索承重桥的索塔，其抗风问题是结构设计的关键[1-4].气动控制措施是一种经济、有效的抗风设计方法，适当的气动措施可有效减小风荷载和风致响应，因而方柱的气动性能优化得到研究者的广泛关注[5-8].

超高层建筑的气动控制措施包括截面的角部处理、截面形状改变、沿高度的形状或尺寸变化等[9-13]. Tanaka等[14-16]对角部处理、截面形状的优化和沿高度的形状或尺寸变化的模型进行了风洞试验，对比分析了不同气动外形柱体的倾覆力矩系数、风力响应等结果.

针对方形截面柱体，常用的平面形态气动优化措施为角部修正，包括圆角、倒角、切角等形式，这类气动措施可有效优化方柱结构的受力与风致振动问题[17]. Tamura及其合作者[18-20]对标准方柱、圆角、切角方柱进行了数值模拟和风洞试验，结果表明，圆角和切角措施可以削弱涡脱强度，减小平均阻力和脉动升力系数. 文献[21-26]中研究了圆角、切角、凹角等角部修正对方柱风致响应的影响，研究表明，3种角部措施中，圆角化方柱能最有效地抑制柱体的横风向响应.

对于圆角化方柱，圆角率会对柱体气动性能产生很大影响. Delany等[27]通过风洞试验研究R/B =

1/50、1/6、1/3（R为模型圆角半径，B为模型名义边长） 3种圆角率方柱气动力的雷诺数效应.  Carassale等[28-29]针对R/B=0、1/15、2/15等3种圆角率的方柱模型，通过风洞试验研究了来流湍流度和雷诺数对其气动性能的影响. Zhang等[30]则通过数值模拟研究了低雷诺数（Re=1 000）下R/B=0、0.125、0.250、0.375、0.5等5种圆角率方柱的流场特性和气动性能. 杜晓庆等[31]对R/B = 1/7的圆角方柱进行大涡模拟研究，解释了圆角化对方柱气动性能影响的流场机理.

從以往研究看，方柱平面形态的修正局限在角部处理气动措施，但随着超高层建筑的大量建设，出现了不少角部修正与凸边、凹边等形状相结合的类方柱结构，而目前尚未见到研究此类方柱气动性能和流场特性的文献. 基于上述考虑，本文采用大涡模拟方法，在雷诺数为2.2×104时，以角（尖角、圆角）、边（直边、凸边和凹边）形状修正的类方柱为研究对象，研究了16种类方柱的气动性能和流场特性，以及凸边圆角柱的气动性能随边部曲率半径的变化规律，分析了角、边形状的改变对柱体气动性能的影响，并从流场角度探讨了形状修正对方柱气动性能的作用机理.

1   数值方法和计算模型

1.1   流动控制方程和亚格子模型

在大涡模拟（LES）方法中，大尺度涡通过滤波后的Navier-Stokes方程直接求解，而小尺度的涡则采用亚格子尺度模型（SGS）模拟. 经过滤波函数的滤波，可得到大尺度涡的不可压缩Navier-Stokes方程：

本文采用的WALE常数Cω是0.325，此为Fluent默认值，该取值下可取得比较满意的结果. 计算采用SIMPLEC格式求解压力速度耦合方程组，空间离散采用中心差分格式，时间离散采用二阶全隐格式.

1.2   计算模型和计算域

图1为以凸边圆角柱为例的计算模型示意图. 均匀来流速度为U，来流湍流度为0，所有模型的截面面积相等，即面积S = B2（其中B为计算模型的名义边长），风向角为0°，雷诺数为2.2×104（根据模型名义边长B和来流风速U计算得到）. 模型的圆角半径定义为R，边部曲率半径定义为K.

1/100和1/7（文中分别称为尖角和圆角）;考虑了3种边部形状，即直边、凸边、凹边. 图2是本文所选取的6种截面类方柱，分别命名为：直边尖角柱（标准方柱）、凹边尖角柱、凸边尖角柱、直边圆角柱、凹边圆角柱和凸边圆角柱，边部曲率半径均为K/B = 2.5. 图3是本文所选取的部分不同曲率半径的凸边圆角柱模型示意图，本文考虑的曲率半径K/B = 0.564（圆柱）、0.7、0.8、1.0、1.2、1.25、1.3、1.35、1.4、1.5、2.5、∞（直边圆角柱）.

计算模型的整体网格和角部局部网格放大图见图4. 由图4（a）可见，入口边界圆弧顶点距离方柱中心为20B，方柱中心距离出口边界距离为30B，横风向计算域宽度为40B，计算域高度（方柱展向长度H）为4B，阻塞率为2.5%. 计算域采用速度入口边界条件，自由出口边界，两侧采用对称边界条件，展向采用周期性边界条件，方柱表面为无滑移壁面边界条件.

2   计算参数和结果验证

本文以标准方柱为研究对象，在文献[31]的基础上进行计算模型的结果验证，结果如表1和图5所示. 文献[31]中，网格验证所采取模型展向长度为2B，为提高网格精度，本文以展向长度4B进行模型验证，且验证计算时间步和展向长度的参数范围更大.

与文献[31]网格验证部分相一致，分别考察了不同周向网格数量、计算时间步、展向长度和展向网格尺度对方柱气动性能计算结果的影响. 本文方柱周向网格数的取值为120、200和320;无量纲时间步长Δt* = ΔtU0 /B（其中Δt为有量纲时间步，U0为来流风速），为0.005、0.010、0.025和0.030;展向长度H为1B、2B、4B和6B;展向网格尺度为0.15B、0.1B和0.075B.

由表1和图5可得，周向网格数量、计算时间步和展向网格尺度对计算结果的影响与文献[31]中的分析相似，故本文所采取的周向网格数量、计算时间步和展向网格尺度与文献[31]相一致，即周向网格数为200，无量纲时间步为0.025，展向网格尺度为0.1B. 在工况A2的基础上，选取展向长度为1B、2B、4B、6B进一步研究展向长度的影响. 比较表1中本文计算的各个工况，方柱展向长度为1B时，得到的气动力系数整体均偏小，对比长度为2B、4B和6B的计算结果，可以发现：随着展向长度的增大，平均阻力系数增大，脉动阻力系数和St数变化不大，其中展向4B和6B时平均阻力系数与文献中的风洞试验结果和大涡模拟结果更为相近.

综合考虑结果精度和计算资源，下文工况所选取的计算模型参数均与工况A3相近，即展向长度为4B，展向网格尺寸为0.1B，周向网格数为200，角部网格加密;近壁面最小网格厚度为0.001B，壁面附近沿径向的网格增长率为1.06，近壁面y+≈1，满足大涡模拟对近壁面的网格精度要求;无量纲时间步Δt*为0.025.

3   角边形状的影响

3.1   气动力系数和St数

圖6给出了6种柱体的气动力系数和斯托罗哈数数的对比情况. 总体上看，与标准方柱相比，只有凹边尖角柱的气动力系数升高，其余柱体的气动力系数均有一定程度的降低，凸边圆角柱的气动力系数均最小.

由图6（a）可知，与标准方柱相比，凹边尖角柱的脉动阻力系数升高，其余柱体的脉动阻力系数降低，且凸边圆角柱脉动阻力系数最小;与角部为尖角的柱体相比，圆角化柱体的脉动阻力系数显著降低;与直边柱体相比，凸边柱体的脉动阻力系数升高，凹边柱体降低. 分析图6（b）和6（c）可得，角、边形状对脉动阻力和脉动升力系数的影响跟平均阻力系数的变化趋势一样.

从图6（d）可见，角、边形状对斯托罗哈数的影响与3种气动力系数的变化趋势相反. 与标准方柱相比，凹边尖角柱的斯托罗哈数降低，其余柱体升高，且凸边圆角柱的斯托罗哈数最大，即6种柱体中凸边圆角柱的涡脱频率最大.

3.2   平均和脉动风压系数

图7分别为6种柱体的表面平均风压和脉动风压系数分布图. 总体上看，各柱体表面风压系数分布曲线整体趋势相近，但数值差异较大，特别是脉动风压系数.

由图7（a）可见，迎风面（ab）上，各柱体平均风压系数基本一致;在上、下侧面（bc和ad）和背风面（cd），与角部为尖角的柱体相比，圆角化柱体平均风压系数绝对值均减小，3种圆角化柱体平均风压系数相差不大;对于3种尖角柱体，不同边部形状下，各柱体的平均风压系数在背风面角点（c点和d点）出现一定的差异，且凹边尖角柱的角点负压绝对值最大，其余位置则较为接近. 分析图7（b）可知，角、边形状对柱体脉动风压系数的影响与平均风压系数的变化趋势基本一致.

3.3   平均风压场

图8给出了0°风向角下6种柱体的时间平均风压场. 由图8可见，与直边柱体相比，凸边柱体上、下侧面及尾流区负压强度减弱，负压极值减小，这是凸边柱体平均阻力系数较小的主要原因;凹边柱体的平均风压场变化趋势与凸边柱体相反;尾流区风压场的差异是造成各柱体平均阻力系数不同的主要原因. 此外，与角部为尖角的柱体相比，圆角化柱体的上、下侧面和尾流区负压强度大大减小.

从总体上看，6种柱体中，凸边圆角柱尾流区域负压绝对值最小，可见凸边修正和圆角措施的组合能最有效地减弱柱体周围平均风压场.

3.4   平均流态结构

图9是6种柱体绕流场的时间平均流线图. 总体上看，气流在各柱体迎风面角点（a点和b点）发生分离，并在柱体上、下侧和尾流中形成6个回流区，侧面和尾流的回流区在方柱的背风面角点附近连通.

由图9可见，与直边柱体相比，凸边柱体上、下侧回流区范围减小，气流分离后的剪切层更贴近柱体侧面，而凹边柱体变化趋势相反;即边部形状的改变主要影响柱体上、下侧边的流场结构. 与角部为尖角的柱体相比，圆角化柱体上、下侧回流区范围显著减小，剪切层更贴近柱体侧面，尾流长度增大，宽度减小.

从总体上看，6种柱体中，凸边圆角柱上、下侧回流区范围最小，尾流长度最大，尾流宽度最小，这导致了凸边圆角柱受到最小的平均阻力系数，尾流的差异是造成各柱体的平均阻力系数不同的主要原因.

4   边部曲率半径的影响

由上文结果可知，凸边修正和圆角措施的组合能最有效改善柱体的气动性能，本节以凸边圆角柱为研究对象，进一步探讨边部曲率半径对凸边圆角柱（角部率R/B=1/7）气动性能和流场特性的影响.

4.1   气动力系数和St数

图10给出了不同边部曲率半径下凸边圆角柱的气动力系数和斯托罗哈数. 注意各柱体气动力系数均采用名义边长B无量纲化.

从图10（a）可见，随着边部曲率半径K/B的增大，柱体的平均阻力系数先减小后增大，在K/B=1.0时取得最小值，当曲率半径大于1.5时，K/B的改变对平均阻力系数的影响较小. 与直边圆角柱和圆柱相比，不同边部曲率半径下凸边圆角柱的平均阻力系数均有不同程度的降低，即圆角措施和凸边修正的组合可以有效降低柱体的平均阻力系数. 分析图10（d）可知，边部曲率半径对斯托罗哈数的影响与平均阻力系数的变化趋势基本相反，在K/B=1.0时斯托罗哈数最大，即涡脱频率最大.

由图10（b）可知，随着边部曲率半径的增大，柱体的脉动阻力系数经历了两次先减小再增大的变化，在K/B=1.0时第1次减小到极小值，在K/B=1.35时第1次增大到极大值，总的来看，在K/B = 1.0时柱体的脉动阻力系数最小. 分析图10（c）可得，边部曲率半径对脉动升力系数的影响与脉动阻力系数的变化趋势基本一致，在K/B = 1.0时柱体脉动升力系数最小.

4.2   平均和脉动风压系数

图11分别为边部曲率半径K/B=0.564（圆柱）、1.0、1.35、2.5、∞（直边圆角柱）时凸边圆角柱的表面平均风压系数和脉动风压系数分布图. 从图中可见，所有柱体中，当K/B=1.0时，柱体的平均和脉动风压系数均小于其余柱体.

由图11（a）可知，迎风面上，各柱体平均风压系数基本一致;迎风面角点（a点和b点）至上、下侧边中点部分，圆柱和K/B=2.5的凸边圆角柱负压较小，K/B=1.35的柱体负压最大;上、下侧边中点至背风面，K/B=1.0的柱体平均风压系数最小，圆柱的负压最大，这导致K/B=1.0的柱体出现最小的平均阻力系数，圆柱平均阻力系数最大. 值得注意的是，K/B=1.35凸边圆角柱的角部分离点附近的负压绝对值较大，从平均风速场（本文未给出）可知，这是由于该柱体气流分离点附近的局部风速更大造成的.

[5]     SOHANKAR A，DAVIDSON L，NORBERG C. Large eddy simulation of flow past a square cylinder：comparison of different subgrid scale models [J]. Journal of Fluids Engineering，2000，122（1）：39—47.

[6]     CAO Y，TAMURA T. Large-eddy simulations of flow past a square cylinder using structured and unstructured grids [J]. Computers and Fluids，2016，137：36—54.

[7]     NISHIMURA H，TANIIKE Y. Fluctuating wind forces on a stationary two dimension square prism[C]// Proceedings of 16th National Symposium on Wind Engineering. Tokyo：Springer，2000：255—260.

[8]   TAMURA T，ONO Y. LES analysis on aeroelastic instability of prisms in turbulent flow [J]. Journal of Wind Engineering and Industrial Aerodynamics，2003，91（12）：1827—1846.

[9]   MOONEGHI M A，KARGARMOAKHAR R. Aerodynamic mitigation and shape optimization of buildings：review [J]. Journal of Building Engineering，2016，6：225—235.

[10]  SHARMA A，MITTAL H，GAIROLA A. Mitigation of wind load on tall buildings through aerodynamic modifications：Review [J]. Journal of Building Engineering，2018，18：180—194.

[11]  趙昕，林祯杉，孙华华. 超高层建筑形态空气动力学优化方法应用综述[J]. 结构工程师，2011，27（3）：133—139.

ZHAO X. LIN Z B，SUN H H. Aerodynamic modification methods for high-rise buildings[J]. Structural Engineers，2011，27（3）：133—139. （In Chinese）

[12]   柯世堂，王浩. 超高层连体建筑风荷载干扰效应大涡模拟研究[J]. 湖南大学学报（自然科学版），2017，44（5）：53—62.

KE S T，WANG H. Large eddy simulation investigation of wind loads and interference effects on ultra high-rise connecting buildings[J].Journal of Hunan University（Natural Sciences），2017，44（5）：53—62. （In Chinese）

[13]   杜晓庆，吴葛菲，王玉梁，等. 尾流干扰下下游圆柱气动性能的流场机理[J]. 湖南大学学报（自然科学版），2019，46（7）：104—112.

DU X Q，WU G F，WANG Y L，et al. Mechanism of wake interference effects on aerodynamic characteristics of a downstream circular cylinder[J]. Journal of Hunan University（Natural Sciences），2019，46（7）：104—112. （In Chinese）

[14]  TAMURA Y，TANAKA H，OHTAKE K，et al. Aerodynamic characteristics of tall building models with various unconventional configurations[C]// ASCE Structures Congress 2010. Orlando，Florida，United States，2010：3104—3113.

[15]  TANAKA H，TAMURA Y，OHTAKE K. Experimental investigation of aerodynamic forces and wind pressures acting on tall buildings with various unconventional configurations[J]. Journal of Wind Engineering & Industrial Aerodynamics，2012，107/108（8）：179—191.

[16]  TANAKA H，TAMURA Y，OHTAKE K，et al. Aerodynamic and flow characteristics of tall buildings with various unconventional configurations[J]. International Journal of High-Rise Buildings，2013，2（3）：213—228.

[17]  顾明，张正维，全涌. 降低超高层建筑横风向响应气动措施研究进展[J]. 同济大学学报（自然科学版），2013，41（3）：317—323.

GU M，ZHANG Z W，QUAN Y. Aerodynamic measures for mitigation of across-wind responses of super tall buildings：State of the art[J]. Journal of Tongji University，2013，41（3）：317—323. （In Chinese）

[18]  TAMURA T，MIYAGI T，KITAGISHI T. Numerical prediction of unsteady pressures on a square cylinder with various corner shapes[J]. Journal of Wind Engineering and Industrial Aerodynamics，1998，74（2）：531—542.

[19]  CAO Y，TAMURA T. Supercritical flows past a square cylinder with rounded corners[J]. Physics of Fluids，2017，29（8）：085110.

[20]  CAO Y，TAMURA T. Shear effects on flows past a square cylinder with rounded corners at Re=2.2×104[J]. Journal of Wind Engineering and Industrial Aerodynamics，2018，174：119—132.

[21]  KWOK K C S，WILHELM P A，WILKIE B G. Effect of edge configuration on wind-induced response of tall buildings[J]. Engineering Structures，1988，10（2）：135—140.

[22]  KAWAI H. Effect of corner modifications on aeroelastic instabilities of tall buildings[J]. Journal of Wind Engineering and Industrial Aerodynamics，1998，74/76（98）：719—729.

[23]  HU J C，ZHOU Y，DALTON C. Effects of the corner radius on the near wake of a square prism[J]. Experiments in Fluids，2006，40（1）：106—118.

[24]  DEY P，DAS AKR. A numerical study on effect of corner radius and Reynolds number on fluid flow over a square cylinder[J]. Sadhana，2017，42（7）：1155—1165.

[25]  ELSHAER A，BITSUAMLAK G，DAMATTY A E. Enhancing wind performance of tall buildings using corner aerodynamic optimization[J]. Engineering Structures，2017，136：133—148.

[26]  HINSBERG N P V，SCHEWE G，JACOBS M. Experiments on the aerodynamic behaviour of square cylinders with rounded corners at Reynolds numbers up to 12 million[J]. Journal of Fluids & Structures，2017，74：214—233.

[27]  DELANY N K，SORENSEN N E. Low-speed drag of cylinders of various shapes[R]. United States：NACA Technical Note，1953：NACA-TN-3038.

[28]  CARASSALE L，FREDA A，MARR？魬-BRUNENGHI M. Effects of free-stream turbulence and corner shape on the galloping instability of square cylinders[J]. Journal of Wind Engineering and Industrial Aerodynamics，2013，123：274—280.

[29]  CARASSALE L，FREDA A，MARR？魬-BRUNENGHI M. Experimental investigation on the aerodynamic behavior of square cylinders with rounded corners[J]. Journal of Fluids and Structures，2014，44（7）：195—204.

[30]  ZHANG W，SAMTANEY R. Low- Re flow past an isolated cylinder with rounded corners[J]. Computers and Fluids，2016，136：384—401.

[31]  杜曉庆，田新新，马文勇，等. 圆角化对方柱气动性能影响的流场机理[J]. 力学学报，2018，50（5）：1013—1023.

DU X Q，TIAN X X，MA W Y，et al. Effects of rounded corner on aerodynamics of square cylinders and its flow mechanisms[J]. Chinses Journal of Theoretical and Applied Mechanics，2018，50（5）：1013—1023. （In Chinese）