壁面展向匀速运动的槽道湍流的直接数值模拟

为了解展向存在匀速流动时湍流槽道流动的特性,应用直接数值模拟方法研究了壁面的展向运动对湍流槽道流动的影响.采用标准谱方法对固定质量流率条件下,上壁面展向运动速度分别等于零、流向流量平均速度以及两倍于流向流量平均速度等3种情况进行了直接数值模拟.结果表明,壁面的展向运动会导致表面摩擦力的增大,二阶湍流统计量诸如均方根速度脉动、雷诺剪切应力和均方根涡量脉动等都会随着壁面展向运动速度的增大而显著增长,同时,槽道截面上的均方根速度脉动也变得更均匀.      

微重力量值对气液两相流相分布及液相湍流统计量影响的数值研究

气液两相流在空间领域具有广阔的应用前景, 深入理解微重力量值(微重力大小)对相分布和液相湍流的影响十分必要. 采用欧拉elax-elax拉格朗日双向耦合模型深入研究了不同微重力量值对相分布和液相湍流的影响. 液相速度场通过直接数值模拟求解, 气泡的运动轨迹由牛顿运动方程跟踪. 研究表明, 气泡分布和液相湍流与微重力量值均具有直接联系. 在低微重力量值下, 气泡近似均匀分布在槽道内, 且对液相湍流统计量几乎没有影响; 然而当微重力量值较高时, 大量气泡聚集在壁面附近, 液相湍流由于气泡的注入受到极大调制.      

端柱面组合密封气膜稳态特性数值仿真分析

针对航空发动机等高速流体机械中流体动密封部位的工况特点,提出一种新型的端柱面组合气膜密封形式.根据气膜密封的结构特点,建立了端柱面螺旋槽组合式密封气膜的数学分析模型.采用CFD分析软件Fluent对气膜三维流场模型以湍流流动形式进行数值模拟计算.研究了密封气膜几何结构参数变化对密封稳态特性（端面承载力F_d 、柱面承载力F_c 、泄漏量Q、摩擦转矩M）的影响.几何结构参数包括端面平均膜厚、端面槽深比、柱面平均膜厚、柱面偏心率、柱面槽深比、柱面槽宽比.      

流向磁场作用下二维磁流体槽道湍流直接模拟

对低磁雷诺数近似下流向磁场作用的二维磁流体槽道湍流进行了直接数值模拟(DNS,Direct Numerical Simulation),给出了 Re=10 000时不同磁相互作用数下二维磁流体槽道湍流的近壁速度分布、湍流脉动速度均方根、雷诺应力等统计量的变化,并与中性流体二维槽道湍流进行了比较.结果表明:流向磁场作用会导致对数区上移;雷诺应力最大值随磁相互作用数呈线性变化;随着磁相互作用数的增大,下壁面平均涡量的时间演化由拟周期性向周期性转变,且脉动周期逐渐增大并当流动层流化后下壁面平均涡量成为常值.      

近壁区湍能耗散率输运方程的理论模型

提出在湍流边界层的近壁区采用三维不稳定波来描述湍流扰动速度,然后根据理论模型对湍能耗散率的生成、粘性破坏以及扩散等特性进行了系统定量的计算和分析.理论计算结果与直接数值模拟结果符合较好,特别是在y<10的区域中明显优于传统湍流模型的计算结果.这不仅在理论上有益于对湍流物理机制的了解,而且可能为湍流模型计算开辟一条新的途径.     

用GAO-YONG湍流模式对分离流动的数值模拟

采用GAO-YONG湍流模式对二维管道内三角形钝体绕流问题以及槽道内后台阶分离流动进行了数值模拟.求解程序基于开源数值计算平台OpenFOAM (Open Field Operation and Manipulation),数值模拟结果很好地预测了钝体绕流问题中漩涡脱落的尾流的流动趋势以及后台阶流动中分离再附的流动结构,同时分析了速度分布以及摩擦系数等参数并与实验值进行了对比,结果符合很好.这表明GAO-YONG湍流模式对大分离流动有较好的预测能力,对工程实践具有指导性作用.     

用GAO-YONG湍流模式对分离流动的数值模拟

采用GAO-YONG湍流模式对二维管道内三角形钝体绕流问题以及槽道内后台阶分离流动进行了数值模拟.求解程序基于开源数值计算平台OpenFOAM (Open Field Operation and Manipulation),数值模拟结果很好地预测了钝体绕流问题中漩涡脱落的尾流的流动趋势以及后台阶流动中分离再附的流动结构,同时分析了速度分布以及摩擦系数等参数并与实验值进行了对比,结果符合很好.这表明GAO-YONG湍流模式对大分离流动有较好的预测能力,对工程实践具有指导性作用.     

