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

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

水利工程论文-沉沙池立面流场数值模拟研究.doc水利工程论文-沉沙池立面流场数值模拟研究.doc -- 2 元

宽屏显示 收藏 分享

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

水利工程论文沉沙池立面流场数值模拟研究摘要采用κε紊流模型,对沉沙池立面流厨行了数值模拟。在验证模型正确的基础上,以大禹渡沉沙池为例,分析了沉沙池中水流的运动规律,为进一步研究泥沙在沉沙池中的运动规律提供了前提,最终为沉沙池的工程设计和运用管理提供科学依据。关键词沉沙池κε模型立面流场大禹渡1前言在多沙河流上修建引水工程,为了减少进入引水渠的泥沙,保证引水质量,往往需要在渠首设置沉沙池沉淀大部分泥沙,防止或减轻引水渠的淤积以及泥沙对水轮机、水泵等过流部件的磨损,防止粗颗粒泥沙进入农田,引起农田沙化。对于沉沙池的研究,我国科研工作者自50年代以来通过物理模型试验进行了大量的研究工作1~3,在沉沙池的结构与形式方面取得了许多成果,但在沉沙池计算理论方面,国内目前大多数还是将沉沙池水流作为一维流或二维均匀流处理,计算泥沙沉降仍采用传统的准静水沉降法、非饱和输沙等经验、半经验公式。而国外自70年代开始就已将各种紊流模型应用到沉沙池计算当中,提出了许多数值方法来模拟沉沙池中的水流泥沙运动规律。国内在这方面的研究尚不多见,因此有必要应用紊流理论来研究沉沙池中的实际水流运动,为沉沙池的工程设计和运用管理提供科学依据和理论指导。2数学模型2.1基本控制方程对沉沙池流场的研究,其最终目的是为了研究泥沙在沉沙池中的垂线分布和沿程淤积情况,在此基础上按照工程要求设计沉沙池的合理尺寸。因此可将沉沙池水流简化为立面二维水流进行研究,水流基本控制方程包括连续性方程1动量方程23其中(4)u、v分别为沉沙池水流方向(x)和水深方向(y)的流速分量,vt为紊动粘性系数,k、ε分别为紊动能及其耗散率,它们通过求解以下输运方程得到k方程(5)ε方程(6)式中,方程(1)~(6)组成了求解沉沙池水流的封闭方程组,模型中五个常用参数的取值见表1。2.2数值方法对以上方程组采用控制体积法进行离散,对流项采用上风差分格式,使计算结果不致发散,划分网格时采用交错网格技术,以避免棋盘格式分布的压力场或流速场,最后,采用TDMA与高斯赛德尔迭代相结合的方法求解非线性方程组。为了避免非线性方程组在求解过程中发散,对各因变量及压力实行欠松弛迭代,以保证解的稳定性和收敛性。表1经验常数值Valueofempiricalfactorscμcε1cε2σkσε0.091.441.921.01.32.3边界条件沉沙池立面流场的边界条件包括(1)进口断面边界条件(2)出口断面边界条件(3)自由水面边界条件和(4)固壁边界条件。分述如下。(1)进口断面边界条件沉沙池进口断面的水流按均匀流计算其行进流速,即水流方向流速U按明渠均匀流的有关公式求解,水深方向流速V0。k、ε分别采用下式计算其参考值,均假设在水深方向均匀分布4。k0.2u278其中,混掺长度lmcμ0.5Hin9Hin为池首进口水深。(2)出口断面边界条件沉沙池的出口往往为一溢流堰,可按堰流计算堰顶平均流速,作为U的出流边界,其余物理量则按自由出流条件给出,即10(3)自由水面边界条件自由水面近似采用刚盖假定,包括①压强采用相对压强,即P0②V0③U、k在自由水面的法向梯度为0,即④ε采用下式计算4εk1.5/0.43H。(4)固壁边界条件沉沙池的固体边界一般包括沉沙池底部、水流经闸门后突然扩大断面的跌坎壁及尾部溢流堰墙三部分,如图3所示。可采用壁面定律模拟粘性底层与紊流区交界面处的流速、紊动能及其耗散率5。①流速壁面条件ures/uτ1/κlnyE11式中ures为平行于壁面的流速分量,uτ为摩阻流速,yyuτ/v是无因次量,y为计算点到壁面的距离,κ为卡门常数,取为κ0.4,E为表征糙率参数,对水力光滑壁面,可取E9.0。②紊动能k的壁面边界条件12③耗散率ε壁面边界条件133模型验证本文计算程序采用FORTRAN77语言编制,对文献6中描述的算例进行了验证计算。该沉沙池有关参数及计算域分别如表2及图1所示。其中L为沉沙池长度,d0为池首闸门开启高度,H为池中水深,q为单宽流量,Uin=q/d0,表示为进口平均流速。图1沉沙池计算域示意图Calculatingsketchofsettlingbasin表2沉沙池尺寸及水力参数DimensionandhydraulicfactorsofsettlingbasinL(cm)H(cm)q(cm2/s)Uin(cm/s)d0(cm)73.011.9109.421.885.0图2流速U垂线分布情况对比ComparisonofverticaldistributionofvelocityU图3大禹渡沉沙池计算域概化剖面图CalculatingsketchofDayudusettingbasin利用上述建立的κε紊流模型对该沉沙池流厨行模拟计算,计算成果与前人成果及实测资料对比如图2所示。图2中,Imam等人采用旋度流函数法模拟沉沙池立面流场7,AbdelGawad等则采取在池首回流区假定流速分布函数进行线性积分及在尾部采用有限元相结合的方法求解水流流场6,与作者采用原始变量法求解水流流场的结果相比,U流速的垂线分布都很相近,且与实测资料吻合较好,这充分证明了以上所建模型的正确性和合理性。4模型应用在验证模型正确的基础上,利用该模型对山西省芮城大禹渡矩形沉沙池立面流厨行模拟计算,得到该沉沙池中流速、紊动能及其耗散率的垂线分布和沿程变化规律8。4.1概况大禹渡沉沙池设在大禹渡电灌站的一级站和二级站之间,总长234m,设计水深5~9.7m,计算时取水深为7m,沉沙池末端设一长192m的溢流堰,概化后的大禹渡沉沙池剖面图如图3所示。4.2计算成果分析取沉沙池的首部、中部及尾部三个典型断面,分析各个断面的流速、紊动能及耗散率的垂线分布,从而揭示沉沙池水流的运动规律。图4为各断面水流方向流速U的垂线分布。分析图4可以看出,在首部断面,沉沙池底部水流方向与池中水流方向相反,说明水流经闸门后断面突然扩大,在沉沙池首部产生局部回流,此处的水流紊动最为剧烈。至沉沙池中部,水流渐趋稳定,流速垂线分布与紊流流速分布规律相符合,即可以用对数律来表示沉沙池中\U\流速的垂线分布。图4流速U垂线分布VerticaldistributionofvelocityU图5流速V垂线分布VerticaldistributionofvelocityV图6紊动能κ垂线分布Verticaldistributionofturbulenceenergyκ
编号:201312122210495978    大小:49.51KB    格式: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人评价)
浏览:11次
奋斗不息上传于2013-12-12

官方联系方式

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

相关资源

相关资源

相关搜索

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