甲烷/空气预混气体火焰的传播特征

利用高速纹影摄像等技术探讨了密闭管道内不同当量比的甲烷/空气预混气体火焰的传播特征。结果表明，当甲烷含量接近当量值时，预混气体火焰传播中会发生火焰阵面由向未燃区弯曲到向已燃区弯曲的转折过程，逐渐由层流燃烧转变成湍流燃烧，并形成Tulip火焰结构；当甲烷含量偏离当量值一定程度时，预混火焰呈现出典型的层流燃烧特征，不会发生火焰阵面由向未燃区弯曲到向已燃区弯曲的转折过程。Tulip火焰结构形成于火焰传播速度迅速降低的区间里，且只有当减速阶段的最大加速度的绝对值大于某一数值时才能形成；Tulip火焰结构是预混火焰由层流燃烧向湍流燃烧转变的一个中间过程。     

燃烧室中悬跨圆柱的火焰绕流

利用高速阴影照相系统,对方形管道内火焰流经悬吊圆柱时所产生的变化及其诱导的流场进行了实验研究,并利用高精度PPM格式对上述现象进行了数值模拟,其中,湍流采用大涡模拟方法,燃烧采用旋涡破碎模型,圆柱边界采用沉浸边界法.实验结果与计算结果比较吻合,在一定程度上揭示了火焰阵面变形和加速的规律.根据研究结果,对火焰与悬吊圆柱相互作用过程中火焰的三维形状、流场中涡量和湍流强度的分布、火焰与流场的相互影响以及火焰传播轨迹等进行了讨论.     

密闭空间可燃气体爆炸超压预测

为避免密闭空间内可燃预混气体爆炸事故造成的伤害，对其进行较为准确的爆炸超压预测是抗爆设计和日常安全管理的关键。结合已有文献实验数据，利用光滑层流火焰传播理论模型建立了爆炸超压模型；对比发现，当体积较大时，光滑层流火焰传播理论模型存在较大的误差。较大体积密闭空间爆炸火焰传播过程中的不稳定性造成火焰前锋面褶皱并引起湍流燃烧，导致火焰前锋面表面积大幅增加，且在火焰传播过程中表现出自相似分形特征。依据褶皱及湍流火焰传播过程中的自相似分形特征，基于分形燃烧理论和相关经验数据，进一步建立了考虑可燃预混气体爆炸火焰褶皱及湍流火焰传播的爆炸超压预测模型，并与实验所得结果进行了对比。结果表明:当密闭空间体积较大时，利用褶皱及湍流火焰传播理论建立的爆炸超压模型进行峰值压力估算时，两种工况下实验所得和理论计算所得相对误差分别为10.4%和11.1%，较光滑层流火焰传播理论爆炸超压模型相比，误差分别减少了72.3%和50.6%。本文所建立理论模型与实验所得结果具有较好的一致性，在一定程度上可满足结构抗爆设计或日常安全管理的需要。     

方形管内楔形障碍物对火焰结构与传播的影响

通过实验与数值模拟方法对CH4/空气预混火焰在有楔形障碍物的卧式燃烧方管内的传播进行了研究。采用多镜头Cranz Schardin高速摄像机和压力传感器等实验设备获得了高清晰度的障碍物诱导火焰失稳的分幅时序照片以及障碍物背风表面压力变化曲线。数值模拟则基于RANS方法与EDU-Arrhenius燃烧模型,计算结果与实验结果基本相符,反映了火焰在管内传播与变形的详细过程。通过综合分析实验与计算结果,得到了由楔形障碍物导致的火焰加速与变形的内在机理,揭示了火焰传播过程中由层流燃烧向湍流燃烧转捩的本质。     

耦合火焰不稳定的爆炸超压预测

基于火焰不稳定和爆炸超压的耦合机制，通过向光滑火焰模型中引入褶皱因子，建立了褶皱火焰模型和湍流火焰模型，对密闭燃烧室内爆炸超压进行理论预测，且对比了绝热压缩和等温压缩对爆炸超压预测的影响规律。结果表明:在增强的流体动力学不稳定作用下，膨胀火焰失稳加剧，且在定容燃烧阶段形成胞状火焰；光滑火焰模型忽略了火焰不稳定，爆炸超压理论预测值比实验值偏低，且等温压缩下超压预测值低于绝热压缩下的预测值；湍流火焰模型高估了火焰褶皱程度，超压预测值远高于实验值；褶皱火焰模型可成功预测丙烷/空气爆炸压力和燃烧室体积V=25.6 m3的甲烷/空气爆炸压力；对于甲烷/空气爆炸，燃烧室体积V≤1.25 m3时，实验压力值介于褶皱火焰模型和绝热光滑火焰模型预测值之间。     

障碍物对预混火焰特性影响的大涡数值模拟

