水流冲击超声速气体射流实验研究

通过水洞实验对有水流速度影响的水下超声速气体射流进行实验研究,通过高速摄像系统记录射流形态演变过程,采用动态测力系统测量射流演变过程中射流周围环境压力的脉动特征.对剪切层涡旋结构进行分析,得到水流冲击射流的剪切涡流动形态演化和压力脉动特征.实验结果表明,射流主体形态的非定常运动依赖于水流速度,无水流速度时,射流主体受到重浮力影响向上弯曲较大,并且可以捕捉到射流的振荡诱导喷管口平面处主频为200 Hz的压力脉动,当存在水流速度时,射流主体向下游发展过程中的偏斜程度较小,射流与水流相互作用形成剪切涡,在水流作用下射流主体向下游发展过程中卷入射流剪切层,与射流主体掺混形成较大尺度的涡结构,喷管口平面处主频消失.     

高速受限混合层结构和效应实验研究

受限混合层的流动主要是喷流与自由来流相互剪切形成的混合层受到壁面的限制而形成的一种流动.文章采用后向台阶平板模型研究了高速高压比条件下的受限混合层的典型流场结构以及冷却效率.实验自由来流Mach数为5, 喷流的Mach数为1.28, 喷流总压为0.2~0.7 MPa, 通过调整冷喷气流的总压, 基于纹影流动显.形成喷口附近波系的欠膨胀流动现象的深刻认识, 提取波系特征与流动参数之间的规律.基于流动显示及实验测量结果, 通过分析流场中大尺度结构的空间演化规律, 揭示流动参数对于冷却效率的影响规律及物理内涵.采用快响应压敏漆(FRPSP)技术在高超声速风洞开展热流分布和冷却效率研究, 获得了平板对受限混合层冷却效率的影响.      

超声速三维空腔流气动噪声被动控制

抑制超声速武器舱空腔流噪声是航空领域中一项重要课题。大量研究表明在空腔前缘采用主/被动控制技术可以在一定程度上抑制腔内噪声水平。利用大涡模拟(large eddy simulation, LES)技术计算分析了Mach 1.4开式矩形方腔及波形、弧形两种前后壁几何修形后空腔的流动及噪声, 探索超声速来流条件下几何修形被动控制技术对开式方腔流噪声的抑制能力。计算结果表明波形和弧形空腔对腔内噪声均具有一定的抑制作用, 且波形空腔噪声控制效果更优。分析认为空腔几何修形能够改变空腔上方剪切层及腔内大尺度涡结构的发展演化, 进而实现对腔内噪声的控制。此外, 还应用LES方法计算分析了增厚的来流边界层条件下超声速方腔流, 发现来流边界层增厚可显著降低腔内噪声水平。      

高精度非定常激波装配法

为正确模拟高超声速绕流中,来流小扰动与弓形激波之间的干扰对流动特征的影响,将弓形激波作为动边界,利用非定常特征关系处理激波处的边界条件.应用五阶精度迎风紧致格式和六阶精度的对称格式与三阶精度的R-K方法相结合,建立高精度非定常激波装配方法.采用该方法数值模拟钝锥高超声速定常流场和二维抛物外形高超声速边界层流动的感受性问题,数值模拟来流小扰动与弓形激波干扰激波后非定常扰动流场,研究扰动波进入边界层产生边界层不稳定波的特征.     

基于NPLS的超声速后台阶流场精细结构及其非定常特性

