强约束球形装药反应裂纹传播和反应烈度表征实验

炸药燃烧的高温高压气体产物可以进入基体裂纹中引发炸药表面热传导燃烧，形成所谓的对流燃烧。在一定约束条件下，不断上升的气体压力反过来又使炸药基体产生更多的裂纹，为对流燃烧提供更多的通道和燃烧表面积，快速生成大量产物气体导致高烈度反应现象的产生。本文中设计了一种新型强约束球形装药中心点火实验，针对一种HMX为基的PBX炸药，对高烈度反应条件下燃烧裂纹传播和反应增长过程进行了观测，实验中采用测得的反应压力和壳体速度历程对反应烈度进行了量化表征。在带窗口结构中，早期炸药中的燃烧裂纹不可见；中期燃烧裂纹扩展到药球表面时，先形成4条沿经线方向近似对称的主裂纹，随后环向贯通并扩展到整个药球表面；最后的剧烈反应造成强烈发光。上述反应演化经历低压增长阶段约为100 μs，之后伴随着壳体变形膨胀产生剧烈的反应，此时产物压力在约10 μs时间内超过1 GPa，并形成约20%相对于裸炸药爆轰的超压输出。在全钢结构中，20 mm厚的壳体膨胀速度最大可达到500 m/s，此时壳体完全破裂。     

基于孔与裂隙网络模型的平行微裂隙对驱油的影响规律研究

认识双重多孔介质中油水两相微观渗流机制是回答形成什么类型的裂隙网络可提高油藏采收率的关键. 微裂隙的分布可以提高多孔介质的绝对渗透率,但对于基质孔隙中的流体介质,微裂隙的存在会引起多孔介质中局部流体压力和流场的变化,导致局部流动以微裂隙流动为主,甚至出现窜流现象,降低驱油效率. 本文基于孔与裂隙双重网络模型,在网络进口设定两条平行等长且具有一定间隔的微裂隙,分析微裂隙的相对间隔(微裂隙之间距离/喉道长度)和微裂隙相对长度(微裂隙长度/喉道长度)对于微观渗流特征的影响. 结果表明:随微裂隙相对长度的增加,出现驱油效率逐渐降低,相对渗透率曲线中的油水共渗区水饱和度和等渗点增加,油水两相的共渗范围减小等现象;随着微裂隙之间相对间隔增大,周围越来越多的基质孔穴间的压力差减小,在毛管压力的限制下,驱替相绕过这些区域,而导致水窜现象.      

对流边界层湍流特性的数值研究

采用大涡模拟方法研究了存在逆温层的情况下大气对流边界层的湍流特性。实际大气边界层中出现逆温层是较常见的，逆温层会导致大气边界层湍流结构的变化，从而影响大气的湍流扩散和输运特性。本文比较了不同逆温梯度的工况，着重分析了逆温层对边界层中热量逆梯度输运（counter gradient heat transportation，CGHT）的影响。计算结果表明：逆温梯度越大，对流边界层的发展越受到抑制；逆温层高度降低会影响整个对流边界层的温度抬升；逆温梯度越大，垂直速度方差越小；在逆温梯度较大的情况下，其逆梯度输运区域要略微低一些，初步分析认为是由于逆温层对热对流的抑制造成的；对于逆温层高度不同的情况，高度越低的逆温层对逆梯度输运的抑制作用更明显。     

浮力对混合对流流动及换热特性的影响

用热线和冷线相结合的技术测量垂直圆管内逆混合对流流体的平均速度、 温度以及它们的脉动. 较详细地研究了浮力对逆混合对流的流动特性和传热特性的影响. 评 估了实验中采用的冷线测量温度补偿速度探头温度敏感的影响. 逆混合对流的传热结果用无 量纲参数Ω (Ω= Grd / Red2 )来表示，其中，基于管道直 径的雷诺数Red变化范围为900~18000, 浮力参数Ω变化范围为 0.004899~0.5047. 研究结果表明，浮力对逆混合对流的换热有强化作用. 随着葛拉晓夫数Grd的增加，温度脉动，流向雷诺正应力和流向温度通量增 大，并且在靠近壁面的流体区域尤其明显. 热线与冷线相结合的技术适合于研究非绝热的流 动测量，可以用于研究浮力对流动和换热特性的影响.     

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

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

