Vortex structure generation on the frontal surface of a cylinder set transversely in a hypersonic flow

The problem of formation of spatially periodic structures on the frontal surface of a cylindrically blunted body set transversely in a hypersonic flow is studied. Within the framework of the model adopted, a possible mechanism of vortex structure generation on the frontal surface of the blunt body is proposed and confirmed by calculations; in this mechanism, the curved bow shock produces a vortex flow, while in its turn the vortex, which persists under weak dissipation, acts on the shock thus maintaining its curved shape. It is shown that the spatially periodic mode of hypersonic flow past a cylinder can exist in the case of a uniform incident flow and under homogeneous boundary conditions on the body surface.     

Experimental modeling of three-dimensional modes of hypersonic flow past cylindrically-blunted bodies

Several sets of experimental studies of the structure of transverse hypersonic flow past blunt bodies (cylinder and truncated wedge) and heat transfer on them are performed in the UT-1M shock tube of the Central Aerohydrodynamics Institute. The purpose of the investigation was to obtain three-dimensional modes of hypersonic flow past the nose surfaces of blunt bodies in an artificially disturbed and nominally uniform flows. The controlled disturbances in the freestream were produced by thin threads pulled over the nozzle exit. In the experiment the flow was visualized using the Töpler method and the heat flux distribution over the cylinder was measured using luminescent temperature transformers. The experiments show that both the flow and the heat transfer in the vicinity of the cylinder nose are very sensitive to vortex disturbances in the oncoming hypersonic flow. In a nominally uniform flow (M_∞ = 8 and Re_∞ = 3160–11670) a steady three-dimensional mode of flow past the nose surface of a blunt wedge could be obtained in the form of a single vortex pair.     

火星科学实验室气动特性数值分析

针对火星着陆探测器进入-下降-着陆过程的高超声速进入阶段,求解三维流体动力学Navier-Stokes 方程与化学反应动力学模型,分析火星科学实验室进入火星大气时的化学非平衡效应、探测器周围的流场结构和气动特性在化学非平衡效应影响下的变化. 结果表明,CO₂ 在激波后大量分解,消耗大量能量;在化学非平衡效应影响下,探测器头部激波脱体距离大幅减小,尾迹旋涡运动减弱;化学非平衡效应影响下探测器升力系数变化不大,阻力系数高于完全气体,升阻比略低,配平攻角小于完全气体.     

Formation and evolution of vortex rings induced by interactions between shock waves and a low-density bubble

The deformation and instability of a low-density spherical bubble induced by an incident and its reflected shock waves are studied by using the large eddy simulation method. The computational model is firstly validated by experimental results from the literature and is further used to examine the effect of incident shock wave strength on the formations and three-dimensional evolutions of the vortex rings. For the weak shock wave case (Ma?=?1.24), the baroclinic effect induced by the reflected shock wave is the key mechanism for the formation of new vortex rings. The vortex rings not only move due to the self-induced effect and the flow field velocity, but also generate azimuthal instability due to the pressure disturbance. For the strong shock wave case (Ma?=?2.2), a boundary layer is formed adjacent to the end wall owing to the approach of vortex ring, and unsteady separation of the boundary layer near the wall results in the ejection and formation of new vortex rings. These vortex rings interact in the vicinity of the end wall and finally collapse to a complicated vortex structure via azimuthal instability. For both shock wave strength cases, the evolutions of vortex rings due to the instability lead to the formation of the complicated structure dominated by the small-scale streamwise vortices.     

EXPERIMENTAL STUDY ON QUANTITATIVE MEASUREMENT OF THREE-DIMENSIONAL STRUCTURE OF HAIRPIN VORTEX BY STEREO-PIV 1)

发卡涡是湍流相干结构研究中最为关注的内容,实现发卡涡三维结构的定量测量并进行流体动力学分析,对深入研究湍流相干结构、实现湍流精准控制等具有重要意义.本研究通过对合成射流装置进行合理控制,使得层流边界层中产生了规则的人造发卡涡结构,进而用体视图像粒子测速仪(Stereo-PIV)锁相实验技术对发卡涡结构所在的三维空间流场进行了定量测量,并得到了一个完整周期内形成的发卡涡三维结构的空间流场. 结果发现,重构所得的三维发卡涡结构质量较高, 实验技术和方案具有可行性.发卡涡结构所在空间流场情况,符合目前人们对于发卡涡、高低速条带、喷射和扫掠事件的常规认识. 此外,对近壁二次流向涡、展向涡量集中区域的展向涡头和强剪切区域、与低速喷射流体相关的汇聚流动和发散流动等有了更细致的认识.同时, 也探讨了"基于二维脉动流场的相关特征去重构发卡涡三维流场"的可行性.为进一步定量探究发卡涡结构的形成演化、不同涡结构的融合及二次诱导等壁湍流相干结构问题提供思路.     

