低柱体雷诺数下柱体上游薄层水流马蹄涡特征研究

马蹄涡是水下柱体基部产生局部侵蚀的主要动力,目前对大雷诺数深水条件的马蹄涡特征研究广泛,而在低柱体雷诺数浅薄层水流条件下因技术条件限制,实验手段难以准确捕获马蹄涡。为准确掌握在该水流条件下马蹄涡的特征,本研究构建了高分辨率高频的粒子图像测速系统,系统测量了6组浅薄层水流条件下柱体前端的瞬时流场。通过分析平均流场特征进而提取出流动分离点,采用旋转强度的方法来识别马蹄涡进而提取马蹄涡位置及强度,使用Oseen涡和纯剪切的控制方程叠加模拟马蹄涡进而计算出马蹄涡的半径。结果表明,在低柱体雷诺数水流条件下(ReDlt;5 000),随着柱体雷诺数的增加流动分离点急剧向下游移动,同时马蹄涡也急剧向柱体端和床面靠近,其半径减小而旋转强度增加。在浅薄层水流条件下,当柱体直径一定时,随着水深的增加,流动分离点向上游运动,同时马蹄涡向远离柱体端和朝水面运动,马蹄涡的半径增加,且其各项参数显著大于明渠水流条件下。随后,结合已有工作,归纳出分离点、马蹄涡特征随柱体雷诺数增加而经历的不同阶段:当5 000lt;ReDlt;8 000后流动分离点仍向下游急剧移动,马蹄涡的各项特征保持稳定;当ReDgt;8 000后,流动分离点向下游移动缓慢移动,马蹄涡的各项参数仍维持稳定。研究结果可为柱体基部科学布设防冲设施提供依据和参考。     

低柱体雷诺数下柱体上游薄层水流马蹄涡特征研究

马蹄涡是水下柱体基部产生局部侵蚀的主要动力,目前对大雷诺数深水条件的马蹄涡特征研究广泛,而在低柱体雷诺数浅薄层水流条件下因技术条件限制,实验手段难以准确捕获马蹄涡。为准确掌握在该水流条件下马蹄涡的特征,本研究构建了高分辨率高频的粒子图像测速系统,系统测量了6组浅薄层水流条件下柱体前端的瞬时流场。通过分析平均流场特征进而提取出流动分离点,采用旋转强度的方法来识别马蹄涡进而提取马蹄涡位置及强度,使用Oseen涡和纯剪切的控制方程叠加模拟马蹄涡进而计算出马蹄涡的半径。结果表明,在低柱体雷诺数水流条件下(ReDlt;5 000),随着柱体雷诺数的增加流动分离点急剧向下游移动,同时马蹄涡也急剧向柱体端和床面靠近,其半径减小而旋转强度增加。在浅薄层水流条件下,当柱体直径一定时,随着水深的增加,流动分离点向上游运动,同时马蹄涡向远离柱体端和朝水面运动,马蹄涡的半径增加,且其各项参数显著大于明渠水流条件下。随后,结合已有工作,归纳出分离点、马蹄涡特征随柱体雷诺数增加而经历的不同阶段:当5 000lt;ReDlt;8 000后流动分离点仍向下游急剧移动,马蹄涡的各项特征保持稳定;当ReDgt;8 000后,流动分离点向下游移动缓慢移动,马蹄涡的各项参数仍维持稳定。研究结果可为柱体基部科学布设防冲设施提供依据和参考。     

低柱体雷诺数下柱体上游薄层水流马蹄涡特征研究

