Schardin问题的数值研究

激波绕过三角楔(Schardin问题)时会产生激波马赫反射与绕射、 三角楔尾涡与涡串等复杂物理现象. 本文利用三阶精度加权基本无振荡(WENO)格式、 结构化矩形网格的自适应加密方法与沉浸边界法对Schardin问题进行了数值模拟. 数值结果清晰地显示了激波与三角楔相互作用, 在楔面发生马赫反射以及在楔角绕射诱导主涡的过程, 并与Schardin等的实验结果及相关数值结果完全符合. 另外, 数值结果还详细反映了先前实验与数值结果没有详细讨论的主涡滑移层上的涡串生成机理, 以及激波与涡串相互作用和产生声波的     

Schardin问题的数值研究

激波绕过三角楔(Schardin问题)时会产生激波马赫反射与绕射、 三角楔尾涡与涡串等复杂物理现象. 本文利用三阶精度加权基本无振荡(WENO)格式、 结构化矩形网格的自适应加密方法与沉浸边界法对Schardin问题进行了数值模拟. 数值结果清晰地显示了激波与三角楔相互作用, 在楔面发生马赫反射以及在楔角绕射诱导主涡的过程, 并与Schardin等的实验结果及相关数值结果完全符合. 另外, 数值结果还详细反映了先前实验与数值结果没有详细讨论的主涡滑移层上的涡串生成机理, 以及激波与涡串相互作用和产生声波的过程.     

高超声速层流绕三维压缩拐角的数值模拟

通过数值模拟研究了高超声速来流绕过压缩拐角的层流分离三维流动特性.数值方法采用三维N-S方程,结合2阶精度Roe格式以及分区结构网格有限体积法进行离散.数值模拟的空间激波结构与实验纹影结果符合较好;激波/边界层干扰区内3条纵向线上的计算压力分布与实验结果进行了对比分析,计算获得在三维楔侧面存在低压力区,与实验结果反映的规律一致,计算结果表明低压力区是由楔体侧缘尖端发起的二次涡的抽吸作用造成的.此外,在楔体后端尾流区的低压沿边界层内的亚声速区往上游传递了一定距离.      

波后滑移线引起的界面变形和R-M不稳定性

实验研究复杂波形结构引起平面界面变形和反射激波冲击下的R-M不稳定性的问题.在竖直激波管中生成稳定的N₂/SF₆平面界面,激波在圆柱绕射后,冲击平面界面,由此研究复杂激波引起的界面变形.平面激波在圆柱绕射后的流场,演化成具有初始入射波、三波点、弯曲反射波、Mach波和Mach反射产生的滑移线等复杂结构.研究复杂结构激波对界面的作用,对认识界面扰动的生成具有较大帮助.绕柱激波冲击后,平面界面仅在两对滑移线内部发生变形.绕柱激波冲击界面后,两对滑移线将界面分成"内界面"和"外界面",界面变形形态同滑移线和界面相交位置相关.反射激波二次冲击下,界面扰动的增长与Jacobs-Sheeley涡量模型较吻合.      

Interaction of shock waves with a jet wake during gas injection into supersonic stream

In the present paper, we report the results of an experimental study of the interaction region of a planar compression shock produced by a wedge in stream with the wake formed behind a cocurrent gas jet (H₂, air, or Ar) injected into the flow. Depending on the gas jet parameters, three modes of interaction could be distinguished: a strong interaction, observed when the flow velocity in the wake was subsonic; a moderate interaction, observed when a subsonic flow region, bounded by a shock of almost conical shape, formed in the vicinity of the compression shock; and a neutral interaction. Three-dimensional non-stationary Euler equations were solved to numerically examine the interaction of an axisymmetric jet with an oblique shock wave. The obtained interaction regimes were found to be in a reasonable agreement with experimental data.     

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

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

Interaction between a composite compression wave and a vortex in a thermodynamically nonideal medium

A numerical analysis is presented of two-dimensional interaction between a transverse vortex and a composite compression wave that can exist in a thermodynamically nonideal medium. It is shown that the interaction of a composite wave involving a “neutrally stable” shock with a vortex generates weakly damped outgoing acoustic waves; i.e., the shock is a source of sound. This phenomenon increases the post-shock acoustic noise level in an initially turbulent flow.     

Anomalous flow deflection at earth's low-Alfvén-Mach-Number bow shock

