| Issue |
JNWPU
Volume 41, Number 1, February 2023
|
|
|---|---|---|
| Page(s) | 188 - 197 | |
| DOI | https://doi.org/10.1051/jnwpu/20234110188 | |
| Published online | 02 June 2023 | |
Hydrodynamic performance analysis of tower wake vortex tidal current energy capture system with swing oscillating hydrofoil
摆式振荡翼塔筒尾涡潮流能捕获系统的水动力性能分析
1
College of Electrical, Energy and Power Engineering, Yangzhou University, Yangzhou 225100, China
2
Smart Energy Internet Research Institute, Yangzhou University, Yangzhou 225127, China
3
School of Marine Science and Technology, Northwestern Polytechnical University, Xi'an 710072, China
Received:
8
April
2022
In order to effectively utilize the tidal current energy in offshore wind farms, a swing oscillating hydrofoil system for tidal current energy capture is proposed. The system is constructed with oscillating hydrofoils as tidal current energy capture device combined with offshore wind turbine underwater tower. Lattice Boltzmann method and large eddy simulation method were combined, and the hydrodynamic characteristics and energy collection efficiency of a NACA0015 airfoil with different motion parameters were numerically analyzed with or without the influence of tower wake. The influence of pitch amplitude, heave amplitude and chord length of the hydrofoil on the energy collection efficiency of the system is analyzed from the aspect of vortex structure in the moving process. The results show that the energy collection efficiency of the oscillating wing is significantly improved by the trailing vortex of the tower, and the vortex shedding behind the tower can reduced by the addition of the oscillating wing, consequently, the fatigue load caused by the flow around the cylindrical body can be reduced. The energy capture efficiency of the oscillating wing can reach a highest efficiently of 30.72% in the wake vortex region of the tower, which has obvious advantages compared with the traditional oscillating wing tidal current energy capture device, and provides an effective way to improve the overall energy utilization rate of offshore wind farms.
摘要
为实现对海上风电场中潮流能有效利用, 提出一种将振荡水翼潮流能捕获装置与海上风力机水下塔筒结合的摆式振荡翼塔筒尾涡潮流能捕获系统。联合格子玻尔兹曼和大涡模拟方法, 对NACA0015翼型在有无塔筒尾流效应时, 不同运动参数对其水动力特性及能量收集效率的影响进行数值分析, 并从水翼运动过程中的旋涡结构分析了俯仰振幅、升沉振幅及振荡翼弦长对系统捕能效率的影响。结果表明, 塔筒尾涡显著提升振荡翼的能量收集效率, 同时振荡翼的加入能减慢塔筒涡脱, 降低由钝体绕流引起的疲劳载荷。摆式振荡翼在塔筒尾涡区域内的捕能效率最高可达30.72%, 与传统振荡翼潮流能捕获装置相比具有明显优势, 为提高海上风场整体能源利用率提供了一种有效的方式。
Key words: tidal current energy / oscillating hydrofoil / lattice Boltzmann method / large eddy simulation
关键字 : 潮流能 / 振荡水翼 / 格子玻尔兹曼方法 / 大涡模拟
© 2023 Journal of Northwestern Polytechnical University. All rights reserved.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
可再生能源的开发利用已然成为全球共识。与其他可再生能源相比，潮流能具有受天气影响小、能量密度高、可预测性强、载荷稳定等优势。潮流能发电装置可根据其原理分为旋转轴式和振荡翼式。与传统旋转轴式水轮机组相比，振荡翼式潮流能收集装置具有以下优点：振荡速度低，对水生生物伤害小，对环境更为友好；结构简单，可靠性高，能量获取效率高。
McKinney[1]最早提出振荡翼可用于能量收集，用正弦函数定义了振荡翼的俯仰和升沉，从理论和实验两方面都证明了振荡翼具有从来流中捕获能量的能力。随后，越来越多的学者开展了对振荡翼式能量收集装置的研究。对于升沉-俯仰耦合运动的传统振荡翼，Kinsey等[2-3]系统地研究了振荡翼运动参数对能量收集效率的影响，发现运动参数对能量收集效率有重要影响。除运动参数外，振荡翼运动过程中的流场结构也对装置的性能有一定影响。Young和Su等[4-5]都指出水翼运动过程中前缘涡的发展与脱落是振荡翼能否达到最佳功率的关键因素。Baik等[6]研究不同雷诺数下振荡翼运动过程中非定常力，发现前缘涡的形成与脱落取决于振荡翼的运动频率和俯仰振幅。Ribeiro等[7]用大涡模拟方法研究了不同雷诺数下振荡翼的能量收集效率，发现前缘涡的生成与脱落主要受振荡翼运动方式影响，与雷诺数的相关度较低。对非均匀流内的振荡水翼，谢玉东等[8]通过数值研究，发现非均匀流可以改善水翼的能量提取性能。
对于多个振荡翼排列方式对整体性能的影响研究，如2个振荡翼前后布置的工况，Akhtar等[9]发现在适当的运动参数下上游振荡翼的尾涡能提高下游振荡翼的推力，使整体的获能效率高于单水翼。对与前后串联布置的振荡翼，Lindsey和Jones等[10-11]进行了实验和数值研究，发现大振幅低频率有利于实现最大能量收集，较大的有效攻角和薄翼型能够实现更高的能量收集效率。Ashraf等[12]对前后布置的双列NACA0014翼型研究发现，尽管双列布置下每个翼型的能量采集量与单翼型相比降低，但总能量收集效率有显著提高。
上述研究都是针对传统简谐运动模型的振荡翼，而乔凯等[13]对水翼的俯仰运动进行了改进，改进后的运动模型能使水翼的升力系数在较长时间段内保持较大值，从而提高水翼的捕能效率。孙光等[14]提出一种带尾缘襟翼的振荡水翼结构，结果表明尾缘襟翼的摆动能增大水翼的攻角，进而提高时均功率系数。Karbasian等[15]首先提出摇臂式振荡翼能量采集运动模式，将水翼的平移与升沉、俯仰耦合，可以在较低的振荡频率下输出较高的功率。进一步地，马鹏磊等[16-17]从摆式振荡翼运动过程中攻角的变化和旋涡结构两方面分析了振荡频率和俯仰振幅对水翼水动力性能的影响。
潮流能和海上风能是可再生海洋能源的重要部分，其中海上风能的开发利用较为成熟，但海上风场中的潮流能却无法有效利用。因此，将摆式运动的振荡翼式潮流能收集装置与海上风电场结合，不仅对提高海上风场效能具有重要意义，还对于消除潜在的塔筒涡激振动具有积极意义。本文提出一种摆式振荡翼塔筒尾涡能量收集装置，运用格子玻尔兹曼方法(lattice Boltzmann method, LBM)和大涡模拟方法(large eddy simulation, LES)分析了不同运动参数下系统的水动力性能及能量收集效率，提出了一种高效收集海上风场潮流能的参考方案。
1 运动模型
图1为摆式振荡翼运动示意图，水翼的俯仰轴绕塔筒圆心做周期性的往复摆动，同时水翼绕其俯仰轴做俯仰运动，其数学模型为：
式中, θ0为振荡翼的俯仰振幅; f为振荡频率; x(t)为水翼俯仰轴x方向的位移; y(t)为水翼俯仰轴y方向的位移; R为摆动半径; γ为摆动角幅值; H0为升沉振幅; 取俯仰轴位置为翼型弦长1/3处, 即XP/c=1/3。
在来流作用下振荡翼产生的水动力可分解为垂直于来流方向的力FY(t)和平行于来流方向的力FX(t)及因俯仰运动而产生的力矩M(t)。振荡翼的瞬时功率系数CP(t)及时均功率系数CP可表示为
式中：ρ为海水密度; U∞为来流速度; c为振荡翼弦长; CPX(t)为分力FX(t)贡献的瞬时功率系数; CPM(t)为俯仰力矩贡献的瞬时功率系数; CPX为分力FX(t)贡献的时均功率系数; CPM为俯仰力矩贡献的时均功率系数; T为振荡翼运动周期; P为时均功率。
攻角是影响水翼性能的重要因素之一, 受振荡翼摆动影响, 等效攻角为
式中, αh(t)为升沉诱导角。
对二维情况而言, 振荡翼的扫掠面积S与其扫掠高度相当, 则振荡翼从来流中捕获能量的效率为
式中：d为振荡翼厚度; H0为扫掠高度。无量纲频率定义为
 |