马蹄涡是水下柱体基部产生局部侵蚀的主要动力,目前对大雷诺数深水条件的马蹄涡特征研究广泛,而在低柱体雷诺数浅薄层水流条件下因技术条件限制,实验手段难以准确捕获马蹄涡。为准确掌握在该水流条件下马蹄涡的特征,本研究构建了高分辨率高频的粒子图像测速系统,系统测量了6组浅薄层水流条件下柱体前端的瞬时流场。通过分析平均流场特征进而提取出流动分离点,采用旋转强度的方法来识别马蹄涡进而提取马蹄涡位置及强度,使用Oseen涡和纯剪切的控制方程叠加模拟马蹄涡进而计算出马蹄涡的半径。结果表明,在低柱体雷诺数水流条件下(ReDlt;5 000),随着柱体雷诺数的增加流动分离点急剧向下游移动,同时马蹄涡也急剧向柱体端和床面靠近,其半径减小而旋转强度增加。在浅薄层水流条件下,当柱体直径一定时,随着水深的增加,流动分离点向上游运动,同时马蹄涡向远离柱体端和朝水面运动,马蹄涡的半径增加,且其各项参数显著大于明渠水流条件下。随后,结合已有工作,归纳出分离点、马蹄涡特征随柱体雷诺数增加而经历的不同阶段:当5 000lt;ReDlt;8 000后流动分离点仍向下游急剧移动,马蹄涡的各项特征保持稳定;当ReDgt;8 000后,流动分离点向下游移动缓慢移动,马蹄涡的各项参数仍维持稳定。研究结果可为柱体基部科学布设防冲设施提供依据和参考。     

基于格子Boltzmann方法后台阶流动特性的研究

格子玻尔兹曼方法（LBM）是近几十年来发展的一种介观尺度的数值模拟方法,具有良好的处理复杂边界的能力和并行计算能力。利用LBM对后台阶流动进行了数值模拟,得到了不同雷诺数下后台阶流动的特性。研究结果表明,在低雷诺数情况下,台阶后再附着点与分离点的距离随着雷诺数的增加而变大;在高雷诺数情况下,台阶后再附着点与分离点的距离随着雷诺数的增加而增加,增长幅度逐渐降低,当Re=8 000时,再附着点与分离点之间的距离趋于稳定,并达到峰值;随着雷诺数的继续增大,这个距离开始减小,最终趋于一个定值。对于高雷诺数下的后台阶流动,采用LES⁃LBM的方法进行模拟,与前人的实验数据吻合得很好,表明LBM⁃LES模型对高雷诺数后台阶流动进行模拟具有较高的精确性。     

旋转圆柱尾涡及受力计算方法研究

为了研究旋转圆柱的三维流场和受力特性,本文采用沉浸边界-格子玻尔兹曼方法(IB-LBM),对均匀来流条件下旋转圆柱进行了数值模拟,分析了不同相对转速和雷诺数的影响。结果表明旋转作用使尾涡与流动方向之间产生一定夹角,相对转速的增大致使停滞点向上移动,并且影响柱体周围流速的损失。相对转速在1～2存在临界点,临界点前后速度波动范围发生改变。雷诺数能改变正负涡的脱落状态,增加柱体尾涡的三维性,并对升阻力产生相反的影响。计算结果验证了IB-LBM在处理旋转问题时的可行性,当前结果与文献中的结果相吻合。     

不同截面形状柱体流激振动能量转换特性分析

为了比较不同截面形状柱体的流激振动特性,探讨适用于流激振动能量转换的柱体截面形式,本文基于雷诺平均N-S方程,结合SST k-ω湍流模型和任意拉格朗日欧拉流固耦合动网格控制方法,对高雷诺数、高阻尼比条件下,圆柱、方柱、类梯形柱、梯形柱及T字形柱等5种不同截面形状的柱体流激振动进行数值模拟,模拟中保持各柱体的特征长度相同,质量比、阻尼比、同一流速下的雷诺数和约化速度等重要无量纲参数也均保持一致。计算结果表明:圆柱的振动表现为典型的涡激振动,其他4种柱体的振动均为驰振,在较低约化流速下T字形柱获得的功率均较其他柱体获得的为大。各柱体能量转换效率均随约化速度的增加先增大后减小,中间存在一个峰值,5种柱体中T字形柱的最大能量转换效率为最大,达到了42. 5%,其次为圆柱,为27%。为了达到最大能量转换效率,圆柱和方柱所需要的流速最小,其次为T字形柱。从分析结果来看,T字形柱和圆柱对于低流速下的能量转换较为有利。     

错列双圆柱的脉动升力特性及其流场机理

