Open Access
Issue
JNWPU
Volume 38, Number 3, June 2020
Page(s) 649 - 656
DOI https://doi.org/10.1051/jnwpu/20203830649
Published online 06 August 2020

© 2019 Journal of Northwestern Polytechnical University. All rights reserved.

Licence Creative Commons
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

混合层广泛存在于亚燃和超燃冲压发动机燃烧室内,实现燃料与空气间的掺混及燃烧组织。其基本结构为平板将2股相互平行的不同或同种类的流体隔开,在隔板下游形成混合层。根据上下2股流流动速度的大小关系,可以分为亚/亚声速混合层、亚/超声速混合层和超/超声速混合层。通常,下标“1”表示高速侧流体参数;下标“2”表示低速侧流体参数。

针对自由混合层已开展大量的深入研究,获得表征自由剪切层的物理量,如剪切层的增长率、卷吸率等[1-5]。随着自由流马赫数的增加,逐渐从亚/亚声速混合层过渡至亚/超声速或超/超声速混合层,二维剪切层增长率降低,压缩性对剪切层的影响逐渐显现。寻找合适的参数表征压缩性的影响成为研究重点。对流马赫数被证实为合适的参数并得到关于该参数的修正关系[6-8]。剪切层内部存在大涡拟序结构,在卷吸和混合过程中起到重要作用[9-11]。

混合层的流动特征主要包括速度分布、压力分布、混合层厚度等。针对不可压低速混合层、射流剪切层、边界层等各类剪切层,已开展较为丰富的研究。研究发现不可压剪切层内的无量纲速度分布具有相似性,各流向速度具有一普遍的无因次分布[12],无量纲速度呈现误差函数分布[13]。对于高速可压剪切层,也发现在火箭轴对称超声速燃气射流基本段内,不同流向位置处速度和动压均落在一无因次分布曲线上[14]。

以往研究多采用纹影、阴影、NPLS等技术对混合层进行可视化,并确定混合层空间发展情况,但是较难定量分析混合层内速度分布规律。传统接触式压力探针可以完成流场的定量分析,但会干扰流场结构。粒子图像测速技术PIV的发展为开展混合层内速度测量提供可能。大量学者开展PIV测量亚声速与超声速气流速度的性能研究,包括PIV粒子要求、粒子特性、投放方法等[15-19]。赵玉新[4]与何霖[20]分别采用PIV技术测量超-超声速混合层和超声速平板边界层,获得速度分布特征。由此说明,PIV对混合层流场的测量具有可行性。然而,对亚-超混合层的PIV测量研究还相对较少。

因此,本文采用粒子图像测速技术PIV对亚-超声速混合层二维流场进行测量,并结合基于雷诺平均方程(RANS)的稳态数值模拟研究,获得了混合层的流场特征,包括速度分布、压力分布和混合层厚度。

1 研究方法

1.1 试验方法

图 1为试验系统示意图,包含2股平行流道。在试验段(编号7)形成混合层(编号9)前,2股流道相互隔开;试验段出口与真空罐(编号11)及真空泵(编号12)相连,模拟试验段内低于大气压力环境。在试验段上游,分别测量两股流体的总压(编号6)、总温(编号5)和流量(编号4)。

在超声速气流侧的壁面沿流向x=20~230 mm布置16个静压测点,同时,在隔板末端上游16 mm处,测量2股流体的静压。压力由精度等级为0.075%的单晶硅压力变送器进行测量;来流温度采用精度为0.1级K型热电偶测量;流量由孔板流量计测量得到。

PIV系统如图 2所示,主要由Nd: YAG双脉冲激光器、片光源组件、同步控制器、CCD相机和计算机组成。

其中,激光器最大激光功率为200 mJ/脉冲(532 nm),每次脉冲时间持续6~8 ns,最大工作频率15 Hz;片光源名义厚度为1 mm;CCD相机的分辨率为2 048×2 048,12 bit灰度等级,每个像素大小为(7.4×7.4)μm。整套系统由MicroVec软件进行控制并对拍摄结果进行后处理。示踪粒子为氧化铝粉末,名义粒子直径10 μm。试验前,将粒子放入110℃的恒温烘箱中烘干。粒子投放分别位于混合层初始形成处上游2 m和1.5 m处。

