泥石流研究述评

泥石流这一特殊的二相流体是山区常见的一种自然灾害,一般可分为稀性泥石流(颗粒碰撞及离散力起主导作用)和粘性泥石流(粘性力起主导作用)两大类。本文概略地评述了泥石流的几个问题,即泥石流形成学、分类学、动力学等,指出急需解决的问题有泥石流汇流理论、泥石流阻力规律等。文中还根据近年来的动态和趋向,推断更加深入的理论研究和以数学模型为核心的泥石流流域管理模式将会成为泥石流学科发展的主要方向。      

一维定常泥石流的数学模型

根据泥石流在山坡斜面的长度（Ｌ）远大于泥石流横截面的特征长度（Ｈ＿０），对连续方程，运动方程和Ｂｉｎｇｈａｍ模型的本构方程进行量级比较，忽略Ｈ＿０／Ｌ量级的量，推导出控制一维定常泥石流运动的微分方程。本文方法可应用于推导非定常一维泥石流的数学模型。     

泥石流预测、预报技术方法的研究现状浅析

泥石流灾害的预测、预报是减灾防灾重要环节和基础。在充分研究泥石流形成条件基础上 ,正确预测、预报泥石流 ,是目前最大限度减少泥石流灾害损失所采取的主要方法。本文综合了当前泥石流预测预报的方法 ,分析了泥石流长期、中期、短期预报及泥石流预报的原理 ,不同时期泥石流预报的特征及存在的问题。我国是受泥石流危害最为严重的国家之一。每年由泥石流造成的直接经济损失约达几十亿元之多 ,减少泥石流灾害造成的损失主要取决于预测预报工作等。     

山区铁路建设中的泥石流灾害与防治对策

在分析已建山区铁路泥石流灾害的基础上 ,指出我国中西部多数在建和规划铁路新线都将穿过泥石流强烈活动区 ,存在着少数泥石流沟难以避饶、部分低频泥石流沟难以判定、青藏高原及周边地区大型冰川和融雪泥石流可能造成铁路毁灭性破坏等重大技术问题 ,提出了改进铁路地质勘察方法、利用高新技术和组织专家对铁路建设中泥石流问题进行专门研究等对策。     

四川凉山州美姑县“6·1”泥石流灾害研究

四川凉山州美姑县“6.1”泥石流灾害实例研究表明,该泥石流约为20年一遇的中小规模的泥石流。流域上游短历时强降雨和冰雹天气过程是这次泥石流暴发的诱因,流域内退化的生态环境和中下两岸不稳定边坡以及沟道内大量的松散堆积物为这次泥石流提供了丰富的固体物质来源。泥石流堆积物具有典型的多峰型粒度特征,且有较高的粘粒含量。巨大的泥石流漂砾、石背石现象、龟裂现象、侧积堤和龙头堆积证实了这次泥石流为粘性泥石流。危险度评价表明,采莫洛沟属于高度危险的泥石流沟,危险度为0.67;乃托沟属于中度危险的泥石流沟,危险度为0.58。风险评估结果可知,两沟都属于泥石流高风险区风险度分别为0.52和0.45。高风险区的泥石流灾害给当地的经济社会造成了严重影响并直接造成了较大的人员伤亡和财产损失。     

四川泸定县泥石流灾害成因、特征与防治建议

