空心活塞运动产生的激波

普通活塞在管道中运动, 不论活塞速度高低都能驱动产生激波.空心活塞在管道中运动与普通活塞不同, 当其运动Mach数低于壅塞Mach数, 活塞前不会形成运动激波; 只有活塞运动Mach数超出壅塞Mach数, 才能驱动产生激波.壅塞Mach数由空心活塞中心孔直径、长度和内壁面摩擦系数以及外直径横截面积决定.火车通过隧道类似空心活塞在管道中运动.将火车车箱与隧道间的空隙换算成有效水力直径后, 可利用空心活塞在管道中运动计算壅塞Mach数的方法, 求出火车通过隧道出现激波的条件.文章提出一种新颖的间接测量列车平均摩擦系数新方案.      

不同界面组分分布对Richtmyer-Meshkov不稳定性的影响

基于二维非定常Euler方程,对平面激波与不同界面组分分布下氦气气柱作用过程所引起的Richtmyer-Meshkov不稳定性现象进行了数值模拟,探讨了激波冲击轻质气柱后气柱界面形态的演变及流场波系结构,定量分析了气柱特征尺度(气柱长度、高度和中轴宽度)和气柱体积压缩率随时间变化.此外,结合流场压强、速度、环量和气体混合率,多角度分析了激波驱动界面气体混合的流动机制,获得了不同界面组分分布对界面不稳定性的影响.结果表明,随着气柱界面从完全扩散界面向间断界面的过渡,界面两侧的声反射系数随之增大,使入射激波与气柱界面的作用由常规透射转变为非常规透射,反射激波逐渐加强,透射激波逐渐减弱,使得Richtmyer-Meshkov不稳定性随之增强;同时,界面两侧阿特伍德数的增大,加强了Rayleigh-Taylor不稳定性和Kelvin-Helmholtz不稳定性的发展.此外,界面不稳定性的加强使得流场环量增大,导致气体混合率的增长速率随之升高.     

斜激波与轴对称边界层干扰的数值模拟

使用两方程Menter-SST模型, 对来流Mach数为3时的斜激波与轴对称边界层的相互干扰现象进行了数值模拟与定性、定量分析。研究了斜楔激波发生器楔角和来流单位Reynolds数变化对干扰区流动的影响, 总结了参数变化引起的流动分离变化规律; 此外, 还计算了与三维计算中心对称面上的入射激波等效的二维情形, 并将三维结果与二维情形进行对比, 对比结果显示中心对称面上的壁面压力系数、分离涡尺寸、涡量分布等与相应的二维情形存在明显差异。      

Laval喷管内激波位置的计算及制冷性能分析

明确激波在Laval喷管内的发生位置,能够为设计及优化Laval喷管线型提供理论依据,从而提高Laval喷管的制冷性能和整个天然气超音速分离器的分离效率。文中对Laval喷管进行了结构设计,对不同背压条件下Laval喷管内激波位置进行了理论分析与计算,并利用FLUENT软件进行了数值模拟,通过研究不同背压条件下Laval喷管内气体马赫数、压力和温度分布,对喷管制冷性能进行了对比分析。结果表明:保持Laval喷管入口压力不变,随出口背压增大,激波位置逐渐从喷管出口向喷管入口方向移动,气体受到激波的影响,在喷管所能达到的最大马赫数不断降低,所能产生的最低温度不断上升,喷管的制冷性能越差;理论计算与数值模拟结果基本一致,数值模拟验证了理论计算的正确性。     

水蒸气跨声速流动中非平衡相变与激波的耦合作用

引入水蒸气非平衡相变的动力学模型和水蒸气真实物性模型,建立了水蒸气跨声速非平衡流动的守恒型数值计算模型,采用Roe-FDS计算格式数值捕捉了水蒸气超声速流动中的非平衡相变与激波效应,在此基础上开展了气动激波与凝结激波的耦合计算,分析了气动激波波锋面与凝结激波波锋面相遇时,气动激波与非平衡相变之间的相互作用规律。研究显示,随着背压的不断升高,气动压缩激波不断向喷管喉口位置推移。当气动激波发生在凝结激波的下游位置时,气动激波的耗散效应使得喷管内的液相质量分数逐渐减小而不会对上游的非平衡相变和凝结激波产生影响。当气动激波随背压继续上行与气动激波交汇时,气动激波强烈的耗散效应使得凝结激波特征迅速减弱,非平衡相变逐渐退化到喷管边界区域,而气动激波由于受到凝结激波的强烈干扰,激波强度迅速减弱,显现出明显的斜弱激波特征.当气动激波上行至喷管喉口附近时,X型凝结激波逐渐消退,非平衡相变在喷管主流区消失。     

