三维化学非平衡超燃整体发动机计算研究.pdfVIP

第十一届全国激波与激波管学术会议 三维化学非平衡超燃整体发动机计算 赵慧勇 乐嘉陵 中国空气动力研究与发展中心吸气式高超声速技术实验室 四川绵阳621000 1、引言 超燃冲压发动机主要由进气道、燃烧室和尾喷管组成。进气道研究的核心问题是进气道 中的边界层特性及激波与边界层的干扰问题。文献【1】指出，激波与边界层的干扰效应虽然已 经研究了近50年，但仍然有一些重要的数值不能很好的模拟；在燃烧室内，燃料通过燃烧把 其中的化学能以热能的形式释放出来，燃烧室的问题是整个高超声速冲压发动机研究的关键 所在。燃烧室内的流动充满激波、边界层、旋涡、非平衡化学反应及其相互作用的三维流动， 流动非常复杂。尾喷管提供了大部份推力，对超燃冲压发动机的整体性能影响很大。 Tanil21通过几何方法研究了进气道和隔离段，修改了外型，认为凸型进气道减弱了唇口激 波强度，但对拐角处的膨胀波影响不大。锯齿壁面的隔离段可以抵消拐角处的膨胀波。 式，湍流模型采用Wilcox的k一功和Menter的k—toSST模型，化学模型采用的7组分8方 式计算不能捕捉到空气泡的横向流动，需要采用半宽计算。侧壁的分离使空气偏向对称面。 他计算的壁面压强和上游的相互作用区明显小于实验值。此外，他的计算结果显示出非对称 的流场，然而这与实验相矛盾。Davidenk014]研究了管道中的氢气燃烧问题。管道包括喷管和 燃烧室。他采用MSD软件，MUSCL插值和Roe通量差分，湍流模型采用带壁面函数的k—E 方程，采用了四种H2／Air化学模型。他指出采用简单的化学模型不能正确地模拟点火延迟和 了半宽和全宽的NAL模型，湍流模型采用RNGk—E方程。他先后采用定常和非定常计算， 仔细分析了上游相互作用，包括旋涡的变化。他指出当上游相互作用区的长度超过特定的值 后，上下对称的假设不能成立。而且这种非对称是物理现象，不是数值模拟带来的。 以上这些文献都没有给出完整的整体发动机的计算和受力的分析。除了Mohiedinl51以外， 它们还缺少流线的细节。本文针对一个超燃整体发动机模型，采用H2作为燃料，对三种油气 比进行了计算，计算结果与实验基本吻合。并且分析了发动机的性能。 2、计算方法 本文针对三维非平衡化学反应N．S方程，时间推进采用LU．SGS隐式方法，无粘通量利 模型采用7组分8方程H2／Air模型。详细的公式可参见文献【6】。湍流模型采用Kokvl的 TNTk一倒双方程湍流模型，该模型比较圆满地解决了k一仞的初值依赖性问题。湍流模型与 流动主控制方程采用解藕计算。非平衡N．S方程计算量相当大，本文采用MPI并行方法，在 江南计算技术研究所的神州大规模并行机上实现了并行计算。 计算边界条件：壁面为绝热非催化壁面；入口边界将来流参数固定；出流边界采用一阶外 推；喷嘴出121的k和仞的法向梯度等于0，其余参数都固定。 3、计算条件 3．1实验模型 计算模型是中国空气动力研究与发展中心超高速研究所提供的。图1是试验的发动机模 型。进气道是具有三波系外压缩段的内外混合压缩式，气流总转角为23．5度，进气道迎风面 232 三维化学非平衡超燃整体发动机计算 扩型燃烧室由台阶、等截面段和扩张段组成，燃烧室扩张比为1．8，燃烧室总长320ram；模型 有一个短的尾喷管，扩张比4．9。模型总长860M。 在50mm长的等截面段末端，设4个燃料喷射支板。在每个支板迎气流方向的端面上， 沿高度均布四个0．5mm直径的燃料喷嘴。可参见图2。 图1超燃整体发动机示意图 图2支板 3．2来流条件 M∞ T∞ (K) P。(pa) P。(kg／m3) 6．OO 235．4543 2978．295 0．04 氧气质量分数 氮气质量分数 水质量分数 分子量(kg／mole) 0．255 0．574 O．17l 0．02634 3．3喷流条件