固体推进剂非稳态燃烧模型数值特性

在修改的KTSS燃烧模型的基础上,研究了热波燃烧模型的数值特性。在大幅度改变压力变化速率参数（β）和热释放参数（H）下,计算结果显示：1）在大β值时,燃速并未发散溢出,直到β=∞,燃速仍然是一个有限值,在一个峰值波动之后,趋向于稳态燃速。2）在增大H值时,燃速的不稳定是逐渐增加的。当H＜1时,燃速振荡后最终回到稳态值；当1＜H＜1．045时,燃速是一系列重复、有限的尖峰；但当H≥1．045时,燃速发散。     

超声速气流中轴向涡强化混合及燃烧的数值研究

通过求解N－S方程，数值计算了扩张圆管内的非平衡流动，并研究了超声速混合和燃烧的强化技术。结果表明，在超声速领域，预旋燃料方式是一种比较有效的强化技术。     

凝相相变过程的热力计算

一般的热力计算软件很少计及凝相物质的相变过程，当燃烧产物温度接近某凝相物质的相变温度时，将出现迭代计算不收敛。根据能量守恒、喉部单位质量流量的流通面积最小和等熵关系，提出了凝相物质在燃烧室、喷喉和喷管出口处于相变过程的热力计算方法，并给出了算例。     

布雷顿和卡诺热机循环性能比较(Ⅱ)定常态流变温热源循环

研究定常态变温热源热机循环性能，导出内可逆卡诺热机和布雷顿热机的最佳功率、效率关系和最大功率及相应的效率界限，并对这两种热机循环的最优性能进行了比较。理论分析表明，在相同的边界条件和热效率下，布雷顿循环的功率可以高于卡诺循环的功率，极限情况下前者是后者的两倍。对于变温热源条件，布雷顿循环主要受益于其工质与热源间有较佳的匹配。所得结果对热机工作参数和工质的最优选择有一定指导意义。     

轴流压气机子午流面计算中三维剪切作用的湍流模型

本文发展了一个轴流压气机子午流面二维计算中模拟三维剪切结构对子午流场作用的湍流模型。试验与这一模型进行了初步的比较。模型表明,在多级轴流压气机子午流场计算中,必须计入实际流场三维剪切结构对子午流场的湍流扩散作用。     

等速上仰翼型绕流结构的数值模拟

通过数值求解涡量-流函数形式的N-S方程研究等速上仰翼型的分离流动结构。用ADI方法解涡量方程;用Poisson方程直接法解流函数方程。在计算得到的流场中除可看到实验中观察到的多涡结构外,还发现一些实验中未观察到的现象。将分离流动结构的计算结果与翼型气动力变化相联系,并结合涡动力学理论,定性地解释了上仰翼型产生高升力的原因。     

转动坐标系中三维跨声速欧拉流的有限体积 TVD格式

在非惯性转动坐标系中,本文采用贴体网格、有限体积法离散和修正数值通量技术,将Harten的一维TVD格式推广到三维。由于转动使方程出现源项,文中通过对源项的巧妙处理,使修改后的格式能用于非齐次双曲守恒律方程组高分辨率的数值计算。为了加速解的收敛,提高显式时间推进的CFL数,本文采用隐式残值光顺技术。三维跨声速带非齐次源项欧拉方程的典型算例表明:捕捉的激波分辨率较高;激波前后没有发现大的数值波动和伪振荡现象;所得的跨声速流场解与实验较接近。     

冷热两股气流在圆管内掺混的速度场、温度场与辐射通量场的数值方法

用数值方法求解冷热两股轴对称平行射流在圆管内掺混的流场、温度场和辐射通量场。计算采用一种通用性的数值方法,其特点是强烈地以物理本质为基础,而不仅仅是单纯的数学推导。这样做,能方便地分析和解释计算结果,使之更真实地接近于实际。 通过计算,分析了当引射系数改变时,掺混的流场,温度场和辐射通量场的变化规律,其结果与实验数据比较一致。     

椭圆截面喷管排气辐射热流数值分析

文中对直升机红外抑制器椭圆截面喷管排气的辐射热流进行数值分析,并预估排气温度分布、速度和燃气成分分布对辐射热流的影响。为了计算辐射热流和考虑尾焰对周围环境的传热,采用形式简单并可用通用方程计算程序计算的六热流辐射模型。由于排气流场内温度和密度变化较大,故紊流模型采用考虑密度变化的修正双方程k-ε模型。本文中所用的比热C_(pm)是随气流温度和成分而变化,温度用迭代法求得。喷管出口截面为椭圆形,排气为三元自由射流,控制方程为三元抛物形偏微分方程,可采用前进积分求解,这样既可满足计算精度要求,又可节省机时。并根据三元喷流周向边界条件周期性变化的特点,采用CTDMA解法。计算结果与试验结果基本一致,说明本计算方法是可行的。     

飞行器完全非定常气动力的一种工程计算方法

本文以谐振运动非定常气动力计算理论为基础,采用偶极子格网技术进行数值计算、有理式拟合、Fourier积分和Duhamel积分等步骤,给出了飞行器任意运动和受任意大气扰动时完全非定常气动力的一种工程计算方法。该方法数学处理方便,得到的气动力表达式简单。算例结果表明:该完全非定常气动力计算方法可行,精度令人满意,具有工程实际意义。