超燃冲压发动机一维稳态跨声速流动奇异初值问题的一个解法

超燃冲压发动机燃烧室一维稳态跨声速流动方程在临界声速点存在奇异初值问题，现有的基于L’Hospital法则的求解方法在原理上存在较大的初值误差，影响一维稳态跨声速计算的精度。为此，本文提出了一种基于变量代换的改进算法，通过定义新的流动变量W=-Ma^2＋2Ma，构造出了非奇异的一维稳态跨声速流动方程。消除了微分方程的奇异性，有效的解决了一维稳态跨声速流动计算过程中存在的奇异初值问题。     

基于随机配点的高超声速飞行偏差演化分析方法

针对不确定性条件下高超声速飞行偏差演化问题,提出一种基于随机配点法的偏差演化分析方法。在分析高超声速飞行偏差主要来源的基础上,综合考虑计算精度和分析效率,基于广义多项式混沌理论,采用随机配点方法进行偏差演化计算,进一步给出随机配点生成的稀疏网格技术。数值仿真表明,随机配点法的引入有效解决了高超声速飞行器系统偏差演化问题,可为工程设计提供有益参考。     

跨声速涡轮叶栅两点优化目标函数研究

针对跨声速涡轮叶栅单点优化方法难以获得整体工况性能提升、多点优化方法难以确定合理目标函数形式的问题，提出了两点优化的方法。为了节约优化时间成本，优化过程采用EIF （Equivalent inviscid flow） 模型进行数值模拟，通过添加惩罚函数保证叶栅满足设计流量和负荷要求，并采用叶栅效率线性平均的目标函数形式进行评价。选择两组跨声速涡轮叶栅进行优化设计，并利用CFD方法分析叶型变化对流场马赫数、激波和损失产生的影响。结果显示，所提出的优化设计方法在保证设计工况性能的同时，能够提升叶栅整体工况性能。通过流场分析，揭示了激波结构变化对不同工况损失影响的定性规律。综合全文研究后，给出了一种适用于跨声速涡轮叶栅两点优化设计的目标函数形式。     

跨声速轴流压气机特性计算用激波损失模型

基于实验结果,发展了一种适用于跨声速叶型的新激波损失模型,该模型考虑了压气机工作状态从堵塞点向喘振点变化时激波系结构的变化。将新激波损失模型应用于某两级跨声速轴流压气机特性预测中,计算与实验结果的比较表明,在跨声速轴流压气机非设计点性能预测中,该模型能够给出更为准确的结果,对于效率曲线的改进尤为明显,并且激波损失系数在全工况下的变化规律更加合理,这表明新激波损失模型能够准确地反映出跨声速压气机内部的流动机理和激波损失大小。     

高超声速飞行器底层性能评价指标的可行性分析

针对控制性能指标是否适用于高超声速飞行任务评价，研究了底层指标的可行性分析方法。首先分析了高超声速飞行器的模型特性，指出底层指标具备高超声速飞行器评价可行性的两点要求。根据上述要求将可行性分析转化为对于任务变量到底层评价指标的映射分析，并提出了一种基于统计实验分析的检验方法。结合高超声速飞行器模型特点，进一步设计了可行性分析流程。然后，分别以静态任务要求和速度、高度跟踪任务要求为例，对应提出频域性能指标和时域性能指标的可行性分析方法。最后针对频域性能指标评价静态任务的可行性设计仿真实验加以验证。仿真结果表明可行性分析的结论与军标观点一致，而且所提出的分析流程具有较好的可实现性。     

跨声速偶极子网格法的改进和在颤振计算中的应用

改进了跨声速偶极子网格法（ＴＤＬＭ）,且在颤振计算中能自动求颤振点速度,从而较大地提高了计算速度,并有利于方法在较大马赫数及较大减缩频率时的应用。对跨声速飞机机翼作了跨声速非定常气动力计算,对跨声速民航机机翼振模型作了颤振计算,与ＴＤＬＭ的比较表明,两者结果有较好的一致性。     

宽包线吸气式高超声速飞行器外形优化研究

针对吸气式宽包线高超声速飞行器的气动优化问题，基于任务要求建立了基于多点权重分配的气动外形优化模型，并采用“CFD+准一维流”方法开展了气动性能分析。为兼顾气动外形优化的效率与精度，通过改进现有的并行加点策略，发展了一套基于代理模型与梯度算法的分层优化框架，并采用函数算例对改进后的加点策略进行了验证。对吸气式高超声速飞行器的气动外形进行了分层优化，在满足各学科约束的情况下使飞行器在各个优化评估点处的气动性能均有所提升。     

