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

1、波浪边界层中细颗粒粘性泥沙的再悬浮和扩散输移摘要：本文用多尺度摄动方法从理论上推导了波浪边界层中细颗粒粘性泥沙的再悬浮和扩散输移的规律,并用几个算例细致地分析了波浪对细颗粒粘性泥沙再悬浮和扩散输移的作用。指出:在近岸及湖泊地区,波浪不仅是细颗粒粘性泥沙再悬浮的主要动力,而且其引起的泥沙扩散输移的强度和风生流相当。关键词：波浪边界层粘性泥沙再悬浮扩散输移1引言底部泥沙的再悬浮和运移是海岸过程的中心问题。由于泥沙往往是重金属和氮、磷等污染物的载体,近年来,泥沙的扩散输移得到了广泛的研究并建立了许多数学模型。Hrikaa(1968)和VanRijn(1993)等讨论了非粘性泥沙扩散输移的数学模型1,

2、20,Sheng(1979,1984,1991a,b)等那么建立了粘性泥沙的数学模型1216。这些模型在计算技术和紊动形式方面比较成熟,但许多特征参数尚赖于野外和实验室资料的率定,这就增加了模型的不确定性,同时这些模型仅仅考虑到波浪对底部切应力的奉献,即对泥沙再悬浮的奉献。事实上,除了波浪对泥沙的侵蚀作用外,对悬浮泥沙的扩散输移也有重要作用。波浪二阶雷诺应力产生的波生流流速随程度方向变化,和风生流流速相当。如美国keehbee湖,水深较深区域的实测流速为510/s,而一波高15,周期3s的代表波,其一阶轨道速度为30/s,产生的流动大于或与风生流流速相当。同深海波生流会改变表层风生Ekan漂移

3、(adsen,1978)一样,此流动也会使近底水质点产生漂移。同时波浪边界层中较大的切应力也将加剧泥沙的扩散输移(Eihian1994)。因此,一个完好的泥沙模型需同时包含波和流。为了更好地理解波浪边界层中泥沙的再悬浮和扩散输移,本文重点从理论上分析细颗粒粘性平坦海床在简谐波边界层中泥沙的再悬浮和扩散输移。2根本方程对波浪边界层的紊动理论已进展了大量的研究并提出了许多半经历半理论模型和数学模型(紊动粘性,混合长度,K、二阶模型等),Sleath(1990)对此进展了回忆和总结指出17,不管紊动粘性系数随时间或空间变化,其给出的速度剖面与常紊动粘性系数的结果相差不大。一些实验说明(Jnssnar

4、lsn1976;Hrikaa,1968)紊动扩散系数随时间变化,而这种变化会对流动产生影响(Trbridgeadsen1984)，。但由于实验数据太少,缺乏以得到可信的形式,为了得到解析解,本文采用常紊动粘性系数和常紊动扩散系数。泥沙的扩散方程为(1)式中为体积含沙量,0为泥沙沉速,Dh,D分别为程度和垂直方向上的紊动扩散系数,i=1,2,(X1,X2)=(X,Y),(U,2)=(u,)分别为程度坐标及流速分量,Z和为垂直坐标及流速分量。对粘性细颗粒泥沙,底部边界条件为(VanRijn,1994)Z=0(2)式中D和通常有如下形式(Patheniades,1965;Krne,1962)(3)(

5、4)其中b、d、分别为底部切应力,泥沙临界淤积切应力和临界冲刷切应力,d,和那么为系数。在大部分河口及湖泊地区,底床表层沉积物呈半固结或非固结状态,很难承受切应力,其d,非常小,伽马仪都很难测到(Patheniades,1965,ehtaPatheniades,1982)6,9。本文只研究泥沙的再悬浮和扩散输移,波浪产生的底部切应力通常大于,因此不失一般性,本文忽略(2)中的沉积项并简化为Z=0(5)在边界层外=0Z(6)假设给定含沙量的初始条件,此问题变成为可冲刷边界条件下泥沙的长时间扩散问题。3量级分析s=D/0(7)另两个特征长度为两种边界层厚度,分别对应于动量和物质扩散。(8)其中波浪