受到上游圆柱尾流作用,下游圆柱在小风向角下易发生尾流激振现象,脉动升力与尾流激振联系紧密,但以往针对脉动升力特性的研究较少,其流场作用机理尚未澄清.针对间距比为P/D=1.5~4的错列双圆柱,采用大涡模拟方法,在高亚临界雷诺数下(Re=1.4×10⁵),重点研究小风向角下（β=0°~30°）圆柱的脉动升力和流场结构之间的内在关系,基于流场特征讨论了上、下游圆柱之间的干扰机理.研究结果表明:在小风向角范围内,上、下游圆柱的脉动升力和流场特性变化剧烈,双圆柱会经历尾流干扰、剪切层干扰、邻近干扰、旋涡撞击及涡的相互作用五种干扰流态;在尾流干扰以及剪切层干扰流态下,下游圆柱的脉动升力远小于单圆柱;在邻近干扰流态下,上、下游圆柱相互间干扰较弱,下游圆柱的脉动升力接近于单圆柱;而在旋涡撞击流态和涡的相互作用流态下,下游圆柱的脉动升力则大于单圆柱,上游圆柱的尾流旋涡对下游圆柱的撞击或者与下游圆柱旋涡的相互作用会导致下游圆柱的脉动升力急剧增加.     

横向磁场作用下Taylor-Couette流动数值模拟

采用Fluent模拟软件对横向磁场作用下液态金属镓的Taylor-Couette流动进行了数值模拟。考虑了有限长和无限长两种模型,模型半径比为1/2,长宽比为4。哈特曼数选取为0～50,雷诺数选取为10和259,分别对应层流Couette流动和Taylor-Couette流动两种工况。结果表明,对于雷诺数等于10工况,流体在横向磁场的作用下产生非轴对称的流型;随磁场强度增加,涡由圆环外侧移动到内侧,并呈现出扁平化趋势。对于雷诺数等于259工况,在磁场作用下,流场内泰勒涡尺寸变小,泰勒涡数量由2对增加至3对,并最终破裂。对于有限长和无限长模型,由于上下壁面边界条件不同,流型变化呈现出不同的特点。     

新型涡流发生器几何结构及换热计算关联式

通过对矩形风道内置有新型涡流发生器一斜截椭圆柱体强化传热性能和压降特性的实验研究,分析了不同雷诺数(Re=800~38000)工况下涡流发生器的几何尺寸和迎流方位对传热和压降的影响.实验结果表明:高宽比h/b=l/2的斜截椭圆柱体较好的几何尺寸和位置参数是斜边倾角α=20°,迎流攻角β=45°.这是马蹄涡和端部涡综合作用的结果.根据77个实验工况回归了斜截椭圆柱体涡流发生器强化换热的准则实验关联式,为工程应用提供了理论依据.实验采用了稳态的气水逆流换热方式,测温数据由HP34970A数据采集系统采集.     

低雷诺数下串列双圆柱涡激振动的数值模拟

针对海洋立管中常发生的流致振动问题,本文采用自主研发的CIP-ZJU数值模型,对雷诺数Re=150条件下串列双圆柱的涡激振动进行模拟。该模型在笛卡尔网格系统下建立,采用具有三阶精度的CIP方法求解N-S(navier-stokes)方程,采用浸入边界法处理流-固耦合问题。本文仅考虑圆柱的横向振动,具体分析不同间距比和折合速度,并分别考虑上游圆柱固定和自由振动两种工况,得到柱体振动响应、受力响应和流场信息,验证了本模型在处理柱体涡激振动问题的有效性。结果表明:串列情况下下游圆柱的最大振幅要明显大于单柱的情况,双圆柱涡激振动的阻力系数普遍比单圆柱涡激振动时要小;上游圆柱固定时,下游圆柱的振动频率几乎不由斯特劳哈尔频率控制;而当上游圆柱自由振动时,在折合速度4≤Ur≤5时,上游圆柱后方产生两列涡,使下游圆柱的运动范围限制在两列涡之间,可对其振动产生了抑制作用。     

串列双圆柱绕流的模态特征及流场重构

