多次激波诱导正弦扰动预混火焰界面失稳的数值研究

激波诱导火焰失稳是实际中常见的现象,为深入研究火焰失稳特性,采用三维单步化学反应的Navier-Stokes方程和9阶weighted essentially non-oscillatory (WENO)的高精度格式,对不同马赫数的入射激波及其反射激波多次诱导正弦型预混火焰界面失稳的现象进行了三维数值模拟,并对计算结果的可靠性进行了验证。研究结果显示,在激波的多次作用下,火焰界面的演变过程主要受Richtmyer-Meshkov (RM)不稳定特性和化学反应特性的双重影响,且随着入射激波强度的增强,上述2种特性均得到进一步强化。为构造体现反应性RM不稳定特性的参数,根据火焰界面混合区平均涡量和化学反应速率,提出了表征界面受不稳定性和化学反应影响的量纲一参数。通过分析发现,在同一入射激波强度下,该参数的对数形式随入射激波和反射激波的多次作用呈基本相同的线性增长趋势;对不同马赫数的入射激波,该参数对数形式的线性增长率也基本一致。这样的变化表明该量纲一参数能够反映反应性RM不稳定过程中火焰界面发展的内在规律。     

球形重质气体物理爆炸特性

基于大涡模拟方法,结合高阶混合格式,对高压重质的SF₆球形气云在空气中爆炸进行了模拟。数值模拟表明,爆炸产生的激波经过气体分界面时分为透射激波以及反射稀疏波,透射激波导致气体分界面处Richtmyer-Meshkov失稳增强,从而加速了2种气体的混合,而反射的稀疏波经过汇聚,在球心处形成二次激波,在该强激波作用下,流场区域基本呈现湍流形态。     

Richtmyer-Meshkov不稳定性扰动增长率冲击模型的微尺度效应

本文通过考虑流体介质输运性质对激波理论的影响,针对微型激波管中的RichtmyerMeshkov(RM)不稳定性,分析了微尺度效应对三种冲击模型扰动增长率的影响。一维气体动力学计算结果表明:流动尺度对于Richtmyer模型、Meyer-Blewett(M-B)模型和Vandenboomgaerde-Mügler-Gauthier(V-M-G)模型中的线性扰动增长率有显著影响。当流动尺度由宏观状态逐渐减小至微尺度时,三种模型的扰动增长率均会经历从少量增长到明显下降,然后迅速上升的过程。微尺度条件下,V-M-G模型的扰动增长率相比宏观尺度有显著提高。与其他两种模型相比,修正后的V-M-G模型更合理地描述了微尺度效应对于线性扰动增长率的影响。此外,对于修正后的V-M-G模型,当入射激波马赫数较低时,扰动增长率受微尺度效应的影响更为明显。     

Air／SF6界面RM不稳定性诱导的物质混合诊断

不同流体伴随激波的作用,在不同密度介质的界面处往往存在激波诱导的界面不稳定性,即RM(Richtmyer-Meshkov)不稳定性,由此不稳定诱导了物质间的相互混合。文中研究了低马赫(1.23)激波作用Air/SF6界面RM不稳定性问题。Air/SF6初始界面由厚度为1～2μm的硝化纤维薄膜相隔得到,利用阴影测试法研究了Air/SF6界面在1.23马赫数激波冲击下,界面混合宽度随时间的发展过程,以及反射激波作用后混合宽度的再增长。实验结果表明混合宽度线性发展前期与经验公式吻合较好。     

柱形汇聚激波冲击球形重气体界面的实验研究

采用高速纹影法实验研究了柱形汇聚激波与球形重气体界面相互作用的 Richtmyer-Meshkov不稳定性问题. 激波管实验段基于激波动力学理论设计, 将马赫数为1.2 的平面激波转化为柱形汇聚激波, 气体界面由肥皂膜分隔六氟化硫(内)和空气(外)得到. 采用高速摄影机在单次实验中拍摄激波运动的全过程, 对柱形激波的形成进行了实验验证, 并进一步观测了汇聚激波与球形气体界面相互作用过程中的波系发展和气体界面变形以及反射激波同已变形界面二次作用的流场演化. 结果表明: 当柱形汇聚激波穿过气泡界面以后, 气泡左侧界面极点沿激波传播方向保持匀速运动, 气泡右侧界面发展成为射流结构, 气泡主体发展成为涡环结构; 在反射激波的二次作用下, 流场中无序运动显著增强并很快进入湍流混合阶段.     

空气自由场中强爆炸冲击波传播二维数值模拟