图1 摆式振荡翼运动模型示意图 |
2 数值方法
2.1 LBM方法
LBM方法是对Boltzmann-BGK方程[18]进行离散, 进而把连续介质处理成位于网格节点的众多离散流体质点粒子, 其中粒子按碰撞及迁移规则在网格上运动, 最后统计各网格流体质点及运动特征得到流体的宏观运动规律。在Boltzmann方程中, 分布函数被用来描述粒子的统计特征
式中：x为粒子的位置坐标；ξ为粒子的速度；a为粒子的加速度；Ω为碰撞项。对碰撞项进行BGK近似(bhatnagar-gross-krook, BGK), 简化碰撞项得到动力学演化方程[19]
式中: fi(x, t)为t时刻，在位置x处, i方向粒子的分布函数; ci为粒子在i方向上的运动速度; fieq是平衡分布函数; τ为无量纲松弛时间; Δt为时间步长, (13)式等号左侧为迁移项, 右侧为碰撞项。
对于离散速度模型, 本文采用单松弛格式的D2Q9模型(2个空间维度, 9个方向的离散速度), 各方向粒子的速度分布如图2所示。其离散速度为
局部平衡态分布函数可表示为
式中：ρ为宏观密度；u为宏观速度；c为格子声速, c=
, RT=1/3。ωi为权系数, 对于D2Q9模型, 有ωi=0=4/9, ωi=1, 2, 3, 4=1/9, ωi=5, 6, 7, 8=1/36。
 |