6、圆频率。一般地,假设三个特征长度相当,也就是(s)=(u)=()(9)那么Shidt数的量级为S=e/Dv=u/=(1)(10)如今我们考虑小振幅简谐波,其圆频率足够大使得=KA1=K1(11)式中K为波数，,A为近底水质点波动半径。不失一般性,假设=(),我们可以引入如下无量纲量Xi=KXi,Z*=Z,t*=t(12)*=0,UI*=UiA,*=KA那么方程(6)保持形式不变,方程(1)成为(13)式中Pe=0D是Pelet数。而底部切应力的量级0为(14)那么可得含沙量量级0为(15)底部边界条件(2)可无量纲化成(16)式中(17)转贴于论文联盟.ll.4扩散方程的多尺度摄动展开通过以上

7、量级分析,我们把量级参数放入无量纲方程(13)和(16)可得(18)Z=0(19)对应垂直和程度长度尺度,此对流扩散过程有两个时间尺度。扩散穿越边界层的时间尺度和波周期相当为(-1)=(2D),而扩散波长的间隔 的时间尺度为(1K2Dh)这两个时间尺度的比为(KDhD)(。考虑到假设=(),本文引进两种时间坐标t和T=t,那么速度和含沙量可摄动展开为Ui=U()i+U()i+()(20)=()+()+(3)(21)=(0)+(1)+()()(22)式中U(n)和(n)为Xi、Z和t的函数,而(n)=(n)(Xi,Z,t,T)。参照Ei和hian(1994)的工作,我们可推得边界层上(垂线平均含

8、沙量)的控制方程(23)式中f为f(z)沿边界层的垂向积分。这是EI和hian(1994)工作的推广,其差异在于ei和hian未考虑泥沙的再悬浮,即没有方程(23)右边的代表底床侵蚀的源项。同样的方程也曾用在定常河流的冲蚀中(VanRijn,1994)8。0上式中产生的对流并产生明显的扩散,这类似于泰勒扩散,第二项那么由紊动产生(eihain,1994)。设边界层外理想流的首阶流速为Ui=R(V0ie-it)(24)式中U0i=U0i(Xj),那么方程(23)可化为(25)(26)(27)(28)式中U和V为欧拉流的垂向加权平均,向量Eij由U()i()得到,系数Hn(n=1,2,3,4)为P

9、e和Se的系数,eihian(1994)曾给出其详细形式和图8。由于H4为复数,那么Eij非对称,也就是EijEji。式(25)右边的源项为一个波周期内的平均侵蚀速率,可以由底部切应力求出(29)下面我们用几个算例研究不同波型和可侵蚀底床几何条件的影响。在所有算例中Dh=D=Ue=D,Pe=1,S=1,s=u=0。由eihian(1994)可知8H1=-0.122058+0.659452I(30)H2=0.033333H3=-0.155391+0.659452IH4=0.023615+0.233866I5前进波作用下泥沙的再悬浮和输移设底床仅0XL,-y可冲刷,一波列垂直此可冲刷无穷条带行进,

10、近底理想流为U(1)=R0(x)e-it=R(Ubeikxe-it)(31)Ub=A=A0/sinhkh(32)式中A0为外表波波高。底部切应力为b=(D/)RU0(X)e-i(t+4)(33)时间平均的侵蚀速率为b=22(34)引入无量纲量X=kx,T=(kb)/T,U0=U0/Ub,=/(35)其中0=2PD/(kb)(36)无量纲系数为U=U/kUb=I(H),D=D/U2b,Exx=ExxU2b=R(H)(37)可得(38)其中(39)此问题的解析解为(40)图1前进波作用下泥沙再悬浮和输移Sedientresuspensinandtransprtatinunderatinsfprgr

11、essiveaves图1为=0.659452,E=0.023615,D=0.001,L=2时的解析解。初始,底部含沙浓度仅在可冲刷条带上方增加,其后部向前传递,而前部那么向波前进方向延展,较长时间以后,在可冲刷区域上形成稳定的含沙浓度。0XL(41)XL(42)在此算例中KU=0.0566?1,因此稳定的含沙量主要由上式第一项决定。当含沙水体前部进入非可冲刷区域,稳定的含沙量也随之进入此区域。当然这只能在可冲刷区域无限可供沙和波浪不衰减这种理想情况下才有可能。假设可冲刷区域为0X2,-Y,那么控制方程为(43)其中(44)此方程可用解析或有限差分法求解。图2为含沙量随时间的变化。含沙水体先从可

12、冲刷区域产生,逐渐成长为一圆角的矩形平台,平台前部随时间沿波传播方向延展,后部同样向前传递。由于紊动粘性较小,横向运移非常缓慢。较长时间后,在平台的后部到达一有限的高度。图2前进波作用下泥沙再悬浮和输移Sedientresuspensinandtransprtatinundertheatinsfprgressiveaves6垂直海墙前入射波浪作用下的泥沙再悬浮和扩散输移设y轴和垂直海墙一致。一列单射波向海墙入射并产生反射,那么外表水质点位移为=2A0s(kxs)e(ikysin-t)(45)近底理想流流速为U0=2iUbssin(kxs)eikysin(46)V0=2Ubsins(kxs)ei

13、kysin(47)同样地用式(37)无量纲化,那么扩散方程中的无量纲系数为(eiandhian,1994)U=2R(H)sR(H)-R(H)sinsin(2)sin(2xs)(48)Exx=4R(H)ssin(Xs)(49)而冲刷项变为(50)其中a=ssin(xs),b=sins(xs)(51)(52)由于输运方程中所有系数和强迫项都与Y无关,的控制方程可简化为一维(53)其中(54)(55)由(48)和(49)可知U和Kxx为X的周期函数,图3(a)为=4的U和KXX。从图3(a)可图3垂直海墙入射波作用下的侵蚀速率Ersinratefavesrefletedfraseaall图4垂直海墙

14、入射波作用下的泥沙再悬浮和输移Sedientresuspensinndtransprtatinfavesrefletedfraseaall知,对流项的零点和扩散项的最小值相应,从而泥沙向直线2Xs=n聚集,在此处扩散最弱。图3(b)为正那么侵蚀速率,图4为可冲刷区域为0X4时的含沙量分布随时间的变化。明显地,泥沙趋向于向对流速度转向处聚集。7结语本文从理论上论证了波浪对泥沙再悬浮和扩散输移的作用,并通过几个算例进展了详细阐述。由于波生流的强度与风生流相当,因此今后工作中应不仅仅考虑波浪对底切力和冲刷速度的奉献。转贴于论文联盟.ll.参考文献1Hrikaa.K.A.atanabe.(1968).

15、Labratrystudynsillatrybundarylayerfl.atalEng.Japan,11,13-28.2Jnssn,I.G.arlsn,N.A.(1976)Experientalandtheretialinvestigatininansillatryturbulentbundarylayer.J.Hydrau.Res.14,45-60.3Krne,R.B.(1962).Fluestudiesfthetransprtfsedientinestuarialshalingpresses.Univ.alif.Hydrau.Eng.SanitaryEng.Res,Lab,Berkley

16、.110pp.4adsen,.S,(1978)asstransprtindeep-ateraves,J.Phy.an.8,1009-1015.5ehta,A.(1984),haraterizatinfhesivesedientprpertiesandtransprtpressesinestuaries.inEstuarinehesivesedientdynais(ed。)A.ehta,Springer-verlag.290-325.6ehta,A.J.Partheniades,E,(1982).Resuspensinfdepsitedhesivesedientbeds.astalEnginee

17、ring.1569-1588.7Ei,.(1989)AppliedDynaisfeanSurfaeaves,rldSientifi.8EI,.hina,.(1994)Dispersinfsallsuspendedpartilesinaavebundarylayer,J.Phy.ean.24,2479-2495.9Partheniades,E.(1965)Ersinanddepsitinfhesivesils.J.HydraulisDiv.ASE,91HYI.105-139.10Partheniades,E.(1971).Ersinanddepsitinfhesiveaterials.River

18、ehanis,H.ShenEd.Vl.。h.25.11Sanfrd,L.P.(1994)ave-induedresuspensinfupperhesapeakeBayuds.Estuaries,17,148-165.12Sheng,Y.P.(1984)delingbttbundarylayerandhesivesedientdynaisinestuarineandastalaters.inEstuarinehesivesedientdynaisTBZ(ed).A.ehta,Springer-verlag.360-400.13Sheng.Y.P.Eliasn,D.E.hiJ.K.hen,X.Le

19、e,H,K.(1991).Effetfsedientresuspensinnnuerialsiulatinfthree-diensinalind-drivenirulatinandsedienttransprtinLakekeehbeeduringspring1989.UniversityfFlridaastalandeangraphiEngineeringLabratryReprt91-019.14Sheng.Y.P.hen,X.(1991),Athree-diensinalnuerialdelfhydrdynais,sedienttransprtandphsphrusdyaisinLakekeehbee:Thery,deldevelpentandduentatin.UniversityfFlridaastalandeangraphiEngineeringLabratryReprt91-2393.15Sheng,Y.P.,k,V.Resuspensinandvertia