基于新型解耦算法的激波诱导燃烧过程数值模拟

对一种模拟化学非平衡流动的时间和空间二阶精度新型解耦算法进行两方面改进,流动算子采用基于Runge-Kutta方法的时间格式以后, 可以推广到更多的空间差分格式,化学反应源项求解算子可以采用梯形公式、拟稳态逼近法和变系数常微分方程求解器. 对H化学非平衡流动; 解耦算法; 计算方法对一种模拟化学非平衡流动的时间和空间二阶精度新型解耦算法进行两方面改进,流动算子采用基于Runge-Kutta方法的时间格式以后,可以推广到更多的空间差分格式,化学反应源项求解算子可以采用梯形公式、拟稳态逼近法和变系数常微分方程求解器.对H_2/Air预混气体中激波诱导振荡燃烧的Lehr试验进行模拟,考察了化学动力学模型、网格尺寸和差分格式耗散大小对计算结果的影响,同时对不同的化学反应源项算子求解算法的计算效率进行了比较.     

高超声速飞行复现风洞理论与方法

针对高超声速飞行伴随的热化学反应流动,本文回顾了郭永怀先生的科研理念和学科布局,综述了他亲手成立的高温气动团队在高超声速飞行风洞实验模拟理论与方法方面的研究进展.高温气体的迅速产生与迅速应用是一种理想的风洞运行方法,而激波管就是这样一种实验装备.论文首先介绍了激波管技术的基本理论与方程,指出将其用于高超声速流动实验模拟时所具有的独特优势.然后讨论了应用激波风洞复现需要的高超声速飞行状态的可行性、基本方程和需要解决的关键问题.针对这些关键问题,进一步介绍了如何应用爆轰现象研发激波风洞驱动技术的理论,并给出了基于爆轰驱动方法的技术发展和工程应用验证.最后,论文介绍了爆轰驱动激波风洞的界面匹配条件,该条件奠定了长实验时间激波风洞运行基础,是其他驱动方法尝试解决而没能完全解决的难题.高温气动团队关于高超声速飞行复现风洞的理论与技术研究,实现了郭永怀先生的战略规划,成就了国际领先的高超声速热化学反应流动研究平台.      

高超声速飞行复现风洞理论与方法

针对高超声速飞行伴随的热化学反应流动,本文回顾了郭永怀先生的科研理念和学科布局,综述了他亲手成立的高温气动团队在高超声速飞行风洞实验模拟理论与方法方面的研究进展.高温气体的迅速产生与迅速应用是一种理想的风洞运行方法,而激波管就是这样一种实验装备.论文首先介绍了激波管技术的基本理论与方程,指出将其用于高超声速流动实验模拟时所具有的独特优势.然后讨论了应用激波风洞复现需要的高超声速飞行状态的可行性、基本方程和需要解决的关键问题.针对这些关键问题,进一步介绍了如何应用爆轰现象研发激波风洞驱动技术的理论,并给出了基于爆轰驱动方法的技术发展和工程应用验证.最后,论文介绍了爆轰驱动激波风洞的界面匹配条件,该条件奠定了长实验时间激波风洞运行基础,是其他驱动方法尝试解决而没能完全解决的难题.高温气动团队关于高超声速飞行复现风洞的理论与技术研究,实现了郭永怀先生的战略规划,成就了国际领先的高超声速热化学反应流动研究平台.     

自由活塞压缩管ALE方法数值模拟

当前国际上实现高焓气体流动的实验手段之一是自由活塞驱动类脉冲设备,包括自由活塞激波风洞和自由活塞膨胀管.采用自由活塞压缩管作为激波风洞和膨胀管的驱动段时,其驱动能力在很大程度上决定了该类设备的性能.本文采用计算流体力学中任意拉格朗日——欧拉方法(arbitrary Lagrangian Eulerian)数值模拟了压缩管内部的自由活塞运动和气体流动特征.采用移动网格技术来适应活塞运动边界,耦合求解网格运动和气体流动过程,并通过双时间步长方法进行流体运动的时间积分.为了满足几何守恒律(geometric conservation law),对移动网格的法向矢量和表面面积计算进行了修正.不同时刻的活塞位置试验测量结果及欧拉方法预测结果,以及基于简单波理论获得的运动活塞底部气体压力、活塞速度与活塞位置都与当前的ALE方法十分一致.该工作为下一步数值模拟自由活塞激波风洞和自由活塞膨胀管中包括压缩管、激波管和喷管等不同部位的耦合流动提供了基础.      

激波风洞高低压段钢膜片破裂特性研究

