基于浸入边界-格子Boltzmann通量求解法的椭圆柱流动特性分析

基于浸入边界-格子Boltzmann通量求解法,开展了雷诺数Re=100不同几何参数下单椭圆柱及串列双椭圆柱绕流流场与受力特性对比研究。结果表明,随长短轴比值的增加,单椭圆柱绕流阻力系数先减小后缓慢上升,最大升力系数则随长短轴比值的增大而减小;尾迹流动状态从周期性脱落涡到稳定对称涡。间距是影响串列圆柱及椭圆柱流场流动状态的主要因素,间距较小时,串列圆柱绕流呈周期性脱落涡状态,而椭圆柱则为稳定流动;随着间距增加,上下游圆柱及椭圆柱尾迹均出现卡门涡街现象,且串列椭圆柱临界间距大于串列圆柱。串列椭圆柱阻力的变化规律与圆柱的基本相同,上游平均阻力大于下游阻力;上游椭圆柱阻力随着间距的变大先减小,下游随间距的变大而增加,当间距达到临界间距时上下游阻力跃升,随后出现小幅度波动再逐渐增加,并趋近于相同长短轴比值下单柱体绕流的阻力。     

三维不等直径串列圆柱绕流的双涡脱落流态频率特征研究

采用大涡模拟方法计算Re=2×10~3三维不等直径串列圆柱(d/D≤1)绕流问题。结果显示,处于双涡脱落流态时,随着串列圆柱间距增加,上游圆柱量纲为一的涡脱频率值St₁总体上升,而下游圆柱量纲为一的涡脱频率值St₂存在先下降后上升的变化规律。在圆柱间距较小的情况下,St₂随着串列圆柱间距的增加而减小,量纲为一的涡脱频率比值、直径比与间距比之间近似满足St₂/St₁∝(L/D)^-1/4(d/D)的幂指数关系;在圆柱间距较大的情况下,圆柱间时均流向速度提高并趋近主流区速度,St₂随间距比增加而上升。在较小直径比串列圆柱情形下,下游圆柱量纲为一的涡脱频率St₂可下降至更低的临界拐点,从而产生“次谐波涡脱锁定”现象。     

串列双方柱体流体动力载荷研究

本文通过流动显示,热线测频和流体动载荷测量在水槽中研究了绕经不同柱间距比S/D(S为双柱间距,D为柱体截面宽)串列双方柱体流动特性。实验雷诺数为Re=6×10~3,柱间距比0.5≤S/D≤10实验测量了涡脱落频率、时间平均阻力、动态阻力和动态升力。通过实验结果综合分析给出临界柱间距范围2.5≤(S/D)_(cr)≤3.0,并将串列双方柱流动随柱间距的变化划分为二种流态区。在临界柱间距,作用于双柱体的流体载荷、涡脱落频率以及流谱都发生跃变。文中分析讨论了两个流态区的特性以及在临界柱间距出现的双稳态特性。     

后柱横向振动串列双圆柱尾迹特性

本文采用电化学显示方法和热线技术就雷诺数为10~2—10~3,后柱横向振动串列双圓柱尾迹和涡街进行了研究。双柱间距为1.5≤s/d≤10,后柱振动振幅为0.1d和0.3d。首先,本文得到了低雷诺数情况下,串列双圆柱绕流旋涡脱落频率随柱间距和雷诺数的变化规律。并且,通过流谱分析,对于出现频率间断现象给予了说明。其次,本文着重研究了振动情况下尾迹特性。当后柱以静止双柱涡脱落频率强迫振动,涡脱落和柱振动之间将发生联锁反应。本文实验揭示了联锁区奇特的涡街流谱——二次涡、涡的聚集和涡街伸长。     

高雷诺数下多柱绕流特性研究