Ti-6Al-4V弹体破坏模式对冲击反应的影响研究

Ti-6Al-4V材料是武器结构轻量化时的重要替代材料，其冲击反应将可能增加战斗部毁伤威力，但目前缺乏对其冲击反应条件及反应机理的研究。本文将采用试验与理论分析方法，研究结构破坏模式对Ti-6Al-4V材料冲击反应的影响，获得其冲击反应条件及反应机理。设计并开展了钛合金弹（头部与壳体均为钛合金）与复合弹（头部碳/碳复合材料、壳体空心钛合金圆柱）正侵彻混凝土试验，撞击速度在222~1008 m/s之间。钛合金弹激发了剧烈的氧化冲击反应，但复合弹未产生冲击反应。破坏模式宏细观分析显示，钛合金弹侵彻后宏观结构基本完整，仅表面发生摩擦磨损，以细观组织剪切变形为主要失效模式，形成尺寸在微米量级至百微米量级的颗粒碎片，碎片个数可高达3×106。复合弹的钛合金空心圆柱被撕裂成块，撕裂面沿剪切带方向发展，碎块尺寸在毫米或以上量级，个数至多百余个。碎片供氧和供热的效率均与碎片尺寸成反比，而特定供氧与供热条件下，碎片尺寸足够小是Ti-6Al-4V材料发生冲击反应的必要条件，这是钛合金弹发生冲击反应而钛合金空心圆柱无法激发冲击反应的本质原因。在具备冲击反应必要条件的前提下，碎片个数越多，冲击反应烈度越高。     

快速热点火熔奥梯铝炸药燃烧转爆轰实验研究

以熔铸型含铝混合炸药熔奥梯铝为对象，研究铸装含铝混合炸药快速热点火后的燃烧转爆轰特性。建立了快速热点火燃烧转爆轰实验平台，由实验装置（加热装置、约束钢管、炸药）、压力测试系统、光纤测速系统组成；加热装置加热15 mm厚45钢钢板，峰值温度大于1 100 ℃，温升速率为85～95 ℃/s。开展了快速热点火带壳熔奥梯铝炸药燃烧转爆轰实验，由加热装置加热约束钢管内熔奥梯铝炸药，炸药化学反应阵面压力和传播速度分别由压电性高压压力传感器和光纤探针测定；实测阵面压力约1 GPa，传播速度最大约2 600 m/s。由光纤数据获得炸药化学反应阵面传播轨迹，通过特征线方法获得冲击形成点，半定量给出冲击形成距离大于850 mm；并比较了管体破片质量实测值与炸药完全爆轰时破片平均质量计算值，实测值远小于计算值。综合实测化学反应阵面传播速度和压力、冲击形成距离分析、破片质量比较，可确定熔奥梯铝炸药没有发生完全爆轰，其化学反应状态为爆燃。另外，采用Adams和Pack模型、CJ燃烧模型，都能够半定量的预估冲击形成距离和燃烧波后压力，为实验设计提供依据，但CJ燃烧模型的计算结果更接近于实测值。     

惰性粉体对油页岩粉尘爆炸火焰的抑制性能和作用机理研究

