射流推力矢量喷管中粒子运动轨迹的数值模拟

用颗粒轨道模型对基于激波控制的二维收缩-扩张(2-DCD)喷管中的气固两相流动进行了数值模拟,研究了1～60 μm不同直径颗粒的运动轨迹.研究表明,颗粒尺寸越小,粒子的随流性越好;颗粒尺寸越大,主流中的粒子与收缩段壁面碰撞越剧烈,无粒子区越大,二次流中的粒子对下壁面的碰撞越剧烈.     

冲压发动机燃烧引起激波运动过程一维分析

分析了冲压发动机喷油燃烧引起内流道内正激波运动的机理,采用一维激波捕捉方法,建立了燃油喷入对正激波运动位置影响的一维仿真模型。通过仿真发现:喷入燃油并逐步增大燃油-空气当量比时,正激波逐步向上游运动;燃油-空气当量比越大,正激波越接近进气道喉道;当燃油-空气当量比增大到一定程度时,正激波距离进气道喉道最近,但并未越过喉道;进一步增大燃油-空气当量比,正激波开始向下游回退进一步分析发现:冲压发动机流道及燃烧组织匹配设计直接影响到正激波在流道内的运动位置,需要在设计中格外重视。燃油-空气当量比与激波位置的关系分析可为冲压发动机设计提供一定的理论参考。     

不同密度弹丸对水冰的超高速撞击成坑实验

为了研究冰冻天体表面撞击坑的形成与演化,开展了水冰的超高速撞击成坑实验。使用二级轻气炮发射1.0 mm直径的球形弹丸,以3 km/s、5 km/s和7 km/s速度对圆柱状冰块进行撞击。弹丸材料包括聚碳酸酯和不锈钢两种,冰块温度为253 K。实验观察到了不同弹丸和不同速度条件下,冰块中撞击坑的形貌特征。对撞击坑直径、深度和剖面形状进行了测量,并与文献中铝弹丸对水冰的撞击坑进行了比较分析。获得了水冰撞击坑特征随撞击参数的变化规律,结果表明：撞击坑直径和深度的主导机制不同,坑深主要由弹丸侵彻作用形成,而坑径主要由冰块的剥落所致;坑深比坑径具有更强的对于弹丸密度的依赖性,高密度弹丸撞击坑直径具有比低密度弹丸更强的对于撞击速度的依赖性;撞击坑体积与撞击能量成正比,高密度弹丸形成的撞击坑直径表现出“能量缩比”行为,而低密度弹丸形成的撞击坑直径表现出“动量缩比”行为。     

高超声速混合流动的粒子模拟耦合算法

通过采用基于三维非结构网格的子网格技术,对适用于高超声速混合流动的自适应时间步长粒子模拟耦合算法(Improved Hybrid Particle Simulation Method,IHPSM)进行了改进,在保证算法计算效率的同时降低数值误差。通过对三维双曲钝锥外形的数值仿真及与DSMC(Direct Simulation Monte Carlo)计算结果的对比分析,证明改进的IHPSM算法具有较高的计算精度,且能较大幅度提高计算效率。基于改进的IHPSM算法,文中针对双曲钝锥外形进行了稀薄气体效应和飞行马赫数对高超声速流动影响规律的研究。结果表明,气体稀薄程度的增加会减缓双曲钝锥前端流场宏观物理量的变化梯度,使流场激波结构变厚;来流马赫数的增大会使激波明显增强,但对激波厚度与结构的影响较小。     

斜激波串与上游激波干扰特性试验研究

为研究斜激波串在背压条件下前移与上游激波相互干扰的流场结构和运动规律,在来流为马赫数 2.7 的直管道内设计一种等宽度斜楔,采用动态压力测量、高速纹影和粒子图像测速(PIV)技术等手段进行了试验。研究结果表明：内置斜楔在管道内产生入射激波、分离激波、膨胀波、再附激波和激波诱导分离等复杂上游激波流场,在分离区附近形成有顺压梯度和逆压梯度的区域。当增大下游压比时,斜激波串逐渐向上游激波流场移动;经过斜楔产生的分离区时,斜激波串的移动速度急剧提升,同时出现非对称分离偏转方向的切换。对比了三种长度尺寸的等楔角斜楔所产生的上游激波流场的差异性,发现在相同的斜楔前缘起始点和楔角时,随着斜楔长度的增加,上游激波流场中激波诱导的分离尺度逐渐变大。     