为探讨不同间距下,串列双圆柱绕流场特征的内在变化规律,在低雷诺数Re=100下,针对间距比为P/D=1.1~5的串列双圆柱,通过动力学模态分解（dynamic mode decomposition,DMD）方法对其绕流场进行模态分解,并基于DMD的主导模态建立降阶模型,重构了串列双圆柱的涡量场。结果表明:串列双圆柱的3种流态,即单一钝体（P/D=1.1~2）、剪切层再附（P/D=3）和双涡脱流态（P/D=4~5）,呈现明显不同的子模态特征;随着间距比的增大,与圆柱涡脱频率对应的模态结构会逐渐由下游圆柱尾流转移至上游圆柱尾流,此主模态表征了3种流态随间距比演变的内在机制;与单一钝体及剪切层再附流态相比,双涡脱流态的下游圆柱近尾流区的高阶模态结构更为复杂,流场重构时需要更多模态才能达到与其他两种流态相似的精度,在尾流旋涡影响区的重构误差最大。     

基于气体动理学格式的平板绕流数值模拟

介绍了基于Bhatnagar-Gross-Krook（BGK）模型的气体动理学格式Gas Kinetic Scheme（GKS）数值算法。给出了平板绕流的物理模型以及边界的处理方式,同时进行了算法准确性验证。对不同雷诺数条件下的单平板绕流进行数值模拟,随雷诺数的增大,涡不断被拉长,逐渐失去对称性,形成卡门涡街;对Re=300的前后平板绕流,当板间距较小时流场处于对称尾流区,没有涡的脱落,当板间距较大时会形成卡门涡街;对Re=300的上下平板绕流,当板间距较小时流场特征与单平板绕流相似,当板间距较大时会形成两列涡,在下游较远处开始逐渐合并。     

湍流模型下双柱绕流升阻力系数数值模拟

本文使用流体软件FLUENT在RNG κ-ε湍流模型下.取雷诺数Re为20 000,间距比s/d等于1.75,模拟了双柱在串列、并列及两柱中心的连线与水流方向成30°和60°夹角错置这四种工况下的绕流升阻力系数.两圆柱阻力系数随夹角增大而增大,以并列双稳态对应的阻力系数最大.双柱串列时,下游圆柱阻力系数为负,两柱升力合力为0.双柱错置和并列时,两圆柱阻力系数都为正值,每个圆柱的升力合力都不为0,以并列发生双稳态流动时的升力合力最大.     

柱体绕流的CIP方法模拟

柱体绕流和流致振动现象是一个复杂的工程问题,利用自主研发的CIP-ZJU模型,对低雷诺数(Re300)圆柱和方柱绕流问题开展了数值模拟。模型在直角坐标系统下建立,采用紧致插值曲线CIP方法作为流场的基本求解器离散了Navier-Stokes方程,基于多相流的理论实现流-固耦合同步求解,利用浸入边界方法处理固体边界。模拟结果与文献结果进行比较,二者吻合情况较好。通过引入压力阻力项和摩擦阻力项,分析了Strouhal数、阻力系数、升力系数等参数随雷诺数的变化情况。结果表明,圆柱绕流与方柱绕流在水动力参数的变化规律上存在多处差异。     

流体横掠圆柱体传热特性数值模拟的流态模型选择

采用层流模型、标准k-ε紊流模型、RNG k-ε紊流模型,对亚临界雷诺数条件下(Re=300~300 000)的二维圆柱扰流进行了数值模拟与仿真,得到了圆柱壁面平均努塞尔数随雷诺数变化的规律。并得出结论:当Re≤3 900时,采用层流模型较合适;当3 900≤Re≤22 000范围内时,采用标准k-ε紊流模型较合适;当22 000≤Re≤300 000时,采用RNG k-ε紊流模型较合适。     

两种管径比的绕流热柱体旋涡信号比较

在考虑壁面热效应的情况下,当雷诺数为0.6×104~5.0×104,分别对管径比为0.4和0.24两种工况下的柱体尾迹旋涡脱落的流动情况进行实验.结果表明:加热使得旋涡流动发生明显变化,其影响表现在旋涡脱落信号傅里叶变换峰值的改变上;这种影响随雷诺数的变化有一临界值ReL,加热在ReL两侧对旋涡所起作用相反,并且ReL随管径比的不同而不同.通过对不同管径比工况的分析比较,发现同等加热强度下小管径比的壁面热效应对旋涡流动影响更明显.     

