台阶式溢洪道滑行水流水面线和消能效果的试验研究

通过模型试验研究了台阶式溢洪道滑行水流的流态、水面线和消能效果。试验得出,台阶式溢洪道的水流流态分为台阶内部的旋滚流态和虚拟底板以上主流区的流态。台阶式溢洪道水面线的变化与来流量、台阶尺寸和掺气浓度有关。当水流进入台阶段时,水面由于失重有一定的降低,然后逐渐抬高至某一高度后再略有降低;在掺气发生点以前,水面线沿程降低,掺气发生点以后,水面线沿程升高;当水流掺气饱和时,水深沿程不再变化,为明渠均匀流。台阶式溢洪道的消能率是逐级台阶的累积效应。水流通过与台阶之间的碰撞、台阶上的水流漩滚以及水流内部的紊动剪切强化了溢洪道的消能效果。在试验范围内,台阶式溢洪道与光滑溢洪道比较,消能率可以提高40% ~70%。通过试验,得出了台阶式溢洪道滑行水流水面线和消能率的计算方法。     

台阶式溢洪道掺气特性的试验研究

通过试验研究了台阶式溢洪道上的掺气现象以及掺气区域的划分，试验表明，台阶式溢洪道沿程掺气可以分为非掺气区、掺气发展区和掺气充分发展区等三个区域；断面掺气分为台阶影响的含气区、楔形清水区、水中气泡区和空中水滴区。试验还测量了台阶式溢洪道上的断面掺气浓度及其沿程分布，计算了断面含水率，提出了初始掺气点的计算方法。与光滑溢洪道相比较，台阶式溢洪道具有初始掺气点前移，掺气发生早，发展快，掺气浓度大等特点。     

90°锥头弹丸不同速度下垂直入水冲击引起的空泡特性

用高速摄像拍摄了90°锥头弹丸低速入水的空泡形态演变过程,全面讨论了不同入水冲击速度下空泡的闭合方式及其演变过程,分析了空泡闭合时间、闭合点水深和弹头空泡长度随入水速度的变化规律以及不同水深位置空泡直径的变化规律;研究了水幕闭合和近液面空泡收缩上升所形成的射流现象及其相互耦合作用过程,探讨了空泡深闭合后其壁面波动规律。结果表明:随着入水速度的增加,空泡分别发生准静态闭合、浅闭合、深闭合和表面闭合,每种闭合方式对应的一个速度区间;弹头产生空泡的临界入水速度为0.657 m/s;不同水深位置的空泡直径呈现非线性变化;随着水深的增加空泡扩张初速增大,空泡最大直径减小,扩张段缩短,收缩段延长;同一时刻水深越大空泡扩张收缩的加速度也越高;水幕闭合后会产生向上和向下两股射流,向下射流速度较大时会对弹丸运动产生影响;近液面空泡收缩上升时会产生强度正比于空泡体积大小和闭合点水深的射流,并与上两股射流相互耦合形成一股更强的向上射流;空泡深闭合后长度缩短和产生的向下射流使弹丸受力发生改变,弹丸速度因受力产生的变化带动了流体质点速度的波动,进而导致空泡壁面发生波动,壁面波动遵循空泡截面独立扩张原理。     

用格子Boltzmann方法模拟液滴撞击固壁动力学行为

首次用格子Boltzmann方法中的伪势模型对液滴撞击固壁的动力学行为进行了数值模拟.详细研究了液滴在壁面上的流动状态以及各种因素对撞击过程的影响.通过数值模拟得到:壁面的可润湿性越小,液滴越容易发生反弹,液滴的回缩速度越快;液滴的撞击速度越大,所得到的相对直径越大,回缩速度越快;液滴的粘性越小,所得到的相对直径越大;液滴的表面张力越大,液滴越容易发生反弹现象.另外,液滴的最大相对直径与We数满足一定的线性关系,这些结果与前人的理论预测和实验结果完全吻合.     

三爆轰管-脉冲爆轰发动机的近场噪声特性

为研究三爆轰管脉冲爆轰发动机噪声的形成机理和特性,设计了正三角形组合的三爆轰管PDE实验系统,分别在与中心轴线成0°、30°、60°及90°方向上,对8个不同距离处的声压进行了测量,结果表明:不同距离处噪声幅值最大值均在30°方向上;爆轰噪声的参考半径为3倍“名义管径”;爆轰噪声的A持续时间随着r(距离三爆轰管名义“中点”的直线距离)的增加近似线性减小,并随着角度的增大而减小。B持续时间与声压值的大小成反比关系,并且随着轴向距离的增加而增加。在不同的方位角上,峰值声压越大,B持续时间越小。     

不同交汇角度明渠交汇口三维水力特性的大涡模拟研究

