超声速进气道起动与不起动的流动特征结构的机理分析.pdf

第 卷 第 期 年 月气 体 物 理 .:./.超声速进气道起动与不起动的流动特征结构的机理分析陈雅倩 胡科琪 王高峰(浙江大学航空航天学院 浙江杭州)()摘 要:基于 等实验的超声速进气道模型 采用 计算方法分析超声速进气道起动和不起动流场特性 通过拓宽计算域和采用来流边界层自由发展方法 准确预测了典型的进气道不起动过程中可能出现的周期性振荡流现象 包括进气道内部高压的产生和降低、下壁面大尺度分离泡的膨胀收缩和迁移 并伴随不起动激波的传播 进气道完全起动状态(末端活门挡板角度)时得到的波系结构、壁面压强和流向速度分布计算结果与实验测量值相吻合 不起动状态()时流场的振荡周期和振幅与实验结果一致 对进气道不起动的非定常流场进行动态模态分解 发现了 个特征频率:主频 .的流场模态揭示了进气道出口的压强振荡最强 而入口及上壁面的速度振荡最强 二倍频.和三倍频.捕捉到的流场模态主要是离散的小尺度高能结构 在不起动状态的振荡过程中 进气道入口外部流场产生了较大的速度和压力脉动 所以对进气道内外流场相互作用的准确描述是预测不起动状态振荡流动的重要因素关键词:进气道 不起动 振荡流动 激波 动态模态分解 中图分类号:.文献标志码:收稿日期:修回日期:基金项目:国家重大科技专项()第一作者简介:陈雅倩()女 硕士 主要研究方向为计算流体力学:.引用格式:陈雅倩 胡科琪 王高峰.超声速进气道起动与不起动的流动特征结构的机理分析.气体物理 ():.:.():.:.:第 期陈雅倩 等:超声速进气道起动与不起动的流动特征结构的机理分析引 言超声速进气道作为冲压发动机的重要组成部分 能够在起动状态下工作是确保推进系统正常运行的基本要求关于进气道不起动问题 自 年 首次观察到超声速进气道入口振荡流动以来 国内外针对这种特殊的振荡现象进行了大量的研究 等对 数为 的超声速进气道进行了风洞实验 在隔离段末端安装挡板来模拟背压通过改变挡板角度 观察到进气道起动、不起动现象 发现不起动状态下进气道内存在 种不同振幅和频率的振荡流动 并且振荡流动中存在激波的运动 等对来流 数为 的矩形进气道的不起动振荡流进行了实验研究 提出以对流、激波串运动和声波传播共同建立反馈回路这一新的振荡机制 使用该机制估计的振幅和频率与实验结果基本一致 等通过实验研究了来流 数为.时进气道的不起动流动 观察到两种新的振荡模式 即混合振荡模式和间歇振荡模式为认识进气道振荡流机理、预测和控制进气道不起动提供了参考 中国科学技术大学的 等在不同出口节流比下对来流 数为.的高超声速进气道进行了实验研究 发现在低节流比下进气道起动 而在近堵塞节流比下进气道不起动 并且出现激波振荡现象 振荡频率随节流比的增加而增加得益于数值仿真技术的发展与应用 超声速进气道数值模拟方法可以对进气道不起动问题进行仿真预测 并提供详细的流场信息分析不起动机理 袁化成等对二维高超声速进气道的不起动和再起动过程进行了数值分析 发现进气道不起动后能否再起动与来流 数有关 常军涛等对高超声速进气道二维流场进行了数值模拟 从流动稳定性的角度解释了进气道不起动/再起动特性形成的原因 等利用大涡模拟()方法对 的实验进行了数值研究 准确预测了进气道起动状态下的流场特征 但对于进气道不起动状态 由于没有考虑风洞实验中实际的内外流耦合情况 等的 结果没有得到实验中观察到的振荡模式 面向流动精细化结构分析时具有较好的求解能力 但是计算代价较高 在使用其计算进气道内外流的耦合问题时对计算能力仍然存在较大的挑战 本文基于 实验中的超声速进气道实验模型 采用 平均()数值计算方法 拓宽计算区域表征内外流场耦合条件和来流湍流边界层发展情况 分析 进气道在起动和不起动状态下的定常和非定常流场 并通过动力学模态分解()方法对进气道不起动振荡流进行模态特征分析 物理模型及数值仿真方法.物理模型本文 数 值 模 拟 分 析 所 采 用 的 研 究 对 象 是 等在德克萨斯大学风洞中开展实验的物理模型 如图 所示 该模型由长度 .的内收缩进气道和长度 .的等截面隔离段组成 其中进气道入口收缩角 入口高度.隔离段高度.模型内部宽度 .在隔离段出口设置了活门挡板装置 实验中通过调整挡板角度()改变通流面积 以此模拟发动机燃烧产生的反压 从而观察进气道从起动状态到不起动状态的流动变化过程 表 给出了实验中提供的自由来流条件 来流 数为.总温 为 、总压为.实验通过纹影成像技术对激波结构进行观察分析 同时采用粒子图像测速()成像技术测量流场 在实验模型的下壁面设置了 个压强监测点 (见图()记录壁面压强随时间的变化情况实验中进气道入口前的来流边界层为湍流边界层并测得来流条件下入口下壁面 位置边界层厚度()为.静压为 ()()图 进气道示意图和实物图.气 体 物 理 年 第 卷表 进气道自由来流参数 /(/).数值仿真方法.计算模型和边界条件本文计算基于 软件 使用有限体积法求解三维可压缩完全气体 方程 空间离散采用迎风格式 并采用 湍流模型对进气道的起动和不起动状态分别进行 和 计算 收敛残差设置为.亚迭代步数为 计算的时间步长为.图 给出了本文根据实验物理模型得到的计算模型 超声速进气道模型的长度为.高度为.为了保证计算边界条件与风洞实验条件更好地保持一致 计算区域拓展到实验件前、后 纵向宽度扩展到 针对进气道不起动状态 除了 工况外 还设置了 的计算模型(见图()分别计算特定 角度下进气道内部的流动 计算入口边界条件设为超声速进口给定来流条件 计算域下边界及进气道表面设为无滑移绝热壁面 上边界及下游出口边界为超声速出口边界条件 计算域展向宽度设为 两侧设为周期性边界()()图 计算模型.计算网格及方法校验针对进气道模型()设计了 套不同的网格方案 以验证数值方法的网格无关性 进气道局部网格如图()所示 表 为具体的网格参数 其中网格 与网格 中第 层网格高度均为.而网格 与网格 中第 层网格高度为.图 给出了 套网格计算得到的壁面压强分布与实验测量值和他人 结果的比较 可以发现除了网格 计算得到的压强曲线有两处明显的凸起变化 网格 和网格 的压强分布基本保持一致 说明网格 的精度已经能够满足 计算需求 因此选择网格 开展进一步计算分析 相关计算结果与参考文献中的 结果也较吻合 与实验结果基本保持一致 因此 本文基于网格()的设置参数完成不同 角度模型的网格划分 模型网格单元总量约为.个表 网格参数 /.图 不同网格的壁面压强分布曲线.结果与讨论.完全起动状态时 进气道完全起动 计算得到的数值仿真结果如下 进气道入口前 壁面处的流向速度沿高度方向的变化曲线如图 所示 进气道入口来流条件与实验中的来流条件吻合较好:处的湍流边界层厚度为.(.)与实验测量值.基本一致 计算得到 点静压为 实验测量值为 相对误差为.第 期陈雅倩 等:超声速进气道起动与不起动的流动特征结构的机理分析图 对称面上 处流向速度分布.如图()所示为实验得到的起动状态下流场纹影图 可以看到其中存在复杂的波系结构图()为计算得到的进气道对称面纹影图和压强云图 比较实验和计算结果可以发现 数值模拟成功捕捉到了实验中观测到的 这 道反射激波、膨胀波 以及上壁面边界层 由于实验测量装置的限制 无法在实验中测量进口段的入射激波 而计算结果可以得到更完整的激波结构 从而分析进气道内的流动特征他人研究指出激波边界层相互作用是影响进气道内部流动的主要因素之一 图 是图()纹影图的红框区域 为数值计算得到的局部流场结构 与激波与边界层作用的理论流场结构基本一致 即入射激波与湍流边界层相互作用后引起湍流边界层气流分离 在分离区的分离点和再附点分别产生一道激波 并在下游相交()(.):()图 进气道对称面纹影图和压强分布.图 激波和湍流边界层相互作用流场结构.气 体 物 理 年 第 卷图 为隔离段/.范围的流向速度分布云图 图 为该范围内不同位置(/.)的流向速度分布曲线 计算结果与实验测量值整体上基本一致 与参考文献中的 结果也比较吻合 说明 对进气道内部流动的预测较好.对进气道不起动特性的影响通过改变 大小 对进气道进行非定常仿真计算 得到 共 个工况的计算结果 图 是计算得到的各工况 点压强随时间的变化曲线 对于 流场完全发展后 点压强已不再随时间变化 分别稳定在.时 计算稳定后 点压强在.范围内波动 振幅只有.考虑到计算误差 可以认为时流动是定常的 而 和 时 点压强随时间呈现周期性大幅振荡模式变化 周期分别为.和.振幅分别为.和.前者的振幅略小而周期较大 相应地 每个周期内 的 点压强比 多一个极大值()()()()()图 流向速度.图 不同位置流向速度分布曲线.图 为不同 对应工况下的进气道对称面的平均压强分布云图 与 的流场相比 活门挡板抬起后 出口气流受阻导致气流压强升高 将造成出口堵塞引起进气道不起动 其中 的进气道内部流场完全发展后保持稳定状态由于三者 相差不大 对应的流场结构基本一致并且由于中压区稳定在下游 并未向上游传播 此时流场中仍存在入射激波、膨胀波、反射激波 和 进气道此时还可以起动 而 时 出口附近压强较高 相比于 高压区向上游发展 阻止了反射激波 的形成 注意到 的流场中 高压区前存在一道激波 将其称为不起动激波 和 的平均流场中 高压气流几乎占满了整个通道 容易导致进气道出口气流反向流动 并且高压区阻断了入射激波的发展 没有形成反射激波 进气道为不起动状态为了定量比较 对进气道内部流动的影响图 给出了不同 对应的进气道下壁面平均压强分布曲线 从图中可知 进气道壁面压强沿程逐渐增加 其中 对应的压强突增点为激波第 期陈雅倩 等:超声速进气道起动与不起动的流动特征结构的机理分析所在位置 比较不同 对应的压强曲线 发现随着 的增大 壁面压强逐渐增大 并且 时 进气道出口附近的压强相比于 显著增加 结合前文分析 可以认为 为进气道进入大幅周期性振荡不起动模式的临界角度图 点压强随时间变化曲线.图 进气道在不同 下的平均压强分布.图 不同 下壁面平均压强分布曲线.图 是 时数值计算和实验的 点压强变化曲线 其中黑色实线为参考文献中的实验结果 红色实线为计算结果 对比实验与计算结果 虽然数值计算和实验结果由于初始设置不同在压力脉动的峰谷值上存在一定差异 但 点压强随时间变化的周期和振幅基本一致 周期为.左右 振幅均在 左右 因此数值计算具有一定的可信度图 时 点压强的实验与计算结果对比.不起动流场及其 特性分析.流场发展过程时进气道内部流动表现为高度非定常性 图 也表明此时流场结构呈周期性变化 是典型的进气道不起动流动模式图 分别为 时一个周期内进气道流场发展过程的压力场和流线图 图中 时刻与图 对应 在 时刻 进气道入口被大尺度分离泡 堵住 高速自由来流被阻挡在进气道外 但由于进气道内部压强较低 自由来流逐渐推动分离泡 向下游移动并逐渐缩小 使得高速气流能够从上部分入口进入进气道()注意到分离泡 前的弱激波也随之向进气道内部移动 当气流持续进入后 由于出口气流受挡板阻挡并且出口面积减小 在 时刻 出口附近开始出现高压区 逆压梯度导致下壁面形成回流 出现分离区 由于高压逐渐扩大 在上游形成不起动激波()分离区向上游延伸至不起动激波后 而主流区主要分布在 和 之间 且沿上壁面向下游流动 这也抑制了上壁面分离泡的发展气 体 物 理 年 第 卷图 的进气道对称面压强分布随时间的变化.图 的进气道对称面流线随时间的变化.分离泡 随自由来流继续向进气道内移动并由于进气道入口主流面积的扩大 分离泡 逐渐缩小 同时也因为大量气流到达下游 高压区进一步扩大 下游的分离区 也推动不起动激波向前传播 直至 融入 前的弱激波也将耗散()高速主流主要沿上壁面分布 此时入口起动 有更多气流进入进气道内 在 时刻 进气道中下游充满高压气流 不起动激波后的分离区膨胀 为其提供前传的能量 此时不起动激波上缘恰好位于进气道入口 随着分离泡继续膨胀、不起动激波向进气道外传播 将入口堵住阻挡了高速自由来流 入口不起动 进气道内的高压气流将从出口排出 压强降低()当进气道下游压强降至来流水平 不起动激波前传至约/后不再前传此时()流场结构与 时一致 不起动激波将随来流向下游移动 进入下一周期上述结果表明 与 相比 时 进气道不起动流场中存在不起动激波的前传播、分离气泡的膨胀收缩和迁移以及下游高压区的产生与降低等流动特征 此时的进气道不起动是高度非定常流动的过程 并且由于计算域向前及纵向拓展 使得不起动激波和分离泡发展至进气道入口前 与自由来流相互作用 获得返回进气道内的能量 从而形成周期性振荡流动由于 时进气道不起动的流动情况与 类似 这里不再赘述 表 分别列出了、时进气道不起动流中不起动激波从下壁面 点传播到 点的平均速度以及流场的振荡周期和主频率 其中频率为周期的倒数 对比可得 时不起动激波的平均传播速度较小 相应的流动周期较长表 不起动激波的平均传播速度 /()/(/)/.特性分析由 提出的动力学模态分析方法提供了一种从均匀采样非定常流场中提取动力学信息的工具 能够用于分析流场复杂的流动特征 与快速 变换()相比 的特点在于从复杂流场中提取单一频率的流动模态 并且可以通过各模态的特征描述流动随时间的演化过程 具有时空耦合的独特优势 能够帮助研究人员关注特定的流动特征和潜在的物理机制 同样可以用于进气道不起动振荡流的后处理分析对 进气道对称面流场进行 处理选取 个样本数据 采样频率为 将 个时间点的流场数据导入 程序进行第 期陈雅倩 等:超声速进气道起动与不起动的流动特征结构的机理分析分析 提取其动态变换过程的特征值和特征模态其中得到的 频谱分布如图 所示 幅值 随频率 的增加呈现下降趋势 即低频模态的幅值较大 其中幅值最大的模态频率为 为静止模态 运动模态最大的频率是.这个频率与.节 流场的振荡周期(.)对应的主频相近 另外 频谱也包含了各个倍频 如二倍频.、三倍频.图 频谱分布.图 分别为 处理后 主频 对应压力场和速度场的幅值和相位图 从压力场幅值分布中可以看到进气道出口附近的流场脉动最强 整体上脉动强度向上游递减 入口上壁面也存在中等脉动强度的流团 相位图表明隔离段内流场动态变化的相位差不大 但与入口上壁面流团的相位差约为/与进气道入口及入口外的流动相位差约为 说明进气道内部和入口压力的动态变化是反向的 从图 速度场模态的幅值分布中可以发现峰值主要集中在上壁面中段(/)入口段(/)上壁面附近以及入口前至下壁面/处的幅值也较大 而从/.到下游 沿下壁面附近的幅值偏小 并且进气道出口附近整个通道内的幅值都较小 说明进气道内部中上游的上壁面以及入口前气流速度的变化程度较大 结合.节 发现这是由于入口大尺度分离泡和上壁面小尺度分离泡的运动造成的流场动态特征 另外 虽然图()相位图中进气道上下壁面附近的颜色差异较大 但由于相位差约为 说明两个区域的速度变化是连续的 而入口下壁面分离涡所在区域与其他流域的相位差约为 即分离涡的运动与进气道内其他流动变化是相反的()()图 频率 对应 模态的幅值和相位(压力场).()()()图 频率 对应 模态的幅值和相位(速度场).相对于 的高频运动模态也有助于理解流场动态变化的规律 图 为 频率对应 模态的压力场分布 其中 的流场特征结构沿上壁面分布 的流场结构除在上壁面外 还包括流道中部和挡板附近 与图 相比 这两个频率下的流场模态主要是离散的小尺度高能结构().气 体 物 理 年 第 卷().图 对应 模态的幅值(压力场).()结论本文基于 等实验的超声速进气道实验模型 采用 计算方法分析超声速进气道起动和不起动流场特性 通过计算域向前和纵向拓展 不仅得到与风洞实验基本一致的来流湍流边界层厚度和进气道入口条件 还能够满足进气道内外流场相互作用的条件 从而在计算分析中得到不起动振荡流动特征针对不同 模型的计算 发现随着 的增大进气道出口和壁面压强逐渐增加 进气道由起动状态进入不起动状态 其中 流场为无振荡稳定模式 流场为大幅低频振荡模式 进气道不起动 在 的进气道完全起动流中 计算得到的壁面压强分布、不同站位流向速度分布均与实验测量值相吻合 不起动工况中的振荡周期和振幅与实验测量值一致对 流场动态发展过程的分析发现进气道不起动是由出口面积减小造成压强增加导致的相比于前人 计算 拓展计算域后实现了进气道内外流的耦合 得到典型的进气道不起动振荡流主要流动特征为出口高压区向上游的发展及压强再次降低、下壁面大尺度分离泡的膨胀收缩和前后移动 并伴随不起动激波的前后传播 对 的流场进行 分析 得到流动主频 .、二倍频.、三倍频.其中频率 的流场模态表明进气道出口附近的压力场振荡最强、速度变化较弱 而出口附近和上壁面的速度场振荡较强 与 相比 高频、捕捉到的流场模态主要是离散的小尺度高能结构参考文献()./.:././.():.():.():.():.袁化成 梁德旺.高超声速进气道再起动特性分析.推进技术 ():.():().常军涛 于达仁 鲍文.攻角引起的高超声速进气道不起动/再起动特性分析.航空动力学报():./.():().():.阎超.计算流体力学方法及应用.北京:北京航空航天大学出版社.:().第 期陈雅倩 等:超声速进气道起动与不起动的流动特征结构的机理分析 .():./:?.():.张悦 谭慧俊 王子运 等.进气道内激波/边界层干扰及控制研究进展.推进技术 ():./.():().:.:.:.寇家庆 张伟伟.动力学模态分解及其在流体力学中的应用.空气动力学学报 ():.():().王建明 王涵 桂琳.压气机叶栅叶顶间隙流的动力学模态分解.推进技术 ():.():()./.:/./.