固体微推力器设计与数值分析

设计了两种基于微光机电技术的固体微推力器结构,其中一种结构具有不连续的非正常拉瓦尔喷管。采用数值计算方法对不同的设计进行了内弹道、燃气流动过程和热损失的数值分析。结果表明:不连续的非正常拉瓦尔喷管同样能够将推进剂燃气加速到超声速,其比冲性能与正常拉瓦尔喷管相当;微光机电技术常用的硅材料会造成推力器比冲性能较大量级的下降,而绝热性能良好的玻璃或陶瓷材料可以有效提高微推力器的比冲性能。     

环形燃烧室性能计算

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

97．5．8华南地区暴雨过程的数值模拟

用非静力平衡的中尺度数值模式MM5V3版，模拟研究了1997年5月8日OOUTC至5月9日OOUTC，发生在华南沿海地区的暴雨过程。结果表明：该模式比较成功地模拟预报了产生该次暴雨的环流形势、急流位置及其强度、暴雨的范围、雨量中心；通过分析数值试验的结果表明：凝结潜热释放的加热作用对暴雨产生和华南低空急流的维持是十分重要的。     

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

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

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

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

腰形柱壳的有限条元解

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

求解大型对称特征值问题的迭代Chebyshev-Lanczos方法

本文研究计算大型对称矩阵极端(几个最大或最小)特征值及相应特征向量的问题,讨论了Chebyshev迭代法对Lanczos方法的应用,提出了Chebyshev-Lanczos方法。计算实践表明迭代Chebyshev-Lanczos方法比迭代Lanczos方法优越。     

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

本文用非定常数值法研究了钝体燃烧器后的紊流参数分布和能量逆转,并对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的高超声速流中压缩角引起的激波／湍流附面层干扰进行了数值研究。计算结果与实验结果吻合较好。