喷管化学反应粘性流场的数值模拟

用时间相关法的显格式求解氢氧发动机喷管的粘性化学反应流动,采用６种组分,８个基元反应有限速率的化学反应模型,湍流模型采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,得到流动参数和产物的质量分数在喷管流场中的分布。计算结果表明,采用显格式方法求解喷管化学反应粘性流动是可行的。     

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

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

液体火箭发动机推力室粘性流场数值模拟和实验验证

运用耦合点隐式方法的MacComak两步差分格式结合κ-ε湍流模型求解肼类燃料发动机推力室中的粘性化学反应流动。化学反应采用17种组分，12个主要反应的有限速率的化学反应模型，得到了流动参数在推力室中的分布情况。其结果同理论分析的相一致，最后用实验验证了用数值模拟方法的正确性，研究结果对发动机的设计、性能和理论研究具有重要的实际意义和参考价值。     

C-H-O-N体系液体火箭发动机喷管流场数值模拟

采用弱耦合点隐式方法的MacCormack格式对C－H－O－N体系的液体火箭发动机喷管粘性流场进行了数值模拟，采用了考虑化学反应和不考虑化学反应的两种模型，其中化学动力学模型为C－H－O－N体系的12组分、14反应有限速率的化学反应模型。数值模拟得到了流场参数在喷管中的分布，结果表明，数值模拟结果与理论分析一致，为发动机的设计及试验提供了参考。     

跨音速喷管化学反应流数值模拟

用时间相关的半隐格式有限差分数值方法求解了化学非平衡反应跨音速喷管流场，在喷管收敛段，流动接近化学平衡状态，控制方程的刚性问题严重，数值积分困难。通过对时间差分项隐式离散、对空间差分项显式离散，流场边界采用参考平面上的特征线法计算等，成功地解决了由于化学反应有限速率带来的数值解不稳定问题。该格式简单、需要计算机存贮空间少。本文完成了一维和轴对称非平衡化学反应喷管流动计算，并与化学平衡流和冻结流的计算结果做了比较。     

粘性跨音速喷管流场计算

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

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

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

双钟型喷管高度补偿特性的数值分析

针对双钟型喷管流场中存在的激波、激波与附面层相互作用、大尺度分离流动等复杂的物理现象，应用RNGκ-ε湍流模型封闭二维粘性可压缩N-S方程组，采用二阶迎风格式及隐式Gauss-Seidel迭代方法进行耦合求解，对双钟型喷管流场进行数值模拟研究。结果表明，双钟型喷管在低空出现壁面可控流动分离，避免了过膨胀现象，性能优于大面积比钟型喷管；在高空下为完全附着流动，性能优于小面积比钟型喷管，实现了高度补偿特性。     

用改进的 SIMPLE 方法求解喷管流动

应用ＳＩＭＰＬＥ方法求解全速度流动问题,这一方法是以速度分量、压力和温度作为原始变量,密度通过状态方程求得,压力对密度的影响隐式处理。利用此方法求解了粘性喷管流动,并与试验结果进行了比较,吻合较好。     

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

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

矢通量分裂隐式有限体积法解超音速喷流强干扰流场

发展了用于轴对称流场,任意非正交曲线网格下的矢通量分裂隐式有限体积法;统一了二维和轴对称流场计算的表达式。该方法对粘性通量项进行了分解和归类,考虑了包括交叉导数项在内的所有项对隐式增量的贡献。对带中心喷流的流动计算表明该方法收敛速度和计算精度均优于ＭａｃＣｏｒｍａｃｋ显式格式。横向喷流强干扰流场的计算初步揭示了姿控发动机喷流产生间接推力的机理,由此,设计中可降低对主发动机额定推力的要求     

斜切反喷管粘性流场数值模

为了在级间分离期间提供反推力，许多固体火箭发动机前端都装有一组斜切反喷管，由于反喷管气动型面具有尖点，并在超音速区有台阶，喷管内存在一系列激波，并伴有流动分离现象。从雷诺平均的非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程出发，利用时间相关法，采用ＭａｃＣｏｒｍａｃｋ两步显格式，结合Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型，数值模拟了斜切反喷管流场。计算得到的壁面压强分布与风洞吹风实验测得的压强分市相当一致。     

