一种模拟边界层泄流与亚声速外流耦合的实验台设计方法

1.本发明涉及亚声速流动实验领域，尤其是一种能够模拟边界层泄流与亚声速外流耦合作用的实验台设计方法。


背景技术：

2.进气道作为吸气式动力飞行器的必备部件，承担着捕获自由来流、向下游发动机提供适量的空气流量和对来流压缩的责任，其性能的好坏，直接影响着整个推进系统的工作效率。根据相关研究表明，进气道总压恢复系数降低1％，动力系统的推力损失1.2％～1.5％。同时，进气道出口流场畸变也直接影响到发动机的工作包线。此外，进气道还是飞机前向三大雷达强散射源之一，其rcs约占整机前向rcs的30～50％。
3.进气道具有内流机理复杂、上下游工况多变、气动性能要求高等特点，所面临的气动设计挑战大，为此流动控制措施在进气道气动设计中得到了广泛使用。较为常见的进气道流动控制措施有边界层吸除、边界层隔道、进口鼓包、涡流发生器等。其中，对于边界层吸除流动控制措施而言，影响其控制效果的关键因素有边界层吸除装置进气效率、内通道流动阻力和泄流阻力。进气效率与内通道流动阻力一般分别由边界层吸除装置入口几何构型与内通道集合构型决定；泄流阻力主要受外流影响，同时边界层泄流与外流具有强耦合作用，因此，泄流阻力对边界层吸除装置流动控制效果影响更为显著。
4.目前针对内外流耦合干扰的研究主要集中于高压气源射流与外流的耦合作用研究，而对于边界层泄流与外流耦合作用(无外接高压气源)的研究较为缺乏。此外，在现有常规风洞中，难以对边界层泄流与外流耦合作用进行详尽的研究。


技术实现要素：

