Schardin问题的数值研究

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

Schardin问题的数值研究

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

Application of the wavenumber jump condition to the normal and oblique interaction of a plane acoustic wave and a plane shock

A kinematic approach is considered whereby the wavenumber jump conditions in conjunction with the appropriate dispersion relations is applied to the investigation of the normal and oblique interaction of a plane acoustic wave with a plane shock wave. For the normal interaction of an acoustic wave with a stationary plane shock a logarithmic shift in the wave spectra is obtained. For the normal interaction with a moving shock front it is shown that for shock Mach numbers above a critical value, the frequency of the transmitted wave becomes negative. This results in the fact that the crests of the transmitted signal arrive at a fixed observer in a reverse order to their generation. Finally, the oblique interaction of an acoustic wave with a stationary shock is considered. The “Snell's Law” for the transmitted wave is derived and two special angles of incidence are identified. The first is a no-refraction angle: i.e., the transmitted wave angle is the same as the incident wave angle. The second is a critical angle such that for incident angles greater than this critical angle there is no transmitted wave. A necessary and sufficient condition for the existence of a transmitted wave is derived in terms of the speed of sound and Mach number of the fluid and the frequency and tangential wavenumber component of the incident wave.The dynamics aspects of the interaction concerning the determination of the frequency independent transmission coefficients and shock displacements are determined for the simple case of the normal interaction with a moving shock as an illustration.     

Differential characteristic of the flow behind a shock wave

Relationships between the derivatives on both sides of a discontinuity in a nonstationary shock wave moving with acceleration in a one-dimensional vortex flow of perfect gas are deduced. The problem of interaction between the shock wave and a weak discontinuity is solved based on these relationships.     

Resonance attenuation of ultrasonic waves in a superconductor with a moving vortex structure

Equations describing the interaction of ultrasonic waves with a moving vortex structure are derived. The addition to attenuation and the relative change in the velocity of longitudinal ultrasonic waves due to this interaction are calculated. It is found that when a longitudinal ultrasonic wave propagates along the direction of motion of the vortex structure and the velocity V of the structure is equal to half the velocity of the wave, then anomalous acoustic attenuation occurs and the contribution from the ultrasound-vortex interaction to the velocity of the ultrasonic wave vanishes. It is shown that if the vortex structure moves at a sufficiently high velocity, then (in contrast to the case of the structure at rest) a weakly damping collective mode propagating with velocity 2V arises in the structure. It is this mode that is responsible for anomalous attenuation of longitudinal ultrasonic waves.     

水下涡流场前向声散射特性分析

研究水下涡声散射特性,在目标探测和流场声成像领域具有重要意义。针对水下低马赫数涡流场前向声散射建立了数值计算方法,探究了其形态函数和指向性。首先,基于摄动声学理论给出了考虑流声耦合作用的涡声散射模型,采用时域有限差分结合完美匹配层构建了数值求解方法;随后,在算法验证的基础上,预报分析了高斯涡涡核尺寸在1~10 m,同时入射平面波无量纲波数在1~10范围内,涡流场强度对前向声散射特性的影响。结果表明,低马赫数下,声散射场具有对称性,且有明显的主瓣和指向性。其前向散射形态函数随入射波波数、涡核尺寸、涡流场强度增加而增大;主瓣方位角随波数增加而趋近入射波传播方向。      

Control of surface plasmon excitation via the scattering of light by a nanoparticle

We analyze the scattering of vortex pairs (the particular case of 2D dark solitons) by a single quantum vortex in a Bose–Einstein condensate with repulsive interaction between atoms. For this purpose, an asymptotic theory describing the dynamics of such 2D soliton-like formations in an arbitrary smoothly nonuniform flow of a ultracold Bose gas is developed. Disregarding the radiation loss associated with acoustic wave emission, we demonstrate that vortex–antivortex pairs can be put in correspondence with quasiparticles, and their behavior can be described by canonical Hamilton equations. For these equations, we determine the integrals of motion that can be used to classify various regimes of scattering of vortex pairs by a single quantum vortex. Theoretical constructions are confirmed by numerical calculations performed directly in terms of the Gross–Pitaevskii equation. We propose a method for estimating the radiation loss in a collision of a soliton-like formation with a phase singularity. It is shown by direct numerical simulation that under certain conditions, the interaction of vortex pairs with a core of a single quantum vortex is accompanied by quite intense acoustic wave emission; as a result, the conditions for applicability of the asymptotic theory developed here are violated. In particular, it is visually demonstrated by a specific example how radiation losses lead to a transformation of a vortex–antivortex pair into a vortex-free 2D dark soliton (i.e., to the annihilation of phase singularities).     

