基于vof模型的振荡水柱气室三维数值模拟

基于vof模型的振荡水柱气室三维数值模拟

作为一个长海的海洋资源大国，中国在利用海洋能源开发，尤其是海上波浪能源方面具有固有的优势。波浪能的分布广泛,能量密度高,波浪能发电多是在沿岸和海上进行,因此不占用土地,不需要迁移人口,且具有综合利用效益,多年来一直是海洋可再生能源中的开发热点。目前,振荡水柱波能发电装置是实用化与商业化开发程度最高的波浪能转换装置之一。该类装置主要由气室、输气管道、空气透平与发电机组构成。气室为一半淹没空仓结构,仓室前墙下部开敞,入射波浪将带动气室内的水柱做上下振荡,并带动水体上部的空气通过输气管道与外界大气之间做往复运动,称为能量一次转换过程。空气驱动透平电机转动,便可实现波浪能向电能的转换,称为能量二次转换过程。其中,气室是实现能量一次转换过程的主要结构。随着对波能发电装置重视程度的日益增加及研究手段的愈加丰富,国内外学者对气室和能量一次转换过程等相关问题进行了广泛而深入的研究。刘月琴等对传统岸式振荡水柱波能装置的气室部分进行了实验研究,重点考察了海岸岸坡之间的效应关系及其规律。当入射波周期较大时,海岸岸坡对振荡水柱的垂向运动响应影响不大,但对吸能效率和单位宽度波能俘获率的影响较为明显。梁贤光等在造波水槽中进行了汕尾100kW波力电站气室模型的性能试验。气室的波能转换效率不仅随波周期变化,而且随波高而异。试验还表明,气室的形状对波能转换率的影响也较大。Brito-Melo用直接计算方法预测振荡水柱装置的水动力特性,为求解气室内压力分布,他们改进了气室内自由水面处的动力学边界条件,同时利用压力辐射法计算气室内往复气流的压力。计算结果证明边界元方法能够有效地计算波浪的绕射和衍射问题。Lewis&Delaure提出一种计算岸式振荡水柱系统的数值模型,利用一阶混合面元法求解在定常边界条件下规则波与波能装置的相互作用。对规则波与不规则波的计算与实验结果的对比,验证了该模型的可靠性。本研究构建了基于计算流体力学软件Fluent的三维数值波浪水槽,用以研究气室前的波浪传播与气室内的水面振荡变化。该数值波浪水槽的控制方程为连续方程和基于雷诺平均的纳维尔斯托克斯方程,自由水面的重构方法为VOF法。通过与试验数据的对比,验证了本研究构建方法的准确性及在相关研究中应用的可能性。13d数值波厅1.1流体体积分数法三维数值波浪水槽造波板位于水槽的左侧,开边界则位于水槽的右侧(图1)。该水槽的数值波浪水槽分别为控制方程和基于雷诺平均的纳维尔斯托克斯(RANS)方程:∂ui∂xj=0(1)∂ui∂xj=0(1)∂ui∂t+ui∂ui∂xi+uj∂uj∂xj=−1ρ∂p∂xi+fxi+∂∂xj(v∂ui∂xj−u′iu′j¯¯¯¯¯¯¯¯)(2)∂ui∂t+ui∂ui∂xi+uj∂uj∂xj=-1ρ∂p∂xi+fxi+∂∂xj(v∂ui∂xj-u′iu′j¯)(2)式中,t为时间;xi,xj和ui,uj为x,z方向的坐标向量和对应的速度分量;ρ,p,ν和fxi为流体的密度、流体的压力、运动粘性常数和体积力。u′iu′j¯¯¯¯¯¯¯¯⋅ρu′iu′j¯⋅ρ定义为雷诺应力,该项也导致了上述方程式不闭合,需引入相应的湍流模型。本研究采用工程中应用较为广泛的标准k-ε湍流模型用于考察波浪运动中的湍动效应。自由水面的追踪与重构采用Hirt&Nichols提出的流体体积分数法,简称VOF方法。该方法假设在某一计算单元内水相所占体积与气相所占体积的总和为定值,水相的计算公式如下:∂aw∂t+∂(awui)∂xi=0(3)∂aw∂t+∂(awui)∂xi=0(3)式中,aw为水相体积分数;t为时间。此外,水气界面重构与通量计算通过由Youngs在1982年提出的假设单元内的线性斜坡的方法获得。1.2辐射边界条件为了更好的揭示振荡水柱波能发电装置的工作机理,主要采用规则入射波进行研究。造波方法为推板式,数值波浪水槽的开边界条件则采用Sommerfield辐射边界条件,该方法已经被证明对规则波的吸收较好。造波与消波边界的功能主要通过Fluent软件的用户自定义(UDF)函数实现,计算程序中的动网格运动则采用边界铺层的方法。计算网格则由Gambit2.0生成,并在水气交界面处进行了加密。控制方程的数值求解方法采用有限体积法,对流项的求解则为二阶迎风格式。速度压力耦合使用“含有分裂算子压力隐式的非迭代时间步进(NITA-PISO)”方法。上述方法与VOF法相兼容,且能够显著降低计算耗时。2气室结构参量振荡水柱波能发电装置前方的波浪传播、气室内的自由水面振荡及输气管内的往复气流的运动都受到气室形状的影响。因此,气室及其相关参量在能量一次转换过程中具有重要作用。开展针对气室形状参量的研究,对揭示振荡水柱波能装置的工作机理、提高气室工作性能将有较高的工程意义和实用价值。本研究所采用的气室主要的气室形状参量包括:气室宽度lf;气室前墙吃水深度ds;气室内静水面上高度lh;气室前墙厚度ls;气室顶部开口宽度ld,底坡水平长度lm;底坡坡角θs;气室墙前平均水深h。若ld=lf,则代表气室为完全开敞,否则为部分开敞。圆柱形输气管安装在气室的顶部中央位置,输气管的直径为ld,高度为hd(图2)。为了考察振荡水柱波能装置结构对气室工作性能的影响、验证所构建三维数值波浪水槽的准确性,研究计算了不同的气室结构参量(表1)。除表中的形状参量之外,输气管的高度hd=13.2m,底坡水平长度lm=22.86m,底坡坡角θs=26°,气室的墙前平均水深为16m。此外,本研究采用的入射波高度为1m,入射波周期为3.5~8.0s。3气室中相对波幅的变化情况工况1~3中,在不同入射波周期条件下气室内波面生成变化特征相似(表1),以工况2入射波周期T=5.5s时气室内的波面升沉变化为例。如图3a所示,当振荡水柱的波面达到波峰时,气室内各点处的波面高程并不相等:气室前墙及后墙附近的波面高程大于中部,但前墙与后墙附近的波面达到最高的时刻不同。前、后墙两位置处的波面高程平均值约等于气室中部的高程值。由图3b可见,气室内振荡水柱达到波谷时,其自由水面的运动特征与其到达波峰时的规律基本相同,即气室内自由水面高程分布沿入射波方向存在区别,而在垂直于入射波方向上波面高程的分布差别不大。相对波幅指气室内的振荡波幅(a)与入射波幅(a0)的比值,反映了气室内水面振荡的幅度和剧烈程度。当气室内的波幅较大时,振荡水柱将推动更多的空气在输气管内做往复运动,意味着更多的波浪能量可以在能量一次转换的过程中变为空气能量。因此,气室内的相对波幅可以较为准确的反应该装置工作性能的优劣。由图4可见工况2条件下,入射波周期为T=5.5s时气室中部自由水面相对波幅的变化情况。气室中自由水面的摆式振荡造成波幅升沉运动不再呈正弦变化。此外,自由水面振荡的周期与入射波周期基本相同,且入射波幅有明显的提高。各工况条件下的振荡水柱波能发电装置气室内中部的自由水面相对振幅三维数值模拟计算结果如图5所示。为更好的验证计算结果,本研究与ZhenLiou等二维数值模拟计算也进行了比较。从总体上看,三维数值模拟结果与试验数据较为吻合;在长周期区,三维数值计算的结果与试验数值更为符合,优于刘月琴中的二维数值模拟计算结果,但仍小于试验值,这是由于相对于二维计算,三维数值模拟能够更好的处理和垂直于入射波向的波面分布情况及其对气室中部的影响。从总体上看,本研究各形状参量条件下的相对波幅沿入射波周期的分布规律基本相同。在短周期区(T=3.5~6.0s),气室内的相对波幅随着周期的增加而增加。在长周期区(T=6.0~8.0s),气室内的相对波幅则随着周期的增加而趋于一个较为稳定的数值,相对波幅在2.0左右。由此可见,在进行波浪能装置运行的选址过程中,应尽量选择实海况平均周期大于6.0s的海域。在上述条件下,波能能量的转换效率较高,且装置的工作性能也较高。4振荡水柱波能发电装置气室水体振荡变化的数值模型本研究构建了基于计算流体力学软件Fluent的三维数值波浪水槽。该数值水槽的控制方程为连续方程和RANS方程,自由水面追踪与重建则采用PLIC-VOF方法,并将上述