Invariance of Charged Particle Time Profiles at Late Stages of SCR Events from the Data of Multisatellite Observations

Based on the data obtained in simultaneous measurements on the Helios-1 and Helios-2 spacecraft, the properties of the spatial and temporal invariance of spectra of protons with energies more than 4 MeV and electrons with energies more than 0.3 MeV are studied during the phase of intensity decay in events initiated by the following associations: a flare—coronal mass ejection—coronal and interplanetary shock waves. Emphasis is placed on research into the fluxes of electrons, whose possibility to be accelerated by shock waves is problematic. It is shown that the spectra and decay times of both protons and electrons are approximately invariant relative to the flare coordinates in a wide angular range. Depending on the western or eastern spacecraft position in relation to the front part of a shock wave, one observes a delay or advance of the onset of the invariant mode relative to the shock wave arrival, which increases with the growth of the corresponding angular distance.     

单喷管火箭自由喷流噪声数值模拟

利用非定常燃气流场数值模拟方法，研究超声速喷流核心区及其外围亚声速区燃气流动态分布和流动特性，在此基础上利用欧拉方程条件的伽辽金有限元方法以及FWH方法，实现并完成喷流噪声传播特性、辐射特性数值模拟。喷流噪声数值模拟结果显示：燃气流推进初期，强喷流噪声区域紧随燃气流前锋；燃气流场相对稳定后，强喷流噪声区域主要位于燃气流等能区末稍，一些小尺度试验的燃气流激波系附近也存在较强喷流噪声。这些强喷流噪声主要由燃气流前锋带动的大涡动态卷吸、燃气流强湍流脉动以及激波扰动引起。受数值模拟网格分辨率影响，当前仅能保持中低频段声压级数值模拟结果与实测结果总体接近。     

Cylindrical heterogeneous detonation waves

Further experimental studies of blast wave initiated cylindrical heterogeneous (liquid fuel drops, gas oxidizer) detonation waves are described. A pie-shaped shock tube, used for these studies, was altered in certain ways so as to improve the modeling of cylindrical waves. These modifications, along with some operational aspects, are briefly discussed. The breech of the facility, where the blast wave is generated by an explosive, became distorted with usage. Results are presented which show that lower detonation velocities are realized with the damaged breech (other conditions being the same). A photographic and pressure switch wave time of arrival study was made to ascertain the wave shape. Photographs are shown which show that the waves, blast as well as detonation, are close to cylindrical. However, in some cases there is appreciable distortion of the wave front by debris ahead of the wave. Presumably this debris comes from the blasting cap used to ignite the condensed explosive. A series of experiments was conducted using kerosene drops of 388 μm diameter dispersed in air through use of a large number of hypodermic needles. Radial fuel void regions were established by cutting off the fuel flow to a number of needles. Preliminary results relating to the effect of the size of the cloud gap on detonation velocity, quenching, and the initiator energy levels required for detonation are discussed.     

Perturbation solutions for variable energy blast waves

The trajectory of and the flow field behind blast waves with time varying energy input is determined. Freeman's (1968) Lagrangean coordinate formulation is modified to include both the geometric factor, α, for plane, cylindrical and spherical shocks and also non-integer values of β, the energy input parameter, in a single computational algorithm. Numerical problems associated with vanishing density at the inner mass boundary or “piston face” are then examined and solved. Second order perturbation solutions about the solution for an infinite strength shock are then obtained in Sakurai's (1965) inverse shock Mach number expansion parameter for 0 β < α + 1. Tables and graphs of significant numerical coefficients are presented for comparison to, and extension of, results of other authors. Graphs of typical shock trajectories and flow field density, pressure and velocity variations are also presented and discussed.     

On flows with vibrational relaxation behind an attached shock wave

The flow conditions are represented by an oblique shock attached to the vertex of a straight edged wedge, placed symmetrically with respect to a uniform supersonic flow in front, so that the flow after the shcok is along the wedge surface. The curvature of the attached shock wave, and the flow variable gradients at the wedge tip caused by non-equilibrium vibrational effects are determined for a steady two-dimensional flow of a vibrationally relaxing gas. A relation between the curvature of the shock and that of the stream line is obtained at the wedge-tip. The variation of entropy and an expression for the vorticity generated by the shock towards down-stream flow are determined. Calculations have been carried out for the shock curvature and flow variable gradients for a wide range of wedge angle and the up-stream flow Mach number.     

