环形燃烧室性能计算

在任意曲线坐标系下对包括扩压器和火焰筒在内的环形燃烧室三维两相反应流进行了数值模拟。燃烧室性能采用多维经验分析法估算所用的数学模型有 κ-ε,紊流模型 ,EBU -Arrhenius紊流燃烧室模型 ,六通量热辐射模型以及颗粒群轨道模型。在非交错网格体系下 ,气相采用 SIMPLE算法和液相采用 PSIC算法求解 ,通过两种工况和三种扩压器进口速度分布计算表明 ,所得的流场各气流参数分布和燃烧室性能较合理 ,本文所建立的程序可为燃烧室优化设计和研制提供有用数据。     

典型几何和流动参数对高超声速进气道性能的影响

用N-S方程模拟了一系列典型二元高超声速进气道内压缩通道及隔离段模型,模拟发现内压缩通道及隔离段增压比、温升比和总压恢复系数等性能参数主要受面积收缩比、内压缩通道收缩角、隔离段长高比等几何参数以及内压缩通道进口马赫数、密度、附面层厚度等流动参数的影响。内收缩比和内收缩通道收缩角的增大都会使压缩增强;隔离段内沿程平均温升比、马赫数和总压恢复系数曲线则均接近平行直线。内压缩通道进口马赫数的增大也会使压缩增强,但较小的进口马赫数可能引起分离,进而增大增压比;而进口密度增大使附面层变薄,对气流的压缩减弱;进口附面层厚度对沿程平均温升比、增压比以及马赫数的影响近似线性。     

高速风洞支架干扰数值修正研究

提出了一种计算高速风洞支架系统对飞行器模型纵向气动力干扰量的数值计算方法 ,从跨声速全位势积分方程出发 ,编制了适用于飞行器全机模型及其带支架情况下的跨声速绕流计算程序。通过对双垂尾模型和GBM 0 3模型两个算例的计算 ,讨论了尾支撑位置及其几何外形参数对模型气动力的影响 ,并对GBM 0 3模型带短支杆情况下的纵向实验结果进行了修正。表明该方法对于分析研究风洞模型支架干扰问题并进行支架干扰修正是可行的、有效的 ,可以作为选择尾支撑位置及其几何外形参数和对跨声速风洞纵向实验结果进行支架干扰修正的工具。     

腰形柱壳的有限条元解

本文用考虑了曲率影响的两结线低阶有限板条元分析了腰形柱壳的力学性能这种条元位移参数少,能较好地体现壳体的几何特性。     

钝体燃烧器后能量逆转和紊流参数的数学模拟

本文用非定常数值法研究了钝体燃烧器后的紊流参数分布和能量逆转,并对Boussinesq的雷诺应力模式和双方程k-ε紊流模型提出了改进方案。这样不仅计算出了能量逆转区,而且紊流参数分布也与实验相符。     

关于Davidson—Lanczos方法的收敛率

本文对文[1]提出的求解大型对称矩阵A的极端(几个最大或最小)特征值及相应特征向量的Davidson—Lanczos方法,用Rayleigh—Ritz逼近理论,研究了该方法的收敛率。证明了由该方法产生的规范正交向量{v_i}_i~m=1是Krylov子空间K_m≡Span(v_1,Av_1,…,A~(m-1)v_1)的一组基。设A的k个最大特征值为又,λ_1>λ_2>…>λ_k,相应的近似特征值为λ_i~(m)(i=1,…,k),得到 这里γ_i(γ_i>1),W_i和W_i~(m)是常数。     

超声速流中激波/湍流附面层干扰数值模拟

采用修正的B／L湍流模型以及多块结构化网格求解了二维N-S方程。分别对超声速流和高超声速流中的激波／湍流附面层干扰进行了数值研究。本文首先研究了进口马赫数为2．96的超声速流。计算结果准确预测了入射斜激波在平直壁面引起湍流附面层分离的流动特征：分离点的反射激波、分离包引起的膨胀扇以及再附点的反射激波。计算的壁面压力分布与实验值吻合较好，计算的分离区长度与实验值比较有一定误差。本文还对进口马赫数为9．22的高超声速流中压缩角引起的激波／湍流附面层干扰进行了数值研究。计算结果与实验结果吻合较好。     

火烧炉内气相燃烧过程的三维数值模拟

本文对某工业火烧炉内的三维湍流气相燃烧场特性进行了数值模拟,湍流模型采用可实现k-ε双方程模型,湍流燃烧采用有限速率/涡耗散模型,辐射换热采用P1模型,数值方法采用SIMPLE算法,实现了炉内火烧试验的数值再现,特别考察了入口压力、燃料/空气比、多烧嘴不同组合工况等因素对炉内温度场和流场的影响,并对该火烧炉的运行提出了一些合理化建议,为其实际运行提供了有益的参考。     

对旋轴流风机非定常流场数值模拟

利用全隐格式的非稳态压力修正求解方法,分析和研究了对旋叶轮级间流场非定常流动情况,上游尾迹与下游叶片位势作用对流场的影响,并与动态实验结果进行了有益的对比.结果表明:非定常流场的计算更能反映叶轮机实际运行性能,如果能很好地组织对旋叶轮的非定常流动,还可以进一步提高其压比和效率.     

人工释能法求解线性代数方程组

本文建立求解线性代数方程组的人工释能法迭代格式,并给出该迭代格式收敛的充分必要条件和最佳释能时机,同时给出实现人工释能法的数值计算方法及其数值稳定性条件。