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

炸药燃烧的高温高压气体产物可以进入基体裂纹中引发炸药表面热传导燃烧，形成所谓的对流燃烧。在一定约束条件下，不断上升的气体压力反过来又使炸药基体产生更多的裂纹，为对流燃烧提供更多的通道和燃烧表面积，快速生成大量产物气体导致高烈度反应现象的产生。本文中设计了一种新型强约束球形装药中心点火实验，针对一种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的增加，温度脉动，流向雷诺正应力和流向温度通量增 大，并且在靠近壁面的流体区域尤其明显. 热线与冷线相结合的技术适合于研究非绝热的流 动测量，可以用于研究浮力对流动和换热特性的影响.     

W/ZrNiAlCu亚稳态合金复合材料破片对RHA靶的侵彻释能特性

为研究W/ZrNiAlCu亚稳态合金复合材料破片对RHA靶板的侵彻释能特性，采用高速摄影弹道枪侵彻实验和能量方程、Avrami-Erofeev方程理论分析的方法，对破片的侵彻释能过程、侵彻规律、释能规律进行了研究。结果表明，破片在撞击并贯穿靶板的过程中激发了材料的燃烧反应，在靶板前方和后方产生了明显火光，随着撞击速度增加，火光范围增加、亮度提高；破片撞击速度、冲塞体速度的关系符合采用能量法推导的包含质量损失的破片侵彻公式，破片理论弹道极限速度为987.1 m·s?1；在实验速度范围内，材料反应效率随着冲击压力的增加而增加，与实验现象吻合。     

铸铁滑动摩擦副表面温度对其磨损机制的影响

在铸铁摩擦融油润滑和线接触的滑动条件下，测定了试样摩擦磨损表面层的温度及其温度分布梯度，考察了表面层的相变特点、组织结构和表貌特征，并且分析了温度对磨损机制转化的影响。可以将摩擦磨损表面的热效应分为３个区：弱热效应区（对应于塑性流动磨损区）、热效应区（对应于氧化磨损区）和强热效应区（对应于严重磨损区）。在临界载荷前，表面层的γ相含量增加，但一临界载荷时，由于表面层的剥落，γ相的含量明显减少。摩擦磨     

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

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

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

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

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.     

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.     

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

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

石英玻璃杆Taylor撞击实验的数值模拟

采用离散元方法模拟石英玻璃杆Taylor撞击问题,再现了其破坏过程:在撞击端,杆以压缩失效波的形式破坏;在自由端,出现了密集的拉伸层裂破坏. 分析表明:层裂是失效波阵面应力快速下降引起的追赶卸载波,与弹性压缩前驱波在自由端反射引起的迎面卸载波相互作用的结果;随着撞击速度的增大,撞击端失效波造成的压缩破坏区域损伤程度增大,反射端层裂破坏损伤区域减小. 进一步对失效波阵面的结构变化及其波速问题进行了研究,发现失效区域随着扩张变成一段裂纹逐渐由密到稀的区段,将此区段分为高损伤区和低损伤区,研究发现由稀疏微裂纹组成的低损伤区的前端面传播速度和弹性前驱波速基本相同,为固定值;而高损伤区前端面的裂纹密度随着传播距离的增加变稀,直至过渡为低损伤区,其传播具有显著的速度衰减、端面模糊直至停止的过程. 高损伤前端面的平均速度随着撞击速度的增大而增大,并逐渐趋近于弹性波速. 最后与已有实验做了对比,发现实验中高速摄影观察到的玻璃中"失效波"阵面实际上是高损伤前端面,而稀疏的低损伤微裂纹很难捕捉.