壁湍流摩阻时空信号无接触测量实验研究

陈纪仲, 田海平, 丁俊飞, 马国祯, 陈怡纯

陈纪仲, 田海平, 丁俊飞, 等. 壁湍流摩阻时空信号无接触测量实验研究[J]. 实验流体力学, 2024, 38(4): 66-74. DOI: 10.11729/syltlx20230109
引用本文: 陈纪仲, 田海平, 丁俊飞, 等. 壁湍流摩阻时空信号无接触测量实验研究[J]. 实验流体力学, 2024, 38(4): 66-74. DOI: 10.11729/syltlx20230109
CHEN J Z, TIAN H P, DING J F, et al. Experimental study on non-contact measurement of spatio-temporal signals in wall turbulence friction[J]. Journal of Experiments in Fluid Mechanics, 2024, 38(4): 66-74. DOI: 10.11729/syltlx20230109
Citation: CHEN J Z, TIAN H P, DING J F, et al. Experimental study on non-contact measurement of spatio-temporal signals in wall turbulence friction[J]. Journal of Experiments in Fluid Mechanics, 2024, 38(4): 66-74. DOI: 10.11729/syltlx20230109

壁湍流摩阻时空信号无接触测量实验研究

基金项目: 国家自然科学基金项目(12172242,12102284,12332017,12302289);山西省基础研究计划项目(20220302122106,20210302124262)
详细信息
    作者简介:

    陈纪仲: (1999—),男,山西运城人,硕士研究生。研究方向:实验流体力学,湍流结构及流动控制。E-mail:chenjizhong2022@163.com

    通讯作者:

    田海平: E-mail:tianhaiping@tyut.edu.cn

  • 中图分类号: O357.5+4

