轴承腔油气两相流数值模拟的试验验证

为探索一种适于模拟轴承腔油气两相流流动特性的数值计算模型,建立了一套简化结构的轴承腔试验装置,通过试验研究获得了不同工况对轴承腔通风口出气比例和回油口滑油温升的影响规律.应用Fluent软件,采用Mixture多相流模型和Eulerian多相流模型分别搭配标准k-ε湍流模型和RNG(renormalization group) k-ε湍流模型对选定试验工况进行数值模拟.对比结果表明:当转速为2200r/min时,使用Mixture多相流模型搭配RNG k-ε湍流模型计算的滑油温升与试验值的相对误差仅为2.72%,体现了该计算模型在高转速时的精确性,并在此基础上获得了油气两相流速度场、压力场和温度场的分布特征.      

二维超音速湍流混合特性研究

着重讨论顺流喷射情况下的超音速混合的流场特征及影响因素。简要介绍了顺流喷射混合的研究现状。同时，探讨了不同的来流参数对流场特征的影响及混合机理。计算采用了一种二阶ＴＶＤ格式有限体积法及Ｂ－Ｌ湍流模型求解二维湍流Ｎ－Ｓ方程。结果表明：所采用的计算方法能成功地抑制非物理振荡，适用于二维超音速内流的湍流混合特性研究。     

等熵侧压诱导的超声速平板边界层二次流研究

为了揭示等熵侧压诱导的超声速平板边界层二次流的形成机理,基于建立的简化流场模型,结合小扰动分析和数值计算,研究了来流边界层马赫数不均匀分布因素对二次流形成的贡献。研究表明,侧压诱导的二次流是流场的平均行为,和瞬时的扰动与涡的相互作用无关。压缩波在主流和边界层中的传播特征不同,导致在侧壁压缩作用下主流和边界层流体压缩不均匀,从而产生了边界层流体相对于主流的横向运动;密度梯度和压力梯度矢量方向不一致,即流动的斜压性是流向涡形成和发展的一个重要因素。     

凝结相变对低温风洞雷诺数试验的影响

为了研究凝结相变对低温风洞雷诺数试验能力的影响,基于Fluent软件建立了气-液两相凝结流动模型,相变模型采用考虑非等温效应修正的经典成核理论和Gyarmathy液滴生长理论;针对不同来流压力和不同试验雷诺数工况条件,对氮气绕流NACA 0012翼型进行了数值模拟。模拟结果表明:随着来流温度的降低,翼型附近区域气体膨胀越过气态饱和线并达到过冷状态,进一步降低来流温度则会在宏观层面上观测到凝结相变对当地流场的改变,与无凝结相变的流场相比,释放的潜热加热气流导致马赫数降低及压力系数的改变;在确保不破坏翼型气动性能试验的前提下,充分利用气体过冷区域来降低来流压力以此减少驱动功率和液氮喷入量是可行的,或者保持来流压力不变提升低温风洞的试验雷诺数。     

高空长航时无人机螺旋桨滑流效应影响研究

推导和分析了以多参考系模型作为螺旋桨计算模型的控制方程。应用数值模拟方法开展高空长航时无人机滑流效应影响的三维数值模拟研究。研究发现,多参考系模型的流动现象能够符合真实螺旋桨的前后流动特征,并且可以较好地模拟螺旋桨滑流对飞机气动性能的干扰。螺旋桨滑流效应使V尾表面流线发生偏转和收缩加速,V尾表面的压力分布明显改变。起飞状态螺旋桨滑流效应对全机气动特性影响最强,爬升状态影响减弱,巡航状态影响最小。滑流效应影响随着推力增加而增大,相同推力不同桨距条件下滑流效应影响基本相同。起飞状态无人机尾部受到螺旋桨大推力滑流效应影响压差阻力急剧增加,导致全机气动性能下降。     

涡轮叶片来流角扰动不确定性量化分析

为了高效精确量化来流条件变化对涡轮叶片性能的影响,发展了一种基于自适应抽样技术的非嵌入式多项式混沌(NIPC)模型,并通过函数试验对自适应NIPC模型进行了验证。在不同尺度来流角变化下,用该模型量化涡轮叶栅能量损失系数的不确定性变化,同时与直接蒙特卡洛模拟和灵敏度分析进行对比。通过流场统计分析揭示来流角变化对气动不确定性的影响机理。结果表明:自适应NIPC在不同尺度来流角扰动下均能快速准确量化能量损失系数变化;能量损失系数不确定性变化对来流角扰动是非线性依赖关系;叶片通道内的激波对来流角变化最为敏感,是引起能量损失系数不确定性变化的主要因素。     

来流条件对SFRJ燃速及自持燃烧性能的影响

