单向流边界层泥沙起动规律

Shields曲线常用于表示泥沙起动的临界条件,基于边界层理论,对Shields曲线各个流区的线型进行了推证;考虑粘结力的作用,对Shields参数及Shields曲线进行了修正,并给出修正Shields曲线表达式;在此基础之上,从边界层角度重新阐述了Shields曲线。结果表明:Shields曲线在光滑紊流及层流区呈直线分布,在过渡区与阻力系数线型保持一致,在粗糙紊流区呈水平直线分布;修正后Shields曲线与原始Shields曲线在形式上保持一致,修正Shields曲线表达式与实测数据吻合较好,适用于粗、细颗粒泥沙起动条件的计算;Shields曲线事实上代表了Shields参数与沙粒周围绕流流态的关系,同一颗粒处于不同流区起动时,其起动切应力不同。     

波浪及波流边界层泥沙起动规律

基于波浪边界层理论及单向流泥沙起动Shields曲线,推证出波浪泥沙起动Shields曲线;基于波流边界层理论,提出表述波流边界层动力特征的波流比因子X及非线性作用因子Y,并建立了Y与X的相关关系;在此基础之上,结合单向流及波浪泥沙起动Shields曲线,推证出波流共同作用下泥沙起动Shields曲线。结果表明:波浪泥沙起动Shields曲线在层流区与单向流光滑紊流区曲线保持一致,粗糙紊流区与单向流粗糙紊流区曲线保持一致,过渡区线型为折线,由层流区及粗糙紊流区曲线延长交汇获得;X及Y能够合理地表征波流边界层动力对比特征及非线性作用特征;波流泥沙起动Shields曲线介于波浪及单向流泥沙起动Shields曲线之间,随着波流比因子X的不同,依据非线性作用因子Y,自动在波浪及单向流泥沙起动Shields曲线之间非线性过渡。建立的波流泥沙起动Shields曲线与试验结果吻合较好,且能够概括单向流、波浪及波流等不同动力及细沙、粗沙等不同粒径的泥沙起动条件。     

波浪及波流边界层泥沙起动规律

基于波浪边界层理论及单向流泥沙起动Shields曲线,推证出波浪泥沙起动Shields曲线;基于波流边界层理论,提出表述波流边界层动力特征的波流比因子X及非线性作用因子Y,并建立了Y与X的相关关系;在此基础之上,结合单向流及波浪泥沙起动Shields曲线,推证出波流共同作用下泥沙起动Shields曲线。结果表明:波浪泥沙起动Shields曲线在层流区与单向流光滑紊流区曲线保持一致,粗糙紊流区与单向流粗糙紊流区曲线保持一致,过渡区线型为折线,由层流区及粗糙紊流区曲线延长交汇获得;X及Y能够合理地表征波流边界层动力对比特征及非线性作用特征;波流泥沙起动Shields曲线介于波浪及单向流泥沙起动Shields曲线之间,随着波流比因子X的不同,依据非线性作用因子Y,自动在波浪及单向流泥沙起动Shields曲线之间非线性过渡。建立的波流泥沙起动Shields曲线与试验结果吻合较好,且能够概括单向流、波浪及波流等不同动力及细沙、粗沙等不同粒径的泥沙起动条件。     

基于不规则钙质砂颗粒发展的拖曳力系数模型及其在细观流固耦合数值模拟中应用

钙质砂广泛分布于我国南海海域,作为岛礁填筑材料,其渗透性很大程度上决定着填筑后土体的固结和沉降。拖曳力系数是表达液体对土体颗粒表面力的参数,也是表征颗粒状土体渗透能力的一个重要指标,目前国内外对具有极不规则形状钙质砂拖曳力系数的研究十分有限。在前期大量钙质砂颗粒沉降试验基础上,建立了能够考虑形状系数的颗粒与液体相互作用拖曳力系数理论公式;在已建立的计算流体动力学-离散元法(CFD-DEM)流固耦合理论框架下,将新发展的拖曳力系数公式嵌入到流固耦合程序中,模拟钙质砂颗粒在液体中沉降过程。通过将模拟结果与室内沉降试验对比,验证了新建立的拖曳力模型及程序实现的准确性。研究结果表明,较CFD-DEM程序中已有的未考虑颗粒形状的拖曳力系数半经验模型,新的拖曳力系数模型能够更准确地模拟不规则颗粒在液体中沉降过程;同时,在数值建模中无需采用异形颗粒即可以反映颗粒形状对拖曳力的影响,这将大大降低数值模拟的计算量,提高计算效率。研究可进一步应用于钙质砂水中堆填后固结沉降以及冲刷等实际工程问题分析中。     