采用改进的延迟分离涡方法数值模拟了高雷诺数下的柱体绕流,包括单圆柱绕流、单方柱绕流、串列双圆柱绕流和串列双方柱绕流,研究了不同雷诺数下圆柱绕流与方柱绕流的水动力特性.计算结果与实验数据及其他文献的数值计算结果吻合良好,研究表明,单方柱绕流在2.0×10~3Re1.0×10~7范围内未出现类似于单圆柱绕流的阻力危机现象,其平均阻力系数C_d、升力系数均方根C'_1及斯特劳哈尔数S t维持在一定范围内波动.串列双圆柱绕流与串列双方柱绕流中,均选取L/D=2.0,2.5,3.0,3.5和4.0这五中间距比进行计算.串列双圆柱绕流中,当Re=2.2×10~4时,在3.0L/D3.5内存在一临界间距比(L_c/D)使得L_c/D前后上下游圆柱的升阻力系数发生跳跃性变化,且当L/DL_c/D时,下游圆柱的阻力系数为负数.而当Re=3.0×10~6时,则不存在临界间距比,且下游圆柱的阻力系数始终为正数.串列双方柱绕流在Re=1.6×10~4和Re=1.0×10~6两种工况下的临界间距比分别处于3.0L/D3.5和3.5L/D4.0区间内,且当L/DL_c/D时,两个雷诺数下的下游方柱阻力系数均为负数.     

分块法研究圆柱绕流升阻力

使用新的分块耦合方法，分别对单圆柱和串列双圆柱绕流进行了数值模拟．对于单圆柱绕流，低Re下计算所得到的定常涡尺寸与实验非常接近．对于串列双圆柱绕流，研究分析了改变双圆柱中心间距对上下游圆柱的升阻力系数和脉动频率所产生的影响，计算结果与实验非常吻合，为进一步研究涡致振动提供了依据．     

分块法研究圆柱绕流升阻力

 使用新的分块耦合方法，分别对单圆柱和串列双圆柱绕流进行了数值 模拟. 对于单圆柱绕流，低$Re$下计算所得到的定常涡尺寸与实验非常接近. 对于 串列双圆柱绕流，研究分析了改变双圆柱中心间距对上下游圆柱的升阻力系数和脉动频率所 产生的影响，计算结果与实验非常吻合，为进一步研究涡致振动提供了依据.     

虚拟边界法研究正交双圆柱及串列双圆球绕流

把Goldstein等人提出的虚拟边界法推广到三维情况,研究了 Re=150时不同间距下正交双圆柱绕流,和Re=250时不同间距下串列双 圆球绕流流场. 对于正交双圆柱绕流,当间距比大于3,下游圆柱对上游圆柱尾流的影响只 限定在下游圆柱的尾流所扫过的范围之内;当间距比小于等于3,下游圆柱对上游圆柱尾流 的影响扩大,下游圆柱尾流扫过区上下出现两排三维流向二次涡结构. 对于串列圆球绕流, 研究发现,在小间距比(L/D≈ 1.5)的情况下,由于上下游圆球尾流区的相互抑 制消除了压力不稳定性,整个流场呈现稳 态轴对称特征;间距比为2.0时,周向压力梯度诱发出流体的周向输运,流场呈现稳态非对 称性,但流场中存在特定的对称面;间距比增大到2.5后,绕流场开始周期振荡,原有的对 称面依旧存在;在间距比3.5时下游圆球下表面的涡结构强度有所减弱,导致占优频率发生 交替;间距比增至7.0时,整个流场恢复稳态特征,两圆球尾部同时出现双线涡,这时流场 对称面的位置发生了变动.     

串列双锥柱绕流流动特性的大涡模拟研究

