利用分区ALE算法数值模拟圆柱湍流涡致振动

用分区算法结合任意拉格朗日-欧拉法(ALE)数值模拟了圆柱湍流涡致振动.求解基于非交错网格系统,压力求解采用压力泊松方程提法,湍流模型采用标准k-ε模型和重整化群(RNG)k-ε模型.计算取中等雷诺数Re=5000、10000、15000、25000、50000等,质量系数M=10.阻尼系数ζ=0.00331,自振频率fn=0.18315、0.1628.计算结果表明:湍流涡致振动下圆柱时均阻力系数大于孤立圆柱绕流,而升力系数(振幅)值都小于孤立圆柱绕流.随着雷诺数增大,湍流粘性系数随之增大,但湍动能和湍流耗散率变化趋势不明显.对孤立圆柱绕流,研究结果与前人实验结果基本相符.     

谱单元法及其在多圆柱绕流分析中的应用

系统阐明谱单元方法,基于谱单元方法对低雷诺数Re=200时不同间距下的顺排两圆柱和Re=150正方形排列的四圆柱绕流及其阻力系数、升力系数等进行数值模拟。研究比较不同间距比L/D(两圆柱圆心距离与圆柱直径之比)对两圆柱和四圆柱绕流的影响,计算分析了涡量图分布、平均阻力系数和斯托罗哈数随间距比的变化。研究表明,间距比对顺排两圆柱和正方形四圆柱绕流影响显著;顺排两圆柱绕流存在临界间距比,在Re=200时临界间距比约为3.6。正方形排列四圆柱存在三种流态。当流场从一种流动形态变成另外一种流动形态时,力学参数发生显著变化,在某些间距比区间内出现骤升或骤降现象。     

多孔介质空心圆柱绕流的摄动解

本文主要讨论了绕多孔介质空心圆柱的有势流动问题。运用摄动方法,给出了速度势的四阶近似,然后得到了圆柱表面的无量纲的压力系数和阻力系数。结果发现在势流的情况下多孔介质圆柱受到附加阻力作用,并且在圆柱的喷射或吸收流体的特征速度与来流特征速度的比值是一个小量的时候,阻力系数与特征速度比成简单的线性关系。     

N-S方程的非定常绕流计算——有迎角Joukowsky翼型的粘性绕流

本文利用N-S方程计算了绕Joukowsky翼型的粘性流动,迎角为0°和42°,对称翼型的厚度为25.8％,无量纲时间间隔为t=0到t=10。计算结果包括:流函数分布,翼面压力分布和升力系数以及阻力系数。     

不可压缩翼型绕流数值研究

使用标准k-ω模型及与色散模型相耦合的k-ω模型分别计算了NACA0012翼型和NACA 4412翼型的低速绕流问题.NACA 0012翼型计算了其来流雷诺数为2.88×106,攻角从0°到15°范围内的流动结构、翼型表面压力分布和升力、阻力特性;NACA 4412翼型计算了临界雷诺数为1.52×106,攻角为13.87°时的流动分离和翼型表面压力系数 ,并与实验数据进行对比.结果表明:在同等条件下,使用与色散模型相耦合的k-ω模型计算得到的NACA 0012翼型的升力和阻力系数比标准k-ω模型提高精度约5%,NACA 4412翼型的表面压力系数精度提高了约3%,进一步验证了其可信性,可将其进一步应用到低速飞行器的气动计算中.      

龙脊风帆模型气动特性风洞实验

通过风洞模型实验，利用六分量高频底座天平技术，测量了称为龙脊风帆的一种单帆模型在均匀来流中各个风向角和俯仰角下的时均升力系数和时均阻力系数以及脉动升力系数和脉动阻力系数，可为工程中类似帆形结构的设计和受力计算提供参考。     

浸入边界法模拟小圆柱对主圆柱涡脱落的抑制

