双排并列三方柱绕流数值模拟

采用有限体积法对雷诺数Re=2.2×104条件下双排并列三方柱绕流问题进行数值模拟,对横向桩距比L/D=1.5、2、2.5、3,纵向桩距比S/D=2、2.5、3共12种工况下绕流流场分别进行计算,得到了每种工况下各方柱绕流阻力系数值及其随桩间距的变化规律,同时分析了绕流流场的演变及方柱前后压差随桩间距的变化趋势,计算结果表明:多方柱绕流较单方柱已明显表现出群桩阻水效应,各方桩绕流参数变化幅度较大,流场演变相对复杂。     

并列双方柱在偏向流流态下的气动性能

小间距并列双方柱存在独特的偏向流现象，为了进一步揭示并列双方柱在偏向流流态下的气动性能，在雷诺数（Re）为8.0×104、间距比（T/B，其中T为方柱中心间距，B为方柱边长）范围为1.25 ~ 5条件下，通过测压风洞试验方法研究了并列双方柱的表面风压、气动力、Strouhal数等气动参数随间距比的变化规律，重点探讨了发生偏向流现象时气动力的时频特性和气动力展向相关性. 并列双方柱的偏向流现象主要发生在T/B = 1.25~2的范围内，气流的偏转方向具有随机性，这一现象会随间距比的增大而逐渐消失；发生偏向流现象时，两个方柱的平均阻力系数、脉动风压系数均小于单方柱，但两个方柱间隙侧壁面受到的负压远强于单方柱；与宽尾流方柱相比，窄尾流方柱的平均阻力系数更大，脉动气动力系数更大，间隙侧壁面受到平均负压更强，气动力系数沿方柱展向的相关性更强.     

全风向角下双方柱脉动气动性能试验研究

为阐明双方柱中上游方柱尾流对下游方柱脉动气动性能的影响,在雷诺数Re=8.0×104均匀来流条件下,完成了间距比P/B=1.25～5(其中P为柱心距,B为方柱边长)、全风向角α=0°～90°的双方柱测压风洞试验.分析了不同P/B条件时下游方柱脉动气动力、升力功率谱、Strouhal数和气动力展向相关性随α的变化规律.发现P/B<3时下游方柱的脉动气动力整体上小于单方柱,且α≥40°时几乎不随风向角变化.但在3≤P/B≤5且0°≤α≤30°时下游方柱脉动阻力远大于单方柱.下游方柱的脉动气动力、升力功率谱和气动力展向相关性随α的变化规律在P/B=3前后差异显著.P/B<3时,下游方柱旋涡脱落受到显著抑制;P/B>3时,下游方柱的脉动气动性能近似于单方柱.     

数值分析串列双圆柱绕流

为研究不可压缩流动中二维的串列双圆柱中不同间距对圆柱间相互作用和尾流特征的影响,选取间距比L/D(L为两圆柱中心间的距离,D为圆柱直径)为2、3、4共3个间距进行了数值分析.计算均在Re=200的非定常条件下进行.计算了圆柱的升阻力系数、尾涡脱落频率等描述绕流问题的主要参量.并与参考文献的数值结果比较,验证了用FLUENT分析的有效性.     

串列双方柱的风压特性及其流场机理

为了研究干扰条件下方柱的风压特性及其流场机理,以串列双方柱为研究对象,采用大涡模拟方法,在雷诺数Re=8×104、间距比P/B=1.1～5的条件下,研究了两个方柱的风压系数、气动力系数、风压非高斯特性、风压相关性随间距比的变化规律,重点探讨了双方柱流场特性及其与风压非高斯特性的内在联系.研究结果表明:随着间距比的增大,串列双方柱依次表现为3种流态,即单一钝体、剪切层再附和双涡脱流态,风压特性与流态密切相关.风压的非高斯特性和风压相关性随流态变化呈现为3种类型:在单一钝体流态下,柱间回流区附近的表面风压呈现明显的非高斯特性且风压相关性较强;在剪切层再附流态下,方柱尾流的涡脱强度低,风压相关性弱,但风压非高斯区域大;在双涡脱流态下,受上游方柱尾流旋涡的作用,方柱侧风面的风压相关性较强,下游方柱的侧风面和背风面出现大范围的风压非高斯区域.     