激波对超音速平板发汗冷却的影响

本文分析了激波影响下,不同激波强度,冷却气体和注入率对发汗冷却效果的影响。马赫数3的超音速主流遇到流道内的楔形激波发生器产生斜激波入射到多孔平板表面。楔块楔角φ分别为0°、4°、8°、12°模拟不同的入射激波强度,冷却介质分别为空气、甲烷与氢气。计算结果表明,激波使多孔区域出口表面静压上升,冷却流体流出受阻而破坏冷却效果,且冷却效率随激波强度增强而下降,随冷却气体分子量增大,冷却效率下降幅度减小;但在激波强度较强时,激波在多孔表面形成逆压梯度,迫使冷却流体流向入射点上游区域而使该区域冷却效果得到恢复。提高冷却剂注入率可以减弱激波对冷却效率的破坏作用,但使壁面温度不均匀性增加。     

激波矢量喷管二次流喷口形态影响研究

基于CFD数值模拟技术,考虑变比热容及温度对传热系数、黏性系数的影响,研究了不同二次流喷口几何构型(矩形及圆孔)在不同落压比、二次流压比等工况下,对激波矢量喷管内的三维流动特性及推力矢量性能的影响。研究表明:随着二次流压比增加,诱导激波角度增加,当诱导激波与上壁面相交时,推力矢量性能恶化;对于矩形喷口无量纲展向长度l1.0时,喷口前分离涡演变为马蹄涡,并在喷口下游诱导尾涡,SPR0.6时随喷口无量纲展向长度增大,推力矢量角度增加;对与多孔喷射模型,在压力缓解机理下,不同圆喷孔数量时,下壁面压力分布及壁面极限流线分布存在差异,然而圆喷孔数量对推力矢量角及推力系数的影响仍待进一步探明。     

弹头激波诱导燃烧特性的数值模拟

为了研究弹头激波诱导燃烧，基于有限体积的考虑化学反应的Navier-Stokes（N-S）方程，对预混氢气-空气化学恰当比时的燃烧流场进行了数值模拟.时间项基于2阶隐式LU-SGS格式，对流项基于Steger-Warming进行离散，化学反应源项采用对角化隐式处理.首先，研究了网格对燃烧爆轰流场结构的影响，并利用Lehr实验结果验证了计算方法的可靠性；其次，研究了弹头的飞行Mach数（Ma=4.18，5.11，6.46）、弹头直径（D=5，10，15 mm）对燃烧流场稳定性的影响.研究表明:计算网格对氢气-空气爆轰流场结构影响很大；弹头直径一定时，氢气-空气燃烧流场稳定性随着飞行Mach数的增大而增强；弹头飞行Mach数一定时，氢气-空气燃烧流场稳定性随着弹头直径减小而增强.      

激波和剪切层相互作用下的超音速射流

欠膨胀超音速射流处于螺旋模式下的中度欠膨胀时,其入射剪切层的激波具有很高强度,激波和剪切层发生了强烈的相互作用,远场辐射的拢动波出现了大间隔、交错的上下行类似螺旋锥面波形图像,该扰动波具有很强的向上游传播的指向性,导致上游噪声高于垂直喷嘴方向的声压级.而在相对压比较低的低度欠膨胀情况下,或高压比下的高度欠膨胀的情形,入射剪切层激波强度相对较弱,远场辐射没有大间隔、交错的上下行远场辐射 关键词： 超音速射流 啸叫 扰动波 激波 剪切层     

入射激波诱导环氧丙烷爆炸特征研究

用自行设计激波管点火测试技术,实验研究了温度范围760-1380K间入射激波诱导下环氧丙烷的点火机理。利用激波管压力传感器测定了H*(486.1) 和O (470.5nm)随激波诱导强度变化的点火时间特征。实验结果表明：在低马赫数下氢氧自由基出现时间较接近,1.5-2.5马赫间随激波诱导强度增大而线性减小;而马赫大于2.5后,氧自由基的出现时间迅速减小,是由于高活化能的氧自由基的点火时间对强激波较敏感,而诱导强度大于3.5马赫后对两者点火影响区别就下明显了。实验数据将有益于含能材料点火时间的研究。     