考虑钙质砂细观颗粒形状影响的液体拖曳力系数试验

钙质砂作为南海岛礁填筑常用的岩土材料,其渗透性很大程度上决定着填筑后土体的固结和沉降。拖曳力系数是表达流体对土体颗粒表面力的参数,也是表征颗粒状土体渗透能力的一个重要参数,目前国内外对钙质砂拖曳力系数的研究十分有限。首先引入一个修正的三维参数 对钙质砂这种天然非规则颗粒材料的形状进行定量描述,然后开展一系列单个钙质砂颗粒在液体中沉降试验,利用高速相机记录颗粒沉降过程,结合图像处理技术获得颗粒沉降平衡速度Ut,进而计算出拖曳力系数CD和雷诺数Re,最后拟合出包含CD、Re及 三个参数的钙质砂拖曳力系数半经验模型。结果发现,在相同雷诺数条件下钙质砂的形状系数 越大,拖曳力系数越小。通过与其他研究结果对比发现,其表面微孔隙越发育,拖曳力系数越小的规律。该模型能够考虑不规则颗粒形状对拖曳力系数的影响,从而提高对土体渗透性预测的精度,对南海岛礁填筑工程中钙质砂固结和沉降的计算也具有重要意义。     

基于不规则钙质砂颗粒发展的拖曳力系数模型及CFD-DEM流固耦合数值方法

钙质砂广泛分布于我国南海海域,作为岛礁填筑材料,其渗透性很大程度上决定着填筑后土体的固结和沉降。拖曳力系数是表达液体对土体颗粒表面力的参数,也是表征颗粒状土体渗透能力的一个重要指标,目前国内外对具有极不规则形状钙质砂拖曳力系数的研究十分有限。在前期大量钙质砂颗粒沉降试验基础上,建立了能够考虑形状系数的颗粒与液体相互作用拖曳力系数理论公式;在已建立的CFD-DEM流固耦合理论框架下,将新发展的拖曳力系数公式嵌入到流固耦合程序中,模拟钙质砂颗粒在液体中沉降过程。通过将模拟结果与室内沉降试验对比,验证了新建立的拖曳力模型及程序实现的准确性。研究结果表明,较CFD-DEM程序中已有的未考虑颗粒形状的拖曳力系数半经验模型,新的拖曳力系数模型能够更准确地模拟不规则颗粒在液体中沉降过程;同时,在数值建模中无需采用异形颗粒即可以反映颗粒形状对拖曳力的影响,这将大大降低数值模拟的计算量,提高计算效率。研究可进一步应用于钙质砂水中堆填后固结沉降以及冲刷等实际工程问题分析中。     

基于不规则钙质砂颗粒发展的拖曳力系数模型及CFD-DEM流固耦合数值方法

钙质砂广泛分布于我国南海海域,作为岛礁填筑材料,其渗透性很大程度上决定着填筑后土体的固结和沉降。拖曳力系数是表达液体对土体颗粒表面力的参数,也是表征颗粒状土体渗透能力的一个重要指标,目前国内外对具有极不规则形状钙质砂拖曳力系数的研究十分有限。在前期大量钙质砂颗粒沉降试验基础上,建立了能够考虑形状系数的颗粒与液体相互作用拖曳力系数理论公式;在已建立的CFD-DEM流固耦合理论框架下,将新发展的拖曳力系数公式嵌入到流固耦合程序中,模拟钙质砂颗粒在液体中沉降过程。通过将模拟结果与室内沉降试验对比,验证了新建立的拖曳力模型及程序实现的准确性。研究结果表明,较CFD-DEM程序中已有的未考虑颗粒形状的拖曳力系数半经验模型,新的拖曳力系数模型能够更准确地模拟不规则颗粒在液体中沉降过程;同时,在数值建模中无需采用异形颗粒即可以反映颗粒形状对拖曳力的影响,这将大大降低数值模拟的计算量,提高计算效率。研究可进一步应用于钙质砂水中堆填后固结沉降以及冲刷等实际工程问题分析中。     

基于不规则钙质砂颗粒发展的拖曳力系数模型及CFD-DEM流固耦合数值方法

