基于五阶WENO格式的燃气在药床中流动过程二维两相流研究

为研究内弹道初始阶段中心点火管燃气在膛内药床中的流动特性和传播规律,设计了可视化点传火实验平台,并进行了膛内假药床的点传火实验。基于加权本质无震荡(weighted essentially non-oscillatory, WENO)格式,构造了膛内轴对称二维内弹道两相流模型,对膛内燃气在假药床中的流动过程进行数值模拟。计算结果与可视化实验结果符合较好,全局压力平均误差为5.35%。表明数值计算准确地描述了燃气流动特性,完整地呈现了点火管燃气在假药床中的发展过程。在点火初始阶段,膛内压力径向效应明显,气相沿径向传播较快,药床药粒基本不会发生运动;随着燃气逐渐在膛内传播,膛内压力呈现径向一致、轴向梯度分布的特征,在压力梯度作用下,气相轴向速度开始占据主导,径向速度在膛底和中部区域减小为零,而固相速度随气相速度变化而变化;气相在到达弹底前,由于固相颗粒的壅塞,会提前出现速度反向波动现象。     

大尺度泄爆构件对室内爆燃压力影响的实验研究

通过在尺寸为2 m×1.2 m×0.6 m的腔体一端安装泄爆板研究不同乙烯浓度下大尺度泄爆构件对腔体内压力变化的影响。选用开启静压不同的2种泄爆板，在乙烯体积分数为4%~11%条件下进行实验，得到3种典型压力时程曲线。实验发现:泄爆结构的实际击穿压力大于静载作用下的击穿压力，且存在使实际击穿压力最大的最佳乙烯体积分数；泄爆构件的开启时间会对腔体内压力变化产生重要影响，低浓度条件下开启时间可达数十毫秒，化学当量比条件下开启时间仅为数毫秒，此时李克山模型针对长方体容器大尺度泄爆仍具有较好适用性，乙烯浓度较高时大尺度泄爆构件会因泄压面积过大造成外部空气大量进入腔体并与未燃气体再次反应发生二次爆炸，高浓度条件下增大泄压面积反而会因二次爆炸对结构造成破坏。     

氮气喷出对管道瓦斯爆炸的阻爆研究

为探究喷出氮气对瓦斯爆炸火焰传播的抑制能力，设置三种氮气喷头布置方式来进行阻爆实验，采用的氮气喷出压力有0.1、0.2、0.3、0.4和0.5 MPa，爆炸开始后喷射氮气，爆炸结束后氮气立刻关闭。结果表明，单喷头距泄压口20 cm时，各压力下喷出的氮气都未能阻爆，但火焰在整个管道内的平均传播速度随氮气压力增大而减小；单喷头距泄压口35 cm时，喷气压力0.5 MPa下成功阻爆，其他喷气压力下未能阻爆；双喷头喷气时，喷气压力0.3、0.4、0.5 MPa情况下都能够阻爆，且喷气压力越大，火焰被阻止的位置越靠前。阻爆的实现，需要氮气在阻爆位置将管道截面上的预混气稀释到可燃极限以下，因而氮气量是影响稀释的重要参数。单喷头时，喷头距离泄压口远更易于实现阻爆。采用双喷头时，氮气区扩大，阻爆所需氮气压力、氮气总量比单喷头时都大为降低。     

Optimal design of gas distributor in fluidized bed for synthesis of silicone monomer

The structure of the gas distributor is closely related to the production efficiency of organosilicon monomers. To improve the production efficiency of organosilicon monomers, this study uses Eulerian-Eulerian two fluid model and proposes a design formula for the gas distributor to optimize the gas distributor. It is proposed that the pressure drop of the gas distributor, the velocity nonuniformity coefficient, the relative standard deviation of the solid holdup, and the solid particle dispersion coefficient are used to evaluate the performance of the gas distributor. The results show that the performance of the gas distributor is significantly improved when the opening ratio Φ = 0.53% is optimized to Φ = 0.18%, in which the relative standard deviation of the solid holdup is reduced by 22%, and the solid particle dispersion coefficient is reduced by 40%. On this basis, this article studies the influence of different arrangements of vent holes on gas-solid fluidization characteristics. The results show that the circular arrangement of vent holes is helpful to the mixing of gas and solid.     

