用隐式近似因子分解法计算喷管跨音速流场

本文用隐式近似因子分解法结合贴体曲线坐标计算了定常、无粘、跨音速喷管流场,选择了参数变化比较剧烈的小喉部曲率半径喷管,陡壁收敛段喷管和潜入喷管三种算例.计算结果表明,该方法收敛快,并且具有良好的精度.Courant数可以取至7.这种隐式法与显式法相比较,花费的机时少得多.本文的方法可望推广到求解N-S方程、PNS方程,和二相有粘喷管流动计算.     

塞式喷管流场计算

采用NND差分格式，对塞式喷管轴对称流场进行了数值模拟。通过流场分析，可以对塞式喷管火箭发动机有清楚的认识。计算表明，NND格式有很好的分辨率；无粘和有粘的情况流场结构和性能相差悬珠；塞式喷管的自适应能力很好。     

粘性跨音速喷管流场计算

本文继用隐式近似因子分解法成功地计算无粘跨音速喷管流场之后,用同样的方法,结合任意曲线坐标系,通过非定常方程在相当长时间后的定常解,针对矩形截面喷管和轴对称喷管粘性跨音速流场,求解了可压缩层流薄层N-S方程,获得了核心流,特别是边界层中的流动参数.对于矩形截面和轴对称两种喷管进行了计算,其结果和实验数据相当一致.将本文的方法应用到粘性两相流动计算,可望较多的节省机时.     

座舱供气喷管计算

利用理论计算方法,对喷管在不同压力下的流量进行计算,得到喷管的流量—阻力特性曲线,并与试验曲线进行对比,验证了理论计算的准确性。该方法可以用来进行环控系统座舱供气喷管设计。     

喷管跨声速流场计算

本文用时间相关法求解定常跨声速喷管流场。对守恒型和非守恒型方程采用全场统一时间步长、当地时间步长和修正超阻尼法进行了计算。计算表明:采用当地时间步长可大大提高收敛速度,采用全场统一时间步长时,数值解非守恒型方程收敛较快;采用修正超阻尼法也能提高收敛速度,但与阻尼系数和时间间隔有关。     

塞式喷管与当量钟型喷管性能的实验测定和分析

为了合理地评价塞式喷管是否具有高度补偿特性，利用空气作为工作介持，对实验塞式喷管和钟型喷管进行了冷流比较实验研究。三种不同的喷管效率计算方法对两种喷管的实验结果进行了处理，得到不同的性能曲线，分析比较两种喷管在不同高度下的效率。与当量实验钟型喷管相比，实验塞式喷管具有高度补偿特性，算法一和算法二分别比较了喷管性能的不同方面，评价喷管性能时可以综合考虑。     

不同环境压强下塞式喷管流场的数值模拟

本文设计了一种环排式塞式喷管,通过求解N-S雷诺平均方程组,对在不同环境压强情况下该塞式喷管的流场结构进行了数值模拟,对采用单方程和双方程两种紊流模型模拟的结果进行了比较。数值模拟的结果表明,用适当选取混合长度的单方程模型与双方程模型计算的结果基本一致;该塞式喷管基本上体现了塞式喷管的优点。      

流体注入的轴对称矢量喷管三维流场计算

采用Roe通量差分分裂格式对基于流体注入控制的轴对称矢量喷管内流场进行了数值模拟。流体注入的位置分别为前孔和后孔,注气压强比为0.75～2.0,注气流量比为2.5％～10.0％,矢量角变化范围为2.8°～7.8°。计算结果表明:随着注气流量和注气压强增加,流体注入所产生的喷管矢量角相应增加;注气位置对喷管矢量角影响较大,注气位置靠近喷管尾沿(后孔注气)比注气位置靠近喷管喉部(前孔注气)所产生的矢量角明显增大。     

斜切反喷管性能计算

本文提出了斜切反喷管流场计算方法,内部点用MacCormack二步显格式数值解控制方程,入口边界和固壁边界的参数用物理边界条件和可应用的沿双特征线的相容性方程计算.固壁边界计算中采用了激波装配技术和开关格式,有效地捕获了流场内的激波系列,计算中得到的激波情况和物理分析一致.本文在流场计算的基础上,进行了反喷管的性能计算.本文提出的方法可应用于反喷管的性能预估.     

变比热容对喷管流场和性能影响的计算

从变比热容角度出发,采用ＣＦＤ方法计算了发动机加力工作时高温状态下喷管内部流场分布和性能,并和定比热容计算结果进行了比较。结果表明,采用ＣＦＤ方法计算加力工作时喷管内部流场和性能,应该考虑变比热容的影响。     

二元矢量喷管流场的N－S方程计算

对二元矢量喷管在不同工况下的流场进行了数值模拟。求解的控制方程为三维非定常Ｎ－Ｓ方程。空间方向用有限体积离散，时间方向用Ｒｕｎｇｅ－Ｋｕｔｔａ法推进。粘性项的计算上采用了一种有限体积框架下新的离散化方法。湍流模型采用了广泛使用的Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数模型。计算表明，本文方法能准确地计算二元矢量喷管、非矢量喷管在设计与大设计落压比下的流场。     

超音速喷管收缩扩张段匹配问题研究