环流喷管的轴对称湍流流场计算

对火箭发动机环流喷管进行了二维数值仿真,用数值方法求解Navier－Stokes方程,数值格式为显格式时间推进格式,湍流模型选择可靠性较好的两层代数湍流模型。由于所研究的流场边界形状复杂,计算中采用了分区的方法。通过多个算例的对比,确定了使环形喷管推力效率较高的几何形状。画出了环流喷管内部的计算网格图、速度矢量图、等马赫线图、等压线图及喷管上壁面及轴线的压力分布图。     

预处理方法在化学非平衡流场数值模拟中的应用

将预处理方法发展应用于化学非平衡流场的数值模拟中。采用有限差分,LU-SGS隐式方法对二维可压化学非平衡控制方程进行了耦合求解,对化学源项和无粘项隐式处理,粘性项则采用显式。通过化学反应剪切层和激波点火的算例表明,所采用的数值方法能够有效求解化学非平衡流场,证明了本算法的可行性,为下一步的工程应用奠定了基础。     

带有分离的二维可压粘性喷管流动计算

采用二维、时变、可压纳维-斯托克斯(N-S)方程计算喷管流动。流动可为无粘、层流或湍流。湍流模型可采用代数混合长模型,一方程和二方程模型。喷管形状为二维(平面或轴对称)。计算网格内点采用MacCormack显式格式。采用局部细网格,快速解法以及时间间隔光滑方法,以加速高雷诺数下收敛。边界点采用参考平面特征线格式计算,粘性项作为源项处理。给出了初步算例结果,与其它计算和实验数据比较,符合良好。     

跨声速轴流风扇转子流场三维粘性数值分析

采用一种快速求解三维粘性流场的计算方法求解跨声速压气机风扇转子内部流场。该方法以ＬＵ型隐式格式和ＭＵＳＣＬ ＴＶＤ迎风格式为基础,结合壁面函数方法和简单的混合长度湍流模型,使用多重网格迭代加速收敛技术,对三维可压缩雷诺平均Ｎ－Ｓ方程进行求解。计算得到了ＮＡＳＡ Ｌｅｗｉｓ ６７跨声速、低展弦比轴流风扇转子的性能曲线,并重点分析了最高效率点工况和近失速工况下的内部流场。计算与实验结果的对比表明,此方     

大迎角振动三角翼跨声速粘性绕流数值分析

本文通过数值求解三维可压缩雷诺平均ＮＳ方程，研究发析了振动三角翼的三维复杂非定常流场，其中对流场项计算格式的空间离散在平行物物理面方向与物面法线方向的空间导数分别采用二阶ＮＮＤ格式和由ＮＮＤ格式高阶插值而得到的三阶迎风格式，时间推进为隐式ＬＵ分解方法，并使用了修正的代数ＢＬ湍流模型模拟高雷诺数的湍流流动。     

高超声速化学非平衡粘性激波层方程的数值解

本文讨论了高超声速粘性激波层方程数值计算时差分格式引起的物理失真问题。具体分析了全隐格式格式粘性的影响,并作了数值试验。为了验证隐式结果的可靠性。在超声速激波风洞中测量了钝锥的表面压力分布,并与计算结果作了比较,两者基本一致。 本文采用隐式有限差分法数值计算了高超声速化学非平衡粘性激波层绕细长球锥的流动。计算时采用连续方程和法向动量方程耦合求解的方法以解决细长体远后身区计算中的问题。应用网格技术和加强系数矩阵主对角元素优势的方法提高了化学非平衡流计算的雷诺数范围。文中给出了高超声速化学非平衡流的计算结果,并与其它文献的结果作了比较。     

亚燃/超燃双重燃烧冷态掺混的数值模拟

通过求解完全的二维Ｎ－Ｓ方程组，对亚燃／超燃双燃烧室冲压发动机超音燃烧室的冷态掺混合性层流流动进行数值计算，采用ＭａｃＣｏｒｍａｃｋ两步显式格式分析一般音速来在股超音速气流间掺混的流动特征，并用代数法生成计算网络。给出的三组算例与文献（１）的实验结果进行对比，结果令人满意。