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

利用计算流体力学软件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值增大。     

不同攻角下2:1矩形截面棱柱绕流的大涡模拟研究

本文研究了高雷诺数下(Re=24000),宽高比为2:1的矩形截面棱柱在不同风攻角下的绕流问题.基于大涡模拟理论和Smagorinsky亚格子应力模型,运用开源CFD软件OpenFOAM对Navier-Stokes控制方程进行求解,对均匀来流作用下的矩形截面棱柱绕流问题进行了三维数值模拟.通过对平均分量、脉动分量的比较,验证了数值模拟结果的准确性与可靠性,并给出了阻力系数、升力系数、斯特罗哈数随攻角的变化规律以及平均流场和湍流流场的流场特征.基于本文研究结果,通过开源CFD可视化软件ParaView给出了不同攻角下的流型变化,结合平均压力系数的分布,详细解释了不同流型的特征,这些结果验证了基于风洞实验所提出的假设流型,相关结论可为结构风工程设计提供参考依据.本文所采用的软件均为开源软件,其开源特性具有良好的科研价值.     

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

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

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

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

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

采用改进的延迟分离涡方法数值模拟了高雷诺数下的柱体绕流,包括单圆柱绕流、单方柱绕流、串列双圆柱绕流和串列双方柱绕流,研究了不同雷诺数下圆柱绕流与方柱绕流的水动力特性.计算结果与实验数据及其他文献的数值计算结果吻合良好,研究表明,单方柱绕流在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分别为26、200、1.4×105的三种典型流态下的单圆柱绕流进行了二维数值模拟研究。Re为26时应用层流模型直接求解N-S方程,而Re分别为200、1.4×105时使用大涡模拟的方法进行计算。数值模拟很好地再现了稳定的涡旋结构、周期性交替脱落的卡门涡街结构、不规则的涡旋结构,在此基础上分析了尾流结构的基本特征及其压强分布规律、平均的流场特性、积分参数(如升力系数、阻力系数、斯特劳哈尔等),并与有关研究成果进行了对比。研究发现,采用不同流动介质时流场特性有所差异,空气为介质时的计算结果更符合实验的成果,而水为介质时计算结果偏差较大,这主要是由尾流涡旋产生的不合理负压造成的。     

基于浸入边界算法的振动钝体绕流模拟研究

传统CFD方法在振动钝体绕流计算中常借助动网格技术,网格再生任务繁重。针对于此,本文利用可在静止网格中计算动边界绕流问题的浸入边界算法(IBM),编写数值模拟程序,分别对竖向强迫正弦振动方柱(Re=UD/v=103、振幅恒定、振动频率变化)以及桥梁断面(Re=UB/v=7.5×103、振幅、振动频率均变化)展开气动特性和流场特征结构分析。初步研究结果表明,振幅恒定为方柱高度的14%时,其涡脱锁定区长度为0.06~0.2,锁定区后端(Stc0.2)振动方柱涡脱频率回归静止涡脱频率;不同振幅下的桥梁断面阻力系数均在静止涡脱频率处产生峰值,桥梁断面升力系数则在此处均出现归零效应,且振幅越大,归零效应愈明显。     

二维后向台阶流流动特性的实验研究

利用PIV技术,通过系统对150相似文献   

上游静止方柱尾流对下游方柱体尾激振动效应影响

基于半隐式特征线分裂算子有限元法,对低雷诺数下串列布置上游静止方柱–下游双自由度运动方柱体结构的尾激振动问题进行了研究.首先与现有文献结果进行对比验证该方法的正确性.然后着重分析了雷诺数(Re)与折减速度(U_r)两个关键参数对下游方柱尾激振动响应的影响,同时将计算结果与单方柱工况进行了对比.数值计算结果表明,雷诺数和折减速度对下游方柱的振幅、振动频率和运动轨迹等动力响应特性的影响较大.随着雷诺数的增大,双柱系统的互扰效应从以涡激效应为主逐渐转变为尾激效应发挥主导作用,从而导致下游方柱的振动响应增强.单方柱工况结构运动轨迹均呈"8"字形.然而,下游方柱的运动轨迹会随着雷诺数的增加而变得复杂.雷诺数较小时(Re=40, 80),下游方柱的运动轨迹基本为"8"字形.雷诺数较大时(Re=120, 160, 200),下游方柱的运动轨迹会出现双"8"字形.同时,下游方柱的尾流场特性主要呈现2S, 2S*, 2P,2T, P+S和稳态6种模式.最后,通过对流场特性进行分析,揭示了串列双方柱系统尾激振动效应的作用机理.     

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

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

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

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

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

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

大湍流度高雷诺数时并列双圆柱的平均和脉动压力分布

