一种爆破片气液两相流临界泄放计算方法

针对爆破片气液两相流超压泄放难以确定的问题,应用分相模型和流体基本方程,在绝热等熵假设下引入滑速比理论,推导出气液两相流临界泄放量的计算模型。在环境压强为0.2MPa~0.6MPa情况下,计算得到:当含气率从0增加到0.1时,质量流量急剧下降;当含气率从0.1增加到0.96时,质量流量下降比较平缓;工况压力越大,质量流量越大。其中在含气率大于0.1时,计算结果与实验结果基本吻合;含气率小于0.1时,随着含气率趋近0,误差越来越大。含气率小于0.1时,通过与文献中Yoon模型、Fauske模型、实验结果进行比较,验证了绝热等熵模型仍然具有较高的准确性。最后将其与GB567-2012方法、等焓膨胀模型、均相计算模型比较,计算结果基本一致,验证了模型实际应用的准确性。     

液态金属熔池中气泡-液体两相湍流的数值模拟研究

应用欧拉-欧拉双流体模型,以及分别描述气泡和液体湍流的两相湍流模型,描述液态金属中气体射流发泡过程的两相流动及气泡分布。与实验结果的对比表明,本模型的预测能力优于前人未考虑气相湍流的模型。考查了液体物性和气泡尺寸对流场特性的影响,结果发现,液体密度越低,气泡尺寸越小,气泡分布就越均匀。     

微重力气液两相流动与池沸腾传热

综述了近年来中国科学院微重力重点实验室(国家微重力实验室)完成的一系列微重力气液两相流动与池沸腾传热方面的地基实验、飞行实验和理论研究等方面获得的主要成果.在微重力气液两相流动方面,提出了半理论Weber数模型用于预测微重力条件下气液两相弹-环状流转换,并采用Monte Carlo方法,针对气泡初始尺寸对泡-弹状流转换的影响进行数值研究.通过俄罗斯"和平号"空间站与IL-76失重飞机实验,获得了微重力下的气液两相流型图,与此同时在地面利用小尺度毛细管模型模拟了微重力气液两相流动特征.实验测量了微重力气液两相流压降,并基于微重力流动特性建立了一个泡状流压降关联模型.在微重力池沸腾传热方面,利用我国返回式卫星完成了两次空间实验,其中,第22颗返回式卫星搭载铂丝表面R113池沸腾实验采用控制温度的稳态加热方式,而实践8号育种卫星搭载平面FC-72池沸腾实验则采用控制加热电压的准稳态加热方式.同时,还进行了地面常重力和落塔短时微重力条件下的对比实验研究.观察到丝状加热表面微重力时轻微的传热强化现象,而平板加热表面微重力核态池沸腾低热流时传热强化、高热流时传热恶化.微重力实验中观察到气泡脱落前存在横向运动现象,据此分析了气泡行为与传热之间关系,并提出了一个预测丝状加热表面气泡脱落直径的半理论模型.旨在对相关领域的进一步发展和空间两相流系统的应用提供数据及理论支持.     

气液两相流压力波传播速度研究

将双流体模型用于绝热无相的管道气液两相流,依据小扰动线化分析原理,导出了压力波波数Ｋ方程通过对不同空隙率下肉体上压力波小随角频率变化的计算,研究了虚拟质量力和狭义相间阻力对压力波波速及其人色散性的影响。对泡状流和弹状流压力波波速的计算结果与前人的测量结果作了比较,两者符合良好。     

通气空泡尾部气泡流仿真和特性分析

基于欧拉-欧拉双流体模型对通气空泡尾部气泡流进行数值仿真, 并采用基于多尺寸分组模型的总体平衡方法预估气泡尺寸分布. 应用改进后湍流耗散系数计算模型, 考虑了气泡体积含量对在湍流作用下气泡扩散现象的影响. 基于上述模型对两种试验工况下流场进行了数值仿真.结果表明模型对空泡尾部回流区特性进行了准确预示, 在回流区高湍流度作用下气泡迅速破碎成小气泡. 并进一步得到试验体尾流区空泡体积分数和速度分布. 尾流区水流速度分布保持了流体经过非流线型对称体时产生的尾流分布规律. 仿真结果与试验数据相一致, 模型适用性得到验证.      

