HLLE++格式在高马赫数空腔流动模拟中的应用

空腔流动存在剪切层运动、涡脱落与破裂,以及激波与激波、激波与剪切层、激波与膨胀波和激波/涡/剪切层相互干扰等现象,流动非常复杂,特别是高马赫数(M>2)时,剪切层和激波更强,激波与激波干扰更严重,对数值格式的要求更高,既需要格式耗散小,对分离涡等有很高的模拟精度,又需要格式在激波附近具有较大的耗散,可以很好地捕捉激波,防止非物理解的出现。Roe和HLLC等近似Riemann解格式在高马赫数强激波处可能会出现红玉现象,而HLLE++格式大大改善了这种缺陷,在捕捉高超声速激波时避免了红玉现象的发生,同时还保持在光滑区域的低数值耗散特性。本文在结构网格下HLLE++格式的基础上,通过改进激波探测的求解,建立了基于非结构混合网格的HLLE++计算方法,通过无粘斜坡算例,验证了HLLE++格式模拟高马赫数流动的能力,并应用于高马赫数空腔流动的数值模拟,开展了网格和湍流模型影响研究,验证了方法模拟高马赫数空腔流动的可靠性和有效性。     

Evolution of the Diffusion Mixing Layer of Two Gases upon Interaction with Shock Waves

A mathematical model of mechanics of a twovelocity twotemperature mixture of gases is developed. Based on this model, evolution of the mixing layer of two gases with different densities under the action of shock and compression waves is considered by methods of mathematical simulation in the onedimensional unsteady approximation. In the asymptotic approximation of the full model, a solution of an initialboundary problem is obtained, which describes the formation of a diffusion layer between two gases. Problems of interaction of shock and compression waves with the diffusion layer are solved numerically in the full formulation. It is shown that the layer is compressed as the shock wave traverses it; the magnitude of compression depends on shockwave intensity. As the shock wave passes from the heavy gas to the light gas, the mixing layer becomes overcompressed and expands after shockwave transition. The wave pattern of the flow is described in detail. The calculated evolution of the mixinglayer width is in good agreement with experimental data.     

Shear layer behavior resulting from shock wave diffraction

All previous studies on shock wave diffraction in shock tubes have spatial and temporal limitations due to the size of the test sections. These limitations result from either the reflection of the expansion wave, generated at the corner, from the top wall and/or of the reflection of the incident diffracted shock from the bottom wall of the test section passing back through the region of interest. This has limited the study of the evolution of the shear layer and its associated vortex, which forms a relatively small region of the flow behind the shock with an extent of only a few centimeters, and yet is a region of significant interest. A special shock tube is used in the current tests which allow evolution of the flow to be examined at a scale about an order of magnitude larger than in previously published results, with shear layer lengths of up to 250 mm being achieved without interference from adjacent walls. Tests are presented for incident shock wave Mach numbers of nominally 1.3–1.5. Studies have been undertaken with wall angles of 10, 20, 30 and 90°. Significant changes are noted as the spatial and temporal scale of the experiment increases. For a given wall angle, the flow behind the incident shock is not self-similar as is usually assumed. Both shear layer instability and the development of turbulent patches become evident, neither of which have been noted in previous tests.     

二维激波与剪切层相互作用的直接数值模拟研究

采用五阶weighed esseritially non-oscillatory (WENO) 格式和三阶total variation diminishing (TVD) Runge-Kutta 格式, 通过求解二维非定常Navier-Stokes 方程, 直接数值模拟了激波与剪切层相互作用, 目的在于揭示激波与剪切层相互作用过程中噪声产生的机理. 研究发现:(1) 当入射激波穿过剪切层时, 剪切层中心位置向下层区域偏移;(2) 入射激波穿过剪切层产生小激波, 在小激波与剪切层接触点处产生声波并向外辐射;(3) 反射激波穿过剪切层后形成了分段弧状激波;(4) 当反射激波穿过剪切层时, 激波在鞍点处泄漏并向外辐射声波, 这是一种激波泄漏机制.     

A calculating method of shock wave oscillating frequency due to turbulent shear layer fluctuations in supersonic flow

