粘性跨音速喷管流场计算

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

旋转固体火箭发动机喷管热结构分析

以某固体火箭发动机喷管为例,研究其旋转情况下的流动、传热及热结构响应性能。利用雷诺应力湍流模型,结合增强型壁面函数,求解N-S方程。以喷管内壁面温度分布为边界条件,求解二维轴对称瞬态热传导方程。将瞬态温度场按时间步长加载,进行瞬态静力学分析,得到不同转速下的流场、温度场及应力场。结果表明,喷管在旋转情况下,喷管内壁面尤其是喉部及扩张段温度分布在转速为100 r/min及600 r/min左右时变化剧烈;旋转内流场与喷管结构的耦合作用加剧了喷管的传热,尤其是喉衬的烧蚀,特征点温度值随转速增大而升高;最大热应力位于喷管最外层尾端,整体热应力在转速低于100 r/min时得到释放,随着转速的不断增大,喷管整体最大热应力及扩张段特征点应力随之增大,而喉衬特征点由于旋转导致了温度梯度降低,其应力值随之减小。     

轴对称Boltzmann模型方程统一算法与空天飞行环境喷管流动模拟

拦截机动飞行器周围大范围区域存在主喷/侧喷流/羽流影响,而传统的纳维-斯托克斯（Navier-Stokes, N-S）方程不能很好模拟发动机喷管扩张段出口附近流动情况,需要一种新的方法来处理这种全流域流动问题。为解决该问题,针对特定轴对称喷管内流动,本文通过数学推导确立描述不同克努森数稀薄环境条件下的轴对称喷管内流动Boltzmann模型方程,初步建立适于该模型方程的数值格式与气体动理论统一算法。通过开展同轴圆筒间的定常/非定常旋转流动以及轴对称喷管内流动数值计算研究,发现统一算法计算流场与其他途径得到的结果吻合较好,验证了统一算法在全局克努森数喷管流动模拟的适应性和可靠性。通过与低密度风洞实验对比,喷管出口核心区羽流结构一致,羽流轴线压力分布一致,表明统一算法可以有效解决喷管入口压缩段到扩张段多流域混合,尤其是出口附近稀薄气体真空低压环境流动问题。     

二维受扰超声速喷管流的动态分析

介绍了STAR -CD软件包所采用的控制方程、湍流模型和数值模拟技术。模拟了有出口旋转挡板的拉瓦尔喷管的超声速流场 ,分析了出口旋转挡板转动对超声速喷管流场速度、压力和Ma等参数的动态影响     

旋转爆震发动机喷管非定常流动特性的数值模拟研究

为分析旋转爆震喷管内非定常流动特性,基于特征线方法和旋转爆震燃烧室出口时均参数设计了一种塞式喷管模型,通过改变塞锥式喷管的进出口压比π(π=15,30,45),分析了喷管内非定常流场结构与喷管工作性能之间的相互关系,并探究了喷管进出口压比变化对喷管工作性能的影响.研究结果表明:旋转爆震燃烧室下游斜激波进入喷管后,在喷管...     

两相喷管流动研究进展

本文综述了两相喷管流动研究的进展,着重介绍了两相喷管流动控制方程的数值解以及两相喷管流动研究的实际应用。     

固液混合火箭发动机喷管流动计算

利用二维轴对称 N-S方程和组分方程对选用液氧 /端羟基聚丁二烯推进剂的固液混合火箭发动机喷管流动进行了计算。计算采用 LU时间隐式格式、MUSCL空间离散和 Van Leer矢通量分裂方法 ,8组分 1 0反应的化学反应模型 ,对化学源相进行了点隐式处理。喷管入口条件通过燃速计算和热力计算得到 ,计算了多个氧化剂流率和装药长度情况下的喷管流场 ,分析了真空比冲和推力与氧化剂流率的关系 ,为固液混合火箭发动机喷管设计提供了依据      

旋转对固体火箭发动机燃烧室燃气流动的影响

采用RNG k-ε湍流模型,计算了某管状装药旋转固体火箭发动机在不同旋转速度情况下的燃烧室喷管统一流场,着重分析了不同转速情况下燃烧室内燃气流动的特点,得到了该发动机燃烧室内旋流切向速度沿径向分布趋势以及旋流涡核半径在不同转速情况下的变化特点,对燃烧室前封头局部区域的流场结构也进行了较为细致的分析.     

旋转飞行器固体火箭发动机内涡核流动的数值分析

为分析旋转飞行器固体火箭发动机内涡核流动,采用数值计算的方法对发动机流场进行模拟。计算结果显示发动机的旋转增强了轴向涡核流动,随着涡核的不断延伸,在喷管入口处,由于潜入式喷管结构以及喷喉直径较小,振荡运动随之加强。涡核在流经喷管时的不断增强,在喷管区域形成了非对称的压力形式,并且产生了实际中的扰动侧向力,最终导致了飞行器的章动不稳定性。     

三轴承旋转喷管矢量偏转规律及流场特性研究

基于短距/垂直起降战斗机用三轴承旋转喷管的特殊设计要求,通过几何约束条件开展了三轴承旋转喷管型面设计方法及运动规律推导与研究,研究了非线性和线性两种喷管矢量角控制规律下的三段筒体随时间的旋转规律。根据小型涡轮喷气发动机的几何尺寸,利用发展的型面设计方法和喷管筒体旋转规律,设计了小尺寸三轴承旋转喷管,并利用CFD数值模拟技术对该喷管的流场特性进行了计算与分析。通过CFD数值模拟技术得到了不同矢量角下喷管的三维流动特性及不同落压比下的气动特性。结果表明:采用喷管矢量角非线性控制规律可以减少非线性控制变量,保证喷管机动性的前提下减小了喷管的设计难度和控制复杂度;基于小型涡轮喷气发动机设计的三轴承旋转喷管0°的推力系数较理想喷管低,90°的推力系数较理想喷管高,喷管在地面最大落压比下0°比90°推力系数高约1%。     

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

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