泸定县地处川西高原与四川盆地的过渡地带,地质环境条件复杂,分布泥石流沟127条,发育密度054条·10km-2。在遥感解译和地面调查基础上,论述了泸定县泥石流灾害成因、发生特征、危害和发展趋势,提出了相关防治建议。(1)成因表现为复杂的断裂构造导致流域内松散物源丰富,较大的相对高差提供了良好势能条件,而汛期集中强降雨则提供了动能条件,且雨量条件中当日雨量与前期有效降雨量之间具有幂函数关系; (2)发生特征表现出典型的空间群发性(包括2005-06-30群发性泥石流、2005-08-11群发性泥石流和2006-07-14群发性泥石流)与时间的夜发性; (3)危害方式主要包括淤埋、冲毁、侵蚀、堵河4种,泥石流已造成74人死亡; (4)发展趋势表现为泥石流沟多处于形成期和发展期(113条,占89%),发生的年际准周期由20a缩短至不足10a,且发生条(次)由4条(次)增加至44条次(2005年); (5)泸定县泥石流防治的突出问题表现为城镇泥石流问题、泥石流堵塞主河问题、风景区泥石流问题与冲沟泥石流问题,建议通过科学规划、风险防御体系构建、土木工程治理、监测预警等措施予以防御。     

泥石流本构模型及动力学模拟研究现状综述

由于泥石流动力学过程的复杂性,其本构模型及数值模拟研究仍然存在众多的难点问题。本文根据大量的文献资料和实例研究,对目前泥石流本构模型和动力学模拟研究现状和存在的问题进行了综合论述,主要包括:(1)已有泥石流动力学本构模型原理及适用性问题;(2)已有数值模拟方法的实现及选择应用;(3)泥石流动力学复杂环境效应研究,复杂地形和水文环境条件下的泥石流应力分布、动能传递等。最后对泥石流动力学G IS模拟的前景进行了讨论。     

金沙江美姑河牛牛坝水电站库区泥石流对工程影响分析

金沙江美姑河牛牛坝水电站库区泥石流沟分布面积广、发生频率高;调查表明库区现有不同类型泥石流沟31条,其中属于高度危险的泥石流沟4条,中度危险的泥石流沟15条;这些泥石流不会造成严重的堵河问题。在施工期泥石流对水电站工程的影响突出,特别是靠近库首的泥石流对工程的安全构成威胁。水库蓄水后,库区泥石流对水电站工程影响有所降低,但位于大坝下游区的泥石流对水电站正常运行仍有较大的影响。     

基于结构两相阻力分析的泥石流堆积形态研究

针对泥石流的结构两相阻力特征,采用结构两相阻力对泥石流的堆积形态进行了静力学分析,推导出在一般条件下泥石流堆积形态的公式,并讨论了在不同条件下泥石流堆积形态的特点。分析表明采用结构两相阻力对泥石流堆积形态分析可以较合理的模化阻力项,能够表现泥石流体复杂的成分结构组成,较真实的反映实际泥石流的堆积形态特征。在特定材料条件下,结构两相阻力项转化为单一材料阻力项,所得到的相应结论与已有对单一材料的分析结论一致。     

汶川地震后四川省绵竹市清平乡文家沟泥石流灾害调查研究

在2008年5月12日汶川地震后的地震灾区暴发了许多泥石流灾害,其中以四川省绵竹市清平乡文家沟泥石流最为显著。文家沟原来不是泥石流沟,在汶川地震时由于滑坡形成的巨大的滑坡-碎屑流堆积体改变了文家沟的泥石流形成条件,在此后的3个雨季内,文家沟先后暴发了5次大规模和特大规模的泥石流灾害,其中以8 ·13文家沟泥石流规模和危害最大。8 ·13文家沟泥石流暴发时的总降雨量为227mm,泥石流持续时间约2.5h,泥石流总量约310×10⁴m³; 泥石流造成7人死亡, 5人失踪, 39人受伤, 479户农房被掩埋,直接经济损失4.3亿元。5次大规模和特大规模的泥石流以及洪水仅带走了16%的可以很容易形成泥石流的滑坡-碎屑流堆积物,文家沟如再遭遇较大降雨还会暴发泥石流。即使在今后的雨季中暴发几次规模如8 ·13泥石流一样大的特大规模泥石流,文家沟在较大降雨下仍然可能暴发泥石流灾害,因此对文家沟泥石流的防治工作将是一个长期的工作。     

The warp and stress transmission in the joint of an open thin wall beam