利用大涡模拟研究了雷诺数Re = 3900下串列双锥柱在间距比L/D_m = 2 ~ 10下的升阻力特性及三维流动结构. 研究发现: 上游锥柱在后方形成的两个展向不对称回流区, 使其后方压力分布不对称. 上游锥柱发展的上洗、下洗和侧面剪切层作用在下游锥柱的附着点位置不同是上游和下游锥柱时均阻力系数和脉动升力系数变化的主要原因, 串列双锥柱间流动结构随间距比变化可分为三种状态: 剪切层包裹状态, 过渡状态及尾流撞击状态. 剪切层包裹状态. 上游锥柱的自由端主导来流在下游锥柱迎风面影响范围广, 上游锥柱剪切层完全包裹住下游锥柱, 从而抑制下游锥柱后方回流区形成, 导致下游锥柱时均阻力系数降低; 尾流撞击状态; 上游锥柱尾流得到充分发展, 其回流区大小随间距比增大不再发生变化, 上游锥柱尾流出现周期性脱落, 撞击在下游锥柱表面, 从而使脉动升力系数大幅增加, 最大脉动升力系数较单直圆柱提升约20.7倍; 过渡状态, 此时时均阻力系数和脉动升力系数均会较剪切层包裹状态增加. 该研究可以为风力俘能结构群列阵布局提供理论支持.      

串列双圆柱体绕流特性与互扰效应研究

基于大涡模拟(LES)方法对亚临界雷诺数(Re=3900)下三维串列双圆柱体绕流问题进行了数值计算。首先,通过求解单圆柱算例来验证计算模型及参数的正确性。然后,着重分析了不同间距比对双圆柱体的流体力系数的影响,并阐述了双圆柱体流场特性变化及其互扰效应内在机理。研究表明:雷诺数Re=3900时,串列双圆柱体绕流临界间距比在3.9～4.0之间;随着间距比的增加,双圆柱体临近流场中二次涡团形成的区域与三维涡结构均会发生变化,导致其结构表面所受的流体力系数在时间与空间上变化的规律性逐渐减弱;达到临界间距比时,流体力系数的变化会呈现出较强的规律性。     

串列双圆柱同步测压实验结果的相关分析

本文在雷诺数2×104下,同步测量了12个不同间距下串列双圆柱的表面压力分布,积分得到脉动升、阻力的时间历程,并对前、后柱之间的脉动升、阻力以及脉动升阻力和圆柱表面的脉动压力进行了相关分析.在本实验中,串列双圆柱流态切换的临界间距比在3.5～4.0之间,在临界间距前后,相关曲线的形态差别很明显.当串列双圆柱之间的距离小于临界间距时,前后柱之间的影响显著,造成脉动升/阻力和脉动压力之间的相关程度下降.串列双圆柱之间的距离大于临界距离时,前、后柱的脉动升力之间的相位差随间距线性增加.分离点的脉动压力和脉动升力之间以及背压和脉动阻力之间都有较好的相关性,可以用分离点的脉动压力和背压来分别表征脉动升力和脉动阻力.     

并列双椭圆柱绕流的格子Boltzmann-虚拟区域方法的模拟研究

采用格子Boltzmann-虚拟区域方法对并列双椭圆柱绕流进行了模拟研究。首先,通过与并列双圆柱的结果进行对比,验证了数值方法的有效性。其次对雷诺数为200时两种间距(g=0.5a和2.0a,g为柱体表面间距,a为椭圆长轴)的情况进行了研究,考察了椭圆长轴与短轴之比,即α值对柱体升阻力系数及涡结构的影响。研究发现,与圆柱相比,对于g=0.5a椭圆柱的升阻力系数可能出现两种变化,一是升阻力随时间演化较规则,接近周期性;二是流场可能长时间偏向于其中一个椭圆柱,这些变化与α的值有关。对于g=2.0a,两个椭圆柱后的某一区域内会出现四列涡街,经过一段时间,四列涡街又会演化成两列向两侧扩张的涡街。     

相同雷诺数下串列双圆柱绕流的仿真分析