Experimental investigation of a strongly shocked gas bubble

A free-falling, spherical, soap-film bubble filled with argon is subjected to a planar M=2.88 shock in atmospheric nitrogen; vorticity is deposited on the surface of the bubble during shock interaction, and the Richtmyer-Meshkov instability ensues. The geometrical development of the shocked bubble is diagnosed with laser sheet imaging and a planar slice showing two cross sections of both the major vortex ring and a secondary vortex ring is revealed experimentally for the first time. Quantitative measurements of the experimental data include the vortex velocity defect, and subsequent circulation calculations, along with a new set of relevant length scales. The shock wave strength, leading to a post-shock compressible regime, allows the study of the instability development in a regime between low Mach number shock tube experiments and high Mach number laser driven experiments that has not been investigated previously.     

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

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

入射和反射激波诱导重气泡变形和失稳的三维数值研究

 采用大涡模拟方法,对入射激波及其反射激波诱导球形重气泡的变形失稳过程进行了三维数值模拟,利用已有实验验证了计算模型的可靠性,重点考察了反射激波与已经失稳的气泡界面的再次作用,讨论了涡环的形成及其三维失稳的过程。研究结果显示：入射和反射激波与球形重气泡作用产生斜压效应,会在流场中产生旋转方向截然相反的多个涡环;反射激波诱导的涡环具有较小的强度,故更加容易失稳,甚至能完全形成具有流向涡量的复杂小尺度涡结构。     

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.     

Visualization of shock-vortex interaction radiating acoustic waves

Unsteady compressible flow fields past a wedge and a cone, evolved by propagation and interaction of shock waves, slip lines, and vortices, are studied by shadowgraphs and holographic interferograms taken during the shock tube experiment. The supplementary numerical calculation also presented time-accurate solution of the shock wave physics which was essential to recognize the similarity and dissimilarity between the wedge and the conical flows. The decelerated shock detained by the vortex interacts with the small vortexlets along the slip layer, producing diverging acoustics: this phenomenon is more distinct in the case of wedge flow for a given shock Mach number. The decelerated shock penetrated through the vortex core constitutes a transmitted shock, which eventually merges with the diaphragm shock that bridges the vortex pair/vortex ring. This phenomenon became remarkably salient in the case of conical flow.     

Nonlinear sound waves upon collapse of a vortex dipole

Acoustic emission by a vortex dipole decreases the dipole energy and changes the vortex velocity distribution. A relative shift of the dipole component occurs. Earlier, such a process of acoustic instability of the vortex was considered as applied to the case of weak acoustic emission described by the wave equation. Features of acoustic instability upon emission of a strongly nonlinear sound by a vortex are considered in this work.     

On whistling of pipes with a corrugated segment: Experiment and theory

Corrugated pipes are commonly used because of their local rigidity combined with global flexibility. The flow through such a pipe can induce strong whistling tones, which is an environmental nuisance and can be a threat to the mechanical integrity of the system. This paper considers the use of a composite pipe: a shorter corrugated pipe segment embedded between smooth pipe segments. Such a pipe retains some flexibility, while the acoustical damping in the smooth pipe reduces whistling tones. Whistling is the result of coherent vortex shedding at the cavities in the wall. This vortex shedding is synchronized by longitudinal acoustic waves traveling along the pipe. The acoustic waves trigger the vortex shedding, which reinforces the acoustic field for a critical range of the Strouhal number values. A linear theory for plane wave propagation and the sound production is proposed, which allows a prediction of the Mach number at the threshold of whistling in such pipes. A semi-empirical approach is chosen to determine the sound source in this model. This source corresponds to a fluctuating force acting on the fluid as a consequence of the vortex shedding. The functional form of the Strouhal number dependency of the dimensionless sound source amplitude is based on numerical simulations. The magnitude of the source and the Strouhal number range in which it can drive whistling are determined by matching the model to results for a specific corrugated pipe segment length. This semi-empirical source model is then applied to composite pipes with different corrugated segment lengths. In addition, the effect of inlet acoustical convective losses due to flow separation is considered. The Mach number at the threshold of whistling is predicted within a factor 2.     

Experimental investigation of primary and secondary features in high-mach-number shock-bubble interaction

