扇环形截面弯曲管道内的粘性流动

本文利用涡量－流函数法系统地研究了扇环形截面弯曲管道内的二次流动，计算结果表明了采用控制体积积分法对涡量－流函数方程进行离散来求解弯曲管道内粘性流动问题是成功的，这样可以避免对压力项的处理；对所得结果的分析揭示了该二次流动不仅受雷诺数和曲率影响，而且扇形中心角、内外径之比对该流动的影响也是显著的。     

颗粒两相流自由剪切层大涡模拟

采用大涡模拟(LES)方法数值仿真了可压缩流自由剪切层内颗粒两相流流场结构.气相控制方程采用二维盒式滤波法滤除小于网格尺度的小涡以得到LES控制方程,而流场中小涡对大涡的影响采用K-ωSST湍流模型模拟.在Lagrange坐标系下用颗粒轨道模型仿真颗粒运动,用双线性插值方法实现Euler坐标系到肌Lagrange坐标系的转换,并考虑了颗粒速度和温度的改变对气相流场的反作用.分析了仿真结果中气相自由剪切层的涡结构特性,不同St数颗粒在流场中的运动特性,颗粒对流场的反作用效果,得到了合理的结论.     

平板湍流转捩过程中Λ涡和环状涡的周围流场研究

采用高精度湍流直接数值模拟技术,对平板流动转捩过程中的涡系结构,尤其是Λ涡和环状涡附近的流场进行了详细的研究.高分辨率的数值模拟给出了实验手段无法得到的全流场信息,以及雷诺平均纳维斯托克斯方程方法(RANS)无法给出的精细涡结构信息.数值模拟在精确捕捉Λ涡和环状涡等转捩过程中典型涡系结构、上喷和下扫等典型流动现象的同时,发现了一些与传统认识所不同的现象:1)Λ涡并不是当地涡量集中的地方;2)环状涡的附近存在一个速度显著大于远场来流速度的区域.这些现象对传统湍流转捩理论提出了挑战,并引发人们对湍流生成与维持机理的进一步研究.     

基于涡限制法的旋翼气动噪声计算

在直升机旋翼流动的数值模拟中,流动求解精度对旋翼气动噪声的预测有着重要影响.采用中心格式有限体积法求解可压缩Euler方程数值模拟旋翼绕流,并采用基于Ffowcs Williams- Hawkings( FW-H)方程的声类比方法计算旋翼气动噪声.为了抑制中心格式的数值耗散,保持机翼翼尖以及旋翼桨尖的尾涡结构,引入了涡...     

翼型近尾迹流中大尺度涡结构的理论模型

根据流动稳定性理论相干结构的相关理论 ,提出了一种二维波模型来描述翼型近尾迹流动中的大尺度涡结构 ,通过采用对近尾迹流动区域分区处理的方法计算了近尾迹区域的大尺度涡结构。计算所得的流线图 ,涡量分布以及雷诺应力分布具有与实验结果相似的分布规律。说明该二维模型能够较好的描述近尾迹流动中这种大尺度涡结构的运动学和动力学特性      

等速上仰翼型绕流结构的数值模拟

通过数值求解涡量-流函数形式的N-S方程研究等速上仰翼型的分离流动结构。用ADI方法解涡量方程;用Poisson方程直接法解流函数方程。在计算得到的流场中除可看到实验中观察到的多涡结构外,还发现一些实验中未观察到的现象。将分离流动结构的计算结果与翼型气动力变化相联系,并结合涡动力学理论,定性地解释了上仰翼型产生高升力的原因。     

细长尖头旋成体大迎角非对称涡系结构

本文采用数值计算方法,对细长尖拱旋成体大迎角背风侧非对称涡系结构以及与沿轴向交变的侧向力分布的关系进行了研究.求解的是N-S方程,采用Jameson中心格式,湍流采用修正的B-L模型.通过数值模拟揭示了在头部开始产生非对称涡,形成二次涡的机理以及在背风侧形成集中涡的过程.从而说明了采用集中涡来模拟这类流场的正确性,表明了二次涡对多个脱体集中涡形成的贡献.     

