复杂旋翼流场的耦合欧拉-拉格朗日数值方法

针对影响旋翼流场求解精度的关键因素“桨叶复杂近体流动”和“尾迹涡畸变”,结合计算流体力学(CFD)方法和黏性涡方法,发展了一套适合于复杂旋翼涡流场分析的耦合欧拉-拉格朗日数值方法:为捕捉桨尖三维效应、激波等细节流场特征,在桨叶近体区域采用CFD方法对其进行求解;针对高雷诺数旋翼流场中桨尖涡的紧凑结构特点,引入黏性涡方法建立了高分辨率的尾迹求解模型;两计算域间的信息交换采用了集中涡源法和边界修正法.应用所建立的计算方法,以旋翼CFD标准验证试验(Caradonna-Tung旋翼)为算例,对尾迹影响明显的悬停状态进行了数值模拟,通过对比耦合边界处流场特征及桨叶表面压力系数分布,验证了方法的有效性.此外,还从旋翼尾迹捕捉精度、涡量耗散特征及计算时间等方面对不同计算方法进行了对比分析,结果表明耦合方法可充分发挥CFD和黏性涡方法各自的优点,在旋翼流场数值模拟方面具有独特的优势.     

液体火箭发动机中频耦合振荡初步研究

为了研究常规液体火箭发动机研制中出现的中频流量型耦合振荡问题,对发动机耦合系统中各环节建立了满足声学分析的线性化频域模型.通过比较供应系统和燃烧室的幅频响应,并分析耦合系统的自由振荡复频率,以判断发动机的耦合稳定性.计算结果符合发动机试车规律,能够有效地反映该类型不稳定的特征.进一步讨论了燃烧时滞、喷注压降、管路长度和安装节流圈对系统耦合不稳定的影响规律,并揭示了该类型不稳定的机理.     

推力室状态对其声学振型及其阻尼特性的影响

为确定工程中冷态条件下获得的推力室声学特性能否表征真实条件下的声学特性,研究了冷态无流动、热态气相流动和湍流两相燃烧三种状态下推力室声学振型及其阻尼特性。在推力室稳态流场中的有限区域施加数值定容弹,激发其具有多模态声学振型的大幅值压力振荡,采用衰减时间和半带宽来定量评价所激发的不同声学振型压力振荡衰减快慢,进而获得其阻尼特性。在相同过载比的数值定容弹激励下,在冷态条件下能激发包含更多声学振型压力振荡,且该振荡衰减时间更长,相同振型压力振荡衰减比热态条件下慢。在冷态条件下,一阶切向振型振幅最大,为最容易被激发声学振型;一阶纵向振型半带宽最小,为最难衰减的振型。在热态条件下,一阶纵向振型为最容易激发声学振型,也为最难衰减声学振型。从所激发的主要振型及其相对衰减的快慢来看,冷态条件下获得的声学特性能够表征真实条件下的推力室的声学特征。     

大机动飞行中孤立旋翼的非定常空气动力模型

为提高旋翼他轴响应的预测,气动模型应当计入机动旋翼的尾迹畸变效应,而准确表示尾迹弯曲参数Kn是关键。为表达机动导致的各种效应,增广Peters-He有限状态尾迹模型可以吸收随飞行状态变化的Kn值。根据烟流实验,该参数随旋翼前飞速度、拉力和机动角速率而变化。旋翼的空气动力实验表明,旋翼轴接受准阶跃输入后,同轴及他轴的空气动力响应有过冲,且俯仰速率愈大,过冲幅值愈大。理论和实验结果的对比显示,为了正确计算机动飞行中的旋翼气动力,特别是操纵的他轴响应,应在模型中采用随状态变化的尾迹弯曲参数Kn。此外,在当前模型基础上联合动态尾迹畸变模型,计算出的旋翼力和力矩更接近试验数据。     

飞行器纵向阻尼动导数直接测量实验研究