串列双方柱气动性能的流场机理研究

采用非定常雷诺平均法和SST k-ω湍流模型,在雷诺数Re=2.2×104时对五种间距串列双方柱的绕流场进行了数值模拟。结果表明：随着间距比的增大,串列方柱的气动性能会发生剧烈变化,其绕流场会经历三种不同流态：间距比为S/B=1.2时,下游方柱完全被上游方柱的分离剪切层包裹,流场呈现单一钝体流态;S/B=1.5和2时,在上游方柱上分离的剪切层会再附在下游方柱侧面,流场呈现剪切层再附流态,并在两个方柱之间形成强烈的回流区;S/B=3和4时,两个方柱的尾流中都会形成规则的涡街,流场呈现双涡脱流态,此时上游方柱的旋涡会与下游方柱发生复杂的相互作用,造成下游方柱受到很大的脉动风压作用。     

串列双方柱气动干扰机理的数模模拟研究

采用非定常雷诺平均法和SST k-ω湍流模型,在雷诺数Re=2.2×104时对五种间距串列双方柱的绕流场进行了数值模拟。结果表明:随着间距比的增大,串列方柱的气动性能会发生剧烈变化;其绕流场会经历三种不同流态:间距比为S/B=1.2时,下游方柱完全被上游方柱的分离剪切层包裹,流场呈现单一钝体流态;S/B=1.5,2时,在上游方柱上分离的剪切层会再附在下游方柱侧面,流场呈现剪切层再附流态;并在两个方柱之间形成强烈的回流区;S/B=3,4时,两个方柱的尾流中都会形成规则的涡街,流场呈现双涡脱流态。此时上游方柱的旋涡会与下游方柱发生复杂的相互作用,造成下游方柱受到很大的脉动风压作用。     

全风向角下小间距双方柱绕流的大涡模拟

采用大涡模拟方法研究了间距比（P/B）为1.5(P为柱心间距,B为方柱边长)、风向角α为0°～90°等条件下双方柱在均匀来流作用下的气动力、流态划分、表面风压和流场特性。研究发现：小间距双方柱流动干扰效应显著,下游柱平均气动力随风向角的变化规律与上游柱和单方柱差异较大,且可能受更大的升力绝对值。将小间距比双方柱绕流分为前角分离流态（α=0°～10°）、分离泡流态（α=20°～30°）、附着流流态（α=40°～60°）及间隙侧分离泡流态（α=70°～90°）四种模式。附着流流态和间隙侧分离泡流态的间隙区出现较强的负压。尾流负压区的强度随风向角先增强后减弱,在α=50°附近达到峰值。前角分离流态时双柱具有一个整体的尾流回流区,而在其他流态下上下游方柱均有独立的回流区。     

顺排管束流动和传热数值模拟

采用数值模拟的方法考察了列管式换热器中顺排管束的流动和传热问题,计算了不同纵向节径比(S_L/D)和横向节径比(S_T/D)对顺排管束传递特征的影响,并综合评价了其换热效果。结果表明,管排结构中大部分圆管均受尾流影响,在前列圆管尾流的作用下,尾涡结构被改变。相比于单圆管,管排结构强化传热的作用较明显,且对管排横向间距变化较为敏感。结果显示,增大S_L/D和S_T/D均能提高努塞尔数(Nu),但同时阻力系数也随之增大。当S_L/D=2.7,S_T/D=2.5时,管束的综合传热能力最优。此外观察到管排结构可以减小涡脱落频率,对涡脱落起一定的抑制作用,并通过线性拟合得到了管排斯特劳哈尔(St)数与雷诺(Re)数的关系。     

BCC_CSM气候模式对中国区域500 hPa气候要素模拟能力的检验