图2 D2Q9速度离散模型 |
2.2 LES方法
LES方法借助滤波器方程完成对“小涡”的滤波处理, 从而形成关于大涡的控制方程, 小尺度的涡则通过亚格子模型来描述。常用的亚格子湍流模型有Smagorinsky模型、Dynamic Smagorinsky模型以及WALE(wall adaping local eddy)模型。WALE模型[20]通过优化应变率张量, 优化了亚格子黏性的计算, 使得对近壁面湍流及转捩的计算结果更准确。因此, 本文选取WALE模型作为亚格子模型。其涡黏度为
式中: Ls=min(κd, Δ1/3); κ为冯卡门系数; d为距壁面最短距离, 过滤尺度为Δ=V1/3; V为计算网格单元的体积。应变率张量为
LBM方法与LES方法的耦合是通过WALE亚格子模型改变LBM方法中碰撞模型的松弛时间来实现的, 简言之滤波尺度之下的小涡对流场的影响是通过加入湍流松弛时间实现的, 具体的方程在此不再赘述, 请参考文献[21], 此文中首次将大涡模拟引入到LBM中。
2.3 可靠性验证
为验证本文数值方法的可靠性, 以NACA0015翼型为研究对象, 并将本文方法与文献[3]在相同工况下的结果进行对比。计算域大小为75c×75c, c为翼型弦长, 格子分布如图3所示, 近壁面区域y+=0.27, 尾流采用动态自适应跟踪优化[22]。时间步长取T/2 000, 根据前期研究及文献[23-24], 确定一个时间步长机构走过振荡翼运动周期的1/2 000。水翼与塔筒的几何边界均为壁面条件, 关于LBM对曲线边界的捕捉方法请见文献[18]。LES采用的WALE亚格子模型常数设置Cw=0.325。
图4a)为一个运动周期内升力系数变化, t/T=0.25时，水翼表面压力系数的对比如图4b)所示, 结果表明本文的计算结果与文献[3]具有较好的一致性。数值计算工况为：NACA0015，Re=1 100，Xp/c=1/3，y0/c=1.0，f*=0.14，θ0=76.33°。数值计算结果对比如表1所示, 其中CX为平均阻力系数, 
为最大升力系数, 
3为最大俯仰力矩系数, 各参数误差都在1%以内, 证明了本文数值计算方法的可靠性。
 |
图3 计算域示意图 |
 |