介绍了新研制的“FL-8风洞大幅升沉和俯仰耦合振荡实验系统”及采用该系统产生纯俯仰运动测量飞行器纵向阻尼动导数的试验技术，给出了部分试验结果。实验结果表明，阻尼导数和时差导数在中、大迎角时存在明显的非线性耦合，采用纯俯仰运动直接测量的阻尼导数在整个迎角范围内都是合理的。     

管内激波串振荡和壁面脉动压力特性

分析了管内激波串振荡引起的壁面压力脉动特性.采用直联风洞实验方式结合动态压力测量技术获得了矩形管道内不同激波串位置下的壁面压力分布数据,对所得压力脉动幅值、功率谱密度曲线等结果进行分析,探讨激波串流动非定常特点与引起压力脉动的机制.     

悬停状态旋翼固定尾迹分析中三角级数的应用

为了在固定尾迹分析中更好地描述桨叶附着涡环量分布,结合叶素理论,采用一个三角级数表达附着涡环量的连续分布,并应用于悬停状态下单旋翼和共轴双旋翼的气动特性分析,通过小量计算即获得连续光滑的附着涡环量分布,将其应用到推导的公式中计算出连续分布的诱导速度,进而分析了共轴双旋翼上下旋翼气动干扰及扭矩平衡条件.计算结果与实验结果进行了比较分析,验证了该方法的有效性,并得出结论:下旋翼总距变化对上旋翼诱导作用影响微小,但对自身诱导作用明显;当上下旋翼扭矩平衡时,下旋翼总距比上旋翼总距约大1°.     

基于N-S方程的尾迹面法翼型气动阻力计算

研究了计算翼型阻力改进的方法——尾迹面法及尾迹面积分的位置和相应的积分技术.在亚跨声速时分别采用了表面积分和尾迹积分求得给定翼型RAE2822的阻力.结果显示,2种方法具有相同的结果,表明尾迹面积分方法是有效的,与实验值吻合较好.尾迹面法中,尾迹面位置应处于离后缘相对弦长距离0.6~1.0之间.尾迹面积分方法中积分结果不依赖于物体的详细几何外形,可以预计对曲率变化大的三维复杂外形,该方法有更大的优势.      

大气准两年振荡对赤道行星波上传的影响

本文讨论了热带大气行星波在QBO风场中向上传播的规律,并运用数值模拟的方法,得到了开尔文波和混合罗斯贝重力波在QBO不同相位下交替上传的结果。试图用QBO对行星波的调制解释QBO在日地相关性中的作用.      

基于WENO-分段线性格式的旋翼流场数值模拟

为了降低二阶JST(Jameson-schmidt-turkel)格式导致的数值耗散,发展了一套适合于格心格式的基于加权本质无振荡(Weighted essentially non-oscillatory,WENO)-分段线性格式的旋翼流场数值模拟方法。采用运动嵌套网格方法生成围绕旋翼的网格系统,主控方程选择Navier-Stokes(N-S)方程。为了更加有效地对桨尖涡等流动细节进行捕捉,在笛卡尔背景网格上采用七阶Roe-WENO格式计算对流通量;在桨叶贴体非结构网格上采用二阶精度的Roe-分段线性格式计算对流通量。时间离散采用了高效的双时间隐式LU-SGS(Lower upper symmetric Gauss-Seidel)方法进行时间推进。最后,应用上述方法对悬停状态的C-T(Caradonna-tung)旋翼和Helishape 7A旋翼进行了数值模拟,将数值计算结果与实验数据进行了对比,计算值与实验值吻合较好;并将桨尖涡模拟效果与二阶JST格式的模拟效果进行了对比。对比结果表明:本文方法能有效对旋翼流场进行计算,且在相同计算条件下,WENO-分段线性格式能够更有效地捕捉旋翼涡流场的流动特性,表明在计算旋翼涡流场时WENO-分段线性格式相比传统二阶JST格式具有更低的数值耗散。