试验参数如表 1所示。超声速气流马赫数Ma1为1.32, 速度U1为396 m/s; 亚声速气流马赫数Ma2为0.11, 速度U2为40 m/s; 2股流体的速度差ΔU为356 m/s; 超声速气流的静压p1为78.8 kPa; 亚声速气流的静压p2为70.6 kPa; 2股流体的静压比p1/p2为1.11, 超声速气流为轻度欠膨胀状态; 2股流体的总温T1tT2t均为303 K; 试验段出口背压pe为75 kPa。

表1

试验参数

thumbnail 图1

试验系统原理

thumbnail 图2

PIV系统及试验段图片

1.2 数值计算方法

1.2.1 计算域及网格

计算域模型如图 3a)所示, 为二维物理模型。隔板厚度5 mm, 亚声速气流进口高度115 mm, 超声速气流进口高度为16 mm, 混合层流向发展长度为900 mm。采用结构化四边形网格并开展网格独立性验证。最终, 计算网格如图 3b)所示, 网格单元数为95 400, 混合层区域最小网格尺寸为0.1 mm。边界条件如图 3b)所示, 进口均采用压力进口, 出口为压力出口; 所有固体壁面采用无滑移壁面。

thumbnail 图3

计算域与网格

1.2.2 湍流模型

采用标准k-ω湍流模型。该模型包含了对低雷诺数影响、压缩性和剪切流扩展等的修正。输运方程为[21]

同时, 考虑剪切流修正和压缩性影响。在kω的耗散项中引入压缩性修正, 压缩性方程为

式中:

1.2.3 数值方法验证

图 4为本文数值计算和Goebel[22]试验无量纲速度分布及Hall和Dimotakis[10]获得的总压分布对比。速度与总压分布均与试验值吻合较好。

thumbnail 图4

数值计算与Goebel和Hall试验结果对比

1.3 参数处理

图 5为混合层示意图。在某流向位置处, 平均速度分别等于U1-0.1ΔUU2+0.1ΔU对应的y值作为混合层的超声速侧和亚声速侧边界点y1y2, 两者间的差值作为混合层的厚度b

无量纲坐标

式中:y为实际坐标; y0为流向位置x处混合层中心位置坐标((y1+y2)/2);=0表示混合层的中心位置;-0.5≤ ≤0.5对应混合层位置, < -0.5为超声速气流部分, >0.5为亚声速气流部分。

无量纲速度Ur和无量纲总压prt

式中, pt1pt2分别为超声速气流和亚声速气流进口平均总压(kPa)。混合层内, 0 < Ur < 1, 0 < prt < 1。亚声速气流无量纲速度Urprt接近0, 超声速气流无量纲速度Urprt接近1。

thumbnail 图5

混合层参数定义示意图

2 结果与分析

2.1 速度分布特征

图 6为混合层x=10~100 mm和x=140~210 mm区域, 时间间隔为1/15 s的2个连续时刻瞬态速度分布云图(流向:从左向右)。上层y>21 mm区域为亚声速气流通道, 下层为超声速气流通道。将2股流体隔开的交界面变化情况反映出2股流体间的混合掺混过程。在x=10~40 mm区域, 交界面近似呈现直线段; 随后, 在x=50~100 mm区间, 交界面出现褶皱和变形, 在混合层的下游区域更为明显, 交界面开始破碎, 并且在亚声速气流通道内出现单独的高速流体微团。对于2个区域, 相邻2个时刻交界面的整体特征相似, 但是形状明显不同, 反映混合层具有瞬态变化特征。

图 7为试验时均速度及数值计算速度分布云图。两者呈现出相似的分布特征, 混合层位于亚声速气流和超声速气流之间, 随着流向发展, 逐渐向两侧流体中扩展, 厚度不断增加; 在x=10~100 mm区域, 试验的流线表明亚声速气流不断流入混合层内, 说明混合层及超声速气流的具有较强卷吸作用。

图 8x=80, 140, 170和190 mm流向位置处, 试验时均速度与数值计算速度对比。对于亚声速气流, 数值计算速度与试验测量速度匹配较好; 在混合层内, 靠近亚声速气流部分, 两者速度吻合也很好; 对于超声速气流, 两者存在差异, 但是随着流向的发展差值逐渐减小。

thumbnail 图6

PIV测量瞬态速度分布云图

thumbnail 图7

PIV时均与数值计算速度分布云图对比

thumbnail 图8

不同流向位置,试验与数值计算速度分布对比

2.2 压力分布特征

图 9为试验和数值计算的超声速侧壁面静压沿流向分布。静压沿流向呈现“波纹”形分布; 波幅随着流动方向不断减小, 波动特征减弱。

