在静止扩压环形叶栅中捕获轴流压气机非定常耦合流型的实验研究

本文研究了在静止扩压环形叶栅中如何实现非定常自然流型(UNFT)向非定常耦合流型(UCFT)的转变，以及耦合流型对轴流压气机能带来多大的时均性能提高。在环形叶栅台上的大量实验研究表明：当以一定方式实现非定常耦合流型后，其分离旋涡缩小，流场中涡量减小，流动损失下降，时均性能有较大幅度提高。实验中将常规的气动感头测量和小惯性热线风速仪测量以及数字式PIV瞬态场测量技术相结合，来获得非定常两代流型的流动图案和时均性能及其相互比较。     

直升机旋翼/机身非定常气动干扰数值分析

基于非定常面元/时间步进全展自由尾迹建立了旋翼/机身非定常气动干扰分析方法。方法中耦合了非定常面元法和时间步进自由尾迹,以准确模拟旋翼非定常气动力、旋翼尾迹及桨叶对机身的非定常干扰效应。为模拟前飞状态下具有升力的机身,将机身离散为汇/偶极子面元,并采用涡线镜像法模拟旋翼尾迹靠近机身表面产生的加速效应。通过计算前飞状态的Maryland、ROBIN(Rotor Body INteraction)旋翼/机身干扰下的非定常压力分布,并与可得到的实验值、CFD计算结果对比,验证方法的准确性。随后分析前飞速度、旋翼与机身距离对旋翼/机身非定气动干扰的影响。计算结果表明机身头部和中部非定常压力主要受桨叶的通过性影响,而机身尾梁主要受尾迹/机身干扰影响,机身非定常气动力频率为桨叶片数的倍频。随前飞速度的增加,机身非定常压力幅值增加,尾梁压力幅值先增加后减小;增加旋翼与机身距离将减小机身和尾梁非定常压力幅值。     

基于动态混合网格的多体分离数值模拟方法

在以往发展的动态混合网格技术的基础上,耦合六自由度运动方程计算和非定常流场计算,建立了一种多体分离问题的非定常数值模拟方法。首先采用刚体运动的六自由度运动方程,建立刚体运动轨迹计算模块,然后以该模块为纽带将以往建立的非定常计算方法与改进的基于Delaunay背景网格插值方法和局部网格重构方法相结合的动态混合网格生成方法耦合起来,形成了精度较高、效率较快的多体分离数值模拟方法。利用该方法对典型多体分离问题进行了数值模拟,计算结果与风洞实验结果吻合较好。     

跨声速无粘流绕振动翼型流动的数值模拟

本文从Euler方程出发,利用LU-ADI格式对绕二元振动翼型跨声速流场进行了数值模拟。文章中给出了绕NACA0012和NACA 64A 010翼型不同状态的几组数值结果,在同某些实验和其它文献结果基本吻合的前提下,研究了振动对翼型表面压力、激波位置的非定常效应以及其它振动参数对流场非定常效应的影响。     

叶轮机械新流型探索雏议

从“流型”的角度, 对比飞行器空气动力学和叶轮机械气体动力学的发展历程, 认为叶轮机械内的流动仍在追求定常附体流型, 由于其流动及结构特性, 使得“有害”分离问题难以避免。本文提出应探索将外流领域内的定常脱体涡流型或下一代非定常脱体涡流型引入到叶轮机械的可能性与技术途径本文并综述了上述探索性研究的有关背景和现状, 以及所设想的基本途径和影响因素。      

基于分离涡模拟的起落架气动噪声研究

用基于SA湍流模式的非定常RANS、分离涡模拟(DES)和延迟分离涡模拟(DDES)分别对四轮基本起落架模型进行了数值模拟,并根据所得的非定常流场计算了表面声压级分布和声压谱。三种方法所用的时间成本大致相同,URANS略低于其它两种,而DDES由于在附着流动区更好地保持RANS特性,故时间成本略低于DES。计算显示,非定常RANS在附着流动区能够得到合理的结果,但不能准确刻画分离流动的流动形态。DES和DDES都能较好地刻画起落架绕流的定常和非定常特性,DDES的结果略好于DES。因为起落架流动属于大分离流动,所以DES也能够得到相对正确的结果。研究结果验证了分离涡模拟在起落架大分离非定常流动预测与噪声预测中的可行性,对减小起落架噪声的方法研究具有一定的意义。     