采用气液两相流混合模型对不同交汇角度下等宽明渠交汇口三维水力特性进行了数值模拟研究。选取大涡模型(Large Eddy Model)封闭两相流时均方程,求解速度与压力耦合方程组时使用半隐式SIMPLE(Semi-Implicit Method for Pressure-Linked Equations)算法,模拟自由水面采用VOF(Volume of Fluid)法。将交汇角度为90°时采用大涡模型计算得到的纵向截面水面线和不同测线上的速度分布与文献中的试验结果相比较,两者吻合良好,且水面线高度误差在4.2%以内,由此可见大涡模型是模拟交汇口水力特性的有效方法。进而将大涡模型用于模拟交汇角度为30°、45°、60°的交汇口水流,得到交汇口处的水深变化及流场的分布规律,并定量分析了交汇口下游各横断面流速不均匀系数的分布规律。结果表明:在整体上交汇角度越大,交汇口各特征横断面流速不均匀系数越大,即水流流速分布越不均匀。     

k-ω-γ模式对转捩影响因素的预测性能研究

高超声速边界层转捩的准确预测对飞行器的防热、减阻至关重要,而影响高超声速边界层转捩的因素众多.从模式角度出发研究边界层转捩的影响因素,采用k-ω-γ 转捩模式对5°圆锥的边界层转捩进行了数值分析,计算了不同头部钝度、来流雷诺数和湍流度情况下的边界层转捩,并与实验结果进行了对比. 研究结果表明:k-ω-γ 转捩模式基本能够反映头部钝度、来流雷诺数、来流湍流度对高超声速圆锥边界层转捩的影响规律,但对转捩后的热流峰值预测不准;从模式构造角度分析发现,雷诺数越高或头部钝度越小,层流区边界层越薄,k-ω-γ 转捩模式中第一、第二模态时间尺度增大,因此转捩起始位置提前;来流湍流度越大,等效脉动动能初值越大,导致层流区发展过程中等效脉动动能越大,因此转捩易于发生.      

考虑多重运移机制耦合页岩气藏压裂水平井数值模拟

页岩作为典型的微纳尺度多孔介质,游离气与吸附气共存,传统的达西定律已无法准确描述气体在页岩微纳尺度的运移规律.基于双重介质模型和离散裂缝模型构建页岩气藏分段压裂水平井模型,其中基岩中考虑气体的黏性流、Knudsen 扩散以及气体在基岩孔隙表面的吸附解吸,吸附采用Langmuir等温吸附方程;裂缝中考虑黏性流和Knudsen扩散,在此基础上建立基岩-裂缝双重介质压裂水平井数学模型并采用有限元方法对模型进行求解.结果表明,基岩固有渗透率越小,表面扩散和Knudsen扩散的影响越大,反之则越小;人工裂缝的性质包括条数、开度、半长以及间距,主要影响压裂水平井生产早期,随着人工裂缝参数值的增加,压裂水平井产能增加,累产气量也越大.其次,页岩气藏压裂诱导缝和天然裂缝的发育程度对页岩气藏的产能有很大的影响,水平井周围只有人工裂缝,周围天然裂缝不开启或不发育时,页岩气藏的水平井的产能较低.      

金属丝网橡胶压缩力学性能研究

对金属丝网橡胶进行了静态压缩试验。利用控制变量法研究了压缩量、相对密度、金属丝丝径、丝材和承压面积对金属丝网橡胶压缩力学性能的影响,并对平均刚度和能量耗散系数随压缩量和相对密度的变化关系进行了研究。试验结果表明:随着压缩量的增加,金属丝网橡胶非线性力学特性逐渐增强;相对密度越大,金属丝网橡胶承压能力越强;金属丝的丝径和丝材主要影响金属丝网橡胶非线性阶段的力学特性,丝径越大,丝材越硬,承压能力越强;承压面积越大,金属丝网橡胶的承压和耗能性能越好;随着压缩量的增加,平均刚度增大,承压能力增强,能量耗散系数减小,减震性能降低;随着相对密度的增加,平均刚度和能量耗散系数均增大,承压能力和减震性能均增强。     

波浪作用下直立结构物附近强湍动掺气流体运动的数值模拟

波浪破碎卷入气体易对建筑物受力产生压力振荡, 了解波浪作用下建筑物附近掺气水流的运动特性是精确计算建筑物受力的前提. 基于OpenFOAM开源程序包和修正速度入口造波方法建立三维数值波浪水槽, 模型采用S-A IDDES湍流模型进行湍流封闭, 并采用修正的VOF 方法捕捉自由液面, 数值模拟了规则波在1:10的光滑斜坡上与直立结构物的相互作用过程, 重点分析了结构物附近的水动力和掺气水流运动特性. 结果表明, 建立的数值模型能精确地捕捉波浪作用下直立结构物附近的自由液面的变化以及气泡输运过程, 较好地描述气体卷入所形成的气腔形态以及多气腔之间的融合、分裂等过程; 波浪与直立结构物相互作用产生强湍动掺气水流, 其运动过程十分复杂; 掺气流体输运过程中水气界面周围一直伴随着涡的存在, 其中, 气泡的分裂与周围正负涡量剪切作用密切相关, 且其输运轨迹主要受周围流场的影响; 研究揭示了结构物附近湍动能与掺气特性的关系, 发现波浪作用下直立结构物附近湍动能的分布与掺气水流特征参数(气泡数量、空隙率)整体呈现一定的线性关系.      

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

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

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

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