钙质砂广泛分布于我国南海海域,作为岛礁填筑材料,其渗透性很大程度上决定着填筑后土体的固结和沉降。拖曳力系数是表达液体对土体颗粒表面力的参数,也是表征颗粒状土体渗透能力的一个重要指标,目前国内外对具有极不规则形状钙质砂拖曳力系数的研究十分有限。在前期大量钙质砂颗粒沉降试验基础上,建立了能够考虑形状系数的颗粒与液体相互作用拖曳力系数理论公式;在已建立的CFD-DEM流固耦合理论框架下,将新发展的拖曳力系数公式嵌入到流固耦合程序中,模拟钙质砂颗粒在液体中沉降过程。通过将模拟结果与室内沉降试验对比,验证了新建立的拖曳力模型及程序实现的准确性。研究结果表明,较CFD-DEM程序中已有的未考虑颗粒形状的拖曳力系数半经验模型,新的拖曳力系数模型能够更准确地模拟不规则颗粒在液体中沉降过程;同时,在数值建模中无需采用异形颗粒即可以反映颗粒形状对拖曳力的影响,这将大大降低数值模拟的计算量,提高计算效率。研究可进一步应用于钙质砂水中堆填后固结沉降以及冲刷等实际工程问题分析中。     

泥沙颗粒形状对非均匀沙起动条件的影响

根据非均匀沙颗粒的不均匀性,引入颗粒形状系数,通过分析泥沙颗粒形状系数和相对暴露度的影响,建立了等效粒径表达式;并通过泥沙起动受力分析,推导得到泥沙起动临界条件表达式.实测资料验证表明,计算值与实测值符合较好.     

泥沙颗粒生长生物膜后起动的实验研究——Ⅰ.起动实验设计及结果分析

天然水体中泥沙颗粒表面容易生长生物膜,生物膜的生长会对泥沙颗粒的运动特性产生显著影响。设计了泥沙颗粒表面生长生物膜的实验,以及生长生物膜后泥沙颗粒的起动水槽实验。在水槽观测段分别放置由去离子水浸泡和天然水培养的两种泥沙,对比观察两种泥沙起动过程的区别,分析起动流速的差异和随培养时间的变化规律。结果表明,天然富营养水体中生物膜生长繁茂,对泥沙起动特性有着显著的影响,能够增强泥沙颗粒的抗冲性,并随时间呈现出先逐步增强到一个最优值后有所削弱的钟形曲线的变化趋势,表明在研究天然河流和湖泊尤其是污染严重的水域的泥沙起动时,需要考虑到生物膜的作用。     

Geometrical mechanism of inverse grading in grain-flow deposits: An experimental revelation

The grain-flow has so far been defined with reference to the distinctive sediment-support mechanism, the dispersive pressure. The role of sediment-support mechanism, however, is required in a multiphase flow to prevent the gravitational settling of the particles through the driving medium during the flow. In a single-phase flow of non-cohesive grains no such secondary mechanism is required to counteract the gravitational pull, the driving force of the flow. So the definition of grain-flow needs a critical revision. This, in turn, involves proper understanding of the grain-flow mechanism, so that the relation between the process and the product can be properly established. The most distinctive feature often demonstrated by a grain-flow deposit is the particle size segregation, which leads to the development of inverse grading. The available explanations for this phenomenon find theoretical constraints. In the present study an attempt was made to understand the mechanism of single-phase non-cohesive granular flow of different flow regime and the particle segregation pattern in the resultant deposit through laboratory simulation. The experimental observations revealed that no sustained granular flow sets in on a slope deviating much from the limiting value of the angle of repose of the granular material. A persistent simple shear flow develops on slopes of this critical value. Each of the grains rolls in response to simple shearing. If the shear stress attains a critical value, theoretically the larger grains can even climb up the adjacent smaller ones towards the down-slope direction. In reality, however, high angle climb is not very common. The larger grains preferably roll over the smaller grains when the common tangent becomes almost horizontal or makes a very low angle with the direction of flow, and by this process gradually reaches the upper surface of the flow causing the development of inverse grading. The upper surface of the resultant deposit remains parallel to the sloping substratum. These properties readily distinguish this variety of granular flow from the other natural flows, and the flow may thus be assigned the distinct status of grain-flow.     

基于物理模型试验的碎屑流流态化运动特征分析