垂直上升圆管内等温泡状流的数值模拟

介绍了描述等温状态下垂直气泡流流动的总体平衡方法(population balanceapproach). 将平均气泡数密度控制方程引入到双流体模型中实现总体平衡. 介绍了描述气泡聚并与破裂机制的Yao和Morel理论模型以及Hibiki等人的实验方法. 利用商业计算流体力学软件ANSYSCFX10, 对圆管特定位置上的5个基本变量:气含率、气泡平均直径、相间表面积浓度以及气体和液体速度沿半径分布进行了模拟计算. 预测结果与实验数据的比较表明,预测结果和实测数据之间有很好的一致性.      

气液两相涡街的数值计算

本文利用有限差分法，建立了模拟气液两相涡街的物理模型和数值计算方法，成功地利用数值方法模拟出了气液两相涡街及涡街中含气率分布等主要特征，该方法简便，理论计算出的结果与实验值相符较好     

垂直向下管内两相流泡状-弹状流型转换研究

建立实验系统,在维持管道出口压力为0.2MPa的条件下,对内径分别为15mm、25mm、40mm、65mm的垂直向下管内空气-水气液两相流动进行了实验研究,获得了两相流泡状-弹状流型分布。实验研究发现:管径对于泡状流与弹状流流型特征有较大影响,并且进一步影响流型转换边界,随着管径增加,泡状流-弹状流的流型转换边界向折算气速减小的方向移动。基于理论推导及实验数据,建立了垂直向下管内气液两相流泡状流-弹状流流型转换预测模型,该模型对本文实验工况条件下的垂直向下管内空气-水气液两相流流型转换具有良好的预测效果,预测模型的计算结果与实验数据之间的误差小于10%。     

气液两相流中旋涡诱发圆柱振动时的脉动升力研究

一定条件下 ,垂直上升矩形截面管内的气液两相流横向冲刷水平布置的圆柱时 ,圆柱在与来流垂直的水平方向上受到了脉动升力的作用。本文在 1.6× 10 4到 6 .0×10 4的 Re数范围和 0到 0 .32的来流含气率范围内 ,进行了实验研究。实验中 ,流体工质为由空气和水混合而成的泡状流 ,采用弹性梁 -电阻应变式小力值传感器来测量圆柱受到的脉动升力 ,利用计算机对传感器的输出信号进行连续采集 ,采集频率为5 0 0 Hz。对由实验得到的脉动升力数据进行自功率谱分析 ,得出在两相流中只有在小于 0 .10～ 0 .12的来流含气率范围内 ,圆柱后部才会形成涡街 ,且此时的 Strouhal数随来流含气率的增大而增大 ;由脉动升力数据的均匀方根值算出脉动升力系数 C′L发现小 Re数时 C′L 随来流含气率的增大而增大 ,而在大 Re时 C′L 随来流含气率的增大先稍下降后再增大。     

不同重力条件下管内冷凝现象研究进展

本文对不同重力环境中水平管内冷凝气液两相流动与传热现象的研究现状进行了全面评介, 重点关注于航天应用中的小管径、低流量和以氨为工质的情形,以及部分重力条件下冷凝现象中的重力效应. 冷凝现象中管壁四周液膜的存在,导致水平管内冷凝气液两相流与沸腾或绝热气液两相流在构型上存在着明显的不同, 表面张力的作用增大,向波状和分层流转换所对应的临界Bond数也随之增大, 进而使小管径、低流量水平管内冷凝气液两相流动与传热现象中的重力效应减弱,甚至可以忽略不计. 实验发现该情形中摩擦压降要小于基于常规气液两相流实验数据的经验关联式的预测结果, 而更接近光滑环状流情形; 对冷凝两相传热系数的预测,只有采用基于空隙率与界面剪切率相互耦合的模型, 才能够给出较好的结果.      

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

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

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

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