不同重力下90°弯管内气液两相流流型及流动特性研究

研究了不同重力条件下90°弯管内气液两相流流型分布形态及流动特性. 通过建立90°弯管内气液两相流流动的三维数学物理模型,采用VOF 方法,对10-6g₀, 10-4g₀, 10-2g₀, 1g₀ (g₀= 9.8m/s2) 重力下的90°弯管内气液两相流流型分布特征、截面空隙率、滑速比及气相尾部最大斜向角进行了比较分析. 研究结果表明:所建立的模型能够正确模拟不同重力条件下90°弯管内气液两相流流型和截面空隙率,并得到气液两相弯管二次流与单相二次流的不同特性. 随着重力水平的提高,90°弯管对气相流型的影响作用减弱,气相整体向弯管内侧积聚靠拢,弯管对尾部的斜向作用减弱.      

快速烤燃条件下B炸药战斗部的临界泄压面积

为了确定战斗部装药在快速烤燃条件下能稳定燃烧的临界泄压面积,基于质量守恒定律和气体状态方程,建立了战斗部壳体内考虑炸药初始温度和排气孔排气的气体压力增长模型。以B炸药圆柱战斗部为研究对象,研究了炸药意外点火后能稳定燃烧的A_V0/S_B(临界泄压面积/炸药外表面积)确定方法,并与实验值进行了比较。结果表明,本文建立的模型能够很好地预测B炸药战斗部的临界泄压面积。研究了战斗部炸药装药表面积、炸药初始温度、空气体积占比和炸药燃速对A_V0/S_B的影响,并将不同温度的模型预测值与实验值进行了比较。结果表明：炸药装药表面积对A_V0/S_B基本没有影响;A_V0/S_B与温度和炸药燃速成正相关,与空气体积占比成负相关;不同温度的模型预测值A_V0/S_B与实验值吻合较好。     

考虑时滞的主动悬架系统控制策略对比研究全文替换

本文以二自由度四分之一汽车悬架系统为研究对象,采用不同控制策略对考虑时滞的悬架系统控制特性进行研究,并对控制效果进行对比分析．首先,采用第二类拉氏方程建立考虑时滞的二自由度悬架控制系统的动力学模型;然后分别基于状态变换法和H_∞控制理论,设计系统的时滞反馈控制律．其中状态变换法主要通过系统状态变量的转换,将系统时滞控制方程转换成不显含时滞的动力学方程,然后采用传统的二次型最优控制方法对系统进行控制．H_∞控制主要通过使用Lyapunov-Krasovskii泛函和自由权矩阵法,设计考虑时滞的H_∞控制律;最后在Matlab/Simulink平台上对控制系统进行仿真分析,并在相同时滞下对两种控制结果进行对比．研究表明,两种控制策略在考虑时滞的情况下均可保证系统的稳定性,且H_∞控制相比于状态变换法有更好的控制效果．     

不同晶型CL-20热分解反应机理计算分析

基于反应力场,采用NPT、NVT系综和Berendsen方法,对ε-、β-和γ-CL-20超晶胞在不同温度下的热分解反应过程进行分子动力学计算。结果表明:3种晶型CL-20的初始分解路径均是五元环和六元环中的N-NO₂键断裂生成·NO₂自由基,生成的主产物包括N₂、H₂O、CO、CO₂、NO₂、NO₃、HNO₂、HNO₃、N₂O₅、N₂O₂和NO等。相同晶型的CL-20发生热分解生成的主产物的反应速率常数随温度的升高而增大。     

2021-08期封面

铝粉反应模型是对悬浮铝粉尘气-固两相爆轰进行数值模拟研究的关键。通过考虑铝粉燃烧产物氧化铝（Al₂O₃）在高温下的分解吸热反应,改进了铝粉的扩散燃烧模型。将该模型嵌入到三维的气-固两相爆轰数值计算程序中,分别对铝粉/空气混合物以及铝粉/氧气混合物的爆轰进行了数值模拟,计算得到的稳定爆轰波速度与实验结果、文献值均吻合较好,误差小于5.5%,表明改进的铝粉反应模型适用于不同氧化气体氛围中铝粉尘爆轰的模拟计算。此外,对两相爆轰参数及爆轰流场的物理量分布进行分析,获得了铝粉反应模型对爆轰波结构的影响规律。     

