Numerical simulation of impulse characteristics on oscillating buoy wave energy converter
-
摘要: 振荡浮子波能装置垂荡获能过程中会受到强烈的波浪冲击作用，因此有必要对浮子的波浪冲击特性开展研究。基于雷诺平均纳维−斯托克斯方程（RANS）建立了波浪与浮子耦合作用的数值模型，通过试验数据验证了所建模型的准确性，采用建立的模型分析了浮子表面的波浪冲击特性。通过数值分析得到以下结论：浮子底部受到的波浪冲击压力近似为正弦脉冲形式，在波浪行进方向上与前侧距离越大，无量纲冲击压力值越大；浮子迎浪侧波浪冲击压力随散射参数的增大而增加，其他位置处的冲击压力受散射参数的影响不明显；浮子上的最大波浪冲击压力受能量摄取（PTO）系统的影响很小。Abstract: The oscillating buoy wave energy converter will be impacted strongly during oscillating and energy capture, so it is necessary to study the wave impact characteristics of the buoy. The numerical model of the coupling effect between wave and buoy is built based on Reynolds-averaged Navier-Stokes equations (RANS). The accuracy of the model is verified by experimental data, and the wave impact characteristics on the surface of the buoy is deduced by the numerical simulation. The numerical analysis result shows that the wave impulse on the bottom of the buoy is approximately sinusoidal, and at the same time, the larger the distance from the front side in the wave traveling direction, the larger the dimensionless impact pressure value. The wave impact pressure on the wave-facing side of the buoy increases with the increase of scattering parameters, while the impact pressure at other positions is not obviously affected by scattering parameters. The maximum wave impact pressure on the buoy is little affected by the power take-off (PTO) system.
-
在“双碳”的战略目标下，海洋可再生能源迎来发展新格局[1]。波浪能因具有开发难度较低、获取方式简单等优势，被认为是品质最高的海洋能[2-3]。振荡浮子波能发电装置[4-5]因结构简单、能量转换效率高等优势是目前常见的三大波浪能利用技术之一。
振荡浮子装置大多是由波能俘获装置、质量体、系泊系统以及发电装置组成，其工作原理是利用俘获装置在波浪的作用下振荡起伏来达到汲取波浪能的目的，之后再由能量摄取（power take-off ，PTO）系统实现波浪能向电能转化。振荡浮子式波浪能发电装置的水动力性能[6]、能量转换效率[7]及性能优化[8]成为目前研究的热点。美国Ocean Power Technologies (OPT) 公司研发了PowerBuoy振荡浮子波能装置，并对其进行了设计优化[9]；Shi[10]和Han等[11]对振荡浮子获能特性进行了理论分析，并进一步提出了获能谱的概念；于通顺等[12-13]探究了振荡浮子前侧波浪爬升特性，为后续浮子竖向尺寸设计提供理论依据；Agamloh等[14]利用有限体积法分别模拟单个浮子和双浮子在不同规则波作用下的运动过程，探究了浮子的能量输出功率。吴必军等[15]和余海涛等[16]均对单浮子波能装置的俘获宽度比进行了试验研究，并指出波浪频率是影响装置俘获宽度比的主要因素。Waters等[17]对驱动线性发电机的浮子进行了试验研究，提出非线性外负载是影响装置俘获宽度比的主要因素。
通过以上研究发现，前人的研究中尚未开展浮子垂荡过程中表面所受的冲击作用的研究。振荡浮子在复杂海况中通过大幅垂荡获能，垂荡过程中浮子与波浪作用过程剧烈，忽略了浮子表面受到的强烈波浪冲击作用，将会增大浮子长期运行的风险。因此本文基于FLOW 3D建立波浪与浮子耦合作用的数值模型，开展浮子所受冲击作用过程的研究，探究浮子所受冲击荷载的作用规律。
1. 数值计算方法
1.1 控制方程
对于波浪和结构物相互作用的问题，波浪被认为是不可压缩黏性流体的运动。考虑湍流流动，将连续性方程和雷诺平均纳维−斯托克斯（Reynolds-averaged Navier-Stokes，RANS）方程作为流体运动的控制方程。假设右手笛卡尔坐标系O–xyz以未受扰动自由面的原点定义，则连续性方程可以写作：
$$ \dfrac{{\partial \left( {{\boldsymbol{u}}{A_x}} \right)}}{{\partial x}} + \dfrac{{\partial \left( {{\boldsymbol{v}}{A_y}} \right)}}{{\partial y}} + \dfrac{{\partial \left( {{\boldsymbol{w}}{A_z}} \right)}}{{\partial z}} = 0 $$ 动量方程为
$$ \begin{gathered} \dfrac{{\partial {\boldsymbol{u}}}}{{\partial t}} + \dfrac{1}{{{V_{\text{F}}}}}\left( {u{A_x}\dfrac{{\partial {\boldsymbol{u}}}}{{\partial x}} + v{A_y}\dfrac{{\partial {\boldsymbol{u}}}}{{\partial y}} + w{A_z}\dfrac{{\partial {\boldsymbol{u}}}}{{\partial z}}} \right) = - \dfrac{1}{\rho }\dfrac{{\partial {\boldsymbol{P}}}}{{\partial x}} + {{\boldsymbol{G}}_x} + {{\boldsymbol{f}}_x} \\ \end{gathered} $$ $$ \begin{gathered} \dfrac{{\partial {\boldsymbol{v}}}}{{\partial t}} + \dfrac{1}{{{V_{\text{F}}}}}\left( {u{A_x}\dfrac{{\partial {\boldsymbol{v}}}}{{\partial x}} + v{A_y}\dfrac{{\partial {\boldsymbol{v}}}}{{\partial y}} + w{A_z}\dfrac{{\partial {\boldsymbol{v}}}}{{\partial z}}} \right) = - \dfrac{1}{\rho }\dfrac{{\partial {\boldsymbol{P}}}}{{\partial y}} + {{\boldsymbol{G}}_y} + {{\boldsymbol{f}}_y} \\ \end{gathered} $$ $$ \begin{gathered} \dfrac{{\partial {\boldsymbol{w}}}}{{\partial t}} + \dfrac{1}{{{V_{\text{F}}}}}\left( {u{A_x}\dfrac{{\partial {\boldsymbol{w}}}}{{\partial x}} + v{A_y}\dfrac{{\partial {\boldsymbol{w}}}}{{\partial y}} + w{A_z}\dfrac{{\partial {\boldsymbol{w}}}}{{\partial z}}} \right) = \\ \qquad\qquad\;\;\; - \frac{1}{\rho }\dfrac{{\partial {\boldsymbol{P}}}}{{\partial z}} + {{\boldsymbol{G}}_z} + {{\boldsymbol{f}}_z} \\ \end{gathered} $$ 式中：
$\rho $ 为流体密度，${\boldsymbol{P}}$ 为流体内部压强，${V_{\text{F}}}$ 为可流动的体积分数，${\boldsymbol{u}}$ 、${\boldsymbol{v}}$ 、${\boldsymbol{w}}$ 分别为在$x$ 、$y$ 、$z$ 这3个方向上的流速分量，${A_x}$ 、${A_y}$ 、${A_z}$ 分别为$x$ 、$y$ 、$z$ 这3个方向可流动的面积分数，${{\boldsymbol{G}}_x}$ 、${{\boldsymbol{G}}_y}$ 、${{\boldsymbol{G}}_z}$ 分别为对应$x$ 、$y$ 、$z$ 这3个方向的重力加速度，${{\boldsymbol{f}}_x}$ 、${{\boldsymbol{f}}_y}$ 、${{\boldsymbol{f}}_z}$ 分别对应3个方向的黏滞力加速度。1.2 紊流模型
本文采用重组化RNG k–ε模型，其表达式为
$$ \begin{gathered} \frac{{\partial {k_{\text{T}}}}}{{\partial t}} + \frac{1}{{{V_F}}}\left( {{\boldsymbol{u}}{A_x}\frac{{\partial {k_{\text{T}}}}}{{\partial x}} + {\boldsymbol{w}}{A_z}\frac{{\partial {k_{\text{T}}}}}{{\partial z}}} \right) = {P_{\text{T}}} + {K_{{\text{T,Diff}}}} - {\varepsilon _{\text{T}}} \\ \end{gathered} $$ $$ \begin{gathered} \frac{{\partial {\varepsilon _{\text{T}}}}}{{\partial t}} + \frac{1}{{{V_{\text{F}}}}}\left( {{\boldsymbol{u}}{A_X}\frac{{\partial {\varepsilon _{\text{T}}}}}{{\partial x}} + {\boldsymbol{w}}{A_Z}\frac{{\partial {\varepsilon _{\text{T}}}}}{{\partial z}}} \right) = \\ \frac{{{C_{{\text{DIS,}}1}} \cdot {\varepsilon _{\text{T}}}}}{{{k_{\text{T}}}}}{P_{\text{T}}} + {\varepsilon _{{\text{T,Diff}}}} - {C_{{\text{DIS,2}}}}\frac{{{\varepsilon _{\text{T}}}^2}}{{{k_{\text{T}}}}} \\ \end{gathered} $$ 式中：PT为由剪切效应引起的湍流动能的生成项，
$ {K_{{\text{T,Diff}}}} $ 和$ {\varepsilon _{{\text{T,Diff}}}} $ 为扩散项，CDIS,1和CDIS,2为模型系数，k为由流速梯度引起的紊动动能的产生项。1.3 运动与碰撞模型
Flow-3D中运动与碰撞模型可以模拟流固耦合下结构物的运动，该运动可以是指定的运动，也可以是与流体动态耦合的运动，移动物体都可以具有6个自由度。该模型还允许存在多个（最多100个）线性弹簧和扭力弹簧、弹性绳索以及系泊缆绳，它们可以附接到移动的物体上，并可对其施加力或扭矩。根据刚体运动学，刚体的一般运动可分为平移和旋转运动。结构物上任何点的速度等于结构物上任意选择的基点速度加上绕该基点旋转引起的速度。选择S为刚体上的基点，R为刚体上任意点，则R点速度表达式为
$$ {{\boldsymbol{V}}_R} = {V_S} + {\boldsymbol{\omega}} \times {{\boldsymbol{r}}_{R/S}} $$ 式中：
${{\boldsymbol{V}}_S}$ 为基点的速度；${{\boldsymbol{r}}_{R/S}}$ 为点R相对于基点S的矢量距离；${\boldsymbol{\omega}}$ 是刚体的角速度矢量，与基点的选择无关。1.4 自由表面追踪的流体体积方法
FLOW-3D中采用了流体体积方法进行自由追踪。定义流体计算域内独个单元的流体体积函数为
${F}（x,y,z,t)$ 。F为流体所占单元体的比值，若单元全部充满流体，则F值为1；若不存在流体时，F值则为0；而当F介于0~1时，则表明该网格内存在自由表面。F函数表示单位体积，并满足方程：$$ \frac{{\partial {F}}}{{\partial t}} + \frac{{\partial {\boldsymbol{u}}{F}}}{{\partial x}} + \frac{{\partial {\boldsymbol{v}}{F}}}{{\partial y}} + \frac{{\partial {\boldsymbol{w}}{F}}}{{\partial z}} = 0 $$ 2. 数值模型的建立
2.1 三维数值波浪水槽的建立
建立数值波浪水槽的侧面示意如图1（a）所示，水槽长度取6倍左右波长，宽度取浮子宽度的5倍。为便于后续进行数值模型的验证，选用Yu等[12]试验中浮子的尺寸进行模型的建立：圆柱型浮子直径为0.4 m，高度为0.4 m，吃水深度为0.12 m。为了合理描述浮子的形状，该水槽由3个网格块组成，中间网格块为浮子运动加密区，浮子放置在距离入口的2倍波长(L)范围内。水槽的边界条件定义为：最左端是可以实现函数造波功能的波浪，最右端为出流，水槽的底部设为壁面，上部为压力，水槽顺波向两侧对称。如图1（b）和图1（c）所示。水槽尾端部分（约1/4水槽长度）设有阻尼消波区用来抑制波浪反射，左端和右端的阻尼系数分别设为1和8。
图 1 三维数值波浪水槽边界设置及网格块划分
下载:
全尺寸图片
2.2 网格的划分及网格独立性验证
本模型采用波高范围内3种不同的网格数量进行了波面网格独立性试验，即设置了网格数量分别为10、20和30的3个不同算例，流域内网格的大小和数量如表1所示。
表 1 无浮子时计算流域网格数量和尺寸波高范围
网格数量网格尺寸/cm x方向总数量 网格总数量 运行时间/min 10 0.02 450 1125000 128.2 20 0.01 900 4100000 436.33 30 0.0067 1350 10125000 1222.3 图2比较了3种网格尺寸下计算区域内相同位置（2L处）自由面高程的历时曲线，结果发现，加密网格数为20与30的计算结果差距较小，加密网格数为10的计算结果与其他2种网格数量的计算结果差距较大，但是加密网格数为30的运行时间明显增加。综合考虑模型计算时间和精度，本文在Z方向的波高(H)范围内网格尺寸建议设置为
$ \Delta z \leqslant {H \mathord{\left/ {\vphantom {H {20}}} \right. } {20}} $ 。
图 2 不同网格尺寸的自由面高程历时曲线下载:
全尺寸图片
为了合理描述浮子的形状，在浮子运动加密区采用3种不同总网格数量，浮子周围网格数量如表2所示，表中∆x、∆y和∆z分别为沿x、y和z方向的网格尺寸。所划分网格可以较好地描述浮子的形状。
表 2 浮子周围的网格数量和尺寸网格信息 网格总数量 ∆x/cm ∆y/cm ∆z/cm 1 676800 0.02 0.02 0.02 2 1664000 0.01 0.01 0.01 3 2995200 0.005 0.01 0.005 图3比较了在不同网格尺寸下采用上述波浪条件捕捉浮子垂荡运动的能力。结果表明，0.01 cm和0.005 cm的网格尺寸下浮子位移相差不大，说明以0.01 cm划分的网格具有足够高的模拟分辨率。因此，浮子周围网格尺寸建议设置为
$\Delta x < L/100$ 和$ \Delta z \leqslant {H \mathord{\left/ {\vphantom {H {20}}} \right. } {20}} $ 。
图 3 不同网格尺寸的浮子垂荡位移历时曲线下载:
全尺寸图片
2.3 三维数值水槽验证
2.3.1 造波能力验证
图4为距波浪入口2L附近（放置浮子处）规则波的数值计算结果与理论值的时程曲线（波高H=0.2 m，周期T=2 s，水深h=1 m）。在入射波传递一段时间到达稳定状态后，很明显波峰处略尖陡，波谷处略平坦，且相对误差低于4%；同时可以发现波浪的反射效应并不明显，这表明水槽的造波能力和消波效果都比较可靠。
图 4 自由波面理论值与计算值的时程曲线下载:
全尺寸图片
2.3.2 波浪对浮子作用的验证
为验证数值模型对浮子运动捕捉的准确性，选用Yu等[12]的试验数据与数值计算结果进行比较。表3给出了连续5个振动周期中浮子运动最高点与最低点的数值结果与试验数据的对比（波高H=0.2 m，周期T=2 s，水深h=1.0 m）。可以看出，数值结果与试验数据的相对误差均在10%以内，所建立数值模型能够较好地模拟波浪对浮子的作用。
表 3 垂荡位移数值结果与试验数据数值结果/m 试验数据/m 相对误差/% 最高点坐标 最低点坐标 振幅 最高点坐标 最低点坐标 振幅 0.1069 −0.1105 0.1087 0.1136 −0.1129 0.1132 3.98 0.1014 −0.1064 0.1039 0.1046 −0.1162 0.1104 5.89 0.0999 −0.1040 0.1019 0.1085 −0.1155 0.1120 9.02 0.1015 −0.1066 0.1041 0.1028 −0.1156 0.1092 4.67 0.0912 −0.1083 0.0997 0.1009 −0.1183 0.1096 9.03 2.3.3 波浪冲击压力观测点设置和定义
根据对称性，在浮子迎浪方向、侧面以及后方（分别为Row A、Row B和Row C）设置3列压力传感器，定义s为浪高仪至浮子底端的垂直距离，d为浮子高度，浮子柱面观测点和压力传感器的布置如图5（a）和表4所示。浮子底部设置了11个压力传感器，θ为所在点与x轴负方向的夹角，相邻2组测点之间的夹角为45°，浮子底面观测点和压力传感器的布置如图5（b）和表5所示。模型中所有虚拟观测点的布置与实际位置误差小于3 mm。
图 5 波浪冲击压力虚拟观测点设置下载:
全尺寸图片
表 4 浮子柱面压力传感器设置s/d Row A Row B Row C 0 A1 B1 C1 0.15 A2 B2 C2 0.30 A3 B3 C3 0.45 A4 B4 C4 0.60 A5 B5 C5 波浪在冲击的过程中涉及到诸多因素，比如水气掺混、瞬时效应等，这些因素导致波浪冲击垂荡浮子的过程具有很强的随机性，因此冲击过程中的压力峰值具有很强的离散性。本文根据文献[18]对每个虚拟的压力传感器采集的压力数据进行了处理，即浮子上的压力测点取样时，将每个周期内的压力峰值由大到小排列，并且取前1/3个峰值的平均值作为各测点的冲击压力峰值特征值，即统计分析的特征值。每个测点的冲击压力峰值的实测数据记为Pt，浮子静止时测点的初始压力记为P0，将实测数据与初始压力无量纲化后记为Pt/P0。浮子所受的水平和竖直方向波浪荷载峰值数据的提取则是将每个周期内的压力峰值提取出来，并且对峰值相差较小的几组数据取平均值为波浪荷载的特征值。水平和垂直波浪荷载分别用Fx和Fz来表示，对波浪荷载无量纲化分析，选取浮子的重力ρgSd为无量纲化参数，其中ρ为浮子的密度，g为重力加速度，S为浮子底面积。
表 5 浮子底面压力传感器设置θ/（°） 测点 0 O~M1~M2 45 O~N1~N2 90 O~P1~P2 135 O~Q1~Q2 180 O~S1~S2 2.3.4 计算工况
数值计算中选用圆柱形浮子，实际浮子高度为1.6 m、质量为3860 kg、吃水深度为0.48 m。计算工况如表6所示。
表 6 计算工况设置工况 波高/m 周期/s 水深/m 半径/m PTO阻尼 阻尼/（N·s·m−1） 刚度/（N·m−1） 1 0.4 4 4 1.6 0 0 2 0.6 4 4 1.6 0 0 3 0.8 4 4 1.6 0 0 4 0.8 4 4 1.2 0 0 5 0.8 4 4 2.0 0 0 6 0.8 4 4 1.6 0 0 7 0.8 4 4 1.6 0 0 8 0.8 4 4 1.6 0 2000 9 0.8 4 4 1.6 0 4000 10 0.8 4 4 1.6 0 6000 11 0.8 4 4 1.6 2000 0 12 0.8 4 4 1.6 4000 0 13 0.8 4 4 1.6 6000 0 14 0.6 6 4 1.6 0 0 3. 波浪对垂荡浮子装置冲击作用
3.1 波浪冲击压力特征
3.1.1 冲击压力峰值在浮子的立侧处分布规律
图6为工况1中A3传感器记录的浮子波浪压力历时曲线，由图可以发现，浮子受到的冲击压力具有明显的周期性，周期与入射波周期保持一致，冲击压力近似为正弦脉冲形式，冲击阶段压力逐渐增大，之后随着波面下降开始逐渐回落，当测点脱离水体，压力维持在零值几乎不变。
图 6 浮子立侧不同压力测点冲击压力历时曲线下载:
全尺寸图片
图7为工况1条件下浮子表面不同测点上的压力分布，1#和2#测点冲击压力在不同位置处差距较小，3#测点的冲击压力差别较大。因为波面高度达不到，故4#和5#测点的冲击压力为0。距离浮子底部越大，无量纲冲击压力值越大，即相较于初始压力提升越大，3#测点的无量纲冲击压力达到最大值。距离浮子底部相同距离的测点中，Row C侧的冲击压力最大，Row A侧的冲击压力最小，表明浮子后侧所受到的冲击压力最大。
图 7 浮子立侧处不同测点压力分布下载:
全尺寸图片
3.1.2 冲击压力峰值在浮子底面分布规律
图8为工况1条件下M1#、O#和S1#测点记录的浮子的底部波浪压力历时曲线，由图可以发现，浮子底端受到的冲击压力具有明显的周期性，周期与入射波周期保持一致；浮子底部中心所受的波浪压力小于波浪传播方向上中心两侧的波浪压力。
图 8 浮子底面不同压力测点冲击压力历时曲线下载:
全尺寸图片
3.2 波浪参数和阻尼参数对浮子上波浪冲击压力垂直分布的影响
3.2.1 波陡参数对浮子周围波浪冲击效应影响
图9给出了不同波陡参数下，最大波浪冲击压力沿的垂直分布。
图 9 不同波陡参数下浮子上最大压力垂直分布下载:
全尺寸图片
从图9中可以发现，最大波浪冲击压力随着波陡参数的增加而增大，距离浮子底部距离越近的位置上波浪冲击压力越大。波陡较大时，浮子靠上位置处的波浪冲击压力最小值出现在Row B侧，这是由于浮子侧面的波浪爬升高度处于最低导致，其他位置的最大冲击压力出现在浮子中间（Row C），最小冲击压力依然发生在迎浪侧（Row A）。
3.2.2 散射参数对浮子周围波浪冲击效应影响
图10给出了不同散射参数下，最大波浪冲击压力沿Row A、Row B和Row C的垂直分布。从图10可以看出，改变浮子半径对于Row B和Row C侧相同压力测点上波浪冲击压力峰值的改变较小。散射参数较小时，Row A侧的相同测点压力峰值随散射参数变化较小；散射参数较大时，Row A侧的相同测点压力峰值有明显提升。总的来说，增加浮子半径对于浮子截面上相同位置的波浪冲击压力影响较小。
图 10 不同散射参数下浮子上最大压力垂直分布下载:
全尺寸图片
3.2.3 线性PTO阻尼参数对浮子周围波浪冲击效应影响
图11给出了不同弹簧刚度系数下，最大波浪冲击压力沿Row A、Row B和Row C的垂直分布。
图 11 不同弹簧刚度系数下浮子上最大压力垂直分布下载:
全尺寸图片
从图11中可以看出，弹簧刚度系数的增加对于浮子上的波浪冲击压力影响较小，浮子上相同压力测点上的压力峰值随弹簧系数的增加几乎保持不变。通过对比发现，施加PTO后对于浮子上的波浪冲击压力改变很小，且浮子上的波浪冲击压力垂直分布规律与浮子自由垂荡运动时保持一致。
图12给出了不同阻尼系数下，最大波浪冲击压力沿Row A、Row B和Row C的垂直分布。从图12中可以看出，增加阻尼系数对浮子截面上波浪冲击压力垂直分布的影响较小；增加PTO阻尼对浮子垂荡运动有影响，但浮子上的波浪冲击由于波浪参数保持一致，导致浮子上的波浪冲击压力受其他因素的影响很小。
图 12 不同阻尼系数下浮子上最大压力垂直分布下载:
全尺寸图片
浮子侧面的波浪冲击压力受波陡参数的影响最大，冲击压力随波陡增大而增大；Row A侧的波浪冲击压力随散射参数增大而增加，其他位置处的冲击压力受散射参数的影响不明显；浮子上的最大波浪冲击压力受PTO系统阻尼的影响很小；浮子侧面的压力分布在s>0.24 m位置上的压力处于最小值，这主要是由于波浪在该位置的爬升较小，冲击强度较低所导致。
3.3 最大波浪荷载与波浪冲击压力关系研究
图13和图14给出了工况3水平和垂直波浪荷载与Row A压力测点时程对照。A1和A2处的波浪压力与水平波浪荷载的历时曲线存在较小的相位差，水平波浪荷载的最大值与最小值出现的时间较早于A1处的波浪压力最值；浮子迎浪侧距离底部较高位置处（A3、A4和A5）的波浪压力与水平波浪荷载的历时曲线几乎不存在相位差，即水平波浪荷载的最大值与最小值出现的时间和该位置处的波浪冲击压力最值同步出现。
图 13 水平波浪荷载与Row A压力测点时程曲线对照下载:
全尺寸图片
图 14 垂直波浪荷载与Row A压力测点时程曲线对照下载:
全尺寸图片
由于水平波浪荷载受力位置与Row A侧相似，都为迎浪面，因此历时曲线存在的相位差较小，在浮子较高位置处的相位差几乎一致。但浮子所受垂直波浪荷载的受力面与浮子Row A侧波浪压力受力面存在较大的差异，从图14中可以发现，垂直波浪荷载的最值出现时间和该位置处的波浪冲击压力最值出现时间存在较大差异。A1和A2处的波浪压力与垂直波浪荷载最值出现的时间相反，即垂直波浪荷载出现最大值时，该位置的波浪压力处于最小位置；浮子迎浪侧距离底部较高位置处（A3、A4和A5）的波浪压力与垂直波浪荷载的最值同步出现，与水平荷载相似。
4. 结论
1）在规则波作用下，浮子所受冲击压力具有明显的周期性，浮子底部波浪冲击压力近似为正弦脉冲形式。
2）在入射波前进方向上，随着浮子底部传感器与浮子前部距离的增大，传感器所受的冲击压力峰值也增大。
3）波浪冲击压力受波陡参数的影响最大，冲击压力随波陡增大而增大。浮子迎浪侧的波浪冲击压力随散射参数的增大有增加的趋势，其他位置处的冲击压力受散射参数的影响不明显。浮子上的最大波浪冲击压力受PTO阻尼系统的影响很小。
4）浮子迎浪侧距离浮子底部较高位置处的波浪冲击压力与波浪荷载最值出现时间几乎保持一致。
5）本研究只针对线性波浪作用下的浮子冲击压力进行了研究，随机波浪作用下浮子荷载特征的研究需进一步开展。
-
图 1 三维数值波浪水槽边界设置及网格块划分
下载:
全尺寸图片
图 2 不同网格尺寸的自由面高程历时曲线
下载:
全尺寸图片
图 3 不同网格尺寸的浮子垂荡位移历时曲线
下载:
全尺寸图片
图 4 自由波面理论值与计算值的时程曲线
下载:
全尺寸图片
图 5 波浪冲击压力虚拟观测点设置
下载:
全尺寸图片
图 6 浮子立侧不同压力测点冲击压力历时曲线
下载:
全尺寸图片
图 7 浮子立侧处不同测点压力分布
下载:
全尺寸图片
图 8 浮子底面不同压力测点冲击压力历时曲线
下载:
全尺寸图片
图 9 不同波陡参数下浮子上最大压力垂直分布
下载:
全尺寸图片
图 10 不同散射参数下浮子上最大压力垂直分布
下载:
全尺寸图片
图 11 不同弹簧刚度系数下浮子上最大压力垂直分布
下载:
全尺寸图片
图 12 不同阻尼系数下浮子上最大压力垂直分布
下载:
全尺寸图片
图 13 水平波浪荷载与Row A压力测点时程曲线对照
下载:
全尺寸图片
图 14 垂直波浪荷载与Row A压力测点时程曲线对照
下载:
全尺寸图片
表 1 无浮子时计算流域网格数量和尺寸
波高范围
网格数量网格尺寸/cm x方向总数量 网格总数量 运行时间/min 10 0.02 450 1125000 128.2 20 0.01 900 4100000 436.33 30 0.0067 1350 10125000 1222.3 表 2 浮子周围的网格数量和尺寸
网格信息 网格总数量 ∆x/cm ∆y/cm ∆z/cm 1 676800 0.02 0.02 0.02 2 1664000 0.01 0.01 0.01 3 2995200 0.005 0.01 0.005 表 3 垂荡位移数值结果与试验数据
数值结果/m 试验数据/m 相对误差/% 最高点坐标 最低点坐标 振幅 最高点坐标 最低点坐标 振幅 0.1069 −0.1105 0.1087 0.1136 −0.1129 0.1132 3.98 0.1014 −0.1064 0.1039 0.1046 −0.1162 0.1104 5.89 0.0999 −0.1040 0.1019 0.1085 −0.1155 0.1120 9.02 0.1015 −0.1066 0.1041 0.1028 −0.1156 0.1092 4.67 0.0912 −0.1083 0.0997 0.1009 −0.1183 0.1096 9.03 表 4 浮子柱面压力传感器设置
s/d Row A Row B Row C 0 A1 B1 C1 0.15 A2 B2 C2 0.30 A3 B3 C3 0.45 A4 B4 C4 0.60 A5 B5 C5 表 5 浮子底面压力传感器设置
θ/（°） 测点 0 O~M1~M2 45 O~N1~N2 90 O~P1~P2 135 O~Q1~Q2 180 O~S1~S2 表 6 计算工况设置
工况 波高/m 周期/s 水深/m 半径/m PTO阻尼 阻尼/（N·s·m−1） 刚度/（N·m−1） 1 0.4 4 4 1.6 0 0 2 0.6 4 4 1.6 0 0 3 0.8 4 4 1.6 0 0 4 0.8 4 4 1.2 0 0 5 0.8 4 4 2.0 0 0 6 0.8 4 4 1.6 0 0 7 0.8 4 4 1.6 0 0 8 0.8 4 4 1.6 0 2000 9 0.8 4 4 1.6 0 4000 10 0.8 4 4 1.6 0 6000 11 0.8 4 4 1.6 2000 0 12 0.8 4 4 1.6 4000 0 13 0.8 4 4 1.6 6000 0 14 0.6 6 4 1.6 0 0 -
[1] 魏文栋. 能源革命: 实现碳达峰和碳中和的必由之路[J]. 探索与争鸣, 2021(9): 23−25. [2] 肖惠民, 于波, 蔡维由. 世界海洋波浪能发电技术的发展现状与前景[J]. 水电与新能源, 2011(1): 67−69. [3] 李永国, 汪振, 王世明, 等. 国外波浪能开发利用技术进展[J]. 工程研究——跨学科视野中的工程, 2014, 6(4): 371−382. [4] FALCÃO A F D O. Wave energy utilization: a review of the technologies[J]. Renewable and sustainable energy reviews, 2010, 14(3): 899−918. doi: 10.1016/j.rser.2009.11.003 [5] MAO Laifeng, ZHANG Jun, QU Quanyou, et al. Investigation and analysis on the performance of the oscillating buoy wave power device under the South China Sea[J]. Applied mechanics and materials, 2013, 321/322/323/324: 1365-1372. [6] 史宏达, 曲娜, 曹飞飞, 等. 振荡浮子波能发电装置浮子运动性能的试验研究[J]. 中国海洋大学学报(自然科学版), 2017, 47(6): 124−130,145. [7] ZHANG Dahai, GEORGE A, WANG Yifei, et al. Wave tank experiments on the power capture of a multi-axis wave energy converter[J]. Journal of marine science and technology, 2015, 20(3): 520−529. doi: 10.1007/s00773-015-0306-5 [8] HALS J, FALNES J, MOAN T. Constrained optimal control of a heaving buoy wave-energy converter[J]. Journal of offshore mechanics and Arctic engineering, 2011, 133(1): 011401. doi: 10.1115/1.4001431 [9] VAN RIJ J, YU Y H, EDWARDS K, et al. Ocean power technology design optimization[J]. International journal of marine energy, 2017, 20: 97−108. doi: 10.1016/j.ijome.2017.07.010 [10] SHI Hongda, CAO Feifei, LIU Zhen, et al. Theoretical study on the power take-off estimation of heaving buoy wave energy converter[J]. Renewable energy, 2016, 86: 441−448. doi: 10.1016/j.renene.2015.08.027 [11] HAN Zhi, CAO Feifei, HANN M, et al. Study on the capture spectrum of wave energy conversion[J]. Applied Ocean research, 2021, 111: 102654. doi: 10.1016/j.apor.2021.102654 [12] YU Tongshun, TANG Yuying, SHI Hongda, et al. Numerical modelling of wave Run-up heights and loads on heaving buoy wave energy converter under the influence of regular waves[J]. Ocean engineering, 2021, 225: 108670. doi: 10.1016/j.oceaneng.2021.108670 [13] 于通顺, 唐渔滢, 黄淑亭. 振荡浮子式波能装置波浪爬升特性试验[J]. 科技导报, 2021, 39(6): 42−46. [14] AGAMLOH E B, WALLACE A K, VON JOUANNE A. Application of fluid-structure interaction simulation of an ocean wave energy extraction device[J]. Renewable energy, 2008, 33(4): 748−757. doi: 10.1016/j.renene.2007.04.010 [15] 吴必军, 刁向红, 王坤林, 等. 10kW漂浮点吸收直线发电波力装置[J]. 海洋技术, 2012, 31(3): 68−73. [16] 余海涛, 刘春元, 费腾, 等. 直驱式波浪发电系统研究[J]. 海洋开发与管理, 2013, 30(S1): 43−46. [17] WATERS R, STÅLBERG M, DANIELSSON O, et al. Experimental results from sea trials of an offshore wave energy system[J]. Applied physics letters, 2007, 90(3): 034105. doi: 10.1063/1.2432168 [18] 丁兆强. 波浪对透空式三维结构物的冲击作用研究[D]. 大连: 大连理工大学, 2009.
