Experimental study on structural characteristics of separation flow induced by 3D wedge in hypersonic laminar flow by oil visualization
-
摘要: 在高超声速风洞中采用油流显示技术开展了三维楔体诱导层流分离现象的实验研究。研究模型为矩形平板/三维楔体、三角形平板/三维楔体。研究结果表明:高超声速(马赫数6.0)局部层流分离结构不仅受楔体压缩角影响,受前体外形的影响也很明显;对于三角形前体,由于楔体上游来流存在横向流动,使得楔体诱导的分离流动结构完全不同于矩形平板模型,且在不同迎角下,分离线形状也存在较大差异;受三角形前体与三维楔体综合作用,诱导层流分离呈现非常显著的三维效应。Abstract: Experimental investigation of hypersonic laminar separation flow over a 3D wedge is carried out in the hypersonic wind tunnel by utilizing oil flow visualization technique. Two test models of rectangular/ triangular flat plate are used. The 3D wedges which have different compression angles are mounted on the afterbody of models. The experimental result shows that the forebody shape has a great influence on the structures of hypersonic laminar separation induced by the 3D wedge. For the triangular plate model, the flow in the upstream of the wedge has obvious crossflow effect, and the separation structure is completely different from that of the rectangular plate model. At different attack angles, the shapes of the separation lines are also very different. The result shows very significant 3D effect of the local laminar separation which is induced by the triangular forebody and the 3D wedge.
-
Keywords:
- hypersonic flow /
- laminar separation /
- 3D wedge /
- oil flow visualization
-
0 引言
0.6m×0.6m连续式跨声速风洞(以下简称0.6m风洞)为国内首座具备负压实验能力的连续式跨声速风洞,试验马赫数范围为0.2~1.6,压力范围为20~250kPa(绝对压力),具有3个试验段,该风洞通过阀门控制系统实现常压、增压和降压运行,其中增压、常压运行通过调节进排气系统实现,降压运行通过调节抽真空系统实现。
风洞的总压调节能力及精度是风洞的一个重要技术指标,对提高试验效率和试验数据精准度具有重 要作用。目前,国内连续式跨声速风洞数量较少,NF-6风洞作为国内首座连续式跨声速风洞,其马赫数范围较窄,试验工况相对较少[1- 2]。目前国内外风洞对马赫数控制方法描述较多,而对连续式风洞总压控制及马赫数和总压的耦合关系描述较少[3-11]。暂冲式风洞压力控制通过主调压阀和引射器调压阀等实现[9-11],而0.6m连续式风洞压力主要通过进排气和抽真空阀门/泵组合来实现,且总压调节受马赫数调节手段影响,可借鉴的内容较少。各种改进的模糊PID控制方法在风洞中尤其是马赫数或风速控制的应用较多[12-13],文献[13]对某暂冲式风洞主气流压力控制原理进行了简单分析,并引入模糊控制概念,总压控制精度达到0.5%,因此本文拟采用模糊PID来进行总压控制,但其控制参数大幅跳动时容易导致控制器输出跳动,进而导致压力波动,因而需要结合0.6m连续式风洞压力特性对控制算法进行优化。
0.6m风洞阀门系统结构见图 1。
不同于现有其他风洞,该风洞压力调节系统具有以下特性:
(1) 总压调节手段多,包含9台阀门和3台泵,不同的压力工况采用不同的阀门组合进行控制。
(2) 压力控制范围广,从绝对压力20kPa至250kPa,涉及常压、增压、负压3种不同的工况,且试验过程中可连续变总压、连续变马赫数和连续变迎角[14]。
(3) 压力控制精度要求高,设计指标为0.2%(常压和增压)[15]。
(4) 作为大型连续式风洞的引导性风洞,风洞运行过程中迎角、二喉道、压缩机转速、驻室抽气流量等随时需要调节,这些因素都对总压有较大影响,从而导致压力场波动较大,需要总压控制策略具备很强的适应能力[14-15]。
本文主要研究内容为:分析0.6m连续式风洞的压力控制特性,摸索多试验工况、宽压力调节范围下具有较好适应性的压力控制策略,并结合控制对象特性通过对模糊PID算法分段变参数的方式进行改进,并实施该风洞的压力控制,为风洞控制策略优化及试验效率的提高奠定基础,同时为后续的大型风洞建设积累经验。
1 风洞总压和马赫数耦合特性
连续式跨声速风洞主要通过调节进气/排气/抽真空阀门的开度来实现总压的调节,由于总压和马赫数在调节时相互耦合,致使调节难度加大,且调节时间较长,经过分析,总压和马赫数调节手段的耦合关系如图 2所示。
其中,G1(s)至G9(s)分别为每种控制量的传递函数,Gab(s)为总压对马赫数的影响函数,Gba(s)为马赫数对总压的影响函数,Ga(s)为总压控制传递函数,Gb(s)为马赫数控制传递函数。
为提高总压控制精度并减少调节稳定时间,理论上需要对总压和马赫数进行解耦控制,但目前总压和马赫数的解耦算法在工程应用中较少,更多是通过控制策略的优化并尽可能减少耦合的方式,因此,本系统设计时将多种马赫数调节手段调节时对总压的影响作为总压控制的扰动对待,同时风洞马赫数控制策略设计时减少总压和马赫数的耦合。0.6m风洞总压控制主要通过调节相应的调节阀开度实现,因此需首先建立调节阀开度和风洞总压的关系。对于可压缩流体,理想情况下,流量Q和压力之间关系为[16]:
(1) 其中,β为缩流处截面等效直径与管道截面等效直径之比; E为与流速有关的系数;A为缩流处截面积;ρ为密度;Q为体积流量; p1为入口静压;p2为出口静压,近似等于大气压;ε为流体膨胀系数。
同时,对于抛物线特性[17]的调节阀,在理想情况下其流量和阀门开度之间满足下式:
(2) 其中,Qmax为流经阀门的最大流量,R为阀门的可调比,为流经阀门的最大流量和最小流量的比值,国内调节阀R一般取为30或50,l/L为阀门开度的百分比。
由式(1)和(2)可得调节阀开度和静压的关系,而风洞总压由静压和动压组成,进而得到调节阀开度和风洞总压之间的关系:
(3) (4) (5) (6) 其中,p0为总压,f1(n)为风速和压缩机转速的函数关系,考虑到稳定段气流速度较低,可忽略不计,f2(α,r,l,m)为马赫数调节对总压的影响量,a为迎角变化梯度,r为压缩机静叶角,l为二喉道位移,m为驻室抽气量。
由此可见,连续式风洞总压影响因素包括:马赫数调节方式、气流介质密度(由风洞内总压和总温确定)、调节阀开度、气源压力等,本文重点对气源压力、阀门开度、压缩机转速和目标总压进行考虑以设计总压控制策略。
2 总压控制策略
为保证总压控制精度并确保风洞运行安全,考虑到PID算法具有结构简单、稳定性好、工作可靠、调整方便的优点,但同时存在自适应能力差的缺点,因此拟采用模糊PID算法,其对压力波动有较强适应性,但由于该连续式风洞压力范围较广,且不同目标压力区间下压力特性不同,若单纯采用模糊PID进行控制,可能导致不同压力工况参数匹配度较差,从而带来震荡或超调,因此该连续式风洞总压控制的核心算法采用分段变参数加模糊PID控制。
结合该连续式跨声速风洞压力变化规律,以及连续式风洞运行时压缩机充气密封系统将向风洞洞体内漏气的特性,为节省气源耗气量,确定该风洞总压控制策略为:
(1) 当工况为常压(控制压力目标在大气压至105kPa)时,由于压缩机充气密封往洞体漏气,因此仅利用排气调节阀进行总压控制;
(2) 当工况为增压(控制压力目标大于105kPa)时,通过预置进气阀开度和调节排气阀进行总压控制;
(3) 当工况为负压(控制压力目标低于94kPa)时,开启真空泵,通过预置进气调节阀开度、开启真空泵和调节抽真空调节阀进行总压控制;
(4) 当压缩机转速大幅上升或下降时,压力会大幅波动,通过计算目标转速和当前转速误差以及转速误差变化率,对参与调节的阀门采用分段限幅跟随调节的方式,以减少积分饱和和超调,同时为减少总压和马赫数耦合,变马赫数时采用多阶梯升压缩机转速的控制方式。
(5) 压力调节过程中,进气阀预置开度与气源压力和目标压力满足调试确定的函数关系,排气主调阀进行压力粗调(0.2%以上),待压力稳定后由排气辅调阀精调(0.2%以下至0.04%)。
3 分段变参数加模糊变参数PID算法
3.1 模糊PID控制算法
模糊PID控制算法以常规PID控制为基础,系统首先计算被控量的误差e和误差变化率ec作为二维模糊控制器的输入变量,然后通过模糊化、模糊控制规则推理、解模糊得到KP,KI和KD的修正量,在控制过程中在线实时整定PID各参数,其原理已经非常成熟,这里不再累述。
3.2 分段变参数加模糊变参数PI算法
由于风洞总压控制范围较宽,在不同的压力区间,风洞的总压特性不同,尤其当总压较高时(p0≥200kPa),风洞排气调节阀开度增大后总压下降明显,同时目前调节阀存在响应时间及执行速度较慢的问题,当压力快速变化时阀门的调节速度跟不上“节奏”,因此若不根据压力控制范围对控制参数进行调整,而仍使用常压时的基准调节参数,必然导致PID算法适应性差,产生较大的超调。因此,以模糊自适应PID为基础,增加目标总压作为输入变量,首先根据目标总压进行分段,基于不同的总压控制目标结合简单PI调试规律初步确定PI参数基准值,尔后由模糊算法进行变参数控制,以提高不同试验状态的自适应能力,原理框图见图 3。
3.3 模糊控制参数及规则表
阀门控制系统中,将计算出的目标压力、实际压力和给定压力的误差e和误差变化率ec作为模糊控制器的输入变量,PI控制器的2个参数KP、KI的修正量ΔKP、ΔKI作为输出变量。误差e基本论域取[-0.75,0.75],其偏差变化率ec的基本论域取[-0.5,0.5],输出变量ΔKP、ΔKI基本论域分别取为[-0.5,0.5]、[-1.5,1.5],E、EC、ΔKP、ΔKI分别是输入变量e、ec、输出变量ΔKP、ΔKI所对应的模糊语言变量。E、EC、ΔKP、ΔKI都分别取7个模糊子集,其中NB、NM、NS、ZO、PS、PM、PB分别代表负大、负中、负小、零(又分成正零ZP、负零ZN)、正小、正中、正大。模糊语言变量选择三角隶属函数。
调试过程先摸索出不同压力区间段的大体P、I参数,同时为防止控制输出突跳,将控制参数进行拟合处理后再作为模块控制的输入,见式7。
(7) 式中KP、KI为分段后基准比例系数和积分系数,p为目标压力。根据不同目标压力分段,能有效缩小适合当前目标压力的基准P、I值范围,基准P、I波动范围缩小,也就缩小了最终的P、I参数波动的绝对值,进而避免了控制参数大幅跳动的问题。
表 1给出了结合风洞压力特性的分段变参数模糊PID控制规则表。 PID的实时控制参数是在基准P、I参数的基础加上模糊算法得到的P、I参数增量。
表 1 阀门控制系统 参数Table 1 parameter of valve control systemE EC NB NM NS 0 PS PM PB NB PB/NL PB/NL PB/NM PB/NM PB/NS PS/0 PS/0 NM PB/NB PB/NB PB/NM PB/NS PB/NS 0/0 0/0 NS PM/NM PM/NM PM/NS PS/NS 0/0 NS/PS NS/PS ZP PS/NS PS/NS PS/NS 0/0 NS/PS NM/PM NM/PM ZN NP/PS NP/PS PS/NS 0/0 NS/PS PM/PS PM/PS PS NS/NS NS/NS 0/0 NS/PS NS/PS PM/PS PM/PS PM 0/0 0/0 NS/PS NM/PS NM/PM PB/NB PB/NB PB PS/0 PS/0 NM/PS NM/PM NM/PM PB/NL PB/NL 3.4 进气阀门开度预置
由于连续式风洞有大量用于气密封的空气往洞内泄露,因此常压时采用气密封漏气加排气阀调节的方式(进气阀关闭);而在增压时尤其高目标压力段,当压力欠压时,若仅通过气密封漏气的方式进行补气则速度过慢,需通过进气阀门对风洞进行补气,调试时为尽量减少气源压力波动对总压控制的影响,结合阀门性能曲线,减少气源消耗,进气阀以一个相对固定的进气量往风洞灌气,进气阀预置开度与气源压力和目标压力满足下列函数关系:
(8) (9) (10) 4 软件流程设计
程序设计时首先基于模块化编程方法,将阀门和泵进行组合以实现不同压力工况控制,风洞运行时,根据工况命令调用不同的模块,以满足多种工况下阀门开关组合的控制逻辑需求,并通过压力控制模块实现闭环控制,压力闭环由阀门系统PLC 315实现,风洞主控系统和阀门控制系统通过Profinet工业以太网进行通信,阀门控制系统软件流程见图 4。
5 试验验证
为测试控制算法性能,对变总压、变马赫数、变迎角等不同工况均进行了调试,迎角范围为-4°~14°,图 5和6列举了不同工况下的典型测试结果,其中图 5为风洞在常压和增压工况下变马赫数变迎角调试结果,图 6为负压20kPa变马赫数变迎角试验结果。
由图 5和6可见,趋于稳态后,常压100kPa时,总压最大值100.071,最小值99.923,平均值100.002,按
[18]计算控制精度为0.079%,增压200kPa时,总压最大值200.152,最小值199.83,平均值200.01,控制精度为0.09%;负压20kPa时,总压最大值20.028,最小值19.972,平均值19.998,控制精度为0.15%;具体测试结果见文献[19]。因此,总压控制精度在常压增压工况时优于0.1%,在负压工况时优于0.2%,达到国军标先进指标[20],图中马赫数上升过程的超调主要由于调节阀的响应速度较慢导致。6 结论
经过试验验证,可以得出以下结论:
(1) 利用分段变参数加模糊PID控制方法实现了0.6m风洞总压精确控制,总压控制精度优于设计指标要求,该控制方法对较宽压力范围或调节阀性能比较敏感的风洞压力控制是有效的;
(2) 如果进一步提高各调节阀的控制精度及响应速度,总压控制精度和稳定速度将进一步提升;
(3) 本文在0.6m风洞探讨的风洞总压控制策略可为后续大型风洞及其他连续式风洞的压力控制方法提供参考。
下一步还将重点对总压和马赫数之间,不同马赫数调节手段与总压控制之间的耦合关系等方面展开更深入的研究。
-
表 1 不同实验条件下的无量纲分离距离L′ 对比
Table 1 Comparison of dimensionless separation distance L´
实验条件 无量纲分离距离L′ δ=5° δ=10° δ=15° δ=20° 矩形平板模型 0.017 0.067 / 0.175 三角形平板模型,α =0° 0 0.067 0.108 0.142 三角形平板模型,α =10° 0 0.225 / 0.367 -
[1] LEWIS J E,KUBOTA T,LEES L. Experimental investigation of supersonic laminar, two-dimensional boundary-layer separation in a compression corner with and without cooling[J]. AIAA Journal,1968,6(1):7-14. doi: 10.2514/3.4434
[2] HOLDEN M S. Shock wave-turbulent boundary layer interaction in hypersonic flow[C]//Proc of the 10th Aerospace Sciences Meeting. 1977. doi: 10.2514/6.1972-74
[3] MALLINSON S G,GAI S L,MUDFORD N R. The interaction of a shock wave with a laminar boundary layer at a compression corner in high-enthalpy flows including real gas effects[J]. Journal of Fluid Mechanics,1997,342:1-35. doi: 10.1017/s0022112097005673
[4] LI S X, CHEN Y K. A study of interaction flow induced by 3D wedge[C]// Proceeding of the 3rd International Conference Fluid-Mechanics. 1997.
[5] CANDLER G, NOMPELIS I, HOLDEN M S. Computational analysis of hypersonic laminar viscous-inviscid interactions[R]. AIAA-2000-0532, 2000. doi: 10.2514/6.2000-532
[6] DOLLING D S. Fifty years of shock-wave/boundary-layer interaction research: What next?[J]. AIAA Journal,2001,39(8):1517-1531. doi: 10.2514/2.1476
[7] BABINSKY H, HARVEY J K. Shock wave-boundary-layer interac-tions[M]. Cambridge: Cambridge University Press, 2011. doi: 10.1017/cbo9780511842757
[8] 余平,段毅,尘军. 高超声速飞行的若干气动问题[J]. 航空学报,2015,36(1):7-23. YU P,DUAN Y,CHEN J. Some aerodynamics issue in hypersonic flight[J]. Acta Aeronautica et Astronautica Sinica,2015,36(1):7-23.
[9] DRUGUET M C,CANDLER G V,NOMPELIS I. Effects of numerics on Navier-Stokes computations of hypersonic double-cone flows[J]. AIAA Journal,2005,43(3):616-623. doi: 10.2514/1.6190
[10] 陈政, 倪招勇, 王军旗, 等. 高超声速层流分离流动模拟影响因素分析[C]// 全国第十五届分离流、旋涡和流动控制会议论文集. 2014. [11] 李艳丽,李素循. 高超声速绕钝舵层流干扰流场特性研究[J]. 宇航学报,2007,28(6):1472-1477. DOI: 10.3321/j.issn:1000-1328.2007.06.007 LI Y L,LI S X. Investigation of interactive hypersonic laminar flow over blunt fin[J]. Journal of Astronautics,2007,28(6):1472-1477. doi: 10.3321/j.issn:1000-1328.2007.06.007
[12] PRINCE S, VANNAHME M, STOLLERY J. Experiments on the hypersonic turbulent shock-wave/boundary layer interaction and the effects of surface roughness[R]. AIAA-99-0147, 1999. doi: 10.2514/6.1999-147
[13] 李素循,倪招勇. 高超声速层流干扰流场研究[J]. 宇航学报,2003,24(6):547-551, 573. DOI: 10.3321/j.issn:1000-1328.2003.06.001 LI S X,NI Z Y. Investigation of laminar interactive flowfield in hypersonic flow[J]. Journal of Astronautics,2003,24(6):547-551, 573. doi: 10.3321/j.issn:1000-1328.2003.06.001
[14] INGER G R. Scaling of incipient separation in high speed laminar flows[J]. AIAA Journal,1995,33(1):178-181. doi: 10.2514/3.12356
-
期刊类型引用(8)
1. 邓海均,熊波,罗新福,洪少尊,李强,王伟仲,刘俊. 跨声速风洞中使用短轴探管测量试验段核心流马赫数影响研究. 南京航空航天大学学报. 2022(02): 225-232 . 百度学术
2. 邓海均,熊波,罗新福,洪少尊,李腾骥,张利波,韩新峰. 连续式跨声速风洞短轴探管设计与试验. 航空动力学报. 2022(06): 1171-1180 . 百度学术
3. 郭良斌,余振东,吴永良,李宁. 暂冲实验气源两级调压阀串联调压的参数设计及测试. 机床与液压. 2021(07): 1-7+21 . 百度学术
4. 陈吉明,吴盛豪,廖达雄,裴海涛,吕金磊,熊波. 0.6 m连续式跨声速风洞流场品质改进试验研究. 南京航空航天大学学报. 2021(02): 236-242 . 百度学术
5. 陈海峰,黄叙辉,凌忠伟,马列波,秦建华. 0.6米暂冲式跨超声速风洞控制系统研制. 自动化与仪器仪表. 2021(04): 35-39 . 百度学术
6. 邓海均,熊波,罗新福,洪少尊,刘奇,蒲泓宇,刘为杰. 跨声速风洞试验段马赫数静压探针移测方法初探. 西北工业大学学报. 2021(04): 786-793 . 百度学术
7. 陈旦,杨孝松,李刚,郭守春,陈天毅. 连续式风洞总压和调节阀相关性研究及其应用. 西北工业大学学报. 2020(02): 325-332 . 百度学术
8. 陈旦,张永双,李刚,郭守春,沈牟. 一种连续式跨声速风洞总压控制方法设计. 实验流体力学. 2019(06): 65-71 . 本站查看
其他类型引用(2)