Experiments to study the compression and unstable evolution of an isolated soap-film bubble containing helium, subjected to a strong planar shock wave (M=2.95) in ambient nitrogen, have been performed in a vertical shock tube of square internal cross section using planar laser diagnostics. The early phase of the interaction process is dominated by the formation of a primary vortex ring due to the baroclinic source of vorticity deposited during the shock-bubble interaction, and the mass transfer from the body of the bubble to the vortex ring. The late time (long after shock interaction) study reveals the presence of a secondary baroclinic source of vorticity at high Mach number which is responsible for the formation of counterrotating secondary and tertiary vortex rings and the subsequent larger rate of elongation of the bubble.     

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.     

Self-sustained oscillations of a shock wave interacting with a boundary layer on a supercritical airfoil

A theory is proposed of the self-sustaining oscillations of a weak shock on an airfoil in steady, transonic flow. The interaction of the shock with the boundary layer on the airfoil produces displacement thickness fluctuations which convect downstream and generate sound by interaction with the trailing edge. A feedback loop is established when this sound impinges on the shock wave, resulting in the production of further fluctuations in the displacement thickness. The details are worked out for an idealized mean boundary layer velocity profile, but strong support for the basic hypotheses of the theory is provided by a comparison with recent experiments involving the generation of acoustic “tone bursts” by a supercritical airfoil section.     

A new approach of a limiting process for multi-dimensional flows

An enhanced Multi-dimensional Limiting Process (e-MLP) is developed for the accurate and efficient computation of multi-dimensional flows based on the Multi-dimensional Limiting Process (MLP). The new limiting process includes a distinguishing step and an enhanced multi-dimensional limiting process. First, the distinguishing step, which is independent of high order interpolation and flux evaluation, is newly introduced. It performs a multi-dimensional search of a discontinuity. The entire computational domain is then divided into continuous, linear discontinuous and nonlinear discontinuous regions. Second, limiting functions are appropriately switched according to the type of each region; in a continuous region, there are no limiting processes and only higher order accurate interpolation is performed. In linear discontinuous and nonlinear discontinuous regions, TVD criterion and MLP limiter are respectively used to remove oscillation. Hence, e-MLP has a number of advantages, as it incorporates useful features of MLP limiter such as multi-dimensional monotonicity and straightforward extensionality to higher order interpolation. It is applicable to local extrema and prevents excessive damping in a linear discontinuous region through application of appropriate limiting criteria. It is efficient because a limiting function is applied only to a discontinuous region. In addition, it is robust against shock instability due to the strict distinction of the computational domain and the use of regional information in a flux scheme as well as a high order interpolation scheme. This new limiting process was applied to numerous test cases including one-dimensional shock/sine wave interaction problem, oblique stationary contact discontinuity, isentropic vortex flow, high speed flow in a blunt body, planar shock/density bubble interaction, shock wave/vortex interaction and, particularly, magnetohydrodynamic (MHD) cloud-shock interaction problems. Through these tests, it was verified that e-MLP substantially enhances the accuracy and efficiency with both continuous and discontinuous multi-dimensional flows.     

Stability and three-wave interaction of disturbances in supersonic boundary layer with mass exchange on the wall

The interaction of disturbances in a boundary layer of the compressible gas is considered in the linear and nonlinear approximation (the weakly nonlinear theory of stability) in the presence of mass exchange (gas blowing or suction) on the surface. The regimes of moderate (the Mach number M = 2) and high (M = 5.35) supersonic velocities of the flow are considered. The suction from the surface is shown to lead to a considerable variation of the linear evolution of disturbances: the vortex disturbances of the first mode and the acoustic disturbances of the second mode are stabilized, the rate of variation is determined by suction intensity. The nonlinear interactions in three-wave systems between the vortex waves in asymmetric triplets at M = 2 and between the waves of different nature (acoustic and vortex waves) ?? in the symmetric ones at M = 5.35 are considered. The planar acoustic wave is the excitation wave in the latter, which excites the three-dimensional subharmonic components of the vortex nature. It is shown that one can delay considerably the transition region with the aid of suction, thereby one can reduce the skin-friction drag. In the gas blowing regime, strong deformations of the mean fields of boundary layers occur, which lead to the destabilization of the vortex and acoustic waves in the linear region and activate the nonlinear processes in transition region. One can expect that this will lead to the acceleration of tripping in supersonic flow.