不同矢量偏角下二元喷管内流特性的数值研究

应用多重网格MacCormack显式格式对不同矢量角下二元喷管的内部流动进行了模拟分析。计算时, 在细网格上采用了雷诺平均N-S方程与k-ε湍流模型方程;在粗网格上则采用Euler方程。数值计算给出了较清晰的流谱, 壁面压力分布的计算也与实验结果相吻合。通过计算结果分析了矢量喷管几何参数对流场的影响, 并对矢量喷管的推力特性进行了数值研究。      

高焓激波风洞喷管流场非平衡特性研究

高焓激波风洞是开展高超声速流动研究的重要地面模拟设备,但其产生的高焓气流在喷管中的膨胀过程是一种典型的热化学非平衡流动,试验段特征参数通过直接实验测量难以完全确定。本文通过求解耦合双温度模型的轴对称Navier-Stokes方程,研究了高焓激波风洞中典型状态下气流的热化学非平衡流动特性,分析了焓值对非平衡特性的影响规律。结果表明,喷管出口自由流均匀区域达到出口截面直径的75%以上,能够为实验提供足够的空间;喷管出口自由流处于热化学非平衡状态,在喷管喉道后约1/5喷管长度处气流即已处于冻结流状态,组分浓度和振动温度随气流流动基本不变;焓值在8.4MJ/kg~19.5MJ/kg之间变化时,非平衡程度随着焓值的增加而增强,但是低焓值时非平衡程度的增强更加剧烈。     

三维数值模拟再生冷却喷管的换热

为了解液体火箭发动机喷管再生冷却的换热特点，采用数值模拟的方法，对内喷管燃气、壁面和冷却液建立不同的三维控制方程，进行流动和传热的耦合计算。在计算中，假定喷管流动为冻结流动，考虑燃气向壁面的对流换热和辐射换热；采用二阶迎风格式离散控制方程，采用DO模型离散求解辐射换热方程，水蒸气的吸收系数根据Leckner公式计算。计算模型采用缩比热试车发动机，数值计算结果与实验结果吻合较好，较准确地模拟出了喷管的壁面热流密度，得到了喷管燃气和冷却液的流场和温度场，对高压再生冷却喷管的设计具有指导意义。     

喷管加装小突片强化射流混合的PIV测量

在常温、超声速流动状态下,采用PIV技术定量测量了收敛喷管和收扩喷管加装小突片时射流纵截面的速度场和涡量场,并对其速度场和涡量场进行了比较。结果表明,对于收扩喷管,虽然其速度边界层厚度明显大于收敛喷管,但加装小突片结构也能够明显的强化射流的混合。并且发现,小突片结构不仅能产生流向涡,而且对周向涡有明显的强化作用。     

任意二元和轴对称细长喷管可压缩有势流动

本文应用摄动法给出了一种简便的计算任意二元和轴对称细长喷管可压缩有势流动的近似解析方法。流场的解是相对于一个表示喷管细长程度的小参数进行渐近展开,而展开式的首项就取喷管当地截面处的一元流的解。对于速度分量,文中给出了展开式中的前两项的表达式。数值算例的结果表明,当喷管壁的斜率大约不超过20°,喷管壁的曲率半径大约不小于喷管截面半高度的4～5倍时,用本方法计算的速度的模的相对误差不超过1％。本方法可用于计算Laval喷管内由亚音速加速到超音速时的流场。     

旋转通道一维非定常计算模型

采用特征线法,对航空发动机二次空气系统的旋转通道结构一维非定常流动及换热的计算模型进行了研究,通过将广义的定截面、有摩擦、考虑对外热交换的旋转通道可压缩非定常流动偏微分方程组转换为常微分方程组,考虑旋转离心力,并采用Colebrook White方程定义其阻力特性,进而获得旋转通道内部沿流动方向压力、质量流量、温度等系统参数的特征线数值解,采用商业软件Fluent对其进行验证.结果表明:特征线法利用压力波传播方式对旋转通道进行求解,能够直观地描述旋转通道一维非定常流动及换热特性,计算误差不超过5%,可以作为二次空气系统旋转通道一维非定常计算研究的有效手段.      

高超声速喷管三维热化学非平衡流场的数值模拟

利用数值求解三维热化学非平衡Navier-Stokes方程的方法和多块网格技术,数值模拟了典型高焓激波风洞锥型喷管的三维非平衡流场,分析了热化学非平衡效应对喷管流动的影响,给出了平动温度、组分N2和O2振动温度、马赫数和组分质量分数在轴对称线上的分布,也给出了喷管出口附近的速度和温度沿径向分布等结果,获得了喷管流场的详细信息,并与已经存在的采用轴对称Navier-Stokes方程计算的结果,进行了比较,二者吻合得很好.计算结果表明:目前的代码是可以模拟多块结构网格下的热化学非平衡流动的.     

一维摩擦喷管流动

本文阐述了一维摩擦喷管流动控制方程的数值求解方法,分析研究了摩擦对喷管流场及发动机性能的影响。确定边界条件时考虑了燃烧室工作过程,并给出了计算例题。     

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

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