应用傅里叶时间谱方法求解二维跨音速流动问题解的精度研究

为了从时间离散精度的角度评估傅里叶时间谱方法的实际计算效率,对用该方法求解二维跨音速流动问题的解进行了精度研究。通过求解Euler方程模拟俯仰振动NACA0012翼型的周期性非定常流动。结果表明:傅里叶时间谱方法可以使用很大的时间间隔较为准确地模拟有激波存在的周期性跨音速流动;对于此类流动问题,相比二阶向后差分公式法,傅里叶时间谱方法在预测翼型表面压力系数积分量特别是升力系数的时间变化规律方面具有明显的计算效率优势;用于模拟翼型表面压力系数自身的时间变化规律时,傅里叶时间谱方法具有不低于二阶向后差分公式法的计算效率。     

不同钝体宽度下的驻涡燃烧室排放试验

为了研究空气流量分配对驻涡燃烧室对排放特性的影响,了解对驻涡燃烧室内污染物生成的过程及其影响因素,设计了一个能够改变中心钝体宽度、仅凹腔供油的驻涡燃烧室.在常压下对该驻涡燃烧室进行了排放特性试验,进口温度保持200℃.试验中,燃烧室进口马赫数为0.15~0.3.影响排放的因素主要包括雾化质量、凹腔当量比以及与进口马赫数相关的驻留时间等.总体来说雾化质量、凹腔当量比的提高对降低CO和HC的排放是有利的,但是这会使NOx排放增加.在低凹腔当量比时,CO排放曲线变化下降比较平缓,甚至出现上升趋势,而HC排放曲线比较陡峭.这是由于HC的消耗速度比CO消耗速度快,随着凹腔当量比的增加,供油压力提高,燃油雾化粒径变小,燃油蒸发时间缩短,使HC排放快速减少,中间产物CO大量产生而来不及消耗.凹腔当量比进一步上升时,由于燃烧温度的提高,使得CO排放快速减少.在燃烧室内燃烧过程中,NOx的形成和消耗是非常复杂的过程,目前只能作一些定性的分析,而CO和HC的反应过程相对简单.通过对不同钝体槽宽下,具有相似凹腔前壁流量的工况的比较,发现CO和HC的形成主要受凹腔内工作状况影响,而NOx的形成过程更复杂,主流也对其产生着重要的影响.      

以柴油为燃料的燃烧室燃烧效率实验及分析

地面动力常常需要改进利用航空燃气轮机技术,为了研究燃气轮机燃烧室在低工况下使用柴油燃料的燃烧效率问题,采用燃气分析试验方法,对燃气轮机燃烧室在不同进口总温、空气流量和油气比等工况下的燃烧效率进行了研究。结果表明:油气比较低时,燃烧区温度低,燃烧不完全导致燃烧效率急剧下降,随着燃烧室油气比的增加,燃烧效率逐渐接近100%;进口空气流量以及进口总温的增加会提高空气雾化喷嘴的雾化能力以及燃烧温度,燃烧更加充分完全,燃烧效率由96%左右提高至99%以上;总结归纳得到了适用于柴油燃料的燃烧效率预估经验关系式。     

针栓式喷注器锥形液膜破碎特性试验

采用高速摄影获得了针栓式喷注器在不同喷注压降和结构参数下的表面波破碎图像,测量了锥形液膜的破碎长度和破碎时间,研究了变工况时液膜破碎长度和破碎时间的变化规律.试验结果表明:在喷注压降不变的条件下,针栓式喷注器能够实现流量的线性调节.针栓式喷注器设计时,在合理的推进剂动量比范围内,狭缝宽度应尽量取小.液膜在低工况时破碎得更快.液膜破碎长度和破碎时间均随喷注压降的增加而减小.      

固体推进剂二维绝热微尺度燃烧模型的数值研究

为实现固体微推力器工作过程的一体化模拟,基于Fluent计算软件的二次开发功能(UDF)和简化化学动力学模型,实现了固体推进剂的二维气-凝相绝热微尺度燃烧模型的建立,该模型针对固体微推力器所用双基推进剂,包含两步凝相反应和五步气相反应,燃速、推进剂表面温度和组分质量分数基于燃面物理特性计算得到,并考虑了粘性作用对气相和凝相反应的影响。针对0.5MPa,1.0MPa,2.0MPa和5.1MPa四种工况进行了计算,结果表明,高压工作环境下出现发光火焰区,且随表面压力增大而逐渐靠近壁面,凝相反应区厚度和嘶嘶区、暗区主要反应物在燃面的质量分数随推进剂表面压力增大而减小。对称面处推进剂燃速,推进剂表面温度和气相火焰结构与实验结果基本一致。由于壁面附近较高的粘性作用,气相火焰在壁面位置更加靠近推进剂燃面,并导致壁面位置推进剂燃速高于对称面位置。该模型实现了二维环境下考虑分步凝相反应的推进剂绝热燃烧模型的一体化计算,较好地拓展了原模型的应用范围。     

