"5 12"汶川地震震后龙门山山前地表塌陷成因探讨

应用NND有限差分格式求解轴对称可压缩N-S方程,研究了不同驻室与 环境压力比条件下欠膨胀超声速射流近场的失稳特性. 计算结果表明欠膨胀超声速射流的失 稳机制根据射流激波结构的特征可分为3种失稳模式：具有规则反射激波结构和单一剪切层 特征的射流不稳定性;带有马赫反射激波结构和双剪切层特征的射流不稳定性;具有弯曲马 赫杆和高度欠膨胀射流的不稳定性. 对于欠膨胀超声速射流,沿射流方向重复出现拟周期性 的射流激波结构是射流稳定发展的特征,这种射流激波结构的消失是射流开始失稳的标志.     

欠膨胀超声速射流不稳定性机理的数值研究

应用NND有限差分格式求解轴对称可压缩N-S方程,研究了不同驻室与环境压力比条件下欠膨胀超声速射流近场的失稳特性.计算结果表明欠膨胀超声速射流的失稳机制根据射流激波结构的特征可分为3种失稳模式:具有规则反射激波结构和单一剪切层特征的射流不稳定性;带有马赫反射激波结构和双剪切层特征的射流不稳定性;具有弯曲马赫杆和高度欠膨胀射流的不稳定性.对于欠膨胀超声速射流,沿射流方向重复出现拟周期性的射流激波结构是射流稳定发展的特征,这种射流激波结构的消失是射流开始失稳的标志.     

欠膨胀高速射流的压力均一化过程

发展了分析欠膨胀高速湍流射流近场特性的数值,计算了不同背压比条件下声速喷管出口附近膨胀区的激波胞格结构,揭示了压力一化过程的的机制。     

超声速中等欠膨胀冲击射流的涡声关系研究

使用大涡模拟方法对冲击面为平面的超声速中等欠膨胀冲击射流进行了数值模拟.利用三阶迎风和四阶对称紧致格式对无量纲化轴对称可压缩滤波N-S方程进行空间离散,时间上推进采用的是三阶精度的TVD型Rugge-Kutta法. 通过与经典的冲击射流实验比较,证实了程序的可靠性. 数值模拟得到了流场中不同尺度的涡结构和激波结构,观察到了上行声波和反射波以及流场中不同位置的声源,分析了冲击区剪切层附近区域的压强和涡旋转强度变化的频率、冲击平板上的压强变化频率以及射流剪切层中不同位置的涡合并出现的频率,发现冲击区剪切层附近区域的压强和涡强度变化以及射流剪切层中的涡合并现象和离散频率的冲击单音有重要关联.      

加权基本无振荡格式研究进展

加权型基本无振荡WENO格式是近十年发展起来的一类高阶、高精度格式,它是在ENO格式的基础上采用加权思想构造的,对流场内的间断和细致结构具有较高的分辨率,适于求解包含激波、膨胀波以及接触间断等复杂结构的流场,目前已发展成为计算流体力学中一种重要的方法。本文针对加权型基本无振荡格式近年来的进展作一简要介绍。     

二维时-空守恒格式及其在激波传播计算中的应用

将改进后的一维时－空守恒格式推广到了二维情形,得到了一个新的一般形式的二维Ｅｕｌｅｒ方程时－空守恒格式,并用该格式对几个具有复杂波系的流场进行了数值模拟。结果表明,该格式保留了一维格式通用性好、结构简单的优点,其计算结果精度高,对激波等间断具有很强的分辨率。     

基于模态分解的轴对称超声速射流啸声产生位置数值分析