The regularities of the warp in the sections of an open thin-wall beam have been discussed using the MODERN ENGINEERING THEORY OF BEAM. Thus, the transmission relationship between the warp displacements and the double-moments in the joints of open thin-wall beams which are connected in different ways has been deduced. Meanwhile, experiments conducted by the author and some conclusions reached from these experiments as well as the tentative idea for further perfection are briefly presented here.Communicated by Chou Li.     

椭圆颗粒在剪切流中旋转特性的数值研究

研究颗粒在流体剪切作用下的运动特性是理解和预测颗粒悬浮流流动行为的关键.当流体的惯性不能忽略时,颗粒的运动往往变得非常复杂.本文采用格子Boltzmann方法对中等雷诺数下椭圆颗粒在剪切流中的旋转运动进行了模拟.首先,研究了雷诺数(0Re 170)的影响,结果表明当雷诺数低于临界值时,颗粒以周期性的方式旋转,角速度最小时对应的长轴方向随着雷诺数的增大而逐渐远离水平方向,而且这一倾角与雷诺数呈分段线性关系;当雷诺数大于临界值时,椭圆形颗粒最终保持静止状态,且静止时的转角与雷诺数呈幂函数关系,雷诺数越大,转角越小,椭圆的长轴越远离水平位置.其次,研究了椭圆颗粒的长短轴之比α(1α10)的影响,结果表明颗粒旋转的周期与α呈幂函数关系,α越大,颗粒旋转周期越小.此外,当α超过临界值时,颗粒也在水平位置附近保持静止状态,此时的转角与α也呈幂函数关系,α越大,转角越小.研究还发现,当雷诺数较大时椭圆颗粒在旋转过程中会产生过冲现象.     

Three-dimensional flow in axial flow fans of non-free vortex design

Three-dimensional laser Doppler anemometer (LDA) measurements were carried out downstream of isolated axial fan rotors of non-free vortex design in order to investigate the role of radial velocity components in design. The structure of secondary flows due to non-free vortex operation was studied in detail. It is pointed out that the tangential gradient of radial velocity at midspan is nearly in direct proportion with the spanwise gradient of ideal total head rise prescribed in design. Design criteria have been established for the neglect of torsion of stream surface segments inside the blading. A linear relationship was proposed in order to estimate the pitch-averaged radial velocities at the rotor exit. Using this relationship, a proposal has been put forward for taking the radial velocity components into account in non-free vortex design with the assumption of conical stream surfaces through the blading.     

On the self-induced motion of vortex sheets

An accurate numerical scheme has been devised to study the self-induced motion of an infinitely thin, free vortex sheet of finite span in an unbounded, inviscid, incompressible fluid. The new numerical scheme has been tested against two vortex sheet problems for which exact solutions have also been obtained. The agreement between the numerical and exact solutions is excellent. The scheme has been further tested against two more examples for which analytical solutions for small times were available. Here too the agreement is excellent.     

Exact solution for entrance region flow between parallel plates

A method of obtaining an exact solution of the complete Navier-Stokes equations for incompressible flow has been developed to study the development of fluid flow between parallel plates. An analytical solution has been developed in the form of a convergent infinite series and a numerical solution obtained for Reynolds numbers of 10 and 500. These results have been compared with data available elsewhere. The method developed here is expected to have general application     

Experimental studies of interfacial instabilities in multilayer pressure-driven flow of polymeric melts

