波浪边界层中细颗粒粘性泥沙的再悬浮和扩散输移.docx

波浪边界层中细颗粒粘性泥沙的再悬浮和扩散输移摘要本文用多尺度摄动方法从理论上推导了波浪边界层中细颗粒粘性泥沙的再悬浮和扩散输移的规律,并用几个算例细致地分析了波浪对细颗粒粘性泥沙再悬浮和扩散输移的作用。指出:在近岸及湖泊地区,波浪不仅是细颗粒粘性泥沙再悬浮的主要动力,而且其引起的泥沙扩散输移的强度和风生流相当。关键词波浪边界层粘性泥沙再悬浮扩散输移1引言底部泥沙的再悬浮和运移是海岸过程的中心问题。由于泥沙往往是重金属和氮、磷等污染物的载体,近年来,泥沙的扩散输移得到了广泛的研究并建立了许多数学模型。Horikawa(1968)和VanRijn(1993)等讨论了非粘性泥沙扩散输移的数学模型[1,20],Sheng(1979,1984,1991a,b)等则建立了粘性泥沙的数学模型[12～16]。这些模型在计算技术和紊动模式方面比较成熟,但许多特征参数尚赖于野外和实验室资料的率定,这就增加了模型的不确定性,同时这些模型仅仅考虑到波浪对底部切应力的贡献,即对泥沙再悬浮的贡献。事实上,除了波浪对泥沙的侵蚀作用外,对悬浮泥沙的扩散输移也有重要作用。波浪二阶雷诺应力产生的波生流流速随水平方向变化,和风生流流速相当。如美国Okeechobee湖,水深较深区域的实测流速为5～10cm/s,而一波高15cm,周期3s的代表波,其一阶轨道速度为30cm/s,产生的流动大于或与风生流流速相当。同深海波生流会改变表层风生Ekman漂移(Madsen,1978)一样[４],此流动也会使近底水质点产生漂移。同时波浪边界层中较大的切应力也将加剧泥沙的扩散输移(Mei&Chian1994)[８]。因此,一个完整的泥沙模型需同时包含波和流。为了更好地理解波浪边界层中泥沙的再悬浮和扩散输移,本文重点从理论上分析细颗粒粘性平坦海床在简谐波边界层中泥沙的再悬浮和扩散输移。2基本方程对波浪边界层的紊动理论已进行了大量的研究并提出了许多半经验半理论模型和数学模型(紊动粘性,混合长度,K～ε、二阶模型等),Sleath(1990)对此进行了回顾和总结指出[17],不论紊动粘性系数随时间或空间变化,其给出的速度剖面与常紊动粘性系数的结果相差不大。一些实验表明(Jonsson&Carlson1976;Horikawa,1968)紊动扩散系数随时间变化,而这种变化会对流动产生影响(Trowbridge&Madsen1984)[１，２]。但由于实验数据太少,不足以得到可信的模式,为了得到解析解,本文采用常紊动粘性系数和常紊动扩散系数。泥沙的扩散方程为式中C为体积含沙量,W0为泥沙沉速,Dh,Dυ分别为水平和垂直方向上的紊动扩散系数,i=1,2,(X1,X2)=(X,Y),(U１,Ｕ2)=(u,υ)分别为水平坐标及流速分量,Z和W为垂直坐标及流速分量。对粘性细颗粒泥沙,底部边界条件为(VanRijn,1994)式中D和ε通常有如下形式(Patheniades,1965;Krone,1962)其中τb、τd、τc分别为底部切应力,泥沙临界淤积切应力和临界冲刷切应力,τd＜τc,m和M则为系数。在大部分河口及湖泊地区,底床表层沉积物呈半固结或非固结状态,很难承受切应力,其τd,τc非常小,伽马仪都很难测到(Patheniades,1965,Mehta&Patheniades,1982)[6,9]。本文只研究泥沙的再悬浮和扩散输移,波浪产生的底部切应力通常大于τc,因此不失一般性,本文忽略(2)中的沉积项并简化为在边界层外如果给定含沙量的初始条件,此问题变成为可冲刷边界条件下泥沙的长时间扩散问题。3量级分析本问题共有三个特征长度。一是稳定含沙量区厚度δs,在此区域泥沙重力和垂直扩散相平衡。另两个特征长度为两种边界层厚度,分别对应于动量和物质扩散。其中ω波浪圆频率。一般地,假设三个特征长度相当,也就是则Schmidt数的量级为现在我们考虑小振幅简谐波,其圆频率足够大使得式中K为波数，,A为近底水质点波动半径。不失一般性,假设∈=O(β),我们可以引入如下无量纲量则方程(6)保持形式不变,方程(1)成为式中Pe=W0δc／Dυ是Peclet数。在边界层中选取长度为dx的微元,微元中周期平均含沙量C的变化由水平方向上波生平流的对流和紊动扩散以及底床的冲刷决定,波生平流流速的量级为O(ωΚΑ２）,则我们可得而底部切应力的量级τ0为则可得含沙量量级C0为底部边界条件(2)可无量纲化成式中(17)4扩散方程的多尺度摄动展开通过以上量级分析,我们把量级参数放入无量纲方程(13)和(16)可得对应垂直和水平长度尺度,此对流扩散过程有两个时间尺度。扩散穿越边界层的时间尺度和波周期相当为O(ω-1)=O(δ2／Dυ),而扩散波长的距离的时间尺度为O(1／K2Dh)这两个时间尺度的比为O(K２δ２≡O(