为了研究惰性粉体对油页岩粉尘爆炸火焰的抑制性能和作用机理,利用粉尘爆炸火焰传播测试系统,选取了五种常用惰性粉体和两种不同油页岩粉尘进行了爆炸火焰抑制实验。通过对爆炸火焰长度、最低惰化比和火焰形态结构的分析,结合惰性粉体的TG-DTG-DSC热特性曲线,系统研究了惰性粉体对油页岩粉尘爆炸火焰的抑制性能和作用机理。研究结果表明,惰性粉体对两种油页岩粉尘爆炸火焰的抑制性能优劣排序为:ABC干粉＞Al(OH)₃＞Mg(OH)₂＞NaHCO₃>岩粉,且两种惰性粉体均对桦甸油页岩（HDOS）的抑爆性能优于龙口油页岩（LKOS）;本文建立了惰性粉体对油页岩粉尘爆炸火焰的抑制机理物理模型,并分析了作用机理,通过作用机理分析表明:高效抑爆粉体应具有热稳定性较好（分解温度在200～400 ℃）,吸热量大,且分解中间态产物能够与燃烧反应活性自由基相结合发挥化学抑制作用等特点。     

非平衡流模拟参数ρ∞L有效性的数值分析

在非平衡流相似律现有研究成果的基础上,通过对化学反应特征时间、具体流场区 域热化学状态、辐射输运方程的分析,并结合若干典型算例,讨论双体碰撞反应的非 平衡尺度效应模拟参数ρ∞L对高空高超声速飞行器全流场适用的条件.结果表明,当模拟参数ρ∞L值低于一定范围时, 在高超声速流场非平衡区中离解非平衡起主导作用,此时模拟参数 ρ∞L不仅对全流场的无量纲参数分布有效,而且对流场电离特性、模型气动热特性、分子光谱范围的驻 点辐射谱强度都是有效的.模拟参数ρ∞L的有效性与自由流速度、不同类物理量和流场不同区域均相关: 当自由流速度增加时, 对全场有效的ρ∞L要求ρ∞L 值减小;与无量纲压力、密度和温度相比, 对化学组元特别是离子组元分布有效的ρ∞L所要求的ρ∞L 值更小; 与模型头部区相比, 在身部区有效的ρ∞L所要求的ρ∞L值更小.     

受限空间内甲醇喷雾液滴形成及其爆炸特性

为防控工业喷雾爆炸和完善喷雾爆炸测试方法,在20 L球形喷雾爆炸测试系统内，实验研究了不同环境压力、喷射压力及浓度下的甲醇喷雾液滴形成及爆炸特性规律。结果表明:增大喷射压力更易致使甲醇破碎成微小液滴，甲醇喷雾液滴爆炸极限范围变宽；环境压力的增大导致甲醇喷雾液滴粒径变大，喷雾液滴爆炸极限范围变窄，一定程度上可以有效抑制甲醇泄露可能导致的次生衍生事故发生。当爆炸容器内环境压力为0.1 MPa、喷射压力为2.1 MPa、甲醇喷雾浓度为356.4 g/m3、甲醇液滴索太尔平均直径为2.5 μm时，爆炸特性参数（最大爆炸压力、最大爆炸压力上升速率及层流燃烧速度）在上述拐点处取得最大值；小粒径（1～15 μm）的液滴在外界能量作用下，更易被点燃，且爆炸过程中瞬态物理化学反应更为迅速和剧烈；较大粒径(22 μm以上)的液滴会出现点火困难现象，然而点火成功后，爆炸特性参数均随甲醇喷雾浓度增加而增加，呈现近似线性规律，此时液滴粒径对上述爆炸特性参数的影响可以忽略。研究结果有助于理解喷雾液滴爆炸规律、完善相应测试方法和安全设计。     

A pore scale study on turbulent combustion in porous media

This paper presents pore scale simulation of turbulent combustion of air/methane mixture in porous media to investigate the effects of multidimensionality and turbulence on the flame within the pores of porous media. In order to investigate combustion in the pores of porous medium, a simple but often used porous medium consisting of a staggered arrangement of square cylinders is considered in the present study. Results of turbulent kinetic energy, turbulent viscosity ratio, temperature, flame speed, convective heat transfer and thermal conductivity are presented and compared for laminar and turbulent simulations. It is shown that the turbulent kinetic energy increases from the inlet of burner, because of turbulence created by the solid matrix with a sudden jump or reduction at the flame front due to increase in temperature and velocity. Also, the pore scale simulation revealed that the laminarization of flow occurs after flame front in the combustion zone and turbulence effects are important mainly in the preheat zone. It is shown that turbulence enhances the diffusion processes in the preheat zone, but it is not enough to affect the maximum flame speed, temperature distribution and convective heat transfer in the porous burner. The dimensionless parameters associated with the Borghi–Peters diagram of turbulent combustion have been analyzed for the case of combustion in porous media and it is found that the combustion in the porous burner considered in the present study concerns the range of well stirred reactor very close to the laminar flame region.     

