用格子Boltzmann方法模拟方柱绕流

基于格子Boltzmann方法的原理和对边界条件处理的便捷性,分别对单方柱和并列双方柱绕流进行了研究.对于单方柱绕流,分别研究了不同雷诺数和阻塞比对它的影响,给出了不同雷诺数下单方柱绕流的流线图和涡量图,以及平均阻力系数和最大升力系数随雷诺数的变化曲线,得到单方柱产生卡门涡街的临界雷诺数,也给出了不同阻塞比下单方柱绕流...     

大型冷却塔双塔组合表面风压干扰效应试验

针对圆截面双曲线冷却塔双塔组合展开刚体模型测压风洞试验,研究双塔在不同间距不同来流方向下的表面动静态风压分布情况.通过与单塔刚体模型测压试验结果进行对比,分析双塔组合风压干扰因子和等效风荷载比例系数随塔距及风向角的变化关系.同时,对双塔在几种典型布置形式下的表面动静态风压分布和绕流状态进行了深入的分析.研究发现:即使塔...     

正三角形布置下三圆柱绕流的数值模拟

利用有限体积法和非结构性网格对正三角形排列下的三圆柱绕流进行了数值模拟.取雷诺数Re=100,并考虑4个不同的间距比(1.5,3,5,7),着重研究了间距比对绕流流场特性的影响,结果表明:在小间距比的情形,下游圆柱发生的偏流现象为单稳态偏流;在大间距比的情形,尾流效应使得上游圆柱脱落的剪切层在下游圆柱间上下波动.     

平板间方柱绕流的格子Boltzmann方法模拟

为研究平板间方柱绕流上下平板对置于其中的方柱绕流所产生的影响,采用格子Boltzmann方法对二维平板间低雷诺数(Re=100)方柱绕流问题进行了数值模拟研究.分析了3种不同阻塞比下,平板边壁对方柱的升、阻力系数、Strouhal数和尾涡流场的影响.结果表明:平板对方柱绕流特性有明显的影响,随着阻塞比的增加,阻力系数和Strouhal数均增大,与无边壁相比阻力系数可增加达30%,而升力系数却随之减小.计算结果与相关实验数据相吻合,验证了格子Boltzmann方法对钝体绕流非定常问题模拟的有效性。     

合成射流对串列双圆柱受力特性及流场特性影响的数值模拟

合成射流能够通过改变单个圆柱尾部的旋涡脱落模态有效抑制其涡致振动,但对多圆柱结构的控制效果尚待研究。采用数值模拟研究了对称布置在上游圆柱尾部的水平合成射流对串列双圆柱的受力特性和尾流特性的影响。在研究中,来流的雷诺数恒定为150,合成射流的位置角度为45°,2个圆柱中心之间的间距比为2～10。研究结果表明,水平合成射流能够同时降低上、下游圆柱的升、阻力系数,但过大的合成射流动量系数反而会导致上、下游圆柱的阻力系数均大于无合成射流控制情况下的值。当动量系数为2和4时,上游圆柱的尾迹完全由水平合成射流主导,变得完全对称;但下游圆柱的尾部始终存在交替脱落的旋涡。     

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

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

串列双圆球的流场转捩研究

本文对虚拟边界法加以改善并推广到三维多连通区域的数值模拟上去,研究了不同雷诺数下串列双圆球的流场转捩现象.结果显示小间距比串列双圆球绕流场的第一分叉点位于雷诺数150～175,第二分叉点位于雷诺数225～250.大间距比绕流场的第一分叉点位于雷诺数170～200.转捩点的位置与单圆球绕流相比有不同程度的提前.给出了切面内流态沿流向的发展变化,捕捉到极限环现象,发现下游尾流场的涡结构分别脱落自下游圆球的上表面和下表面,两个来源处的涡脱频率为倍频的关系,随着雷诺数的增加存在频率竞争和占优频率交替的现象.     

基于数值模拟与风洞试验的多分裂导线阻力系数干扰效应

为研究分裂导线阻力系数随风攻角与分裂间距的变化规律,采用CFD数值模拟方法研究了不同分裂数与分裂间距下的低风压导线阻力系数随风攻角的变化规律.通过在边界层风洞中开展90 cm长的导线节段模型试验,获得了低风压导线不同风攻角条件下导线阻力系数随雷诺数、尾流干扰的变化规律,并与数值模拟结果进行对比.研究结果表明:多分裂导线...     

高雷诺数下并列双圆柱绕流的大涡模拟

