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

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

   首页 人人文库网 > 资源分类 > DOC文档下载

水利工程论文-坡面降雨入渗产流规律的数值模拟研究.doc

  • 资源星级:
  • 资源大小:21.52KB   全文页数:16页
  • 资源格式: DOC        下载权限:注册会员/VIP会员
您还没有登陆,请先登录。登陆后即可下载此文档。
  合作网站登录: 微信快捷登录 支付宝快捷登录   QQ登录   微博登录
友情提示
2:本站资源不支持迅雷下载,请使用浏览器直接下载(不支持QQ浏览器)
3:本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰   

水利工程论文-坡面降雨入渗产流规律的数值模拟研究.doc

水利工程论文坡面降雨入渗产流规律的数值模拟研究摘要本文采用运动波理论和两次改进后的GreenAmpt入渗模型建立了坡面降雨入渗产流的动力学模型,并得到了实验资料的良好验证。运用该模型研究了简单坡面上降雨入渗产流的动力学规律,分析了雨强、土壤初始含水量、渗透系数、坡面阻力,以及坡长、坡度等因素对坡面产流过程的影响规律,得出了一些有益的结论。关键词入渗产流坡面动力学1概述雨水降落在坡面上将产生雨水的聚集并形成坡面水流。坡面水流是土壤水蚀过程的主要动力,搞清产流的动力学特点是进一步研究侵蚀过程规律的基础。坡面水流不同于一般明渠流动,其水深极浅一般只有几毫米,沿程不断有质量源和动量源加入,使其随时间和空间有较大的变化。且坡面流的坡度较一般河渠陡得多,边界条件也更为复杂。这些特点使得对坡面水流的研究有相当的难度。坡面产流研究已有很长历史,但对它的数学求解还只有三十多年。60年代后期Woolhiser和Ligget1967将运动波模型引入坡面水流研究,大大简化了计算工作,促进了研究的发展。运动波模型是从一维圣维南方程简化而来,其基本假设是水流的能坡和底坡相等,并借助Chezy阻力公式得到流量和水深的关系。Woolhiser和Ligget的研究结果表明在运动波波数k10时,运动波模型可以很好地描述坡面水流运动。而实际坡面流的运动波波数一般远大于10沈冰等,1996。因此,运动波近似是一种较好的数学描述方式。其后,又有对运动波理论的修正Ponce,1978,Govindaraju,1988,保留了水深的沿程变化项,相当于压力梯度,被称为扩散波模型。该模型扩展了适用的参数范围,但并无实质性改进,因此实际应用仍以运动波为主。也有使用完整圣维南方程求解实际问题的戚隆溪,1997。土壤入渗过程的研究也有很长历史,从1911年提出概念明确形式简单的GreenAmpt积水入渗模型开始,相继有Horton1940,Philip1957等模型出现,但GA模型仍以其简单的形式,明晰的物理概念,良好的扩展性和可信的应用效果受到广泛重视,特别是经过MeinLarson1973和Chu1978的两次改进,使其可应用于不均匀降雨的入渗计算,更使它成为最有效和应用最广泛的模型。在国内,GA模型尚未受到重视,Horton模型曾得到相当广泛的运用,但其参数的物理意义明显不如GA模型明晰。也有研究者使用更基本的土壤水分运动微分方程,但所需的参数更加难于获取,计算也更为复杂。本文工作旨在建立物理概念明晰的降雨入渗产流综合计算模式,并用以研究简单坡面的产流过程,分析各主要因素的影响和各主要因素的影响和各主要参量的变化规律。以期对坡面产流的动力学规律有清楚的认识。2计算模式坡面流运动十分复杂,目前主要采用运动波理论、扩散波或完整圣维南方程进行描述。正如前文所述,运动波近似理论在大多数情况下可以很好地描述坡面流运动过程,且计算简单。因此本文仍采用一维运动波理论,即坡面流基本方程为1此处第二式直接使用了水力学中熟知的Chezy公式和Manning公式。其中,x为沿坡面向下的坐标,t为时间s,h为水深m,q为单宽流量m2/s,p为降雨强度m/s,此处假设降雨方向垂直向下,i为入渗率m/s,S0为坡面坡度,S0sinθ,θ为坡面倾角,n为Manning糙率系数。土壤的入渗过程对坡面流的形成和流动过程影响很大,本文采用形式简单、物理概念明晰的GA入渗模型,其计算方程为idI/dtK1θSθiS/IIKtSθSθiln1I/SθSθi2其中K为土壤饱和导水率渗透系数m/s,θS为土壤饱和含水率,即有效孔隙率,θi为土壤初始含水率,S为土壤吸力m,I为累积入渗量m。经典的GreenAmpt模型是干土积水入渗模型,其前提是在整个入渗过程中地表始终有积水。MeinLarson1973年将其推广应用至降雨入渗的情况。设有稳定的雨强p,只有p大于土壤的入渗能力时,地表才能形成积水。而在降雨的初始阶段,全部降雨都渗入地下。由GA模型知,入渗率是随累积入渗量的增加而减小的。设想当累积入渗量达到某一值时,ip,此时开始积水,称此累积入渗量为Ip。因此由GA模型入渗公式可以导出开始积水时的Ip值IpθSθiS/p/K13开始积水时间由tpIp/p给出。因此整个过程的入渗率可表示为ipt≤tpiK1θSθiS/Ittp4式中的I为积水开始后的累积入渗量包含未积水时段的入渗量在内。由于不是由t0开始积水,I的计算须采用修正后的公式KttptSISθSθiln1I/SθSθittp5tS表示假设由t0开始积水,到入渗量IIp或ip时所需时间,可理解为一个虚拟时间,可计算如下KtSIpSθSθiln1Ip/SθSθi6改进的主要思想是将整个过程假设为从一开始就是积水入渗,这样该曲线在积水后部分相对于实际入渗曲线将向左平移tptS,将这条曲线向右平移tptS,再加上积水前的入渗强度等于降雨强度的关系,就得到真实的入渗过程。但稳定降雨在实际应用中远不能满足需要,Chu1978将MeinLarson1973改进的GA模型再作推广,使其可应用于变化的降雨过程。基本作法是,对每个计算时段将地表状态分为四种情况1.开始无积水,结束无积水2.开始无积水,结束有积水3.开始有积水,结束有积水4.开始有积水,结束无积水在每一时段开始,已知降雨总量与入渗总量,剩余总量。根据两个因子判断时段结束时是否有积水。若时段开始无积水,使用因子cu,若时段开始有积水,使用因子cp,其表达式为cuPtnRtn1KSM/iKcu0时段结束将积水,cu0时段结束仍有积水,cu0积水消失7其中M代表θSθi,Ptn代表tn时刻降雨总量,Rtn1代表tn1时刻剩余总量。可以证明,时段结束时积水与否与此两因子的正负等价。当iK时,始终无积水,不用此两因子判断。3模型求解模型的求解也分为两部分。运动波方程是一个非线性的对流型方程,求解采用了简单的一阶显式迎风格式。经过比较,一阶显式迎风格式对此问题的计算结果在幅值、相位、守恒性几个方面的综合效果较其他一些格式为优。我们曾使用过Preissmann四点偏心格式,它在计算此类坡面薄层水流运动波方程时耗散过大,表现为达到平衡产流的时间偏快。我们也曾试过中心差分格式和高阶迎风格式,结果发现标准的中心差分格式仍然和在普通的线性对流型方程中的表现一样

注意事项

本文(水利工程论文-坡面降雨入渗产流规律的数值模拟研究.doc)为本站会员(奋斗不息)主动上传,人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知人人文库网([email protected]),我们立即给予删除!

温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。

copyright@ 2015-2017 人人文库网网站版权所有
苏ICP备12009002号-5