One of the more severe fluctuating pressure environments encountered in supersonic orhypersonic flows is the shock wave oscillation driven by interaction of a shock wave withboundary layer.The high intensity oscillating shock wave may induce structure resonanceof a high speed vehicle.The research for the shock oscillation used to adopt empirical orsemiempirical methods because the phenomenon is very complex.In this paper atheoretical solution on shock oscillating frequency due to turbulent shear layer fluctuationshas been obtained from basic conservation equations.Moreover,we have attained theregularity of the frequency of oscillating shock varying with incoming flow Mach numbersM_∞and turning angleθ.The calculating results indicate excellent agreement withmeasurements.This paper has supplied a valuable analytical method to study aeroelasticproblems produced by shock wave oscillation.     

Particle image velocimetry measurements of a shock wave/turbulent boundary layer interaction

Particle image velocimetry is used to investigate the interaction between an incident shock wave and a turbulent boundary layer at Mach 2.1. A particle response assessment establishes the fidelity of the tracer particles. The undisturbed boundary layer is characterized in detail. The mean velocity field of the interaction shows the incident and reflected shock wave pattern, as well as the boundary layer distortion. Significant reversed flow is measured instantaneously, although, on average no reversed flow is observed. The interaction instantaneously exhibits a multi-layered structure, namely, a high-velocity outer region and a low-velocity inner region. Flow turbulence shows the highest intensity in the region beneath the impingement of the incident shock wave. The turbulent fluctuations are found to be highly anisotropic, with the streamwise component dominating. A distinct streamwise-oriented region of relatively large kinematic Reynolds shear stress magnitude appears within the lower half of the redeveloping boundary layer. Boundary layer recovery towards initial equilibrium conditions appears to be a gradual process.     

Analytical solution of a problem of a collision of two hypersonic gas flows from symmetric sources

An asymptotic solution of the Euler equations that describe stationary interaction of two hypersonic gas flows from two identical spherically symmetric sources and an integral equation determining the shock wave shape are obtained with the use of a modified method of expansion of the sought functions with respect to a small parameter, which is the ratio of gas densities in the incoming flow and behind the shock wave. The solution of this equation near the axis of symmetry allows the shock wave stand-off distance from the contact plane and the radius of its curvature to be found. It is shown that the solution obtained agrees well with the known numerical solutions.     

Supersonic low-concentration dusty-gas flow past a plane cylinder in the presence of an oblique shock wave interacting with the bow shock

The motion of an inertial dispersed admixture near a plane cylinder immersed in a steady-state hypersonic dusty flow in the presence of an oblique shock wave interacting with the bow shock is considered. It is assumed that the free-stream particle mass concentration is small and the particles do not affect the carrier flow. The III and IV shock wave interaction regimes are considered. The gas flow parameters in the shock layer are calculated from the numerical solution of the full Navier-Stokes equations for the perfect gas. A TVD second-order finite-difference scheme constructed on the basis of a finite volume method is used. For calculating the dispersed-phase parameters, including the concentration, the full Lagrangian method is used. On a wide range of variation of the particle inertia parameters, the patterns of the particle trajectories, velocity, concentration, and temperature in the shock layer are studied. The possibility of aerodynamic focusing of the particles behind the shock wave intersection point and the formation of narrow beams with a high particle concentration is revealed. These beams impinge on the cylinder surface and result in a sharp increase in the local heat fluxes. The maximal possible increase in the heat fluxes caused by the particles colliding with the cylinder surface is estimated for the flows with and without the incident oblique shock wave.     

Interaction of a high-speed combustion front with a closely packed bed of spheres