图 10为数值计算总压与静压分布云图。亚声速气流总、静压保持不变; 随着流向发展, 混合层内的总压不断增加, 超声速气流的总压及高压区域逐渐减小。超声速气流的静压出现离散的相互交替的低、高静压区, 区域形状由初始的三角形逐渐变为半圆形。

研究中超声速气流与亚声速气流静压比为1.1, 超声速气流轻度欠膨胀, 在隔板下游将经历一系列膨胀加速和压缩减速过程。膨胀时, 静压降低, 压缩时, 静压升高。因此, 超声速壁面静压呈现波纹分布, 静压分布云图中出现离散交替的低、高静压区。

thumbnail 图9

试验与数值计算超声侧壁面静压分布

thumbnail 图10

数值计算总压与静压分布分布

2.3 相似性

混合层发展中, 当其进入充分发展状态后, 各流向位置的发展具有相同的特征。当对各流向位置采用无量纲量统一标度后, 各流动参数分布曲线相互重合, 即流动存在相似性。

图 11为不同流向位置处, 数值计算无量纲速度Ur、无量纲总压prt与无量纲坐标分布规律。各流向位置处无量纲速度和无量纲总压分布重合。Ur=0.5和prt=0.5分别对应=-0.11和=-0.45, 偏向超声速气流侧。其中速度分布在-0.11≤ ≤0.5区间, 误差函数与各流向位置处的速度吻合较好; 而在-0.5≤ ≤-0.11区间, 速度变化率较大, 高斯分布函数较为接近。总压分布在prt < 0.5时, 近似呈现双曲正切函数分布, 而在prt>0.5时, 仍然呈现出高斯函数分布。试验测量的无量纲速度与拟合曲线的对比如图 12所示。各流向位置处速度分布基本相重合且分布型式与拟合的误差函数及高斯函数分布型相近。整体而言, 高斯函数可以表征混合层内速度与总压分布。

thumbnail 图11

不同流向位置处,数值计算Urprt分布规律

thumbnail 图12

PIV测量无量纲速度与拟合曲线对比

2.4 混合层厚度发展

图 13为试验与数值计算的混合层边界线及混合层厚度发展。随着流向发展, 混合层边界向两侧扩展, 混合层厚度不断增加。x=10~100 mm区间, 试验测量的亚声速侧边界数值较大, 可能原因为该位置混合层初始形成, 速度差较大, 混合较为强烈。对试验数据进行最小二乘法线性拟合, 混合层增长速率为0.135 mm/mm。

thumbnail 图13

混合层边界线及厚度随流向发展

3 结论

采用粒子图像测速技术PIV对亚-超声速混合层流场进行测量, 同时, 采用基于雷诺平均方程RANS的标准k-ω湍流模型开展稳态数值模拟, 考虑亚超声速混合层压缩性影响, 获得亚-超混合层基本流场特征, 包括速度分布、压力分布和混合层厚度发展。超声速气流马赫数Ma1为1.32, 亚声速气流马赫数Ma2为0.11, 混合层内静压约为70 kPa。研究得到的主要结论如下:

1) 亚-超声速混合层具有非稳态特征, 2股流体分界面从近似线性分布, 随后发生褶皱变形, 最后发生破碎并出现离散的流体微团;

2) 混合层的时均特征为沿流向近似呈现线性增长, 增长率为0.135 mm/mm。

3) 混合层内速度和总压分布具有相似性特征, 曲线分布呈现高斯函数分布。