在马赫数3.4的低噪声超声速风洞中,分别采用纹影技术及基于纳米示踪的平面激光散射技术(NPLS),对台阶高度为3 mm的层流后台阶绕流结构进行了流动显示实验,旨在研究后台阶流场再附区的平均流动规律及再附边界层的瞬态流场精细结构,并通过比较不同时刻的NPLS数据,发现后台阶再附边界层涡结构的非定常演化特性;并且与马赫数4.2条件下的纹影结果和前期研究成果进行了比较.结果表明,两种流动显示方法均可以捕捉激波膨胀波的位置,但NPLS技术在揭示流场某一截面内的瞬态精细结构方面更具优势,其时间分辨率可达6 ns,空间分辨率微米量级;该流场条件下的回流剪切层的倾角为-3.1?,再附边界层的平均增长斜率为0.07519;再附边界层中发卡涡脱落的特征时间尺度为大约10μs.相同膨胀比条件下,再附随着马赫数增大而推迟;膨胀比越大,再附越提前,并且流动折转角变大.     

耗散紧致格式的频谱特性研究与应用

采用Fourier分析方法,给出线性耗散紧致格式和基于m级Runge-Kutta时间积分方法的全离散格式的频谱特性,并应用五阶耗散紧致格式模拟典型高频波传播和超声速平面Couette流动的特征值问题及其稳定性边值问题,展示耗散紧致差分格式良好的频谱特性.     

带喷流超声速后台阶流场精细结构及其运动特性研究

采用基于纳米示踪的平面激光散射技术(NPLS)对带超声速喷流的后台阶流动精细结构进行了研究. 来流马赫数为3.4, 喷流实测马赫数为2.45, 而名义马赫数为2.5. 结果清晰地揭示了激波、剪切层、混合层、Kelvin-Helmholtz涡、羊角涡及湍流大尺度结构等大量典型流场结构. 基于大量流场精细结构图像, 对典型位置处的结构进行了空间两点相关性分析, 在喷流混合层前端涡结构小于湍流充分发展的尾端, 结构角相对也小. 喷流工作时, 模型台阶下游表面由一薄层气膜覆盖. 获得了模型流向和不同高度展向平面内的流场结构, 对照纹影试验结果, 分析了流动特点及时间演化规律. 采用微型压力扫描系统测试了模型表面的压力系数分布, 靠近喷流下游处压力系数区域0.0146. 针对NPLS图像做了流动的分形维数的分析, 发现在流动初始阶段分形维数接近于1, 越靠下游分形维数越高.     

基于三角波瓣混合器的超声速流场精细结构和掺混特性

在超声速吸气式混合层风洞中,采用基于纳米粒子的平面激光散射(NPLS)技术对平板混合层和三角波瓣混合器诱导的混合层流场精细结构进行了对比实验研究.上下两层来流的实测马赫数分别为1.98和2.84,对流马赫数为0.2.NPLS图像清晰地展示了Kelvin-Helmholtz涡、流向涡、波系结构以及大尺度涡结构的配对合并过程.通过对比分析时间相关的NPLS流场图像,发现了大尺度拟序结构随时间发展演化的非定常特性.基于流动显示结果,采用分形维数和间歇因子指标对流场结构和混合特性进行了定量分析.实验研究表明,三角波瓣混合器诱导的流向涡结构显著提高了上下两层来流的掺混效率,其流动远场的分形维数突破了平板混合层中完全湍流区的分形维数值,达到了1.88,流场结构表现出明显的破碎性,有利于流动在标量层面的扩散和掺混.流动间歇性分析表明,流向涡与展向涡的相互剪切作用主导着混合层的掺混特性,同时由于流向涡的卷吸作用,三角波瓣混合器诱导的混合层混合区域更大,更多的流质被卷入混合区完成混合.     

超声速层流/湍流压缩拐角流动结构的实验研究