不完全膨胀超声速射流的势核中会产生准周期的激波栅格结构, 其与剪切层内拟序结构的相互作用会产生激波噪声. 啸声是主要向上游方向传播的、具有离散频率的高强度激波噪声, 其产生是受一种非线性的声反馈环机制驱动. 精确定位啸声的声源位置是定量理解啸声反馈环机制和发展准确的啸声预测模型的一个关键所在. 为了分析近场啸声, 本文采用高精度数值方法直接求解轴对称可压缩Navier-Stokes方程, 数值模拟了完全膨胀射流马赫数为1.10和1.15的圆形声速喷管欠膨胀超声速冷射流, 得到了A₁和A₂两种轴对称模态啸声. 通过傅里叶模态分解、本征模态分解和动态模态分解, 分析了射流时序压力场和速度场, 研究了啸声关联拟序流动结构的空间演化, 精确定位了轴对称模态啸声的声源位置. 研究表明: 啸声关联拟序流动结构存在饱和态区域, 啸声声波是在其饱和态区域产生并向外传播; 在本文所涉及的射流马赫数范围内, A₁和A₂两种轴对称模态啸声的有效声源位置分别是在第4和第3个激波栅格结构的尾缘.      

有限谱ENO格式及其应用

首先对王健平提出的有限谱法^[1-3]做了进一步的理论研究，发现了一些新的有限谱法的插值基函数组，并互将有限谱法应用在ENO格式中，构造了有限谱ENO计算格式，然后通过对一维Euler方程的几个经典的模型问题和二维湍流与弱激波相互作用问题的数值计算，并且与理论解或准精确解进行比较分析，从而表明了此格式对于激波和其他间断具有较高的分辨率，在激波附近基本上没有明显的数值振荡，而且对于流场中的细致结构也具有相当高的精度。     

A STUDY ON JET PHENOMENON OF R22 GAS CYLINDER UNDER THE IMPACT OF SHOCK WAVE

重物质气泡在平面激波作用下出现射流是一个非常有趣的物理现象, 它与气体属性、激波强度密切相关, 气泡内外波系结构复杂, 值得进行深入探讨. 论文采用水平集方法结合真实虚拟流体方法捕捉物质界面, 对R₂₂圆柱气泡在平面激波作用下的射流现象进行了数值模拟, 详细阐述了波系结构在气泡内外的演变过程, 总结了导致气泡射流现象形成的典型波系结构特征, 从波系结构角度解释了射流形成的原因.论文也对不同强度的激波冲击R₂₂圆柱气泡的情况进行了归纳, 发现激波强度对射流现象的出现以及射流结构的尺寸特征有重要影响, 得出了 R₂₂圆柱气泡能否出现射流结构的临界激波强度, 以期对实验结果起到指导与验证作用.     

无粘可压缩流动的改进高精度方法

为更准确捕捉复杂流场的流动细节,通过对WENO格式的光滑因子进行改进,发展了一种新的五阶WENO格式。对三阶ENO格式进行加权可以得到五阶WENO格式,但是不同的加权处理,WENO格式在极值处保持加权基本无振荡的效果不同,本文构造了二阶精度的局部光滑因子,及不含一阶二阶导数的高阶全局光滑因子,从而实现WENO格式在极值处有五阶精度。基于改进五阶WENO格式,对一维对流方程、一维和二维可压缩无粘问题进行算例验证,并与传统WENO-JS格式和WENO-Z格式进行比较。计算结果表明,改进五阶WENO格式有较高的精度和收敛速度,有较低的数值耗散,能有效捕捉间断、激波和涡等复杂流动。     

Shock wave configurations and flow structures in non-axisymmetric underexpanded sonic jets

An experimental and numerical study of underexpanded free sonic jet flows issuing from rectangular, elliptical and slot nozzles has been undertaken. Aspect ratios (AR) of 1, 2, and 4 are described at pressure ratios (exit plane pressure to ambient pressure), of 2 and 3. There is good qualitative agreement between the experimental observations and the numerical predictions. In the case of rectangular jets, a complex system of shock waves forming the incident shock system is identified. This shock wave system originates at the corners of the nozzle exits, and proceeds downstream. Mach reflections are found to occur on the incident shock wave surface as well as the presence of a Mach disk terminating the first jet cell. This Mach disk has the shape of a square, a hexagon, or an octagon depending on the nozzle shape. For slot and elliptical jets, the formation of the incident shock wave was not observed along the minor axis plane of the nozzle for AR > 2. The incident shock wave was observed to originate downstream of the nozzle exit in the major axis plane. This wave system undergoes a transition to Mach reflection as it propagates downstream of the nozzle exit. In all cases tested, the shape of the jet boundary is significantly distorted. In rectangular jets, the narrowing of the jet boundary along the diagonal axis of the nozzle exit is observed, and in the case of the elliptical and slot jets axis switching is noted.     