水中斜激波的传播

 利用水的简化物态方程讨论了水中斜激波波后流动马赫数、密度比、压强增量及冲击角与来流马赫数的关系,给出了一个无量纲化的激波极线方程;对附体激波转变为脱体激波及波后流动为等声速流的条件进行了讨论。     

Large-eddy simulation of shock-turbulence interaction in supersonic diffuser flows

LES based on explicit filtering is used to study the shock train phenomenon in turbulent supersonic diffuser flows with circular cross-section and isothermal wall with an incoming pipe flow at friction Reynolds number 245 and centerline Mach number 1.7. Alternate regions of compression and expansion are found in the shock train which is followed by a shock-free ‘mixing’ region as observed in experiments and simulations in the literature. Turbulence amplification and local peaks in pressure-dilatation correlation are observed in the vicinity of the shocks. Low-frequency oscillations of the shock train are also observed.     

强梯度复杂流场中的粒子动力学响应试验研究

示踪粒子在(高)超声速流场中的动力学响应是粒子成像测速等粒子示踪测量技术的关键问题之一.现有文献对粒子动力学响应的试验测量往往是通过单个斜激波响应的测量方法.然而,当示踪粒子用于测量高速飞行器发动机内部复杂的激波串流场时,粒子将经历由多道激波导致的速度、压力、黏性等剧烈变化.本文结合目前(高)超声速飞行器的研究热潮,重点关注示踪粒子在应用于发动机内部具有连续激波的复杂流场测量中存在的跟随性评估方面,开展了一系列的相关试验研究.包括测量超声速风洞的喷管出口速度分布以验证测试系统的性能,在马赫4.2和3.0流场中测量了粒子对二维10°和15°单斜劈绕流中的斜激波动力响应,并测量了模拟发动机内部连续梯度的双斜劈粒子斜激波动力响应.结合粒子动力学的理论模型,得到了各状态的粒子弛豫时间、弛豫距离、Stokes数.基于图像方法、统计学规律分析了激波非定常抖动对测量结果的影响,并对测量结果进行了修正.结果显示,相同斜劈角度下,马赫数越高,粒子的弛豫时间、弛豫距离就越大.但是在相同的来流马赫数下,斜劈角度越大,粒子的弛豫时间、弛豫距离反而减小.在强梯度之后由于流场的雷诺数和黏性系数变化剧烈,粒子的跟随性降低了大约5.7%,stokes数增加了约1%.虽然在本文条件下Stokes数仍满足超声速流场对粒子跟随性的要求,但粒子响应的降低无疑是值得关注的,尤其是当其被应用于具有更多连续梯度的复杂流场测量中.     

激光引致激波演化数值研究

在高重复频率激光推进的研究中,激波的合并是发生在激波演化后期的,同时由于脉冲间隔短,脉冲宽度对流场演化的影响也需要详细研究。考虑了激光辐照过程对流场演化的影响,通过数值计算对激波演化特性进行了研究。结果表明,初期波阵面的椭球形状逐渐转化为一个球形,球心与击穿点的距离随时间逐渐减小并最终趋于稳定。基于激波合并的应用,当激波马赫数在1~2之间时,给出了激波波阵面半径随时间的变化规律,以及激波高压区长度和波峰压强随激波波阵面半径变化规律的经验公式。     

On the Global Existence and Stability of a Multi-Dimensional Supersonic Conic Shock Wave

We establish the global existence and stability of a three-dimensional supersonic conic shock wave for a compactly perturbed steady supersonic flow past an infinitely long circular cone with a sharp angle. The flow is described by a 3-D steady potential equation, which is multi-dimensional, quasilinear, and hyperbolic with respect to the supersonic direction. Making use of the geometric properties of the pointed shock surface together with the Rankine–Hugoniot conditions on the conic shock surface and the boundary condition on the surface of the cone, we obtain a global uniform weighted energy estimate for the nonlinear problem by finding an appropriate multiplier and establishing a new Hardy-type inequality on the shock surface. Based on this, we prove that a multi-dimensional conic shock attached to the vertex of the cone exists globally when the Mach number of the incoming supersonic flow is sufficiently large. Moreover, the asymptotic behavior of the 3-D supersonic conic shock solution, which is shown to approach the corresponding background shock solution in the downstream domain for the uniform supersonic constant flow past the sharp cone, is also explicitly given.     

考虑熵层效应的扫掠激波/湍流边界层干扰特性研究