Heating zone dynamics in in situ combustion

The integrodifferential equation of the quasisteady regime of a moving in situ combustion front is obtained and its exact solution is constructed in a particular case; the possibility of the heat generated at the combustion front being projected into the region ahead of the front is analyzed and the heating zone dynamics in the reservoir and the surrounding rock are investigated. In a number of studies of in situ combustion it is assumed that an increase in the water-air factor or, what amounts to the same thing, an increase in convection velocity in the reservoir leads to the total transfer of the heat into the region ahead of the combustion front [1–3]. In [3] the area of the heating zone ahead of the combustion front was calculated in accordance with the Marx-Longenheim model [4]. Below, on the basis of exact solutions of model problems it is shown that in the case of quasisteady Newtonian heat transfer between the surrounding medium, when the latter is assumed to be a thermal reservoir, i.e., maintain a constant temperature, this projection of heat is possible if the convection velocity exceeds the velocity of the combustion front. In the case of unsteady heat transfer in accordance with the Leverrier model there is no total projection of heat into the region in question; in the steady-state regime a limited heating zone, proportional in depth to the square of the difference of the convection and combustion front velocities, is formed ahead of the front.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 166–172, July–August, 1987.The author wishes to thank V. M. Entov for his valuable advice and useful discussions.     

Natural convection heat transfer on surfaces of copper micro-wires

The natural convection heat transfer characteristics and mechanism for copper micro-wires in water and air were investigated experimentally and numerically. The wires with diameters of 39.9, 65.8 and 119.1 μm were placed horizontally in water inside of a sealed tube and in air of a large room, respectively. Using Joule heating, the heat transfer coefficients and Nusselt numbers of natural convection for micro-wires in ultra pure water and air were obtained. A three dimensional incompressible numerical model was used to investigate the natural convection, and the prediction with this model was in reasonable accordance with the experimental results. With the decrease of micro-wire diameter, the heat transfer coefficient of natural convection on the surface of micro-wire becomes larger, while the Nu number of natural convection decreases in water and air. Besides, the change rate of Nu number in water decreases apparently with the increase of heat flux and the decrease of wire diameter, which is larger than that in air. The thickness of boundary layer on the wall of micro-wire becomes thinner with the decrease of diameter in both water and air, but the ratio of boundary layer thickness in water to the diameter increases. However, there is almost no change of this ratio for natural convection in air. As a result, the proportion of conduction in total heat transfer of natural convection in water increases, while the convective heat transfer decreases. The velocity distribution, temperature field and the boundary layer in the natural convection were compared with those of tube with conventional dimension. It was found that the boundary layer around the micro-wire is an oval-shaped film on the surface, which was different from that around the conventional tube. This apparently reduces the convection strength in the natural convection, thus the heat transfer presents a conduction characteristic.     

Diffusive Effects on Recovery of Light Oil by Medium Temperature Oxidation