流态化运动是高速远程滑坡的主要运动形式,是揭示高速远程滑坡运动机理的重要基础。基于粒子图像测速（PIV）分析方法,采用物理模型试验对不同粒径组成条件下的颗粒流内部的速度分布、剪切变形及流态特征进行了研究,并对高速远程滑坡流态化运动特征进行了讨论分析。结果表明:碎屑流流态化运动特征与颗粒粒径呈显著的相关性,随着粒径的减小或细颗粒含量的增加,颗粒流底部相对于边界的滑动速度以及整体的运动速度均呈逐渐减小的趋势,颗粒流内部剪切变形程度增加,颗粒的运动形式由“滑动”向“流动”转变;当颗粒粒径较小或细颗粒含量较高时,颗粒流内部剪切速率增大的趋势在颗粒流底部更加显著,反映了粒径减小有助于促进颗粒流内部剪切向底部的集中;在同一颗粒流的不同运动阶段及不同纵向深度,其流态特征具有显著差别,颗粒流前缘及尾部主要呈惯性态,颗粒间以碰撞作用为主,而主体部分则主要呈密集态,颗粒间以摩擦接触作用为主;在颗粒流表面及底部,颗粒间相互作用方式主要是碰撞作用,中间部分则以摩擦作用为主;对于不同粒径的颗粒流,随着粒径的增大或粗颗粒含量的增加,颗粒流内部颗粒的碰撞作用加强,颗粒流整体趋于向惯性态转变。     

考虑颗粒棱角影响的直剪试验的离散元模拟

为了研究颗粒棱角对颗粒材料力学行为的影响,建立了具有不同棱角度的对称多面体颗粒,采用了一种简单并适合任意颗粒形状的接触本构模型,对三维离散元开源程序YADE进行了修改,研究了颗粒棱角度在模拟直剪试验中的影响以及接触力各向异性在剪切过程中的演化规律。研究结果表明,颗粒棱角度越小,颗粒间相互咬合自锁的作用越小,颗粒受剪更易转动,致使颗粒体系的剪切强度和剪胀性下降;竖向加载力越大,颗粒棱角度的影响越明显;法向接触力的各向异性在剪切过程中表现为先增后减最后趋向稳定的趋势;法向接触力的各向异性变化程度随颗粒棱角度的增大而增大。     

Critical length of force chains and shear band thickness in dense granular materials

This paper investigates the existence of the critical force chain length and the buckling of unconfined grain columns in dense granular materials. Tests on assemblies of flat pentagon photoelastic particles were first carried out to demonstrate the maximum length of force chains. Then, the theoretical buckling analysis and distinct element method (DEM) simulations for grain columns composed of mono-sized elliptical particles were performed. The results revealed the existence of critical column length, which is generally affected by the particle shapes, the rotational resistance at particle contact points and the end constraints to the grain columns. The interparticle friction does not have explicit effect on the critical force chain length, but it has significant influence on the grain column’s curvature when collapse takes place. The thickness of shear band in granular soils can be determined as the critical length of grain columns by appropriately imposing the constraints on the boundaries, as confirmed by DEM simulations and experimental results.     

基于水下摄影的床面泥沙运动特性试验研究

为提高推移质试验的观测精度与效率,基于水下摄影和粒子跟踪（UP/PTV）技术,实现明槽流床面泥沙运动状态的精细试验观测。通过剔除床面颗粒震颤干扰、设置颗粒临界运动阈值和多重滤波筛选程序,提高采集样本数据的有效性。利用多组低强度推移质试验,提取床面颗粒运动轨迹、速度、单步时长等数据,并进行粒子运动的Lagrange过程分析和概率密度分布（PDF）研究。研究表明:粒间碰撞和近底紊流扫荡的影响使粒子速度在单步步长内呈现先急剧增加再缓慢衰减的变化特征;速度PDF曲线显示细尾Gamma函数特性,同时受粒子震颤效应影响,加速度PDF曲线则具有拉普拉斯分布特征;粒子速度与摩阻流速之间关系密切,保持3.4~3.5的比值;粒子单步时长与步长的联合分布呈幂函数变化规律,拟合曲线指数一般为1.25~1.3。幂律指数大小受推移质输沙强度与床面粒子异质性的影响。     

Role of particle crushing on particle kinematics and shear banding in granular materials