References

  1. Abramowich G N. The Theory of Turbulent Jets[M]. Cambridge, MA: MIT Press, 1963 [Google Scholar]
  2. Sabin C M. An Analytical and Experimental Study of the Plane, Incompressible, Turbulent Free-Shear Layer with Arbitrary Velocity Ratio and Pressure Gradient[J]. Journal of Fluids Engineering, 1965, 87 (2): 421– 428 [Google Scholar]
  3. Dimotakis P E. Two-Dimensional Shear-Layer Entrainment[J]. AIAA Journal, 1986, 24 (11): 1791– 1796 [Article] [Google Scholar]
  4. Zhao Yuxin. Experimental Investigation of Spatiotemporal Structures of Supersonic Mixing Layer[D]. Changsha: National University of Defense Technology, 2008(in Chinese) [Google Scholar]
  5. Konrad J H. An Experimental Investigation of Mixing in Two-dimensional Turbulent Shear Flows with Applications to Diffusion-limited Chemical Reactions[D]. Pasadena, California: California Institute of Technology, 1977 [Google Scholar]
  6. Bogdanoff D W. Compressibility Effects in Turbulent Shear Layers[J]. AIAA Journal, 1983, 21 (6): 926– 927 [Article] [Google Scholar]
  7. Papamoschou D, Roshko A. The Compressible Turbulent Shear Layer an Experimental Study[J]. Journal of Fluid Mechanics, 1988, 197 (1): 453 [Article] [Google Scholar]
  8. Dimotakis P E. On the Convection Velocity of Turbulent Structures in Supersonic Shear Layer[C]//AIAA 22nd Fluid Dynamics, Plasma Dynamics and Lasers Conference, 1991 [Google Scholar]
  9. Brown G L, Roshko A. On Density Effects and Large Structure in Turbulent Mixing Layers[J]. Journal of Fluid Mechanics, 1974, 64: 775– 816 [Article] [Google Scholar]
  10. Hall J L. An Experimental Investigation of Structure, Mixing and Combustion in Compressible Turbulent Shear Layers[D]. Pasadena, California: California Institute of Technology, 1991 [Google Scholar]
  11. Zhang C X, Liu Y, Fu B S, et al. Direct Numerical Simulation of Subsonic-Supersonic Mixing Layer[J]. Acta Astronautica, 2018, 153: 50– 59 [Article] [Google Scholar]
  12. Ping Jun. Theory and Application of Jet[M]. Beijing: Astronautics Press, 1995 (in Chinese) [Google Scholar]
  13. Smits A J, Dussauge J P. Turbulent Shear Layers in Supersonic Flow[M]. Liu Hong, Zhang Bin, Translator. Beijing: Aviation industry press, 2016(in Chinese) [Google Scholar]
  14. Zhang Fuxiang. Dynamics of Rocket Jet[M]. Harbin: Harbin Engineering University Press, 2004 (in Chinese) [Google Scholar]
  15. Xu Jinglei. The Development of the PIV Experimental study of the Super/Hypersonic Flow Field[J]. Advance in Mechanics, 2012, 42 (1): 85– 94 [Article](in Chinese) [Google Scholar]
  16. Rong Zhen. Research on the PIV Technique and Application in Supersonic Flow[D]. Shanghai: Shanghai Jiaotong University, 2012(in Chinese) [Google Scholar]
  17. Goyne C P, Mcdaniel J C, Krauss R H, et al. Velocity Measurement in a Dual-Mode Supersonic Combustor Using Particle Image Velocimetry[R]. AIAA-2001-1761 [Google Scholar]
  18. Humphreys W M, Bartram S M, Blackshire J L. A Survey of Particle Image Velocimetry Applications in Langley Aerospace Facilities[R]. AIAA-1993-0411 [Google Scholar]
  19. Haertig J, Havermann M, Rey C, et al. Particle Image Velocimetry in Mach 3.5 and 4.5 Shock-Tunnel Flows[J]. AIAA Journal, 2002, 40 (6): 1056– 1060 [Article] [Google Scholar]
  20. He Lin. Experimental Investigation of Supersonic Boundary Layer and Shock Wave/Boundary Layer Interaction[D]. Changsha: National University of Defense Technology, 2011(in Chinese) [Google Scholar]
  21. Chen Maozhang. Fundamentals of Viscous Fluid Dynamics[M]. Beijing: Higher Education Press, 2002 (in Chinese) [Google Scholar]
  22. Goebel S G, Dutton J C. Experimental study of Compressible Turbulent Mixing Layers[J]. AIAA Journal, 1991, 29 (4): 538– 546 [Article] [Google Scholar]

All Tables

表1

试验参数

All Figures

thumbnail 图1

试验系统原理

In the text
thumbnail 图2

PIV系统及试验段图片

In the text
thumbnail 图3

计算域与网格

In the text
thumbnail 图4

数值计算与Goebel和Hall试验结果对比

In the text
thumbnail 图5

混合层参数定义示意图

In the text
thumbnail 图6

PIV测量瞬态速度分布云图

In the text
thumbnail 图7

PIV时均与数值计算速度分布云图对比

In the text
thumbnail 图8

不同流向位置,试验与数值计算速度分布对比

In the text
thumbnail 图9

试验与数值计算超声侧壁面静压分布

In the text
thumbnail 图10

数值计算总压与静压分布分布

In the text
thumbnail 图11

不同流向位置处,数值计算Urprt分布规律

In the text
thumbnail 图12

PIV测量无量纲速度与拟合曲线对比

In the text
thumbnail 图13

混合层边界线及厚度随流向发展

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.