HPAM溶液在石英微圆管中的流动实验

针对HPAM溶液在油藏孔隙中的复杂渗流特性问题,在内径为10～350μm微圆管中进行了部分水解聚丙烯酰胺(HPAM)溶液的流动实验.实验结果表明:在本文实验条件下,HAPM溶液在管径尺度较小的微圆管内的流动规律明显偏离常规尺度下的非牛顿流体力学流动规律,其截面平均流速高于按照常规尺度流动推算出的结果,且管径越小,偏离程度越大,微尺度效应越强.     

微尺度气体流动

了解微尺度气体流动特点是微机电系统设计和优化的基础.有关的 研究可以上溯到20世纪初Knudsen的平面槽道流动质量流量的测量和 Millikan的小球阻力系数的测量，实验结果揭示了稀薄气体效应即尺度效 应对气体运动的重要影响.由于流动特征长度很小，微尺度气流经常处于 滑流区甚至过渡领域，流动的相似参数为Knudsen数和Mach数.因此可以考 虑利用相似准则，通过增大几何尺寸、减小压力的途径，解决微机电系统 实验观测遇到的困难.为解决直接模拟Monte Carlo方法分析微机电系统中 低速稀薄气流遇到的统计涨落困难，我们提出了信息保存法(IP)，该方法 能够有效克服统计散布，并已成功用于多种微尺度气流.     

HPAM稀溶液在微圆管中流动特性的实验研究

实验研究了部分水解聚丙烯酰胺(HPAM)稀溶液(水溶液)在内径 为10.1 ~325 \mu m石英微管中的高剪切速率(1 215 ~23 120 s ^{ - 1} )流动. 结果表明, 聚合物溶液在管径小于100 \mu m微管中的流动具有 明显的微尺度效应, 实测流速高于由同种溶液在常规管径圆管中的流动 规律所预测的流速. 实测流速与预测流速之间的偏离程度与管径和剪切速率有关: 相同剪切 速率下, 管径越小, 偏离越明显; 在小管径微管中, 偏离随剪切速率的增加而减小, 而当管 径超过30.7 \mu m后, 偏离不再随剪切速率的变化而改变.     

页岩气和致密砂岩气藏微裂缝气体传输特性

页岩气和致密砂岩气藏发育微裂缝,其开度多在纳米级和微米级尺度且变化大,因此微裂缝气体传输机理异常复杂.本文基于滑脱流动和努森扩散模型,分别以分子之间碰撞频率和分子与壁面碰撞频率占总碰撞频率的比值作为滑脱流动和努森扩散的权重系数,耦合这两种传输机理,建立了微裂缝气体传输模型. 该模型考虑微裂缝形状和尺度对气体传输的影响. 模型可靠性用分子模拟数据验证.结果表明:(1)模型能够合理描述微裂缝中所有气体传输机理,包括连续流动,滑脱流动和过渡流动;(2)模型能够描述不同开发阶段,微裂缝中各气体传输机理对传输贡献的逐渐变化过程;(3)微裂缝形状和尺度影响气体传输,相同开度且宽度越大的微裂缝,气体传输能力越强,且在高压和微裂缝大开度的情况下表现更明显.     

CO2微气泡溶解动力学及提高采收率机理研究

CO₂微气泡是一种具有潜力的提高采收率与碳埋存方法,本文在自主设计的CO₂微气泡发泡装置的基础上,表征了高温高压条件下微气泡形态,进一步研究了微气泡的溶解特征,研究结果表明:10 MPa下制备出的微气泡直径10～70μm,平均直径34.43μm; 15 MPa下制备的微气泡直径更小,平均直径25.03μm;地层水高矿化度条件下,平均气泡直径277.17μm,且气泡稳定性降低.微气泡的溶解实验结果表明CO₂微气泡的溶解速率较高,但是未溶解的CO₂仍以气泡的形式在地层中运移,微气泡注入地层后将形成“碳化水+微气泡”的运移模式.采用可视化微流控平台,首次研究了高温高压条件下无化学剂辅助CO₂微气泡的提高采收率机理:(1)提高微观洗油效率;(2)通过体积膨胀、溶解携带作用将油滴带出盲端,采出盲端中的剩余油;(3)打破油滴的毛管压力平衡状态,采出柱状残余油;(4)在流动中产生“贾敏效应”,封堵大孔隙、提高波及效率.本文研究可为CO₂微气泡提高油藏采收率与碳封存提供指...     