障碍物在预混气体火焰传播过程中对其结构及传播特性造成较大影响,对火焰的加速和爆燃转爆轰过程(deflagration-to-detonation transition, DDT)起到直接的促进作用。通过障碍物条件下可视管道中甲烷/空气预混火焰传播实验,捕获其火焰微观结构变化。采用三维物理模型,采用壁面自适应局部涡黏模型(wall-adapting local eddy-viscosity, WALE)的大涡模拟(large eddy simulation, LES),并用火焰增厚化学反应模型(thickened flame model, TFM)对实验过程进行重现。分析开口管道中预混火焰翻越障碍物后的复杂流场变化,并分析层流向湍流转变过程的特点。揭示了在障碍物影响下预混火焰扰动失稳现象的直接原因,是由障碍物引发的3个气流涡团同时作用而形成Kelvin-Helmholtz不稳定及Rayleigh-Taylor不稳定现象耦合作用所导致。     

多孔障碍物对预混火焰传播的影响

以甲烷为代表性气体，研究了半封闭管道中设置多孔障碍物对可燃气体爆炸火焰传播的影响，基于大涡模拟对实验进行了重现，对比了实验与模拟中火焰传播过程的形状、位置及速度，分析了模拟结果中火焰穿过障碍物前后的流场和表面积变化，给出了衡量火焰褶皱程度的指标及算法。结果表明:大涡模拟结果与实验结果有较好的一致性；火焰在存在障碍物的管道内传播，经历层流快速膨胀、受阻回流、湍流快速发展和脉动减速4个阶段，各阶段火焰依次分别呈现加速、减速、二次加速、二次减速的波动变化；当可燃气体在开口与点火位置同端的管道内爆炸，火焰在接近障碍物时，受管道封闭端和障碍物约束显著，而出现脉动回流现象；火焰穿过多孔障碍物后，传播速度骤升至峰值，较未穿过障碍物前的最大速度可增加58.7%；障碍物是导致火焰面破碎以及面积褶皱率增大的直接原因，火焰褶皱率最大可达44.8%，比未穿过障碍物前的最大褶皱率增大39.27%。     

射流扩散火焰闪烁模态及频率特性的实验研究

层流扩散火焰的闪烁是一种经典的火焰不稳定性现象,然而现阶段对其不稳定性模态及频率特性的研究尚不充分.文章利用高速摄影和粒子图像测速同步测量技术,对宽工况下准静态环境中的浮力主控和动量主控圆口射流扩散火焰的不稳定性及其频率进行了实验研究.实验发现,燃料射流流量增大会导致火焰失稳,所引起的闪烁现象可分为varicose模态和sinuous模态.流场测量结果表明,闪烁火焰的主要流动结构表现为分别位于火焰面内外的两个剪切层.火焰外剪切层卷起形成的大尺度环形涡的周期性产生、增长和脱落是导致火焰面周期性形变(即火焰闪烁)的原因.闪烁火焰具有准周期性,其流场具有统一的主频率,且该频率与火焰脉动频率一致,说明闪烁火焰在本质上是种整体流动不稳定性的体现.浮力主控火焰sinuous模态的频率比varicose模态高3 Hz左右,并且在燃料流量较大的工况下存在varicose和sinuous模态间的切换.浮力主控火焰的频率符合经典的1/2标度律,但不同模态对应着不同的标度律系数.动量主控火焰的频率显著偏离1/2标度率,且偏离程度随射流动量增加而增强.研究表明了流动不稳定性模态对于射流扩散火焰频率特性具有不可...     

湍流加速火焰的三维数值模拟

火焰在设有障碍物的管内传播时会自身加速，并可能导致爆炸。本文基于湍流κ-ε模型和改进的EBU—Arrhenius反应模型，对该现象进行了三维空间的数值模拟。计算结果反映了障碍物、湍流和火焰之间相互作用的正反馈机理，描绘了火焰在管内加速传播的三维图像。     

湍流的诱导及其对瓦斯爆炸过程中火焰和爆炸波的作用

在实验的基础上,研究了管内瓦斯爆炸过程中湍流的诱导及其对瓦斯爆炸过程中火焰和爆炸波的影响作用.研究结果表明,管道面积突变对瓦斯爆炸过程中湍流的产生具有重要影响.管道面积突变(变大、变小)时,产生附加湍流,并使下游火焰气流的湍流度增加,瓦斯爆炸过程中火焰的传播速度迅速提高,并可诱导激波的产生.在80×80mm等截面直管中(瓦斯浓度为理论上最猛烈的爆炸浓度9.5%),瓦斯爆炸最大火焰传播速度为40.8m/s,管内各点均为压力波信号,当管道加装一Φ300mm圆管形成面积突扩11倍和突缩11倍两断面后,面积突扩处(L/D=22)火焰速度增大5.05倍,达到64.4m/s,面积突缩处(L/D=28)火焰速度为156.0m/s, 增大4.55倍,并在L/D=48倍处形成激波(超压1.6976atm、波速416.7m/s),在L/D=98倍处,激波强度最大.在面积突变管内加装加速环可使瓦斯爆炸过程中湍流度加剧,火焰的传播速度更高,激波生成的位置(L/D=28)、最强点位置(L/D=70)均前移,激波强度增大.研究结果对指导现场如何防治瓦斯爆炸,减轻瓦斯爆炸的威力具有一定的指导意义.     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

