On nonequilibrium of pyrolysis process in the manufacture of ethene

Nonequilibrium process for cracking ethane and n-buthane in the manufacture of ethene has been analytically and numerically investigated in a Heavi-side function temperature field and through a normal shock wave. The results demonstrate that, while the reaction temperature increases, the maximum value of ethene yield is increased, and the optimal reaction duration is sharply shortened. For the identical initial reaction temperature, the maximum value of ethene yield through a stationary normal shock wave is less than that in a Heavi-side function temperature field. However, the ethene consumption after the maximum value in the former case is less than that in the latter. Higher ethene yield will be obtained by using the gasdynamic heating method than by using the current methods.

激波在带扩容室的变截面管中传播和消波效果的数值分析

在利用实验结果验证程序的基础上，对激波在具有复杂结构管道内的传播特性进行了数值模拟，比较了不同入射激波条件和特征尺寸对扩容室消波效应的影响，并考虑了二维(轴对称、平面)模型和三维模型的影响。计算结果表明，同一尺寸的扩容室对不同马赫数人射激波的消波效果是相近的，扩容室特征尺寸不同则消波效果差别很大。

The Mach reflection of a detonation based on soot track measurements

This paper presents a series of soot tracks formed by gaseous detonation waves diffracting around wedges with different wedge angles. These cellular structure patterns describe the Mach-reflection processes of a detonation and reveal some unique characteristics. They can be used to analyze the relationship between the trajectory angle of the triple point, wedge angle, and initial pressure in Mach reflection. Compared to the Mach-reflected one-dimensional shock wave in nonreactive air, all these unique characteristics for a Mach-reflected detonation should be attributed to the transverse-wave structure of the detonation front; meanwhile, the precursor shock wave and transverse wave influence the Mach-reflected detonation, respectively. The experimental results support the recently published numerical simulation of this complex phenomenon.

用环形激波聚焦实现爆轰波直接起爆的数值模拟

利用基元反应模型和有限体积法对环形激波在可燃气体中聚焦实现爆轰波直接起爆进行了数值模拟.研究结果表明,标准状态下的氢气-空气混合气体在马赫数为3.1以上的环形激波聚焦产生的高温高压区作用下会诱发可燃气体的直接起爆形成爆轰波,爆轰波与激波和接触间断相互作用产生了复杂的波系结构;爆轰波爆点位置在对称轴上并不是固定的点,而是随着初始激波马赫数的变化而发生移动;可燃气体初始温度和压力对起爆临界马赫数都有影响,但是初始温度的影响大得多.

NS方程激波计算的摄动有限差分方法

摄动有限差分（PFD）方法从一阶迎风差分格式出发，将差分系数展开为网格步长的幂级数，通过提高修正微分方程的逼近精度来获得更高精度的差分格式。由于格式基于一阶迎风格式，因此具有迎风效应、网格节点少等特点。本文首先通过对Burgers方程的摄动差分格式的推导，将摄动有限差分格式引入时间相关法的计算，并构造了守恒形式的摄动有限差分格式，然后推广到一维Navier-Stokes方程组的计算。数值比较研究表明：本文构造的NS方程摄动有限差分格式具有比一阶迎风较高的精度和分辨率，而且保持了一阶迎风格式的无振荡性质。

Shocked Flows Induced by Supersonic Projectiles Moving in Tubes

A numerical study on shocked flows induced by a supersonic projectile moving in tubes is described in this paper. The dispersion-controlled scheme was adopted to solve the Euler equations implemented with moving boundary conditions. Four test cases were carried out in the present study: the first two cases are for validation of numerical algorithms and verification of moving boundary conditions, and the last two cases are for investigation into wave dynamic processes induced by the projectile moving at Mach numbers of M-p = 2.0 and 2.4, respectively, in a short time duration after the projectile was released from a shock tube into a big chamber. It was found that complex shock phenomena exist in the shocked flow, resulting from shock-wave/projectile interaction, shock-wave focusing, shock-wave reflection and shock-wave/contact-surface interactions, from which turbulence and vortices may be generated. This is a fundamental study on complex shock phenomena, and is also a useful investigation for understanding on shocked flows in the ram accelerator that may provide a highly efficient facility for launching hypersonic projectiles.

Direct Numerical Simulation Of Three-Dimensional Richtmyer-Meshkov Instability

Direct numerical simulation (DNS) is used to study flow characteristics after interaction of a planar shock with a spherical media interface in each side of which the density is different. This interfacial instability is known as the Richtmyer-Meshkov (R-M) instability. The compressible Navier-Stoke equations are discretized with group velocity control (GVC) modified fourth order accurate compact difference scheme. Three-dimensional numerical simulations are performed for R-M instability installed passing a shock through a spherical interface. Based on numerical results the characteristics of 3D R-M instability are analysed. The evaluation for distortion of the interface, the deformation of the incident shock wave and effects of refraction, reflection and diffraction are presented. The effects of the interfacial instability on produced vorticity and mixing is discussed.

Effects of Mach number on turbulent separation behaviours induced by blunt fin