Numerical experiments for supersonic and hypersonic flows

This paper contains the results of numerical experiments carried out by means of a modified non-stationary “large particle” method developed for the investigation of properties of both steady and non-steady flows associated with complicated gasdynamic problems. In the solution obtained by numerical analysis one may observe the formation of shock waves, local supersonic zones, separation regions with recirculation, and similar effects. Such phenomena are all of significance to a great number of practical problems.     

Sur des conditions d''amorçage par onde de choc de la détonation du nitrométhane

Results concerning the thermodynamic and mechanical properties of nitromethane (N.M.) at pressures up to 11 GPa are presented. The pressure is generated either by an incident shock wave or by a shock wave reflected within a medium prepressurized by a first shock wave. Calculation of the temperature behind these shock waves, based on the Walsh-Christian model, calls, in particular, upon the knowledge of the N.M. shock polar relative to normal temperature and pressure conditions, but also of those corresponding to prepressurized states.

Taking advantage of the phenomena of N.M. electric polarization under shock, we determine, on the one hand, the relation between pressure and particle velocity and, on the other hand, the influence of pressure conditions on induction delays of the explosive.

According to calculation, for a same pressure level the N.M. temperature behind a single shock is higher than that obtained by two successive compressions.

Experimentally, we observe that N.M. compressed at 11 GPa by means of two shock waves does not detonate (during the observation time of about 0,5 μs), while in the case of a single shock wave of the same amplitude the induction delay is lower than 0.1 μs. These results show the important role of temperature, as opposed to that of pressure.     

超音速流在正规反射激波系后的参数优化特性

研究了超音速气流中正规反射波加一正激波结构下的参数优化特性，给出了求解这种优化结构的方法，通过计算得到了优化条件存在的区域和优化解的特性，它用可于指导超音速飞行器的进气道及相关装置的设计。     

Numerical validation and parametric investigation on the cold flow field of a typical cavity-based scramjet combustor

The three-dimensional coupled implicit Reynolds Averaged Navier–Stokes (RANS) equations and the two equation standard k–ε turbulence model has been employed to numerically simulate the cold flow field in a typical cavity-based scramjet combustor. The numerical results show reasonable agreement with the schlieren photograph and the pressure distribution available in the open literature. The pressure distribution after the first pressure rise is under-predicted. There are five shock waves existing in the cold flow field of the referenced combustor. The first and second pressure rises on the upper wall of the combustor are predicted accurately with the medium grid. The other three shock waves occur in the core flow of the combustor. The location of the pressure rise due to these three shock waves could not be predicted accurately due to the presence of recirculation zone downstream of the small step. Further, the effect of length-to-depth ratio of the cavity and the back pressure on the wave structure in the combustor has been investigated. The obtained results show that there is an optimal length-to-depth ratio for the cavity to restrict the movement of the shock wave train in the flow field of the scramjet combustor. The low velocity region in the cavity affects the downstream flow field for low back pressure. The intensity of the shock wave generated at the exit of the isolator depends on the back pressure at the exit of the combustor and this in turn affects the pressure distribution on the upper wall of the combustor.     

Shock waves in water foams

Propagation of shock waves in tubes filled with water foams is studied using pressure gauges. Low amplitude shock waves consist of a precursor which propagates at a velocity slightly less than the acoustic velocity in the gas, and of a main compression wave which propagates slower than the precursor. Stronger shock waves have a single front. Maximum pressure rise in the incident and reflected shock waves cannot be calculated using one-dimensional conservation equations at the shock front. It is suggested that the flow of the liquid in foam cells has to be taken into account in order to predict the behavior of shock waves in foams. The nature of the gas which fills the cells is shown to have a strong effect on the quenching of blast waves in foams.     

Growth and decay of weak discontinuities in radiative gasdynamics