图4 升力系数与压力系数对比 |
数值计算结果对比
3 结果与分析
本节的数值计算均以NACA0015翼型为研究对象。与传统运动方式的振荡翼不同, 本文研究的摆式振荡翼置于塔筒尾涡内, 振荡翼的运动频率须与塔筒涡脱频率保持一致。
3.1 塔筒与振荡翼的耦合影响
以塔筒直径D为参照, 取Re=400进行塔筒尾涡与翼型绕流耦合特性分析。取振荡翼弦长c=D, 摆动半径R=6D, 俯仰振幅θ0=80°, 升沉振幅H0=0.8D。表2为有无塔筒工况下的振荡翼水动力性能, 通过对比可以发现塔筒尾涡能够提升水翼的能量捕获效率。水翼对塔筒载荷及涡脱频率的影响见表3, 其中Cd为塔筒平均阻力系数, 
为塔筒最大升力系数, T为塔筒涡脱周期, 发现在塔筒尾涡区域加入振荡翼后, 可以降低塔筒的水动力载荷, 减慢塔筒涡脱, 降低由钝体绕流引起的疲劳载荷。
由(1)~(4)式定义的运动模型, 振荡翼在(0~0.5)T为下行程, (0.5~1.0T)为上行程, 2个行程的运动方向相反, 振荡翼所受载荷也完全相反。为便于分析, 以下行程为例, 结合图5~8对振荡翼的水动力特性进行分析。
在(0~0.1)T期间, 上个周期的上行程形成的前缘涡在水翼上表面的尾缘脱落, 使上表面附近压力逐渐从负压转变为正压; 又由于水翼的逆时针俯仰以及向下的升沉, 使次级前缘涡无法附着在翼面进一步发展, 在上表面形成狭长的涡面并逐渐脱落, 上表面压力逐渐从负压向正压转变; 同时水翼做逆时针的俯仰运动, 且俯仰轴离前缘较近, 尾缘上表面附近的流体受挤压, 压力增大, 尾缘下表面边界层加速分离, 压力减小, 水翼上下表面压差快速转变, 升沉力逐渐从做负功转变为做正功。在(0.1~0.2)T期间, 水翼前缘逐渐进入尾涡涡团, 在涡团与逆时针俯仰的共同作用下, 水翼上表面压力进一步升高, 下表面压力进一步降低, 表面压差增大, 时均功率系数逐渐达到第一个峰值(见图8b))。
在(0.2~0.4)T期间, 受下表面附近塔筒尾涡的涡团影响, 延缓水翼下表面前缘涡的生成与发展, 水翼发生失速延迟, 升力系数变化不大(见图8a))。图7中在0.25T时刻水翼的负压区更大, 0.4T时刻水翼下表面尾缘反向流减弱, 前缘涡发展更充分, 使水翼两侧压差稳定, 升力系数CY维持在较高值。在(0.4~0.5)T期间, 水翼逐渐远离塔筒尾涡, 前缘涡进一步发展并逐渐脱落, 使CY达到第二个峰值(见图8a))。随着水翼的进一步运动, 前缘涡发展至尾缘附近, 在上行程开始时逐渐脱落。
振荡翼置于塔筒尾涡区域内和无塔筒尾涡作用下振荡翼的升力系数曲线和功率系数曲线如图8所示。塔筒尾涡延缓了振荡翼前缘涡的发展, 使升力系数CY与升沉速度VY的同步性更好, 一个运动周期内升沉力大部分时间对翼型做正功, 因而能够获得更高的效率, 从功率系数曲线中也体现了这一点; 尤其在(0.25~0.5)T期间, 塔筒尾涡对振荡翼的功率提升有着明显的积极作用。
有无塔筒情况下的振荡翼水动力性能
有无水翼情况下的塔筒载荷与涡脱周期
 |
图5 下行程振荡翼压力云图 |
 |
图6 下行程振荡翼涡量云图 |
 |
图7 0.25T和0.4T时刻振荡翼压力云图 |
 |
图8 振荡翼的升力系数CY和时均功率系数CP曲线 |
3.2 俯仰振幅
为进一步分析俯仰振幅θ0对摆式振荡翼的影响, 保持其他运动参数相同(Re=400, c=D, R=6D, H0=0.8D), 有无塔筒情况下振荡翼能量收集效率η随俯仰振幅θ0的变化如图9所示。无论是否有塔筒尾涡作用, 当θ0=80°时, 振荡翼的能量收集效率η都达到最大值, 分别为30.72%和27.31%, 在塔筒尾涡的作用下η相对提高了12.49%。塔筒尾涡对能量收集效率有着明显的提升作用, 尤其当θ0=50°时候, 在塔筒尾涡的作用下, 能量收集效率由负值转为正值, 即振荡翼由推进模式转变为能量收集模式。随着θ0的进一步增大, 振荡翼能量收集效率又逐渐降低。
振荡翼运动过程中受力的变化是水动力性能的直接体现, 图10显示了有塔筒情况下不同俯仰振幅下振荡翼力系数在一个运动周期内的变化情况。随着俯仰振幅的增大, 振荡翼的升力系数Cy和阻力系数Cx均有所增大, 当θ0超过80°后, 升力系数不再增大, 而阻力系数继续增大, 使振荡翼的能量吸收效率达到峰值后开始降低。
为进一步分析有塔筒情况下俯仰振幅θ0对振荡翼能量收集效率的影响, 保持其他运动参数相同(Re=400, R=6D, H0=0.8D), 其能量收集效率η随俯仰振幅的变化如图11所示, η不随俯仰振幅或水翼弦长c单调增加或减少, 在参数域内存在高效区间, 在θ0=80°, c=0.9D附近取得。
 |