An experimental study of the interaction between shock wave and turbulent boundary layer induced by blunt fin has been carried out at M-infinity = 7.8 using oil flow visualization and simultaneous measurements of fluctuating wall pressure and heat transfer. This paper presents the effects of Mach number on turbulent separation behaviours induced by blunt fin.

Novel experimental phenomena of fine-particle fluidized beds

The dynamic response of bed height and concentration waves in liquid-solid fluidized beds to a step change in the fluidization velocity is considered. We experimentally study the liquid-solid fluidized beds, spherical beadings, with sizes ranging from 230 to 270 mesh and the inner diameter of columns made from glass is 2.4 mm. Experimental results find that under certain conditions, fine particles with large Richardson-Zaki exponent n display different dynamic behavior from usual particles with smaller n during expansion and collapse of the fluidized state. (c) 2007 Elsevier Inc. All rights reserved.

Numerical-simulation of low supersonic-flow around a sphere

The controlled equations defined in a physical plane are changed into those in a computational plane with coordinate transformations suitable for different Mach number M(infinity). The computational area is limited in the body surface and in the vicinities of detached shock wave and sonic line. Thus the area can be greatly cut down when the shock wave moves away from the body surface as M(infinity) --> 1. Highly accurate, total variation diminishing (TVD) finite-difference schemes are used to calculate the low supersonic flowfield around a sphere. The stand-off distance, location of sonic line, etc. are well comparable with experimental data. The long pending problem concerning a flow passing a sphere at 1.3 greater-than-or-equal-to M(infinity) > 1 has been settled, and some new results on M(infinity) = 1.05 have been presented.

Upwind compact scheme with dispersion control for aerodynamic equations

A compact upwind scheme with dispersion control is developed using a dissipation analogy of the dispersion term. The term is important in reducing the unphysical fluctuations in numerical solutions. The scheme depends on three free parameters that may be used to regulate the size of dissipation as well as the size and direction of dispersion. A coefficient to coordinate the dispersion is given. The scheme has high accuracy, the method is simple, and the amount of computation is small. It also has a good capability of capturing shock waves. Numerical experiments are carried out with two-dimensional shock wave reflections and the results are very satisfactory.

The effect of variable currents on internal solitary waves

The effect of variable currents on internal solitary waves is described within the context of a variable coefficient Korteweg-de Vries (KdV) equation, and the approximate slowly varying, solitary-wave solution of this equation. The general theory which leads to the variable coefficient KdV equation is described; a derivation for the special case when the solitary wave and the current are aligned in the same direction is given in the Appendix. Using further simplifications and approximations, a number of analytical expressions are obtained for the variation in the solitary wave amplitude resulting from variable shear in the basic current or from when the basic current is a depth-independent flow which is a simple representation of a geostrophic current, tidal flow or inertial wave.

Phase relaxation model for 2-phase flows

A phase relaxation model (PRM) for 2-phase flows is presented in this paper on the basis of three principal assumptions. The basic equations for PRM arc derived from the Boltzmann equations for gas-partlcle mixture, The general characteristics and solving process of the PRM's basic equations are also presented and discussed. Many terms in the PRM's basic equations contain a factor ε= ρgρp/ρg+ρp2 which is an intrinsic small parameter for 2-phase mixture, with ρg and ρp being respectively the densities of gas and particle phases.This makes it possible to simplify the computation of the PRM's basic equations. The model is applied to for example, studying file steady propagation of shock waves in gas-particle mixture. The analysis shows that with an increase of shock wave strength the relaxation process behind a gasdynamics shock front becomes a kind of dynamics relaxation instead of the standard exponential relaxation process. A method of determining experimentally the velocity and tem...更多perature relaxation rates (or times) of gas-particle flows is suggested and analyzed.

可燃气体中激波聚焦的点火特性

数值模拟了二维平面激波从抛物面上反射在可燃气体中聚焦的过程，研究了形成爆轰波的点火特性,对理想化学当量比氢气／空气混合气体，在初始压强20kPa的条件下，马赫数2．6-2．8的激波聚焦能产生两个点火区：第1个点火区是反射激波会聚引起的，第2个点火区是由入射激波在抛物面上发生马赫反射引起的．这种条件下流场中会出现爆燃转爆轰，起爆点分别分布在管道壁面、抛物反射面和第2点火区附近．起爆机理分别为激波管道壁面反射、点火诱导激波的抛物面反射和点火诱导的激波与第2点火区产生的爆燃波的相互作用,不同的点火和起爆过程导致了不同的流场波系结构，同时影响了爆轰波传播的波动力学过程．

可燃气体中激波与障碍物作用在下游形成爆轰波的数值研究

对平面激波和单个矩形障碍物作用的过程进行了数值模拟,研究了反射产生的上行爆轰波在下游可燃气体中形成爆轰波的过程。数值结果表明,下游爆轰波形成过程主要有2种模式:爆轰波直接绕射和绕射波在上壁面反射,这和已有的实验结果是一致的。通过研究下游爆轰波的形成过程受入射激波马赫数、混合气体的压力和管道尺度的影响,分析了上游爆轰波向下游传播的波动力学过程,讨论了2种形成过程的作用规律和控制因素,阐明了下游爆轰波的形成规律。