北京首都机场畅通滑行时间研究

随着机场规模不断扩大和流量持续增长,准确把握拥挤态势成为运行管理中的难题。美国联邦航空局(FAA)和欧洲空管局(EUROCONTROL)采用畅通滑行时间(unimpeded taxi time)作为量化分析机场运行效能的标杆。国内笼统地规定了少数机场的平均滑行时间作为机场运行指标,但无法反映不同机场、不同停机位、不同起飞跑道造成的运行差异。因此可借鉴机场畅通滑行时间这个指标来考核机场效能,但是欧美机场的该指标数值并不完全适合国内机场,因此在欧美计算方法的基础上提出了改进的畅通滑行时间计算方法,并构建多元线性回归模型,运用北京首都机场实际运行数据进行检验,证明了该方法的可行性。     

燃气轮机燃烧室网络式流量分配设计

采用一维网络式方法对某型燃烧室进行流量分配设计,并作三维数值模拟验证。网络式方法将燃烧室分割为若干个独立的单元,各单元互相关联构成整体网络。整体网络通过压力关系式求解,得到各元件压力损失、质量流量等性能参数。对某型燃烧室结构进行网络化建模并作优化设计,目标包括:1各排主燃孔及气膜冷却孔流量分配;2燃烧室进出口压降4.2%。对优化后的燃烧室结构作三维数值模拟并与一维网络优化结果对比验证。结果显示:二者吻合很好,各排孔的流量分配以及燃烧室总压损失目标均实现良好,显示一维网络方法在燃烧室复杂结构设计中的精确性和可靠性。     

一种FMCW SAR地面运动目标谱图域参数估计方法

对调频连续波(FMCW)合成孔径雷达(SAR)地面运动目标的参数估计方法进行了研究,采用相位中心偏置天线(DPCA)技术对地杂波进行抑制,分析了载机及地面运动目标连续运动对回波信号的影响,推导了采用DPCA技术引入的回波慢时间包络(STE)项与地面运动目标参数之间的关系。在此基础上,提出了一种地面运动目标谱图域参数估计方法,该方法首先利用Radon变换在谱图域估计导致回波信号距离走动的目标等效径向速度,并对距离走动进行校正;其次,在谱图域中提取运动目标回波幅度,根据STE项引起的回波幅度变化与目标方位向速度之间的关系,估计目标的方位向速度,并进一步求解相应的目标径向速度。所提方法能够在谱图域完成地面运动目标二维速度估计,最后的仿真实验验证了所提方法的有效性和可行性。     

一种适用于三维混合网格的GMRES加速收敛新方法

为提高流场计算收敛效率,发展了一套适用于三维混合网格Naiver-Stokes方程求解的并行广义最小残差（GMRES）隐式时间推进方法。该方法由科学计算可移植扩展工具包（PETSc）中的Krylov子空间求解器实现,线性方程系统中的系数矩阵直接以显式给出以提高算法的稳定性。为进一步提高GMRES方法的收敛速度,对非结构网格的序号进行了重排序,使得系数矩阵的非零元素尽量向主对角线靠近。利用所发展的GMRES方法,完成了对ONERA-M6机翼、AIAA阻力预测会议通用研究模型（CRM）等算例的计算,计算结果与试验结果吻合良好。通过与其他隐式推进方法进行比较,对算法的收敛特性进行了研究。结果表明,所发展的GMRES方法计算更加稳定,残差下降速度相对LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法更快,尤其是气动力系数向着收敛解逼近的速度更加明显,提高了计算效率。     

低马赫数下多凹腔燃烧室非稳态燃烧过程

作为稳定火焰的有效手段之一,凹腔构型在冲压发动机燃烧室研究中占有重要地位。在对以煤油为燃料的多凹腔燃烧室冷/热态流动特性分析的基础上,重点研究低进口马赫数条件下燃烧室点火起动初期非稳态过程。结果表明：上游凹腔内大涡结构有助于提高燃料的驻留时间,未燃混气被高速主流带入下游凹腔内继续反应,进一步提高燃烧效率;燃油喷射速度决定被卷吸进回流区的燃油质量分数的大小,进而影响燃烧效率高低;燃烧室点火起动初期出现了主流熄火、火焰逆流传播以及主流再着火等复杂现象,火焰逆流传播现象是在上游凹腔内燃料自燃与下游燃烧释热压缩来流两种机制共同作用下完成的。