编写了适用于模拟具有高密度比、高压力比的强激波问题的二维柱对称多介质流体计算程序。利用有限体积方法求解流体的Euler方程组,采用level set方法捕捉爆炸产物与空气的运动界面,并通过求解物质界面两侧Riemann问题的精确解来计算爆炸产物与空气之间的数值通量。研制了三角形网格自适应技术来实现网格的自动加密和粗化,在保证捕捉激波峰值的前提下有效地提高了计算效率。利用计算程序对1 kt TNT当量的空气自由场强爆炸问题进行数值模拟,计算得到的峰值超压、冲击波到达时间等物理参数与点爆炸理论结果基本一致。     

用改进的ghost fluid方法模拟强激波与界面的相互作用

为了解决原来的ghost fluid方法在计算强激波和界面相互作用时界面附近出现的速度和压力振荡问题,对原来的ghost fluid方法进行了改进,通过在界面处构造Riemann问题并求出界面的压力和速度,ghost fluid流体的压力和速度分别用界面的压力和速度代替,ghost流体的密度通过熵常数外推得到。改进的ghost fluid保持了原来的ghost fluid的简单性,对一维强激波与气-气、气-液界面的相互作用问题以及射流问题进行了数值计算,得到了分辨率较高的计算结果。     

再冲击载荷作用下流动混合的数值模拟

基于多介质的流体体积分数VOF（volume of fluid）方法和PPM（piecewise parabolic method）方法,给出和发展了可用于多介质粘性流体动力学的数值计算方法和计算代码MVPPM（multi-viscosity-fluid piecewise parabolic method）。为了检验和验证此计算代码,对某激波管实验再冲击载荷作用下的流体动力学不稳定性及其导致的混合过程进行了数值模拟,计算结果与实验结果一致。同时还研究了激波反射冲击作用下流体混合区的演化情况,在反射激波和混合区相互作用的瞬间,混合区的宽度明显减小,之后又迅速增大;另外,混合增长率与初始扰动的频谱有很大关联。通过对有粘性（分子动力学粘性）和无粘性结果的对比,发现粘性对混合区的影响很小。     

A STUDY ON JET PHENOMENON OF R22 GAS CYLINDER UNDER THE IMPACT OF SHOCK WAVE

重物质气泡在平面激波作用下出现射流是一个非常有趣的物理现象, 它与气体属性、激波强度密切相关, 气泡内外波系结构复杂, 值得进行深入探讨. 论文采用水平集方法结合真实虚拟流体方法捕捉物质界面, 对R₂₂圆柱气泡在平面激波作用下的射流现象进行了数值模拟, 详细阐述了波系结构在气泡内外的演变过程, 总结了导致气泡射流现象形成的典型波系结构特征, 从波系结构角度解释了射流形成的原因.论文也对不同强度的激波冲击R₂₂圆柱气泡的情况进行了归纳, 发现激波强度对射流现象的出现以及射流结构的尺寸特征有重要影响, 得出了 R₂₂圆柱气泡能否出现射流结构的临界激波强度, 以期对实验结果起到指导与验证作用.     

Richtmyer-Meshkov不稳定性流体混合区发展的实验研究

使用矩形激波管，在马赫数分别为$M=1.5$和1.7的条件下实验研究了气/液界面 上(即Atwood数$A$接近1时)由Richtmyer-Meshkov不稳定性引起的流体混合现象. 得到 了气/液界面上Richtmyer-Meshkov不稳定性后期流体混合区域宽度随时间的发展呈现出线 性关系的结果，即$h \propto t$. 比较了不同马赫数和初始扰动下的发展情况，发现当 马赫数增加时，同一时间混合区域 宽度随之增加，混合区域宽度增长变快；而相比于波长差别不大的弱多模态初始扰动(无人 为干扰界面), 当界面初始扰动获得随机外界干扰时，界面混合区域具有较大的宽度以及增 长速度. 并且增加激波马赫数和初始扰动多模态性，流体混合程度更为剧烈.     

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

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

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


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).     

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).     

Dynamics and Breakup of Viscoelastic Liquids (A Review)

The phenomena of hydrodynamic breakup of liquid jets, drops, films, bridges, and filaments are reviewed for liquids with viscoelastic properties. The reasons for breakup are capillary instabilities, collisions with rigid obstacles, and other forms of dynamic action. The relationship between the properties of the liquids and the features of the breakup process is discussed.     

Flow stability analysis and excitation using pulsating jets

Classical flow stability applied to transition from laminar to turbulent flow may also describe the behavior of vorticity fluctuations created by a pulsating jet placed along a solid boundary. A numerical laminar flow experiment involving a pulsating jet placed along the surface of a duct with flow separation downstream, resulted in eliminating most part of the separated flow region. Applying the same approach to a turbulent flow, it was possible to develop a turbulent stability flow formulation and apply successfully turbulent pulsating jet flow separation control. To cite this article: D. Skamnakis, K. Papailiou, C. R. Mecanique 333 (2005).