Volatile oil recovery by means of air injection is studied as a method to improve recovery from low permeable reservoirs. We consider the case in which the oil is directly combusted into small products, for which we use the term medium temperature oil combustion. The two-phase model considers evaporation, condensation and reaction with oxygen. In the absence of thermal, molecular and capillary diffusion, the relevant transport equations can be solved analytically. The solution consists of three waves, i.e., a thermal wave, a medium temperature oxidation (MTO) wave and a saturation wave separated by constant state regions. A striking feature is that evaporation occurs upstream of the combustion reaction in the MTO wave. The purpose of this paper is to show the effect of diffusion mechanisms on the MTO process. We used a finite element package (COMSOL) to obtain a numerical solution; the package uses fifth-order Lagrangian base functions, combined with a central difference scheme. This makes it possible to model situations at realistic diffusion coefficients. The qualitative behavior of the numerical solution is similar to the analytical solution. Molecular diffusion lowers the temperature of the MTO wave, but creates a small peak near the vaporization region. The effect of thermal diffusion smoothes the thermal wave and widens the MTO region. Capillary diffusion increases the temperature in the upstream part of the MTO region and decreases the efficiency of oil recovery. At increasing capillary diffusion the recovery by gas displacement gradually becomes higher, leaving less oil to be recovered by combustion. Consequently, the analytical solution with no diffusion and numerical solutions at a high capillary diffusion coefficient become different. Therefore high numerical diffusion, significant in numerical simulations especially in coarse gridded simulations, may conceal the importance of combustion in recovering oil.     

甲烷/煤尘复合爆炸火焰的传播特性

为揭示甲烷/煤尘复合爆炸火焰的传播机理,利用气粉两相混合爆炸实验系统,在低于甲烷爆炸下限条件下,采用高速摄影机记录火焰传播图像,通过热电偶采集火焰温度,研究了煤尘种类以及甲烷体积分数对甲烷/煤尘复合火焰传播特性的影响。结果表明:挥发分是衡量煤尘燃烧特性的主导因素;随着煤尘挥发分的升高,燃烧反应增强,火焰传播速度升高,火焰温度升高;挥发分含量差异较小时,水分含量越低,燃烧反应越剧烈;在相同条件下,焦煤的燃烧反应强度最高,其次为长焰煤,最后为褐煤;随着甲烷体积分数的增加,煤尘颗粒的燃烧可由释放挥发分的扩散燃烧转变为气相预混燃烧,燃烧反应增强,火焰传播速度和火焰温度显著升高;热辐射和热对流作用促进煤尘颗粒热解,释放挥发分进行燃烧反应,维持复合火焰的持续传播;随着混合体系中甲烷体积分数的增加,混合爆炸机制由粉尘驱动型爆炸转为气体驱动型爆炸,燃烧反应增强;甲烷/煤尘复合爆炸火焰可由未燃区、预热区、气相燃烧区、多相燃烧区和焦炭燃烧区5部分组成,湍流扰动导致燃烧介质空间分布存在差异,使得燃烧区无规则交错分布。     

Wave propagation in the combustion of a low-concentration aerial suspension

In the general case the convective combustion of aerial suspensions is described by the equations of mechanics of multiphase media [1]. If the volume particle content is neglected and it is assumed that in the initial stage of convective front propagation the particles are stationary, and that during combustion their temperature is constant, then the equations for describing the combustion process reduce to the equations of gas dynamics for a distributed supply of heat and mass [2, 3]. The equations and model constant mass burning rate kinetics are used to solve the plane one-dimensional problem of the combustion of an aerial suspension in part of a region bounded on one side by a fixed wall. A small parameter proportional to the mass concentration and the heat value of the fuel is introduced. The method of matched asymptotic expansions [4] is used to construct a uniformly applicable first approximation. The solution obtained describes the wave propagation in aerial suspension combustion processes. The resulting pattern includes an inclined compression wave propagated with the speed of sound followed by a convective hot reaction product front whose propagation velocity is much less (in conformity with the small parameter introduced) than the speed of sound.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 63–73, March–April, 1986.     

Numerical investigation of transient natural convection heat transfer of freezing water in a square cavity

