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

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

水利工程论文-耦合算法在幕墙式消浪结构性能研究中的应用.doc水利工程论文-耦合算法在幕墙式消浪结构性能研究中的应用.doc -- 2 元

宽屏显示 收藏 分享

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

水利工程论文耦合算法在幕墙式消浪结构性能研究中的应用摘要作者针对近岸波浪与结构物相互作用问题提出了一种耦合数值方法,即用时均化的二维雷诺平均的Navierstokes方程流体体积法模型表达内域流动,用一维Boussinesq方程表达外域流动,通过速度、压力和波面匹配边界条件实现两种数值模型的同步求解。耦合模型中的二维子模型能够较好地表达结构物附近流动的细部特征,包括漩涡结构一维子模型的计算效率很高,可通过延长其计算域以达到有效地避免二次反射波的影响。所建立的耦合数值模型被证实可应用于幕墙式消浪结构防波性能的研究。关键词幕墙式消浪结构耦合数值方法VOF方法Boussinesq方程以往的研究成果表明,迎浪面开孔的沉箱直立堤可以有效地减小反射波,但消浪室的宽度即开孔前墙和不透水后墙之间的宽度一般应达到当地波长的四分之一1。如果入射波为涌浪或者其他类型的长波,这意味着理想的消浪室宽度在实际工程上可能无法实现。最近,日本学者提出了一种能有效消减直立堤前反射波的新型结构幕墙式消浪结构2curtainwalleddissipater,其断面如图1所示。设在直立墙前的垂直屏障称为幕墙。幕墙至直立墙的距离B为消浪室的宽度。图1幕墙式消浪结构幕墙吃水深度用c表示。由于入射波引起的消浪室内水体振荡运动和幕墙下面的涡旋运动使得波能大量耗散,从而实现消减反射波浪的目的。这种新型消浪结构的主要优点是能够有效地减小消浪室的设计宽度。数值模拟是揭示幕墙式消浪结构水力学性能和消浪机理的有效手段。由于在幕墙和直立墙处产生的反射波在造波边界处可能形成二次反射,通常的方法需要在二次反射波传播到结构物之前停止计算。这意味着计算域的长度必须足够大。然而,在一个很长的立面二维计算域上全部采用粘流波浪数值模型做精细模拟一方面计算工作量很大,另一方面必要性也不充分。为此,本文提出了一个耦合求解策略,即将二维RANSVOF模型与一维Boussinesq方程模型耦合起来解决问题。在耦合模型中,一维子模型由于其计算效率很高,可以考虑足够长的计算域二维子模型则能够较好地反映流场的细部,包括粘性对流动的影响。1耦合模型的原理如图2所示,耦合模型是将整个计算域Ω划分为Ω1和Ω2两个子域。这两个子域通过一条公共的重叠带衔接起来。Ω1为包含幕墙和直立墙的近场,流动以二维紊流运动方程,即Reynolds方程RANS为控制方程,采用标准kε紊流模型封闭,并在近壁区应用壁面函数理论6自由水面的描述采用定义流体体积函数的方法5。图2耦合模型区域划分Ω2域内流动的控制方程采用色散性改进的Boussinesq方程4的一维形式,经差分离散后得到系数矩阵为三对角矩阵形式的代数方程组,采用追赶法快速求解。在耦合模型中,两个子模型RANS和Boussinesq各自独立求解,耦合的实现体现在重叠带上流动信息的匹配。为了便于耦合处理,Boussinesq方程和RANS方程均采用交错网格进行差分离散。其中,Boussinesq方程的离散参考了Madsen和Sorensen所用的格式4。RANS动量方程中时间项的离散格式为向前差分,粘性项的离散格式为二阶中心差分。为消除数值粘性的影响,动量方程中对流项的离散格式采用了三阶迎风差分格式3。差分方程的求解采用了SOLAVOF方法5。其基本思想是首先用前一时刻的流场计算结果代入动量方程的显式差分格式,求出当前时刻流场的近似值再通过对压力厨行迭代修正,使得连续方程在一定的精度条件下得以满足,对表面单元要求满足自由表面的动力学边界条件,即通过线性插值确定表面单元中心处的压力值在完成压力迭代后,再对速度进行校正,然后用校正后的速度值代入kε方程相应的差分格式求解紊动动能和紊动动能耗散率最后,应用施主与受主单元模型计算当前时刻的流体体积函数,确定流体自由表面的位置。由于动量方程、紊动动能方程和紊动动能耗散率方程对近壁区网格细密程度的要求不同,耗散率方程的要求最严,动量方程和动能方程的要求基本一致,为了既保证解的精度而又不致使网格划分太密,本文在固壁区附近采用了壁面函数方法6进行处理。即在壁面附近引入以下关系其中β0是常数,与壁面粗糙度有关,本文取β00.0005L是特征长度,计算中取为近壁区网格中心到壁面的距离。合理地设置匹配边界条件,使得内域流动和外域流动在匹配边界处光滑而连续地过渡,是保证子模型耦合的关键。本文的做法如图3所示。匹配边界ΓB是Ω2域的出流边界,Boussinesq模型执行每一时间步的计算之前需要预知该边界上的速度和波面值,由于边界ΓB同时又在Ω1域的计算节点上,于是ΓB边界的速度条件ΓB可利用Ω1域得到的流场信息表达如下,图3匹配边界附近的差分网格单元1式中u为水平速度,F为流体体积函数5,δy为垂向的网格步长,j为垂向网格节点编号,Jmax表示垂向网格节点的最大编号。给ΓB边界的波面赋予匹配条件时数值试验表明,需要利用连续方程反映的水位流量关系给出匹配条件其效果好于直接给定水位过程条件。这是因为通过水位和流量的相互调整,计算域内的反射波可得以减弱。因而,在实际计算中ΓB边界的波面可表示为利用Boussinesq方程推导过程中引入的近似展开关系,在波面函数η和深度平均速度导的结果,ΓV边界上水平速度u\,垂向速度v和压力p的匹配条件可按下面表达式给出耦合模型同步求解过程可简单概括如下。首先,考虑Ω2域左边界处的入射波条件,在Ω2域内执行Boussinesq模型,当Ω2域中接近匹配边界的节点上的水平速度值第一次达到103m/s量级时,开始在Ω1域内执行RANSVOF模型。在某一时间步n,执行Boussinesq模型所需要的匹配边界ΓB处的速度和波面条件按式1和式2给出Boussinesq模型在当前时间步的计算完成后,随即利用式3~5计算出匹配边界ΓV处的波面、速度和压力边界条件,并启动RANSVOF模型RANSVOF模型在当前时间步的计算完成后,即按照式1和2计算出匹配边界ΓB所在位置的速度和波面匹配条件。然后进入下一时间步的计算。2计算结果为了验证计算结果,参考已知的模型试验条件,取耦合模型为等水深d0.42m,计算域全长18.36m.其中,子区域Ω2和Ω1的长度分别为12.58m和5.78m.Ω1域中幕墙厚度q0.042m,幕墙下端为一45°尖角。幕墙吃水深度按相对吃水c/d0.3~0.7计算了多个工况。幕墙和直立墙均作为全反射边界处理。消浪室宽度B0.29m取为固定值。通过改变入射波的周期,消浪室的相对宽度在B/L0.08~0.16L为波长的范围内变化。计算网格的划分情况如下。Ω2域按等步长Δx0.042m划分为300个单元格Ω1域的垂向也按等步长Δz0.021m划分,但其水平方向是按变步长划分的,共划分了150个单元格。Ω1域中单元编号i1~120的各单元为Δxi0.042m,i121~125的各单元为Δxi0.035m,i126~130的各单元为Δxi0.028m,i131~150的各单元为Δxi0.021m.幕墙附近和消浪室区域的网格划分得比较细的目的是为了更有效地分辨幕墙下面涡
编号:201312122226366474    大小:26.92KB    格式: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人评价)
浏览:5次
奋斗不息上传于2013-12-12

官方联系方式

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

相关资源

相关资源

相关搜索

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