CN103940580B - 一种计算波浪作用下悬沙输移相位滞后效应的方法 - Google Patents
一种计算波浪作用下悬沙输移相位滞后效应的方法 Download PDFInfo
- Publication number
- CN103940580B CN103940580B CN201410155863.6A CN201410155863A CN103940580B CN 103940580 B CN103940580 B CN 103940580B CN 201410155863 A CN201410155863 A CN 201410155863A CN 103940580 B CN103940580 B CN 103940580B
- Authority
- CN
- China
- Prior art keywords
- wave
- granule
- half period
- particle
- speed
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Underground Or Underwater Handling Of Building Materials (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明提供了一种计算波浪作用下悬沙输移相位滞后效应的方法,包括将颗粒动能转换为势能确定在波浪半周期内颗粒的悬浮高度;通过对指数分布沿深度进行积分平均,获得颗粒悬浮层内水质点的深度平均的垂向流速;根据颗粒在沉降过程中受到颗粒之间的相互碰撞作用后的受阻碍沉速和颗粒悬浮层内水质点的深度平均的垂向流速修正获得波浪边界层内悬沙颗粒的沉速;根据波浪半周期内悬浮颗粒的沉降时段、波浪边界层内悬沙颗粒的沉速和颗粒的悬浮高度计算获得波浪半周期的相位滞后效应参数。本发明能够反应水质点垂向流速对悬浮颗粒在波峰半周期和波谷半周期内沉降的影响,考虑了颗粒之间的相互碰撞作用造成的能量耗散比例,因此计算结果更为准确。
Description
技术领域
本发明涉及一种计算波浪作用下悬沙输移相位滞后效应的方法,属于海岸工程泥沙输移计算技术领域。
背景技术
近岸地区波浪作用下的泥沙输移是一种非恒定输沙过程。相位滞后效应反映的是水体中悬沙的含量与近底水动力变化不同步的现象,其本质是在波浪前半个周期内被悬浮的泥沙不能在该半周期内及时沉降而是保持悬浮直至下半个周期内向相反方向继续输移沉降的过程。泥沙颗粒粒径越细、波浪周期越短、近底流速越大,波浪的非对称性越强、相位滞后效应就越明显。相位滞后效应的强弱能够影响波浪作用下净输沙率的大小和方向。
现有的技术方案中对相位滞后效应的定量刻画均是采用参数化处理,即悬浮颗粒沉降到床面上所需要的时长tfall与波浪周期内可供悬沙颗粒沉降时段t的比值P。
因此,现有的技术方案就集中体现在如何确定波浪作用下悬浮颗粒沉降到床面上所需要的时长tfall和波浪周期内可供悬沙颗粒沉降的时段t两个方面。现有的技术方案在确定二者时一般采用半周期的方式,即针对波浪波峰通过时的半周期Tc和波谷通过时的半周期Tt分别确定。具体的,在确定tfall时,现有的技术方案处理为悬沙颗粒被悬浮的高度Δs与颗粒沉速ωs的比值,即:
但是如何确定悬沙颗粒被悬浮的高度Δs和反映波浪边界层内水质点的垂向流速对颗粒沉速ωs以及颗粒在沉降过程中受到颗粒之间相互碰撞阻碍的影响等方面,特别是在确定制约颗粒沉降历时的因素时均采用经验性的处理,存在着较大的误差,具体包括:
(1)在确定悬沙颗粒被悬浮的高度Δs时,现有的技术方案或者假定是沙纹波高或者临底高含沙层厚度的数倍而具体的倍数需要依据实验数据进行率定,或者假定悬沙颗粒运动的动能全部转化为势能来计算颗粒的悬浮高度,而没有考虑颗粒的部分动能由于颗粒之间的碰撞作用被耗散掉了。这些方面的不足都造成了计算结果的误差。
(2)在考虑颗粒的沉速受到波浪边界层内水质点的垂向运动速度的影响时,如何计算波浪边界层内水质点的垂向运动速度是一个关键的技术环节。现有的技术方案计算的是水质点在沙纹波高处或者临底高含沙层顶部的垂向运动速度,没有考虑水质点在波浪边界层内的速度分布,存在着较大的误差。
(3)在确定波浪半周期内可供悬沙颗粒沉降时段t时,现有的技术方案或者采用整个波浪半周期或者假定为半周期内水质点减速运动的时段,均为经验性的处理,存在较大的误差。
发明内容
本发明为解决现有的相位滞后效应的计算方法存在的没有考虑颗粒的部分动能由于颗粒之间的碰撞作用被消耗、没有考虑水质点在波浪边界层内的速度分布以及在确定波浪半周期内可供悬沙颗粒沉降时段时均采用经验性处理而导致存在着较大误差的问题,进而提供了一种计算波浪作用下悬沙输移相位滞后效应的方法。为此,本发明提出了如下的技术方案:
一种计算波浪作用下悬沙输移相位滞后效应的方法,包括:
通过将预定百分比的颗粒动能转换为势能确定在波浪半周期内所述颗粒的悬浮高度;
当波浪边界层内水质点时均垂向运动速度呈指数分布时,通过对所述指数分布沿深度进行积分平均,获得颗粒悬浮层内水质点的深度平均的垂向流速;
根据颗粒在沉降过程中受到颗粒之间的相互碰撞作用后的受阻碍沉速和所述颗粒悬浮层内水质点的深度平均的垂向流速计算获得波浪边界层内悬沙颗粒的沉速;
根据波浪半周期内悬浮颗粒的沉降时段、所述波浪边界层内悬沙颗粒的沉速和所述颗粒的悬浮高度计算获得波浪半周期的相位滞后效应参数。
由上述本发明提供的技术方案可以看出,通过假定波浪边界层内水质点时均垂向运动速度呈指数分布,再通过沿深度进行积分平均得到颗粒悬浮层内水质点的深度平均垂向流速,从而能够更为充分地反应水质点垂向流速对悬浮颗粒在波峰半周期和波谷半周期内沉降的影响,因此对颗粒沉速的修正就更为准确;通过将颗粒的动能转化成势能来计算其被悬浮的高度,由于考虑了颗粒之间的相互碰撞作用造成的能量耗散比例,因而计算结果更为准确;根据悬沙颗粒在波浪半周期内何时沉降、沉速大小和波浪边界层内水质点运动速度等因素的影响计算获得沉降时段,进一步提高了计算结果的准确性。
附图说明
图1为本发明的具体实施方式提供的计算波浪作用下悬沙输移相位滞后效应的方法的流程示意图。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
本具体实施方式提供了一种计算波浪作用下悬沙输移相位滞后效应的方法,如图1所示,包括:
步骤1,通过将预定百分比的颗粒动能转换为势能确定在波浪半周期内所述颗粒的悬浮高度。
具体的,在本具体实施方式中假设98%的颗粒动能转变成势能来确定颗粒的悬浮高度,即:
其中,Δi表示波浪半周期内颗粒被悬浮的高度,下角标i=c、t分别表示波峰半周期和波谷半周期,ui表示波浪波峰速度和波谷速度,s表示泥沙颗粒的密度与水流密度之比且等于2.65,g表示重力加速度。
并且,本具体实施方式采用如下公式计算波浪半周期内悬浮颗粒的沉降时段:
ti=Ti 0.74 (4)
步骤2,当波浪边界层内水质点时均垂向运动速度呈指数分布时,通过对所述指数分布沿深度进行积分平均,获得颗粒悬浮层内水质点的深度平均的垂向流速。
具体的,悬沙颗粒在减速时段的沉降会受到波浪边界层内水质点垂向流速的影响,在波峰半周期内垂向流速向下,会加强颗粒的沉降;相反,在波谷半周期内,垂向流速向上,会减缓颗粒的沉降,鉴于水质点垂向流速是一个随时间和高度变化的过程,本具体实施方式假设波浪边界层内水质点时均垂向流速的分布为:
其中,Wfi,max表示水质点在颗粒悬浮高度Δi处的最大垂向流速,z表示距床面的距离。通过对积分式(5)取深度平均值得到:
其中,Wf,max采用Kos’yan等(2007)提出的公式进行计算,h表示静止水深,H表示波高,Ti表示波浪半周期,λ表示波长。
步骤3,根据颗粒在沉降过程中受到颗粒之间的相互碰撞作用后的受阻碍沉速和所述颗粒悬浮层内水质点的深度平均的垂向流速计算获得波浪边界层内悬沙颗粒的沉速。
具体的,由于颗粒的沉降过程中还会受到颗粒之间的相互碰撞作用产生的阻碍作用,相应的沉速采用Richardson和Zaki(1954)公式计算:
其中,Whs表示受阻碍后的沉速,Ws表示无阻碍时的沉速,利用斯托克斯公式计算,c表示悬沙内平均体积比浓度,取值为0.365,cs表示底床最大体积比浓度,取值为0.65。
综合公式(6)和公式(7),则波浪边界层内悬沙颗粒的沉速计算公式为:
其中,ωsc、ωst分别表示波峰半周期和波谷半周期内颗粒的沉速。
步骤4,根据波浪半周期内悬浮颗粒的沉降时段、所述波浪边界层内悬沙颗粒的沉速和所述颗粒的悬浮高度计算获得波浪半周期的相位滞后效应参数。
具体的,波浪半周期的相位滞后效应参数Pi的计算表达式为:
其中,ti=Ti 0.74表示波浪半周期内悬浮颗粒的沉降时段。如果P值小于或等于1,则表示颗粒能够在半周期内及时沉降到床面上,如果P值大于1,则表示颗粒沉降所需要的时间要长于可供沉降的时段,因而会有部分悬沙不能及时沉降到床面上而是保持悬浮到下一个半周期内继续输移沉降,产生相位滞后效应。
下面通过具体的计算示例对本发明提出的计算波浪作用下悬沙输移相位滞后效应的方法进行说明:
有如下的波浪要素和泥沙粒径:H=1.50m,T=6.5s,Tc=2.51s,Tt=3.99s,h=3.5m,D50=0.138mm,Uc=1.55m/s,Ut=-0.83m/s。
(1)利用公式(3)计算可以得到波峰、波谷半周期内悬沙颗粒被悬浮的高度:Δc=0.073m、Δt=0.021m;
(2)利用公式(4)计算可以得到波峰、波谷半周期内悬浮颗粒的沉降时段:tc=1.98s、tt=2.78s;
(3)利用公式(6)计算可以得到波峰、波谷半周期边界层内水质点垂向运动平均流速:
(4)利用公式(7)计算可以得到颗粒的受阻碍沉速:Whs=2.77×10-4m/s;
(5)利用公式(8)计算可以得到波峰、波谷半周期内悬浮颗粒的沉速:ωsc=0.023m/s、ωst=0.022m/s;
(6)最后利用公式(9)计算可以得到波峰、波谷半周期内的相位滞后效应参数:Pc=1.63、Pt=0.35。
上述计算结果表明,波峰半周期内的泥沙输移具有相位滞后效应,因而部分泥沙会保持悬浮状态继续向波谷半周期内输移,造成向岸的净输沙减少。
本具体实施方式提供的计算波浪作用下悬沙输移相位滞后效应的方法具有如下的技术效果:
1.计算波浪边界层内水质点垂向运动速度对计算相位滞后相应参数十分关键。现有的技术方案仅是简单地计算沙纹波高处或者临底高含沙层顶部水质点的垂向运动速度。因而在对悬浮颗粒的沉速进行修正时就会产生较大的误差。
本具体实施方式提供的方法首先假定波浪边界层内水质点时均垂向运动速度呈指数分布(公式5),通过沿深度进行积分平均得到颗粒悬浮层内水质点的深度平均垂向运动速度(公式6),该时均、深度平均的速度能够更为充分地反映水质点垂向运动速度对悬浮颗粒在波峰半周期和波谷半周期内沉降的影响,并且计算简单。因而,对颗粒沉速的修正就更为准确,在准确计算相位滞后效应参数中起到关键作用。
2.如何准确计算悬沙颗粒的被悬浮高度对计算相位滞后相应参数也相当关键。现有的技术方案或者只是简单地假定该悬浮高度是沙纹波高或者临底高含沙层厚度的数倍,而具体的倍数需要依据实验数据进行率定,因而误差较大;或者假定颗粒运动的全部动能转化成势能来计算颗粒的被悬浮高度,因而也存在着一定的误差。
本具体实施方式提供的方法也采用将颗粒的动能转化成势能来计算其被悬浮的高度,但是考虑了颗粒之间的相互碰撞作用造成的能量耗散比例(公式3),因而计算结果更为准确,在准确计算相位滞后效应参数中起到重要作用。
3.如何确定在波浪半周期内悬沙颗粒沉降的时长对于准确计算相位滞后效应参数很关键。现有的技术方案一般都是假定整个波峰半周期和波谷半周期内悬沙颗粒都在沉降或者是经验性地给出某一时长作为沉降时段,因而误差较大。
本具体实施方式提供的方法通过充分考虑悬沙颗粒在波浪半周期内何时沉降、沉速大小和波浪边界层内水质点运动速度等因素的影响,给出了计算该沉降时段的公式(公式4),因而计算结果更为准确,在准确计算相位滞后效应参数中也起到重要作用。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明实施例揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (5)
1.一种计算波浪作用下悬沙输移相位滞后效应的方法,其特征在于,包括:
通过将预定百分比的颗粒动能转换为势能以确定在波浪半周期内所述颗粒的悬浮高度;
当波浪边界层内水质点时均垂向运动速度呈指数分布时,通过对所述指数分布沿深度进行积分平均,获得颗粒悬浮层内水质点的深度平均的垂向流速;
根据颗粒在沉降过程中受到颗粒之间的相互碰撞作用后的受阻碍沉速和所述颗粒悬浮层内水质点的深度平均的垂向流速修正获得波浪边界层内悬沙颗粒的沉速;
根据波浪半周期内悬浮颗粒的沉降时段、所述波浪边界层内悬沙颗粒的沉速和所述颗粒的悬浮高度计算获得波浪半周期的相位滞后效应参数。
2.根据权利要求1所述的计算波浪作用下悬沙输移相位滞后效应的方法,其特征在于,所述颗粒的悬浮高度通过以下公式计算获得:
其中,Δi表示波浪半周期内颗粒被悬浮的高度,下角标i=c、t分别表示波峰半周期和波谷半周期,ui表示波浪波峰速度和波谷速度,s表示泥沙颗粒的密度与水流密度之比,g表示重力加速度。
3.根据权利要求2所述的计算波浪作用下悬沙输移相位滞后效应的方法,其特征在于,所述波浪边界层内水质点时均垂向运动速度呈指数分布满足以下公式:
其中,Wfi,max表示在颗粒悬浮层内水质点在颗粒悬浮高度Δi处的最大垂向流速,z表示距床面的距离。
4.根据权利要求3所述的计算波浪作用下悬沙输移相位滞后效应的方法,其特征在于,所述在颗粒悬浮层内水质点的深度平均的垂向流速通过以下公式计算获得:
其中,h表示静止水深,H表示波高,Ti表示波浪半周期,λ表示波长。
5.根据权利要求4所述的计算波浪作用下悬沙输移相位滞后效应的方法,其特征在于,所述颗粒在沉降过程中受到颗粒之间的相互碰撞作用后的受阻碍沉速通过以下公式计算获得:
其中,Whs表示受阻碍后的沉速,Ws表示无阻碍时的沉速,能够利用斯托克斯公式计算,c表示悬沙内平均体积比浓度,取值为0.365,cs表示底床最大体积比浓度,取值为0.65。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410155863.6A CN103940580B (zh) | 2014-04-17 | 2014-04-17 | 一种计算波浪作用下悬沙输移相位滞后效应的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410155863.6A CN103940580B (zh) | 2014-04-17 | 2014-04-17 | 一种计算波浪作用下悬沙输移相位滞后效应的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103940580A CN103940580A (zh) | 2014-07-23 |
CN103940580B true CN103940580B (zh) | 2016-10-05 |
Family
ID=51188344
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410155863.6A Expired - Fee Related CN103940580B (zh) | 2014-04-17 | 2014-04-17 | 一种计算波浪作用下悬沙输移相位滞后效应的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103940580B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107702885B (zh) * | 2017-09-29 | 2019-08-30 | 中国水利水电科学研究院 | 床面颗粒起悬时间的提取方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4436049A (en) * | 1982-05-20 | 1984-03-13 | The United States Of America As Represented By The Secretary Of The Navy | Simulator for testing sonobuoys |
US4495809A (en) * | 1983-03-29 | 1985-01-29 | The United States Of America As Represented By The Secretary Of The Navy | Deep depth undex simulator |
CN1664252A (zh) * | 2005-04-21 | 2005-09-07 | 中国科学院力学研究所 | 一种海底输沙模拟方法 |
CN102359862A (zh) * | 2011-08-12 | 2012-02-22 | 河海大学 | 粉沙质和淤泥质海岸泥沙运动数值模拟方法 |
CN102841188A (zh) * | 2012-09-28 | 2012-12-26 | 河海大学 | 淤泥质海岸风天波能功率损耗与含沙量分布关系评判方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100890563B1 (ko) * | 2007-11-16 | 2009-03-25 | 삼성중공업 주식회사 | 선박의 조종 및 내항 성능 실험 장치 |
-
2014
- 2014-04-17 CN CN201410155863.6A patent/CN103940580B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4436049A (en) * | 1982-05-20 | 1984-03-13 | The United States Of America As Represented By The Secretary Of The Navy | Simulator for testing sonobuoys |
US4495809A (en) * | 1983-03-29 | 1985-01-29 | The United States Of America As Represented By The Secretary Of The Navy | Deep depth undex simulator |
CN1664252A (zh) * | 2005-04-21 | 2005-09-07 | 中国科学院力学研究所 | 一种海底输沙模拟方法 |
CN102359862A (zh) * | 2011-08-12 | 2012-02-22 | 河海大学 | 粉沙质和淤泥质海岸泥沙运动数值模拟方法 |
CN102841188A (zh) * | 2012-09-28 | 2012-12-26 | 河海大学 | 淤泥质海岸风天波能功率损耗与含沙量分布关系评判方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103940580A (zh) | 2014-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Garcia et al. | Experiments on the entrainment of sediment into suspension by a dense bottom current | |
Powell | Glacimarine processes at grounding-line fans and their growth to ice-contact deltas | |
Patacci et al. | Flow behavior of ponded turbidity currents | |
KNELLER et al. | Modelling the effects of salt-induced topography on deposition from turbidity currents | |
Omid et al. | Effects of bed-load movement on flow resistance over bed forms | |
Gladstone et al. | On the application of box models to particle-driven gravity currents | |
Börzsönyi et al. | Two scenarios for avalanche dynamics in inclined granular layers | |
CN103940580B (zh) | 一种计算波浪作用下悬沙输移相位滞后效应的方法 | |
Morris et al. | Changes in flow direction at a point caused by obstacles during passage of a density current | |
Smyth et al. | Numerical modelling of turbulent flow structures in a trough blowout | |
Zhong et al. | A bed-load function based on kinetic theory | |
Chang et al. | Numerical simulation of sediment particles released at the edge of the viscous sublayer in steady and oscillating turbulent boundary layers | |
Udo et al. | Application of Rouse’s Sediment Concentration Profile to Aeolian Transport: Is the suspension system for sand transport in air the same as that in water? | |
Hadeler et al. | Granular matter in a silo | |
Jensen et al. | Sediment transport and backfilling of trenches in oscillatory flow | |
Arora et al. | Resistance to Flow and Velocity Distribution in Rigid Boundary Channels Carrying Sediment‐Laden Flow | |
Guan et al. | Translational momentum of deformable submarine landslides off a slope | |
CN110188310B (zh) | 梯级水库多流态复合的排沙预测方法及装置 | |
Hurzeler et al. | Spreading model for a turbidity current with reversing buoyancy from a constant-volume release | |
Hibino et al. | Effect of a capping inversion on the stability of an Ekman boundary layer | |
Saxton et al. | Density segregation of minerals during high-velocity transport over a rough bed: implications for the formation of placers | |
Armitage et al. | A unit stream power model for the prediction of local scour in rivers | |
Takakuwa et al. | Structures of flow with particles saltating over the permeable fixed rough bed | |
Cheng et al. | Dynamic response of water flow and sediment transport off the Yellow River mouth to tides and waves in winter | |
Emms | A streamtube model of rotating turbidity currents |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20161005 Termination date: 20200417 |
|
CF01 | Termination of patent right due to non-payment of annual fee |