基于黏性涡模型的旋翼流场数值方法

 建立了一种适用于旋翼非定常流场特性分析的黏性涡数值方法。在该方法中:流场中的大尺度涡被离散为若干微小的涡元,通过求解涡量-速度形式的Navier-Stokes方程模拟涡元的输运等过程;黏性扩散效应采用高精度的粒子强度交换法进行计算,而桨叶附着涡以及新生涡环量采用了Weissinger-L升力面理论进行求解;为显著提高计算效率,在诱导速度及其梯度的计算中还引入了快速多极子算法(FMM)。应用上述方法,对悬停和前飞状态下的多个旋翼流场算例进行了计算,通过对比旋翼尾迹涡量特征和诱导速度分布等,验证了该方法的有效性。此外,还将本方法与旋翼计算流体力学(CFD)方法及传统的自由尾迹方法进行了比较,结果表明黏性涡方法在兼顾效率的同时,还能够更好地捕捉旋翼尾迹运动。     

基于黏性涡域法的二维钝体绕流模拟

针对传统计算流体力学方法需要划分网格且计算耗时长等问题，一种确定性涡方法(黏性涡域法)逐渐发展并日臻成熟。本文利用该方法计算模拟了二维钝体的流场、气动系数及斯特劳哈尔数，并与其他成熟CFD软件计算结果对比来验证模拟精度。黏性涡域法的模拟过程为生成涡量、涡量演变及气动计算。文章首先推导了该方法各过程的基本方程，然后基于Lua语言编制了圆柱和方柱绕流计算程序，最后利用Gnuplot实现流场可视化。利用黏性涡域法模拟不同雷诺数下圆柱绕流、不同无量纲频率和无量纲振幅的振动圆柱气动力时程、不同风向角下的方柱绕流，结果表明：雷诺数为2×10²～1×10⁵时，基于黏性涡域法的圆柱绕流气动计算结果与其他方法的计算结果基本一致；基于该方法的不同风向角下方柱绕流气动计算结果与其他文献结果趋势相同。     

三角翼大迎角绕流特性数值模拟的网格处理技术研究

研究网格对三角翼大迎角绕流特性数值模拟的影响.结果表明,网格的生成需要结合流动现象,对网格拓扑结构和网格点分布进行选择与搭配.C-H型网格适宜模拟尖前缘分离涡流态,法向网格在一定范围内应等距增长,沿流向逆压梯度较大的区域内适当增加网格点,尾迹区网格则应做上翘处理.Euler方程具有模拟三角翼旋涡及预测涡破裂特性的能力,但对二次涡等粘性引起的流动细节把握能力不足.利用层流假设的N-S方程,通过合适的网格,也可得到满意的计算结果,但对涡破裂后的强烈非定常湍流流动模拟能力不足.采用旋涡螺旋度可准确反映主涡与二次涡流动,描述旋涡的破裂现象.用轴向速度迅速减小并小于来流速度的点作为涡破裂判据似应更合理.     

波瓣混合流场流向涡的数值分析

借助N-S方程对波瓣喷管混合装置的流场进行了数值分析。采用三维贴体曲线坐标,非交错网格法和K－ε两方程湍流模型,对波瓣喷管内外两侧的气流通道及波瓣喷管下游的混合段流场统一进行了数值计算。得出了速度场和流向涡及正交涡的分布,并对流向涡的发展规律及其对正交涡分布规律的影响作了分析。并对混合管壁面的压力分布和出口速度型与实验值进行了比较。     

高压比离心压气机二次流旋涡结构研究

为揭示高压比离心压气机的流动特性,采用数值方法对高压比离心压气机的旋涡结构和流动损失的产生及演变规律进行了研究。根据不同类型旋涡的具体特征,给出了分别适用于受迫涡和自由涡的二次流识别方法,包括截面旋线法和拟定主流的截面流线法。应用给出的二次流识别方法并结合耗散函数,探讨了压气机内旋涡的形成机理以及旋涡与损失的关联性。研究表明:当涡量与截面法矢量夹角的余弦值大于零时,旋线方向与实际气流方向定性一致,否则相反;旋线显示的涡轴方向与截面法矢量夹角大于90°时,识别出的旋涡不存在;刮削涡和泄漏涡既是低能流体的聚集区也是能量的耗散区,是影响离心压气机损失产生及分布的关键因素;诱导轮尾迹会抑制导风轮流道内叶表通道涡的形成。     

预测轴流压气机叶栅通道端壁附面层及叶片力亏损

本文将二次流理论与端壁附面层理论相关联,提出一种用于预测轴流压气机叶栅通道内的端壁附面层及叶片力亏损的新方法。附面层主流流动采用包含叶片力亏损的动量积分方程,横流流动采用双层速度分布模型,使用二次流分析和有限差分计算方法预测横流外层速度分布,主流流动和横流流动交替计算。用于两种高负荷压气机叶栅的预测结果表明,不仅端壁附面层发展的预测结果与实验结果吻合较好,而且叶片力亏损发展的预测值也与实验值一致性较好;此外,该方法不仅自动,而且还能较准确地预测S形横向流动的发展。     

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.     