以NCEP/NCAR的1948?2005年逐月再分析资料为实况场,利用标准化均方根误差、相关系数等方法对BCC_CSM模式模拟中国区域500 hPa不同季节的位势高度场、温度场、相对湿度场、风场等气象要素的能力进行检验.结果表明:1)位势高度场的标准化均方误差在中国区域的不同季节均小于1,对应相关系数均通过显著性水平α=0.05的显著性检验;2)温度场的标准化均方误差在中国区域不同季节均小于1,对应相关系数除了夏季的中国南方外均通过显著性水平α=0.05的显著性检验;3)相对湿度、经向风与纬向风在中国区域夏、秋、冬三个季节的标准化均方根误差都比较小,实况纬向风场与模拟纬向风场全年相关性都比较好,相关系数在中国绝大部分地区能通过显著性水平α=0.05的显著性检验;4)模式模拟出了中国不同地区各个气象要素的年内变化,除中国西南地区的经向风外,各气象要素实况场与模拟场的相关系数均通过水平α=0.05的显著性检验.     

Richtmyer-Meshkov不稳定的双椭圆气柱相互作用的PIV研究

激波冲击双气柱既包含了激波与界面的相互作用又包含了界面之间的相互作用,是探索实际应用中Richtmyer-Meshkov(RM)不稳定性演化的很好的桥梁.利用高分辨率的粒子图像测速(PIV)技术获得了4种不同初始中心间距的双椭圆气柱在平面激波冲击下的速度场、涡量场和环量,定量表征了界面相互作用对RM不稳定性演化的影响.界面相互作用较弱时,双椭圆气柱的演化由两对涡对结构主导,内涡弱于外涡,且随着中心间距的减小内涡的弱化愈加显著;界面相互作用较强时,双椭圆气柱演化模式发生转变,内涡消失,只由一对涡对结构主导,外涡随中心间距的减小而增强.界面相互作用导致双椭圆气柱的演化行为与单椭圆气柱的演化行为不同,典型的特征是内涡的削弱及涡对结构的旋转运动.双椭圆气柱演化行为与双圆气柱演化行为也存在差别,表现为强界面相互作用时外涡随中心间距的减小而增强.     

基于声扰动方程的矩形柱绕流噪声数值模拟

为了研究低雷诺数下钝体绕流噪声,提出一种混合声学预报方法.采用声扰动方程(APEs)结合浸入边界算法(IBM)对矩形柱绕流噪声进行模拟研究.首先通过方柱绕流模型对混合方法和边界条件进行验证,然后讨论不同迎角下方柱的流噪声特性及0°迎角下矩形柱长宽比对声场的影响规律.结果表明:方柱绕流仿真结果与文献计算值符合较好,说明本文方法有效;随着迎角增加,声压级和声传播角相应变化,30°迎角下声压级较大;矩形柱的辐射噪声低于方柱辐射噪声,并随着长宽比增大逐渐减小,主要是因为尾流区形成卡门涡街的位置向下游移动,降低了尾涡强度和升力变化;矩形柱尾涡脱落和升力变化是主要的噪声源.     

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

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

串列双桥墩紊流特性研究

采用概化水槽试验和数值模拟相结合的方法,建立了水流与桥墩的二维概化水槽模型。通过对迹线、尾涡强度及脱落频率的研究,分析了桥墩间距对桥墩周围紊流特性的影响。研究结果表明,串列双桥墩间距较小(L/D≤2.0)时,上游桥墩没有形成涡的脱落,上游桥墩的剪切层交替的附着在下游桥墩上,导致剪切层在运动之前变厚,尾涡脱落频率减小,紊动强度增大;当3.0≤L/D5.0时,上、下游桥墩均有涡的脱落,且桥墩间距越大涡脱落现象越明显,尾涡的紊动强度随间距的增大而减小;当桥墩间距增大到5.0后,上、下游尾涡的紊动及尾涡脱落类似单个桥墩。     

方柱绕流的大涡模拟

采用大涡模拟(LES)方法对雷诺数为2.2×104的方柱绕流流场进行了数值模拟.使用非交错网格的有限差分法,分别对准三维物理模型和真实三维物理模型求解不可压N-S方程,将沿流向方向方柱水平中心线上的时均速度的计算结果与实验数据进行了比较,结果表明,三维模型的模拟结果优于准三维模型的模拟结果.比较升力系数和阻力系数发现,与二维模拟(RNG)方法相比, 三维模拟的结果更加接近实验测试数据.     