螺旋槽端面微间隙高速气流润滑密封特性

考虑入口气流压力损失和出口阻塞效应,建立了微间隙端面高速气体润滑密封分析数学模型,对螺旋槽端面微间隙高速气流润滑密封特性进行研究.重点分析了不同密封间隙、密封压力和转速等工况条件下,入口压力损失和出口阻塞效应对开启力、泄漏率及气膜刚度等密封特性参数的影响规律.结果表明:高速气体阻塞效应使出口压力高于环境压力,压力损失使入口气膜压力下降,导致泄漏率和气膜刚度明显下降,并使开启力增加.随着密封压力和密封间隙的增加,阻塞效应增强,导致泄漏率和气膜刚度显著降低.密封压力10 MPa时,泄漏率降低可达20%,气膜刚度的下降可达30%以上.     

Effect of distributed gas injection on aerodynamic characteristics of a body of revolution in a supersonic flow

Results of experimental and numerical investigations of the effect of gas injection through a permeable porous surface on the drag coefficient of a cone-cylinder body of revolution in a supersonic flow with the Mach number range M _h = 3–6 are presented. It is demonstrated that gas injection through a porous nose cone with gas flow rates being 6–8% of the free-stream flow rate in the mid-section leads to a decrease in the drag coefficient approximately by 5–7%. The contributions of the decrease in the drag force acting on the model forebody and of the increase in the base pressure to the total drag reduction are approximately identical. Gas injection through a porous base surface with the flow rate approximately equal to 1% leads to a threefold increase in the base pressure and to a decrease in the drag coefficient. Gas injection through a porous base surface with the flow rate approximately equal to 5% gives rise to a supersonic flow zone in the base region.     

考虑活塞回复的侧向后喷武器两相流数值模拟

为研究活塞回复运动对火药燃气流动的影响,基于两相流理论对活塞控制侧向后喷武器的发射过程进行了数值模拟研究。考虑控制侧向后喷通道开闭的活塞-弹簧系统的往复运动,建立了结合膛内气固两相流、活塞腔内流固耦合和侧向排气管内气体瞬态流动的武器发射过程数学模型,并将数值模拟结果与相关文献进行了比较验证。得到了该武器发射过程中膛内流场分布与稀疏波传播特性,并与普通武器的膛内流场进行了对比分析。进一步研究了活塞回复运动对火药燃气流动和减后坐效率的影响。结果表明:相对于不考虑活塞的回复运动,在弹丸初速都降低1.52%的情况下,因为活塞回复关闭后喷通道,其减后坐效率由38.86%下降到32.88%,说明在此类武器研究中,不可忽视活塞回复运动。     

瓦斯/空气预混气体爆燃流场测试系统多参数耦合同步控制实验方法

为了更精确地获得爆炸激波管内瓦斯/空气预混气体爆燃过程中,激波形成过程、压力和火焰传播速度以及火焰与惰性阻燃剂相互作用的流场演化图像。通过分析激波管测试系统中多个目标的时间响应特征及控制方式,利用超高速相机、光电倍增管、时间延时器、固态继电器、电荷放大器和数据采集系统等设备,设计实验方案,分别对激波管中瓦斯/空气预混气体爆燃高压点火系统的响应时间和惰性介质阻燃剂喷射系统的响应时间进行测试。实验结果表明电火花点火的响应时间为微秒量级,而阻燃剂喷射系统的响应时间为毫秒量级,以响应时间为依据,通过设置精确的延迟时间实现多目标同步控制,为完成激波管内瓦斯/空气预混气体爆燃过程的微观流场显示奠定基础。     

Influence of gas injection on phase inversion in an oil–water flow through a vertical tube