为了研究来流条件对固体燃料燃速以及自持燃烧性能的影响规律,采用3阶单调迎风格式(MUSCL)重构方法,AUSMPW⁺通量分裂格式,k-w剪切应力输运(SST)湍流模型,7组分3反应有限速率化学反应模型以及2阶矩湍流燃烧模型,编制了二维轴对称湍流燃烧仿真程序。仿真结果显示:来流质量流率和总温的增加会提高燃速,温度和质量通量的影响指数大约为0.95~1.00和0.67;来流条件影响燃速的主要机理是壁面附近温度和湍流黏性的变化,其中湍流黏性是主要影响因素;随着质量流率的增加,特征速度逐渐降低,当入口马赫数接近1时,特征速度大约为850 m/s,固体燃料将无法自持燃烧。     

来流含水对航空发动机风扇/压气机特性的影响

结合大小水滴模型和由大量试验数据得到的损失和落后角模型,采用两相流理论与数值模拟的手段,在来流含水的状态下计算了某型发动机压缩系统的特性;分别研究了来流含水量、水滴尺寸等对压气机性能的影响。结果表明:进口气流中水的含量以及它存在形式均对压缩系统性能有不可忽视的影响,同时也验证了本文方法有效性和应用价值。      

侧流效应对发动机稳定性的影响

本文是在以二阶滞后模型研究周向稳态压力畸变对发动机稳定性影响基础上，增加侧流影响因素的，因侧流影响在本质上与压力畸变流影响相同，都是改变了压气机进口面的攻角分布，故使问题深入研究得以简化。     

格子气模型在机场出港大厅旅客疏散过程应用

利用偏向随机行走格子气模型中方向选择概率的计算方法和该模型中的串行更新规则,建立了偏向随机行走格子气行人流模型,将其应用到机场候机楼出港大厅,并对出港大厅内的出港旅客流的形成、平衡、缓解和消散四个阶段进行了模拟和研究,给出了模拟过程,并得到了出港旅客流的流动规律和运动性能参数.     

内转式进气道/冯·卡门乘波体一体化设计方法

应用特征线理论设计了内转式轴对称基准流场以及外压缩轴对称基准流场，利用激波交线、流线追踪方法等相关技术提出了一种头部进气式的高超声速飞行器内转式进气道/冯·卡门乘波体一体化设计方法，并对生成的一体化构型进行了数值模拟及分析，数值结果验证了该方法的正确性和有效性。该一体化设计方法基本保留了内转式进气道的优良特性，并以高升阻比乘波体为原型构建较高升阻比的一体化构型，从流场耦合的角度出发为减弱机体与进气道之间复杂的波系干扰，实现飞行器内外流的完全耦合进行了探索。     

弹用内乘波式进气道起动性能研究

 为弄清内乘波式进气道在低马赫数状态下的流动特征,分析影响内乘波式进气道起动能力的因素,研究与弹体匹配设计的内乘波式进气道的起动问题。首先基于一种有利于出口均匀性的基本流场,采用流线追踪技术,设计了来流马赫数为4.0且进出口形状适应弹体安装要求的双模块弹用内乘波式进气道;此后,采用计算流体力学(CFD)方法获得了低马赫数下进气道的三维波系结构和流动特征。研究表明,进气道溢流口位置是影响内乘波进气道起动能力的重要因素:在溢流口位置由两侧改至最下端后,起动马赫数由3.6下降为3.3;采用单模块方案,溢流口设置在下端后,起动马赫数下降为3.25。此外,设计内乘波式进气道基本流场也对起动性能有影响:设计出口马赫数不变,双模块方案下,入口气流偏转角每增大2°,起动马赫数约下降0.1;单模块方案下,提高入口气流偏转角最大可使起动马赫数下降为3.1;进气道内收缩比对起动能力的影响体现在入口气流偏转角不变时,进气道起动能力仅取决于内收缩比,设计出口马赫数每增加0.2,起动马赫数约减小0.2。研究所分析的各个弹用内乘波式进气道在设计条件下均可捕获99%的来流,在扩大了工作马赫数范围的同时,保持了高流量捕获性能和高总压恢复系数的优势。     

分形树状通道换热器内的流动换热特性

建立了分形树状通道换热器中层流流动与传热的三维稳态模型,采用流固耦合计算方法对入口水力直径为4mm的矩形截面树状通道内流动换热进行了数值模拟,重点研究了分叉效应对传热的强化机理和换热器受热面的温度分布。研究结果表明:分叉处形成的二次流能有效地强化换热;与传统的蛇形通道相比,分形树状通道换热器具有温度均匀性好、压降小的明显优势。在相同入口雷诺数时,分形树状通道换热器受热面的最大温差远小于蛇形通道换热器,另外,分形树状通道的层流流动压降较之蛇形流道可减小50%以上。同时,加工了分形树状通道换热器及蛇形通道换热器各一套,对数值模拟结果进行了实验验证。实验值与模拟值能较好地吻合,证明了所建流动换热三维数值模型正确可信。     

高性能流场并行粒子追踪数据管理系统