超音速喷管由收缩和扩张两段组成,是高速风洞中获得超音速流动的重要部件。对于超音速喷管收缩段设计方法的研究较多,但尚未见到喷管收缩段的几种设计型面曲线如何与扩张段合理匹配问题的研究。针对以上的问题,利用CFD数值模拟方法对几种收缩曲线与扩张段匹配问题进行了计算。通过对计算结果进行对比分析,发现不同收缩曲线对超音速喷管内流有一定的影响,分析了流场特征,展示了马赫数在喷管内部的分布。     

化学非平衡喷管流计算方法的研究

本文提供了考虑化学反应非平衡和振动非平衡的喷管流动计算方法。在忽略扩散、粘性和热传导及辐射损失的假设下，得到完全气体一维定常流动方程，再加上振动能松驰方程和化学反应速率方程即可数值积分求解。采用扰动法得到接近平衡流的一组积分初值。采用“逆积分方法”求解通过喉道。计算结果是令人满意的。     

复杂几何形状喷管内外三维流场的数值模拟

研究了球面收敛矩形扩张的二元喷管（ＳＣＦＮ）流场的数值模拟,为了有效地捕获激波与滑流面及较好地模拟分离,采用了Ｒｏｅ格式、ＭＵＳＣＬ方法及Ｖｅｎｄａｔａｋｒｉｓｈｎａｎ限制器、Ｓｐａｌａｒｔ－Ａｌｌｍａｒａｓ方程紊流模型。研究中给出了两个喷管模型四个流动状态的数值模拟结果。结果表明,除了有严重分离的流动情况以外,该算法能得到相当满意的计算结果。     

高精度ENO格式在喷管流动模拟中的应用

本文分析了ＥＮＯ格式的特点,并应用于Ｅｕｌｅｒ方程和全ＮＳ方程的迁移项和压力项,模拟了各种喷管流动。首先计算了ＪＰＬ喷管流动,取得了和实验一致的结果,分析了喷管内产生弱激波的原因;其次计算了非定常二维喷管流动,给出了相应的计算结果,并分析了壁面附近流场的拓扑结构;最后计算了摆动喷管的流动及推力矢量计算,取得了推力随摆角变化的规律,和实验结果表明基本一致。     

轴对称矢量喷管三维传热计算研究

采用封闭腔理论算法，利用矢量喷管的三维内流场数值计算结果，对轴对称矢量喷管的壁温分布进行了数值计算。计算结果表明，在非矢量状态，喷管同一轴向位置处的周向壁温基本相同，壁温分布可按一维处理。在矢量状态，喷管同一轴向位置处的周向壁温相对差值达15％，壁温分布呈现复杂的形式，必须采用三维计算方法计算壁温分布；喷管扩张段壁温明显高于非矢量状态，所开发的计算方法和程序既可用于轴对称矢量喷管的工程设计，也可为试验验证提供理论依据。     

轴对称喷管流场分析的有限体积法

本应用有限体积法对轴对称喷管的跨音速流场进行了分析，采用显式预校两步MacCormack格式进行计算。运用隐式残差平均技术及当地时间步长来加速计算的收敛，提高格式的稳定性。对三个不同几何参数的轴对称喷管的数值分析表明，计算结果与试验数据吻合良好。     

固体火箭发动机内流场分区耦合数值计算

采用分区耦合方法计算固体火箭轴对称燃烧室与喷管流场。对于低速的燃烧室流场，选用不可压流的Ｎ－Ｓ方程描述并用ＳＩＭＰＬＥＣ方法数值求解；对于高Ｒｅ数的喷管流场，则采用Ｅｕ－ｌｅｒ方程描述并用ＳＣＭ方法求解。计算时用燃烧室出流为喷管流场提供入口参数，同时用喷管流场压强分布反馈影响燃烧室流动状况。对耦合边界条件处理方法进行了探讨。对典型的侧壁加质燃烧室与喷管流场进行了计算，计算结果揭示了单独喷管流场计算难以反映的喷管收敛段近壁区的低速区域，与已有的燃烧室流场实验结果一致并反映了燃烧室与喷管流场之间的联系，较好地模拟了流动中的物理现象。     

固体火箭喷管两相粘性跨音速流场计算

本文进行了固体火箭喷管两相粘性湍流跨音速流场计算.粘性的气相控制方程用隐式近似因子分解法求解,粒子方程采用跟踪粒子轨迹的特征线法求解,粘性湍流选用代数模型,气相和凝相充分地偶合.CFL数可以取500左右,收敛速度快,使粘性两相跨音速喷管流场计算耗费的机时达到工程计算可以接受的程度.通过计算,获得了在粘性和粒子同时作用下的流场参数分布.这对固体火箭发动机流场需要同时考虑粒子和粘性作用的专题研究大有裨益.     

潜入喷管背壁流场模拟显示

为了研究固体火箭发动机潜入喷管背壁流场的特点,设计了一种冷模拟实验装置.利用引入外加烟雾以及片光照明的显示技术进行背壁区域的流动显示,并用照相机和摄象机记录显示结果.结果表明:在雷诺数较高、潜入度较大和喷管潜入部分型面为凹弧面条件下,背壁区域容易形成旋涡;同时,针对实验条件进行内流场的数值计算,实验和计算结果相比较,二者在定性上具有良好的一致性.