Experimental study on non-contact measurement of spatio-temporal signals in wall turbulence friction

  • 摘要:

    壁湍流摩擦阻力的产生机理与实验测量是壁湍流基础研究中的重要课题,关系到流动减阻发展和工程技术变革。在壁湍流复杂相干结构迁移演化过程中,会引起内外区强烈的动量交换,使黏性底层速度梯度呈现时空变化,即局部摩擦阻力与相干结构迁移演化呈现时空关联特性。目前鲜有实验测量技术能同时对摩擦阻力时空信号及与其关联的湍流结构发展演化过程进行精准、同步的测量。作为一种无接触的流场测量方法, PIV在摩擦阻力测量方面不断被寄予厚望。本文基于高时空分辨率2D2C TRPIV时序流场的粒子图像数据,开发了一种可调时空分辨率互相关算法(Adjustable Spatio-Temporal Resolution Cross Correlation, ASTRCC)。利用该算法计算了具有高时空分辨率的壁面摩擦阻力,探究了摩擦阻力的空间波动性与时间延迟性,并从时空演化角度探索了壁湍流相干结构与局部流场摩擦阻力的时空关联特性,验证了该算法的正确性。

    Abstract:

    The generation mechanism and experimental measurement of wall turbulence frictional resistance are critical topics in the fundamental research of wall turbulence, which are related to flow drag reduction and transformative engineering technologies. The migration and evolution of complex coherent structures in wall turbulence cause intense momentum exchange between the inner and outer regions. This leads to spatiotemporal variations in the velocity gradient of the viscous sublayer. Local frictional resistance is correlated with the spatiotemporal changes in these coherent structures. Few experimental measurement techniques currently exist that can accurately and simultaneously measure the spatio-temporal signals of the friction drag and the associated evolution of turbulent structures. As a non-contact flow field measurement method, the Particle Image Velocimetry (PIV) has been increasingly expected to excel in frictional resistance measurement. This paper developed an Adjustable Spatio-Temporal Resolution Cross Correlation (ASTRCC) algorithm, based on high spatiotemporal resolution 2D2C TRPIV time-series flow field particle image data. Wall frictional resistance with high spatiotemporal resolution was calculated using this algorithm. The spatial fluctuations and temporal delays of frictional resistance were investigated. From the perspective of spatiotemporal evolution, the spatiotemporal correlation between coherent structures in wall turbulence and local flow field frictional resistance was explored. The correctness of this method was validated.

  • 高超声速飞行器是指巡航Ma>5.0,以吸气式、组合动力发动机或火箭等为推进系统,在大气层及跨大气层中实现高超声速巡航的飞行器。超燃冲压发动机作为高超声速飞行器推进系统的最佳解决方案之一,逐渐发展为各个国家的重点研究项目[1-2]。然而超燃冲压发动机必须达到一定的飞行速度才能启动工作,否则将无法满足从起飞到高超声速巡航的宽速域飞行任务。于是,有研究者提出3种组合发动机推进系统的概念——涡轮基组合循环发动机(TBCC)、火箭基组合循环发动机(RBCC)及空气涡轮火箭/冲压组合循环发动机(ATR)。相较于火箭基组合循环发动机(RBCC),涡轮基组合循环发动机(TBCC)可以产生大一个数量级的比冲,且能使用普通燃料实现常规起降,大大提高了经济性。因此,TBCC组合发动机推进系统具有更好的应用前景,各国都相继开展了TBCC的相关研究计划,如美国的RTA计划[3]、日本的HYPR计划[4]和欧洲的LAPCAT计划[5]等。有关研究表明:在飞行马赫数Ma=6.0时,非对称膨胀喷管产生的推力占到推进系统总推力的70%[6]。另有研究表明:尾喷管的推力系数下降1%会造成安装净推力4%的下降[7]。因此,喷管的设计非常重要。

    近年来,TBCC组合发动机推进系统的研究热点在发动机总体和进气道方面,国内外对TBCC推进系统喷管的设计及性能研究较少。由于要满足飞行器和发动机一体化的要求,飞行器的后体表面一般情况下要与动力系统的喷管相融合,于是就产生了非对称形式的喷管膨胀型面,即所谓的单边膨胀喷管(SERN)[8]。Dusa[9]最先对TBCC尾喷管进行研究,提出喷管必须采用几何可调结构,以满足喷管在十分大的落压比范围内的正常工作,喷管在冲压模态工作时需要采取有效的冷却措施,以减小热负荷带来的负面影响。Gamble等[10]使用simMechanic和MATLAB/Simulink几何设计程序整合了并联式TBCC喷管几何结构和运动子系统,使用滑动及转动结构的4个机械部件,实现喷管的调节控制,但整个喷管可动部件较大,调节方式比较复杂。我国对TBCC发动机的研究起步较晚,近年来,在国家的大力支持下,北京航空航天大学[11-13]、南京航空航天大学[14-16]、西北工业大学[17-19]等院校和研究所都开展了相关工作,并取得了一些重要的进展。

    大多数研究者所设计的喷管长度与进口高度的比值(D/Hi)即长高比均在10左右。对排气系统而言,在满足进气道及燃烧室出口位置的限制下,要在飞行器后体的有限空间内完成满足性能要求的喷管气动方案设计,是一个较为严格且具有挑战性的命题。本文的研究对象是一种强几何约束下的组合发动机喷管,且飞行器总体设计要求喷管的长高比为5.28,设计难度明显增大。在借鉴上述科研成果,且在满足总体强几何约束条件下,对长高比为5.28的小长高比组合发动机喷管的气动设计展开了研究。采用喷管特征线法开展喷管二维型线设计,对喷管的设计点选取、三维侧向膨胀角、喷管双通道相对位置等影响因素进行了研究,给出了组合发动机喷管气动方案,并对其气动性能开展了数值计算及分析。

    本文所设计及研究的组合发动机喷管飞行马赫数范围为Ma=0~6.0,图 1为组合发动机喷管模态转换过程和冲压模态过程中的几何示意图。为满足几何尺寸要求,冲压喷管与涡轮喷管均为圆形进口,涡轮发动机喷管从进口到转轴处为面积逐渐减小的圆转方收缩段;为满足喉道可调,转轴到喉道位置处设计为二元等宽的收缩段,喉道后设计为二元等宽的扩张段。由于冲压喷管进口为超声速气流,相应地在冲压喷管的二元等宽扩张段前设计了一段具有等面积变化规律的圆转方过渡段,其中涡轮喷管收缩段型线皆采用五次方曲线法[20]设计。基于文献[21]用Fluent软件对非对称喷管内外流场的计算与验证,本文采用商用CFD软件Fluent对模型进行数值计算选取基于密度的隐式求解器,选择的计算模型为RNG、k-ε湍流模型,壁面函数设置为非平衡壁面函数,并采用变比热的理想气体模型,分子黏性系数采用Sutherland公式计算。图 2给出了喷管对称面上的网格图,整个流场采用结构化网格划分,总网格数约为300万。如图 3所示,喷管进口、出口以及延伸面上的边界条件分别取为压力进口、压力出口。为保证计算准确性,在喷管喉道、壁面附近、出现激波和剪切层等的流动区域内进行了网格加密。边界层第一层高度为0.05 mm,取12层,且壁面y+小于10。下面对网格无关性进行研究,选取3种数量的网格进行比较,数量分别为200万(粗网格)、300万(中等网格)和400万(细网格)。

    图  1  排气系统示意图
    Fig.  1  Geometry of the exhaust system
    图  2  典型工作状态时喷管网格图
    Fig.  2  Mesh distribution of the TBCC exhaust system
    图  3  排气系统网格图
    Fig.  3  Computational domain and mesh of the exhaust system

    表 1中可见,3种网格计算的结果非常一致,推力系数Cfx最大相差为0.4%,升力L最大相差为1%,俯仰力矩M最大相差0.57%,反应出其网格无关的特性。为保证求解精度与计算效率,将网格量控制在300万左右。

    表  1  不同网格量的喷管性能参数
    Table  1  Nozzle performance parameters with different mesh numbers
    Mesh numbers L/N M/(N·m) Cfx
    200万 2575.8 -2069.9 0.9132
    300万 2583.4 -2088.3 0.9163
    400万 2610.1 -2100.1 0.9200
    下载: 导出CSV 
    | 显示表格

    为了验证本文数值计算方法的准确性,选取NASA Spaid等[22]发表的喷管实验报告,进行相同状态下的数值计算,并与实验结果进行对比。验证算例模型具体的几何尺寸如图 4所示。喷管进口高H0=15.24 mm,喷管进口马赫数1.665,总压172 kPa,总温478 K;外部自由来流马赫数6.0,总压2520 kPa,总温478 K。图 5给出了按文献[22]的数值计算方法,对该喷管算例进行数值计算得到的喷管沿程压力分布曲线对比。从图中可以看出,数值计算得到的喷管沿程压力分布趋势与实验结果变化趋势一致,吻合较好,表明本文所采用的计算方法可以较为准确地模拟此类喷管流动,数值计算结果可信。

    图  4  Computational domain and mesh of the exhaust system
    Fig.  4  Scheme of the nozzle model
    图  5  壁面压力的实验与计算结果比较
    Fig.  5  Comparison between the experiment and numerical calculation of the wall pressure

    本文评价喷管的性能参数包括:轴向推力系数Cfx、升力L和俯仰力矩M

    轴向推力系数定义为轴向推力与等熵推力之比:

    (1)

    (2)

    (3)

    (4)

    其中, i为进口参数下标,mi为进口质量流量,k为比热容比,R为气体常数,Tt4为喷管进口总温,Pt4为喷管进口总压,NPR为落压比,Vix为进口轴向速度,p4为进口静压,Aix为进口面积在y方向的投影,fx为作用在喷管壁面上的相对压力积分在x方向的投影。如图 6所示,升力则是出口AB截面冲量在y负方向的分量,为整个流通通道壁面上的相对压力积分在y方向的投影,记为L。喷管壁面所受力对冲压发动机喷管进口中心点的力矩记为M

    图  6  喷管性能参数示意图
    Fig.  6  Schematic diagram of nozzle performance parameters

    图 7给出了推进系统总体对喷管的几何尺寸及位置约束条件。如图所示,冲压发动机喷管进口形状为圆形,直径H,涡轮发动机喷管进口形状为圆形,直径为1.44H。上、下喷管通道中心线相对距离为1.73H。涡轮及冲压发动机喷管进口水平相对距离为0.52H。喷管长度为5.28H

    图  7  喷管几何尺寸及位置约束示意图
    Fig.  7  The constraint of size and position

    在上述几何尺寸及相对位置的约束下,采用课题组前期研发的喷管特征线法设计程序[22]开展喷管二维型线设计,图 89给出了特征线法设计的Ma=3.0二元非对称喷管型面及流场图。常规的喷管长高比都在10左右,而本文研究的喷管的长高比仅为5.28。对于排气系统而言,在满足进气道及燃烧室出口位置的限制下,要在飞行器后体的有限空间内完成满足性能要求的喷管气动方案设计,是一个较为严格且具有挑战性的命题,针对上述小长高比喷管设计问题,本文给出了一套喷管设计流程,先后对喷管的设计点选取、三维侧向膨胀角、喷管双通道相对位置开展参数化研究。先通过喷管设计点选取的研究,确定冲压发动机喷管的二维气动方案,接着在此基础上,通过三维侧向膨胀角与下壁面截短的研究,确定冲压发动机喷管三维气动方案,最后通过喷管双通道相对位置研究确定组合发动机喷管的三维气动方案。此设计流程能够在几何受限情况下同时兼顾空间的有效利用与喷管的气动性能。

    图  8  特征线法设计的二元非对称喷管型面
    Fig.  8  The nozzle designed by MOC
    图  9  二元非对称喷管流场马赫数等值图
    Fig.  9  Contours of the nozzle

    当飞行器在低马赫数飞行时,排气系统的设计膨胀面积比较大,而工作压比很小,内部流动处于严重的过膨胀状态,推力系数普遍较低。由于低压比下喷管性能恶化程度与喷管设计点有关,因此可以考虑减小喷管的设计点来改善其在低马赫数飞行时的性能;当喷管的设计点降低后,设计面积比也相应地减小,在低落压比工作时,过膨胀程度就相应减小。基于这种考虑,选取了Ma=3.0、3.5、4.0和5.0的冲压发动机喷管进口参数(见表 2)作为喷管的设计点参数,开展喷管气动型面设计。图 10给出了各设计点下原始长度和上下壁面等比例截短后的二维喷管型线。

    表  2  冲压发动机喷管进口参数
    Table  2  The parameters of ramjet nozzle inlet
    Ma h/km pt/Pa Tt/K
    3.0 18.0 138 319 1432
    3.5 20.0 177 797 1496
    4.0 21.9 164 240 2195
    5.0 24.8 202 909 2665
    下载: 导出CSV 
    | 显示表格
    图  10  不同设计点下喷管上下型线对比
    Fig.  10  The nozzle model with different design Mach numbers

    图 1112给出了喷管的推力系数和升力随飞行马赫数的变化曲线,其中横坐标为飞行马赫数,纵坐标为喷管的推力系数和升力,不同的曲线代表按不同设计点获取的喷管方案。可以看出,在相同的设计马赫数下,随着飞行马赫数的增加,喷管推力系数先增加后缓慢减小,升力则逐渐增大。在相同的飞行马赫数下,随着设计点马赫数的增大,推力系数和升力均减小,且不同的飞行马赫数下减小幅度不同,在飞行Ma=3.0时,设计点Ma=3.0的喷管比设计点Ma=5.0喷管的推力系数提高了约10%,而改变设计点马赫数对喷管在巡航点Ma=6.0的轴向推力系数影响不大,设计点为Ma=3.0和5.0时,轴向推力系数相差最大2.5%。设计点Ma=3.0和3.5下的喷管推力性能非常接近,且设计点Ma=3.0下的喷管性能相对更为优越。综合考虑,确定Ma=3.0作为组合发动机喷管的设计点。

    图  11  不同设计点的喷管推力系数变化图
    Fig.  11  The axial thrust coefficient with different design Mach numbers
    图  12  不同设计点喷管升力变化图
    Fig.  12  The lift force with different design Mach numbers

    在喷管设计点二维气动方案基础上,本节针对冲压发动机喷管的三维气动方案进行探究,重点研究侧向膨胀角度对喷管气动性能的影响。由于二维气动方案的确定实则是对喷管出口高度形成了限制,喷管的出口面积并未确定,而喷管的进出口面积膨胀比是影响喷管推力特性的重要因素,因此研究过程中,在保证喷管出口高度不变的前提下,首先将二维喷管型线直接侧向拉伸,形成进出口等宽的二元喷管作为基准构型,如表 3所示,并再选取2个侧向膨胀角,保证喷管上下壁面长度相同,进行数值计算及性能分析。不同侧向膨胀角的喷管构型的俯视图如图 13所示。

    表  3  喷管构型的侧向膨胀角
    Table  3  The lateral expansion angle of the nozzle
    Model Angle Width
    A 1.50 Ht
    B 2.23° 2.06 Ht
    C 6.06° 3.03 Ht
    下载: 导出CSV 
    | 显示表格
    图  13  喷管的侧向膨胀角示意图
    Fig.  13  The contour of the lateral expansion angle

    图 14为3种喷管构型的推力系数随飞行马赫数的变化曲线,其中横坐标为飞行马赫数,纵坐标为喷管的推力系数。从图中可以看出,在相同的侧向膨胀角度下,随着飞行马赫数的增加,喷管推力系数先增加后缓慢减小,在低马赫数飞行条件下,随着侧向膨胀角度的增加,喷管的推力系数逐渐降低,反之,在高马赫数飞行条件下,喷管的推力系数随侧向膨胀角度增加而增加。构型A在Ma=3.0~4.0时具有较高的推力系数,在Ma=3.0时比构型C推力系数高约8.7%,而在Ma=5.0~6.0时推力系数有所下降,在Ma=6.0时相较于构型C降低了约2%,这主要是由于喷管进出口面积比小导致欠膨胀,从而推力性能下降。相对而言,带有一定侧向膨胀角度的构型B和C,在Ma=5.0~6.0高马赫数飞行条件下具有较高的推力系数,但低马赫数Ma=3.0~4.0性能下降明显,对于构型C而言,Ma=3.0时推力系数相较于Ma=6.0时下降幅度超10%,喷管过膨胀的问题仍十分突出。由此可见,在保持喷管出口高度不变的情况下,随着喷管侧向膨胀角的增大,高马赫数性能有所提高,低马赫数性能下降明显,尤其在飞行条件Ma=3.0时,构型C相较构型A的推力系数下降幅度约8.7%,但高马赫数时提升幅度仅为2%。从总体考虑,并优先保证Ma=3.0~4.0的喷管性能,故Model A为最佳方案。将其作为组合发动机喷管的基准型面。

    图  14  不同构型喷管的推力系数
    Fig.  14  The axial thrust coefficient of the models

    图 15所示,为实现巡航点Ma=6.0正升力,现初步将下壁面长度由D截短为0.38 D图 16给出了2.4节中构型A喷管的推力系数和升力随飞行马赫数的变化曲线,其中横坐标为飞行马赫数,左右纵坐标分别为喷管的推力系数和升力。由图 16可知,喷管的升力为负。喷管的升力主要来自于喷管上下壁面的压力差,其值为负的主要原因是喷管上壁面的压力在竖直方向上小于下壁面的压力。为实现正升力,需减小下壁面带来的负压力,最为直接的方式就是将下壁面进一步截短。

    图  15  喷管下壁面截短对比
    Fig.  15  The nozzle after cutting wall below
    图  16  截短前后喷管性能对比
    Fig.  16  The nozzle performance parameters after cutting wall below

    图 16中截短后的喷管性能曲线可以看出,对下壁面进行截短可有效提高喷管的升力。下壁面的进一步截短实则改变了喷管出口面积,影响了喷管的面积膨胀比,下壁面的截短同时也损失了壁面上的压力,对喷管的性能产生了直接的影响,使推力系数略有下降,在Ma=6.0时,下降幅度达到最大,约为1.4%。但从工程应用的角度来看,升力提升的效果是十分可观的,在Ma=3.5~6.0时,喷管实现了正升力,即便低马赫数(Ma=3.0)仍有负升力现象,但对于巡航点Ma=6.0而言,升力提高5倍,已达到要求。

    总体要求涡轮与冲压喷管的进口形状为圆形,且2个喷管进口中心线相对距离为1.73H。为保证模态转换过程喷管几何可调,将喷管前段设计为圆转方过渡段,在中间档板起始位置过渡为二元等宽喷管以便进行几何调节。

    图 17所示,定义中间调节板的起始位置在a点。由于本文设计总体要求喷管长度与高度的比值为5.28,远低于常规的喷管长高比几何约束,导致喷管气动性能降低。为使涡轮发动机喷管获得更大的膨胀比,即增加共用的喷管膨胀壁面长度,现将涡轮发动机喷管圆转方收缩段出口的中心线竖直向下偏置一定距离,使中间档板的转轴前移,即改变Ha的大小,所以Ha是喷管设计中的一个重要参数。由于组合发动机在转级点附近需要较高的推力性能来保证转级的顺利,因此可以通过获得转级点最优的推力性能来寻求最佳的a点位置。选取5种不同的偏置距离Ha,并选取相同的分流板角度,保持流量相同,研究涡轮发动机喷管相对位置对组合发动机喷管性能的影响。选取模态转换Ma=3.5进行性能评估及对比分析。

    图  17  组合发动机喷管偏置距离示意图
    Fig.  17  The schematic diagram of the offset distance

    图 18给出了冲压发动机喷管、涡轮发动机喷管、组合发动机喷管的推力系数随偏置距离的变化曲线,横坐标为偏置距离,纵坐标为喷管的推力系数。可以看出,冲压发动机喷管推力系数随着偏置距离的增加而降低,在Ha=0.3Dt时达到最低。而随着偏置距离的增加,涡轮发动机及组合发动机喷管总的推力系数呈现先增后减的趋势,在Ha=0.267Dt时,推力系数达到最大值,相较于Ha=0.3Dt时组合发动机喷管推力系数提升了约1%。由于喷管受几何尺寸限制,偏置距离改变范围较小,相应地造成性能参数的变化幅度较小,但仍可得出喷管推力性能在喷管双通道不同的位置布局情况下的变化趋势。可见偏置距离的增加,造成了分流板转轴的前移,增加了涡轮发动机喷管的面积膨胀比,在一定范围内可以提高涡轮发动机喷管的推力性能,同时势必会造成冲压发动机出口面积的减小,冲压喷管推力性能下降。而就整个组合发动机喷管的推力性能而言,存在一个合适的偏移距离Ha,使得喷管总体推力性能最优。对于本次评估的5个偏移方案,Ha=0.267Dt方案相对较优。因此,按此几何尺寸确定组合发动机可调喷管的三维气动方案, 图 19给出了最终确定的组合发动机喷管的三维气动方案。

    图  18  不同偏置距离下喷管推力系数
    Fig.  18  The axial thrust coefficient of offset distances
    图  19  组合发动机喷管三维气动方案示意图
    Fig.  19  Three-dimensional of the TBCC exhaust system

    2.6节确定了组合发动机喷管的最终气动方案,本节通过数值计算方法对其开展详细的流场计算和性能分析,讨论了推力系数、升力、俯仰力矩随飞行马赫数的变化关系,并对飞行器在冲压模态和模态转换时组合发动机喷管的性能进行分析。

    表 4给出了Ma=3.0~6.0时组合发动机喷管的推力系数、升力和俯仰力矩。可以看出,随着飞行马赫数的增加,升力和俯仰力矩逐渐增加,而喷管推力系数先增加后缓慢减小。由于设计点为Ma=3.0,巡航点Ma=6.0的推力性能略有下降,但Ma=3.0~5.0时推力性能都能达到0.920以上。

    表  4  冲压模态时喷管性能参数
    Table  4  Nozzle performance parameters at supersonic mode
    Ma h/km L/N M/(N·m) Cfx
    3.0 18.0 -452.9 -136.1 0.9204
    3.5 20.0 120.2 561.0 0.9508
    4.0 21.9 223.2 571.4 0.9416
    5.0 24.8 746.6 1214.4 0.9259
    6.0 26.0 1131.7 1761.2 0.9184
    下载: 导出CSV 
    | 显示表格

    图 20给出了在冲压模态过程中马赫数为3.0、3.5、4.0、5.0和6.0这5个典型工作状态下喷管对称面的马赫数分布图。从图中可以看出, 冲压模态过程中排气系统的流场结构较为相似,相较于转换模态过程中的流场分布,冲压模态过程中流场结构简单,在低马赫数下冲压发动机喷管处于过膨胀状态。由于中间分流板的结构设计,在喷管上壁面有一处向上偏转,由于偏转处面积比的突增,在此处形成膨胀波,气流经过膨胀波后偏转,在下游又形成一道压缩波。气流经过激波后马赫数下降0.2左右。在冲压模态过程中,排气系统流场变化均匀,没有出现明显的流动分离现象。

    图  20  冲压模态时喷管对称面马赫数等值线图
    Fig.  20  Contours of nozzle Mach number at supersonic mode

    根据飞行器总体要求,喷管模态转换工作点为Ma=3.0及3.5,要求喷管能够兼顾2个模态转换点的性能。当飞行高度为20 km、飞行马赫数为3.5时,组合发动机处于共同工作状态,开始从涡轮模态向冲压模态转变,转换过程中来流马赫数保持不变,涡轮发动机逐渐关闭,冲压发动机开启工作。本节选取Ma=3.5模态转换过程进行气动性能分析,针对模态转换过程中的5个典型工作点进行数值计算。表 5给出了模态转换过程中,组合发动机喷管的推力系数、升力和俯仰力矩。可以看出,随着中间分流板的偏转,涡轮喷管通道逐渐关闭,推力系数逐渐增加,俯仰力矩逐渐增加,而升力逐渐降低。这是由于模态转换过程中,2个喷管出口气流参数差异较大,2股气流相互干扰严重,最终影响到喷管的性能参数。

    表  5  模态转换时喷管性能参数
    Table  5  Nozzle performance parameters at transition mode
    Transition mode L/N M/(N·m) Cfx
    1 2583.5 -2088.3 0.9163
    2 2413.6 -1861.9 0.9190
    3 2178.1 -1739.3 0.9229
    4 1599.6 -1693.3 0.9344
    5 801.8 -337.1 0.9515
    下载: 导出CSV 
    | 显示表格

    图 21为涡轮发动机和冲压发动机喷管共同工作的马赫数等值图。模态转换过程中,中间分流板逐渐向上转动关闭涡轮通道,喉道高度逐渐减小。涡轮喷管的进口总压逐渐降低;而冲压通道的入口气流条件保持不变。

    图  21  模态转换时喷管对称面马赫数等值图
    Fig.  21  Contours of nozzle Mach number at transition mode

    模态转换过程中整个流场中存在复杂的激波、压缩波和滑移面现象。涡轮喷管通道内的气流在喷管出口上壁面受到远方来流的作用,在涡轮通道气流一侧产生了一道向下偏转的激波。由于2股气流的速度、密度等参数不同,两者之间形成一道滑移面(Shear layer 1)。由于中间分流板末尾起着涡轮发动机喷管喉道的作用,高温高压燃气经过分流板后继续向下膨胀,并与冲压通道的气流相互作用,形成一道滑移面(Shear layer 2)。由于气流在滑移面处产生偏转,由此形成一道向上的压缩波,并与上壁面形成的激波交汇后继续向上延伸。

    由于冲压喷管通道出口处气流与远方来流相互作用,在冲压通道下壁面出口处形成了膨胀波;2股气流之间也形成一道滑移面(Shear layer 3)。由于分流板作用,在冲压通道内形成一道向下的弱激波,并与下方的膨胀波相交后,继续向下偏转。自由来流在下唇板尾缘与冲压通道内气流相遇,在滑移面下侧产生两道激波。

    随着分流板不断向上偏转,涡轮通道出口气流角度不断增加,涡轮气流逐渐向上偏转,滑移面(Shear layer 2)也随之向涡轮通道侧偏移,而冲压通道固定不变,滑移面(Shear layer 3)变化不明显。由于分流板的向上偏转,在分流板末端又产生一道向上的激波,并与涡轮喷管上壁面相交,产生一道反射激波。

    本文对强几何约束下的小长高比组合发动机可调喷管的气动设计进行了研究,并对Ma=3.0~6.0范围内的喷管气动性能开展研究,得出以下结论:

    (1) 对于强几何约束下的小长高比喷管设计,可通过降低设计点马赫数的方法改善组合发动机喷管在低马赫数飞行时的气动性能,可避免喷管由于设计膨胀面积比较大而工作压比较小,而出现严重过膨胀现象。

    (2) 就本文研究的组合发动机喷管而言,在出口高度相同的情况下,侧向膨胀角度为0°时气动性能相对较优。

    (3) 涡轮/冲压发动机喷管出口相对位置对并联布局的组合动力喷管转级点气动性能影响较大,就本文研究而言,可通过收缩段的偏移实现2个喷管流道设计膨胀比的调节,当Ha=0.267Dt时,组合发动机喷管气动性能相对较优。

    (4) 并联式组合动力喷管的双通道气流在共同工作阶段的涡轮出口处相互干扰严重,整个过程中推力系数约为0.920;而冲压模态过程中,并无流动分离现象,流场结构相对较好,推力系数不低于0.918。

  • 图  1   实验布置示意图

    Fig.  1   Schematic of the experimental setup

    图  2   ASTRCC方法示意图

    Fig.  2   Schematic diagram of the ASTRCC method

    图  3   瞬时流向速度的线性拟合及残差情况

    Fig.  3   Linear fitting and residual of instantaneous flow velocity

    图  4   速度剖面示意图

    Fig.  4   Schematic diagram of velocity profile

    图  5   互相关分析图

    Fig.  5   Cross correlation analysis diagram

    图  6   近壁流场扫掠与喷射示意图

    Fig.  6   A schematic sketch to show the "sweeping" and "ejection" in the near-wall flow field

    图  7   摩阻与结构关联示意图

    Fig.  7   Schematic diagram of the correlation between friction stress and structure

    图  8   u′v′散点图

    Fig.  8   u′v′ scatter diagram

    图  9   流向空间摩阻时间信号频谱图

    Fig.  9   Flow direction space friction stress and time signal spectrum diagram

    图  10   不同流向空间摩阻时间序列

    Fig.  10   Time series of spatial friction stress in different flow directions

    图  11   摩阻时间信号相关性分析示意图

    Fig.  11   Friction stress time signal correlation analysis diagram

    图  12   流向速度与摩阻互相关结果对比图

    Fig.  12   Flow velocity and friction stress cross-correlation diagram

    图  13   时空FTLE与摩阻相关性分析

    Fig.  13   Analysis of the correlation between spatiotemporal FTLE and frictional stress

    表  1   湍流边界层基本流动相关参数

    Table  1   Basic flow related parameters of turbulent boundary layer

    参数
    u/(mm﹒s−1) 158
    边界层厚度δ/mm 33.576
    动量损失厚度θ/mm 3.26
    统计平均的壁面摩擦速度uτ/(mm﹒s−1) 8.6
    基于动量损失的雷诺数${{Re} _\theta }$ 428
    基于边界层位移的雷诺数${{Re} _{{\delta ^ * }}}$ 4411
    基于摩擦速度的雷诺数${{Re} _\tau }$ 240
    下载: 导出CSV
  • [1] 许春晓. 壁湍流相干结构和减阻控制机理[J]. 力学进展, 2015, 45(1): 111–140. DOI: 10.6052/1000-0992-15-006

    XU C X. Coherent structures and drag-reduction mechanism in wall turbulence[J]. Advances in Mechanics, 2015, 45(1): 111–140. doi: 10.6052/1000-0992-15-006

    [2] 王轩, 范子椰, 陈乐天, 等. 流向凹曲率壁面湍流边界层的TRPIV实验研究[J]. 实验流体力学, 2022, 36(6): 1–9. DOI: 10.11729/syltlx20210084

    WANG X, FAN Z Y, CHEN L T, et al. Experimental study of TRPIV for turbulent boundary layer of longitudinal concave curvature wall[J]. Journal of Experiments in Fluid Mechanics, 2022, 36(6): 1–9. doi: 10.11729/syltlx20210084

    [3]

    TANG Z Q, MA X Y, JIANG N, et al. Local dynamic perturbation effects on the scale interactions in wall turbulence[J]. Journal of Turbulence, 2021, 22(3): 208–230. doi: 10.1080/14685248.2020.1864388

    [4]

    TANG Z Q, JIANG N. The effect of a synthetic input on small-scale intermittent bursting events in near-wall turbulence[J]. Physics of Fluids, 2020, 32(1): 015110. doi: 10.1063/1.5129042

    [5] 高南, 刘玄鹤. 实用化壁面切应力测量技术的综述与展望[J]. 空气动力学学报, 2023, 41(3): 1–24. DOI: 10.7638/kqdlxxb-2021.0450

    GAO N, LIU X H. A review of wall-shear-stress measurement techniques for practical applictions[J]. Acta Aerodynamica Sinica, 2023, 41(3): 1–24. doi: 10.7638/kqdlxxb-2021.0450

    [6] 严宇超, 姜澄宇, 马炳和, 等. 壁面剪应力标定方法研究综述[J]. 实验流体力学, 2017, 31(2): 20–25. DOI: 10.11729/syltlx20170007

    YAN Y C, JIANG C Y, MA B H, et al. Review of the calibration methods and devices for wall shear stress[J]. Journal of Experiments in Fluid Mechanics, 2017, 31(2): 20–25. doi: 10.11729/syltlx20170007

    [7] 赵荣娟, 吕治国, 黄军, 等. 基于压电敏感元件的摩阻天平设计[J]. 空气动力学学报, 2018, 36(4): 555–560. DOI: 10.7638/kqdlxxb-2016.0112

    ZHAO R J, LÜ Z G, HUANG J, et al. Design of skin friction balance based on piezoelectric ceramics[J]. Acta Aerodynamica Sinica, 2018, 36(4): 555–560. doi: 10.7638/kqdlxxb-2016.0112

    [8]

    BAARS W J, SQUIRE D T, TALLURU K M, et al. Wall-drag measurements of smooth- and rough-wall turbulent boundary layers using a floating element[J]. Experiments in Fluids, 2016, 57(5): 90. doi: 10.1007/s00348-016-2168-y

    [9]

    DRIVER D M. Application of oil-film interferometry skin-friction measurement to large wind tunnels[J]. Experiments in Fluids, 2003, 34(6): 717–725. doi: 10.1007/s00348-003-0613-1

    [10]

    GUBIAN P A, STOKER J, MEDVESCEK J, et al. Evolution of wall shear stress with Reynolds number in fully developed turbulent channel flow experiments[J]. Physical Review Fluids, 2019, 4(7): 074606. doi: 10.1103/PhysRevFluids.4.074606

    [11] 张伟, 葛耀君, 杨泳昕. 粒子图像测速技术互相关算法研究进展. [J] 力学进展, 2007, 37(3): 443-452.

    ZHANG W, GE Y J, YANG Y X. Research progress of cross correlation algorithms in particle image velocimetry[J]. Advances in Mechanics, 2007, 37(3): 443-452.

    [12] 王康俊, 白建侠, 唐湛棋, 等. 用平均速度剖面法测量湍流边界层壁面摩擦速度的对比研究[J]. 实验力学, 2019, 34(2): 209–216. DOI: 10.7520/1001-4888-17-190

    WANG K J, BAI J X, TANG Z Q, et al. Comparative study of turbulent boundary layer wall friction velocity measured by average velocity profile method[J]. Journal of Experimental Mechanics, 2019, 34(2): 209–216. doi: 10.7520/1001-4888-17-190

    [13]

    ZHU Y D, JIANG X Y, ZHANG Y C, et al. Iterative PIV interrogation for complex wall-bounded flows[J]. Measurement Science and Technology, 2019, 30(9): 095302. doi: 10.1088/1361-6501/ab20f4

    [14]

    NGUYEN C V, WELLS J C. Direct measurement of fluid velocity gradients at a wall by PIV image processing with stereo reconstruction[J]. Journal of Visualization, 2006, 9(2): 199–208. doi: 10.1007/BF03181763

    [15]

    KÄHLER C J, SCHOLZ U, ORTMANNS J. Wall-shear-stress and near-wall turbulence measurements up to single pixel resolution by means of long-distance micro-PIV[J]. Experiments in Fluids, 2006, 41(2): 327–341. doi: 10.1007/s00348-006-0167-0

    [16] 许德辰, 张悦, 刘欣乐, 等. 基于粒子追踪测速的壁面摩擦应力测量[J]. 实验流体力学, 2022, 36(2): 131–138. DOI: 10.11729/syltlx20210156

    XU D C, ZHANG Y, LIU X L , et al. Measurement of wall-shear stress via micro-particle tracking velocimetry[J]. Journal of Experiments in Fluid Mechanics, 2022, 36(2): 131–138. doi: 10.11729/syltlx20210156

    [17]

    WILLERT C E. High-speed particle image velocimetry for the efficient measurement of turbulence statistics[J]. Experiments in Fluids, 2015, 56(1): 17. doi: 10.1007/s00348-014-1892-4

    [18]

    SHEN J Q, PAN C, WANG J J. Accurate measurement of wall skin friction by single-pixel ensemble correlation[J]. Science China Physics, Mechanics & Astronomy, 2014, 57(7): 1352-1362. doi: 10.1007/s11433-014-5462-9

    [19]

    OKAMOTO K, NISHIO S, SAGA T, et al. Standard images for particle-image velocimetry[J]. Measurement Science and Technology, 2000, 11(6): 685–691. doi: 10.1088/0957-0233/11/6/311

    [20]

    NOBACH H, HONKANEN M. Two-dimensional Gaussian regression for sub-pixel displacement estimation in particle image velocimetry or particle position estimation in particle tracking velocimetry[J]. Experiments in Fluids, 2005, 38(4): 511–515. doi: 10.1007/s00348-005-0942-3

    [21] 申俊琦, 王建杰, 潘翀. 平板湍流边界层瞬时摩擦阻力的光学测量和统计分析[J]. 气体物理, 2020, 5(5): 13–23. DOI: 10.19527/j.cnki.2096-1642.0873

    SHEN J Q, WANG J J, PAN C. Optical measurement and statistical analysis of instantaneous wall-shear stress in a turbulent boundary layer[J]. Physics of Gases, 2020, 5(5): 13–23. doi: 10.19527/j.cnki.2096-1642.0873

    [22]

    MARUSIC I, ADRIAN R J. The eddies and scales of wall turbulence[M]//Ten Chapters in Turbulence. Cambridge: Cambridge University Press, 2012: 176-220. doi: 10.1017/cbo9781139032810.006

    [23] 樊星, 姜楠. 用平均速度剖面法测量壁湍流摩擦阻力[J]. 力学与实践, 2005, 27(1): 28–30. DOI: 10.3969/j.issn.1000-0879.2005.01.007

    FAN X, JIANG N. Skin friction measurement in turbulent boundary layer by mean velocity profile method[J]. Mechanics in Engineering, 2005, 27(1): 28–30. doi: 10.3969/j.issn.1000-0879.2005.01.007

    [24]

    WANG J J, PAN C, WANG J J. Characteristics of fluctuating wall-shear stress in a turbulent boundary layer at low-to-moderate Reynolds number[J]. Physical Review Fluids, 2020, 5(7): 074605. doi: 10.1103/PhysRevFluids.5.074605

    [25]

    PAN C, KWON Y. Extremely high wall-shear stress events in a turbulent boundary layer[J]. Journal of Physics: Conference Series, 2018, 1001: 012004. doi: 10.1088/1742-6596/1001/1/012004

    [26]

    HERPIN S, STANISLAS M, FOUCAUT J M, et al. Influence of the Reynolds number on the vortical structures in the logarithmic region of turbulent boundary layers[J]. Journal of Fluid Mechanics, 2013, 716: 5–50. doi: 10.1017/jfm.2012.491

    [27] 王康俊. 湍流边界层相干结构特性及其主动控制减阻机理的实验研究[D]. 天津: 天津大学, 2021.

    WANG K J. Experimental study on coherent structure characteristics of turbulent boundary layer and its active control drag reduction mechanism[D]. Tianjin: Tianjin University, 2021.

    [28]

    TIAN H P, YI X R, XU F, et al. Lagrangian-based spatial-temporal topological study on the evolution and migration of coherent structures in wall turbulence[J]. Acta Mechanica Sinica, 2022, 38(1): 321465. doi: 10.1007/s10409-021-09006-1

    [29] 李鹏, 单桂华, 迟学斌. 时变流场的有限时间李雅普诺夫指数(FTLE)并行算法研究[J]. 科研信息化技术与应用, 2014, 5(2): 43–50. DOI: 10.11871/j.issn.1674-9480.2014.02.005

    LI P, SHAN G H, CHI X B. Parallel algorithm of FTLE computation for time-varying flow field[J]. e-Science Technology & Application, 2014, 5(2): 43–50. doi: 10.11871/j.issn.1674-9480.2014.02.005

图(13)  /  表(1)
计量
  • 文章访问数:  415
  • HTML全文浏览量:  17
  • PDF下载量:  69
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-08-20
  • 修回日期:  2023-10-02
  • 录用日期:  2023-10-30
  • 刊出日期:  2024-08-24

目录

/

返回文章
返回
x 关闭 永久关闭