Earth's magnetosphere is an obstacle to the supersonic solar wind and the bow shock is formed in the front side of it. In ordinary hydrodynamics, the flow decelerated at the shock is diverted around the obstacle symmetrically about the Earth-Sun line, which is indeed observed in the magnetosheath most of the time. Here we show a case under a very low-density solar wind in which duskward flow was observed in the dawnside magnetosheath. A Rankine-Hugoniot test shows that the magnetic effect is crucial for this "wrong flow" to appear. A full three-dimensional magnetohydrodynamics (MHD) simulation of the situation confirming this interpretation and earlier simulations is also performed. It is illustrated that in addition to the "wrong flow" feature, various peculiar characteristics appear in the global picture of the MHD flow interaction with the obstacle.     

Investigation of magnetohydrodymamic interaction in a supersonic air flow at <Emphasis Type="Italic">M</Emphasis> = 8

The influence of magnetohydrodynamic interaction localized before a model on the position of a shock wave attached to a wedge is experimentally and numerically investigated. The investigation is carried out in an air flow with a Mach number of 8. It is shown that, for a hydromagnetic interaction parameter on the order of 0.1, the slope angle of the shock wave can be increased by 10°. Experiments are conducted for the case when the flow is ionized by an electron beam or by a pulsed electric discharge. Good agreement between experimental and numerical results is obtained for both ways of ionization if the Joule heating of the gas is insignificant. The conclusion is drawn that the way of providing a nonequilibrium conductivity of the flow has a minor effect on the position of the oblique shock wave near the wedge with the hydromagnetic interaction parameter being the same.     

Vortex reconnection in a swirling flow

Processes of vortex reconnection on a helical vortex, which is formed in a swirling flow in a conical diffuser, have been studied experimentally. It has been shown that reconnection can result in the formation of both an isolated vortex ring and a vortex ring linked with the main helical vortex. A number of features of vortex reconnection, including the effects of asymmetry, generation of Kelvin waves, and formation of various bridges, have been described.     

Clustering of the low-inertia particle number density field in random divergence-free hydrodynamic flows

We consider the diffusion of the low-inertia particle number density field in random divergence-free hydrodynamic flows. The principal feature of this diffusion is the divergence of the particle velocity field, which results in clustering of the particle number density field. This phenomenon is coherent, occurs with a unit probability, and must show up in almost all realizations of the process dynamics. We calculate the statistical parameters that characterize clustering in three-dimensional and two-dimensional random fluid flows and in a rapidly rotating two-dimensional random flow. In the former case, the vortex component of the random divergence-free flow generates a vortex component of the low-inertia particle velocity field, which generates a potential component of the velocity field through advection. By contrast, in the case of rapid rotation, a potential component of the velocity field is generated directly by the vortex component of the random divergence-free flow (linear problem).     

Global visualization and quantification of compressible vortex loops

The physics of compressible vortex loops generated due to the rolling up of the shear layer upon the diffraction of a shock wave from a shock tube is far from being understood, especially when shock-vortex interactions are involved. This is mainly due to the lack of global quantitative data available which characterizes the flow. The present study involves the usage of the PIV technique to characterize the velocity and vorticity of compressible vortex loops formed at incident shock Mach numbers ofM=1.54 and1.66. Another perk of the PIV technique over purely qualitative methods, which has been demonstrated in the current study, is that at the same time the results also provide a clear image of the various flow features. Techniques such as schlieren and shadowgraph rely on density gradients present in the flow and fail to capture regions of the flow influenced by the primary flow structure which would have relatively lower pressure and density. Various vortex loops, namely, square, elliptic and circular, were generated using different shape adaptors fitted to the end of the shock tube. The formation of a coaxial vortex loop with opposite circulation along with the generation of a third stronger vortex loop ahead of the primary with same circulation direction are of the interesting findings of the current study.     

Wavenumber-extended high-order oscillation control finite volume schemes for multi-dimensional aeroacoustic computations

A new numerical method toward accurate and efficient aeroacoustic computations of multi-dimensional compressible flows has been developed. The core idea of the developed scheme is to unite the advantages of the wavenumber-extended optimized scheme and M-AUSMPW+/MLP schemes by predicting a physical distribution of flow variables more accurately in multi-space dimensions. The wavenumber-extended optimization procedure for the finite volume approach based on the conservative requirement is newly proposed for accuracy enhancement, which is required to capture the acoustic portion of the solution in the smooth region. Furthermore, the new distinguishing mechanism which is based on the Gibbs phenomenon in discontinuity, between continuous and discontinuous regions is introduced to eliminate the excessive numerical dissipation in the continuous region by the restricted application of MLP according to the decision of the distinguishing function. To investigate the effectiveness of the developed method, a sequence of benchmark simulations such as spherical wave propagation, nonlinear wave propagation, shock tube problem and vortex preservation test problem are executed. Also, throughout more realistic shock–vortex interaction and muzzle blast flow problems, the utility of the new method for aeroacoustic applications is verified by comparing with the previous numerical or experimental results.     