In this study, a transient heat transfer process of freezing water inside a two-dimensional square cavity has been investigated numerically. Water was used as a phase-change medium, and the numerical model has been created with control volume approach by using C++ programming language. To be able to accelerate the numerical calculations, CUT (Consistent-Update-Technique) algorithm has been implemented in the numerical code. Span-wise variations of the vertical component of the velocity have been represented in comparison with the experimental measurements from the literature at various vertical positions to examine the accuracy of the numerical scheme. The influence of natural convection has been considered by comparing the conduction and convection dominated solidification under same boundary conditions. Comparative results have been obtained regarding time-wise variations of the cold wall temperature and the dimensionless effectiveness. Moreover, the streamlines and isotherms have been represented to understand the differences between the conduction and convection driven phase change processes.Results indicate that natural convection becomes remarkable and has different forms at the initial periods of the phase change process. Increasing the effect of natural convection in the cavity increases the cooling rate of water. Near the density inversion temperature of water (4°C), temperature variations fluctuate and counter currents observed in the domain.     

气体辅助注射成型过程动力学分析

基于气体穿透机理和Hele-Shaw流动模型,对管状模腔中气体冲破熔体前沿形成中空制品气辅成型过程进行了研究,推导出反映充模流动压力梯度比、非牛顿幂率指数等影响因素与计算表层熔体厚度比之间关系的数学公式,建立了熔体前沿和气体前沿速度与位移演化关系的数学模型,分别得出了牛顿流体和非牛顿流体选取不同影响参数时熔体前沿和气体前沿速度、位移的演化曲线.模拟结果表明,所建数学模型能较好的反映熔体前沿和气体前沿速度、位移演化关系.在气体冲破熔体前沿以前,气体接近匀加速运动,前沿位移梯度逐渐增加;熔体前沿的速度几乎保持不变,位移随时间接近线性增长.当气体冲破熔体前沿时,熔体和气体前沿的速度和位移均急速上升.     

Magnetic control of natural convection in the horizontal Bridgman configuration using a spectral method: transversal plan

The present study investigates the electromagnetic braking of buoyancy convective flows occurring in differentially heated cavities, filled with low Prandtl, dilute, incompressible and electrically conducting alloys, and subjected to a constant horizontal temperature gradient. In practice, such flows known as ‘Hadley circulation’ are relevant in material processing technologies, such as the horizontal Bridgman configuration. A collocation spectral numerical method is developed to solve the two-dimensional Navier–Stokes equations, modelling the flow phenomena occurring in such configurations, using a vorticity–stream function formulation. The two components of the velocity are deduced from the stream function and the temperature distribution is obtained through the resolution of the energy conservation equation. The results in terms of velocity and temperature distributions for a given Grashof number are obtained for various Hartmann numbers and show that as the Hartmann number increases, the electromagnetic braking of the flow is observed. Moreover, the results illustrate the changes affecting the flow structure which becomes quasi-parallel in the core region of the cavity for sufficiently high values of Ha and the onset of the Hartmann and parallel layers along the boundaries. Also, with increasing Ha, the isotherms are less affected by the convective flow and become parallel to the vertical walls indicating that heat transfer is mainly achieved by conduction.     

油页岩原位注热开采储层有效热解区变化规律分析

针对油页岩原位注热开采过程中储层有效热解区变化规律不清,实际热解效果无法准确判断难题,采用数值模拟方法,以抚顺油页岩储层为研究对象,建立了油页岩原位注热开采热流固耦合力学模型,与前人结果对比,验证了模型可靠性。重点考察水力压裂裂缝通道短路问题,分析得到了油页岩原位注热开采过程中储层有效热解区、储层有效热解区中地应力、注汽压力及沉降量随注热时间变化规律。结果表明,过热蒸汽沿水力压裂裂缝流动不会出现裂缝通道短路现象,过热蒸汽可通过水力压裂裂缝加速油页岩储层热解;采用过热蒸汽对流加热油页岩储层效率高,只需1年能使96%的油页岩储层达到热解所需温度;油页岩储层有效热解区中部形成应力集中区,最大地应力为21.6 MPa;热解后靠近注热井处岩层发生沉降,热解2年后最大沉降量达0.85 m。所得结论对现场油页岩原位注热开采有参考意义。