在Ma=3.0的超声速风洞中, 分别对上游边界层为超声速层流和湍流, 压缩角度为25°和28°的压缩拐角流动进行了实验研究. 采用纳米粒子示踪平面激光散射(NPLS)技术获得了流场整体和局部区域的精细结构, 边界层、剪切层、分离激波、回流区和再附激波等典型结构清晰可见, 测量了超声速层流压缩拐角壁面的压力系数. 从时间平均的流场结构中测量出分离激波、再附激波的角度和再附后重新发展的边界层的增长情况, 通过分析时间相关的流场NPLS图像, 可以发现流场结构随时间的演化特性. 实验结果表明: 在25°的压缩角度下, 超声速层流压缩拐角流动发生了典型的分离, 边界层迅速增长失稳转捩, 并引起一道诱导激波, 流场中出现了K-H涡、剪切层和微弱压缩波结构, 而超声速湍流压缩拐角流动没有出现分离, 湍流边界层始终表现为附着状态; 在28° 的压缩角度下, 超声速层流压缩拐角流动进一步分离, 回流区范围明显扩大, 诱导激波、分离激波向上游移动, 再附激波向下游移动, 分离区流动结构复杂, 相比之下, 超声速湍流压缩拐角流动的回流区范围明显较小, 边界层增长缓慢, 流场中没有出现诱导激波、K-H涡和压缩波, 流动分离区域的结构也相对简单, 但分离激波的强度则明显更强. 关键词： 压缩拐角 层流 湍流 流动结构     

求解Navier-Stokes方程组的组合紧致迎风格式

给出一种新的至少有四阶精度的组合紧致迎风(CCU)格式,该格式有较高的逼近解率,利用该组合迎风格式,提出一种新的适合于在交错网格系统下求解Navier-Stokes方程组的高精度紧致差分投影算法.用组合紧致迎风格式离散对流项,粘性项、压力梯度项以及压力Poisson方程均采用四阶对称型紧致差分格式逼近,算法的整体精度不低于四阶.通过对Taylor涡列、对流占优扩散问题和双周期双剪切层流动问题的计算表明,该算法适合于对复杂流体流动问题的数值模拟.     

Numerical simulations of shock-wave interaction with a boundary layer in the plane supersonic flows with jet injection

A supersonic air flow in a plane channel with a transverse turbulent jet of hydrogen injected through a slot on the bottom wall is simulated. The algorithm for solving the Favre-averaged Navier-Stokes equations for the flow of a perfect multispecies gas on the basis of the WENO scheme is proposed. The main attention is paid to the interaction of the shock-wave structure with the boundary layers on the upper and lower duct walls under the conditions of an internal turbulent flow. Namely, a detailed study of the structure of the flow is done, and separation and mixing depending on the jet slot width are investigated. It is found that in addition to well-known shock-wave structures produced by the interaction of the free stream with the transverse jet and the bow shock interaction with the boundary layers near the walls, an additional system of shock waves and the flow separation arise on the bottom wall downstream at some distance from the jet. The comparison with the experimental data is performed.     

Computation of the normal injection of a hydrogen jet into a supersonic air flow

The interaction of a planar supersonic air flow with the hydrogen jet injected perpendicularly across the slot from the duct walls is studied numerically. An algorithm is constructed for solving the Favre-averaged Navier — Stokes equations for the flow of a thermally perfect multicomponent gas based on the ENO scheme. The influence of the jet Mach number and the ratio of the jet and flow pressures on the shockwave structure of the flow and the jet penetration depth is shown.     

激波/边界层干扰对等离子体合成射流的响应特性

利用高速纹影系统和数值模拟方法研究了激波/边界层干扰对逆流喷射的等离子体合成射流的响应特性,并揭示了流动控制机理.实验在来流马赫数Ma=3.1的风洞中进行,测试模型采用钝头体和压缩斜坡的组合模型,等离子体合成射流激励器安装在钝头体头部.纹影系统捕捉了放电频率为f=1 kHz和f=3 kHz的激励对附体激波形态和分离激波运动的控制效果.等离子体合成射流使压缩斜坡激波/边界层干扰区域的起始点向下游移动,分离泡尺寸减小,附体激波强度减弱,发生弯曲,再附点移向上游,与此同时分离激波向附体激波逼近.与f=3 kHz激励相比,f=1 kHz激励的射流流量更大,对激波/边界层干扰的影响范围更广、控制效果更好.通过数值模拟,揭示了射流与来流相互作用对下游流场的作用机理:射流与来流相互作用诱导出大尺度旋涡,大尺度旋涡耗散发展增强了近壁面流场的湍流度;压缩斜坡上游近壁面的流场性质发生变化,进而导致了压缩斜坡激波/边界层干扰区域流动的变化.     