跨音压气机转子叶尖流场试验与分析

本文利用高频响动态压力测试技术测量了某跨音压气机转子叶尖间隙流场,其中包括对跨音转子叶尖漏流和激波／漏流干扰的细致流场特征．测试结果表明激波在机匣壁面处受漏流干扰的影响相当大,这一干扰分别来自于前缘漏流涡的干扰和叶尖第二次漏流对通道激波的直接冲击干扰,使激波结构呈“S”形.在低背压高攻角条件下,转子存在十分明显的尾涡．     

Peculiarities of evolution of shock waves generated by boiling coolant

Simulation of compression wave generation and evolution at the disk target was performed for the case of explosive-type boiling of coolant; the boiling is initiated by endwall rupture of a high-pressure pipeline. The calculations were performed for shock wave amplitude at different times and modes of pipe rupture. The simulated pressure of a target-reflected shock wave is different from the theoretical value for ideal gas; this discrepancy between simulation and theory becomes lower at higher distances of flow from the nozzle exit. Comparative simulation study was performed for flow of two-phase coolant with account for slip flow effect and for different sizes of droplets. Simulation gave the limiting droplet size when the single-velocity homogeneous flow model is valid, i.e., the slip flow effect is insignificant.     

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

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

ILES for mechanism of ramp-type MVG reducing shock induced flow separation

A high order implicit large eddy simulation (ILES) is carried out to study the mechanism of shock induced flow separation reduction under ramp-type MVG control. The mechanism was originally considered as that MVG can generate streamwise vortices which strongly mix boundary layer and the boundary layer becomes more capable to resist strong adverse pressure gradient caused by shock and to keep the boundary layer attached. However, according to our ILES, a chain of ring-like vortices is generated behind the ramp-type MVG and goes further to interact with the shock. When the ring-like vortices pass through the shock, the shock wave is weakened and altered while the vortex structures are quite stable. The instantaneous simulation shows that the spanwise ring-like vortex, not the streamwise vortex, plays a key role to weaken the shock and reduce the shock-induced separation. Detailed investigation on ring-like vortices and shock interaction will be given in this paper.     

Classification of shock wave reflections from a wedge. Part 1: Boundaries and domains of existence for different types of reflections

Classification of various types of the reflections of a shock wave over straight wedge is supposed. The idea about entire reflection phenomenon as a result of interaction of two processes—the shock wave refection process and the flow deflection process—serves as a basis for the classification. To recognize the types of reflection, changes in the shapes of the reflected wave, Mach stem and contact surface are taken into account. The boundaries and domains of existence for various types of reflection are reported. New names for some types of reflection are proposed.     

A hybrid CFD/characteristics method for fast characterization of hypersonic blunt forebody/inlet flow

A hybrid CFD/characteristic method(CCM) was proposed for fast design and evaluation of hypersonic inlet flow with nose bluntness, which targets the combined advantages of CFD and method of characteristics. Both the accuracy and efficiency of the developed CCM were verified reliably, and it was well demonstrated for the external surfaces design of a hypersonic forebody/inlet with nose bluntness. With the help of CCM method, effects of nose bluntness on forebody shock shapes and the flowfield qualities which dominate inlet performance were examined and analyzed on the two-dimensional and axisymmetric configurations. The results showed that blunt effects of a wedge forebody are more substantial than that of related cone cases. For a conical forebody with a properly blunted nose, a recovery of the shock front back to that of corresponding sharp nose is exhibited, accompanied with a gradually fading out of entropy layer effects. Consequently a simplification is thought to be reasonable for an axisymmetric inlet with a proper compression angle, and a blunt nose of limited radius can be idealized as a sharp nose, as the spillage and flow variations at the entrance are negligible, even though the nose scale increases to 10% cowl lip radius. Whereas for two-dimensional inlets, the blunt effects are substantial since not only the inlet capturing/starting capabilities, but also the flow uniformities are obviously degraded.