图9 俯仰振幅θ0对能量收集效率η的影响 |
 |
图10 有塔筒情况不同俯仰振幅θ0振荡翼力系数的变化 |
 |
图11 有塔筒情况不同俯仰振幅振荡翼能量收集效率图 |
3.3 升沉振幅
由(9)式可知, 振荡翼的能量收集效率与时均功率系数CP呈正相关, 与升沉振幅H0呈负相关。图12为有无塔筒情况下振荡翼时均功率系数CP和能量收集效率η随升沉振幅H0的变化曲线(Re=400, c=D, R=6D, θ0=80°), 随着升沉幅度的增大, (8)式中升沉诱导角αh(t)增大, 使最大等效攻角逐渐小于失速角, 时均功率系数先增后减, 在H0=D时达到最大值; 但时均功率系数不与升沉振幅呈相应比例增加, 能量收集效率反而随升沉振幅增大而降低。从图12b)也可以看出, 振荡翼的能量收集效率η随着升沉振幅H0的增大而降低, 且随着振荡翼升沉振幅的增大, 有无塔筒情况下的能量收集效率逐渐接近。
图13为有塔筒情况0.4T时刻不同升沉幅度下振荡翼的压力和涡量云图, 对比发现, 随着H0的增大, 水翼逐渐远离塔筒尾涡涡团, 使涡团对前缘涡发展的延缓作用逐渐减弱; 同时振荡翼摆动速度也随着H0的增大而增大, 前缘涡更靠近下表面且距尾缘更近, 使下表面尾缘附近的反向流增强, 在下表面尾缘附近形成正压, 减小了水翼两面的压差, 从而使水翼的能量收集效率降低。图13表明, 随着振荡翼升沉振幅的增大, 塔筒尾涡区域在振荡翼摆动范围内的占比逐渐降低, 塔筒尾涡对振荡翼的影响逐渐降低, 因此能量收集效率随着H0的增大逐渐与无塔筒情况对应的值接近。
 |
图12 升沉振幅H0对时均功率系数CP和能量收集效率η的影响 |
 |
图13 有塔筒情况0.4T时刻不同升沉振幅H0振荡翼的压力和涡量云图 |
3.4 振荡翼弦长的相关性
由于振荡翼置于塔筒尾涡区域内, 振荡翼运动频率与塔筒涡脱频率保持一致, 因此折算频率f*仅与振荡翼弦长c有关。图14为有无塔筒情况下振荡翼时均功率系数和能量收集效率随弦长的变化(Re=400, R=6D, θ0=80°, H0=0.8D), 从中可以看出, 随着翼型弦长的增大, 时均功率系数和能量收集效率都先增后减。由(9)式可知, 振荡翼的能量收集效率与时均功率系数CP呈正相关, 与振荡翼弦长也呈正相关。无塔筒工况下, 当c=0.9D时η达到最大值30.14%, 对应的折算频率为0.190 9;有塔筒的工况下, 最高效率30.72%在c=D时获得, 此时的折算频率为0.211 5。
在塔筒尾涡的影响下, 振荡翼的最佳折算频率增大, 能量收集性能有显著提升; 且随着弦长的增大, 塔筒尾涡对能量收集效率的积极作用越发明显。在0.4T时刻不同弦长振荡翼的压力云图如图15所示, 对比发现, 随着振荡翼弦长的增大, 前缘涡逐渐靠近前缘且与下表面接近, 尾缘反向流随之增强, 使振荡翼性能下降。有塔筒情况下, 振荡翼下表面的尾缘反向流弱于无塔筒情况, 且随着弦长的增大, 塔筒尾涡对尾缘反向流的抑制作用越明显, 从而使振荡翼能保持较高水动力性能。
 |