为澄清并列双圆柱结构发生偏向流现象的流场机理,采用大涡模拟(LES)方法,在高雷诺数下(Re=1.4×10~5)研究了并列双圆柱的气动性能及其流场特性随圆柱间距比P/D(P为圆心间距,D为圆柱直径)的变化规律,重点探讨了小间距并列双圆柱的偏向流现象及其对圆柱气动性能的作用机理.研究结果表明:大涡模拟方法得到的气动力结果与文献风洞试验值吻合良好;随着并列双圆柱间距的增大,绕流场会呈现单一钝体、偏向流和平行涡街等多种流态结构;当P/D=1.1时,绕流场会间歇性地出现单一钝体和偏向流流态,两种流态的气动性能和流场特性有很大差异,圆柱的气动力会随时间发生剧烈变化,呈现非稳态特征;当P/D=1.2～1.5时,绕流场呈现偏向流流态,两个圆柱的气动力和尾流呈现不对称现象,偏向流的偏转方向会出现间歇性地变化,尾流涡脱强度弱,气动力脉动小;当P/D=2～4时,绕流场总体呈现平行涡街流态,尾流涡脱强度强,气动力脉动大,气动干扰减弱.     

来流湍流度对圆柱绕流特性的影响

采用大涡模拟方法,在雷诺数为1.4?05时,研究来流湍流度对圆柱绕流特性的影响规律。研究表明：湍流度会使得自由剪切层间的相互作用更加剧烈,使得流场的三维特性更为强烈;同时加强了尾流区湍流能量的相互传递,减少能量的耗散,从而延长了回流区长度。随着回流长度的增加,圆柱背风面的平均风压系数得以增大,最终导致圆柱阻力系数随之减小。在尾流湍流得到增强的同时,尾流区会形成多尺度的湍流,从而促使漩涡形成区形成大涡量的漩涡,导致尾流区的漩涡呈现多频率的脱落。     

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

为探讨不同间距下,串列双圆柱绕流场特征的内在变化规律,在低雷诺数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种流态随间距比演变的内在机制;与单一钝体及剪切层再附流态相比,双涡脱流态的下游圆柱近尾流区的高阶模态结构更为复杂,流场重构时需要更多模态才能达到与其他两种流态相似的精度,在尾流旋涡影响区的重构误差最大。     

中等雷诺数圆球绕流的数值研究

采用三维黏性不可压缩流场分块耦合计算方法计算圆球绕流问题.动量方程分别采用高精度和常规精度离散,压力泊松方程采用二阶中心差分格式离散.出口边界条件为无反射边界条件;计算网格为全流场封闭的H-O型网格.计算了雷诺数为8.9～1 000.0时圆球绕流受到的阻力,分析尾涡形状及长度随雷诺数的变化规律.计算得到的阻力系数和尾涡长度与实验结果符合较好.用高精度算法研究了尾涡结构问题,在雷诺数为1 000.0时成功地获得了尾流发夹涡现象.     

三维错列双波浪锥柱绕流流动特性数值仿真

利用大涡模拟分析了雷诺数Re=3 900下间距比为L/Dm=4和5、交错角为α=0～15°的错列双波浪锥柱升阻力特性、流场结构及尾迹干涉效应。研究发现:对于特定间距比L/Dm=4、5,下游波浪锥柱脉动升力系数得到显著提高,在α=10°时较单直圆柱分别提升20.1倍和21.4倍,这主要是由上游波浪锥柱尾涡卷起撞击在下游波浪锥柱的一侧,下游波浪锥柱表面产生周期性的耦合力导致;受上游波浪锥柱两个自由端及其侧面来流与下游波浪柱作用产生的回流区的影响,下游波浪锥柱时均阻力系数显著降低。随着交错角增加,对于间距比L/Dm=4、5,下游波浪锥柱时均阻力系数逐渐接近上游柱,且在α=10°时,时均阻力系数较单圆柱分别降低34.9%和18.8%。涡量图展示了错列双波浪锥柱之间的完全撞击状态,侧面撞击状态和尾流干扰状态;且由于波浪锥柱表面形状的影响,使得流经波浪锥柱后方的尾涡存在明显分层现象,侧面撞击状态相对于其他两种状态能提供更大的脉动升力系数。本文研究结果可为风力俘能结构列阵的布局提供理论支持。     

地震作用下双塔弱连接连体结构受力性能研究

对单塔楼和双塔连体结构的抗震性能进行了比较分析和研究,采用三维有限元模型,通过3个算例,分别计算了单塔楼、连接体与主体结构固接的对称双塔楼和与主体结构弱接的对称双塔楼结构,从结构振型、位移、剪力3个方面分析了弱连接体对双塔连体结构抗震性能的影响,得出了一些有益的结论,可为同类工程设计参考.     

