谐波平衡法在非定常流场中的应用

谐波平衡法（HBM）适用于模拟周期性非定常流场。求解这类问题时,只需要一个周期中几个等距时刻的流场解,即可重建整个周期流动的时间历程,计算效率较高;同时还采用谱算子精确求解流场变量对物理时间的导数项,保障了算法的计算精度。鉴于此,详细介绍了谐波平衡法的原理和实现过程,并给出了如何在现有计算程序下实现谐波平衡法的流程。取跨声速翼型和超声速钝锥两个周期性强迫俯仰振荡绕流为例,以双时间步方法为参考,对谐波平衡法的效率和精度作了详细的考核和分析。计算结果表明,谐波平衡法的内存消耗较大,谐波数N_H=2时约为双时间步方法的4倍,但谐波平衡法取得与双时间步方法相近结果的计算时间仅为后者的1/5。因此,谐波平衡法是一种有应用前景的工程实用方法。对于长周期问题,谐波平衡法的优势更加明显。     

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

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

旋翼翼型俯仰沉浮运动非定常气动特性实验研究

在南航非定常风洞内研制、开发了一套两自由度电控液压动态实验台,利用该设备,完成了翼型俯仰运动、沉浮运动以及俯仰沉浮耦合运动测压实验,并积分计算了翼型动态运动过程中的升力系数.翼型俯仰运动结果表明,翼型上的压力分布和升力系数与翼型振动频率和振动幅值有很大关系,随振动频率和振动幅值的增加,升力系数迟滞包线相应增加,而这一过程与翼型上前缘涡的脱落过程有关.翼型沉浮运动研究了翼型在不同迎角、不同风速下做沉浮运动时的压力分布和升力系数,结果表明在20m/s风速下,模型做沉浮运动时只对翼型前缘的压力分布产生影响,而在10m/s风速下,沉浮运动不仅使前缘压力分布产生变化,而且影响翼型下游的压力分布,这与翼型前缘形成的分离气泡脱落有关.另外给出了翼型做俯仰沉浮耦合运动时的压力分布初步结果.     

外挂俯仰连接刚度对机翼颤振的影响

本文使用v-g法、Niblett法和模型风洞试验3种方法,分析了机翼和外挂间俯仰连接刚度对机翼颤振的影响。并提出了“气动压力中心线”概念。用振型节线与压力中心线相对位置的变化,解释机翼/外挂颤振类型变化的机理。     

基于CFD方法的俯仰静、动导数数值计算

应用时空二阶精度的隐式迭代NND算法,数值模拟了尖锥、钝锥、弹道外形和飞船返回舱等典型再入飞行器作强迫俯仰振荡的粘性动态流场,给出动态气动力和力矩系数的时间历程,在此基础上,运用积分法和奇异分解线性最小二乘法辨识静、动态气动稳定性参数,并与试验结果及半经验理论预测方法进行了比较,表明本文方法具有准确度较高、效率较高及适应性广等优点.     

大迎角下俯仰振荡对细长体气动特性的影响

通过对头部无扰动和有扰动的两个细长体模型进行数值计算,分析了各个流场中涡系结构和气动力变化历程,并对比两者之间的差别.发现在高频率小振幅的俯仰振荡作用下,头部无扰动和有扰动的细长体的流场随时间的发展慢慢由非对称转变为对称.相对于头部有扰动的细长体,在第一个振荡周期内非定常俯仰运动对头部无扰动的细长体作用更明显.     

细长旋成体大迎角非定常运动数值研究

利用DES方法对处于非对称流态下的细长旋成体进行了俯仰和偏航振荡的数值模拟,观察细长体气动特性尤其是背风面非对称分离涡的变化。计算结果表明在大迎角非对称多涡系情况下,固定频率和振幅的非定常运动可以改变流场结构和气动力,对细长体背风面流态有巨大的影响。俯仰振荡对非对称分离涡有明显的控制作用,抑制流场的非对称性,使分离涡趋于对称;而固定频率的偏航振荡则破坏背风面分离涡的稳定性,激励非对称背涡产生随时间变化的周期性脱落。所进行的非定常运动与细长体流场耦合作用研究在国内外研究尚少。     

网格变形在俯仰振荡翼型中的应用

为解决航空器结构变形问题,在研究相关非定常流场基础上,建立了包含结构变形的非定常流场数值求解方法,运用无限插值法（TFI）生成任意时刻的非定常动态网格,通过与刚性动网格比较,证明了方法的可行性,并以NACA0012翼型为算例,对翼型俯仰振荡运动中表面柔性变形问题进行了数值模拟,发现俯仰振荡中叠加变形形式为一阶振型时,可以明显地提高升力、增加动稳定、提高升阻比。     

某型飞机俯仰方位控制系统设计与仿真

为提高大迎角飞行状态下的战斗机等作战飞机控制系统的稳态精度,改善其控制系统的动态性能,对某型飞机俯仰方位控制系统结构进行改进设计。通过引入积分环节,将其俯仰方位控制系统由0型系统变为I型系统,提高了控制精度及抗干扰能力。通过合理配置俯仰方位控制系统控制规律传动比,对其动态性能也进行了有效改进。仿真结果显示：所设计的飞机俯仰方位控制系统结构合理,控制规律传动比配置恰当,对大迎角飞行状态下的飞机,其控制系统稳定性得到提高,其动、静态性能得到改善。     

考虑多种物理效应的钝锥俯仰稳定性参数影响分析

首先采用二阶ROE格式计算了非定常NS方程下高超声速钝锥体的俯仰阻尼导数,通过与文献中内伏牛顿法和实验结果对比验证了本文所采用非定常计算方法的正确性.同时采用二阶ROE格式及双时间步长、刚性旋转动网格技术等,数值分析了NS方程计算条件下高空高马赫数钝锥体模型的强迫振动非定常流场特征,并运用最小二乘法辨识出其俯仰静稳定性导数与俯仰阻尼导数.通过计算结果对比分析了不同高度(含稀薄气体效应)、Maxwell滑移边界条件及五组元化学非平衡模型对钝锥体模型非定常流动特征及静、动导数的影响.结果分析表明,本文所计算的钝锥模型在所假设强迫振动下的俯仰稳定性参数随高度变化剧烈,且考虑滑移边界条件与化学反应会使得俯仰力矩迟滞曲线形态发生较大程度改变.准确计算高空、高马赫数条件下飞行器俯仰阻尼导数应考虑化学非平衡效应和稀薄气体效应的影响,是否考虑这些物理效应会关系到飞行器静、动稳定的程度.