复杂微通道气体流动的格子Boltzmann模拟

构建了一个模拟复杂微通道内气体流动的多松弛格子Boltzmann模型。该模型采用动力学曲面滑移边界,考虑了微尺度效应和努森层影响。此外,为了更准确地描述微通道内气体的滑移速度,在模型中引入孔隙局部Kn数来代替平均Kn数。之后采用Poiseuille流对模型进行验证,模拟结果与用直接模拟蒙特卡洛方法和分子模拟结果吻合较好,证明了该模型模拟微通道内处于滑移区和过渡区气体流动的有效性。最后,采用该模型模拟多孔介质内气体渗流过程。结果表明,随着孔隙平均Kn数的增加,多孔介质内的高渗区域增加,且优先从小孔隙中开始增加,这是由于小孔隙中微尺度效应更加明显,相对大孔隙流动阻力更小所致。     

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模型,当入射激波马赫数较低时,扰动增长率受微尺度效应的影响更为明显。     

基于微尺度造型的干气密封流动有序性数值分析

基于干气密封微尺度流动特性,提出一种有序微造型的干气密封模型,可在提升密封性能的同时为激光开槽提供新的借鉴和参考.选择螺旋槽和T型槽两种干气密封经典槽型为研究对象进行仿真分析,通过文献对比验证了计算方法的正确性,在有序微造型的基础上进行了有无微造型性能分析和对比,最后对微尺度下微造型的密封性能开展了系统研究.结果表明:同工况下,微造型结构的开启力较传统结构有明显提升,在高速高压及微尺度时的提升量愈加显著;T槽型的对称性使得其微造型结构兼具一定的减漏效果,且随膜厚增大减漏效果越明显;微造型深度、数量和面积对密封性能影响较大,存在一个使开启力较大同时泄漏量较小的最优槽深(螺旋槽和T型槽分别为5.5和2.5μm)及微造型深度区间(螺旋槽和T型槽皆为0.9～1.2μm).具微造型结构良好的增压性能,有助于进一步提高干气密封的稳定性.     

考虑稠密气体效应的微尺度气体轴承模拟

当气体流动区的局部压力升高到一定程度时,气体稠密效应对其输运特性的影响不可忽略,其影响可以用Enskog理论来描述。利用直接模拟蒙特卡罗(DSMC)方法以及考虑气体分子稠密效应的衍生方法对气体粘度随特征尺度的变化规律和气体粘度随密度变化的规律进行了计算分析,在此基础上研究了考虑稠密效应的微尺度气体轴承的承载能力。计算结果表明,稠密效应在抵消掉气体稀薄效应的同时还改变了输运系数和气体薄膜的润滑特性,提高了气体轴承的承载能力。     

Navier-Stokes方程二阶速度滑移边界条件的检验

对微尺度气体流动,Navier-Stokes方程和一阶速度滑移边界条件的结果与实验数据相比,在滑移区相互符合,在过渡领域则显著偏离.为改善Navier-Stokes方程在过渡领域的表现,有些研究者尝试引入二阶速度滑移边界条件,如Cercignani模型,Deissler模型和Beskok-Karniadakis模型.以微槽道气体流动为例,将Navier-Stokes方程在不同的二阶速度滑移模型下的结果与动理论的直接模拟Monte Carlo(DSMC)方法和信息保存(IP)方法以及实验数据进行比较.在所考察的3种具有代表性的二阶速度滑移模型中,Cercignani模型表现最好,其所给出的质量流率在Knudsen数为0.4时仍与DSMC和IP结果相符;然而,细致比较表明,Cercignani模型给出的物面滑移速度及其附近的速度分布在滑流区和过渡领域的分界处(Kn=0.1)已明显偏离DSMC和IP的结果.     

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

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

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

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