汶川地震诱发文家沟巨型滑坡 碎屑流基本特征及成因机制初步分析


2008年5月12日,汶川M80地震在四川省绵竹市清平乡文家沟内诱发一巨型滑坡。通过现场调查得知,滑坡前后缘高差455m,厚度20～30m,滑面为基岩层面,初始方量2750×107m3。滑体在运动中转化为碎屑流。滑坡-碎屑流总的水平运动距离为4022m,垂直运动距离为1443m,遗留的堆积物体积达5×107m3。滑坡距映秀—北川断裂仅36km,位于其下盘,地震烈度达XI度。滑坡导致文家沟中48人遇害,并形成一条完整的地震次生地质灾害链。初步分析表明滑坡启动速度快,滑坡向碎屑流转化过程明显、地点明确。碎屑流运动过程复杂,伴有强烈的“气垫效应”和“前缘气浪冲击效应”。作者认为,文家沟滑坡的高启动速度是长持时强烈地震动作用的结果,与山体的猛烈碰撞是导致滑体解体并转化为碎屑流的原因。     

Preface: Orbits around asteroid

One of the core issues in modern celestial mechanics is the orbital dynamics in the near-regime gravitational field of as- teroids, which provides deep insights into the mathematical nature of a class of nonlinear systems, and plays as a critical basis for in situ explorations of different science goals. Lots of efforts have been made to reveal the characteristics of orbital motion in the vicinity of asteroids, and to improve the skills of asteroid research in methodology.     

Effet de l'anisotropie élastique cristalline sur la distribution des facteurs de Schmid à la surface des polycristauxInfluence of crystalline elasticity anisotropy on Schmid factor distribution at the free surface of polycrystals

Experimental studies of the plasticity mechanisms of polycrystals are usually based on the Schmid factor distribution supposing crystalline elasticity isotropy. A numerical evaluation of the effect of crystalline elasticity anisotropy on the apparent Schmid factor distribution at the free surface of polycrystals is presented. Cubic elasticity is considered. Order II stresses (averaged on all grains with the same crystallographic orientation) as well as variations between averages computed on grains with the same crystallographic orientation but with different neighbour grains are computed. The Finite Element Method is used. Commonly studied metals presenting an increasing anisotropy degree are considered (aluminium, nickel, austenite, copper). Concerning order II stresses in strongly anisotropic metals, the apparent Schmid factor distribution is drifted towards small Schmid factor values (the maximum Schmid factor is equal to 0.43 instead of 0.5) and the slip activation order between characteristic orientations of the crystallographic standard triangle is modified. The computed square deviations of the stresses averaged on grains with the same crystallographic orientation but with different neighbour grains are a bit higher than the second order ones (inter-orientation scatter). Our numerical evaluations agree quantitatively with several observations and measures of the literature concerning stress and strain distribution in copper and austenite polycrystals submitted to low amplitude loadings. Hopefully, the given apparent Schmid factor distributions could help to better understand the observations of the plasticity mechanisms taking place at the free surface of polycrystals. To cite this article: M. Sauzay, C. R. Mecanique 334 (2006).     

Preface: Current research progress on mechanics of high speed rail

The construction and operation of high-speed rail（HSR）grid within the past two decades in China,in terms of the scale,may be possibly comparable to any national-wide construction in the history of China,even the Great Wall.By counting railways with commercial train service at the speed of200 km/h,China has the world’s longest HSR network with over 19 369.8 km of track in service today and this number     

ACTA MECHANICA SINICA 2014 ANNUAL SUBJECT INDEX

正~~     

Nonlinear behavior of matrix-inclusion composites under high confining pressure: application to concrete and mortar

This paper is devoted to a micromechanics-based simulation of the response of concrete to hydrostatic and oedometric compressions. Concrete is described as a composite made up of a cement matrix in which rigid inclusions are embedded. The focus is put on the role of the interface between matrix and inclusion which represent the interfacial transition zone (ITZ). A plastic behavior is considered for both the matrix and the interfaces. The effective response of the composite is derived from the modified secant method adapted to the situation of imperfect interfaces. To cite this article: T.H. Le et al., C. R. Mecanique 336 (2008).