5.为解决上述内容，本发明的目的在于提供一种便于研究边界层泄流与亚声速外流耦合作用的实验台设计方法，该实验台可用于研究不同来流马赫数以及不同来流边界层厚度条件下边界层泄流与亚声速外流的耦合作用。此外，该实验台结构简单、尺寸小、便于维护，对实验场地和气源能力要求较低。
6.为了达到上述目的，本发明采用的技术方案如下：
7.一种模拟边界层泄流与亚声速外流耦合作用的实验台设计方法，其特征在于，包括以下步骤：
8.(1).设置两路管道，根据实验模型确定两路管道基本构型；每路管道均包括依次延伸的入口段、整流段、收缩段、边界层发展段、试验段、转接段；
9.(2).为使试验段入口上游流场发展一致，两路管道试验段入口截面上游管道几何尺寸一致；
10.(3).根据实验室已有真空气源抽吸能力、实验时长以及试验段马赫数确定管道喉道面积；
11.(4).将喉道形状与尺寸确定后，在保证两路管道互不相干的约束条件下，确定实
验台入口段、整流段和收缩段几何形状与尺寸；入口段曲线方程为：
12.ρ2＝a2cos(2θ)
13.其中，ρ为曲线在极坐标系中的半径，θ为曲线在极坐标系中的角度，a为曲线特征参数，点a与点b为当θ取值范围为(7π)/8～(5π)/4时曲线的端点；入口段bc段曲线方程为：
14.y＝-1
×
10-5
x
5-4
×
10-7
x
4-0.0005x
3-0.3628x
2-121.79x-15871x∈[-872.4,-577.0]
[0015]
其中，x为流向坐标，y为法向坐标，下同；当x＝-872.4、x＝-577.0时，由bc段曲线方程确定点b、c法向坐标y值；整流段为平直段，收缩段de段曲线方程为：
[0016]
y＝-1.447
×
10-9
x
7-4.32
×
10-6
x
6-5.522
×
10-3
x
5-3.919x
4-1668x
3-4.257
×
105x
2-6.032
×
107x-[0017]
3.662
×
109x∈[-459,-400]；
[0018]
当x＝-459、x＝-400时，由de段曲线方程确定点d、e法向坐标y值，由此得到入口段、整流段和收缩段对称面上壁面型线；
[0019]
(5).步骤(4)将各段曲线函数确定后，将收缩段de曲线以喉道截面轮廓线为引导线进行扫掠形成管道曲面，由此得到收缩段对称面下壁面型线lm；
[0020]
(6).以步骤(5)得到的收缩段入口轮廓为基准，沿曲线dcb进行扫掠，随后以扫掠得到的管道终截面轮廓为基准，沿此终截面进行扫掠，最终得到入口段和整流段型面，进而得到入口段和整流段对称面下壁面型线ijkl；
[0021]
(7).为消除壁面黏性效应的影响，调整边界层发展段上壁面沿流向扩张角以保证主流区域内流向马赫数分布一致，根据实验具体要求设置边界层发展段沿流向长度，由此得到边界层发展段对称面型线，边界层发展段对称面型线为上壁面ef及下壁面mn，随后以喉道截面轮廓为基准，沿边界层发展段对称面型线进行扫掠，得到边界层发展段管道型面；
[0022]
(8).两路管道中，部分气流经过实验模型从一路管道流入另一路管道，因此实验模型入口下游一路管道气流流量有所减小，另一路管道试验段出口流量有所增加；根据上述流量变化规律，对实验模型入口下游一路管道试验段进行型面迭代优化设计；
[0023]
(9).另一路管道试验段有流量注入，且管道内流动为亚声速流动，因此为避免对实验模型出口上游流场产生影响，该另一路管道试验段上壁面在某点x'处增大上壁面扩张角度，在另一路管道试验段与实验模型出口相交的下游某点y'处壁面扩张角转为0.3
°
，随后在点x'和点y'处对上壁面进行倒圆处理；上壁面倒圆处理后，在实验模型出口上游1.33h处扩张角由0.3
°
逐渐增大至5
°
，在点p1'处扩张角逐渐由5
°
减小至0.3
°
，点w'下游壁面扩张角保持为0.3
°
，p1'为点p'沿-z方向在上壁面的投影；
[0024]
(10).将边界层发展段出口截面轮廓线以步骤(8)和(9)得到的试验段对称面型线为引导线进行扫掠，得到试验段型面；试验段对称面型线为上壁面fg/f'g'及下壁面nq/n'q'；转接段型面根据试验段出口截面gq(g'q')光滑过渡至流量调节阀入口截面st(s't')得到，转接段扩张角不大于8
°
；最终得到的实验台两路管道型面。
[0025]
进一步的，实验台入口来流条件为当地大气环境，由试验段下游真空气源提供负压条件，驱动气流在试验段达到实验所要求的速度。
[0026]
进一步的，实验台实验时长、试验段尺寸和试验段气流马赫数受真空气源能力决
定。
[0027]
进一步的，实验台整流段长度不小于2倍收缩段出口高度，收缩段收缩比不小于5。
[0028]
进一步的，试验段型面在保证与试验段上下游型面光滑过渡的条件下根据具体实验模型进行更换。
[0029]
进一步的，另一路管道试验段上壁面在点x'处增大上壁面扩张角度至5
°
。
[0030]
有益效果：
[0031]
与常规亚声速风洞相比，本发明实现的实验台可精细化地模拟边界层泄流与外流，进而可对边界层泄流与外流的耦合作用进行详尽地试验研究，为可实现高效费比流动控制的边界层吸除流动控制装置(依靠边界层内气流自身动能驱动，无外接高压气源)的设计提供数据支撑。此外，该发明实现的实验台模拟来流马赫数范围较广，且对实验场地及气源要求较低，为研究边界层泄流与亚声速外流的耦合作用提供了一种切实可行的实验台设计方法。
附图说明
[0032]
图1是实验台整体部件结构示意图；
[0033]
图2是本文实验台设计案例中所使用的实验模型(边界层吸除装置)；
[0034]
图3是实验模型与实验台主/副管道连接结构示意图；
[0035]
图4是双通道抽吸实验台主管道对称面型线；
[0036]
图5是双通道抽吸实验台副管道对称面型线；
[0037]
图6是实验台喉道下游管道横截面构型；
[0038]
图7是收缩段管道半截面构型；
[0039]
图8是收缩段入口(整流段)截面轮廓尺寸；
[0040]
图9是入口段、整流段和收缩段半截面构型；
[0041]
图10是主管道试验段对称面型线；
[0042]
图11是副管道试验段对称面型线；
[0043]
图12是实验台主/副管道整体型面构型；
[0044]
图13是试验段入口气流马赫数为0.54条件下主/副管道对称面距离下壁面0.26h、0.50h和0.74h度处流向马赫数分布；
[0045]
图14是试验段入口气流马赫数为0.75条件下主/副管道对称面距离下壁面0.26h、0.50h和0.74h处流向马赫数分布。
具体实施方式
[0046]
本发明公开一种模拟边界层泄流与亚声速外流耦合作用的实验台设计方法。请参阅图4、图5、图6、图7、图8、图9、图10、图11、图12所示，下面对采用本发明方法设计的试验段马赫数为0.50～0.75、两通道流量之和最大为1.5kg/s的实验台设计方法的详细实施步骤进行叙述(需要说明的是下述实验台设计方法已考虑了试验段实验模型，实验模型为埋入式微小通道，其构型如图2所示，实验模型入口连接主管道、出口连接副管道，连接结构如图3所示)。
[0047]
(1).主管道和副管道对称面型线如图4和图5所示。在实验台初始设计阶段，首先
根据实验模型(图2)确定管道基本构型，由于本案例中实验模型入口与出口均为平面，因此主/副管道实验模型安装面(边界层发展段和试验段下壁面mnq段/m'n'q'段)均为平面(实验模型入口安装位置为图4所示op段，实验模型出口安装位置为图5所示o'p'段)，且沿x方向(流向)保持水平；边界层发展段管道上壁面(ef段/e'f'段)沿流向存在0.3
°
扩张，以修正边界层黏性效应的影响；管道横截面构型如图6所示，其中，管道展向尺寸w沿流向保持不变。
[0048]
(2).为保证试验段入口上游流场发展一致，主管道与副管道试验段入口截面(截面fn/截面f'n')上游管道几何尺寸一致，此外，由于实验模型入口平面法向与实验模型出口平面法向垂直，为满足实验模型的安装条件，主/副管道边界层发展段和试验段下壁面(mnq段/m'n'q'段)法向也垂直。
[0049]
(3).根据实验室已有真空气源抽吸能力、实验时长以及试验段马赫数确定实验台管道喉道面积(喉道截面位置为截面em和e'm')，喉道截面构型如图6所示，其中，w＝94mm，h＝38mm，r1＝20mm，r2＝5mm，喉道截面面积为3389.6mm2。
[0050]
(4).步骤(1)～(3)将喉道形状与尺寸确定后，在保证两路管道互不相干的约束条件下，确定实验台入口段、整流段和收缩段几何形状与尺寸。本案例中入口段曲线方程为：
[0051]
ρ2＝a2cos(2θ)
[0052]
其中，ρ为曲线在极坐标系中的半径，θ为曲线在极坐标系中的角度，a为曲线特征参数，点a与点b为当θ取值范围为(7π)/8～(5π)/4时曲线的端点；本案例中a＝63mm，θ取值范围为(7π)/8～(5π)/4。入口段bc段曲线方程为：
[0053]
y＝-1
×
10-5
x
5-4
×
10-7
x
4-0.0005x
3-0.3628x
2-121.79x-15871x∈[-872.4,-577.0]
[0054]
其中，x为流向坐标，y为法向坐标，下同；当x＝-872.4、x＝-577.0时，由bc段曲线方程确定点b、c法向坐标y值。整流段为平直段，长度为3.1h(h为喉道高度)，收缩段de段曲线方程为：
[0055]
y＝-1.447
×
10-9
x
7-4.32
×
10-6
x
6-5.522
×
10-3
x
5-3.919x
4-1668x
3-4.257
×
105x
2-6.032
×
107x-3.662
×
109x∈[-459,-400]
[0056]
当x＝-459、x＝-400时，由de段曲线方程确定点d、e法向坐标y值。由此得到入口段、整流段和收缩段对称面上壁面型线。
[0057]
(5).步骤(4)将各段曲线函数确定后，将收缩段de曲线以喉道截面轮廓线为引导线进行扫掠形成管道曲面，图7给出收缩段管道半截面构型，图8所示为收缩段入口(整流段出口)截面轮廓尺寸，收缩段收缩比为5.35。
[0058]
(6).以步骤(5)得到的收缩段入口轮廓为基准，沿曲线dcb进行扫掠，随后以扫掠得到的管道终截面轮廓为基准，沿此终截面进行扫掠，最终得到入口段和整流段型面，入口段、整流段和收缩段半截面构型如图9所示。
[0059]
(7).为消除壁面黏性效应的影响，边界层发展段上壁面沿流向扩张角为0.3
°
，以保证主流区域内流向马赫数分布一致，边界层发展段沿流向长度为9.71h，由此可得到边界层发展段对称面型线(上壁面ef，下壁面mn)，随后以喉道截面轮廓为基准，沿边界层发展段
对称面型线(ef和mn)进行扫掠，得到边界层发展段管道型面，边界层发展段出口高度为1.05h。
[0060]
(8).试验段型面构型需根据具体实验模型而定。本案例中实验模型构型如图2所示，主管道部分气流会经过实验模型流入副管道，因此实验模型入口下游主管道气流流量有所减小，副管道试验段出口流量有所增加。根据上述流量变化规律，对主管道试验段进行型面迭代优化设计，对称面型线如图10所示。其中，fn截面为试验段入口，gq截面为试验段出口，oo1截面为实验模型入口前缘所在截面，pp1截面为实验模型入口后缘所在截面，试验段总长为10h。为修正试验段入口下游气流流量减小而带来的影响，o1p1段存在0.5
°
收缩，p1g段存在0.3
°
扩张，同时对o1p1段和p1g段在交点p1处进行倒圆处理倒圆半径r3＝200mm。
[0061]
(9).副管道试验段对称面型线如图11所示。其中，f'n'截面为试验段入口，g'q'截面为试验段出口，o'o1'截面为实验模型出口前缘所在截面，p'p1'截面为实验模型出口后缘所在截面，试验段总长为10h。考虑到副管道试验段有流量注入，且管道内流动为亚声速流动，因此为避免对实验模型出口上游流场产生影响，副管道试验段上壁面在点x'处增大上壁面扩张角度至5
°
，在实验模型出口下游点y'处壁面扩张角转为0.3
°
，随后在点x'和点y'处对上壁面进行倒圆处理(和)，倒圆半径r4和r5分别为900mm和250mm。上壁面倒圆处理后，在实验模型出口上游1.33h处扩张角由0.3
°
逐渐增大至5
°
，在点p1'处扩张角逐渐由5
°
减小至0.3
°
，点w'下游壁面扩张角保持为0.3
°
。
[0062]
(10).将边界层发展段出口截面轮廓线以步骤(8)和(9)得到的试验段对称面型线为引导线进行扫掠，即可得到试验段型面。实验台转接段型面可根据试验段出口截面gq(g'q')光滑过渡至流量调节阀入口截面st(s't')得到，但转接段扩张角应不大于8
°
。最终得到的实验台主/副管道型面如图12所示。
[0063]
针对本发明，采用数值仿真方法对上述技术方案的效果进行验证，模拟了试验段入口马赫数为0.54和0.75条件下实验台的流场品质，图13和图14给出了主/副管道对称面距离下壁面0.26h、0.50h和0.74h高度处流向马赫数分布(注：实验中主要关注实验模型上游流场)。由图可知，试验段入口上游主流区气流马赫数基本保持不变，试验段入口流场品质满足实验要求。
[0064]
另外，本发明的具体实现方法和途径很多，以上所述仅是本发明的优选实施方式。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。