Growth and decay properties of weak discontinuities headed by wave fronts of arbitrary shape in three dimensions are investigated in a thermally radiating inviscid gas flow. The effects of radiative transfer are treated by the use of a general differential approximation for a grey gas of arbitrary opacity including effects of radiative flux, pressure and energy density. The transport equations representing the rate of change of discontinuities in the normal derivatives of the flow variables are obtained, and it is found that the nonlinearity in the governing equations does not contribute anything to the radiation induced waves. In contrast to the radiation induced waves, the nonlinearity in the governing equations plays an important role in the interplay of damping and steepening tendencies of a modified gasdynamic wave. An explicit criterion for the growth and decay of a modified gasdynamic wave along bicharacteristics curve in the characteristic manifold of the governing differential equations is given and the special reference is made of diverging and converging waves. It is shown that there is a special case of a compressive converging wave for which the stabilizing influences of thermal radiation and the wave front curvature are not strong enough to overcome the tendency of the wave to grow into a shock.     

High speed gas flows in explosions of cavitated explosives

A possibility of attaining steady flow of detonation products with specific energy much larger than the specific chemical energy of explosive is demonstrated in the case when a cylindrical charge of explosive is fitted with an evacuated cavity. Simple estimates and results of numerical analysis of the process are presented. Steady process may be considered to occur under the following assumptions: (1) effects arising due to jet interaction with cavity walls are negligible; (2) the detonation process is steady. In the case of limited explosive lengths these assumptions have been shown to be correct.When the cavity is filled with gas or liquid, a variety of steady and non-steady flow regimes is possible, depending on the density of the filling medium. One well-known case is that of flow with irregular reflection of shock waves at the cavity axis accompanied with the formation of Mach intersections. Another interesting flow regime is observed to occur in the case of low density filling medium (liquid hydrogen, for example). In this case the filling medium is driven by a “detonation piston” at constant velocity, equal to the velocity of detonation, forming a uniform growing column of hot shock-compressed matter, specific energy of which exceeds by one order of magnitude the specific energy of the explosive. Obviously, the walls of the vessel containing hydrogen must be able to withstand radial loads for a sufficiently long time (20 μ sec).The relative merits of these methods in comparison to others in high speed gas-dynamics is discussed.     

Rankine涡作用下激波变形和分叉的数值研究

数值模拟了激波和Rankine涡相互作用，重点研究了波涡相互作用引起的激波结构的变化过程。应用非结构化适应网格下的二阶精度Godunov型的PLM格式来求解Elder方程，对一平面激波与不同涡强Rankine涡相互作用下瞬态激波的运动、发展进行了研究。漩涡与激波间相对强度的大小对激波结构的影响明显，弱相互作用诱导激波变形，强相互作用诱导激波变形并产生分叉。计算结果表明，应用非结构化适应网格的PLM格式所得到的数值解较好地反映了瞬态激波结构的有关特征和信息。     

火箭垂直回收中发动机布局与喷流壁面效应影响研究

针对重复使用火箭垂直着陆过程的喷流流场问题开展研究,利用计算流体力学(Computational Fluid Dynamics, CFD)方法研究了壁面效应和发动机布局对超声速喷流的影响。研究表明,着陆距离(L)在2.24D～11.2D(D为喷管出口的直径)的范围内,地面效应对喷管出口中心处的温度分布影响较小;在当前计算条件下,当L<2.24D时,超声速喷流撞击地面会形成强烈的激波,随着离地高度的降低,该激波位置往喉部方向移动,由于壁面效应,喷管内部形成斜激波,导致中心喷管壁面处的温度升高;中心喷管相对外侧喷管往外突出增大了壁面流动速度,导致外侧喷管出口的温度降低;研究还表明子级火箭底部端面的喷管数量增加后,会导致喷管的温度升高。研究结果将为火箭发射及回收方案选取提供参考。     

喷管尺寸对电弧风洞流场和热环境的影响分析

为获得电弧风洞喷管尺寸对试验流场以及模型表面热流的影响规律,针对某特定模拟参数试验状态,采用高焓流动数值模拟方法对不同尺寸锥形喷管下的球柱校核模型试验流场进行了模拟和比较分析。研究发现,在模拟气流焓值和模型驻点热流的条件下,采用出口尺寸小的喷管所需电弧加热功率更低,同时单位流向截面上气流能量转化为模型驻点气动热的比例更低。不同喷管出口尺寸下,试验流场喷管出口区域热力学非平衡程度、波后氧原子质量分数、模型驻点区域压力以及表面传导热流和扩散热流占比都比较接近,但相较飞行状态存在明显差异;不同喷管出口尺寸下来流速度、激波脱体距离以及驻点线上平动温度之间的差异明显,喷管出口尺寸越大,其与飞行状态越接近。