旋转爆震燃烧室轴向和周向长度对其出口流场压力和温度的影响

为了研究旋转爆震燃烧室与涡轮的匹配特性,利用二维欧拉方程数值研究了基于当量H₂/Air燃烧的旋转爆震燃烧室出口流场特性,对比分析了不同燃烧室轴向长度和周向长度出口总压脉动、总压畸变以及出口温度分布规律。结果表明:旋转爆震燃烧室在稳定工作状态下,其出口总压的脉动值会呈现周期性振荡;燃烧室尺度对发动机出口流场的不均匀性有很大影响,随着燃烧室轴向长度的增大或周向尺寸的减小,其出口总压脉动均值、畸变指数和出口温度分布系数均会减小,其出口流场均匀性提高。此外,爆震波高度随着周向尺寸的增大而增大;轴向尺寸对爆震波高度几乎不产生影响。     

火箭发动机燃烧过程的鲁棒非脆弱H∞控制

针对某液体火箭发动机燃烧室的燃烧过程,设计了鲁棒非脆弱 H_∞ 状态反馈控制器。首先,基于一种新型的时滞分割法和交互式凸组合技术,借助于构造一个包含四重积分项的Lyapunov-Krasovskii泛函(LKF),并利用新的积分不等式方法给出了LMI形式的时滞相关有界实判据;其次,在有界实判据的基础上,采用矩阵的合同变换以及变量替代技巧将燃烧过程非线性矩阵不等式线性化,通过求解线性矩阵不等式得到相应的非脆弱H_∞控制器的可行解。模拟结果验证了本文设计方法的有效性。     

大速差射流燃烧室中烟煤与贫煤燃烧的数值模拟

本文基于颗粒相的轨道模型,对大速差射流燃烧室中烟煤粉与贫煤粉的二维流动,混合及燃烧进行了数值模拟,模拟结果从两相耦合的角度,阐明了煤粉颗粒在燃烧室中运动的规律,煤粉与大速差射流诱导的中心气体逆流之间的混合及其对煤粉火焰稳定的影响,指出此种燃烧室中煤粉火焰稳定的回流区燃烧机理,气相流场及回流区的预报结果与实验符合良好。     

点火条件对密闭管道内预混氢气/空气燃爆特性的影响

基于流体动力学软件Fluent,开展数值模拟,研究点火位置（距管左端壁面100、200和500 mm）、点火温度（1 000、1 500和2 000 K）和点火面积（管左端壁面处半径为50、35和20 mm的点火域）等点火条件对1 000 mm密闭管道中预混氢气/空气（H₂/air）燃爆特性的影响。研究表明:点火位置距管左端壁面越远,中间节点处温度越高,温升越快;不同点火温度下管内最高温升速率基本同步,且提高点火温度,使得燃烧反应更剧烈,能提高管内气体温升速率,但却降低管内的压力峰值;点火面积越小,预混H₂/air燃烧前期温升越快。当采用半径为35 mm的点火域和点火位置距管左端壁面100 mm的点火方式时,预混H₂/air燃爆的各项参数相对较高。不同点火条件对密闭管内气体的动能和内能的影响规律类似于其对管内气体的流速和温度的影响规律,而对涡量的影响不明显。     

钝头体高超声速绕流底部失稳特征数值模拟

利用数值模拟方法对高超声速钝锥及Apollo返回舱底部尾迹流场进行了研究, 分析尾迹流动的失稳过程. 对钝锥模型, 在M_∞=6, Re=1.71× 10⁶(Re以球头半径为参考长度)条件下观察到了底部流动的不稳定性. 不添加任何扰动, 数值模拟首先得到的流动是稳定解, 在底部发展出一个主分离区和一个二次分离区, 流动是轴对称状态. 继续进行计算, 发现二次分离线率先变形, 底部流场发展出非定常周期流动. 对Apollo返回舱模型, 在相同条件下 (Re以前面圆弧半径为参考长度), 数值模拟首先得到的流动同样是稳定解, 出现以二次分离线率先变形为起始的结构失稳, 演化出周期性过程, 但持续时间较短, 很快出现了非周期非对称状态. 研究表明, 高超声速钝锥及Apollo返回舱底部流场均存在不稳定性问题, Apollo返回舱的底部流场更加不稳定.     

The pressure and temperature behind bodies with flat bases in a supersonic flow when inert and active gases are introduced into the base region