The interfacial deformation and stability of two-(A-B) as well as three-layer symmetric (A-B-A) and asymmetric (A-B-C) pressure-driven flow of viscoelastic fluids has been investigated. Flow visualization in conjunction with digital image processing has been used to observe and measure the rate of encapsulation and interfacial stability/instability of the flow. Specifically, the encapsulation behavior as well as stability/instability of the interface and the corresponding growth or decay rate of disturbances as a function of various important parameters, namely, number of layers and their arrangement, layer depth ratio, viscosity and elasticity ratio as well as disturbance frequency, have been investigated. Based on these experiments, we have shown that the encapsulation phenomena occurs irrespective of the stability/instability of the interface and in cases when both encapsulation and instability occur simultaneously their coupling leads to highly complex and three-dimensional interfacial wave patterns. Moreover, it has been shown that the simple notion that less viscous fluids encapsulate more viscous fluids is incorrect and depending on the wetting properties of the fluid as well as their first and second normal stresses the reverse could occur. Additionally, in two- and three-layer flows it has been shown that by placing a thin, less viscous layer adjacent to the wall longwave disturbances can be stabilized while short and intermediate wavelength disturbances are stabilized when the more elastic fluid is the majority component. Furthermore, in three-layer flows it has been demonstrated that in the linear instability regime no dynamic interaction between the two interfaces is possible for short and intermediate wavenumber disturbances. However, in the nonlinear stability regime dynamic interactions between interfaces have been observed in this range of disturbance wavenumbers leading to highly chaotic flows. Finally, in the parameter space of this study no subcritical bifurcations were observed while supercritical bifurcations resulting in waves with a pointed front and a gradual tail were observed.     

Numerical study of magnetohydrodynamic laminar flow separation in a channel with smooth expansion

The flow of an electrically conducting incompressible viscous fluid in a plane channel with smooth expansion in the presence of a uniform transverse magnetic field has been analysed. A solution technique for the governing magnetohydrodynamic equations in primitive variable formulation has been developed. A co‐ordinate transformation has been employed to map the infinite irregular domain into a finite regular computational domain. The governing equations are discretized using finite‐difference approximations in staggered grid. Pressure Poisson equation and pressure correction formulae are derived and solved numerically. It is found that with increase in the magnetic field, the size of the flow separation zone diminishes and for sufficiently large magnetic field, the separation zone disappears completely. The peak u‐velocity decreases with increase in the magnetic field. It is also found that the asymmetric flow in a symmetric geometry, which occurs at moderate Reynolds numbers, becomes symmetric with sufficient increase in the transverse magnetic field. Thus, a transverse magnetic field of suitable strength has a stabilizing effect in controlling flow separation, as also in delaying the transition to turbulence. Copyright © 2008 John Wiley & Sons, Ltd.     

Relationship between mass discharge rate and granular temperature of silo flow with variance of outlets

A lower granular temperature zone near the orifice in silos was found by speckle visibility spectroscopy technology. Surprisingly, the values of granular temperature in this area are inversely proportional to the mass discharge rate (MDR). This finding that links the macroscopic MDR with mesoscopic dynamics temperature of the particle system, is obtained by looking at the effect that different outlet sizes of silos on MDR and granular temperature field. Apart from evidencing the relevance of a parameter that has been traditionally overlooked in silos, this relationship supposes a benchmark with which to explore the influence of small orifice sizes on MDR. Among these, we found that the frequency of intermittent clogging and their dependence on the size of the outlet.     

Numerical computation of turbulent flow through a square-sectioned 90° bend

Finite-volume, semi-elliptic computations are reported of the three-dimensional flow around a 90° square-sectioned bend for which detailed laser-Doppler measurements have been reported by Taylor et al.¹ While the standard k- eddy-viscosity model has been used in the main flow region, in place of the usual “wall functiond”, the mixing-length hypothesis has been employed to resolve the flow in the layer immediately adjacent to the wall. The scheme is successful in predicting the details of the primary and secondary flow fields both within the bend and downstream thereof.     

Scattered light streaming birefringence in colloidal solutions

Summary The light scattered from particles or molecules in a fluid flow is often used to perform local measurements in a three-dimensional flow field. The scattered light technique has been adopted by McAfee and Pih to the field of streaming birefringence. Their empirical relationship between the observed optical pattern and the state of the flow is extended here to describe more general flow situations. The new set of equations which includes all possible shear and dilatation terms is tested with a number of experiments in the flow of a colloidal solution. The range of validity of these equations is restricted to low-Reynolds number flows.With 14 figures