一种跨声速轴流压气机性能预测的数值方法

采用二维流线曲率法数学方法,考虑到不同因素对落后角和损失的影响,发展了一种适应于跨声速轴流压气机的落后角和损失模型.对轴流跨声速压气机NASA(National Aeronautics and Space Administration)Stage 37进行了数值计算,得到了设计点与非设计点的特性曲线,并与试验数据进行了...     

跨声速轴流压气机非轴对称端壁造型优化设计

为探索非轴对称端壁提高跨声速轴流压气机性能的潜力,基于三维数值模拟优化平台,针对跨声速轴流压气机转子进行了非轴对称轮毂优化造型,深入分析了非轴对称轮毂造型对跨声速压气机流场结构及性能的影响机理。结果表明:优化造型后的非轴对称轮毂同时提高了该跨声速压气机的效率和压比,最高效率点效率提高0.31%,压比提高0.31%,并改善了压气机的变工况性能。     

跨声速多级轴流压气机非设计性能预测

基于公开发表的研究成果, 发展了一种适用于跨声速轴流压气机的新型损失和落后角模型, 并进行了3-D的修正, 运用改进的流线曲率法对某两级跨声速轴流压气机进行了数值模拟, 得到了设计点和非设计点的特性曲线.通过计算结果与实验值的比较, 验证了本模型和方法在设计点和非设计点均具较高的精度, 可用于工程计算.      

全隐式TVD格式求解二维Euler方程的对角化算法

本文采用有限体积法,Chakravarthy提出的基于近似Riemann解的全隐式TVD差分格式求解二维非定常Euler方程。采用牛顿法对半离散化后的全隐式差分方程进行线性化,然后对线性化后的方程采用近似因式分解法处理,得到了两组三对角块矩阵方程组。在求解该方程组时,采用了矩阵对角化思想使原来的4×4三对角块矩阵方程组转化为4组分离的三对角代数方程组,因而大大节省了求解方程组所耗费的计算时间,提高了计算效率。通过对翼型的跨声速绕流问题及圆柱超声速流动问题的计算,证实了本文算法不但简便可靠,而且还具有较强的通用性。采用TVD格式捕捉到的激波分辨率高,上、下游无任何波动。     

二维高亚声速空腔流激振荡的数值模拟研究

采用基于P .L .Roe的近似Riemann解的修正Osher&Chakravarthy (MOC)三阶TVD有限差分格式 ,数值求解二维雷诺平均全Navier Stokes方程 ,并用Cebeci Smith代数湍流模型 (对腔内区域作修正 )来模拟湍流效应 ,时间方向的积分采用四阶Runge Kutta方法 ,通过对GAMM超音速前台阶绕流的计算 ,验证了格式及程序的有效性 ,对高亚声速来流下的空腔流动作了数值模拟研究 ,取得了较好的效果     

TVD格式计算一维Euler方程间断解存在的几个理论问题

本文讨论了TVD格式按本文第二节A段意义推广到拟线性方程组的几个理论问题。指出:1.对于Riemann间断分解,Euler方程的解和Riemann不变量都没有TVD性质,因此不能将TVD格式自然推广;2.TVD格式是提高分辨率的一种途径,它不能保证得到物理解,能得到物理解的TVD格式一定是耗散守恒格式,因此将TVD格式推广时,必须保证格式是耗散的,否则仍会得到非物理解。并给出了一阶迎风格式,TVD2,UNO2得到的若干非物理解算例;3.非物理解现象只是发生在具有剧烈膨胀的单波区的计算上,而某些TVD格式若会得到非物理解,则一定发生在特征值变号的单波区的计算上。     

叶轮机转子内流动的三维跨声Navier-Stokes方程组数值计算