为研究均匀水流场中串列排布的柱群之间的干涉影响,本文以三维串列双圆柱为例,通过计算流体力学(CFD)软件FLUENT15.0中双方程k-ε模型,分析模拟了双圆柱所受平均阻力、平均升力、后柱周向压力、斯特劳哈尔数等水动力特性。结果表明:在雷诺数为Re=2×10~4的串列双圆柱绕流中,两圆柱中心间距L与圆柱直径D的比值为L/D=4时,后柱受前柱绕流尾流影响大,明显高于单圆柱绕流的平均阻力系数,后柱的周向压力值也随前柱尾流的摆动呈现显著的不对称性;当L/D=8时,前柱绕流尾流对后柱影响逐渐减弱;当L/D=12时,两圆柱之间的相互干扰几乎可以忽略,可以看作是相互独立的单圆柱绕流。最后,计算的斯特劳哈尔数与单圆柱绕流对应的斯特劳哈尔数相近且仿真数值在计算数值范围之内,验证了整个仿真分析的准确性,也进一步说明了双圆柱绕流的柱群的干涉影响。双圆柱间距越大,前、后柱之间的干涉影响越弱。     

流向磁场作用下圆柱绕流的直接数值模拟

绕流是托卡马克装置中液态包层内常见的流动形态,对流场与热量分布有着重要的影响.本文通过直接数值模拟(DNS),研究了不同磁场强度下$Re=3900$的圆柱绕流,分析了磁场强度对于湍流尾迹的影响.无磁场情况下,直接数值模拟的结果与前人的实验及模拟结果吻合很好.圆柱下游的尾迹中,随着流向距离的增大, 流向速度剖面逐渐从U型进化呈V型, 并慢慢趋于平缓,这表明尾迹中的流动结构受圆柱影响逐渐减小.圆柱后方两侧的剪切层中,由于Kelvin-Helmholtz不稳定性的影响,可以清晰地看到小尺度剪切层涡的脱落.通过对无磁场的计算结果施加流向磁场,本文计算了哈特曼数($Ha$)分别为20, 40和80的工况,以研究磁场效应对于湍流的影响.结果表明磁场较弱时,流动依然呈三维湍流状态.随着磁场增强, 近圆柱尾流区受磁场抑制明显,回流区被拉长,剪切层失稳位置向下游转移.圆柱后方的涡结构由于受到竖直方向洛伦兹力的挤压作用,随着哈特曼数的增加尾迹区域逐渐变窄.相比于无磁场情况的涡结构,由于磁场的耗散作用,相应的涡结构尺度变小.该研究不仅扩展了现有磁场下湍流运动的参数范围,对于液态包层的设计及安全运行同样具有重要的理论指导意义和工程应用价值.      

DIRECT NUMERICAL SIMULATIONS ON THE TURBULENT FLOW PAST A CONFINED CIRCULAR CYLINDER WITH THE INFLUENCE OF THE STREAMWISE MAGNETIC FIELDS1)

绕流是托卡马克装置中液态包层内常见的流动形态,对流场与热量分布有着重要的影响.本文通过直接数值模拟(DNS),研究了不同磁场强度下$Re=3900$的圆柱绕流,分析了磁场强度对于湍流尾迹的影响.无磁场情况下,直接数值模拟的结果与前人的实验及模拟结果吻合很好.圆柱下游的尾迹中,随着流向距离的增大, 流向速度剖面逐渐从U型进化呈V型, 并慢慢趋于平缓,这表明尾迹中的流动结构受圆柱影响逐渐减小.圆柱后方两侧的剪切层中,由于Kelvin-Helmholtz不稳定性的影响,可以清晰地看到小尺度剪切层涡的脱落.通过对无磁场的计算结果施加流向磁场,本文计算了哈特曼数($Ha$)分别为20, 40和80的工况,以研究磁场效应对于湍流的影响.结果表明磁场较弱时,流动依然呈三维湍流状态.随着磁场增强, 近圆柱尾流区受磁场抑制明显,回流区被拉长,剪切层失稳位置向下游转移.圆柱后方的涡结构由于受到竖直方向洛伦兹力的挤压作用,随着哈特曼数的增加尾迹区域逐渐变窄.相比于无磁场情况的涡结构,由于磁场的耗散作用,相应的涡结构尺度变小.该研究不仅扩展了现有磁场下湍流运动的参数范围,对于液态包层的设计及安全运行同样具有重要的理论指导意义和工程应用价值.     

