薄翼失速翼型前缘分离泡对失速特性的影响

DES方法结合了RANS(Reynolds-averaged Navier-Stokes)和LES(Large Eddy Simulation approaches)的优点。在近壁面它体现为RANS模型的特点而在远离物面处又起到LES的亚格子模型的特性。论文应用DES(Detached-EddySimulation)方法讨论了影响薄翼失速的分离泡对翼型的升力特性影响。     

基于当地流动结构的RANS/LES混合模型

针对传统雷诺平均Navier-Stokes/大涡模拟(RANS/LES)混合模型固有的RANS和LES区域衔接难题,借鉴Vreman提出的亚格子(SGS)模型,提取其中描述流动状态的判断函数,利用不可压缩槽道湍流标定模型系数,构造了一个新的RANS/LES混合模型。评估算例包括不可压缩槽道湍流、稳态超声速平板边界层、NACA4412翼型绕流和亚临界雷诺数圆柱绕流。该模型作为壁面模化LES,解决了传统RANS/LES方法中速度型偏离对数律(LLM)的缺陷,且拓宽了原Vreman模型在粗糙网格上的适用性,具有较高精度。该模型同时解决了模化应力折损问题,且对大尺度分离流动的模拟精度相较传统的RANS/LES混合模型如脱体涡模拟(DES)有进一步提升。     

一种改进的类DES湍流模拟方法

构造了一种基于一方程S-A(Spalart-Allmaras)模型和一方程Yoshizawa亚格子模型的混合RANS/LES(Reynolds-averaged Navier-Stokes/large eddy simulation)湍流模拟方法.在涡黏假设的基础上,将Yoshizawa亚格子湍动能方程转化为等效的亚格子湍流涡黏性输运方程,并采用混合函数将其与S-A模型方程进行混合,从而改进了DES(detached eddy simulation)模型的亚格子行为,同时克服了其依靠网格控制模型转换的缺点.模拟了超声速的带斜坡凹腔流动,并与相同网格下的LES及DES结果进行了比较,结果表明该混合RANS/LES方法在远离壁面的自由剪切流区域与LES行为一致,而在附着边界层区域表现优于LES和DES方法.      

翼型结冰状态复杂分离流动数值模拟综述

结冰触发的复杂分离流动将导致翼型气动性能特别是失速特性全面恶化。结冰状态气动特性的准确预测和流动机理的深入剖析依赖于分离流场结构的精确求解。随着计算流体力学特别是湍流模拟方法的不断完善，数值模拟能够更为清晰和完备地反映非定常分离流场的细节特征及物理本质、提供更加翔实和丰富的气动力数据。从雷诺平均（RANS）、大涡模拟（LES）和RANS/LES这3类典型湍流模拟方法的应用层面出发，综合评述了近年来数值模拟研究在翼型结冰状态失速特性预测与分离特征描述等方面取得的主要进展，并从高精度冰形构造、新型湍流模拟方法、深层次非定常特性、实时耦合分析等方面对现阶段研究发展的相关趋势进行总结和展望。     

多段翼型大迎角分离流动的DELAYED RANS/LES混合算法

多段翼型的大迎角绕流发生大范围附面层分离,具有明显的三维与非定常流动特性。RANS/LES混合算法继承了LES对流动分离区大尺度漩涡准确模拟的优点,避免了纯LES算法需求网格量巨大与亚格子模型壁面函数不成熟等问题,对分离流动的模拟效果优于RANS算法。以S-A湍流模型与Smagorinsky亚格子模型为基础,借鉴DDES的附面层延迟控制思想,构造了可用于对接网格、重叠网格的DELAYED RANS/LES混合算法。研究了GA（W）-1多段翼型的大迎角分离流动及其气动特性。     

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

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

宽速域RANS-LES混合方法的发展及应用

传统的雷诺平均方法(RANS)已经不能满足大范围分离、激波振荡、压力脉动、动载荷等极端工况下的流动预测需求;大涡模拟(LES)、直接数值模拟(DNS)等方法资源耗费多、效率低,离工程湍流问题仍较为遥远。RANSLES混合方法结合了RANS高效率和LES高精度的特点,近期有望大规模应用到工程湍流问题中。首先对现有的RANS-LES混合方法进行了归类,对各自的构造思想、特点进行了分析。然后报告了脱体涡模拟(DES)类方法的发展历程和现状,讨论了使用DES类方法计算分离流动时,对流项离散格式对分离特性、小尺度结构及频谱特性等的影响,并构造了自适应耗散函数。最后介绍了近年来国内外RANS-LES混合方法在宽马赫数范围(马赫数从0.1到20)内的机理研究和工程应用。现有的以DES类方法为代表的RANS-LES混合方法能够较为精细地模拟非定常大分离流动中的复杂现象,但在计算效率等方面还有较大的改进空间;植入式DES方法在模拟全机带部件流动上具有较高的效率和模拟精度,是重要的发展方向。RANS-LES混合方法在动态失速、燃烧、气动弹性、气动噪声、气动光学等与非定常流动密切相关的方面也有广阔的应用前景。     

