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

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

水利工程论文-异重流潜入运动的剖面二维数值模拟.doc水利工程论文-异重流潜入运动的剖面二维数值模拟.doc -- 2 元

宽屏显示 收藏 分享

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

水利工程论文异重流潜入运动的剖面二维数值模拟摘要本文采用变密度流基本方程和混合有限分析法,求解了突扩边界下的异重流潜入运动。通过对数值计算结果分析,讨论了进口密度弗汝德数Fre对潜入运动的影响,加深了对异重流潜入规律的了解。关键词异重流潜入运动数值计算混合有限分析法1前言异重流是自然界中常见的一种流动现象,如水库中的浑水异重流,河口地区出现的盐水楔,工业中遇到的温差异重流,底部突然下沉的小型沉沙池中的泥沙异重流等,在水库、船闸引航道、河港、河口、沉沙池等实际工程中,浑水异重流的淤积问题往往是必须考虑的问题1。由于异重流潜入发生处的边界条件复杂,具有与一般明渠流不同的水流泥沙运动规律,且潜入处的水力泥沙因素与异重流的发生、发展及运行状况密切相关2,历来广为研究者们所关注。用一维的方法,再加上一些经验关系解异重流运动,无疑具有其实际意义,但异重流运动具有较强的三维性,不是一维流动理论基础所能较好解决的,并且在用一维方法求解异重流运动时,常常碰到诸如上下水层交界面位置及交界面上水体之间的动量和质量交换系数等难以确定的困难。近年来,逐步采用数值模拟技术对异重流潜入运动进行研究,如YoonT.H。采用流函数涡量法对异重流的潜入运动进行了模拟和研究3方春明、韩其为用立面二维对异重流潜入流动进行数值模拟,探讨了影响异重流潜入条件的因素4。但由于问题本身的复杂性,许多问题仍有待进一步研究。考虑到异重流潜入运动的过程中,不仅具有初始动量,而且还受到因密度差引起的浮力作用,与静止环境中的浮力射流运动有相似之处,不同之处在于异重流潜入发生处的边界条件复杂,因此,本文将异重流的潜入运动视为一种具有复杂边界的浮力射流运动,从NavierStokes方程出发,引入满足布辛涅斯克假定的一组变密度方程,并采用混合有限分析法对该运动进行一系列的数值模拟,通过对数值计算结果进行讨论和分析,进一步加深对该流动的认识。2数学模型图1二维异重流潜入运动的示意图Schematicdiagramof2Ddensitycurrentplungingflow图1所示密度为ρ0Δρ的浑水从水平河道流入密度为ρ0的清水矩形水库中,由于密度分布不均匀,流场中就出现了浮力,在几何形状突然扩大的流道中,形成了有浮力回流的流动。图中x、z分别表示水平和竖直方向,进口浑水流速u0满足抛物线分布。2.1控制方程该模型遵循以下假定忽略泥沙与水流之间的速度差异根据布辛涅斯克近似,仅重力项计及密度变化认为泥沙在水平和竖直方向上扩散系数相同仅对异重流潜入段进行模拟,此时浮力作用较大,且计算区间不太长,因此在这里只考虑异重流的对流和扩散,不考虑其沿程沉降4。为使计算和讨论具有一般性,以水库初始水深H、进口平均流速u0为特征尺度,引入无量纲量Xx/H,Zz/H,Uu/u0,Ww/u0,Ttu0/H,Pp/ρu20则控制方程的无量纲形式为连续方程1运动方程23浓度方程4式中Γρρ0/ρ0为浑水的有效密度,ReU0H/v为雷诺数,FrU0/为弗汝德数,Scv/D为斯密特数,v、D分别为流体的运动粘性系数和泥沙扩散系数。2.2初始和边界条件2.2.1初始条件初始时刻,水库水体静止,进口断面的浑水流速满足抛物线分布,且认为此断面的密度分布均匀初始压力按静水压力给出PX,Z1Z。2.2.2边界条件研究的流动区域共有进口、出口、自由表面、上游固壁及底部五个不同的边界,如图1所示。进口浑水流速满足抛物线分布,断面密度均匀分布,即U7.565ZZHH0/H,W0,ΓΔρ/ρ0出口认为流动已充分发展,则壁面给定无滑移边界条件,即U0,W0,上游或底部自由表面对自由面作刚盖假定,并采用对称平面条件,即。3计算方法3.1求解区域及网格划分该流动问题是开边界问题中的一种,当出口采用充分发展假定时,求解区域应保持足够的长度,以便在我们感兴趣的区域准确地模拟出流体的运动规律。依文献5,本文计算区域长度取为L20。因水库进口附近的流动是主要的,远离进口的流动是渐变的,故在数值计算中,X方向上采用非均匀网格,进口附近的网格密,越远离进口,网格越疏在Z方向上,除了0、1附近的网格密些(与边界条件有关),其余部分采用均匀网格。依控制方程及边界区域等特点,采用原始变量法求解为便,故采用交错网格来离散并求解整个流场和密度场6。计算时,动量方程和浓度方程的控制精度为107,连续方程的控制精度为104。3.2计算格式采用混合有限分析法的五点格式进行离散。在局部单元格上,将非线性方程(2)~(4)线性化,写成下列一般形式5则各方程的系数可归纳为表1的形式。方程5的混合有限分析解为6表1控制方程系数CoefficientsofgoverningequationsφE2A2BSUReReUReWWReReUReWRe1/Fr2ΓΓScReScReUScReW0式中Ci1,j、Ci1,j、Ci,j1、Ci,j1和Cp分别为混合有限分析法系数,具体表达见文献7。其中,SEφtn为第n次的迭代值,Sn用中心差分计算,φnt用向前差分表示,即npφijn1φnij/Δt,将i1,j、i1,j、i,j1、i,j1依此用1、2、3、4表示,则式(6)可写为7对整个计算区域,由式(2)~(4)形成代数方程组,它们都是对角占优的,可用追赶法求解。4计算结果和分析4.1Fre对潜入形态的影响异重流在潜入水库的过程中,不仅具有初始动量,而且还受到因密度差而引起的浮力作用,流动是非定常和变密度的,此时影响流动的一个重要的特征数是进口密度弗汝德数Fre,具体表达形式为8式中u0为进口处浑水的特征流速,L为特征长度,在本文中,u0等于进口处浑水的平均流速,L等于进口水深H0,ρ为进口处浑水的密度,ρ0为水库的清水密度。Fre反映了进口水流惯性力作用与浮力作用之比当其它条件不变时,Fre越小,说明浑水惯性力越小,浮力作用越大反之亦然。为研究异重流潜入形态随Fre的变化规律,本文在进口单宽流量一定的情况下,通过改变进口处浑水的密度ρ,分别计算了Fre0.25,0.501,1.00,1.58时异重流各特征量的变化情况,并绘出了进口附近区域在T0.8(T为无量纲化时间,Ttu0/h)时刻相应的有效密度等值线图(Γ/Γ0)和速度矢量图(图2和图3),由图可看出,Fre对异重流潜入流动特性的影响是较大的。图2是不同Fre的相对有效密度等值线图当Fre0.25时,负浮力占主导地位,浑水入库后不久便下沉,在进口附近形成一股潜流当Fre0.501,1.00,1.58时,入库浑水的初始动量作用逐渐增大,负浮力失去其主导作用,浑水继续向前水平运动,直到水库清水对其阻力与维持其运动的力平衡为止,此时因负浮力作用浑水开始潜入,该位置通常也称为潜入点。从该图还可以看出,Fre越大,浑水与水库水体的混合越充分,在水库中的扩散程度越大,潜入点距进口的距离越长,越不容易形成异重流反之,则越容易形成异重流。图3是相应的速度矢量随Fre的变化图。该图反映了在进口下游一定区域内会形成较大的回流旋涡区,随着Fre的改变,回流旋涡区的形状、大小和中心位置相应发生改变。图2异重流潜入时的密度场Γ/Γn随Fre的变化T0.8DensityfieldΓ/ΓnwiththevariationofFreT0。8inthecaseofdensitycurrentplunging当Fre较小时,浑水在进口附近区域很快吸附近岸,在水库中形成一个逆时针方向的旋涡,且该回流旋涡区范围较大,密度场的范围相应地也较大随着Fre的增大,在进口附近将出现一顺时针方向的旋涡,同时逆时针方向的回流旋涡区范围将随之减小,密度场的范围相应地也减小当Fre增大到一定程度后,进口附近区域仅出现顺时针方向的回流旋涡,进口附近水流的掺混作用剧烈。图3异重流潜入时的流场随Fre的变化T0.8VelocityfieldwiththevariationofFreT0.8inthecaseofdensitycurrentplunging图4Xp/H与Fre之间的关系ReleationbetweenXpandFre为进一步研究异重流潜入形态与Fre之间的关系,我们还利用数值计算的结果,在图4中点绘了潜入点与进口之间的无量纲距离Xp/H与Fre之间的关系。该图清楚地反映了两者之间的关系Fre1.5时,随着Fre增大,Xp/H迅速增大,进口附近的流动主要受浑水在进口处的初始动量所控制,此时流动为射流型。在卷流型和射流型之间存在一个过渡区,过渡区的流动不仅
编号:201312122157385541    大小:49.55KB    格式: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人评价)
浏览:1次
奋斗不息上传于2013-12-12

官方联系方式

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

相关资源

相关资源

相关搜索

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