The head-on collision of a combustion front with a closely packed bed of ceramic-oxide spheres was investigated in a vertical 76.2 mm diameter tube containing a nitrogen diluted stoichiometric ethylene–oxygen mixture. A layer of spherical beads in the diameter range of 3–12.7 mm was placed at the bottom of the tube and a flame was ignited at the top endplate. Four orifice plates spaced at one tube diameter were placed at the ignition end of the tube in order to accelerate the flame to either a “fast-flame” or a detonation wave before the bead layer face. The mixture reactivity was adjusted by varying the initial mixture pressure between 10 and 100 kPa absolute. The pressure before and within the bead layer was measured by flush wall-mounted pressure transducers. For initial pressures where a fast-flame interacts with the bead layer peak pressures recorded at the bead layer face were as high as five times the reflected Chapman–Jouget detonation pressure. The explosion resulting from the interaction developed by two distinct mechanisms; one due to the shock reflection off the bead layer face, and the other due to shock transmission and mixing of burned and unburned gas inside the bead layer. The measured explosion delay time (time after shock reflection from the bead layer face) was found to be independent of the incident shock velocity. As a result, the explosion initiation is not the direct result of the shock reflection process but instead is more likely due to the interaction of the reflected shock wave and the trailing flame. The bead layer was found to be very effective in attenuating the explosion front transmitted through the bead layer and thus isolating the tube endplate. This paper is based on work that was presented at the 21th International Colloquium on the Dynamics of Explosions and Reactive Systems, Poitiers, France, July 23–27, 2007.     

The shock–vortex interaction patterns affected by vortex flow regime and vortex models

We have used a third-order essentially non-oscillatory method to obtain numerical shadowgraphs for investigation of shock–vortex interaction patterns. To search different interaction patterns, we have tested two vortex models (the composite vortex model and the Taylor vortex model) and as many as 47 parametric data sets. By shock–vortex interaction, the impinging shock is deformed to a S-shape with leading and lagging parts of the shock. The vortex flow is locally accelerated by the leading shock and locally decelerated by the lagging shock, having a severely elongated vortex core with two vertices. When the leading shock escapes the vortex, implosion effect creates a high pressure in the vertex area where the flow had been most expanded. This compressed region spreads in time with two frontal waves, an induced expansion wave and an induced compression wave. They are subsonic waves when the shock–vortex interaction is weak but become supersonic waves for strong interactions. Under a intermediate interaction, however, an induced shock wave is first developed where flow speed is supersonic but is dissipated where the incoming flow is subsonic. We have identified three different interaction patterns that depend on the vortex flow regime characterized by the shock–vortex interaction.      

Wave induced thermal boundary layers in a compressible fluid: analysis and numerical simulations

The response of a semi-infinite compressible fluid to a step-wise change in temperature of its boundary is investigated analytically and numerically. Numerical results of the boundary layer structure are compared with Clarke’s analytical solution for a gas with thermal conductivity proportional to temperature. To avoid unwanted numerical dissipation in the numerical analysis, the space-time conservation element and solution element (CESE) method has been adopted to solve the unsteady 1-D Navier-Stokes equations. Good agreement between analytical and numerical results has been found for the development of the thermal boundary layer on a long time scale. Weak shock waves and expansion waves induced by the thermal boundary layer due to its compressibility, are observed in the numerical simulation. Finally, the numerical method has been applied to the reflection of a non-linear expansion wave and to a shock wave from an isothermal wall, thereby illustrating the effect of the boundary layer on the external flow field.     

改进虚拟边界算法在超声速流动问题求解中的应用

提出了一种改进的虚拟单元浸没边界法, 并与一种高阶格式的有限差分算法相结合, 运用于求解超声速复杂几何绕流问题.该算法的核心思想在于在固体边界的内部和外部分别施加满足边界关系的作用点, 使得几何边界离散更加细化, 起到了壁面附近网格局部加密的作用.采用源空间内流体点作为反距离插值算法的重构点, 有效避免了插值点数目过少而与作用点相重合情况.通过对二维激波反射现象 (马赫数为 2.81) 和三维超声速球体绕流问题 (马赫数为 1.2) 的数值模拟, 与实验结果对比表明, 本文改进算法相对一般的虚拟边界法来说能显著提高数值精度, 减小计算误差.计算结果揭示了球体绕流中剪切层、压缩波系和尾迹的相互作用导致自由剪切层失稳的机理.剪切层厚度和湍流雷诺脉动经历了线性增长、大幅度震荡和小幅度波动三个阶段, 导致剪切层表面褶皱因子变化呈指数规律增长.其湍流结构表现出明显的各向异性, 具体在流向雷诺正应力在湍流脉动中占主导地位, 激波的压缩作用对不同方向雷诺正应力的影响存在空间迟滞效应.      

高超声速激波湍流边界层干扰直接数值模拟研究