旋翼翼型非定常动态失速特性的CFD模拟及参数分析

构建了一套基于运动嵌套网格技术和可压缩RANS方程的旋翼翼型非定常流动特性模拟的高效、高精度的CFD方法。首先,发展了基于Poisson方程求解的围绕翼型的粘性贴体正交网格生成方法,并提出了基于最小距离法(MDM)改进策略的运动嵌套网格生成方法,克服了弹簧法可能导致网格畸变的不足;其次,为准确模拟由湍流分离和气流再附引起的气动力的迟滞效应,基于RANS方程、双时间方法和高阶插值格式,建立了旋翼翼型非定常气动特性分析的高精度数值方法,并采用能够较好捕捉气流分离现象的S-A湍流模型;再次,针对旋翼后行桨叶动态失速时桨叶剖面来流速度较低、迎角较大的特点,为解决低来流速度时L-B半经验模型在旋翼翼型非定常动态失速计算中的局限性,并克服可压缩方程对低速流场计算收敛困难和精度低的问题,建立了基于Pletcher-Chen低速预处理方法、FAS多重网格法和隐式LU-SGS方法相结合的高效数值方法。应用发展的方法,分别针对NACA0012、SC1095旋翼翼型静态和轻度、深度动态失速进行计算,精确捕捉了气动力迟滞效应以及翼型前缘脱体涡的产生、对流和脱落过程,验证了本文方法的有效性;最后,着重针对NACA0012动态失速状态,开展了振荡参数对旋翼翼型非定常动态失速特性影响的分析,研究结果表明翼型迎角平均值、振幅及减缩频率的变化均能引起迟滞效应的改变并使得气动力峰值发生有规律的前、后移现象等。     

混合型三阶格式及关于两种DES算法的比较计算研究

为开展RANS/LES混合模拟,在传统三阶迎风偏置格式的基础上,提出一类混合型的三阶计算格式,格式在中高波数范围具有可调的耗散水平.在此基础上,文中通过计算比较研究两种RANS/LES混合算法:一种是基于Spalart-Allmaras一方程模型的DES模型(DES-SA),另一种是基于混合长度模型的DES模型(DES-ML),使用 发展的三阶格式对圆柱绕流进行模拟,并将得到的结果与文献比较,对计算格式、DES-ML算法进行了初步探讨.     

现代超临界翼型设计及其风洞试验

开展了现代超临界翼型的设计研究,对现役飞机的压力分布形态进行了分析,针对现役飞机在巡航状态和阻力发散点的压力分布进行对比,提取了现役飞机超临界剖面设计的要点。采用类函数/型函数变换(CST)参数化方法、基于二阶震荡及自然选择的随机权重混合粒子群算法(RwSecSelPSO)、雷诺平均Navier-Stokes(RANS)方程、Kriging代理模型结合定期望值型的目标函数建立了优化设计系统。针对提高阻力发散马赫数和降低巡航低头力矩的设计指标,利用优化设计系统通过调整目标期望值设计了一系列满足设计指标但阻力发散马赫数不同的超临界翼型,并选择了其中具有典型特性的翼型进行了对比分析,验证了提高阻力发散马赫数和低速失速特性的设计方法,指出了在阻力发散点形成平顶形压力分布的超临界翼型具有较好的综合性能。对设计的超临界翼型进行了高、低速风洞试验验证,试验结果表明:设计结果达到了设计指标要求,提出的低速改进方案有效,层流对超临界翼型失速特性影响较大。     

一种减少失速保护虚警率的混合滤波方法

针对飞机失速保护设计过程中易出现虚警的情况,提出一种基于FIR和传统滤波器相结合的混合滤波方法。结合飞机本体的运动特征,对参与失速保护逻辑运算的两个重要参数,即侧向加速度和迎角信号进行频率分析,基于信号的特征设计混合滤波器,然后采用Simulink滤波工具箱设计滤波器对信号进行滤波处理。仿真结果表明该方法能够减少虚警率的发生,特别是大背景噪声下的虚警率能够得到较好抑制,验证了滤波器设计正确性,从而提高失速保护功能对适航条款符合性。该滤波器设计方法,不但可以有效解决本文提到的问题,对于其他涉及信号处理的系统设计同样具有借鉴意义。     

带执行装置的压气机系统建模及仿真

针对压气机主动稳定控制对模型的要求,以Moore-Greitzer模型为基础,首先将模型进行空间离散化,用平均分布在压气机周向环面的离散位置点来描述流量扰动,建立了反映压气机周向动态特性的多维状态空间模型;其次以喷气装置作为主动稳定控制系统的执行装置,通过分析加入喷气装置后对流动区域产生的质量和动量影响,同时考虑气流从喷气装置出口到压气机起始平面存在的时间滞后,建立了带执行装置的压气机系统动态模型.以某压气机为例,对压气机失速行为及喷气状态下对失速的抑制作用进行仿真,结果表明:所建立的模型可以反映压气机旋转失速动态过程,通过控制喷气流量能拓展压气机的稳定工作范围,有效抑制压气机失速,模型可用于压气机主动稳定控制系统分析与设计.      