Influence of the opening angle of a conical supersonic nozzle on the structure of initial interval of non-isobaric jet

The ideal gas exhaustion from an infinite volume into a gas at rest through a supersonic conical Laval nozzle is considered. The problem was solved numerically by steadying in time in a unified formulation for the regions inside the nozzle and in the ambient environment. In such a statement, the nozzle outlet section is no internal boundary of the region under consideration, and there is no need of specifying the boundary conditions here. Local subsonic zones arising in the flow lie inside the region under consideration, which eliminates the possibility of using a marching technique along one of the coordinates. The numerical solution is constructed by a unified algorithm for the entire flow region, which gives a possibility of obtaining a higher accuracy. The computations are carried out in the jet initial interval, where, according to monograph [1], the wave phenomena predominate over the viscous effects. The exhaustion process is described by the system of gas dynamics equations. Their solution is constructed with the aid of a finite difference Harten’s TVD (Total Variation Diminishing) scheme [2], which has the second approximation order in space. The second approximation order in time is achieved with the aid of a five-stage Runge-Kutta method. The solution algorithm has been parallelized in space and implemented on the multi-processor computer systems of the ITAM SB RAS and the MVS-128 of the Siberian Supercomputer Center of SB RAS. The influence of the semi-apex angle of the nozzle supersonic part and the pressure jump between the nozzle outlet section and the ambient environment on the flow in the initial interval of a non-isobaric jet is investigated in the work. A comparison with experimental data is presented. The computations are carried out for the semi-apex angles of the nozzle supersonic part from 0 (parallel flow) to 20 degrees. For all considered nozzles, the Mach number in the nozzle outlet section, which was computed from the one-dimensional theory, equaled three. Computations showed that in the case of flow acceleration in a conical supersonic nozzle, its geometry is one of the main factors determining the formation of the jet initial interval in ambient environment.     

合成冷/热射流控制超声速边界层流动稳定性

为了探究超声速边界层流动稳定性及其转捩控制机理,提出基于合成冷/热射流的边界层速度-温度耦合控制方法,并通过数值模拟研究了Ma=4.5超声速平板边界层不稳定波的传播,采用线性稳定性理论中的时间模式分析了壁面吹吸、射流温度、扰动频率、扰动振幅等对不稳定波控制效果的影响.结果表明:无射流控制时,边界层内同时存在不稳定的第一模态扰动波和第二模态扰动波,且二维波形式的第二模态占主导地位;壁面吹吸作用下,仅出现更加不稳定的第二模态,第一模态被抑制;速度-温度耦合控制下,射流温度对扰动模态的不稳定区域大小及扰动增长率影响显著,射流温度与来流温度不同时,温度的脉动使得流动转捩为湍流的速度加快,边界层速度型更加饱满,抗干扰能力增强,流动稳定性提高;高频的吹吸扰动对流场的控制效果优于低频扰动,扰动频率超过400 Hz时,第二模态扰动波时间增长率降低,扰动分量对边界层速度剖面和温度剖面的修正加快,第二模态更加稳定;扰动振幅减小为主流速度的1%时,仅出现时间增长率较小的第二模态,控制效果较好,进一步减小时,第一模态重新出现,并且波数范围与第二模态先重合后分离,对应的时间增长率先增加后减小.研究结果为边界层转捩控制技术提供了新的思路.     

等离子体激励器的静特性与高速特性仿真研究