飞行器大迎角下俯仰静、动导数的数值计算

 采用一种适用于非定常流动计算的逐次超松弛（SOR）时间推进方法耦合求解非定常Navier-Stokes方程和强迫运动方程,数值模拟了NACA0012翼型、弹道外形(HBS)和有翼导弹等飞行器在不同初始迎角下做强迫俯仰振荡的黏性动态流场,给出动态气动力及力矩系数的时间历程;采用积分法获得飞行器的静导数和动导数,将计算结果与实验结果及半经验理论预测方法进行比较并分析飞行器在大迎角下动态特性的非线性特征;计算结果与实验值相吻合,表明本文采用的非定常数值模拟方法在大迎角下能够较准确地模拟飞行器的动态特性,具有工程应用价值。     

基于激励盘理论的螺旋桨滑流数值模拟研究

针对螺旋桨滑流流动复杂、非定常数值模拟难度大的问题,对激励盘模型进行深入研究,提高其对滑流模拟的精度。首先根据非定常数值模拟单独螺旋桨流动结果确定激励盘载荷。分析了考虑旋向载荷分布对滑流数值模拟的重要性。采用激励盘模型得到的定常滑流速度分布与螺旋桨非定常滑流的时间平均结果非常吻合,且采用激励盘理论能够明显地缩短模拟时间,在工程应用中具有很大的优势。     

座舱盖抛放过程中的非定常气动特性研究

通过对N-S方程的数值求解模拟了座舱盖在静止和旋转两种状态下的气动力分布,分析了来流马赫数、旋转角速度以及迎角对气动系数的影响,研究了座舱盖自由抛放过程的非定常效应,并最终建立了预测座舱盖自由抛放过程气动系数的气动系数图。计算表明:(1)非定常效应是影响抛放过程座舱盖气动力分布的主要因素;(2)外形左右对称的座舱盖,均有非定常气动系数平均值与无量纲角速度成线性关系的规律,而且变化斜率几乎与舱盖形状尺寸等无关。     

进口畸变下非定常数值模拟进口边界数学模型

进口畸变下压气机内存在的不同空间和时间尺度的非定常流动之间的相互作用是进口畸变对压气机性能影响的重要流动机理.介绍了一种用于周向进口总压畸变下单转子通道内非定常数值模拟的进口边界简化模型,并基于该模型进行了三维非定常数值模拟.数值模拟结果显示了进口周向总压畸变的大尺度总压脉动和压气机转子通道非定常分离流场的相互作用过程.该模型为周向进口畸变的三维非定常数值模拟提出了一种简化方法.      

旋翼干扰流场数值模拟中的新型Blend低速预处理方法

针对定常低速预处理在双时间步长内收敛性差、旋翼干扰流场数值模拟中非定常低速预处理自动退化为无预处理的问题,在非定常低速预处理方法的基础上,采用量级分析方法研究了对流项和黏性项量级,通过组合定常和非定常预处理特征值矩阵,建立了一种新型Blend低速预处理方法。通过圆柱非定常扰流、前飞直升机旋翼/机身干扰非定常流场和直升机/舰船干扰非定常流场的数值模拟,对各低速预处理方法的收敛性和预测精度进行了细致研究。结果表明,本文方法解决了在旋翼干扰流场数值模拟时非定常低速预处理自动退化为无预处理的问题,相比定常低速预处理大幅提高了收敛效率。     

欧拉方程计算非定常流动的一种数值方法

提出了利用欧拉方程计算非定常流动的一种数值方法。在利用中心差分的有限体积法和双时间推进法求解非定常欧拉方程过程中，通过在固定物面边界上满足完全的非定常边界条件，计算出非定常流场。算例结果表明，本文方法具有高效、简便的特点，而且能描述激波的非定常运动和变化情况。     

三维可压缩非定常流的壁面分离判据及其分离线附近的流动形态

通过引入依赖于密度的物面法向速度变换wr=-1/ρh1h2∫ozh1h2(e)ρ/(e)tdz,描述物理速度空间(ua,va,wa=u,v,w+wr)具有无滑移壁面条件的三维可压缩非定常连续方程可转换成变换速度空间(u,v,w)内具有无滑移条件定常连续方程.因此,采用定常壁面分离的分析方法和结论,再通过变换和研究wr的贡献,给出了三维可压缩非定常壁面分离的判则以及分离线附近的流动形态.研究指出,二维和三维情况下,都出现伴有壁外附着的壁面分离情况.数值模拟证实了理论和结论.     

跨声速轴流压气机近失速状态的间隙泄漏流流动特性

