会员注册 | 登录 | 微信快捷登录 支付宝快捷登录 QQ登录 微博登录 | 帮助中心 人人文库renrendoc.com美如初恋!
站内搜索 百度文库

热门搜索: 直缝焊接机 矿井提升机 循环球式转向器图纸 机器人手爪发展史 管道机器人dwg 动平衡试验台设计

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

宽屏显示 收藏 分享

页面加载中... ... 广告 0 秒后退出

资源预览需要最新版本的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方向为例,它们是Ashidaetal.1975Ewx0.0015/RixParkeretal.1987Ewx0.075/1718Rix2.42.5Fukushimetal.1985Ewx0.0015/0.0204Rix以及俞维升1991Ewx0.0034/Rix2.9等。在方程2、3中,剪切摩擦力τx、τy是作用在异重流上下界面摩擦力的总和,它们由下式所示。其中cdcf1α′,被称为总阻力系数,范家骅等1980的cd值约0.003,Parkeretal.1987的cd值在0.001~0.058之间,俞维升1991的cd值在0.003~0.004之间,可见cd值有很大的不确定性。在方程4中,泥沙饱和浓度Es是描述泥沙异重流流动特征参数的函数,对黄河口可用下式表示中国水科院1997公式2,3s≤15kg/m3,Es123V2/gh0.36V/ω0.331Δh/h0.2s15kg/m3,Es9.7V2/gh0.01V/ω0.161Δh/h0.22V为合流速,h为水深,ω为泥沙沉速,Δh为潮差。张青玉经验式Es9.83V2/h0.23定解条件边界条件陆地边界Vn=0n为岸线的法线方向,s′0′表示陆地水边界htht,ss表示水边界河口边界u1.5m/ss50kg/m3初始条件u0v0s0h0.5m采用有限差分法中的ADI法求解方程组1、2、3、4。事实上,我们应当同时求解描述泥沙异重流流动的方程,因为它们是一组相互耦合的方程组。然而同时求解相当复杂,我们便采取一种叫冻结系数法4的近似耦合法来求解。计算黄河口泥沙异重流时,泥沙异重流起点在黄河入海处。数值计算区域如图1所示,x、y方向空间步长均取1/3千米,时间步长取24秒计算区海底地形如图2所示柯氏参数f=2Ωsinφ,Ω取7.2722105s1,φ取36.7°泥沙密度ρs取2.65g/cm3,水体密度ρw取1.015g/cm3,泥沙异重流密度ρf取1.040g/cm3。图1泥沙异重流数值计算区域Areaofnumericalsimulationofsedimentdensitycurrent图2计算区海底地形图Mapofareaofnumericalsimulation2平面二维泥沙异重流的参数试验为了定量地评价基本控制参数对泥沙异重流流动的影响,本文设计了一系列数值试验或称敏感性试验来展示各种参数与泥沙异重流流动特征之间的定量关系。所选取的基本控制参数包括水下斜坡坡角θx、θy,总阻力系数cd,水力挟带系数Ew,泥沙饱和浓度Es和泥沙沉速ω,进行数值试验时它们各自的输入值见表1。以试验13为标准进行对比,试验13计算结果如图3。表1数值试验中基本控制参数的输入值Inputvaluesofdifferentbasiccontrollingparameters试验序号斜坡坡角θx~θy°总阻力系数cd水力挟带系数Ew饱和浓度Eskg/m3泥沙沉速ωcm/s10.1~0.00.010.00.03620.5~0.00.010.00.03631.0~0.00.010.00.0364真实地形0.0010.00.0365真实地形0.0050.00.0366真实地形0.010.0015/Ri0.0367真实地形0.050.00.0368真实地形0.010.00034/Ri2.90.0369真实地形0.010.0公式10.03610真实地形0.010.00.45711真实地形0.010.00.08612真实地形0.010.00.02013真实地形0.010.00.036指张青玉经验式,公式1指水科院经验式图3数值试验13的计算结果图本图的坐标与图1的坐标一致CalculatingresultsofnumericaltestNo.132.1水下斜坡坡角的影响水下斜坡坡角越大,有效重力作用越强,泥沙异重流加速越快,可以运行更远的距离,而且可以一直保持较高的含沙量。通过数值试验13、1、2、3可揭示坡角大小对泥沙异重流流动产生的显著影响。这些数值试验的控制参数除坡角不同外,其余参数均相同,见表1。流速变化图4显示坡角越大,流动速度越大。在河口以东2km处图4,试验13的流速为0.80m/s,试验1、2、3的流速分别为1.41m/s,1.82m/s和2.67m/s。数值试验13如图3a,由于其输入真实海底地形,泥沙异重流从黄河口流出后,便向周围扩散,这是由于实际上口门处地形最高,无论向北向南,还是向东均有一定的坡度。流动厚度变化图5同样显示,坡角越大,其流动厚度越大。试验13,1,2,3的0.8m等厚度线所括范围占整个计算区域的百分比分别为25.0,66.6,75.0,83.3。悬沙含量变化图6显示,坡角越大,愈易保持高的悬沙含量,在河口以东2km处,试验13,1,2,3的悬浮含量分别为18.5,22.3,31.6,37.9g/l。2.2总阻力系数cd的影响总阻力系数cd的大小主要与边界条件有关,对于不同边界不同性质的泥沙异重流,其总阻力系数cd的值很不相同。cd值大小对泥沙异重流影响很大,cd越大,泥沙异重流在运动过程中动量损失越大,越易发生沿程沉降而导致快速消亡。本文设计了数值试验4、5、13、7来定量分析cd对泥沙异重流流动的影响。cd值分别是0.001、0.005、0.01和0.05,依次增大。图4显示随着cd值的增大,流速逐渐减小,在河口以东2km处,试验4的流速为1.36m/s,试验5,13的流速分别为1.05m/s,0.80m/s试验7的cd值高达0.05,几乎是所量测到的cd值中的最大值,在如此大的阻力下,流速快速衰减,在2km处流速已降至0.32m/s。流动厚度变化图5表明,流动厚度的衰减速度随cd值的增大而迅速增大,以0.8m等深线为例,在试验4中,0.8m等深线所括范围很大,占据整个计算区的3/4,而在试验13中,它的范围已大大缩减,仅占整个计算区域的1/4左右,试验7的范围更小,占整个区域不到1/8。流速和流动厚度的减小,必然导致悬沙含量随之减小,如图6所示,试验4、5、13、7在河口以东2km处的悬沙含量分别是32.6,25.3,18.5和16.4kg/m3。2.3水力挟带系数Ew的影响当泥沙异重流流动处于超临界流状态,水力挟带系数Ew大于零,这时有一定水量的卷入当泥沙异重流处于亚临界流状态,水力挟带系数Ew等于零,甚至小于零出现负挟带。黄河口的泥沙异重流一般情况下,其流速小于1.5m/s,悬沙浓度小于60g/l,属低浓度泥沙异重流,通常处于亚临界状态,因此在本文的数值试验中,水力挟带系数Ew大部分取零值,只有试验6和试验8中的Ew分别取0.0015/Ri和0.00034/Ri2.9。三个试验的流速和悬沙含量基本类似,说明有无水力挟带项以及水力挟带系数的大小对流速和悬沙含量的分布影响不大然而三个试验的流动厚度大相径庭图7,试验13的流动厚度沿程减小,而试验6,试验8的流动厚度沿程增大,Ew值越大,卷入的水量便越多,试验8的Ew略大于试验6的Ew,因此,试验8流动厚度的增加略快于试验6。图4泥沙异重流流速m/s变化图河口以东2km处的流速Variationofvelocityofsedimentdensitycurrent图5泥沙异重流流动厚度变化图0.8m等厚度线所括范围与整个计算区域之比Variationofflowdepthofsedimentdensitycurrent图6泥沙异重流悬沙浓度(g/l)变化图河口以东2km处的浓度Variationofsedimentconcentrationofsedimentdensitycurrent图7Ew值大小对泥沙异重流流动厚度m的影响注图中阿拉伯数字代表数值试验序号VariationofdensitycurrentdepthinfluencedbyEw2.4泥沙饱和浓度Es的影响泥沙饱和浓度Es是影响泥沙异重流数值模拟的关键因素,Es越大,泥沙异重流越易保持而不消亡。本文设计了试验9与试验13来定量分析Es不同对泥沙异重流所造成的影响。试验9的Es用水科院关于黄河口输沙力的公式来表达,试验13的Es用张青玉公式表达,除Es不同外,两个试验的其余参数均相同。流速变化图4清楚地显示,当Es值增大时,流速增大,口门以东2km处的流速从0.56m/s增到0.80m/s。在流动厚度变化图5中,试验9的0.8m等值线范围仅占整个区域的1/10,而试验13为1/4左右,即流动厚度增加了2.5倍左右如悬沙含量图6所示,试验9的悬沙含量在河口以东不到2km处已降至15kg/m3,而试验13的悬沙含量在河口以东2km处仍为18.5kg/m3,说明泥沙饱和浓度越小,泥沙异重流越易发生沿程沉降而导致其快速消亡。
编号:201312122239036847    大小:58.49KB    格式:DOC    上传时间:2013-12-12
  【编辑】
2
关 键 词:
专业文献 建筑水利 精品文档 水利工程
温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 人人文库网仅提供交流平台,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
  人人文库网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

暂无评论,赶快抢占沙发吧。

当前资源信息

4.0
 
(2人评价)
浏览:6次
奋斗不息上传于2013-12-12

官方联系方式

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

相关资源

相关资源

相关搜索

专业文献   建筑水利   精品文档   水利工程  
关于我们 - 网站声明 - 网站地图 - 友情链接 - 网站客服客服 - 联系我们
copyright@ 2015-2017 人人文库网网站版权所有
苏ICP备12009002号-5