基于Stereo-PIV技术的三维发卡涡结构定量测量研究

发卡涡是湍流相干结构研究中最为关注的内容,实现发卡涡三维结构的定量测量并进行流体动力学分析,对深入研究湍流相干结构、实现湍流精准控制等具有重要意义.本研究通过对合成射流装置进行合理控制,使得层流边界层中产生了规则的人造发卡涡结构,进而用体视图像粒子测速仪(Stereo-PIV)锁相实验技术对发卡涡结构所在的三维空间流场进行了定量测量,并得到了一个完整周期内形成的发卡涡三维结构的空间流场. 结果发现,重构所得的三维发卡涡结构质量较高, 实验技术和方案具有可行性.发卡涡结构所在空间流场情况,符合目前人们对于发卡涡、高低速条带、喷射和扫掠事件的常规认识. 此外,对近壁二次流向涡、展向涡量集中区域的展向涡头和强剪切区域、与低速喷射流体相关的汇聚流动和发散流动等有了更细致的认识.同时, 也探讨了"基于二维脉动流场的相关特征去重构发卡涡三维流场"的可行性.为进一步定量探究发卡涡结构的形成演化、不同涡结构的融合及二次诱导等壁湍流相干结构问题提供思路.      

Three dimensional shock wave/boundary layer interactions

An investigation into a three-dimensional, curved shock wave interacting with a three-dimensional, curved boundary layer on a slender body is presented. Three different nose profiles mounted on a cylindrical body were tested in a supersonic wind tunnel and numerically simulated by solving the Navier–Stokes equations. The conical and hemispherical nose profiles tested were found to generate shock waves of sufficient strength to separate the boundary layer on the cylinder, while the shock wave generated by the ogival profile did not separate the boundary layer. For the separated flow, separation was found to occur predominantly on the windward side of the cylinder with the lee-side remaining shielded from the direct impact of the incident shock wave. A thickening of the boundary layer on the lee-side of all the profiles was observed, and in the conical and hemispherical cases this leads to the re-formation of the incident shock wave some distance away from the surface of the cylinder. A complex reflection pattern off the shock wave/boundary layer interaction (SWBLI) was also identified for the separated flow cases. For comparative purposes, an inviscid simulation was performed using the hemispherical profile. Significant differences between the viscous and inviscid results were noted including the absence of a boundary layer leading to a simplified shock wave reflection pattern forming. The behaviour of the incident shock wave on the lee-side of the cylinder was also affected with the shock wave amalgamating on the surface of the cylinder instead of away from the surface as per the viscous case. Test data from the wind tunnel identified two separation lines present on the cylindrical surface of the hemispherical SWBLI generator. The pair of lines were not explicitly evident in the original CFD simulations run, but were later identified in a high-resolution simulation.     

Characteristics of unsteady type IV shock/shock interaction

Characteristics of the unsteady type IV shock/shock interaction of hypersonic blunt body flows are investigated by solving the Navier–Stokes equations with high-order numerical methods. The intrinsic relations of flow structures to shear, compression, and heating processes are studied and the physical mechanisms of the unsteady flow evolution are revealed. It is found that the instantaneous surface-heating peak is caused by the fluid in the “hot spot” generated by an oscillating and deforming jet bow shock (JBS) just ahead of the body surface. The features of local shock/boundary layer interaction and vortex/boundary layer interaction are clarified. Based on the analysis of flow evolution, it is identified that the upstream-propagating compression waves are associated with the interaction of the JBS and the shear layers formed by a supersonic impinging jet, and then the interaction of the freestream bow shocks and the compression waves results in entropy and vortical waves propagating to the body surface. Further, the feedback mechanism of the inherent unsteadiness of the flow field is revealed to be related to the impinging jet. A feedback model is proposed to reliably predict the dominant frequency of flow evolution. The results obtained in this study provide physical insight into the understanding of the mechanisms relevant to this complex flow.     

Computational study of curved shock–vortex interaction

The interaction between a curved shock wave and a compressible vortex is numerically investigated. The investigation concentrates on the local deformation of the shock structure due to the shock–vortex interaction. The essentially non‐oscillatory (ENO) scheme is used to solve the unsteady two‐dimensional Euler equations. A curved shock wave is obtained by the diffraction of an initially planar shock wave around a right‐angled corner and then allowed to interact with a strong compressible vortex superimposed on the flow. The same vortex affects the shock wave differently depending on the placement of the vortex because of the varying strength of the shock wave. This effect could range from a non‐symmetric deformation of the shock wave to a local disruption in the shock structure depending on the strength of the shock wave in the interaction region. This process leading to a local disruption in the shock structure is analyzed in detail. It is shown that such a disruption in the shock structure can be predicted by simple one‐dimensional considerations. Copyright © 1999 John Wiley & Sons, Ltd.     