可变转速工况的压气机分布式动态模型与仿真

针对目前压气机动态模型用于主动稳定控制时,存在既不反映压气机变转速工况又不能模拟压气机内部周向位置失速信号发展过程的不足,在一维状态变量描述的压气机失稳动态模型基础上,以喷气装置作为执行装置,通过傅里叶变换,采用多维状态变量描述压气机动态过程,建立了含喷气执行装置的压气机系统变转速工况下分布式动态模型。仿真结果表明,所建模型可实现定转速和变转速工况下的压气机不稳定行为动态仿真,采用该模型进行压气机主动稳定控制,可使压气机失速模态幅值由0.28减小为1×10~(-5),能有效抑制压气机失速的发生,增强了模型的适用性。     

飞机结冰后非线性气动力建模及动态响应特性研究

当前结冰后的动力学仿真主要基于线性气动力模型及结冰影响模型进行,对于结冰后非线性气动力模型的研究较少,然而这对于飞机结冰后在大迎角等非线性区域的动力学分析极为重要.研究了飞机结冰后的非线性气动力模型,将结冰后线性气动力模型的估算方法扩展到非线性气动力模型上,并在结冰后大迎角区域风洞试验数据的基础之上,对失速及过失速区的气动参数变化进行修正.研究了飞机在不同升降舵脉冲信号作用下的动态响应.仿真结果表明,计算结果能够较好地反映飞机在结冰状态下的动力学特性,提出的非线性结冰影响模型能够为结冰条件下研究飞机失速及过失速区域的动力学特性提供理论支撑.     

压气机旋转失速的一个气动声学模型

 根据广义Lighthill方程建立了压气机旋转失速的气动声学模型。提出了一种应用声学方法获取压气机旋转失速特征参数(旋转失速传播速度、数目)的新途径并结合试验结果进行了具体的分析。此外,在简化假设下,用数值法分析了旋转失速特征参数对声音幅值变化的影响。     

小型电动无人机深失速回收仿真研究

对几种典型的小型无人机回收方式进行了比较,针对小型电动无人机翼载荷较小的特点,深失速回收方式具有明显优势,并给出了具体的深失速回收方案.通过分析大迎角情况下无人机的动力学模型,推导了在深失速情况下无人机的状态方程,建立了无人机深失速回收仿真模型,并以某型近程电动无人机为对象,在Matlab/Simulink中对深失速模...     

发动机轴流压缩系统逐级模拟技术

基于轴流压气机逐级过失速特性,建立了一个能够显示压缩系统喘振和旋转失速等过失速行为的一维非定常逐级模拟技术,发展了分析轴流压缩系统过失速响应的动态滞后方法。并将该模型应用于国内某型号的高转速、双转子压缩系统的过失速行为的分析,讨论了时间滞止的常数、燃烧室对压缩系统失速行为的影响,得到了合理的定性结果。     

新一代战斗机非定常流动数值研究综述

新一代战斗机强调超机动能力和强隐身性，其中大攻角下的静态失速、动态失速及内埋弹仓绕流是与高机动和强隐身密切相关的、极具挑战的几类典型的非定常流动，它们对数值方法提出了极高的要求。为了高精度地仿真流场、清楚地揭示流动机理，有效地控制非定常流动，非常有必要发展高精度且高效率的RANS-LES混合方法体系，包含RANS-LES混合方法本身、与RANS-LES混合方法匹配的高精度自适应耗散格式、基准湍流模式、高质量计算网格、高精度时间推进方法、非定常量的统计方法等，具有极强的紧迫性。提出、发展、验证并应用该类方法数值仿真新一代战斗机（包括单独部件、组合体、甚至全机）的非定常流动，数值预测结果与风洞实验数据吻合良好；此类方法可为新型战斗机设计提供理论依据和分析手段。     

Study on effects of thickness on airfoil-stall at low Reynolds numbers by cusp-catastrophic model based on GA(W)-1 airfoil

The phenomena of an airfoil stall present the behaviors of catastrophe and hysteresis at low Reynolds numbers. Numerical simulation results of two-dimensional airfoil GA(W)-1 show that the width of the hysteresis loop of airfoil stall will gradually decrease and even disappear with the decrease of thickness ratio. These nonlinear characteristics are in accordance with the topological features of the cusp catastrophic model. According to the topological invariant principle, a novel topological mapping method is developed to establish the mapping relationship between cusp catastrophic model and stall characteristics of the airfoil, then the effect of thickness ratio on airfoil stall is successfully described quantitatively by cusp catastrophic model. Further, based on the established topological mapping relationship, combined with the mean flow field of the airfoil stall, potential function approach of cusp catastrophic model is first introduced to interpret the catastrophe and hysteresis of the airfoil stall, and it is found that as the thickness ratio decreases, the system's maximal potential energy gradually disappears, and the short separation bubble at the leading edge of the airfoil changes to long separation bubble, so the airfoil stall changes from a bistable system to a monostable system.