小振辐振荡圆柱诱导的二次旋涡流动

本文给出充满粘性不可压缩流体时同心圆柱间内柱作小振幅振荡所诱导的二次旋涡流动的解析解,对典型的柱径比给出容器外边界对二次旋涡流动流线、涡量分布和动能分布的影响。给出了等涡量线和等动能线分布的鞍点。采用数值求解非定常Stokes流动考察了二次旋涡流动的发展过程。二次定常旋涡流动的显示实验表明,解析解、数值计算和实验显示所示流谱基本一致。     

有限空间内低旋流流动与燃烧特性

为了对低旋流流动特性及其对火焰传播与稳定影响进行深入理解,通过实验与数值模拟相结合的方法,对有限空间内低旋流流动与燃烧特性进行研究.结果表明:改变中心孔板小孔直径,可以形成不同旋流数的旋流喷射器;在冷态条件下,随着旋流数增加,在展向平面上逐渐形成“W”形速度矢量分布;在铅垂截面内,速度矢量分布存在一个环形高速区域、中心低速区和外侧低速区.由于燃烧气体吸热膨胀、涡量诱导和温度斜压效果,使得冷态高旋回流流动转变为了燃烧条件下低旋无回流流动;在有限空间内,低旋流燃烧没有形成脱体火焰,而是形成类似于扩散燃烧火焰;随着旋流数增大,温度分布边界和化学反应边界在展向平面逐渐缩短,而在铅垂截面逐渐扩大.      

垂直于物面的横截面上流态的拓扑

对于旋成体有迎角绕流，本文研究了垂直于物面的横截面上的流态结果。利用解析方法，给出了物面轮廓线上、迎风和背风对称线上以及全横截面上截面流线方程奇点的分布规律以及奇点沿流动方向的演变规律，从而可定性地确定截面流态的结构。通过求解ＮＳ方程数值模拟了钝锥有迎角高超声速绕流流场，计算获得的垂直于物面的横截面流线的拓扑结果和理论分析完全一致。     

基于横向二次射流的水下推力矢量方法

提出一种基于横向二次射流的水下推力矢量技术,通过二次射流的横向速度场诱导主流发生偏转,建立了推力矢量偏角与流速偏角的数学关系,证明了通过主流偏转实现推力矢量偏转的有效性。通过数值计算方法分析了不同二次射流深度、不同二次射流/主流体积比及不同二次射流/主流速度比条件下主流偏转角度变化。结果显示:随着二次射流深度的增加,主流受到壁面阻碍作用增强,因而偏转角度减小。随着二次射流/主流体积比的减小,出口负压区所占比例减小,主流偏转角度增加,且当体积比减小到一定值后,负压影响可以忽略,主流不再随体积比而变化。主流偏转角度随速度比增加而增加,且在速度比一定的条件下,速度数值的变化对主流偏转没有影响。设计了一种主流为圆形射流的水下矢量推进器,对其数值分析结果揭示:当位于射流中剖面同侧的二次射流全部作动时,主流可以取得最大的偏转角度,且主流的偏转方向可以通过使不同的二次射流组合处于作动状态进行控制。      

不同进口条件对OGV/ 前置扩压器流场大涡 模拟结果的影响

为了研究不同进口条件对压气机出口导向叶片(Outlet Guide Vane,OGV)/前置扩压器流场性能的影响,分别采用合成涡方法与白噪声方法生成进口条件,研究不同进口条件对流场数值模拟结果的影响,并与文献给出的试验结果进行对比。结果表明:采用合成涡方法生成的进口速度场满足试验测得的时均速度分布与脉动均方根分布,且其生成的进口湍流结构可传播到下游较远处;合成涡方法能够较好地预测OGV内的流场,与白噪声方法相比,所预测的OGV吸力面分离区较小;在OGV出口截面、扩压器内和扩压器出口截面,采用2种方法得到的速度分布相差不大,与试验相比分布趋势相同。     

应用PIV研究二维后向台阶起动涡流动结构特性

采用适合于低速水流的白光气泡 PIV技术以及氢泡时间线流动显示技术在水槽中对绕二维后向台阶的起动涡瞬时速度场、涡量场进行了定量测量和流谱观察。得到了起动涡涡量的最大值、总环量、涡的尺度及涡心位置随时间变化的特性曲线 ;并观察到起动涡发展的 4个阶段