应用修正湍流模型对旋转通道内换热的模拟

针对旋转光滑矩形通道分别应用针对旋转状态修正的 k-ε、标准k-ω以及提出的针对旋转状态修正的k-ω湍流模型进行流动和换热的数值模拟,通过与实验结论的对比,讨论了采用不同湍流模型对计算结果的影响.计算工况为旋转数Ro=0.24,流体进口雷诺数Re=25000.计算结果表明:采用所提出的针对旋转状态修正的k-ω湍流模型的计算结果要比采用针对旋转状态修正的k-ε以及标准k-ω湍流模型的计算结果更接近实验结论.     

二维微槽道电渗流的增强混合

微尺度流动中的混合一直是微流系统研制中的一个大问题,对壁面zeta势非均匀分布情况下的二维微槽道电渗流进行了数值分析,重点分析了这种情况下的微流动混合增强.基于电渗流的Helmholtz-Smoluchowski滑移速度模型,应用SIMPLE方法计算了4种滑移速度在整个壁面上的分布不均匀一致时的微槽道电渗流动的流场及被动标量在其中的混合.计算结果表明,滑移速度的非均匀一致导致微槽道电渗流中的流线发生剧烈的转折,从而使得微槽道流动的混合效率有很大的提高.流场中流线的转折越剧烈、转折的次数越多,微槽道流动混合效率的提高就越大.     

基于共振三波的大涡模拟壁面模型

大涡模拟的一个重要发展方向是处理工程流动问题,但高雷诺数下近壁湍流边界层对计算网格的过高要求成为制约其应用的主要瓶颈;要减少计算量,构建合理且能够正确反映近壁湍流动力学特性的壁面模型是一条有效途径.在此背景下,构造了一种计算中便于实现的基于共振三波的壁面模型.通过对充分发展槽道湍流算例的研究,完成了对该模型构造方法、方案的论证和验证计算.结果表明使用壁面模型后,在减少计算量的同时,定性和定量的结果均较为合理,能够得到重要的瞬态特征结构和正确的一阶、二阶统计量,从而部分验证了该方法的有效性.     

沟槽面湍流减阻研究综述

对近20年来沟槽面湍流边界层特性、湍流拟序结构、湍流减阻及其机理的研究进展进行了综述.内容涉及沟槽面平板、旋成体、机翼等在压、跨、超音速流动情况下的实验研究成果;压力梯度、攻角、侧滑角等的影响;湍流猝发特性、紊动特性、近壁区带条结构的特征及减阻机理等方面的工作.为更有效地减少表面摩阻,必须深入开展对沟槽面湍流边界层特性、湍流拟序结构及湍流减阻机理等方面的研究.      

几何和流动参数对上翘后体阻力系数的影响

利用N-S方程对9个上翘后体模型进行了气动力计算,主要研究后体几何参数和流动参数对上翘后体阻力系数的影响.研究结果表明, 后体的压差阻力系数分别随上翘角、收缩比的增加及迎角的减小而明显增加;后体摩擦阻力系数分别随后体的长细比的增加和雷诺数的减小而增加;后体越扁平,其压差阻力系数越大;在跨音速时,波阻系数也与上翘角有关,上翘角增加会导致波阻系数进一步加大.     

基于多岛遗传算法的二维翼型吸气减阻优化

针对吸气减阻技术吸气质量较大的缺点,采用多岛遗传算法与计算流体力学(CFD)相结合的方法,对吸气减阻技术进行了优化.数值计算采用E387二维翼型,将吸气开口位置和吸气质量通量作为优化参数,通过寻找上翼面最佳的吸气开口位置,使翼型阻力和吸气质量最小.结果表明:多岛遗传算法能够有效地优化吸气开口位置和吸气质量;当吸气开口位于最佳位置时,吸气质量达到了寻优区间下限,而翼型减阻效果可达8.3%;吸气没有延缓流动转捩的发生,翼型阻力减小的原因主要是由于压差阻力的减小.      

底凹结构减阻效应数值分析

为研究弹丸底凹结构的减阻机理,使用三维定常CFD方法对M910弹丸的流场特性进行了数值模拟。给出了零升阻力系数随马赫数的变化规律,所得结果与实验数据符合很好。在此基础上,为M910弹丸引入底凹结构并进行数值模拟。对比了不同弹底结构的底部流场特性,对底凹结构减阻效应的产生机理进行了分析。结果表明:亚声速下,底凹结构在底凹腔体内引入了高压“死水区”,并以“屈从”的流体边界代替了原固体底面,从而改变了尾部涡街的形成位置、形状和强度,最终增大底部压力,减小弹丸阻力;跨声速下,由于尾部涡街远离弹丸底面,固体底面与流体边界面的作用相同,使得底凹不再具有减阻效果;超声速下,底凹结构的减阻机理与底排弹丸减阻机理类似,即底凹结构中的流体为弹丸底部回流区添加质量从而达到减阻作用。      