为了探究熵层对扫掠激波/湍流边界层干扰特性的影响规律,采用仿真方法对尖鳍/钝板物理模型进行研究。结果表明:扫掠激波上游的熵层厚度随着平板前缘钝化半径的增大而增加,同时边界层厚度也随着熵层厚度的增加而增加。熵层的引入并不改变扫掠激波/湍流边界层干扰固有的准锥形相似特性,也不会改变拟锥原点（virtual conical origin,VCO）的位置,仅会改变干扰形成的上游影响线和分离线的角度。扫掠激波/湍流边界层干扰形成的锥形主旋涡和角涡的尺度随着熵层厚度的增加而增大。上游熵层的引入增大了下游扫掠激波/湍流边界层干扰区的总压损失,但扫掠激波/湍流边界层干扰自身造成的相对总压损失并不受上游熵层的影响。      

Direct numerical simulation of shock/turbulent boundary layer interaction in a supersonic compression ramp

A direct numerical simulation of the shock/turbulent boundary layer interaction flow in a supersonic 24-degree compression ramp is conducted with the free stream Mach number 2.9.The blow-and-suction disturbance in the upstream wall boundary is used to trigger the transition.Both the mean wall pressure and the velocity profiles agree with those of the experimental data,which validates the simulation.The turbulent kinetic energy budget in the separation region is analyzed.Results show that the turbulent production term increases fast in the separation region,while the turbulent dissipation term reaches its peak in the near-wall region.The turbulent transport term contributes to the balance of the turbulent conduction and turbulent dissipation.Based on the analysis of instantaneous pressure in the downstream region of the mean shock and that in the separation bubble,the authors suggest that the low frequency oscillation of the shock is not caused by the upstream turbulent disturbance,but rather the instability of separation bubble.     

Numerical Study of Air Flow Induced by Shock Impact on an Array of Perforated Plates

This study is focused on the propagation behavior and attenuation characteristics of a planar incident shock wave when propagating through an array of perforated plates. Based on a density-based coupled explicit algorithm, combined with a third-order MUSCL scheme and the Roe averaged flux difference splitting method, the Navier–Stokes equations and the realizable k-ε turbulence model equations describing the air flow are numerically solved. The evolution of the dynamic wave and ring vortex systems is effectively captured and analyzed. The influence of incident shock Mach number, perforated-plate porosity, and plate number on the propagation and attenuation of the shock wave was studied by using pressure- and entropy-based attenuation rates. The results indicate that the reflection, diffraction, transmission, and interference behaviors of the leading shock wave and the superimposed effects due to the trailing secondary shock wave are the main reasons that cause the intensity of the leading shock wave to experience a complex process consisting of attenuation, local enhancement, attenuation, enhancement, and attenuation. The reflected shock interactions with transmitted shock induced ring vortices and jets lead to the deformation and local intensification of the shock wave. The formation of nearly steady jets following the array of perforated plates is attributed to the generation of an oscillation chamber for the inside dynamic wave system between two perforated plates. The vorticity diffusion, merging and splitting of vortex cores dissipate the wave energy. Furthermore, the leading transmitted shock wave attenuates more significantly whereas the reflected shock wave from the first plate of the array attenuates less significantly as the shock Mach number increases. The increase in the porosity weakens the suppression effects on the leading shock wave while increases the attenuation rate of the reflected shock wave. The first perforated plate in the array plays a major role in the attenuation of the shock wave.     

Flow characteristics of supersonic gas passing through a circular micro-channel under different inflow conditions

Gas flow in a micro-channel usually has a high Knudsen number. The predominant predictive tool for such a microflow is the direct simulation Monte Carlo(DSMC) method, which is used in this paper to investigate primary flow properties of supersonic gas in a circular micro-channel for different inflow conditions, such as free stream at different altitudes, with different incoming Mach numbers, and with different angles of attack. Simulation results indicate that the altitude and free stream incoming Mach number have a significant effect on the whole micro-channel flow field, whereas the angle of attack mainly affects the entrance part of micro-channel flow field. The fundamental mechanism behind the simulation results is also presented. With the increase of altitude, thr free stream would be partly prevented from entering into micro-channel.Meanwhile, the gas flow in micro-channel is decelerated, and the increase in the angle of attack also decelerates the gas flow. In contrast, gas flow in micro-channel is accelerated as free stream incoming Mach number increases. A noteworthy finding is that the rarefaction effects can become very dominant when the free stream incoming Mach number is low. In other words, a free stream with a larger incoming velocity is able to reduce the influence of the rarefaction effects on gas flow in the micro-channel.