An experimental study has been made of the influence of gas injection on the phase inversion between oil and water flowing through a vertical tube. Particular attention was paid to the influence on the critical concentration of oil and water where phase inversion occurs and on the pressure drop increase over the tube during phase inversion. By using different types of gas injectors also the influence of the bubble size of the injected gas on the phase inversion was studied. It was found that gas injection does not significantly change the critical concentration, but the influence on the pressure drop is considerable. For mixture velocities larger than 1 m/s, the pressure drop over the tube increases with decreasing bubble size and at inversion can become even larger than the pressure drop during the flow of oil and water without gas injection.     

连通容器泄爆过程的影响因素

以甲烷/空气混合物为研究对象,开展了连通容器气体泄爆影响因素的实验研究。结果表明:连通容器泄爆片泄爆时,随着破膜压力和量纲一泄压比的减小,大、小球容器的最大泄爆压力均增大;在等量纲一泄压比条件下,随着连接管道长度的增加,传爆容器的最大泄爆压力增大。连通容器无膜泄爆时,大球点火条件下,无论管长如何,起爆容器和传爆容器均比单个容器最大泄爆压力大。小球点火条件下,当管道长度为0.45 m时,起爆容器和传爆容器的最大泄爆压力均小于单个容器。连通容器无膜泄爆且量纲一泄压比相同时,当管道长度为0.45 m时,大、小容器内的最大泄爆压力基本相等;当管道长度为2.45 m时,大容器点火时,传爆容器最大爆炸压力大于起爆容器,但小容器点火时,起爆容器最大泄爆压力大于传爆容器;当管道长度为4.45和6.45 m时,传爆容器最大泄爆压力均大于起爆容器。     

Effects of a countercurrent gas flow on falling-film mode transitions between horizontal tubes

A liquid film falling between horizontal tubes is known to take the form of droplets, jets or sheets, depending on the liquid flow rate; the form of the flow is the so-called “falling-film mode”. Although previously neglected in studies of mode transition, a countercurrent gas flow often exists in falling-film heat exchangers, and its effect on the liquid flow might be important: it could impact the flow regime, lead to local “dryout,” and decrease the heat transfer rate. Experiments are conducted to explore the effects of a countercurrent gas flow and liquid feeding length on falling-film mode transitions for a liquid flowing over horizontal tubes. The effects on mode transition are shown to depend on fluid properties and are explained in terms of unsteadiness and film thickness. In general, transition hysteresis is reduced with an increasing gas velocity. A correlation is developed to predict the countercurrent gas flow effects on falling-film mode transitions. The liquid feeding length can affect mode transitions in quiescent surroundings and when a countercurrent gas flow imposed.     

均压槽结构形状对静压干气密封性能影响分析

静压干气密封(DGS)中的均压槽起着均布压力和二次节流的作用,开展了典型形状均压槽静压DGS的性能对比和结构优选.基于静压气体润滑理论,建立了圆形、椭圆形、扇形和环形等四种典型均压槽静压型DGS的几何模型和数学模型,采用有限差分法求解获得了四种均压槽静压型DGS端面的膜压分布和稳态密封性能参数,分析了径向开槽比和周向开槽比对四种均压槽静压型DGS密封性能的影响规律;以获得较大的密封开启力和气膜刚度为目标,计算得到了均压槽径向开槽比和周向开槽比的优选值范围.结果表明:当均压槽径向开槽比0.15Wd0.45时,环形均压槽静压型DGS可获得较大的开启力和气膜刚度,其他三种均压槽静压型DGS的径向开槽比优选值范围为0.3Wd0.45;扇形和椭圆形均压槽静压型DGS具有相似的密封性能,其密封性能仅次于环形均压槽静压型DGS;当均压槽周向开槽比0.6Lθ1.0时,扇形均压槽静压型DGS可获得较大的开启力和气膜刚度.     

Radial diffusion in the poiseuille flow of a gas mixture

