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

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