旋转爆轰胞格结构的实验和数值研究

对爆轰波在环形圆管(预混气体为2H2/O2/Ar)内的传播分别进行了实验和数值研究。实验研究 采用烟迹板记录了环形圆管内爆轰波的胞格结构。数值计算利用二阶附加半隐的Runge-Kutta法和五阶 WENO格式分别离散欧拉方程的时间和空间导数项,采用基元反应简化模型描述化学反应过程,得到了旋转 爆轰的流场及数值胞格结构。实验和数值模拟结果表明:爆轰波在圆环管中传播时,由于圆环的内壁发散、外 壁收敛,圆环内侧爆轰强度小于外侧,胞格尺寸较大;内侧OH 的分布区域大于外侧,浓度较低。旋转爆轰的 这种性质,使爆轰波能以稳定的角速度绕轴旋转。     

弯管内爆轰波传播的流场显示和数值模拟

采用激光纹影系统拍摄了爆轰波在不同位置的流场照片. 用二阶附加半隐的龙格- 库塔法和五阶WENO格式 分别离散欧拉方程时间和空间导数项，用基元反应来描述爆轰化学反应过程，获得了压力、 温度、典型组元质量分数分布及数值胞格结构和爆轰波平均速度. 结果表明：受壁面稀疏波 和压缩波影响，爆轰波阵面发生畸变. 但由于弯管曲率半径较大，未出现爆轰波熄灭. 靠近 凹壁面的激波强度大于凸壁面侧，且凹壁面侧的反应区宽度较凸壁面侧要窄. 弯管出口处的 三波点数目较入口处减少，爆轰波衰减. 在出口直段，受扰动的爆轰波可恢复为自持爆轰波. 爆轰波流场、胞格结构、平均爆轰波速度的计算和实验结果定性一致.     

粘性及热传导对于爆轰波的影响

对在满足化学当量比的氢氧混合气体中传播的一维和二维连续旋转爆轰波进行了数值模拟,以此检验粘性和热传导对爆轰波发展和结构的影响。模拟分别基于NS和Euler控制方程,采用二步化学反应模型,对流项采用5阶MPWENO格式求解,时间方向采用3阶TVD Runge-Kutta法,粘性项采用中心差分格式进行离散。结果表明:粘性和热传导不会对爆轰波流场的基本流场结构产生影响;在具体数值上,粘性和热传导的影响值在爆轰波、斜激波、接触间断等速度或温度剧烈变化处相对较大,但总体上其影响量均比爆轰波流场的步进值小三个量级。因此,在没有壁面效应的一维爆轰和二维连续旋转爆轰波流场中,粘性和热传导项作为很小的扰动存在,对爆轰波的流场结构和数值大小基本不会产生影响。     

不稳定性对爆轰波楔面马赫反射的影响规律研究

采用一种两步化学反应模型对胞格爆轰波的楔面马赫反射过程进行了数值研究,从而澄清和解释胞格不稳定性对马赫反射发展模式和自相似性的影响。考虑到反应欧拉方程源项的刚性问题,本文采用附加RungeKutta方法耦合非刚性对流项和刚性反应源项,对流项的离散采用五阶精度的WENO格式。计算结果表明,对于稳定胞格爆轰波而言,其马赫反射过程本质上与ZND爆轰波的马赫反射是一致的,整体上不存在自相似性,胞格不稳定性只是造成了三波点轨迹线局部小振幅的波动。在楔面顶点附近,由于马赫杆是强过驱的,爆轰波的马赫反射过程是自相似的。在远场,爆轰波马赫反射的三波点轨迹线渐近的趋向于一条直线,说明重新获得了自相似性。对于不稳定的爆轰波,由于自身的不稳定性可以与马赫反射的强度相匹配,定义其三波点的轨迹是困难的,进行自相似性分析没有意义。     

直管内胞格爆轰的基元反应数值研究

基于基元反应和二维欧拉方程，对直管内胞格爆轰进行了数值模拟。采用5阶WENO(weighted essentially nonoscillatory scheme)求解对流项，采用2阶附加半隐的龙格-库塔法处理化学反应源相引起的刚性。获得了密度、压力、温度和典型组元质量分数流场及数值胞格结构等。结果表明:网格精度的差异明显影响胞格的规则性和爆轰的平衡模数，随着网格尺度的减小，胞格由不规则变为规则。预混气组成、初压、初温及管道宽度给定，三波点数收敛为确定值。足够强度的初始扰动可再现胞格爆轰，最终形成的自持胞格爆轰模数与初始扰动的形状、大小、位置均无关。沿胞格中心线，爆轰波速度变化范围为0.88DCJ~1.5DCJ，爆轰波平均速度与CJ爆轰速度仅偏差0.88％。峰值压力与初压之比为14~50。计算爆轰波平均速度、胞格宽长比与实验值基本一致，但计算胞格宽度比实验值略小。数值模拟加深了对横波的产生和发展、未反应气囊、爆轰胞格的二次起爆等胞格爆轰特性的认识。     