使用浸入边界法研究了小圆柱对主圆柱涡脱落的抑制.方法中使用非贴体笛儿尔网格,易于处理包含复杂边界的流动问题.采用离散附加力直接加入边界条件方法对虚单元进行重构,使边界条件在浸入边界上精确满足.使用隐式分步法解二维非定常不可压Navier-Stokes方程,通过速度和压力解耦提高计算效率.数值模拟单圆柱绕流及不同位置小圆柱和主圆柱的流动干扰,通过分析流场涡结构和升、阻力系数,得到小圆柱对主圆柱涡脱落的延迟和抑制作用.计算结果与已有实验结论和数值结果对比,计算误差不超过5%,说明浸入边界法可以简单有效地处理圆柱涡脱落抑制这类流动干扰问题.     

桥梁断面雷诺数效应

雷诺数效应是流体力学的传统基础性问题,随着桥梁跨度的增加,桥梁断面雷诺数效应研究日显重要.在近两个数量级的雷诺数范围内测量了两类桥梁断面的三分力,以及表面压力.研究了桥梁裸体断面三分力系数及表面压力系数的雷诺数效应,宽高比对裸体断面三分力系数,以及表面压力系数雷诺数效应的影响.近流线型桥梁断面的阻力系数、升力系数都有明显的雷诺数效应,Ⅱ型断面阻力系数雷诺数效应与近流线型断面规律相反,雷诺数对两类断面的表面压力分布具有较大影响.宽高比对三分力系数雷诺数效应有明显的影响.     

基于N-S方程的非定常地面效应研究

研究了二维翼型在移动地面和固定地面条件下的地面效应及流动特性问题.流场求解采用雷诺平均N-S方程,湍流模型采用k-ω和SST模型.对NACA0015对称翼型在离地速度为1 m/s情况下的非定常绕流进行了计算,选取移动地面条件下不同高度的定常解与非定常计算结果进行了比较.结果表明,考虑非定常效应时升力和阻力的变化明显小于定常解的结果,相对而言力矩的差别较小,固定地面条件下的升力系数和阻力系数略大于移动地面条件下的计算结果.     

舰载机折叠翼气动特性分析

折叠翼技术是舰载机与航母相匹配的关键,研究折叠翼的气动特性对舰载机飞行动态及安全性具有重要意义。基于离散化的思想对折叠翼进行建模,并利用计算流体力学CAE软件Fluent对其进行数值计算,从升力系数和阻力系数变化及压力分布分析舰载机折叠翼的气动特性。结果表明:当折叠角为75°时,外翼的升力系数及侧向力系数将达到最大,阻力系数最小;一旦折叠机构失效,外翼的折叠将导致机翼整体升力系数及阻力系数大幅度下降,气动特性变差;与折叠机构失效前相比,折叠后的机翼失速迎角变小,且在失速后升力系数有缓慢上升趋势。     

Experimental Investigation of the Effects of Various Plasma Actuator Configurations on Lift and Drag Coefficients of a Circular Cylinder Including the Effects of Electrodes

In this paper, the effects of the existence of plasma actuator electrodes and also various configurations of the actuator for controlling the flow field around a circular cylinder are experimentally investigated. The cylinder is made of PVC (Polyvinyl Chloride) and considered as a dielectric barrier. Two electrodes are flush-mounted on the surface of the cylinder and are connected to a DC high voltage power supply for generation of electrical discharge. Pressure distribution results show that the existence of the electrodes and also the plasma are able to change the pressure distribution around the cylinder and consequently the lift and drag coefficients. It is found that the effect of the existence of the electrodes is comparable with the effect of plasma actuator in controlling the flow field around the cylinder and this effect is not reported by other researchers. Eventually it is concluded that the existence of the electrodes or any extra objects on the cylinder and also the existence of the plasma are capable of changing the flow field structure around the cylinder so that the behavior of the lift and drag coefficients of the cylinder will be changed significantly.     

Viscous Crossflow Around a Rotating Circular Cylinder

Unsteady viscous crossflow around a rotating circular cylinder is studied in the laminar and turbulent flow regimes. The dependences of drag and lift coefficients of the cylinder on the rotation speed are obtained.     

Two-dimensional viscous vortex flow around a circular cylinder

