论文摘要
钝体绕流时的旋涡脱落是自然界中最普遍的一种物理现象，包含着许多复杂的流动机理。旋涡脱落在钝体表面形成周期性脉动作用力，当钝体固定方式为弹性支撑或允许发生形变时，脉动力将引发钝体产生周期性振动，钝体振动又进一步影响尾流旋涡脱落形态，这种流体与结构物之间相互作用的问题即为流致振动（Flow-induced motion, FIM）。在许多与流体有关的机械与工程中，流致振动是一个涉及安全性的重大问题。涡致振动水生清洁能源（VIVACE）系统将流致振动与能源利用相结合，使这种潜在的破坏性现象在清洁能源领域得到有效利用。为了提高VIVACE系统的能量密度，本文采用被动湍流控制（PTC）和多柱体结构改善系统能量收集性能。另外，随着近年来海洋工程、风工程、航天航空和核工程的飞速发展，多柱体结构流致振动问题受到广泛关注，已成为流体力学理论研究和应用领域的重要课题之一，尚存在大量的问题需要进一步深入研究。因此，对于多柱体流致振动的研究具有非常重要的学术意义与应用价值。本文针对被动湍流控制下的多圆柱流致振动问题进行数值分析和实验研究，详细分析了多圆柱流致振动特性和尾流旋涡形态，探讨了各种因素对多柱体旋涡脱落的影响，揭示了多柱体流致振动产生机理，探索改善多柱体流致振动和提高涡致振动水生清洁能源系统能量密度的方法。采用Spalart-Allmaras湍流模型，求解非稳态雷诺平均纳维-斯托克斯方程组（RANS）获得多圆柱流致振动的数值解，并在美国密歇根大学海洋可再生能源实验室（MRELab）完成相应实验测试，得到了以下结论。首先，针对PTC单圆柱流致振动系统展开数值计算和实验研究，获得不同来流速度下圆柱振幅、频率和近尾迹流场特征，建立了描述柱体流致振动能量转换装置发电功率的数学模型，分析了最优能量收集速度范围。当雷诺数（Re）在3×104≤Re≤1.3×105范围内时，PTC单圆柱流致振动RANS数值结果与实验测试极为吻合，能够清晰观察到涡致振动（VIV）初始分支、VIV上部分支和驰振。柱体每个振动周期产生的脱体旋涡数量随着Re数的上升而增加。PTC单圆柱VIVACE系统在整个测试Re数范围内都能有效进行能量收集，模拟结果表明，在VIV上支和驰振起始区域具有最优能量收集效率。其次，采用GGI（General Grid Interface）界面结合拓扑网格变形技术，首次对串行排列的弹性支撑PTC双圆柱结构涡致振动和驰振进行数值计算，并在MRELab完成了实验测试。研究表明PTC双圆柱的流致振动振幅和频率响应数值模拟与实验结果趋势一致，两个圆柱均清晰观察到VIV初支和上支、VIV-驰振过渡分支及驰振。当Re=30000时，上游圆柱尾流旋涡形态为典型的2S模式，下游圆柱受到上游圆柱脱体旋涡的影响，运动几乎受到抑制；当Re=59229时，上游圆柱尾涡形态在2P和2P+2S之间切换，下游圆柱尾涡形态为2P模式。在PTC双圆柱发生驰振时，数值模拟获得高达3.5D的振幅（D为圆柱直径），实验测试中的圆柱撞击水槽安全停止阀，达到实验水槽所允许的振幅极限。然后，针对串行排列的三个PTC圆柱，研究了各个圆柱流致振动规律，首次通过数值模拟可视化分析了涡致振动过渡到驰振的原因及驰振产生机理。数值模拟和实验研究均表明，在光滑圆柱表面采取PTC措施，能有效刺激边界层分离和旋涡脱落，扩大柱体流致振动同步范围，使高Re数下数值结果与实验测试相吻合。在圆柱振动进入VIV上部分支以后，随着来流速度增加，线性粘滞阻尼模型能有效模拟高Re数PTC圆柱流致振动。圆柱流致振动各个分支间的转移往往伴随尾流旋涡形态的改变，在PTC圆柱由VIV过渡到驰振时，流致振动产生机理也随之发生变化。PTC圆柱驰振的产生是由升力不稳定性引起的，在脱体旋涡形成过程中，圆柱上下侧的分离剪切层被拉伸到几乎垂直于来流方向的位置，造成圆柱上下侧出现较大压差，从而使圆柱受到的升力急剧增大产生驰振。最后，通过数值模拟和实验分析对串行排列的PTC四圆柱的流致振动特性开展了研究。研究发现，在串行排列的四圆柱系统中引入PTC，是实现数值结果与实验测试相互匹配的关键因素。在3×104≤Re≤1.05×105范围内，四个圆柱均观察到了VIV初支、VIV上支、VIV-驰振过渡分支和驰振。PTC圆柱振动由VIV上支向驰振过渡时，实验发现水槽自由液面会造成圆柱振幅突然降低，使柱体驰振出现振幅波动。在Re=50000时，来自上游第1个圆柱的脱体旋涡使第2个圆柱受到的流体作用力变得不稳定，极大削弱圆柱横向升力，从而导致第2个圆柱出现低振幅和低频率振动现象。由于上游PTC圆柱产生的脱体旋涡对下游圆柱运动产生明显影响，使各个圆柱振动频率相互关联，第2/3/4圆柱的振动频率随流速的变化趋势相似。上游PTC圆柱每个振动周期旋涡产生数量往往多于下游柱体，尾流旋涡形态也更为复杂。