气相爆轰在T形管中传播新现象的实验研究

对2H2/O2/Ar系统爆轰波在T形管(截面为40mm×40mm)中传播现象进行了实验研究.用烟迹片记录了T形管中爆轰波的胞格结构,用压电传感器记录了分叉口附近指定点压力时间曲线,得到了爆轰波在分叉口附近的平均速度和胞格图案演变.结果表明:初压P0≥2.67kPa,在水平和垂直支管下游区域(距离分叉口约3.5—6倍方管截面边长),分叉口影响消失,爆轰波恢复稳定,且强度基本保持不变.在分叉口绕射过程中,爆轰波在膨胀区中衰减,诱导激波阵面弯曲.两个支管中发生马赫反射,三波点迹线清晰可见.该传播特性是爆轰波的诱导激波和横波共同作用的结果.分叉口附近的胞格结构先消失再恢复,在无胞格和平衡胞格之间的区域存在细密胞格的过渡区,表征了在诱导激波与化学反应阵面分离后的区域中出现二次点火.P0=2.00kPa,水平支管中稳定自持爆轰能重建,垂直支管中爆轰熄灭.P0<2.00kPa,分叉口上游已不能形成稳定爆轰.还对胞格结构中的几个特征参数进行了测量,并初步分析了P0对这些参数的影响.     

气相爆轰波传播过程中的自点火效应

基于基元反应模型和单步反应模型,对直管道中H₂-air混合气体中爆轰波的传播过程进行了数值模拟,揭示了气相爆轰波传播过程中的自点火效应。利用数值模拟方法计算了不同爆轰模型的点火延迟时间,并得到了爆轰波三波点的传播过程以及所形成胞格结构的尺寸。结果表明,胞格宽度与点火延迟时间成正比;爆轰波诱导区内气体的点火延迟时间与三波点的运动周期基本一致。进一步对结果分析可知,爆轰波的自维持传播取决于点火延迟时间(表征化学反应的特征时间)和三波点的运动周期(表征流动的特征时间)的匹配;当二者相匹配时,经过前导激波压缩后形成的高温高压爆轰气体,在短时间内实现了自点火,同时释放出大量的能量推动了爆轰波的前进,即爆轰波的稳定自维持传播依靠其自点火机制。     

气相爆轰波在半圆形弯管中传播现象的实验研究

对气相(2H2/O2/Ar系统)爆轰波在半圆形弯管中的传播现象进行实验研究。用烟迹膜记录了弯管中爆轰波的胞格结构,采用压电传感器测量了沿弯管内外母线指定点的压力时间曲线,得到了爆轰波沿弯管内、外母线的平均速度和胞格尺寸的变化。结果表明:当平面爆轰波进入弯管后,受壁面的几何形状作用,诱导激波阵面发生弯曲。沿诱导激波阵面,自内母线到外母线方向,激波强度逐渐增大。同时,爆轰波后的化学反应区也受到影响,胞格尺寸发生较明显的变化。在本文条件下,当初压p08.00kPa,受扰动的爆轰波在弯管出口下游仍恢复为强度不变的稳定爆轰。胞格记录的三波点迹线表明:受扰动的爆轰波在出口段发生了马赫反射。实验结果还表明:当p0降至5.33kPa,平面稳定爆轰波经过半圆形弯管后,其强度发生衰减并直至出现熄灭。     

T型管内流动气体中爆轰绕射过程的数值模拟

基于带化学反应的二维Euler方程,对H₂、O₂、Ar体积比为2:1:1的混合气体系统在T型管内的爆轰绕射进行了数值模拟。用二阶附加半隐的Runge-Kutta法和五阶WENO格式分别离散欧拉方程的时间和空间导数项,采用9组分48步基元反应简化模型描述爆轰波在静止系统和流动系统中的传播过程,得到了温度、压力、典型组元H质量分数的分布及数值胞格结构。结果表明:在流动系统中,迎风面上波阵面为斜爆轰结构,静止系统两侧和顺风面上的波阵面为完全解耦的前导激波;在水平管中,波阵面与上下壁面经历一系列马赫碰撞后,最终形成正爆轰;在流动系统中,胞格结构明显向下游偏移;横向爆轰波的产生对爆轰波的再生起到了关键作用。     

