基于fluen软件的二维数值波式生成方法.docxVIP

基于fluen软件的二维数值波式生成方法 0 维数值波浪浚合模型 fluent是一种应用于商业软件的数值造波和伪波模拟的软件。作为一种通用的计算水流软件，fluent没有模拟波浪的模块。建立值波水槽，用计算机模拟物模，进行试验，降低了非接触点波测量的成本，并减少了测量结果。避免由于传感器尺寸和模型变形等因素对流场的影响，从而获得更详细的流动信息。中国的波浪与水工建筑物之间的相互作用研究仍停留在建模过程中。 本文利用FLUENT的二次开发功能建立二维数值波浪水槽模型,可有效模拟微小振幅线性波.通过模拟不同参数的线性波,对两种造波方法的优劣性进行分析,给出两种造波方法的优缺点及适用性,为工程应用中合理选用造波方法提供理论依据. 1 数学模型 1.1 流体力学模型 二维数值波浪水槽的控制方程采用以速度和压力为变量的不可压黏性流体的二维N-S方程和连续性方程.坐标系为欧拉直角坐标系,水平方向为X轴,右向为正;垂直方向为Z轴,向上为正.u、w分别为X、Z方向的速度分量,ρ为流体密度,p为流体压力,υ为流体的运动学黏性系数. 连续方程: ?u?x+?w?z=0(1)?u?x+?w?z=0(1) 流体区域和阻尼层消波段内的统一形式的动量方程为: ?u?t+u?u?x+w?u?z=-1ρ?p?x+υ(?2u?x2+?2u?z2)-μ(x)u?w?t+u?w?x+w?w?z=g-1ρ?p?z+υ(?2w?x2+?2w?z2)-μ(x)w(2)?u?t+u?u?x+w?u?z=?1ρ?p?x+υ(?2u?x2+?2u?z2)?μ(x)u?w?t+u?w?x+w?w?z=g?1ρ?p?z+υ(?2w?x2+?2w?z2)?μ(x)w(2) 其中,μ(x)为波浪消波系数.对流体区域μ(x)=0;对阻尼层消波段,μ(x)=K(x-x0)/(xl-x0).K=4.0. 1.2 波浪场的速度场 对于平衡位置在原点,造波板冲程为X0,角频率为ω的推板式造波机,其推板作简谐运动的速度: u(t)=X0ω2cosωt(3)u(t)=X0ω2cosωt(3) 推板运动产生一定周期和波高的波浪, 并向外传播.根据势流理论,在水深为h的波浪水槽中,距离造波板x处产生的波面可表示为 η(x,t)=2X0ω2coshkh?sinhkhkg(sinh2kh+2kh)?cos(kx-ωt)(4) 另外,由线性波浪理论可知,波浪场速度势函数为 ?(x,z,t)=gAω?coshk(z+h)coshkh?sin(kx-ωt)(5) 则波浪的速度场 u=Aωcoshk(z+h)sinhkhcos(kx-ωt)w=Aωsinhk(z+h)sinhkhsin(kx-ωt)(6) 其中:k为波数;ω为圆频率.在FLUENT软件中,动边界由墙体边界和动网格来实现.由冲程和波面方程关系求出冲程X0.同样,可以根据波浪周期、水深计算得到相关系数k、ω.推板在计算开始时做简谐运动从而生成波浪,可完成波浪在FLUENT软件中的实现. 1.3 水质点速度定义 定义水质点速度造波是指在水槽左端一定范围内,定义水质点的速度. 其目的是强迫造波区水质点的水平速度,使之达到模拟相应波浪的目的.定义水质点速度 U(c,t)=X0ωsin(kx-ωt)2.0(7) 2 动边界造波模型 二维数值波浪水槽模型总高0.75 m,长14 m,水深0.5 m,上部空气0.25 m.由Gambit划分为结构化网格,网格划分后单元格长0.08 m,高0.025 m.x方向网格数为175,z方向网格数为30,总网格数为5250.模型如图1(a)、(b)所示. 将动边界造波模型定义为模型1,定义水质点速度造波模型为模型2.模型1和模型2的边界条件设置为:模型1左边界AB取墙体边界条件并设置为动边界,模型2左边界AB取固壁边界条件,右边界CD均取为固壁边界条件,上边界AD均取为压力入口边界条件,底边界BC均取为固壁边界条件. 初始条件:流场中初始速度u=0,w=0,迭代控制精度ε=0.001,每步计算最多迭代30次. 利用FLUENT的前处理软件Gambit进行网格的划分和已知模型边界的定义.采用FLUENT中VOF模型、Segregated求解器、标准k-ω模型、压力速度耦合的PISO算法求解非定常状态下的紊流问题.应用DEFINESOURCE(mom-source,cell,thread,dS,eqn)宏来添加动量源项,实现水槽尾端消波. 3 数值波浪特征 在H0/h=0.1-0.3,h/λ=0.124-0.219,选择周期T0=1.3、1.5、1.7 s三种情况进行数值计算模拟.每种情况下选取3个波高,分别为H0=0.09、0.12、0.15 m.通过模拟计算得到的波浪周期和波长的统计结果分别如表1—2及图2—5所示. 从