基于格子Boltzmann方法的低雷诺数四圆柱绕流流动模式与受力特性模拟

利用格子Boltzmann方法模拟了雷诺数为100时,均匀来流条件下的二维菱形排布的四柱绕流现象,得到了不同柱间距比下的绕流流动模式及阻力变化规律。结果表明:圆柱互扰效应与柱间距比有关,当L/D≤1.2时为单钝体模式,圆柱互扰效应以临近效应为主;当1.2相似文献   

不同雷诺数下二维椭圆柱绕流的数值模拟研究

利用计算流体力学软件Fluent对不同雷诺数(Re=100,3900,3.5×10~6)下二维椭圆柱绕流进行了数值模拟研究,分析了不同轴长比(2b/2a=cosθ,θ=0°,15°,30°,45°,60°)下椭圆柱绕流的特性。通过对比尾部涡流情况、升力系数C_L、阻力系数C_D以及斯特劳哈尔数St初步发现:由于椭圆形截面偏流线型的特点,在三种雷诺数下随着θ的增大椭圆柱绕流尾涡强度减小,流场的变化使圆柱表面的压力系数减小,最终导致圆柱的升力系数幅值与阻力系数均值减小。而斯特劳哈尔数St在三种雷诺数下的变化不同,随着θ的增大,层流雷诺数(Re=100)下St值减小;亚临界雷诺数(Re=3900)下St值在45°处轻微上扬,在60°处明显减小;超临界雷诺数(Re=3.5×10~6)下St值增大。     

湍流涡致振动频率特性

用数值模拟方法对固定圆柱湍流涡脱落频率与弹性圆柱湍流涡致振动频率特性进行了研究,湍流计算模型采用标准κ-ε模型,压力泊松方程提法基于非交错网格系统.研究结果表明:固定圆柱湍流绕流涡脱落频率基本不随雷诺数而变,对于同一固有频率弹性圆柱,涡振频率基本不随雷诺数而变;对于某一固定雷诺数流动涡振频率在一定范围内与系统固有频率有关.     

正方形排列四圆柱体绕流特性及互扰效应研究

基于计算流体动力学理论,运用大涡模拟方法对雷诺数Re=3900三维正方形排列四圆柱体结构群的绕流问题进行数值计算,主要分析来流攻角与间距比两个参数对四圆柱体结构群流体参数及流场模态的影响。结果表明:来流攻角与间距比均对四圆柱体结构群绕流特性有较强的影响;来流攻角θ=0°、22.5°、45°下,临界间距比分别为3.5、4.0、3.0;间距比的变化会导致下游圆柱表面压力系数分布发生改变;另一方面,间距比较小时,四圆柱体结构之间的互扰作用均以临近效应为主;随间距比增大,上游圆柱尾流对下游圆柱有显著影响,其互扰作用会转变尾激效应。     

串列双圆柱绕流问题的数值模拟

本文运用有限体积方法,对绕串列放置的双圆柱的二维不可压缩流动进行了数值计算。为研究两圆柱不同间距对圆柱相互作用和尾流特征的影响,选取间距比Ｌ／Ｄ（Ｌ为两圆柱中心间的距离,Ｄ为圆柱直径）在１．５～５．０之间每隔０．５共八个有代表性的间距进行了计算模拟。计算均在Ｒｅ＝２００条件下进行。计算结果表明：对该绕流问题,流动特征在很大程度上取决于间距的大小。且间距存在一临界值,间距比从小于临界值变化到大于临界