本文通过风洞实验研究了来流湍流度,Iu=10％雷诺数分别为Re=1.95×10~9和Re=6.5×10~5时单个圆柱和不同间距比下并列双圆柱的平均和脉动压力分布。结果表明:在Re=1.95×10~5时单个圆柱的平均压力分布类似于低湍流度高超临界雷诺数时的压力分布;当雷诺数增大至6.5×10~5时,绕圆柱表面流动的分离点前移和背压绝对值提高,总的阻力系数随之增加。并列双圆柱的间距比变化对圆柱表面压力分布影响很大,在极小间距比(N/d=1.05)时,双圆柱间的缝隙流使附近柱面产生高达-5的压力系数峰值(Re=6.5×10~5),同时脉动压力也大为增加;在较小间距比时(1.5     

两不等直径圆柱高雷诺数时均风力场的数值模拟

圆柱间气动干扰研究具有重要的理论和现实意义.尽管国内外开展了圆柱组合风效应的风洞试验研究,但主要针对等直径圆柱,并且雷诺数多为10~5以下.考虑到工程结构风场的高雷诺数特征,采用数值模拟方法,模拟两不等直径圆柱在串列、并列及交错排列下的高雷诺数(Re=4.5×10~5)时均绕流场.通过改变组合的间距和风向,分析探讨两柱阻力、升力及总风力的变化规律.     

地面效应影响下的高速火箭橇气动特性数值与风洞试验研究

精确预示地面效应下高速火箭橇的气动特性及流场规律对高速火箭橇的设计和评估具有重要意义。本文应用有限体积方法,研究了湍流模型对火箭橇气动特性计算精度的影响,建立了基于realizable k-ε湍流模型的火箭橇气动特性的高精度数值计算方法;结合风洞试验方法,研究了雷诺数和地面效应对高速火箭橇流场流动规律的影响,分析了火箭橇气动特性。结果表明,火箭橇阻力系数随雷诺数增大而减小,升力系数和俯仰力矩系数随雷诺数增大而增大,但雷诺数对高速火箭橇气动特性的影响较小;地面效应会使火箭橇流场发生激波-激波干扰、激波-边界层干扰和激波反复反射等复杂气动现象,大幅提升了火箭橇的升力系数和俯仰力矩系数,但对阻力系数的影响较小。研究为高速火箭橇气动外形的设计及运动稳定性的评估提供依据。     

圆角化对受迫振动方柱绕流特性影响机理分析

为分析圆角化对低雷诺数下受迫振动方柱绕流特性的影响机理, 对Ansys Fluent软件进行二次开发, 即通过用户自定义函数中的DEFINE_ CG_MOTION宏对柱体周期性受迫振动的函数进行编程, 并对流场计算域进行区域划分以便利用动网格技术中动态层法实现柱体受迫振动, 从而实现对受迫振动柱体绕流流场的流固耦合模拟.在雷诺数Re = 200时, 考虑方柱截面不同圆角的影响, 对均匀流作用下5种圆角化r/D = 1/2, 1/4, 1/5, 1/8和0受迫振动方柱的绕流进行数值模拟, 分析了这5种参数下受迫振动方柱的升阻力系数、尾流涡量和锁定区间的变化规律, 澄清了圆角化对受迫振动方柱稳定性的影响机理.研究表明: 与尖角方柱相比, 圆角化方柱升阻力系数有了明显的减小, 且升力、阻力系数随圆角增大而减小; 低振幅比下圆角方柱的涡旋脱落模式均为2S模态, 涡旋尾迹变窄; 锁定区间范围基本关于F = 1对称, 锁定区间的变化趋势与圆柱类似.      

电磁极宽度对圆柱绕流场影响的数值分析

在雷诺数Re=200的情况,利用Maxwell方程直接数值计算表面包覆电极与磁极圆柱体产生的电磁力分布,将其加入到动量方程中,在各种电磁力作用参数和电磁极宽度的组合下,对表面覆盖电磁极圆柱体在弱电解质中的绕流场结构及其升阻力特性进行了数值模拟与分析.结果表明,当电磁极宽度较小时,圆柱体绕流场的分离点越容易接近后驻点,而电磁力对总阻力的影响并不明显,但对压差和摩擦阻力均有明显影响.当电磁极宽度较大时,圆柱体尾部区域越容易产生射流现象,而且总阻力随电磁力作用参数和电磁极宽度增大而减小.在电磁力尚不足以完全抑制周期性涡脱落的情况下,升力幅值随电磁力作用参数增大而减小,但随电磁极宽度则先减小后略有增加,升力脉动频率则均随电磁力作用参数和电磁极宽度增大而增加.研究表明,电磁力可以有效地改善圆柱体绕流场结构,达到减小圆柱体阻力并抑制其脉动升力之目的,因此是圆柱型结构的一种有效流动控制手段.     

基于雨燕翅膀的仿生三角翼气动特性计算研究

针对低雷诺数微型飞行器的气动布局,设计出类似雨燕翅膀的一组具有不同前缘钝度的中等后掠(Λ=50?)仿生三角翼.为了定量对比研究三角翼后缘收缩产生的气动效应,设计了一组具有同等后掠的普通三角翼.为了深入研究仿生三角翼布局的前缘涡演化特性以及总体气动特性,采用数值模拟方法详细地探索了低雷诺数(Re=1.58×104)流动条...     

棱柱绕流流动的数值模拟

利用数值方法对长宽比为1／3, 1和3的棱柱绕流在雷诺数为100的非稳态流动特性进行了分析和研究。采用有限体积法对棱柱绕流的二维流动N-S方程进行离散求解,分析和研究了非稳态的棱柱绕流流场,升力系数,阻力系数和涡动特性,数值模拟的结果与相关文献的数据比较吻合。通过上述研究能够为了解棱柱绕流的非稳态流动特性提供有力的帮助。而对棱柱三维流动的模拟分析和对雷诺数的变化对棱柱流动特性的影响进行研究,将为掌握棱柱绕流的工程特性打下基础。     

对称翼型低雷诺数小攻角升力系数非线性现象研究

采用Rogers发展的三阶Roe格式,求解非定常不可压N-S方程,时间方向为二阶精度双时间步方法, 数值模拟了对称翼型SD8020低雷诺数(Re=40000,100000)条件下,流场层流分离涡结构和升力系数随攻角的变化．同试验比较证明了数值模拟的正确性．通过对数值模拟时均化流场结果的详细分析,发现对称翼型在小雷诺数0°攻角附近出现的层流分离泡,其内部结构和演化规律都不同于经典层流分离泡模型,从而提出了一种后缘层流分离泡模型．并应用该模型对对称翼型小攻角低雷诺数流场特性以及升力系数非线性效应的形成机理进行了研究和解释．