随着当下计算能力和存储性能的提升,流场数据产出的规模越来越庞大,针对流场数据的可视化应用对于硬件及软件算法的要求也随之提高。基于国家数值风洞（NNW）工程支持,主导设计并开发了高性能流场并行粒子追踪数据管理系统,帮助用户探索和分析大规模流场数据。该系统针对流场数据提供多种高效的数据管理方法,在超算集群上针对并行粒子追踪过程进行了数据预取优化与负载均衡优化。对于粒子追踪过程中产生的流线（或迹线）及进程工作记录数据,该系统支持用户在本地平台上进行性能诊断和分析。使用不同流场数据集开展的两个应用实例验证了该系统的有效性。     

基于离散等价方程的非结构网格有限差分法

激波装配法把解析关系式嵌入流场避免了间断引起的理论问题，使全场一致高精度的实现成为可能，有望解决超声速流动转捩研究中的感受性模拟难题。但装配复杂激波结构以后，被分割的流场空间常出现不规则的几何形状，这给常规的基于结构网格的有限差分法应用带来困难。基于离散等价方程理论，提出了一种新的基于非结构网格的有限差分法，在空间二维离散点附近仅用3条网格线就可以构造出一阶迎风格式。数值算例表明收敛过程对网格质量不敏感，解决了激波装配法模拟激波相交出现小夹角后使用结构网格进行计算存在的难题。根据这种方法的特点展望了未来的应用前景。     

流体运动中的非Fourier传热

 分析了运动流体中非Ｆｏｕｒｉｅｒ传热的数理基础。运用特征一次特征方法阐明基本方程组的影响域、决定域等概念及其物理内涵；强调热Ｍａｃｈ波的存在。以流动中固定热源问题证明基本关系正确；并给出当流动速度不小于热传播速度时,流汤中发生热激波。     

有翼飞行器大攻角分离涡模拟

基于势流方程,采用涡格法和三维离散涡相结合,模拟具有升力面分离涡的大攻角弹冀、弹体、尾翼组合体绕流。该方法适用于高亚音速、小攻角或大攻角有翼飞行器的压力分布、法向力和压力中心等气动特性的数值计算,并具有显示三维涡流场的能力。     

高精度ENO格式在射流数值模拟中的应用

本文分析了ENO格式的特点，并应用于全NS方程的迁移项和压力顶，模拟了轴对称和三维射流。首先对轴对称自由射流为算例比较了高阶TVD格式和ENO格式，证实了ENO格式捕获复杂激波波系的能力；其次计算了四喷管湍射流干扰流场，给出了相应的计算结果；最后计算了四股射流流入圆管内的流场，给出了物理量沿壁面和截面的分布。     

基于时间谱方法的振荡翼型和机翼非定常黏性绕流数值模拟

 传统的双时间方法在非定常计算中长时间的过渡迭代推进求解导致其计算效率相对较低,针对周期性非定常流动问题的流动特征,发展了一种基于离散傅里叶变换的高效时间谱方法,用于求解振荡翼型和机翼的非定常黏性绕流。在时空耦合的雷诺平均Navier-Stokes(RANS)方程的求解中,对流项的离散应用了Roe的通量差分格式,物理时间项的离散方法为时间谱方法,伪时间推进采用了隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式。考虑到湍流的时空耦合效应,时空耦合的Spalart-Allmaras一方程湍流模型的物理时间项同样采用时间谱方法进行离散。为了进一步提高计算效率,当地时间步长和多重网格技术等加速收敛的措施均被采用。算例对俯仰振荡NACA0012翼型和Lann机翼的周期性非定常流场进行了数值计算。结果表明:对于周期性非定常流场的数值模拟,相比于传统的双时间方法,用时间谱方法近似物理时间项,不仅能够提高流场的计算精度,而且更能够大幅度提高计算效率。     

Aeroelastic dynamic response of elastic aircraft with consideration of two-dimensional discrete gust excitation

Design loads generally require a one-dimensional discrete gust profile without consideration of the spanwise effect, and this profile cannot represent the true gust field exactly. For a high aspect ratio aircraft, two-dimensional gusts may cause critical load conditions, and approaches for calculating dynamic responses under two-dimensional discrete gust excitation are rarely presented. In this paper, a spanwise non-uniform vertical discrete gust field is established based on a one-dimensional ‘1-cos’ gust profile in reference to a DARPA proposal, while frequency and hybrid approaches to the dynamic response analysis of flexible aircraft under this two-dimensional gust excitation are presented. Solution techniques have been applied to a high aspect ratio aircraft to assess the different response characteristics with a comparison between one-dimensional and two-dimensional discrete gust field conditions. The results show that the two-dimensional discrete gust model produces a higher bending moment than that of the one-dimensional condition. Therefore, the critical load conditions that are derived from the two-dimensional discrete gust for high aspect ratio aircraft should be seriously considered. According to the analysis, an active control scheme to alleviate the bending loads caused by the two-dimensional gust is designed, and alleviation effects in different gust conditions are compared.