固体火箭发动机零维两相燃烧室压强计算方法研究

为了更准确地预估含金属燃料固体火箭发动机的燃烧室压强,在压强计算中考虑两相流的影响,从一维两相喷管流动的求解出发,通过两相平衡流模型、两相常滞后模型、两相等温流模型、颗粒定温模型等模型的简化,分别推导不同模型下喷管中两相混合物的流量计算公式,再把流量公式应用到发动机零维内弹道理论中,推导并简化得到零维燃烧室平衡压强的计算公式。把压强公式用于HTPB推进剂固体火箭发动机和铝冰固体火箭发动机的燃烧室压强计算,结果表明,当固体推进剂中金属含量较高时(如铝含量为21%的HTPB推进剂发动机),用传统零维燃烧室压强公式预估的压强与实验误差较大,而使用合适的两相流模型和对应的零维燃烧室压强计算方法,在HTPB发动机中,能把压强预估结果与实验的误差降低到6%以内。如果使用多维内流场计算的方法,燃烧室压强预测结果的误差将下降到2.5%以内。结论发现在含金属固体火箭发动机的燃烧室压强计算中,考虑两相流的影响是必要的,而使用两相流修正后的零维燃烧室压强计算公式能够快速、较准确地预估这些发动机的燃烧室压强。     

基于大涡模拟的圆转矩喷管尾喷流强化掺混机理研究

为揭示2元圆转矩喷管尾喷流强化掺混的内在机制,应用大涡模拟(LES)方程对2种相同进、出口直径的喷管模型(轴对称、2元圆转矩)在Ma=0.8、高雷诺数(2×10~5)条件下进行了数值模拟计算。结果表明:与轴对称喷管相比,圆转矩喷管射流掺混效应增强,速度衰减快,核心区长度和高温区域面积减小。同时尾喷流拟序结构变化说明:2种喷管主要拟序结构均包含涡环、涡辫、发卡涡、螺旋涡等相似结构;但圆转矩喷管在射流近场诱导出的涡旋更丰富,边角剪切涡发展更快,形成明显的CVP结构,导致其射流柱失稳时刻提前、距离缩短;同时,喷管形式的改变使得射流剪切层内雷诺应力增大,速度脉动增强。拟序结构发展及雷诺剪切应力变化说明在射流流场中涡旋发展耗散速度增大、速度边界层脉动增强、射流柱易失稳是导致射流掺混增强的本质因素,为异形喷管的强化掺混机理提供了依据。     

基于CFD和CSM耦合的通用静气弹分析方法

提出了一种适用于有限元精细化建模的流固耦合插值点选择方法,通过RBF(径向基函数)方法实现流固耦合面的数据交换,实现了基于CFD/CSM(computational fluid dynamics/computational structural mechanics)耦合的通用非线性静气弹分析方法。以HIRENASD(high Reynolds number aero-structural dynamics)风洞试验模型为验证对象,数值结果很好地与风洞试验结构变形、气动压力分布吻合,验证了所发展非线性CFD/CSM耦合静气弹求解器的精度。详细研究了HIRENASD模型在大迎角(AOA)流动下的静气动弹性特性,以及该模型弹性变形对机翼气动特性影响规律。研究表明:HIRENASD弹性模型变形后其升力小于刚性模型;在小迎角范围内刚性、弹性模型升力差随迎角增大呈线性增长;当迎角大于4°后,升力差先减小后基本保持不变,呈非线性关系。      

基于双eN方法的翼身组合体流动转捩自动判断

发展翼身组合体复杂外形流动转捩自动判断方法,对高亚声速民机自然层流（NLF）机翼设计具有重要意义。使用多块结构化网格和三维雷诺平均Navier-Stokes（RANS）方程求解器,耦合边界层方程求解和基于线性稳定性理论（LST）的完全双e^N方法,发展了一套可同时计及Tollmien-Schlichting波和横流不稳定性扰动诱导转捩的翼身组合体流动转捩自动判断方法。对DLR-F4翼身组合体绕流进行了转捩自动判断,将得到的转捩位置与试验结果进行比较,验证了所发展方法的正确性。使用上述方法对配置自然层流机翼的中短程民机翼身组合体外形进行了数值模拟,并将结果与单独机翼的转捩位置进行了对比,结果表明机身三维位移效应增强了自然层流后掠机翼边界层的横流不稳定性强度,导致翼根转捩位置提前至前缘区。     

透平机械动密封动力特性数值方法研究