The paper provides an in-depth exploration of the role of particle crushing on particle kinematics and shear banding in sheared granular materials. As a two-dimensional approximation, a crushable granular material may be represented by an assembly of irregularly shaped polygons to include shape diversity of realistic granular materials. Particle assemblies are subjected to biaxial shearing under flexible boundary conditions. With increasing percentage of crushed particles, mesoscale deformation becomes increasingly unstable. Fragmented deformation patterns within the granular assemblies are unable to form stable and distinct shear bands. This is confirmed by the sparsity of large fluctuating velocities in highly crushable assemblies. Without generating distinct shear bands, deformation patterns and failure modes of a highly crushable assembly are similar to those of loose particle assemblies, which are regarded as diffuse deformation. High degrees of spatial association amongst the kinematical quantities confirm the key role that non-affine deformation and particle rotation play in the generation of shear bands. Therefore, particle kinematical quantities can be used to predict the onset and subsequent development of shear zones, which are generally marked by increased particle kinematic activity, such as intense particle rotation and high granular temperature. Our results indicate that shear band thickness increases, and its speed of development slows down, with increasing percentage of crushed particles. As particles crush, spatial force correlation becomes weaker, indicating a more diffuse nature of force transmission across particle contacts.     

Threshold of motion of bivalve and gastropod shells under oscillatory flow in flume experiments

The threshold of motion of non-fragmented mollusc shells was studied for the first time under oscillatory flow. In this regard, flume experiments were used to investigate the threshold of motion of three bivalve and three gastropod species, two typical mollusc classes of coastal coquina deposits. The sieve diameters ranged from 2·0 to 15·9 mm. These experiments were performed on a flat-bottom setup under regular non-breaking waves (swell) produced by a flap-type wave generator. The critical Shields values for each species of mollusc were plotted against the sieve and nominal diameter. Moreover, the dimensionless Corey shape factor of the shells was evaluated in order to investigate the effect of mollusc shell shapes on the threshold of motion. According to their critical Shields parameter, the mollusc threshold data under oscillatory flow present smaller values than the siliciclastic sediments when considering their sieve diameter. In addition, the mollusc datasets are below the empirical curves built from siliciclastic grain data under current and waves. When considering the nominal diameter, the critical Shields parameter increases and the mollusc data are closer to siliciclastic sediments. Bivalves, which have a flat-concave shape (form factor: 0·27 to 0·37), have a higher critical Shields parameter for smaller particles and more uniform datasets than the gastropod scattered data, which have a rounded shape (form factor: 0·58 to 0·62) and have varied morphologies (ellipsoidal, conical and cubic). The comparison between previous current-driven threshold data of bioclastic sediment motion and the data of mollusc whole shells under oscillatory flow shows a fair correlation on the Shields diagram, in which all datasets are below the mean empirical curves for siliciclastic sediments. These findings indicate that the shape effect on the transport initiation is predominant for smaller shells. The use of the nominal diameter is satisfactory to improve the bioclastic and siliciclastic data correlation.     

Experimental study of the dynamic behavior and segregation of density-bidisperse granular sliding masses at the laboratory scale

From the understanding of dynamics and processes of rapid granular flows and the granular-segregation mechanism in gravity-driven flow, we can clarify the particle-composition structure in the downstream areas of avalanches in geophysical contexts, such as landslides, rock falls, and snow-slab avalanches. Such dynamics also provide a basis for geophysical studies. This study experimentally investigates the dynamic behavior and segregation phenomena of a density-bidisperse, rapid, granular flow down a quasi-2D, rough, inclined rectangular chute. Particles with two density ratios are used to investigate the mechanism of density-induced segregation, and four chute-inclination angles are tested to examine the influence of driving forces. The dynamics of the mixture flow—which includes the flow-depth evolution, stream-wise and depth-wise velocity profiles, shear rate, and granular temperature in the upper high-shear band of the flow—are obtained from particle image velocimetry (PIV) measurements. The two-dimensional concentration distributions of the particles in the stream-wise direction are also obtained using 2D image processing to determine the segregation state. In the upstream region, the variation in the concentration of heavier particles is defined as the strength of the density-induced segregation state, S_d. Our results indicate that the mixture-flow parameter—particularly the shear rate and the granular temperature in the upper high-shear band—crucially influence the strength of particle segregation in granular avalanches. In the upstream region, a higher shear rate and a higher granular temperature in the upper high-velocity band result in a smaller drag force in the mixture flow, causing stronger density-induced particle segregation. These results well describe the entire processes of dense granular flows, from upstream initiation to the downstream steady state. Therefore, they reveal the structure of the mixed flow in the depth direction and are expected to explain various gravity-driven mixture granular flows.