基于动态亚格子模型的方柱绕流大涡模拟

目的 对高雷诺数下的方柱绕流进行数值模拟,得到方柱下流场时间平均结果以及湍流流场的流场特征,用于验证亚格子模型对大涡模拟(LES)计算精确性的影响.方法 基于有限体积法(FVM)的流体计算开源软件OpenFOAM,采用两种不同动态亚格子模型下的大涡模拟方法,模拟高雷诺数(Re =2.2 ×104)下的方柱绕流.结果 通过数值计算得到流场完整的平均和脉动结构,并得到阻力系数、升力系数、斯托拉哈数等重要的气动力参数,经对比与试验实际情况吻合.结论 高雷诺数复杂流场下亚格子应力的封闭形式对大涡模拟的计算结果有重要的影响,采用动态模拟技术一方程模型计算结果优于常规的Smagorinsky模型,在高雷诺数复杂流场下的计算是高效的,为结构抗风的数值研究提供了参考.     

刚性薄平板绕流特性的数值模拟

为了揭示雷诺数、来流剪切参数对平板尾流涡街演化过程以及平板阻力特性的影响规律,本文基于浸没边界-格子Boltzmann方法建立了刚性薄平板绕流模型,并在雷诺数为100～1 200内,对不同雷诺数与来流剪切参数下的刚性薄平板绕流问题进行数值模拟。分别研究了平板尾流区沿流向的2次失稳跃迁过程以及平板流向载荷随来流剪切参数的变化及其对雷诺数的依赖性,并对其机理进行了讨论。研究表明：随着流场剪切参数的增加,平板近尾流经过数次准周期分岔过程进入了由大尺度涡结构所主导的混沌状态,进一步诱发了低频、混沌的结构阻力;雷诺数越大,向混沌特征流向载荷转变所需的临界剪切参数越小。     

水下钝体空化绕流数值模拟

基于均质平衡流理论,利用Fluent 6.3对水下钝体空化绕流进行数值模拟,研究钝体空泡形态的变化规律,分析钝体的减阻特性。结果表明:空化数较小时,钝体从水绕流发展成空化绕流,随着空化数减小,空泡的无量纲直径和无量纲长度都增大;压差阻力系数减小,摩擦阻力系数减小并趋于0,总阻力系数比水绕流时降低,基本等于压差阻力系数。     

圆柱绕流涡致振动的平面湍流数值模拟

为了研究工程中存在的湍流涡致振动问题,应用平面湍流k-ε模型和重整化群k-ε模型,结合分块耦合方法及任意拉格朗日欧拉法（ALE）数值模拟了圆柱绕流涡致振动,求解基于非交错网格系统, 其中包含满足连续性方程的压力泊松方程解法.计算了湍动能、湍流耗散率、湍流黏性系数及柱面涡量的变化情况.结果表明,当雷诺数为5 000～2×105时计算所得的阻力系数与实验数据吻合良好.研究发现,在涡致振动情况下的湍流耗散率比固定圆柱涡脱落情况下小得多,湍流黏性系数与雷诺数呈对数律增长.     

带襟翼的机翼尾涡合并数值计算

为获得带有襟翼的机翼尾涡的合并动力学过程,在验证数值方法的基础上,数值模拟带有襟翼的机翼绕流尾流场,根据涡合并特征将合并过程划分为诱导共转阶段、合并阶段和轴对称化阶段.采用涡间距量纲—化尾流区域描述二涡诱导合并.变换弦长雷诺数、襟翼翼梢与机翼翼梢的间距、襟翼角度,改变襟翼翼梢涡与机翼翼梢涡的强度比,得到尾涡合并差异的特征.计算结果表明:随着雷诺数的增加涡的强度增加,涡量的扩散程度减低,涡合并过程被推迟,空间诱导运动过程得到延长,涡系空间诱导运动增强,涡合并的雷诺数效应随着雷诺数增加而减弱,缩短涡间距,加速涡合并过程,但合并后的远场涡尺寸和形态没有显著改变;襟翼涡随着襟翼角的减小而减弱,强度逐渐减弱的襟翼涡逐步被翼梢涡拉伸卷吸,微弱的襟翼涡系在近场中的合并过程完全改变,翼梢涡的运动轨迹并未受到诱导运动.翼梢涡合并的雷诺数效应表现为诱导运动过程的增强,尾流中强度小的涡系起不到明显的诱导作用.     

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

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