流道内分隔板长度对方柱绕流特性的影响

为了揭示流道内方柱-分隔板钝体的绕流特性及分隔板长度对尾迹脱涡的影响,对低雷诺数(Re=150)下阻塞比为1/6的流道内的方柱-分隔板结构进行数值研 究,讨论方柱和分隔板的受力及自由剪切层的变化趋势.结果表明：在有限流道内,分隔板对方柱的尾流结构产生了明显的影响,抑制了尾流区剪切层之间的相互作 用,并在一定范围内使方柱阻力减小;随着分隔板长度L的变化(0D～6D,D为方柱边长),流动特性呈现出3个不同阶段,即当0≤L/D≤1.00时,随着板长度的增加,自 由剪切层沿顺流方向逐渐下移;当1.25≤L/D≤4.50时,自由剪切层在板尾端卷起并在分隔板尾迹边缘清晰地出现二次涡;当L/D≥4.75时,分隔板完全阻断剪切层 之间的相互作用,使自由剪切层再次附着在分隔板上.     

明渠湍流横向涡旋的尺度与环量特征

为了分析明渠湍流中横向涡旋的尺度与动力学特征,利用自主研发的高频粒子图像测速系统测量了6种水流条件下的明渠恒定均匀流时间序列流场,使用涡旋提取方法识别涡旋,提出一套简洁客观的涡旋尺度及环量确定方法,分析了明渠湍流横向涡半径与环量的特征。涡旋的无量纲尺度与强度仅与所在位置有关,与雷诺数等其它因素均无关系。从壁面附近至对数区的上边界,涡旋的平均半径增大。受水面影响,半水深以上的涡旋半径有逐渐减小的趋势。在粘性作用下,涡旋的平均环量由壁面至水面总体呈减小趋势。时均剪切对逆向涡旋的尺度与强度均有抑制作用,其平均尺度与环量在全水深均小于正向涡。     

排桩式波浪能发电装置附近流场特性研究

基于雷诺时均Navier-Stokes方程和k?ω湍流模型,研究单排桩式振荡水柱式波浪能发电装置在规则波作用下的桩基附近流场特性. 通过物理模型实验验证所建数值波浪水槽的准确性. 模拟规则波作用下装置附近的流场特性,分析装置附近的涡特征以及其对装置附近泥沙冲刷情况的潜在影响. 结果表明,测试工况下,在桩式振荡水柱装置桩柱体的桩基附近观察到马蹄涡和尾涡现象,马蹄涡强度随着Keulegan–Carpenter（KC）数的增大而减小,尾涡强度随着KC数的增大而增大. 马蹄涡因强度过小,在模拟的KC数范围内不是装置桩基泥沙冲刷的主要因素,高强度的尾涡很可能成为装置桩基附近泥沙起动、输运和冲刷的重要影响因素.     

完全湍流剪切层对圆柱涡激振动特性的影响

为研究完全湍流剪切层(2×10~4~4×10~4Re1×10~5~2×10~5)下雷诺数对涡激振动动态响应的影响,在OpenFOAM平台下采用有限体积法求解雷诺平均Navier-Stokes方程和Spalart-Allmaras一方程湍流模型,采用二阶范德波尔方程描述二维圆柱体涡激振动,对完全湍流剪切层下的圆柱绕流涡激振动进行数值计算,分析涡激振动的振幅、频率特性,并在Williamson-Roshko(W-R)图谱中分析旋涡尾迹脱落模式.结果表明:在相同折减速度下提升雷诺数时,圆柱振幅会显著增大,最大振动振幅值高于自由剪切层过渡旋涡区结果及修正Griffin图谱结果.雷诺数对提高圆柱的涡激振动振幅和同步区的范围具有决定性作用,在完全湍流剪切层中,当提升阻尼参数(α≥0.38)时,涡激振动的上部分支消失,振动频率出现"阶跃".