激波风洞是用于高超声速飞行器气动外形设计和优化的常用地面试验装置,基于爆轰驱动技术,激波风洞能够在短时间(毫秒级)内产生高温、高压的驱动气体来模拟高超声速试验气流.主膜片位于激波风洞中的爆轰驱动段和激波管段之间,试验时膜片在爆轰脉冲压力下打开,膜片的打开状态和脱落情况对激波风洞气流品质有很大的影响. 同时,膜片也是形成激波的先决条件. 传统的风洞采用铝质膜片进行试验,在激波风洞中需要承压能力更强的膜片, 此时铝质膜片不再适用, 需要采用钢质膜片.因此, 对激波风洞中的钢膜片破裂特性进行研究很有必要.将数值计算结果与试验结果进行比较, 发现数值计算结果与试验结果吻合得比较理想,计算结果具有可靠性. 基于膜片的应力-应变模型, 建立了膜片打开的动力学模型,根据CJ爆轰理论, 采用有限元软件计算模拟了膜片破裂的过程,分析总结了膜片破裂的机制和力学特性规律.采用控制变量法对不同厚度和凹槽长度的膜片进行分析研究,得到了膜片破膜压力和有效破膜时间的变化规律. 在激波风洞试验中,根据膜片总破膜时间设计了适用于JF-12复现风洞的膜片参数.      

喷管内水蒸气非平衡凝结流动数值模拟研究

基于水蒸气非平衡凝结动力学模型和水蒸气真实热力学性质模型,建立了水蒸气跨声速非平衡凝结流动数值模型,对喷管内水蒸气非平衡凝结流动进行了数值模拟,模型成功地捕捉到了凝结激波,并且模拟结果与实验数据吻合较好,证明了所建立数值模型的准确性和可靠性。在此基础上研究了水蒸气在跨声速流动过程中发生的非平衡凝结现象以及对流场产生的影响,发现水蒸气凝结现象对流场的影响非常显著;并对凝结激波的产生过程和原因进行了详细分析。     

基于立方型气体状态方程的激波风洞准一维流动数值研究

采用基于立方型气体状态方程的准一维流动数值模拟方法研究了反射式高焓激波风洞的真实气体流动,重点关注了高压真实气体效应对风洞全场流动时空结构和驻室区气流参数的影响,并以理论分析揭示了高压真实气体效应对激波管内流动的作用机理。研究表明：对于以冷高压气体驱动的激波风洞,使用考虑分子体积和分子间作用力的真实气体状态方程能够更准确地描述气体的状态和风洞内的流动状况。高压真实气体效应主要在冷驱动气体中发生作用,其作用效果主要是使当地声速增大,从而使得入射稀疏波和反射稀疏波的传播速度加快;另一方面,高压气体效应在高温气体效应较显著的被驱动气体中作用微弱,且对激波管产生激波的强度和激波后的流动状态影响甚微。稀疏波的加快传播改变了激波管波系的相干时空关系。提前抵达的稀疏波可在一定情况下侵蚀激波风洞的有效试验时间。对于所测试的激波风洞构型,在150 MPa氢气驱动110 k Pa氮气的工况下,高压效应导致的有效试验时间缩短约38%。适当加长驱动段长度和采用高温气体驱动均可有效减弱高压真实气体效应的影响。     

用于研究双波迎面相互作用的电控双驱动激波管实验技术

本文介绍了用于研究平面运动激波与头激波迎面相互作用的电控双驱动激波管技术。研究成功的用电容高压放电破膜(聚脂薄膜)技术,破膜时间可控制在100微秒左右,破膜重复性为±50微秒;提出了同时测量和记录管内前后两道运动激波的速度,以及两激波间时间间隔的方法;并以第二激波作为纹影照相的同步信弓,实现了±2微秒的同步精度,从而用单火花和双火花纹影系统获得了双波迎面相互作用的流场照片。     

含化学反应膛口流场的无网格数值模拟

基于无网格方法,对包含大位移运动边界和非平衡化学反应的膛口流场进行了数值模拟。所发展算法是基于线性基函数最小二乘显式无网格方法,忽略黏性及湍流的影响,对流场采用ALE(arbitrary Lagrangian-Eulerian)形式的Euler方程描述,对流通量和化学反应源项采用多组分HLLC(Harten-Lax-van Leer-Contact)格式和有限速率反应模型计算,对于运动边界造成的点云畸形采用局部点云重构方法处理,重构过程中采用虚拟边阵面推进。对圆柱绕流和激波诱导燃烧流场进行了数值模拟,验证了重构方法和化学反应计算的有效性。最后对12.7 mm口径机枪膛口流场进行了模拟,结果同实验照片、非结构网格方法结果吻合较好,数值结果清晰地再现了膛口初始冲击波、膛口冲击波、欠膨胀射流波系结构的动力学发展过程,以及膛口焰的时间、空间分布特征。     

高超声速三维化学非平衡流动简化NS方程数值解

本文研究用时间稳定法求解高超声速三维化学非平衡流动简化ＮＳ方程，把整个非平衡流体动力学问题分裂为流体力学问题和化学问题。流体力学与化学问题之间的耦合通过迭代方式加以处理。考虑高温空气分解与电离，七个化学组元，并认为化学反应是以有限速率进行的．差分格式采用无自由参数的混合反扩散格式，激波及壁面处理应用特征线理论．对来流攻角为１０°，马赫数为２６，飞行高度为８５ｋｍ的钝锥体，给出了流场参数和物面加热率的数值计算结果．     

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

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

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

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