螺旋线对斜拉索涡激振动特性影响的试验研究

为研究不同参数的螺旋线对斜拉索涡激振动特性和气动力的影响规律,进行了针对表面无螺旋线斜拉索和缠绕不同参数螺旋线斜拉索的节段模型风洞试验,螺旋线参数包括根数、直径和缠绕间距. 结果表明：斜拉索表面缠绕1根螺旋线的涡激振动振幅远大于2根或3根螺旋线;在亚临界雷诺数区,缠绕螺旋线斜拉索的平均气动阻力系数均小于无螺旋线斜拉索,最多可减小10.9%;缠绕2根螺旋线的斜拉索平均阻力系数整体比3根螺旋线小. 斜拉索表面缠绕间距s=8D（D为斜拉索直径）的2根螺旋线时,缠绕直径d<15 mm螺旋线的斜拉索涡激振动现象明显;d ≥15 mm的螺旋线可显著降低斜拉索涡激振动响应幅值,其脉动升力系数相较于无螺旋线斜拉索有明显降低,证明缠绕大直径螺旋线对涡激振动响应控制的有效性. 螺旋线的缠绕间距对斜拉索涡激振动的影响与螺旋线直径有关,当螺旋线直径较小时,缠绕间距对涡激振动的影响并不明显;当螺旋线直径较大时,间距的大小直接影响涡激振动的抑制效果.     

方柱非定常绕流与涡激振动的数值模拟

为了研究二维刚体方柱的涡激振动和非定常绕流,利用计算流体力学软件Fluent6.0,引入雷诺应力模型求解不可压粘性流体Navier-Stokes方程,将Newmark-β方法的代码嵌入用户自定义函数(UDF)同软件连接来求解方柱的动力响应.用分块的方法建立初始网格,方柱和流场的耦合作用通过Fluent的动网格技术来实现,进行雷诺数从7000到42000范围内的模拟.得到静止方柱的升阻力系数和涡脱频率,其中雷诺数为22000结果与试验一致;同时得到方柱振动的位移和速度时程,以及"频率锁定"、"拍"和锁定段的位移"失谐"现象.     

周期性环肋通道内的流动和换热

用非稳态数学模型对包含入口段的带肋环形通道内的流动和换热进行了数值模拟.在肋片参数Sp/D=0.187 5、Sp/D=1.125,雷诺数在400～800范围内,流动和换热都出现了自维持振荡特性,且入口段后,各几何周期的速度场、无因次过余温度场基本相同,随时间的变化规律也基本相同.肋片高度增加,平均努塞尔数和摩擦阻力系数增加.肋片间距增大,对摩擦阻力系数影响较小,对努塞尔数影响较大.雷诺数增大,努塞尔数增大,阻力系数减小.     

横肋管中的流动阻力

以三种不同直径的铜管和有机玻璃管构成实验管路系统,研究了横肋管中空气流动阻力特性。在峰高为(0.01～0.11)d,间距为(7～40)d,Re=400～100000范围内,试验提供了从层流到湍流范围内完整的阻力系数曲线。从阻力形成机理,分析了横肋的峰高k,间距s,厚度w等因素对阻力系数的影响。并得到了简明的数据关联式f_r/f_s=1+128(k/d)~(1.36)(s/k)~(-0.41)Re~(0.25) 式中f_r为横肋管阻力系数;f_s为光滑管阻力系数;d为管径。     

风场中两串列球形沙粒的阻力及流场特性的数值研究

利用数值方法对雷诺数Re≤300的风场中沿风向串列的两球形沙粒的阻力系数及流场特性进行了研究.结果表明:两串列球形沙粒,后沙粒对前沙粒的阻力系数影响不大,而前沙粒对后沙粒阻力系数的影响很大,这一影响随着两沙粒间距的增大而减小,并且两沙粒的阻力系数均逐渐趋近单颗沙粒的阻力系数;两串列沙粒的阻力系数均随流场雷诺数Re的增大...