氧化剂含量和粒度对NEPE推进剂燃速影响的模型化

以高能固体推进剂热分解特性和燃烧模型的研究成果为基础,建立了由化学结构参数计算NEPE推进剂的燃速和压强指数的公式,计算了氧化剂组分含量和粒度对燃烧性能的影响。经验证,计算结果与实测燃速值的偏差全部在±20%以内,其中80%的偏差在±10%以内。这说明所建立的模型基本合理,编制的NEPE类推进剂燃速计算程序基本可行。     

发动机短舱内外流场与结构温度场耦合计算

采用CFX软件,利用SST k-ω湍流模型,同时引入了辐射模型考虑热辐射的影响.采用CFX的共轭传热方法,在流体域求解N-S方程,固体域求解热传导方程,然后在流固耦合交界处交换数据,实现了流场与结构温度场的耦合数值分析.为了验证此方法的正确性,对文献中轴对称喷管进行了数值模拟,所得喷管壁温与试验数据吻合良好;另外进行了高超声速钝体流场与结构温度场的耦合计算,计算所得激波位置与试验结果一致,固壁非定常温度场分布与文献结果吻合.在此基础上,进行了航空发动机短舱内外流场与结构温度场耦合的数值研究,并得出了几个短舱外表面蒙皮冷却的可行性措施,为飞行器红外隐身设计提供参考.     

桥梁断面静风荷载的格子Boltzmann方法数值计算

通过引入反映湍流涡粘性的湍流松弛,得到了模拟高雷诺数湍流的BGK方程.在速度相空间、物理空间和时间上对BGK方程进行离散得到了三维十九速离散速度模型;结合分区计算技术,设计了格子Boltzmann并行算法;根据亚格子Smagorinsky模型,提出了直接从粒子分布函数计算湍流松弛时间的方法.用开发的并行计算程序对分体双箱截面和闭口箱梁截面的静风荷载进行了数值识别,得到的静力三分力系数和流场压力分布与风洞试验结果及CFD宏观方法计算结果吻合,并从表面压力分布入手分析了两种桥梁截面的绕流特点.     

马赫数2.5超燃室中超/亚声流流向涡掺混冷态流场数值研究

结合实验,对一种带亚声速预燃室和流向涡掺混器的超声速燃烧模型燃烧室,在其进口马赫数为2.5的来流条件下,进行了冷态流场的数值研究.计算与实验得到的燃烧室沿程压力分布相当一致,计算与实验得到的激波结构也基本吻合.计算结果表明:首先,在马赫数2.5的来流条件下,亚燃预燃室易于达到启动状态;其次,流向涡掺混器增强超/亚声速流之间的掺混的效果明显,但其掺混深度尚有限;最后,流向涡掺混的超燃室掺混段有着复杂的激波膨胀波波系,波涡干涉和激波附面层干涉结构.     

新概念燃烧室应用高温空气燃烧技术的数值模拟

为了将高温空气燃烧技术(HiTAC)应用于航空发动机燃烧室设计,本文设计了一个全新概念的燃烧室.该燃烧室将小型预燃室燃烧产生的高温贫氧烟气直接喷入与其连接的火焰简内组织成高温贫氧空气燃烧.本文对该燃烧室内的燃烧状态进行了数值分析.结果表明:该新概念燃烧室具有温度均匀分布,壁面有效冷却,NOx低排放,节约燃油等优势.     

巡飞器无舵偏情况下的气动数值分析

采用SSTk-ω二方程湍流模型,通过解耦求解雷诺平均N-S方程实现对巡飞器无舵偏情况下湍流流场的数值模拟。对巡飞器低雷诺数绕流流场进行了计算,分析了其在无舵偏情况下气动特性随迎角、侧滑角的变化情况。仿真结果表明,巡飞器采用充气式机翼后,具有很好的失速性能与较好的大迎角稳定性。     

30P-30N多段翼型复杂流场数值模拟技术研究

采用"亚跨超CFD软件平台"(TRIP2.0)数值模拟了30P-30N多段翼型的复杂流场,主要目的是考核湍流模型、转捩位置对多段翼型压力分布和典型站位速度型的影响。本文通过求解任意坐标系下的雷诺平均的N-S方程,采用多块对接结构网格技术,在与相应试验结果对比的基础上,详细研究了SA一方程湍流模型、SST两方程湍流模型、不同的转捩位置对该翼型压力分布和典型站位速度型的影响。本文的研究结果表明,采用全湍流模拟方式可以较好地模拟该多段翼型的压力分布,但对速度型的模拟精度较差;模拟试验的转捩位置可以改善主翼附面层与前缘缝翼边界层尾迹区的模拟精度;采用微吸气技术推迟前缘缝翼的转捩位置,可以进一步提高缝翼尾迹区的数值模拟精度。     

一种剪切修正的涡黏性亚格子模型及其在槽道湍流中的应用

提出了一种利用壁面附近强剪切对涡黏系数进行修正的亚格子模型,称为剪切修正的涡黏模型.该模型能使涡黏系数在壁面附近自动满足壁面修正,在脉动为零的地方,涡黏系数自动为零,保证了壁面附近涡黏系数趋于零的特性.用该模型对槽道湍流壁面无吹吸、有吹吸,较低和较高雷诺数的情况进行了数值模拟,与直接数值模拟的结果进行了对比,结果吻合较好.      

基于需用速度增益曲面的大气层外超远程拦截导引方法

考虑拦截器使用耗尽关机固体推进系统的情况,提出需用速度增益曲面概念,设计了基于该概念的大气层外超远程拦截导引方法。根据Lambert导引,给定拦截时间就有唯一的指令推力方向与之对应,导引过程分为两个阶段:零控拦截到达阶段,选择最优拦截时间,需用速度增益曲面迅速与助推时间-拦截时间平面相切;零控拦截保持阶段,拦截时间不断减小,保持需用速度增益曲面在助推时间-拦截时间平面上滑动,消耗多余的推进剂。利用拟零控拦截概念使两个导引阶段平缓过渡,避免了导引方法切换时推力方向的跳变。导引律计及了J2项摄动对滑行段弹道的影响,导引精度对推进系统参数偏差的鲁棒性强。仿真结果表明本文制导方法用于大气层外超远程目标拦截,脱靶量仅为km量级,推进系统参数偏差对导引精度的影响很小。     

适于黏性计算的自适应笛卡儿网格生成及其应用

复杂外形的黏性网格生成已经成为计算流体力学(CFD)在工程应用中的主要瓶颈,高效而鲁棒的网格生成技术研究显得尤为迫切。使用基于物面外形自适应加密的笛卡儿网格方法生成计算网格,投影方法拟合壁面边界,将投影得到的柱形单元沿法向分层获得适于黏性计算的网格。应用特征恢复技术恢复凹面特征并改善凸面处的网格质量。流场解算中应用最小二乘法重构获得2阶精度,Roe格式计算无黏通量,Holmes-Connell方法计算黏性通量。加入Spalart-Allmaras一方程湍流模型以获得湍流解。上下对称高斯 赛德尔迭代(LU-SGS)格式隐式时间推进,对网格重排序保证格式的稳定性和平衡性。算例显示网格生成方法是快速的和鲁棒的,流场计算精度满足工程需求,因此该方法非常适合应用于实际工程计算。