Three-dimensional supersonic flow past blunt bodies with contour discontinuities,with account for equilibrium and frozen state of the gas in the shock layer

An algorithm is devised for calculating by the finite difference method the supersonic flow region for three-dimensional steady-state flow of a viscous gas past a blunted body with many contour discontinuities. The state of this gas at high hypersonic flight speeds can be characterized by equilibrium or frozen physicochemical processes. Generally speaking, any arbitrary number and sequence of either compression or expansion discontinuities is permitted. The computational scheme adopted provides identification of the vortex layers, regions with different equations of state, and so on. We use a flow model that is either frozen throughout the entire shock layer or only in the portion of the layer adjacent to the body surface. The pressure at certain points on the surface of spherically blunted cones with half-angles θ?10° may differ by a factor of 2 or more in equilibrium and frozen flows. Example calculations are presented, and the results are analyzed.     

Two-dimensional self-similar flows generated by the interaction between a shock and low-density gas regions

A plane time-dependent flow generated by the interaction between a normal shock and a low-density gas region occupying a quarter of the plane is theoretically investigated. Numerical simulation is performed on the basis of the Euler equations. It is established that after the shock has come in contact with the low-density region two-dimensional self-similar flows of different type can develop. On regular interaction the original shock is refracted on the low-density region with the matching of the accelerated and original shock and the refracted contact discontinuity at a common point. On irregular interaction a complicated flow occurs; it includes curved and oblique shocks, a contact discontinuity with points of inflection, multiple matching points, a high-pressure jet, and a layered vortex. The jet and vortex structures are investigated in detail. The tendency of the gasdynamic structure development with variation in the control parameters of the problem is determined. A simplified, near-analytical technique for estimating the slopes of the main shocks and the gas parameters behind them is proposed.     

Method of mean-mass parameters for a three-dimensional boundary layer in a flow with vorticity

When a gas flows with hypersonic velocity over a slender blunt body, the bow shock induces large entropy gradients and vorticity near the wall in the disturbed flow region (in the high-entropy layer) [1]. The boundary layer on the body develops in an essentially inhomogeneous inviscid flow, so that it is necessary to take into account the difference between the values of the gas parameters on the outer edge of the boundary layer and their values on the wall in the inviscid flow. This vortex interaction is usually accompanied by a growth in the frictional stress and heat flux at the wall [2, 3]. In three-dimensional flows in which the spreading of the gas on the windward sections of the body causes the high-entropy layer to become narrower, the vortex interaction can be expected to be particularly important. The first investigations in this direction [4–6] studied the attachment lines of a three-dimensional boundary layer. The method proposed in the present paper for calculating the heat transfer generalizes the approach realized in [5] for the attachment lines and makes it possible to take into account this effect on the complete surface of a blunt body for three-dimensional laminar, transition, or turbulent flow regime in the boundary layer.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 80–87, January–February, 1981.     

高超声速稀薄流中横向喷流干扰特性实验研究

喷流干扰是高超声速飞行高精度控制的一种有效手段,研究者们以往大部分都主要集中于连续流条件下喷流干扰效应的机理研究,并给出了喷流干扰流场的典型结构,而稀薄流条件下喷流干扰特性的实验数据还十分匮乏.本文利用JFX爆轰激波风洞产生高超声速稀薄自由流,基于平板模型开展不同喷流压力和自由来流参数对横向喷流干扰特性影响的实验研究,采用高速纹影成像及图像处理技术,获得稀薄流条件下喷流干扰流场演化过程及流场结构的变化规律.相比于无喷流条件形成的流场,横向喷流与稀薄自由流相互作用形成的流场结构更为复杂,喷流压力由于受到稀薄来流的扰动,斜激波会短暂穿透喷流干扰流场并延伸至楔形体上部.喷流干扰流场内桶状激波的影响范围随着喷流压力的升高而逐渐变宽,位于三波点上游的斜激波空间位置不会随喷流压力的变化而改变,而位于三波点下游的弓形激波则向上游移动,当喷流压力过低时,桶状激波不会与其他两种激波交汇形成三波点.高超声速稀薄来流压力的降低同样会使桶状激波的影响范围变宽,弓形激波同样也会向上游移动,但基本不会对斜激波空间位置产生任何影响.     

Nonequilibrium hypersonic flow of gas past a wing of small aspect ratio