为了研究间隙泄漏流以及泄漏涡自身的非定常性对轴流压气机的旋转失速的影响.对跨声速轴流压气机NASA转子37进行全三维定常、非定常的数值模拟,对比了最大效率和近失速工况的实验和数值模拟结果,定常计算所获得的总性能与试验结果符合良好.对于非定常计算,详细分析了NASA转子37近失速工况下流场结果,揭示了NASA转子37在近失速工况点,间隙泄漏流存在较明显的非定常性,这种非定常性表现为间隙泄漏流激波干涉引起间隙泄漏涡的周期性破碎.      

叶型偏差对涡轮性能影响的非定常数值模拟研究

 叶型的加工和装配误差难以避免,这将对涡轮的气动性能产生一定的影响。采用定常和非定常数值模拟结合整机实验的方法,研究了叶型偏差对涡轮气动性能及内部非定常流动细节的影响,并比较了定常和非定常计算结果的差别。研究结果表明：在整个工作范围内,叶型偏差都会造成涡轮性能的明显下降,通道涡和泄漏流的发展、前缘吸力峰的强度、激波的形态和强度等流动结构以及叶排间的非定常相互作用都会随叶型的变化而显著改变,而非定常数值模拟对涡轮气动性能和流动结构受叶型偏差影响的捕捉更为准确。     

引入DMD方法研究有/无控气流分离的动态结构

为分析非定常流动控制技术抑制分离流的机理,对弯曲扩压通道的试验模型进行了数值模拟,针对扩压通道在无控和采用最佳射流频率状态下的计算结果引入了动力模态分解(DMD)技术进行分析。通过DMD技术能够将包含时空信息的扩压通道复杂流场进行分解,捕获流场包含的动力信息和对应的拟序流动结构。将无控和有控流场分解的结果进行对比分析后表明:采用有效激励措施时,和脱落涡频率一致的涡系对流场的影响更加突显,流场整体上表现得更加有序;非定常控制抑制了一部分涡的增长,使得各模态整体上更加稳定;而有控流场占主导地位的涡系结构相比无控流场较为有序,且对主流区未形成明显的直接影响。     

单级跨音速压气机三维非定常粘性流动计算

本文以数值求解雷诺平均非定常N-S方程方法完成了单级跨音速压气机动静叶相干三维非定常粘性流动计算,给出了高分辨率NND格式的时间推进有限差分解。对某跨音速单级压气机进行了定常和非定常流的详细数值计算,设计转速下定常粘性流场计算所得到的性能参数与实验结果吻合较好。在此基础上,对最高效率状态进行了动/静叶相互干扰的非定常流动计算,并对结果进行了分析。      

基于RANS/LES混合方法的分离流动模拟

飞行器在大迎角、快速俯仰机动时,流场中含有大尺度、非定常的涡结构,传统雷诺平均Navier-Stokes （RANS）模型不能准确模拟流场结构,根据国际上相关研究的发展趋势,需要采用混合RANS/大涡模拟（LES）模型来对复杂分离流动进行准确模拟。本文对基于分区混合与湍流尺度混合的双重RANS/LES混合计算模型进行发展与应用。通过典型简化模型的静、动态湍流大分离流动,测试和验证所采用的脱体涡模拟（DES）类方法,重点研究改进的延迟DES （IDDES）模型在动态问题应用中的正确性和有效性,并对所采用的数值模拟方法和相应的计算软件的可靠性、鲁棒性以及精度进行了考核验证。典型算例包括超声速圆柱底部流动、跨声速方腔流动、NACA0015机翼深失速分离涡模拟等。计算表明：发展的IDDES类混合计算模型可有效解决对数层不匹配的问题;对于定态非定常分离流动,DES、DDES、IDDES等模型计算结果差别不大,随着流动的非定常特性增强,IDDES模型的优势逐渐显现;对于动态非定常分离流动,则需要采用IDDES类模型。     

汽轮机末级三维非定常流动数值模拟

采用三维粘性非定常数值模拟方法,对某汽轮机设计状态末级非定常流动进行了模拟,并对设计状态下各因素引起末级动叶表面产生非定常力的大小和作用机理进行了分析。结果表明,本文中的汽轮机末级在设计条件下,上游静叶尾缘激波和位势场对下游转子的影响比尾迹作用的影响更大,静叶尾缘激波和位势场导致下游转子叶片表面非定常力大小呈周期性变化;此外上游静子尾迹以及动叶尾缘出口膨胀波等也是产生转子叶片表面非定常力的因素。