In a number of experiments (see [1], in which experimental papers are listed), diffusion has been observed in the radial direction in the process of flow of a mixture along tubes at low pressures. The heavier molecules accumulate near the tube axis. The attempt made in [1] to explain this phenomenon by the influence of the Burnett contribution to the diffusion did not lead to success, and the Burnett terms in the radial diffusion velocity indicate a motion of heavy molecules away from the tube axis. In the present paper, a complete analysis is given of this phenomenon. We consider the problem of the flow of a mixture along a cylindrical tube of finite length for given pressure difference p between its ends. On the basis of the hydrodynamic equations of the Burnett and super-Burnett approximations, a consistent asymptotic (with respect to the small parameter ) solution is given; = (p/p)R/L is the relative change in the pressure along the tube at a distance of order R (R and L are the radii and length of the tube). Radial diffusion occurs in the quadratic approximation in . It is shown that the radial diffusion velocity contains new terms not present in [1]; these are due to the inhomogeneity of the temperature and the pressure over the tube section, the expansion of the gas, and the super-Burnett correction to the diffusion velocity. The most important is the thermodiffusion term, which is determined by the hydrodynamic equations of the Navier-Stokes approximation. The remaining terms have order relative to it of Kn² (Kn = 1 /R is the Knudsen number, and 1 is the mean free path of the molecules). The expression obtained for the diffusion velocity agrees in sign with the experiment.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 91–96, September–October, 1979.I am grateful to G. E. Skvortsov, who drew my attention to this problem, and Yu. N. Grigor'ev for discussing the results.     

An aerodynamic wake in a hypersonic free-molecular gas flow

An experimental investigation was made of the form of an aerodynamic wake and the flow in it for the case where the flow is set up using a point source and is propagated in a motionless residual gas. It is shown that the geometry of the wake is determined by the mutual arrangement of the source and model, that the ground pressure is negligibly small, and that the flow in the wake is part of the oncoming flow, scattered by the residual gas. The value of the transverse velocity ratio exceeds by several orders of magnitude the value of the longitudinal ratio.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 183–186, January–February, 1976.     

Investigations of gas and particle dynamics in first generation needle-free drug delivery devices

Abstract. Transdermal powdered drug delivery involves the propulsion of solid drug particles into the skin by means of high-speed gas-particle flow. The fluid dynamics of this technology have been investigated in devices consisting of a convergent-divergent nozzle located downstream of a bursting membrane, which serves both to initiate gas flow (functioning as the diaphragm of a shock tube) and to retain the drug particles before actuation. Pressure surveys of flow in devices with contoured nozzles of relatively low exit-to-throat area ratio and a conical nozzle of higher area ratio have indicated a starting process of approximately 200 s typical duration, followed by a quasi-steady supersonic flow. The velocity of drug particles exiting the contoured nozzles was measured at up to 1050 m/s, indicating that particle acceleration took place primarily in the quasi-steady flow. In the conical nozzle, which had larger exit area ratio, the quasi-steady nozzle flow was found to be overexpanded, resulting in a shock system within the nozzle. Particles were typically delivered by these nozzles at 400 m/s, suggesting that the starting process and the quasi-steady shock processed flow are both responsible for acceleration of the particle payload. The larger exit area of the conical nozzle tested enables drug delivery over a larger target disc, which may be advantageous. Received 12 March 2000 / Accepted 8 June 2000     

Modeling of three-phase heavy oil–water–gas bubbly flow in upward vertical pipes

A bubbly gas–bubbly oil flow pattern may occur when water, heavy oil and gas flow simultaneously in vertical pipes in such a way that water is the continuous phase. In this work, a one-dimensional, thermal, transient two-fluid mathematical model, for such flow, is presented. The model consists of mass, momentum and energy conservation equations for every phase whose numerical solution is based on the finite difference technique in the implicit scheme. The model is able to predict pressure, temperature, volumetric fraction and velocity profiles. For accurate modeling of multiphase flows, the key issue is to specify the adequate closure relationships, thus drag and virtual mass forces for the gas and oil phases were taken into account and special attention was paid on the gas–oil drag force. When this force was included into the model it was found that: (1) such force had the same order of magnitude than the oil drag force and both forces were smaller than the gas drag force, (2) the pressure, gas and oil velocities and gas and oil volume fraction profiles were affected, (3) the numerical stability was increased. The model predictions are in agreement with experimental data reported in literature.     