Influence of nozzle geometry on the near-field structure of a highly underexpanded sonic jet

Detailed near-field structures of highly underexpanded sonic free jets have been investigated with the help of computational fluid dynamics. Two-dimensional, axisymmetric Euler equations have been chosen to predict the underexpanded jets, and the third-order total variation diminishing finite-difference scheme has been applied to solve the system of governing equations numerically. Several different nozzles have been employed to investigate the influence of the nozzle geometry on the near-field structures of highly underexpanded sonic free jets. The results obtained show that the distance from the nozzle exit to the Mach disk is an increasing function of the jet–pressure ratio, which also significantly influences the shape of the jet boundary. The diameter of the Mach disk increases with the jet–pressure ratio, and it is further significantly influenced by the nozzle geometry, unlike the distance of the Mach disk from the nozzle exit. However, such a dependence on the nozzle geometry is no longer found when an effective-diameter concept is taken into account for the flow from a sharp-edged orifice. A good correlation in the diameters of the Mach disk is obtained, so that the near-field structure of highly underexpanded sonic free jets is a unique function of the pressure ratio, regardless of the nozzle geometry.     

Starting of an axisymmetric convergent-divergent nozzle in a hypersonic flow

The starting of an axisymmetric convergent-divergent nozzle, with the result that supersonic flow is formed within almost the entire channel, is modeled, as applied to the hypersonic aerodynamic setup of the Institute of Mechanics of Moscow State University. A successful starting is realized when the nozzle is thrown in a uniform supersonic air flow at a fairly high Mach number. The steady flow structure is studied. It is numerically shown that in the convergent section of the channel there arises an oblique shock wave whose interaction with the nozzle axis leads to the formation of a reflected shock and a curvilinear Mach disk with a region of unsteady subsonic flow in the vicinity of the throat. The mathematical model is based on the two-dimensional Euler equations for axisymmetric gas flows.     

Experimental Investigation of Nozzle Exit Reflector Effect on Supersonic Jet

The present study describes an experimental work to investigate the effect of a nozzle exit reflector on a supersonic jet that is discharged from a convergent–divergent nozzle with a design Mach number of 2.0. An annular reflector is installed at the nozzle exit and its diameter is varied. A high-quality spark schlieren optical system is used to visualize detailed jet structures with and without the reflector. Impact pressure measurement using a pitot probe is also carried out to quantify the reflector’s effect on the supersonic jet which is in the range from an over-expanded to a moderately under-expanded state. The results obtained show that for over-expanded jets, the reflector substantially increases the jet spreading rate and reduces the supersonic length of the jet, compared with moderately under-expanded jets. The reflector’s effect appears more significant in imperfectly expanded jets that have strong shock cell structures, but is negligible in correctly expanded jet.     

Pitot pressures of correctly-expanded and underexpanded free jets from axisymmetric supersonic nozzles

The structures of the axisymmetric free jets from supersonic nozzles with the exit Mach numbers of 1.5 and 2.0 are studied with special attention to the decay of the Pitot pressures downstream of the Mach disk. The Pitot pressure probe and schlieren method are used in the experiments to diagnose the flowfield. A TVD numerical method is also applied to the Euler equations, and the computed jet structures are compared with experiments. In the underexpanded jet, the experimentally obtained Pitot pressure near the jet centerline is found to substantially recover downstream of the Mach disk. By comparing the numerical computation, this phenomenon is thought to be caused by the turbulent momentum transfer to the central region from the region outside the slip line where the stagnation pressure loss is small.     