It is well known [1] that nonequilibrium physicochemical processes taking place in gases at high temperature influence the gas-dynamic parameters and aerodynamic characteristics of bodies in hypersonic flight. In the present paper, the thin shock layer method [2–4] is used to consider the problem of nonequilibrium hypersonic flow of gas past a wing of small aspect ratio at an angle of attack. It is shown that the flow component of the vorticity is conserved along the streamlines, and this property is exploited to obtain an analytic solution of the equations of the three-dimensional nonequilibrium shock layer. The influence of the disequilibrium on the thickness of the shock layer and the pressure distribution is investigated.     

可压缩燃烧反应转捩混合层直接数值模拟

针对三维时间发展可压缩氢/氧非预混燃烧反应平面自由剪切混合层,采用5阶迎风/6阶对称紧致混合差分格式以及3阶显式Runge-Kutta时间推进方法,直接数值模拟了伴随燃烧产物生成和反应能量释放, 流动受扰动激发失稳并转捩的演化过程. 在转捩初期, 获得了${\it\Lambda}$涡、马蹄涡等典型的大尺度拟序结构,观察到了流动失稳后发生双马蹄涡三维对并的现象, 大尺度结构呈较好的对称性.在流动演化后期, 大尺度结构逐次破碎形成小尺度结构, 混合层进入转捩末期,呈明显的不对称性.      

关于吸气式高超声速推进技术研究的思考

回顾了吸气式高超声速推进技术的研究进展, 分析了超燃冲压发动机研制面临的关键科学问题, 并从不同角度探讨了增大超燃冲压发动机推力的可能方法.这些方法包括: 能够降低总压损失的高超声速来流压缩方法、生成三维涡流的超声速混合增强技术、碳氢燃料的预热喷射、可以控制燃烧过程的燃烧室设计优化方法、通过减小发动机流道湿面积来降低摩擦阻力和催化复合解离的燃气降低高温气体效应.考虑到等压热力学循环的热效率,还建议研究在高超声速推进系统中应用热效率高的爆轰过程, 并探讨了爆轰推进方法研究的进展与问题.吸气式高超声速推进技术是高超声速飞行器发展的关键技术, 认真思考和探索其发展方向是非常必要的.      

Hypersonic flow of nonequilibrium air over blunt bodies

A numerical method is described for computing nonequilibrium three-dimensional supersonic flow of a gas in the shock layer over the forward surface of blunt bodies with discontinuities of shape. The basic idea is to divide the original system of differential equations into two subsystems, which are solved in succession: first for the gasdynamic variables, the velocity components and the pressure, and then for the relaxation parameters and the enthalpy. To calculate the velocity components and the pressure we use the iterative marching method [1, 2] in the form given in [3]. The relaxation equations and the enthalpy equation are integrated numerically along the stream lines. A discussion is given of the effect of nonequilibrium of physical and chemical reactions on the distribution of parameters in the inviscid shock layer and on the aerodynamic coefficients of blunt bodies in hypersonic air flow. The unsteady aerodynamic coefficients are calculated by the curved body method [4]. The computational algorithm takes the form of a program in “ALGOL-60” for the BéSM-6 computer.     

Baroclinic vortex sheet production by shocks and expansion waves

Vortex sheet production by shocks and expansion waves refracting at a density discontinuity was examined and compared using an analytical solution and numerical simulations. The analytical solution showed that with a small exception, vortex sheet strength is generally stronger in fast/slow shock refractions. In contrast, expansion waves generated a stronger vortex sheet in slow/fast refractions. This difference results in larger vorticity deposited by shocks in fast/slow refractions and by expansion waves in slow/fast refractions. Shock refractions become irregular and the analytical solution fails when either incident, transmitted or reflected shock, exceeded the angle limit for an attached shock. To investigate vortex sheet production outside the range of analytical solutions and to verify the applicability of the planar-interface analytical solution to a curved interface, shock refraction through a sinusoidal interface was numerically simulated in the shock frame of reference. It is found that variation in the local incidence angle along the curved interface creates pressure waves that affect the level of deposited vorticity. This contributes to the difference between predictions from local analysis and numerical computation. Furthermore, an interesting behavior of the shock and expansion wave-deposited vorticity in supersonic ramp flow was discovered. When the high- and low-density streams were swapped, while keeping the incident flow Mach numbers constant, a vortex sheet of equal magnitude but of opposite sign was generated.     

Solution of the Chaplygin equation in plane-parallel transonic flow with a normal shock

In [1], for the case in which the exact Chaplygin equation of motion is replaced by the Tricomi equation, Frankl constructed an example of plane-parallel transonic flow with a normal shock terminating within the flow (later this example was extended by Biibosunov to obtain a curved shock [2]). This flow is an analog of the point vortex flow and has concentric circles as streamlines, and a unique dependence of the velocity vector on the Cartesian coordinates is provided everywhere. In the present paper this example is extended to the case of the Chaplygin equation, and the shock is considered to be normal.The author wishes to thank S. V. Fal'kovich for discussion of the study results and helpful advice.