多重网格法求二元跨声速流在一种半自适应网格上的欧拉解

本文采用MacCormack两步显式格式,用有限体积法求解了二元跨声速欧拉流。推导了物面边界条件,采用了特征远场边界条件及远场环量修正。利用保角变换方法生成O型贴体网格,并修改得到了一种在激波处局部加密的半自适应贴体网格。采用多重网格及焓阻尼加速收敛技术计算了NACA0012翼型的跨声速气动特性,得到了十分满意的收敛过程和计算结果。     

建筑物风环境三维定常数值仿真

针对建筑物的风环境问题，采用计算流体动力学CFD软件对两幢建筑物进行了数值模拟，得出了建筑物的风场分布，以及风载的分布规律，结果表明，气流在来流方向第一幢建筑物的顶部、两侧的前沿产生分离，建筑物之间产生强烈的漩涡，顶部的局部负压最大。     

环形燃烧室性能计算

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

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

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

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

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

腰形柱壳的有限条元解

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

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

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

一种采用解析积分的直线单元的边界元法

本文提出一种采用解析的直线单元计算受有离心力的二维问题的边界元法,推导了相应直线单元的积分解析式。算例表明,本文方法能获得较高的数值精度,尤其在计算靠近边界的内点应力时,效果更好。     

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

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

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

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