超声速进气道可压及不可压流动数值模拟

采用Chio-Merkle预处理矩阵对可压NS(Navier-Stokes)方程组进行时间预处理,分析了预处理方法的物理和数学背景.用有限体积方法,结合LU-SGS(Lower-Upper-Symmetric-Gauss-Seidel)隐式时间积分和AUSM+(P)(Precondition Advection Upstream Splitting Method)格式、中心差分2种空间离散格式,求解该预处理NS方程组.通过圆弧凸包流动和二维方腔顶盖驱动流动的数值试验表明,该方法克服了传统时间迭代方法模拟低速流动时的刚性问题,加速了收敛过程,可以同时有效地模拟可压及不可压流动.用此方法与块结构化网格技术相结合,进行了在不同飞行马赫数、攻角和出口反压条件下,混压式轴对称超声速进气道三维流场数值模拟.计算结果表明:该方法能准确地捕获复杂激波系,清楚地揭示进气道三维流动现象;进气道性能随工作条件的变化规律,与理论分析一致.      

仿生减阻表面的进展与挑战

减阻表面由于其在海洋船舶、管道运输、飞机和国防军事等各个领域展现出的巨大潜力和广泛的应用前景而受到越来越多的关注。自然界中的动植物进化出了许多独特的形态或表面结构，均展现出了卓越的减阻性能。对近年来受海豚表面启发的顺服表面减阻、鲨鱼皮表面启发的微结构减阻及受荷叶表面启发的超疏水减阻等进行了总结，并介绍了这些技术的研究进展、减阻机理和挑战，对目前仿生减阻表面的研究进行了总结与展望。      

边界层减阻研究综述

阐述了目前边界层减阻研究中的一些主要方法、手段及其原理,介绍了国内外对边界层减阻研究和应用的状况,指出了其中存在的一些问题,分析了今后减阻研究和应用的趋势。文中也包括了作者的一些观点和看法。     

临近空间表面波等离子体减阻性能分析

为研究临近空间表面波等离子体减阻效果,基于流体宏观模型的基本特征,分析了表面波等离子体流动控制机理的基础,以飞艇为模型,对其在 0°攻角情况下的流场进行仿真计算,比较了不同激励器控制方案的减阻效果,研究了飞艇尾部区域的等离子体流动控制效果. 结果表明,表面波等离子体具有增加飞艇升力、减小飞艇阻力的效果;单侧控制方案最大减阻效果达7%左右,对称控制方案减阻效果明显优于单侧控制方案,最大减阻效果可达32%左右;表面波等离子体对飞艇尾部的流动分离具有很好的消除抑制作用.      

尾翼对民机后体流动特性的影响

在1.5 m低速风洞中采用表面油流、烟线/激光片光显示和部件测力等实验方法对2种民机布局后体的绕流流型和气动特性进行了研究.实验的风速为40 m/s,雷诺数为1.191×105,迎角范围为-10°~20°.实验结果表明:随着迎角由负到正的变化,2种布局模型后体绕流都经历了上分离—无分离—下分离流型的变化,分离流型属开式分离并由尾部起始向机身发展;小迎角下2种布局模型后体阻力的差别主要由尾翼的摩擦阻力引起,而较大迎角下尾翼绕流分离引起的压差阻力起主要作用;垂尾分离随迎角增加而不断减弱,从而导致正常式布局后体阻力随迎角变化表现为非对称.      

机翼后缘连续变弯度对客机气动特性影响

后缘连续变弯度机翼在提高民用客机气动特性方面有较大的潜力,近年来被广泛关注。基于建立的全局优化设计系统,研究了机翼后缘连续变弯度对宽体客机翼身组合体气动特性的影响。首先,采用自由型面变形(FFD)技术建立了后缘连续变弯度的参数化方法。然后,采用RANS方程作为流场评估方法,针对翼身组合体构型设计点附近升力系数开展了机翼后缘连续变弯度气动减阻优化设计。最后,探索了仅外翼段后缘连续变弯度和内外翼后缘均连续变弯度优化设计结果的异同。优化结果表明,升力系数小于设计升力系数时,在只考虑外翼段后缘连续变弯度的设计中,不易实现激波阻力和诱导阻力同时降低,考虑内翼段后缘连续变弯度后,减阻量较前者更为明显;升力系数大于设计升力系数时,外翼段和内外翼的后缘偏转均可实现诱导阻力和激波阻力的同时降低,且减阻量相差不大。