密封动力特性系数是评价透平机械转子稳定性的重要参数。对现有密封动力特性数值方法进行了综述,为了比较分析不同的密封动力特性数值方法,分别建立了基于控制体方法的双控体求解模型和基于CFD方法的稳态旋转坐标系与瞬态转子平动、转子单/多频椭圆涡动求解模型,综合比较分析了控制体法、CFD稳态旋转坐标系法和CFD瞬态平动法、单/多频法之间的差异以及适用范围,研究了转子涡动频率对密封气流力和动力特性系数的影响。研究结果表明:控制体法求解速度较快,但求解精度较低,稳态旋转坐标系法只适用于转子做小轨迹同心涡动的求解模型,瞬态平动法求解速度较慢;瞬态单/多频法考虑了转子的涡动频率,更符合密封的实际工作情况;应用瞬态单/多频法得到的密封气流力变化频率与转子涡动频率相同,但气流力的相位滞后于转子涡动位移的相位;密封的直接刚度系数随着转子涡动频率的增加而增大,交叉刚度系数和阻尼系数的绝对值随着转子涡动频率的增加而减小。     

Arcadia建模方法与SysML建模方法比较研究

基于模型的系统工程与传统系统工程相比,有着形象化、具体化、沟通效率高等多种优势。基于系统建模(System Modeling Language,简称SysML)语言进行建模用途广泛,但是通过实践,针对系统架构和系统运行场景分析,基于Sys ML的建模方法特点并不突出。而Arcadia的基于模型系统工程建模有着简单、清晰等特点,对于系统架构和系统定义均有着Sys ML所不具备的特点。通过系统设计的目标与特点,对基于Sys ML的建模方法和Arcadia方法进行比较研究和实践,总结两种方法在基于模型的系统工程领域中的优劣。     

对航空货运运价的制定方法研究

AResearchonthePriciugMethodoftheAirCaryoTransport自中国的改革开放政策实施以来,我国国民经济的迅速发展,进出口贸易的大幅度增加和人民生活消费水平的提高,我国航空运输业得到了迅速的发展。特别是航空货运的发展已大大超过航空客运,近10年以来,我国航空货运量年均增长20％,为世界肮空货运发展速度的3倍。正是由于货运业的高速发展,使得许多航空公司开始从日益加剧的客运市场慢慢转向市场前景看好的航空货运市场。航空货运业有望成为在新世纪推动中国航空运输业发展的新的经济增长点。虽然我国航空货运的发展速度令世人瞩目…     

导热-辐射耦合传热的多尺度分析和数值模型

复合材料高温传热特性的准确预测对飞行器热防护结构的设计有重要意义,相关模型也是国家数值风洞工程中多相多介质计算模型的重要组成部分。针对周期性结构复合材料的高温传热问题,利用多尺度渐进分析方法,对包含导热方程和辐射传输方程的导热-辐射耦合传热模型开展了研究。建立了表征单元尺度模型及宏观平均导热-辐射耦合传热模型,获得了材料宏观等效导热系数与表征单元模型间的关系,发现宏观等效辐射吸收和散射等系数可通过表征单元内的体积平均求取。结合有限容积方法与格子Boltzmann方法,建立了复合材料导热-辐射耦合传热多尺度数值模型。采用二维常物性材料传热过程的模拟验证了多尺度模型的有效性,结果表明所建立的多尺度模型能够对温度场给出准确高效的计算结果。该方法有助于为复合材料传热特性的预测提供数值手段。     

考虑燃料雾化的气液两相连续旋转爆轰数值模拟

为研究气液两相连续旋转爆轰发动机燃烧室内爆轰波的传播特性,以汽油为燃料,富氧空气为氧化剂,建立了欧拉-拉格朗日模型进行二维数值仿真,其中气相方程采用时空守恒元与求解元方法求解,液相方程采用标准四阶龙格库塔法求解。在两相旋转爆轰模型还考虑了液滴雾化破碎过程。计算结果表明：起爆后形成的初始爆轰波经过初始燃料填充区域后逐渐衰减,随后入口附近新生成的压力波经过一系列发展形成了自持稳定传播的旋转爆轰波;旋转爆轰波的传播模态受燃料与氧化剂的喷注压力和氧化剂填充比影响,在不同工况下旋转爆轰波呈现出4种传播模态,即稳定单波模态、稳定双波模态、不稳定双波模态和不稳定单波模态;在双波模态工况下,燃烧室内初始只形成1个爆轰波,后由入口附近局部爆炸产生的压力波发展为新的爆轰波,转化为双波模态后爆轰波的强度略有下降,但燃烧室整体推力更加平稳。     

带凹坑圆孔气膜冷却流动传热特性

利用大涡模拟方法对平板带凹坑气膜孔冷气射流与主流的相互作用展开研究,并结合红外热成像壁温测试和粒子图像测速结果进行实验验证,揭示了凹坑对气膜冷却的强化机理.结果表明:凹坑在气膜孔下游诱发反肾型涡对,将冷气向壁面卷吸,改善了气膜贴壁性;与圆孔气膜冷却相似,凹坑气膜冷却流场也出现马蹄涡、发卡涡等大尺度拟序结构,但与圆孔下游...