Numerical investigation of inviscid shock wave dynamics in an expansion tube

Highly complicated shock wave dynamics has been numerically calculated by solving the Euler equations for a circular shock tube suddenly expanded three times of the original tube diameter atx=0. Shock waves of different shock Mach number,M _s =1.5 and 2.0, have produced remarkably distinct blast jet structures. A planar shock wave took its final form after the blast by repeated Mach reflections of the blast wave: the first one at the wall and the second one at the central axis. The central Mach disc overtook and merged with the annular Mach stem before the planar shock wave was formed. In contrast to the blast wave which would propagate spherically in an open space, the present blast wave undergoes complex morphological transformation in the restricted flow passage, resulting in an unstable and oscillatory blast jet structure of highly rotational nature. The slipstream tube emanating from the shock tube exit corner decomposed into a chain of small vortex rings that interacted with the barrel shock of the jet, which caused periodic collapse of the jet structure. The finite volume-FCT formulation equipped with the time-dependenth-refinement adaptive unstructured triangular mesh technique in the present paper has contributed to resolution of the intricate physical discontinuities developing in the blast flow fields.     

An experimental study of the under-expanded swirling jet with secondary coaxial stream

The present study addresses experimental results for investigating the details of the near field flow characteristics produced in an under-expanded, dual, coaxial, swirling jet. The under-expanded swirling jet is discharged from a sonic inner nozzle. An outer annular nozzle produces co- and counter-swirling streams relative to the inner primary swirling jet. The interaction between both the outer annular swirling stream and inner under-expanded swirling jet is quantified by impact and static pressure measurements, and visualized by using the shadowgraph method. Experiments are performed for several different pressure ratios. The results show that the outer secondary co-swirling jet significantly changes the structure of the inner under-expanded swirling jet, such as the shock structures and the recirculation region generated at the jet axis. The effect of the outer secondary stream on the major structures of the inner primary swirling jet is strongly dependent on the pressure ratio of the inner swirling jet, regardless of the swirl direction of the outer stream.Received: 17 May 2004, Accepted: 27 September 2004, Published online: 26 November 2004[/PUBLISHED]H.D. Kim: Correspondence to     

Numerical study of spherical blast-wave propagation and reflection

The objective of this study is to understand the flow structures of weak and strong spherical blast waves either propagating in a free field or interacting with a flat plate. A 5th-order weighted essentially non-oscillatory scheme with a 4th-order Runge-Kutta method is employed to solve the compressible Euler/Navier-Stokes equations in a finite volume approach. The real-gas effects are taken into account when high temperature occurs. A shock-tube problem with the real-gas effect is first tested in order to verify the solver accuracy. Moreover, unsteady shock waves moving over a stationary wedge with various wedge angles, resulting in different types of shock wave reflections, are also tested. It is found that the computed results agreed well with the existing data. Second, the propagation of a weak spherical blast wave, created by rupture of a high-pressure isothermal sphere, in a free field is studied. It is found that there are three minor shock waves moving behind the main shock. Third, the problem of a strong blast wave interacting with a flat plate is investigated. The flow structures associated with single and double Mach reflections are reported in detail. It is found that there are at least three local high-pressure regions near the flat plate. Received 27 July 2000 / Accepted 25 January 2002 – Published online 17 June 2002     

Flow of a mixture of gas and solid particles in supersonic underexpanded jets

The presence of a second phase in a gas jet flowing out of a nozzle leads to considerable changes in the flow pattern [1–3]. Thus, as the particle concentration increases, the central jump in compression [shock wave] moves in the direction of the nozzle cutoff, while the Mach number on the axis of the jet in front of the forward jump diminishes. In this paper we shall consider the numerical solution of the problem of an axisymmetrical, two-phase, underexpanded jet flowing out of a straight nozzle into a submerged space. It is assumed that the distribution of the flow parameters is uniform over the jet cross section and that no thermal or dynamic retardation of the particles occurs.