A parametric study of high-frequency plasma jet actuator was carried out, using the experiment- tally measured energy distribution law of arc discharge as an ideal heat source. The influence of the exit angle of the actuator on the flow field was explored. The jet flow field characteristics of the spark discharge actuator under supersonic flow (Ma0=2.0) were investigated. The results show that the energy density of heat flux increases and the jet front and forward shock wave moves faster with the decrease of discharge region, and the smaller the exit angle of the jet is, the stronger the momentum injection ability of the actuator along the flow direction is. The rule still applies under high-speed air flow conditions. Compared to the static condition, the momentum injection capability of the jet is stronger and the influence domain is larger under supersonic flow conditions.     

Interaction of an Optical Discharge with a Shock Wave

The results of an experimental study of the impact of the focused pulsed-periodic radiation from a CO₂ laser on a gas-dynamic structure in a supersonic jet are presented. The radiation of the CO₂ laser is propagated across the stream and focused by a lens on the axis of the supersonic jet. To register the flow structure, a shadow scheme with a slit and a flat knife located along the flow is used. The image is fixed by a speed camera with an exposure time of 1.5 μs and a frame rate of 1000 1/s. In the flow, the plasma initiated by the pulsedperiodic laser is visualized in order to identify and determine the period of plasma development, as well as the motion of the initial front of the shock wave. It is shown that at the transverse input of laser radiation into the stream the periodic structure of the thermal trace is created with the formation of an unsteady shock wave from the energy release zone. At small repetition rates of laser radiation pulses, the thermal spot interacts with the flow in the pulsed mode. It is shown that elliptic nonstationary shock waves are formed only at low subsonic flow velocities and in a stationary atmosphere. The process of nonstationary ignition by an optical discharge of a methane–air mixture during a subsonic outflow into a motionless atmosphere is shown experimentally. The results of optical visualization indicate burning in the trace behind the optical discharge region.     

Control of spray spot in cold spray technology. Part 1. Gas dynamic aspects

This paper presents a study of supersonic jets formed by approaches that are new for cold spray technique: the main flow is swirled, the nozzles with permeable profiles and with exit slots on the supersonic section are engineered. The flow swirling achieved in the nozzle prechamber retains downstream to substrate surface. The system of vortices created within the permeable nozzles changes the shock wave features of the overexpanded jet and the geometry of the bow shock wave ahead of the substrate surface. These new features of flow may affect particle motion and particlesubstrate interaction under the conditions of cold spray process; this offers tools for obtaining the necessary shape of a spray spot.     

进口条件对超燃冲压发动机氢气燃烧流场特性的影响分析

对不同进口条件下的超燃冲压发动机燃烧室内氢气喷流超声速燃烧流动特性进行了数值模拟与分析.宽范围超燃冲压发动机是吸气式高超声速飞行器推进系统设计中的热点问题之一,受实验设备硬件条件及实验技术限制,数值模拟技术仍然是超燃冲压发动机燃烧室内燃气燃烧特性及流场特性的主要研究手段.采用基于混合网格技术的多组元N-S方程有限体积方法求解器,在不同进口Mach数及压强条件下,对带楔板/凹腔结构的燃烧室模型氢气喷流燃烧流场进行了数值模拟,对比分析了氢气喷流穿透深度、喷口前后回流区结构、掺混效率及燃烧效率等流场结构与典型流场参数的变化特性及影响规律.研究成果可为宽范围超燃冲压发动机喷流燃烧流动特性分析提供参考.      

弹体尾部轴对称干扰流场数值模拟

绕弹体的超声速气流与发动机喷流相互作用,在尾部形成复杂的干扰流场。本文用有限差分法求解全N-S方程,对这一复杂流场进行了数值模拟,得到了实验观察到的各种流场结构及其随喷口压力比的变化规律。外流M_∞=1.94,R_e∞=2.2×10⁵,喷流M_j=3.0,喷口压力比p_j/p_∞分别为1.03,0.527,0.15三种。差分算法为一种改进的Beam-Warming格式。计算底部压力和激波在喷流中心的反射位置与实验数据进行了比较,吻合较好。