Propagation of pressure waves in two-phase flow

We deal with a pressure wave of finite amplitude propagating in a gas and liquid medium or in the fluid in an elastic tube. We study the effects of pipe elasticity on the propagation velocity of the pressure wave. Pressure waves of finite amplitude progressing in the two-phase flow are treated considering the void fraction change due to pressure rise. The propagation velocity of the two-phase shock wave is also investigated, and the behavior of the reflection of the pressure wave at the rigid wall is analyzed and compared to that in a pure gas or liquid. The results are compared to experimental data of a pressure wave propagating in the two-phase flow in a vertical shock tube.     

Numerical simulation of flow in Hartmann resonance tube and flow in ultrasonic gas atomizer

The gas flow in the Hartmann resonance tube is numerically investigated by the finite volume method based on the Roe solver.The oscillation of the flow is studied with the presence of a needle actuator set along the nozzle axis.Numerical results agree well with the theoretical and experimental results available.Numerical results indicate that the resonance mode of the resonance tube will switch by means of removing or adding the actuator.The gas flow in the ultrasonic gas atomization (USGA) nozzle is also studied by the same numerical methods.Oscillation caused by the Hartmann resonance tube structure,coupled with a secondary resonator,in the USGA nozzle is investigated.Effects of the variation of parameters on the oscillation are studied.The mechanism of the transition of subsonic flow to supersonic flow in the USGA nozzle is also discussed based on numerical results.     

Experimental analysis of the swirling flow in a Francis turbine draft tube: Focus on radial velocity component determination

The draft tube of a hydraulic turbine is the component where the flow exiting the runner is decelerated, thereby converting the excess of kinetic energy into static pressure. In the case of machine refurbishment of an existing power plant, most of the time only the runner and the guide vanes are currently modified. For financial and safety reasons, the spiral casing and the draft tube are seldom redesigned, even if these components present some undesirable behaviour. In some cases, the installation of an upgraded runner leads to a peculiar and undesirable efficiency drop as the discharge is increased above the best efficiency point value. It is found to be related to a corresponding sudden variation in the draft tube pressure recovery coefficient at the same discharge.The swirling flow exiting the runner is complex and highly turbulent. The radial velocity is rarely measured because a quite complicated measurement setup is needed. However, this velocity component is greatly needed in order to properly initialize the numerical simulations, and its influence is important in spite of its small magnitude. Velocity measurements downstream of the runner include radial component made at CREMHyG (Grenoble) by LDV, and PIV techniques are presented. An analytical formulation for this velocity component based on the formulation for the conical diffuser and on the three vortices structure is proposed and compared with measurements.     

CFD-DEM simulation of spouting of corn-shaped particles

Three dimensionally coupled computational fluid dynamics (CFD) and discrete element method (DEM) were used to investigate the flow of corn-shaped particles in a cylindrical spouted bed with a conical base. The particle motion was modeled by the DEM, and the gas motion by the k-? two-equation turbulent model. A two-way coupling numerical iterative scheme was used to incorporate the effects of gas–particle interactions in terms of momentum exchange. The corn-shaped particles were constructed by a multi-sphere method. Drag force, contact force, Saffman lift force, Magnus lift force, and gravitational force acting on each individual particle were considered in establishing the mathematical modeling. Calculations were carried out in a cylindrical spouted bed with an inside diameter of 200 mm, a height of 700 mm, and a conical base of 60°. Comparison of simulations with experiments showed the availability of the multi-sphere method in simulating spouting action with corn-shaped particles, but it depended strongly on the number and the arrangement of the spherical elements. Gas–solid flow patterns, pressure drop, particle velocity and particle concentration at various spouting gas velocity were discussed. The results showed that particle velocity reaches a maximum at the axis and then decreases gradually along the radial direction in the whole bed. Particle concentration increases along the radial direction in the spout region but decreases in the fountain region, while it is nearly constant in the annulus region. Increasing spouting gas velocity leads to larger pressure drop, remarkably increased speed of particle moving upward or downward, but decreased particle concentration.