The effect of introducing active and inert gases into the base region behind two-dimensional and axisymmetric bodies on the pressure and temperature in that region is studied in the absence of heat transfer between the body and the base region. The flow in the mixing zone is assumed to be turbulent. To investigate the effect of introducing active gases on the parameters in the base region, a model of diffusion combustion is used. The base pressure is calculated by Korst's method [1]. For the velocity profile in the mixing zone an expression is used which results from integrating the equation of motion in von Mises variables. The temperature and concentrations of the components in the base region are determined from integral equations for the conservation of enthalpy and concentration, which make it possible to calculate the parameters in the base region for arbitrary flow rates of the active gas, including arbitrarily small rates. The results are given of calculating the base pressure when hydrogen and argon are blown into the base region. In the latter case computational and experimental results are compared.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 48–57, March–April, 1971.     

点火具点火射流蚀剥作用对底排药柱燃烧特性的影响

为了研究底排点火具射流蚀剥作用对底排药柱碎块脱落情况及其燃烧性能的影响,采用半密闭爆 发器实验装置模拟底排弹出膛口瞬态卸压工况,借助高速录像记录点火具点火与底排药柱燃烧的序列图像。 建立底排药柱在半密闭爆发器燃烧的数学模型,计算分析了不同碎块脱落质量引起的平均压力、质量流率和 燃烧时间的变化情况。研究结果表明:蚀剥作用发生在点火具点火初期,导致脱落的底排药剂碎块来不及燃 烧;蚀剥作用主要是强点火射流对底排药柱侵蚀与冲击造成的;碎块总量约占底排药柱初始质量的7%~ 9%,严重影响点火和燃烧的一致性。计算结果与实验值吻合较好。     

Study of gas recirculation in the base region of a model launcher

The technique and the results of a study of gas recirculation in the base region of a launcher are described. The investigation is based on the measurement of the concentrations of chemical substances on a model in a wind tunnel. The engine jets are modeled by the combustion products of a solid fuel.     

Functional mathematical model of a hydrogen-driven combustion chamber for a scramjet

A functional mathematical model of a hydrogen-driven combustion chamber for a scramjet is described. The model is constructed with the use of one-dimensional steady gas-dynamic equations and parametrization of the channel configuration and the governing parameters (fuel injection into the flow, fuel burnout along the channel, dissipation of kinetic energy, removal of some part of energy generated by gases for modeling cooling of channel walls by the fuel) with allowance for real thermophysical properties of the gases. Through parametric calculations, it is found that fuel injection in three cross sections of the channel consisting of segments with weak and strong expansion ensures a supersonic velocity of combustion products in the range of free-stream Mach numbers M_∞ = 6–12. It is demonstrated that the angle between the velocity vectors of the gaseous hydrogen flow and the main gas flow can be fairly large in the case of distributed injection of the fuel. This allows effective control of the mixing process. It is proposed to use the exergy of combustion products as a criterion of the efficiency of heat supply in the combustion chamber. Based on the calculated values of exergy, the critical free-stream Mach number that still allows scramjet operation is estimated.     

横向气流穿过垂直下落液滴场的二维流动数值计算

对Cl2/He混合气体横向穿过垂直下落的BHP(按重量25%的KOH,25%的H2O2及50%的H2O)液滴场的化学反应流动作了数值计算.模拟的流场是气体/液滴两相流场,在气相方程中,考虑了液滴与Cl2反应产生及释放O2(1Δ)的质量源项及表示液滴对气流阻碍作用的动量源项.由于气相动量小,液滴在下落过程横向偏移小,下落...     

Turbulence energy balance in reacting turbulent flows