高超声速激波与湍流边界层干扰会导致飞行器表面出现局部热流峰值,严重影响飞行器气动性能和飞行安全. 针对高马赫数激波干扰问题,以往数值研究多采用雷诺平均方法,而在直接数值模拟方面的相关工作较为少见. 开展高超声速激波与湍流边界层干扰的直接数值模拟研究,有助于进一步提升对其复杂流动机理认识和理解,同时也将为现有湍流模型和亚格子应力模型的改进提供理论依据. 采用直接数值模拟方法对来流马赫数6.0,34°压缩拐角内激波与湍流边界层的干扰问题进行了研究. 基于雷诺应力各向异性张量,分析了高超声速湍流边界层在压缩拐角内的演化特性. 通过对湍动能输运方程的逐项分析,系统地研究了可压缩效应对湍动能及其输运的影响机制. 采用动态模态分解方法,探讨了干扰流场的非定常运动历程. 研究结果表明,随着湍流边界层往下游发展,近壁湍流的雷诺应力状态由两组元轴对称状态逐渐演化为两组元状态,外层区域则由轴对称膨胀趋近于各向同性. 干扰流场内存在强内在压缩性效应(声效应),其对湍动能输运的影响主要体现在压力--膨胀项,而对膨胀--耗散项影响较小. 高超声速下压缩拐角内的非定常运动仍存在以分离泡膨胀/收缩为特征的低频振荡特性,其物理机制与分离泡剪切层密切相关.      

Time-resolved stereo PIV measurements of shock–boundary layer interaction on a supercritical airfoil

Time-resolved stereo particle-image velocimetry (TR-SPIV) and unsteady pressure measurements are used to analyze the unsteady flow over a supercritical DRA-2303 airfoil in transonic flow. The dynamic shock wave–boundary layer interaction is one of the most essential features of this unsteady flow causing a distinct oscillation of the flow field. Results from wind-tunnel experiments with a variation of the freestream Mach number at Reynolds numbers ranging from 2.55 to 2.79 × 10⁶ are analyzed regarding the origin and nature of the unsteady shock–boundary layer interaction. Therefore, the TR-SPIV results are analyzed for three buffet flows. One flow exhibits a sinusoidal streamwise oscillation of the shock wave only due to an acoustic feedback loop formed by the shock wave and the trailing-edge noise. The other two buffet flows have been intentionally influenced by an artificial acoustic source installed downstream of the test section to investigate the behavior of the interaction to upstream-propagating disturbances generated by a defined source of noise. The results show that such upstream-propagating disturbances could be identified to be responsible for the upstream displacement of the shock wave and that the feedback loop is formed by a pulsating separation of the boundary layer dependent on the shock position and the sound pressure level at the shock position. Thereby, the pulsation of the separation could be determined to be a reaction to the shock motion and not vice versa.     

The flow structure of dilute gas-particle suspensions behind a shock wave moving along a flat surface

The asymptotic and numerical investigations of shock-induced boundary layers in gas-particle mixtures are presented.The Saffman lift force acting on a particle in a shear flow istaken into account.It is shown that particle migration across the boundary layer leads tointersections of particle trajectories.The corresponding modification of dusty gas model isproposed in this paper.The equations of two-phase sidewall boundary layer behind a shock wave moving at aconstant speed are obtained by using the method of matched asymptotic expansions.Themethod of the calculation of particle phase parameters in Lagrangian coordinates isdescribed in detail.Some numerical results for the case of small particle concentration aregiven.     

On the Propagation of Long‐Wave Perturbations in a Two‐Layer Free‐Boundary Rotational Fluid

A mathematical model for the propagation of longwave perturbations in a freeboundary shear flow of an ideal stratified twolayer fluid is considered. The characteristic equation defining the velocity of perturbation propagation in the fluid is obtained and studied. The necessary hyperbolicity conditions for the equations of motion are formulated for flows with a monotonic velocity profile over depth, and the characteristic form of the system is calculated. It is shown that the problem of deriving the sufficient hyperbolicity conditions is equivalent to solving a system of singular integral equations. The limiting cases of weak and strong stratification are studied. For these models, the necessary and sufficient hyperbolicity conditions are formulated, and the equations of motion are reduced to the Riemann integral invariants conserved along the characteristics.     

激波与转捩边界层干扰非定常特性数值分析