旋转爆轰的三维结构和侧向稀疏波的影响

基于带化学反应的三维Euler方程,采用氢气-空气的9组分19步基元反应简化模型,对圆环形燃 烧室内的旋转爆轰进行了数值模拟,讨论了旋转爆轰波的三维结构及侧向稀疏波对旋转爆轰波阵面的影响。 数值结果表明,爆轰波能够以旋转方式沿预混气层稳定传播。在侧向稀疏波作用下,爆轰波阵面发生变形。 与理想的C-J爆轰相比,爆轰波强度和爆轰参数都有所下降。     

氢气/丙烷/空气预混气体爆轰性能的实验研究

通过采用压力传感器和烟灰板两种测试设备,开展了常温常压下氢气/丙烷和空气混合气体爆轰性能的实验研究。实验过程中观察到自持爆轰波,爆轰速度比值在0.99~1之间,爆轰压力比值在0.8~1.2之间。爆轰胞格尺寸在10~50 mm范围内,建立了爆轰胞格尺寸和化学诱导长度的关系式。随着丙烷不断添加,爆轰速度减小,而爆轰压力和胞格尺寸增加。这种变化趋势起初较快,而后变缓。因为起初氢气摩尔分数较大,混合气体趋向于氢气/空气的爆轰性能;而后因丙烷摩尔质量较大,丙烷逐渐起主要作用,混合气体表现出丙烷/空气的爆轰性能。     

Experiments on thermal characteristics of a natural circulation loop with latent heat energy storage under cyclic pulsed heat load

 An experimental investigation has been made of thermal characteristics of a rectangular, annular single-phase natural circulation loop with the inner tube filled with a solid–liquid phase change material (PCM) under cyclic pulsating heat load. A rectangular, annular loop of 150 cm in height and 75 cm in width was constructed with an annular gap of 0.6 cm, within which water was filled. The inner tube of the annular loop was filled with a PCM (n-Eicosene) or air. Under the cyclic pulsating heat load, temperature field within the water-filled annular loop with PCM- or air-filled inner tube was found to evolve into a steady periodic variation for the range of parameters considered. The water temperature and/or its fluctuating amplitude along the heated or cooled sections of the loop with the PCM-filled inner tube were found to be markedly lower than those measured in the loop with the air-filled inner tube under the identical conditions. On the other hand, along the insulated sections of the loop a somewhat minute difference in temporal variations of the water temperatures exists between the loops with PCM- and air-filled inner tube. In addition, at the outer wall along the cooled section, a time-periodic variation of temperature was detected in synchronizing with the pulsating heat load. Parametric effects of varying amplitude and time-period of the pulsating heat input, as well as of varying the inlet coolant temperature of the cooling jacket were investigated. Received on 30 June 2000 The authors are grateful for the support for this study from National Science Council of Republic of China in Taiwan under the Project Nos. NSC87-2212-E006-054 and NSC88-2212-E006-022.     

Effect of elastic wall motion on oscillatory flow of micropolar fluid in an annular tube

Summary Oscillatory flow of a micropolar fluid in an annular tube is investigated. The outer wall of the tube is taken to be elastic and the variation in the diameter of the elastic wall due to pulsatile nature of pressure gradient is assumed to be small. The wall motion is governed by a tube law. The nonlinear equations governing the fluid flow and the tube law are solved using perturbation analysis. The steady-streaming phenomenon due to the interaction of convected inertia with viscous effects is studied. The analysis, is carried out for zero mean flow rate. It presents the effects of the elastic nature of the wall combined with micropolar fluid parameters on the mean pressure gradient and wall shear stress for different catheter sizes and frequency parameters. It is found that the effect of micropolarity is of considerable importance for small steady-streaming Reynolds number. Also, it is observed that the relationship between mean pressure gradient and the flow rate depends on the amplitude of the diameter variation, flow rate waveforms and the phase difference between them.     

Comparison of critical conditions for DDT in regular and irregular cellular detonation systems

