首页 理论教育推移质颗粒单步运动的统计规律分析

推移质颗粒单步运动的统计规律分析

【摘要】:分析具有代表性的6组试验的推移质颗粒的单步运动,对大量粒子的运动规律进行统计分析。图1-34典型跃移轨迹xz平面图1-33、图1-34中L 为粒子跃移的单步步长,L1为颗粒起跃至达到最高点的长度,H 为跃高,d 是颗粒单步运动横向摆动幅度,uo、vo、wo分别为起跃点的纵向、垂向和横向速度,ue、ve、we则为降落点的纵向、垂向和横向速度。

分析具有代表性的6组试验的推移质颗粒的单步运动,对大量粒子的运动规律进行统计分析。从试验结果来看,颗粒单步步长以及跳跃高度均随水流速度增大而增大,随颗粒粒径和比重增大而减小,较粗颗粒的单步运动距离并没有Einstein假定的100D 这么大。从统计规律来看,推移质单步步长的概率密度服从Γ 分布,颗粒起跃段步长与总步长之比接近一个常数。

1.2.5.1 基础资料

从表1-6中选取的6组资料的颗粒跃移运动的各统计参数见表1-10。

表1-10 颗粒跃移统计参数

在计算分析中,当颗粒在一步运动中的最高点离开床面超过1D 高度后,即认为其进行跃移运动,选取典型轨迹见图1-33和图1-34。

图1-33 典型跃移轨迹xy 平面

图1-34 典型跃移轨迹xz平面

图1-33、图1-34中L 为粒子跃移的单步步长,L1为颗粒起跃至达到最高点的长度(以后简称起跃段长),H 为跃高,d 是颗粒单步运动横向摆动幅度,uo、vo、wo分别为起跃点的纵向、垂向和横向速度,ue、ve、we则为降落点的纵向、垂向和横向速度。定义起跃角η=arctg(vo/uo)、降落角φ=arctg(ve/ue),可以描述粒子跃移和降落时在垂向上的运动幅度,定义颗粒横向摆动幅角δ=arctg(d/L),用来描述颗粒在一次单步运动过程中的横向运动幅度。

表中各参数定义同前,带上划线的量为该量在一组试验中的平均值。实测资料表明:颗粒单步步长最小值在10D 以下,最大则超过100D,步长值与颗粒的本身性质及水流条件有关。根据表1-10的资料可知,步长随颗粒粒径和比重的增大而减小,随水流速度(比降)的增大而增大,颗粒的跃移高度与单步步长有相同的变化趋势。起跃角的变化范围是0.9°~24°,降落角则为0.4°~9.5°,从平均值来看,水流速度越快,起跃角和降落角越小;颗粒粒径和比重越大,颗粒跳跃的角度越小。从横向上的运动来看,颗粒在一次跳跃的过程中左右摆动的幅度最大不超过10D,平均值则基本接近于零,横向摆动幅角的计算结果也体现了同样的规律,说明颗粒跃移时横向上没有运动趋势。

图1-35 与Θ 的关系

颗粒起跃段长L1的平均值与平均步长的比值基本不变,平均为0.368左右,也即L/L1=2.72,见图1-35,这与胡春宏的研究结果是一致的(L/L1=2.79)。

1.2.5.2 单步运动参数的概率密度分布

关于跃高和跃长究竟符合什么分布,从最初Einstein认为的单步距离符合定参数的指数分布,发展到目前普遍承认的Γ 分布,前人已经进行了比较多的讨论和研究。试验结果表明,颗粒单步步长及跃高的概率密度均服从Γ 分布,以粒径0.13cm 组的试验数据为例,图1-36和图1-37分别画出了跃长和跃高的概率密度分布拟合曲线。

图1-36 颗粒跃长的概率分布

图1-37 颗粒跃高的概率分布

图1-36和图1-37中点据即为实测各区的概率密度,曲线则为拟合的Γ分布曲线:

其中A 为调整系数,α和β 是分布参数,当α=1时为指数分布,x0为分布曲线的起点参数。

其他几组试验的统计规律基本相同,具体参数值见表1-11和表1-12。第六组试验由于水流作用力较大,所获得的典型跃移颗粒数较少,故不作分析。

表1-11 步长Γ分布的统计参数

表1-12 跃高Γ分布的统计参数

为了使不同水流和泥沙条件下的实测资料能相互比较,采用相对步长L/L来对五组试验颗粒的单步运动距离进行统计,如图1-38所示。从图可见,相对步长的概率密度分布依然服从Γ 分布,并且不同粒径系列的三组试验数据基本重合,说明颗粒的相对步长的概率密度分布不受粒径的影响;但随比重的减小或水流速度的增大,概率密度分布的峰值变小,也即在轻质沙或高流速的条件下,颗粒的步长分布范围增大,说明颗粒的比重和水流速度对相对步长的概率密度分布有一定的影响。

1.2.5.3 跃移速度的概率密度分布

颗粒单步跳跃过程中在起跃点和降落点处三个方向上的平均速度分量见表1-13。

图1-38 相对步长的概率密度分布

表1-13 起跃速度和降落速度分量

其中uo、 vo、 wo分 别代表颗粒起跃时在纵向、垂向和横向上速度分量的平均值,ue、ve和 we则是颗粒降落时的相应速度平均值。不难发现,颗粒在跃移运动过程中,在距离床面相同高度处,下落过程中的粒子纵向速度总大于上升阶段的速度;横向平均速度的绝对值维持在0.1cm/s以下,可以认为粒子在横向上没有趋势性的运动。

图1-39~图1-44别给出了五组试验中各速度分量的概率密度分布,由图1-39和图1-40可知,无论是相对起跃速度还是相对降落速度,在纵向上的概率密度均接近正态分布:

式中:M 为修正系数;μ 和σ 为分布参数,具体统计数值见表1-14。

图1-39 纵向相对起跃速度的概率分布

图1-40 纵向相对降落速度的概率分布

图1-41 垂向相对起跃速度的概率分布

图1-42 垂向相对降落速度的概率分布

图1-43 横向起跃速度的概率分布

图1-44 横向降落速度的概率分布

表1-14 纵向相对速度概率密度分布参数

表1-15 垂向相对速度概率密度分布参数

图1-41和图1-42表明,在垂向上,由于颗粒与床面之间的碰撞,起跃速度和降落速度出现低速运动的几率大大减小,统计的概率密度分布偏离了正态分布,而更倾向于Γ分布。根据式(1-51),分布参数见表1-15。

表1-16 横向速度概率密度分布参数

图1-43和图1-44进一步表明,横向上的颗粒运动并不受边壁的干扰,速度平均值很小,可以认为粒子只是随机摆动,各统计参数见表1-16。