激波与边界层干扰的非定常问题是高速飞行器气动设计中基础研究内容之一.以往研究主要针对层流和湍流干扰,在分离激波低频振荡及其内在机理方面存在着上游机制和下游机制两类截然不同的理论解释.分析激波与转捩边界层干扰下非定常运动现象有助于进一步加深理解边界层状态以及分离泡结构对低频振荡特性的影响规律,为揭示其产生机理指出新的方向.采用直接数值模拟方法对来流马赫数2.9,24?压缩拐角内激波与转捩边界层干扰下激波的非定常运动特性进行了数值分析.通过在拐角上游平板特定的流向位置添加吹吸扰动激发流动转捩,使得进入拐角的边界层处于转捩初期阶段.在验证了计算程序可靠性的基础上,详细分析了转捩干扰下激波运动的间歇性和振荡特征,着重研究了分离泡展向三维结构对激波振荡特性的影响规律,最后还初步探索了转捩干扰下激波低频振荡产生的物理机制.研究结果表明:分离激波的非定常运动仍存在强间歇性和低频振荡特征,其时间尺度约为上游无干扰区内脉动信号特征尺度的10倍量级;分离泡展向三维结构不会对分离激波的低频振荡特征产生实质影响.依据瞬态脉动流场的低通滤波结果,转捩干扰下激波低频振荡的诱因来源于拐角干扰区下游,与流场中分离泡的收缩/膨胀运动存在一定的关联.     

A Newton/upwind method and numerical study of shock wave/boundary layer interactions

The objective of the paper is twofold. First we describe an upwind/central differencing method for solving the steady Navier–Stokes equations. The symmetric line relaxation method is used to solve the resulting algebraic system to achieve high computational efficiency. The grid spacings used in the calculations are determined from the triple-deck theory, in terms of Mach and Reynolds numbers and other flow parameters. Thus the accuracy of the numerical solutions is improved by comparing them with experimental, analytical and other computational results. Secondly we proceed to study numerically the shock wave/boundary layer interactions in detail, with special attention given to the flow separation. The concept of free interaction is confirmed. Although the separated region varies with Mach and Reynolds numbers, we find that the transverse velocity component behind the incident shock, which has not been identified heretofore, is also an important parameter. A small change of this quantity is sufficient to eliminate the flow separation entirely.     

The Multi-point Optimization of Shock Control Bump with Constant-Lift Constraint Enhanced with Suction and Blowing for a Supercritical Airfoil

Both shock control bump (SCB) and suction and blowing are flow control methods used to control the shock wave/boundary layer interaction (SWBLI) in order to reduce the resulting wave drag in transonic flows. A SCB uses a small local surface deformation to reduce the shock-wave strength, while suction decreases the boundary-layer thickness and blowing delays the flow separation. Here a multi-point optimization method under a constant-lift-coefficient constraint is used to find the optimum design of SCB and suction and blowing. These flow control methods are used separately or together on a RAE-2822 supercritical airfoil for a wide range of off-design transonic Mach numbers. The RANS flow equations are solved using Roe’s averages scheme and a gradient-based adjoint algorithm is used to find the optimum location and shape of all devices. It is shown that the simultaneous application of blowing and SCB (hybrid blowing/SCB) improves the average aerodynamic efficiency at off-design conditions by 18.2 % in comparison with the clean airfoil, while this increase is only 16.9 % for the hybrid suction/SCB. We have also studied the SWBLI and how the optimization algorithm makes the flow wave structure and interactions of the shock wave with the boundary layer favorable.     

Stability of Steady‐State Shear Jet Flows of an Ideal Fluid with a Free Boundary in an Azimuthal Magnetic Field Against Small Long‐Wave Perturbations

The problem of the linear stability of steadystate axisymmetric shear jet flows of a perfectly conducting inviscid incompressible fluid with a free surface in an azimuthal magnetic field is studied. The necessary and sufficient condition for the stability of these flows against small axisymmetric longwave perturbations of special form is obtained by the direct Lyapunov method. It is shown that if this stability condition is not satisfied, the steadystate flows considered are unstable to arbitrary small axisymmetric longwave perturbations. A priori exponential estimates are obtained for the growth of small perturbations. Examples are given of the steadystate flows and small perturbations imposed on them which evolve in time according to the estimates obtained.