图14 振荡翼弦长c对时均功率系数CP和能量收集效率η的影响 |
 |
图15 0.4T时刻不同弦长振荡翼的压力云图 |
4 结论
本文采用LBM-LES方法, 对振荡翼塔筒尾涡潮流能捕获系统进行二维数值模拟, 分析不同运动参数对能量捕获效率的影响。在验证数值方法可靠性的前提下, 分析了有无塔筒情况下, 俯仰振幅、升沉振幅和弦长变化对振荡翼塔筒尾涡潮流能量捕获系统的能量收集效率及水动力性能的影响。二维数值模拟具有一定的局限性, 但是其结果表明:
1) 在塔筒尾涡区域加入振荡水翼后, 可以降低塔筒的水动力载荷, 减慢塔筒涡脱, 降低由钝体绕流引起的疲劳载荷; 在合适的运动参数下, 塔筒尾涡能改善振荡翼的水动力性能, 显著提高能量收集效率。
2) 塔筒尾涡对振荡翼水动力性能的提升作用随着俯仰振幅的增大逐渐减弱, 当θ0=50°时, 在塔筒尾涡的作用下, 振荡翼由推进模式转变为能量收集模式。
3) 不同的升沉振幅下, 塔筒尾涡都能显著提高振荡翼的能量收集效率, 且随着振荡翼升沉振幅的增大, 塔筒尾涡对振荡翼的影响逐渐降低, 有无塔筒情况下的能量收集效率逐渐接近。
4) 塔筒尾涡使振荡翼的最佳折算频率增大, 能量收集性显著提升; 且随着弦长的增大, 塔筒尾涡对能量收集效率的积极作用越发明显, 在c=D时获得最高效率30.72%, 此时的折算频率为0.211 5。
References
- MCKINNEY W, DELAURIER J. Wingmill: an oscillating-wing windmill[J]. Journal of Energy, 1981, 5(2): 80–87 [Google Scholar]
- KINSEY T, DUMAS G. Computational fluid dynamics analysis of a hydrokinetic turbine based on oscillating hydrofoils[J]. Journal of Fluids Engineering, 2012, 134(2): 1–16 [Google Scholar]
- KINSEY T, DUMAS G. Parametric study of an oscillating airfoil in a power-extraction regime[J]. AIAA Journal, 2008, 46(6): 1318–1330. [Article] [NASA ADS] [CrossRef] [Google Scholar]
- YOUNG J, ASHRAF M A, LAI J, et al. Numerical simulation of fully passive flapping foil power generation[J]. AIAA Journal, 2013, 51(11): 2727–2739. [Article] [NASA ADS] [CrossRef] [Google Scholar]
- SU Y, BREUER K S. Resonant response and optimal energy harvesting of an elastically mounted pitching and heaving hydrofoil[J]. Physical Review Fluids, 2019, 4(6): 1–18 [Google Scholar]
- BAIK Y S, BERNAL L P, GRANLUND K, et al. Unsteady force generation and vortex dynamics of pitching and plunging aerofoils[J]. Journal of Fluid Mechanics, 2012, 709: 37–68. [Article] [NASA ADS] [CrossRef] [Google Scholar]
- RIBEIRO B L R, FRANK S L, FRANCK J A. Vortex dynamics and Reynolds number effects of an oscillating hydrofoil in energy harvesting mode[J]. Journal of Fluids and Structures, 2020, 94: 1–15 [Google Scholar]
- XIE Yudong, MA Penglei, WANG Yong, et al. Analysis of energy extraction performance of oscillating foil in non-uniform flow[J]. Acta Energiae Solaris Sinica, 2018, 39(4): 886–891. [Article] (in Chinese) [Google Scholar]
- AKHTAR I, MITTAL R, LAUDER G V, et al. Hydrodynamics of a biologically inspired tandem flapping foil configuration[J]. Theoretical & Computational Fluid Dynamics, 2007, 21(3): 155–170 [CrossRef] [Google Scholar]
- LINDSEY K. Feasibility study of oscillating-wing power generators[D]. Monterey: Naval Postgraduate School, 2002 [Google Scholar]
- JONES K D, LINDSEY K, PLATZER M F. An investigation of the fluid-structure interaction in an oscillating-wing micro-hydropower generator[J]. Physics of Atomic Nuclei, 2003, 71: 73–82 [Google Scholar]
- ASHRAF M A, YOUNG J, LAI J, et al. Numerical analysis of an oscillating-wing wind and hydropower generator[J]. AIAA Journal, 2011, 49(7): 1374–1386 [NASA ADS] [CrossRef] [Google Scholar]
- QIAO Kai, WANG Qixian, WANG Yong, et al. Energy harvesting performance analysis on improved motion model of oscillating hydro-foil[J]. Journal of Shandong University, 2020, 50(6): 40–47. [Article] (in Chinese) [Google Scholar]
- SUN Guang, WANG Yong, XIE Yudong, et al. Hydrodynamic characteristics of oscillating hydrofoil with trailing edge[J]. Journal of Shandong University, 2020, 50(6): 23–29. [Article] (in Chinese) [Google Scholar]
- KARBASIAN H R, ESFAHANI J A, BARATI E. The power extraction by flapping foil hydrokinetic turbine in swing arm mode[J]. Renewable Energy, 2016, 88: 130–142 [CrossRef] [Google Scholar]
- MA Penglei, YANG Zhihong, WANG Yong, et al. Energy extraction and hydrodynamic behavior analysis by an oscillating hydrofoil device[J]. Renewable Energy, 2017, 113: 648–659 [CrossRef] [Google Scholar]
- MA Penglei, WANG Yong, LIU Haibin, et al. Analysis on hydrodynamic performance of flapping foil with swing motion[J]. Acta Energiae Solaris Sinica, 2018, 39(3): 665–672. [Article] (in Chinese) [Google Scholar]
- HE Yaling, WANG Yong, LI Qing. Theory and application of lattice Boltzmann method[M]. Beijing: Science Press, 2009 (in Chinese) [Google Scholar]
- QIAN Y H, D'HUMIRES D, LALLEMAND P. Lattice BGK model for Navier-Stokes equation[J]. Europhysics Letters, 1992, 17(6): 479–484 [CrossRef] [Google Scholar]
- NICOUD F, DUCROS F. Subgrid-scale stress modelling based on the square of the velocity gradient tensor[J]. Flow Turbulence and Combustion, 1999, 62(3): 183–200 [Google Scholar]
- HOU S, STERLING J, CHEN S, et al. A lattice Boltzmann subgrid model for high Reynolds number flows[J]. Fields Institute Communications, 1994, 6: 1–18 [Google Scholar]
- XU Chang, HUANG Haiqin, SHI Chen, et al. Numerical simulation of wind turbine wakes in typical complex terrains based on LBM-LES Method[J]. Proceedings of the CSEE, 2020, 40(13): 4236–4243. [Article] (in Chinese) [Google Scholar]
- RIBEIRO A, CASALINO D, FARES E. Lattice-Boltzmann simulations of an oscillating NACA0012 airfoil in dynamic stall[M]. Mykonos: Springer International Publishing, 2016 [Google Scholar]
- POVITSKY A. Modeling of pitching and plunging airfoils in proximity for thrust generation[C]//46th AIAA Fluid Dynamics Conference, Washington, 2016 [Google Scholar]
All Tables
All Figures
|
图1 摆式振荡翼运动模型示意图 |
| In the text | |
|
图2 D2Q9速度离散模型 |
| In the text | |
|
图3 计算域示意图 |
| In the text | |
|
图4 升力系数与压力系数对比 |
| In the text | |
|
图5 下行程振荡翼压力云图 |
| In the text | |
|
图6 下行程振荡翼涡量云图 |
| In the text | |
|
图7 0.25T和0.4T时刻振荡翼压力云图 |
| In the text | |
|
图8 振荡翼的升力系数CY和时均功率系数CP曲线 |
| In the text | |
|
图9 俯仰振幅θ0对能量收集效率η的影响 |
| In the text | |
|
图10 有塔筒情况不同俯仰振幅θ0振荡翼力系数的变化 |
| In the text | |
|
图11 有塔筒情况不同俯仰振幅振荡翼能量收集效率图 |
| In the text | |
|
图12 升沉振幅H0对时均功率系数CP和能量收集效率η的影响 |
| In the text | |
|
图13 有塔筒情况0.4T时刻不同升沉振幅H0振荡翼的压力和涡量云图 |
| In the text | |
|
图14 振荡翼弦长c对时均功率系数CP和能量收集效率η的影响 |
| In the text | |
|
图15 0.4T时刻不同弦长振荡翼的压力云图 |
| In the text | |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.