The turbulence accompanying combustion and the propagation of detonation waves in gases has been studied theoretically and experimentally in many papers [1–8]. The attention of researchers has been concentrated on essential questions like how the turbulent flow field interacts with the kinetics of the chemical reaction and to what extent the process of chemical change is intensified, and how the turbulence itself is deformed by the heat released and the accompanying expansion of the gases. The various mechanisms proposed for these phenomena are based on various hypotheses concerning the structure of the combusion zone and the determinative stage of the interaction of the turbulence with the chemical-reaction kinetics. The mechanism of turbulence generation by combustion proposed in a number of papers [3–6] is based on the observation in turbulent flow of a weakly curved flickering laminar flame. This gives rise to a nonuniform flow field of the gas, part of the energy of which goes over into the energy of turbulent fluctuations. Other authors [7, 8] considered the turbulence field to interact with the chemical-reaction kinetics via a volume mechanism and suggested a criterion of turbulence intensification based on certain physical considerations, e.g., the condition for the intensification of thermogaskinetic oscillations proposed by Rayleigh [9]. In the present paper the problem is analyzed by introducing Kolmogorov's general equation for the turbulence energy balance in reacting turbulent flows [10]. In accordance with, this equation the turbulence energy can vary due to energy exchange between the turbulent motion and the mean gas flow as a result of the work on turbulent mass transport in the acceleration field of the mean flow, and due to the effect of pressure fluctuations on the rate of thermal expansion from the chemical reaction. Each of these effects is considered and analyzed.     

Experimental study on the unusual turbulence intensity distributions in rod-to-wall gap regions

The phenomenon of unusual turbulence intensity distributions in closely spaced rod bundles was experimentally studied. Fully developed turbulence intensities in rod-to-wall gap regions, formed by a single rod regularly mounted in a trapezoidal duct, were measured with a hot-wire anemometer. Starting from a rod setting giving a small bottom rod-to-wall gap, a series of gap spacings were achieved by moving the rod upward, resulting in larger bottom gaps and at the same time smaller top gaps. A region of high turbulence kinetic energy occurred first in the bottom gap region corresponding to the small bottom gap spacing and then shifted to the top gap region for the cases of small top gap spacings. There was no evidence of a high kinetic energy patch when the top gap size was extremely small. To unveil the mechanism behind this phenomenon, measurements on secondary velocities and energy spectra were conducted for one rod setting, and these results are described in detail. It was found that the measured secondary flow patterns could not explain the abnormal turbulence intensity distributions. However, energy density spectra revealed certain peak frequencies that support the argument that cross-gap large-scale eddy motion is of importance as both a transport process and production mechanism for subchannels formed by closely spaced rod bundles.     

开区注氮采空区自燃温度场的数值模拟研究

为了研究开区注氮采空区遗煤自燃和被抑制的力学过程,用非均质多孔介质中的渗流连续性方程、气体弥散方程和综合传热方程的联立,建立了采空区注氮非定常数值模型．结合实例,用迎风格式有限元方法求解．给出了采空区的漏风流态、氮气流态,描绘了采空区自燃过程中氧、CO浓度和温度分布的变化过程．计算考虑了瓦斯涌出和工作面推进的影响．得到注氮条件下采空区高温区形状变窄;随着注氮量的提高,使自然发火期逐渐变长,直至不自燃．理论计算与实际情况吻合．     

激波诱导火焰失稳与爆轰的条件研究

采用九阶WENO和十阶中心差分格式数值求解激波与火焰作用过程,考察了激波强度、火焰尺寸对激波与球形火焰作用过程的影响。结果表明,增大激波强度或火焰尺寸均可在流场中引发爆轰,但激波强度的影响更大,并且其引发的爆轰可使火焰迅速膨胀,放热率提高,从而影响燃烧特性;此外,爆轰波传播过程中会迅速消耗可燃预混气,合并原有的反射激波,并在流场中形成局部高压区,极大地改变流场结构。     

Laminar mixed convection in horizontal concentric annuli with non-uniform circumferential heating

 Steady, laminar, mixed convection in the fully developed region of horizontal concentric annuli has been investigated numerically for the case of non-uniform circumferential heating. Two heating conditions were studied, one in which a 180^∘ arc encompassing the top half of inner surface of the inner cylinder is uniformly heated while the bottom half is kept insulated, and the other in which the heated and the insulated surfaces were reversed. The fluid flow and heat transfer characteristics were found to be affected by the heating conditions. For the investigated range of the governing buoyancy parameter, the modified Grashof number (Gr^*), it was found that bottom heating arrangement gives rise to a vigorous secondary flow, with the result that the average Nusselt numbers are much higher than those for pure forced convection. On the other hand, the local Nusselt numbers are nearly circumferentially uniform. In the case of top heating arrangement, a less vigorous secondary flow is induced because of temperature stratification, with average Nusselt numbers that are substantially lower than those for bottom heating and with large circumferential variation of the local Nusselt number. Received on 15 March 2000