Abstract. The results of an experimental study of DDT in mixtures with regular and irregular detonation cellular structures are presented. Experiments were carried out in a tube 174 mm i. d. with obstacles (blockage ratios were 0.1, 0.3, and 0.6). Mixtures used were hydrogen–air and stoichiometric hydrogen–oxygen diluted with , Ar, and He. The critical conditions for DDT are shown to depend on the regularity of the cellular structure of test mixtures. The critical values of the cell sizes in Ar- and He-diluted mixtures are shown to be significantly smaller than those in -diluted mixtures. This means that systems with a highly regular detonation cellular structure have far less capacity for undergoing DDT compared to irregular ones with the same values of detonation cell sizes. Received 18 November 1999 / Accepted 15 May 2000     

Three-dimensional cellular structure of detonations in suspensions of aluminium particles

Recently, we have used scarce available data on the detonation cell size in suspensions of aluminium particles in air and oxygen to adjust the kinetic parameters of our two-phase model of detonations in these mixtures. The calculated detonation cell width was derived by means of two-dimensional (2D) unsteady simulations using an assumption of cylindrical symmetry of the flow in the tube. However, in reality, the detonation cells are three-dimensional (3D). In this work, we have applied the same detonation model which is based on the continuous mechanics of two-phase flows, for 3D numerical simulations of cellular detonation structures in aluminium particle suspensions in oxygen. Reasonable agreement on the detonation cell width was obtained with the aforementioned 2D results. The range of tube diameters where detonations in $\text{ Al/O}_2$ mixture at a given particle size and concentration would propagate in the spinning mode has been estimated (these results make a complement to our previous analysis of spinning detonations in Al/air mixtures). Coupling these results with the dependencies of detonation cell size on the mean particle diameter is of great interest for the understanding of fundamental mechanisms of detonation propagation in solid particle suspensions in gas and can help to better guide the experimental studies of detonations in aluminium suspensions. It is shown that the part of detonation wave energy used for transverse kinetic energy of both gas and particles is quite small, which explains why the propagation velocity of spinning and multi-headed detonations reasonably agrees with the ideal CJ values.     

爆震室壁面条件对脉冲爆震发动机爆震性能的影响

在爆震室内快速形成稳定传播的爆轰波是脉冲爆震发动机的关键．本文利用有限速率化学反应模型,考虑粘性、热对流,基于N-S方程对氢气与空气/氧气为反应混合物的爆震发动机爆震室内流场进行计算．从流场压力、速度、涡量、湍流动能等方面研究爆震室壁面条件对燃烧爆轰性能的影响,分析流场爆轰波压力与流场湍动能的关系,讨论可燃气体燃烧转爆轰的机理．结果表明：爆震室内燃烧爆轰机理受到化学反应能量释放、壁面摩擦效应、壁面与外界热交换的影响．在文中讨论的范围内,相比于半圆形和三角形的爆震室装置,矩形的爆震室增强装置能在更短的时间内得到较高的爆轰波压力和湍动能峰值．壁面粗糙层高度（粗糙度）影响爆震室的燃烧爆轰性质．当壁面粗糙度为0.15mm时,粗糙度对爆轰的激励作用大于抑制作用,能较快形成稳定的爆轰波,且推力为35.5N;随着壁面对流换热系数的增大,爆震室壁面的散热加剧．当壁面对流换热系数大于临界值2.6W/(m2·K)时,爆震室内不能形成稳定的爆震波．     

Heat transfer and two-phase flow characteristics of refrigerants flowing under varied heat flux in a double-pipe evaporator

A mathematical model based on the annular flow pattern is developed to simulate the evaporation of refrigerants flowing under varied heat flux in a double tube evaporator. The finite difference form of governing equations of this present model is derived from the conservation of mass, energy and momentum. The experimental set-up is designed and constructed to provide the experimental data for verifying the simulation results. The test section is a 2.5 m long counterflow double tube heat exchanger with a refrigerant flowing in the inner tube and heating water flowing in the annulus. The inner tube is made from smooth horizontal copper tubing of 9.53 mm outer diameter and 7.1 mm inner diameter. The agreement of the model with the experimental data is satisfactory. The present model can be used to investigate the axial distributions of the temperature, heat transfer coefficient and pressure drop of various refrigerants. Moreover, the evaporation rate or the other relevant parameters that is difficult to measure in the experiment are predicted and presented here. The results from the present mathematical model show that the saturation pressure and temperature of refrigerant decrease along the tube due to the tube wall friction and the flow acceleration of refrigerant. The liquid heat transfer coefficient increases with the axial length due to reducing the thickness of the liquid refrigerant film. Due to increase of the liquid heat transfer coefficient, increasing wall heat flux is obtained.Finally, the evaporation rate of refrigerant increases with increasing wall heat flux.