应用有限体积法及Baldwin-Lomax代数紊流模型求解叶轮机转子内流动的三维可压缩非定常薄层Navier-Stokes方程组,计算一个绕X轴等转速旋转的转子内流场问题,用当地时间步长方法来加速收敛。采用隐式算法对对流通量进行求解,而扩散通量采用显式算法。采用高阶TVD下的Roe格式计算对流通量,采用中心差分格式计算扩散通量。关于守恒变量Q的非线性方程组,采用牛顿迭代方法求解,线化后的方程用对称的Gauss-Seidle方法求其渐近解。应用该算法模拟了转子NASARotor67流场,与试验做了大量的对比研究,包括性能、等马赫线以及总压和总温等等。在最高效率点附近计算与试验比较接近,但在靠近失速点时计算与试验略有差别。      

应用TVD格式分离求解不可压N-S流动

提出一种新的将连续方程与动量方程分离求解的TVD引用方法。对动量方程中的对流项采用三阶精度、迎风TVD格式离散, 粘性项和压力项采用中心差分格式离散。压力采用压力替代法求解。计算所用的网格是具有交错网格主要特点的修正非交错网格, 并采用一种压力修正法和G-S中心对称迭代法加速收敛。给出了二维管流和二维空穴流动在不同状态下的计算结果, 并与其它解法的计算结果或精确解进行比较。结果表明, TVD格式不仅能用于不可压流数值计算, 而且还能得到比非TVD差分法精度更高的结果。      

固体火箭发动机喷管气固两相流动的数值模拟

对颗粒相采用颗粒轨道模型,气相求解可压缩N-S方程组,计算方法采用显式Runge-Kutta时间推进法与有总变差衰减(TVD)性质的高精度MUSCL-Roe格式;自主开发了曲线坐标系下二维轴对称可压缩N-S方程组的解算器Solve2D,研究了固体火箭发动机喷管中颗粒相对流场的影响以及不同尺寸颗粒运动规律.结果表明:颗粒相对流场的影响主要表现在喷管喉部以及扩张段,和单相流场相比,沿轴线马赫数减小,且颗粒尺寸越小减少得越多;沿轴线气相温度升高,且颗粒尺寸越小温度升高越多;颗粒尺寸越小,无粒子区越小;颗粒越大与收缩段壁面碰撞越剧烈,无粒子区越大.      

一种求解透平叶栅三维流场的高精度TVD格式

本文提出了一种应用三阶精度Chakravarthy-OsherTVD格式求解环形透平叶栅三维定常流场的隐式近似因子分解方法。控制方程为柱座标速度分量表示的任意曲线座标系下的三维欧拉方程。采用Beam-Warming近似分解, 对隐式项进行迎风差分, 得到了时间精度为一阶、空间精度为5点三阶的隐式格式。算例对5点三阶和5点二阶精度的格式进行了比较。表明空间三阶精度隐格式对计算结果没有明显改进, 但收敛残差有较大降低。      

用CSCM格式求解激波反射问题

本文对Euler方程采用CSCM格式,数值求解了激波反射问题,与目前流行的交替方向隐式(ADI)格式以及总变差递减(TVD)格式的计算结果进行了比较,并作了简要分析讨论。     

一种修正的Osher通量在高阶WCNS中的特性

将一种基于多维修正的Osher通量运用于高阶加权紧致非线性格式(WCNS)中,该修正通量主要在垂直于激波面的界面上增加耗散,能够改善Osher通量的激波捕捉稳定性,同时对边界层和接触间断的分辨率影响非常小。对修正的Osher通量在高阶WCNS中的特性进行研究,通过数值模拟测试了基于Osher通量的WCNS的激波稳定性、热流预测精度、边界层模拟能力、激波边界层干扰模拟能力,并与Steger-Warming通量和Roe通量进行了对比。结果表明修正后的Osher通量比Harten修正的Roe通量具有更好的激波捕捉鲁棒性,而边界层、驻点热流值和激波边界层干扰的模拟则明显优于Steger-Warming通量。上述结果说明了基于修正的Osher通量的高阶WCNS具有较好的激波捕捉特性、热流预测精度和边界层计算能力。     

高超声速电离绕流的数值模拟

本文从三维N－S方程出发，采用稳式L－U分解算法和7组分15个反应的化学模型，数值模拟高超声速电离空气绕流。首先采用对称TVD格式、AUSMPW＋格式和Van Leer的矢通量分裂格式计算了高超声速球头绕流，并对它们的计算结果做了对比分析。然后用前两种格式，对RAM－C飞行试验模型三个再入高度（81km、71km、61km）的流场进行了数值模拟，计算的流场电子数密度值和试验测量数据符合较好。