圆柱亚临界临界超临界流动特性的大涡亚格子模型比较研究.pdfVIP

第１４卷第４期２０１６年８月 动力 学 与控 制 学报 Ｖｏｌ．１４Ｎｏ．４ １６７２６５５３／２０１６／１４ ／３７１８ Ａｕｇ．２０１６ ⑷ ＪＯＵＲＮＡＬＯＦＤＹＮＡＭＩＣＳＡＮＤＣＯＮＴＲＯＬ 圆柱亚临界临界超临界流动特性的大涡亚格子 模型比较研究  祝志文　邓燕华 （湖南大学土木工程学院，长沙　４１００８２） 摘要　为获取圆柱在亚临界、临界和超临界区内的气动特征，基于大涡模拟并分别采用Ｓｍａｇｏｒｉｎｓｋｙ亚格子 ４ ５ 尺度模型和动力亚格子尺度模型计算了Ｒｅ＝４．１×１０～８．２×１０内的绕流场．将获得的圆柱定常与非定常 气动特性与相关文献进行了对比．研究表明，两种亚格子模型均能给出圆柱升力系数平均值、ＲＭＳ值和涡脱 Ｓｔ数的合理估计，能预测阻力系数的下降．但两种亚格子模型均无法准确预测临界区阻力系数的快速下降、 不同Ｒｅ数下压力系数平均和ＲＭＳ值在圆柱表面的分布．在临界和超临界区，Ｓｍａｇｏｒｉｎｓｋｙ亚格子尺度模型 可能高估了亚格子湍流粘性；相反，动态亚格子模型可能在流场部分区域明显低估了亚格子湍流粘性． 关键词　Ｓｍａｇｏｒｉｎｓｋｙ模型，　动态亚格子模型，　气动特性，　临界区 ＤＯＩ：　１０．６０５２／１６７２６５５３２０１５８０ 且有时脉动升力根方差（ＲＭＳ）很大．但在临界 Ｒｅ 引言 数区及附近，文献［１］和［３］等报道了较大非零平 土木工程结构如主缆、桥墩、拉索和吊杆，以及 均升力的出现，引起了土木工程领域的广泛关注． 电视塔、桅杆、烟囱和冷却塔，甚至一些超高层建筑 都具有圆形截面特征．这些工程结构不仅会受到流 体的静荷载作用，产生内力和变形，甚至可能会在 周期性流体动力作用出现大幅振动现象．虽然圆柱 形截面外形简单，但其绕流形态如分离点位置（ ） θ ｓ 等随Ｒｅ数的变化而不断变化（见图１），其气动力 ［１］ 和漩涡脱落表现出很强的Ｒｅ数效应 ，其流体动 力学特性远比具有尖锐棱角的钝体截面复杂．一般 ５ ５ 图１　圆柱流动特性的Ｒｅ数分区 认为，在２×１０＜Ｒｅ＜４～５×１０的临界区内，圆 柱的阻力系数会出现快速的大幅下降，称之为阻力 Ｆｉｇ．１　Ｒｅｎｕｍｂｅｒｐａｒｔｉｔｉｏｎｏｎｆｌｏｗｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆｃｉｒｃｕｌａｒｃｙｌｉｎｄｅｒ 危机（ｄｒａｇｃｒｉｓｉｓ，图１），此后阻力系数又缓慢增 大，因而会导致作用在圆形结构上的顺风向荷载随 与试验相比，数值方法能非常方便地设定和处 来流风速急剧变化．虽然此时Ｒｅ数比较高，但对大 理来流和边界条件，不存在模型支架和测量系统的 尺度的工程结构是完全可能达到的．目前试验研究 干扰，因堵塞度可显著减小，可明显减少洞壁的干 报道的Ｒｅ数已覆盖了从低Ｒｅ数、亚临界、临界到 扰，因而在土木工程中得到了越来越多的应 超临界和高超临界的广泛范围，但由于不同研究的 用［４－５］．本文数值方法采用大涡模拟（ＬＥＳ），该方