A code based on a vortex method for simulating viscous vortex flow around a circular cylinder has been developed. The original method to obtain the no-through and no-slip condition introduced here is described in detail. The diffusive part of the vorticity transport equation is treated using a deterministic method based on the solution of the heat equation. Results have been obtained and are presented concerning the vorticity field, the velocity field, the evolution of drag and lift coefficients and the Strouhal number.     

吸气条件对圆柱非定常分离流影响的数值研究

本文根据二维可压非定常Ｎ－Ｓ方程，基于Ｂｅａｍ ａｎｄ Ｗａｒｍｉｇｎ及Ｓｔｅｇｅｒ的方法，在圆柱绕流出现周期分离后，施加吸气条件，研究气对于圆柱非定常绕流的影响，计算结果表明，在恰当地区域施加恰当的吸气条件，可以希望通过减小压差阻力来减小总阻力，吸气还有抑制尾迹中非对称流动的作用，从而使圆柱所受的垂直于来流方向的升力的变化幅度减小，变化频率减小，但是吸气会增加磨擦阻力，过份的吸气会导致摩阻增大太     

垂直交叉双圆柱绕流强迫振动数值模拟

采用虚拟体法数值模拟了平面单圆柱绕流强迫振动,以及空间垂直交叉双圆柱绕流下游圆柱强迫振动流场。验证了单圆柱强迫振动中的锁定状态以及相位突变现象,从而证实了该数值方法模拟振动流场的可靠性。研究了在雷诺数Re=150、间距为5倍圆柱直径、下游圆柱按正弦曲线振动时,对尾流场的影响。下游圆柱两端尾流场在振动的作用下,涡的横向间距增大,而中心尾流场由于受上游圆柱尾流的影响而保持原先的状态。下游圆柱在锁定区的振动使尾流场变得稳定,表现在流向二次涡结构的减少,以及尾流场的速度场频谱趋向单一化。     

趋于临界马赫数的圆柱跨声速绕流特性分析

通过深化认识趋于临界马赫数Ma_cr的圆柱跨声速绕流特性,明确新型飞行器增升减阻设计的空气动力学理论依据。采用大涡模拟方法数值研究了来流马赫数Ma_∞为0.75和0.85、雷诺数Re为2×10⁵的圆柱跨声速绕流。结果表明:当Ma_∞趋于临界马赫数(Ma_cr≈0.9)时,圆柱的阻力下降且升力系数振荡被抑制;通过力的分解,得知圆柱的阻力减小来自旋涡力的影响,而非可压缩性;圆柱的阻力减小与其背压上升有关;剪切层初始阶段的对流马赫数Ma_c随Ma_∞的增加而增大,而增长率相反,这使得剪切层更为稳定、柱体背压更高。此外,由于Ma_∞=0.85时边界层分离点处的激波和尾迹处的激波向下游推移,使得近尾迹处的湍流脉动减弱,进而导致柱体的表面压力振荡和升力系数振荡被抑制。     

临近空间螺旋桨气动性能数值分析

在实验条件不成熟和计算性能不够好的情况下,提出了一种研究螺旋桨气动性能的方法.首先对叶素周围的流场进行数值模拟,分析叶素的气动性能并得出叶素升阻系数,然后差值得到叶素升阻系数沿桨径的连续分布,再利用叶素理论对螺旋桨的整体气动性能进行分析.研究表明,所提出的研究方法不仅可以对叶素周围流场进行详细的分析,还能够较好地分析螺...     

各种排列双圆柱绕流相干性研究

本文讨论了各种排列双柱绕流的相干性问题。低湍流度高Re数的实验表明升力是由尾涡诱导引起的,而双稳态现象是由尾涡不稳定性引起的。此外,还对高湍流度状态作了探索性的工作,发现由于尾涡的可分性和狭缝流方向的不稳定性所共同引起的双稳态现象与低湍流度的情况有很大不同,升力和阻力也是如此。     

Simulation of the flow around a stationary circular cylinder taking into account laminar-turbulent transition

We consider a problem of a stationary incompressible viscous fluid flow around a flat circular cylinder. In the vicinity of the critical Reynolds number Re _cr a stepwise drop of the cylinder drag takes place, which is called the drag crisis.