水利工程论文-黄河口泥沙异重流基本控制参数的数值试验.doc水利工程论文-黄河口泥沙异重流基本控制参数的数值试验.doc

收藏 分享

资源预览需要最新版本的Flash Player支持。
您尚未安装或版本过低,建议您

水利工程论文黄河口泥沙异重流基本控制参数的数值试验摘要为了定量地评价泥沙异重流基本控制参数对泥沙异重流流动的影响,本文设计了一系列数值试验来展示各种参数与泥沙异重流流动特征之间的定量关系。采用ADI法求解三维方程经垂向积分后的平面二维泥沙异重流方程。所选取的基本控制参数包括水下斜坡坡角,总阻力系数,水力挟带系数,泥沙饱和浓度和泥沙沉速。通过数值试验得出水下斜坡坡角Θ、泥沙饱和浓度ES越大,总阻力系数CD、泥沙沉速Ω越小,相应地泥沙异重流流速、流动厚度以及悬沙浓度越大。其中斜坡坡角对泥沙异重流影响最大。当泥沙异重流流动处于超临界流状态,水力挟带系数EW大于零,这时由于有一定水量的卷入,所以流动厚度沿程增加。关键词黄河口泥沙异重流控制参数数值试验ADI法泥沙异重流及形成的地质体在沉积学和石油地质中有重要意义。在黄河口,泥沙异重流是泥沙搬运的重要方式,是河口区沉积动力的重要过程,是认识河口三角洲发育演化的重要途径。为使研究深入一步,我们进行了泥沙异重流的数值模拟。近年来,国内对异重流数模研究较少,仅见方春明等1997关于泥沙异重流潜入时的立面二维数值模拟以及邱晨霞1995对盐水异重流的二维两层数值计算。1平面二维泥沙异重流控制方程及数值模拟目前,包含三个基本控制方程流体质量、动量守恒方程和泥沙质量守恒方程的泥沙异重流数值模拟被广泛应用,本文把一维方程推广到平面二维的情形,成为四方程模型,如下所示11234其中U,V分别为X,Y方向速度,ΘX,ΘY分别为X,Y方向坡度,EWX与EWY分别为X、Y方向的水力挟带系数,ΤX、ΤY为X方向与Y方向的剪切摩擦力,KX、KY分别为X方向与Y方向的悬沙扩散系数,T为时间,H为异重流厚度,S为深度平均的悬沙含量,ΡS,ΡW,ΡF分别为泥沙颗粒、水体和泥沙异重流的密度,ES为泥沙饱和浓度,Ω为泥沙颗粒的沉降速度,Α为泥沙运动饱和系数,G为重力加速度,F为柯氏力参量。在方程1中,水力挟带系数EWX与EWY是流动RICHARDSON数RIX、RIY的函数,可用若干经验式表达5。以X方向为例,它们是ASHIDAETAL1975EWX00015/RIX;PARKERETAL1987EWX0075/1718RIX2425;FUKUSHIMETAL1985EWX00015/00204RIX;以及俞维升1991EWX00034/RIX29等。在方程2、3中,剪切摩擦力ΤX、ΤY是作用在异重流上下界面摩擦力的总和,它们由下式所示。其中CDCF1Α′,被称为总阻力系数,范家骅等1980的CD值约0003,PARKERETAL1987的CD值在0001~0058之间,俞维升1991的CD值在0003~0004之间,可见CD值有很大的不确定性。在方程4中,泥沙饱和浓度ES是描述泥沙异重流流动特征参数的函数,对黄河口可用下式表示中国水科院1997公式2,3S≤15KG/M3,ES123V2/GH036V/Ω0331ΔH/H02S15KG/M3,ES97V2/GH001V/Ω0161ΔH/H022V为合流速,H为水深,Ω为泥沙沉速,ΔH为潮差。张青玉经验式ES983V2/H023定解条件边界条件陆地边界VN=0N为岸线的法线方向,S′0′表示陆地水边界HTHT,SS表示水边界;河口边界U15M/SS50KG/M3初始条件U0V0S0H05M采用有限差分法中的ADI法求解方程组1、2、3、4。事实上,我们应当同时求解描述泥沙异重流流动的方程,因为它们是一组相互耦合的方程组。然而同时求解相当复杂,我们便采取一种叫冻结系数法4的近似耦合法来求解。计算黄河口泥沙异重流时,泥沙异重流起点在黄河入海处。数值计算区域如图1所示,X、Y方向空间步长均取1/3千米,时间步长取24秒;计算区海底地形如图2所示;柯氏参数F=2ΩSINΦ,Ω取72722105S1,Φ取367;泥沙密度ΡS取265G/CM3,水体密度ΡW取1015G/CM3,泥沙异重流密度ΡF取1040G/CM3。图1泥沙异重流数值计算区域AREAOFNUMERICALSIMULATIONOFSEDIMENTDENSITYCURRENT图2计算区海底地形图MAPOFAREAOFNUMERICALSIMULATION2平面二维泥沙异重流的参数试验为了定量地评价基本控制参数对泥沙异重流流动的影响,本文设计了一系列数值试验或称敏感性试验来展示各种参数与泥沙异重流流动特征之间的定量关系。所选取的基本控制参数包括水下斜坡坡角ΘX、ΘY,总阻力系数CD,水力挟带系数EW,泥沙饱和浓度ES和泥沙沉速Ω,进行数值试验时它们各自的输入值见表1。以试验13为标准进行对比,试验13计算结果如图3。表1数值试验中基本控制参数的输入值INPUTVALUESOFDIFFERENTBASICCONTROLLINGPARAMETERS试验序号斜坡坡角ΘX~ΘY总阻力系数CD水力挟带系数EW饱和浓度ESKG/M3泥沙沉速ΩCM/S101~00001000036205~00001000036310~000010000364真实地形00010000365真实地形00050000366真实地形00100015/RI00367真实地形0050000368真实地形001000034/RI2900369真实地形00100公式1003610真实地形00100045711真实地形00100008612真实地形00100002013真实地形001000036指张青玉经验式,公式1指水科院经验式图3数值试验13的计算结果图本图的坐标与图1的坐标一致CALCULATINGRESULTSOFNUMERICALTESTNO1321水下斜坡坡角的影响水下斜坡坡角越大,有效重力作用越强,泥沙异重流加速越快,可以运行更远的距离,而且可以一直保持较高的含沙量。通过数值试验13、1、2、3可揭示坡角大小对泥沙异重流流动产生的显著影响。这些数值试验的控制参数除坡角不同外,其余参数均相同,见表1。流速变化图4显示坡角越大,流动速度越大。在河口以东2KM处图4,试验13的流速为080M/S,试验1、2、3的流速分别为141M/S,182M/S和267M/S。数值试验13如图3A,由于其输入真实海底地形,泥沙异重流从黄河口流出后,便向周围扩散,这是由于实际上口门处地形最高,无论向北向南,还是向东均有一定的坡度。流动厚度变化图5同样显示,坡角越大,其流动厚度越大。试验13,1,2,3的08M等厚度线所括范围占整个计算区域的百分比分别为250,666,750,833。悬沙含量变化图6显示,坡角越大,愈易保持高的悬沙含量,在河口以东2KM处,试验13,1,2,3的悬浮含量分别为185,223,316,379G/L。22总阻力系数CD的影响总阻力系数CD的大小主要与边界条件有关,对于不同边界不同性质的泥沙异重流,其总阻力系数CD的值很不相同。CD值大小对泥沙异重流影响很大,CD越大,泥沙异重流在运动过程中动量损失越大,越易发生沿程沉降而导致快速消亡。本文设计了数值试验4、5、13、7来定量分析CD对泥沙异重流流动的影响。CD值分别是0001、0005、001和005,依次增大。图4显示随着CD值的增大,流速逐渐减小,在河口以东2KM处,试验4的流速为136M/S,试验5,13的流速分别为105M/S,080M/S;试验7的CD值高达005,几乎是所量测到的CD值中的最大值,在如此大的阻力下,流速快速衰减,在2KM处流速已降至032M/S。流动厚度变化图5表明,流动厚度的衰减速度随CD值的增大而迅速增大,以08M等深线为例,在试验4中,08M等深线所括范围很大,占据整个计算区的3/4,而在试验13中,它的范围已大大缩减,仅占整个计算区域的1/4左右,试验7的范围更小,占整个区域不到1/8。流速和流动厚度的减小,必然导致悬沙含量随之减小,如图6所示,试验4、5、13、7在河口以东2KM处的悬沙含量分别是326,253,185和164KG/M3。23水力挟带系数EW的影响当泥沙异重流流动处于超临界流状态,水力挟带系数EW大于零,这时有一定水量的卷入;当泥沙异重流处于亚临界流状态,水力挟带系数EW等于零,甚至小于零出现负挟带。黄河口的泥沙异重流一般情况下,其流速小于15M/S,悬沙浓度小于60G/L,属低浓度泥沙异重流,通常处于亚临界状态,因此在本文的数值试验中,水力挟带系数EW大部分取零值,只有试验6和试验8中的EW分别取00015/RI和000034/RI29。三个试验的流速和悬沙含量基本类似,说明有无水力挟带项以及水力挟带系数的大小对流速和悬沙含量的分布影响不大;然而三个试验的流动厚度大相径庭图7,试验13的流动厚度沿程减小,而试验6,试验8的流动厚度沿程增大,EW值越大,卷入的水量便越多,试验8的EW略大于试验6的EW,因此,试验8流动厚度的增加略快于试验6。图4泥沙异重流流速M/S变化图河口以东2KM处的流速VARIATIONOFVELOCITYOFSEDIMENTDENSITYCURRENT图5泥沙异重流流动厚度变化图08M等厚度线所括范围与整个计算区域之比VARIATIONOFFLOWDEPTHOFSEDIMENTDENSITYCURRENT图6泥沙异重流悬沙浓度(G/L)变化图河口以东2KM处的浓度VARIATIONOFSEDIMENTCONCENTRATIONOFSEDIMENTDENSITYCURRENT图7EW值大小对泥沙异重流流动厚度M的影响注图中阿拉伯数字代表数值试验序号VARIATIONOFDENSITYCURRENTDEPTHINFLUENCEDBYEW24泥沙饱和浓度ES的影响泥沙饱和浓度ES是影响泥沙异重流数值模拟的关键因素,ES越大,泥沙异重流越易保持而不消亡。本文设计了试验9与试验13来定量分析ES不同对泥沙异重流所造成的影响。试验9的ES用水科院关于黄河口输沙力的公式来表达,试验13的ES用张青玉公式表达,除ES不同外,两个试验的其余参数均相同。流速变化图4清楚地显示,当ES值增大时,流速增大,口门以东2KM处的流速从056M/S增到080M/S。在流动厚度变化图5中,试验9的08M等值线范围仅占整个区域的1/10,而试验13为1/4左右,即流动厚度增加了25倍左右;如悬沙含量图6所示,试验9的悬沙含量在河口以东不到2KM处已降至15KG/M3,而试验13的悬沙含量在河口以东2KM处仍为185KG/M3,说明泥沙饱和浓度越小,泥沙异重流越易发生沿程沉降而导致其快速消亡。
编号:201312122239036847    类型:共享资源    大小:58.49KB    格式:DOC    上传时间:2013-12-12
  
2
关 键 词:
专业文献 建筑水利 精品文档 水利工程
  人人文库网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
关于本文
本文标题:水利工程论文-黄河口泥沙异重流基本控制参数的数值试验.doc
链接地址:http://www.renrendoc.com/p-206847.html

当前资源信息

4.0
 
(2人评价)
浏览:6次
wenku上传于2013-12-12

官方联系方式

客服手机:17625900360   
2:不支持迅雷下载,请使用浏览器下载   
3:不支持QQ浏览器下载,请用其他浏览器   
4:下载后的文档和图纸-无水印   
5:文档经过压缩,下载后原文更清晰   

精品推荐

相关阅读

人人文库
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们

网站客服QQ:2846424093    人人文库上传用户QQ群:460291265   

[email protected] 2016-2018  renrendoc.com 网站版权所有   南天在线技术支持

经营许可证编号:苏ICP备12009002号-5