大分工業高等専門学校紀要

第 48 号(平成 23 年 11 月)

風レンズ風車の非定常特性に関する研究―正弦的に風速変動する風況下での電圧負荷制御による出力特性―
利光

和彦1・菊川裕規1・佐藤浩平2・松崎純平1・三浦剛1

1機械工学科， 2九州大学大学院

The wind turbines with a flanged-diffuser shroud -so called “wind-lens turbine”- are developed as one ofhigh performance wind turbines by Ohya et al. The performances of the conventional wind turbine and thewind-lens turbines – the compact-type and long-type are experimentally and numerically investigated for
unsteady winds. Experimental and numerical examples are presented to demonstrate the dependences of thewind velocity frequency on the power coefficients of the wind turbines in the harmonic oscillating velocity
wind. Consequently, using the present designed rotor, the wind-lens turbinesshowthe higher performance thanthe propeller type wind turbinewithout the diffuser in the oscillating velocity wind. Furthermore, the numericalestimation can predict the power coefficients for the fluctuating flows in 94 to 102% of accuracy.
キーワード : Wind Turbine, Unsteady Flow, Wind Energy, Power Coefficients, FlangedDiffuser Shroud

１．緒 言る．また，烏谷らの通常プロペラ風車理論出力予測値と比地球温暖化防止のため，自然エネルギーを利用する風車較し予測精度を検討する．の普及が急務の工学的課題となっている．このような時代の要請に答える形で，九州大学の大屋・烏谷ら 1)によって，２．風車モデルと風車トルク測定装置

プロペラ風車につば付の拡大管を設置した「つば付きディ

本実験で用いたコンパクトタイプおよびロングタイプ

フューザ風車（通称，風レンズ風車）」と呼ばれる高効率風レンズ風車の実験モデルおよび概寸を図 1 に，設計諸な風車の開発が進められている．この風レンズ風車を小型

元を表 1, 2 示す．翼車は，設計周速比λ を 3.7 および 5.0

分散型の高効率・低騒音な風力発電装置として市街地などの 2 種 類とし，翼 断面形状は 翼根から 翼端にかけてに設置すれば，風力エネルギー利用の促進が期待される．

NACA63218（翼根）～NACA63212 （翼端）で変化する．

しかしながら風車を設置する市街地はある程度風況の良

風車性能実験装置の概略を図 2 に示す．風洞は測定断面い場所を選定しても風速や風向が時間的に変化する不安

1m×1m で 66 個（横 6 個×縦 11 個）の複数ファンにより定な風況となることが多く，風速および風向が変動する風

変動風や剪断風をアクティブにコントロールできるマルに対する風車性能や特性を把握した上で風車開発を行う

チファン風洞である．翼車の回転軸はトルクメータ（小野ことが重要となる．この点に着目した研究として，例えば，測器製，SS-005 定格 0.5N･m，MP-981）にカップリング

烏谷ら 2), 3)は正弦変動する風速に対して通常の固定ピッチを介して接続され，発電機負荷として直流モータ（Maxon

翼風車の出力について，定常風に対する空力特性と風車翼社製，RE25，定格 20W，40W）を使用する．このモータの動的回転運動応答の特性を組み合わせて，変動風下での負荷トルクを直流電源を用いて調整することで，任意の周風車の理論的な出力予測式を提案している．また，同グル速比に対する翼車回転数および翼車トルクを計測する．

ープは大気乱流が風車性能に与える影響を風レンズ風車フィールド試験によって計測し，的確に変動風を捉えるこ

３．風車周りの定常流れ場の出力特性
風レンズ風車の出力係数 Cw を式(1)で定義する．図 3とができれば，通常の安定した風より出力が大きくなることを示している．4) 本報では，特に風速が正弦変動する場合に着目し，従来型プロペラ風車および風レンズ風車の非定常風速変動に対する出力特性を実験的に明らかにすに流入風速 U =5.0m/s の場合の定常風における Cw と周速比  の関係を示す．この Cw は，式(1)において翼車の直径を基準とした円の面積（ A  2.96  10 2 m2)を受風面積―4―大分工業高等専門学校紀要

第 48 号((平成23 年 11 月)

出力係数は，翼車単体，A とした．定常風における風車の出力係数の追随特性に指針を与えて
f r [Hz]を次式で定義し，翼車のコンパクトタイプ，ロングタイプの順にに大きくなり，翼車いる．

設計周速比λ =5.0 では，各タイプでそれぞれそれぞれ最高出力の周速比は，λ=2.5, 3.7, 5.0，翼車設計周速比λ
翼車設計周速比 =3.7 では，Table 1 Dimensions of the wind turbine of theflanged-diffuser
diffuser shroud.
Tip clearance
3mm (ht/D=1.5%)
ht
Diffuser typeLong
Compact

λ=2.0, 3.0, 4.7 である．なお，この出力係数
出力係数 Cw から後述
する第 5 章での変動風における風車の出力特性値

C w ( w ) を計算する．

Diffuser angle

PCw 
1 / 2AU 3

DiffuserlengthFlange height(1)４．正弦的風速変動風における風車出力特性
風車出力特性12 deg.

－L/D

1.250.225

h/D

0.50.1Table 2Dimensions
imensions of the rotor blades.

次に，上流風速を正弦波状に変化させたさせた変動風U  U  u~ sin( 2f  t )Rotor radius
Hub to tip ratior
Dh / D

Designed tip-speedratiorU(2)に対する風車性能実験について述べる．
．表3に実験条件を示す．平均風速は U =5m/s，変動風速のの振幅は u~ =1m/s

0.2
5NACA63218 (root)～
NACA63212(tip)Blade profiles（±20%），変動周波数 f = 0.033Hz，0.05Hz
0.05Hz，0.083Hz，97 mm

0.25Hzとした．
烏谷ら(4)は，プロペラ風車の定常風におけるにおける静的な空力

Pitot TubeRotor

Flanged-Diffuserr

Torque Detector(ONOSOKKI SS-005)

DC Motor(Maxon RE25)

特性を基に，翼周りの流れが風速の変化変化に応じた状態にな
Steady FlowandUnsteady Flow

□1000

るために必要な時間（風が翼を通過するする時間のオーダー）る準静的近似が成立する条件下で，風車風車の応答周波数

1400500は，翼車が 1 回転する時間に比べて十分小十分小さいと仮定でき
1000z

ht=3

Wind Tunnel

Torquue Converter(ONOSO
OKKI TS-3200A)

Variable Resistor

PC

DC Power Supply

Fig. 2 Schematics of measuring
easuring system of windturbineperformance.

h=20

Dh=40

φ206

27
L=45

D=200
φ226
DF =φ266

x

θ=12°
L=250

φ306
DF =φ506
h=100

Dh=40

D=200
ht=3(a) Compact-type

Fig. 3 Power coeffcients based on rotor swept area for(b) Long-type
Fig. 1 Schematics and photograph of the wind turbines
turbinesteadywind velocity,5m/s.

with the flanged-diffuser shroud.

―5―大分工業高等専門学校紀要

第 48 号(平成 23 年 11 月)

Table 3 Conditions of fluctuating winds and turbine

U (t )  U  u(t )characters.
Uu~ [m/s]

Oscillation Frequency, f [Hz][m/s]

( u~ / U %)(Oscillation Period, T  1 / f [s] )1.0 (20%)

0.25

0.083

0.05

0.033

(4)

(12)

(20)

(30)

w

Wind turbines
5.0変動風下での時間平均出力 PA は，平均する時間を T ，翼車トルクおよび角速度を T w ， とすると次式で表され，風速の 3 乗に比例する．

frD

3.7

5.0

3.7

5.0

Rotor only

2.09

2.47

0.65

0.64

Compact-type

3.06

3.51

1.01

0.87

Long-type

4.86

4.74

0.98

1.03

3r 4Cw (w )U
fr 
2 I2w(4)

PA 

1 T
1 T
1Tw dt   Cw ( (t )) A[U (t)]3 dt (5)T 0T 02式(4)を(5)に代入して整理し，変動風速 u (t ) と出力係数C w (  ( t )) は無相関を仮定すると，非定常風における出力係数 Cw は式(6)で近似される．

C w  C wA 1  3 2

 (6)

(3)

2 

ここで，  は流体密度， r は翼車半径， C w (  w ) はトルク係数が周速比のみに依存すると仮定した出力係数，  w は1 T
C w (  ( t )) dt
T 0
1

C w ( ) d (8)
2 g ( f / f r )  C wA 

定常風下において最高効率を示す周速比， U は時間平均風速， I は翼車の慣性モーメントを表す．なお，本実験風車では I =1.52×10-4 kg  m 2 である．実験条件（表 3）に風車の応答性の参考値として f r 値を示す．なお，本実験  w  g ( f / f r )  h( f / f r ) ,での準静的仮定成立条件の詳細は紙面の都合で省略するが，翼端通過時間および軸方向翼車通過時間が共に翼車が一回転する時間より十分小さいことが計算され，理論上で   w 2 ,は準静的仮定が成立していると考えられる．

比 λD =5.0 では，それぞれの風車タイプにおいて平均周速領域に渡って変動風のほうが定常風より出力係数が大きい．一方，ロングタイプ風レンズ風車では  <5.0 の場合

  w 2 2 ,(10), (11)h  f / f r   1.371  exp  3.76 f / f r  , (13)イプおよびロングタイプの出力係数を示す．翼車設計周速

よびコンパクトタイプでは， f ≦0.083 の場合に，ほぼ全(9)g  f / f r   1 .04 1  exp  7.49 f / f r  , (12)

図 4～9 に時間平均出力係数の一例として U =5m/s で20% 変動風振幅(1m/s)における翼車単体，コンパクトタ比  が約 2.5，3.5 および 4.7 で最大となる．翼車単体お221 T  u (t ) 1  u~ (7)
dt


T 0  U 2 U 

ここで， w は定常風での出力係数が最大となる周速比，

 ，  は風速の分散に対して式(10), (11)で関係づけることができる．また， g  f / f r  と h  f / f r  は，烏谷ら 3)は定常風の場合が変動風の場合より大きく，  >5.0 の高によると風車の慣性モーメント，翼半径，出力係数や制御回転域では，全周波数において変動風の方が大きくなる．

間隔にかかわらず，式(12)，(13)の近似式で表すことがで

但し， =4.3 付近では翼車回転共振点のために変動風（ f= 0.05Hz，0.083Hz，0.25Hz）の C w が不連続に小さくなっているため，ロングタイプは全周速比において該当する
変動風の場合が定常風の場合より C w が大きい傾向があると考えられる．

きる．なお， f は風速の変動周波数，fr は式(3)の応答振動数である．以上より，正弦変動する場合は，推定した CwAより出力係数が式(6)右辺第 2 項分 3 2 だけ大きくなる．
本実験値と理論値を比較するため式(14)で予測精度を検討する．最高効率点での理論値と実験値との比較を図
10 にまとめる．同色の棒グラフは左から，翼車単体，コンパクト，ロングタイプを表す．結果より烏谷らの理論式

５．変動風での最高出力係数における理論推定値と実験値との比較による予測は，本実験条件では 94%～102%の精度で実験値と一致することが分かる．

一般的に，変動風風速 U (t ) は平均風速 U と変動風速u (t ) の線形結合として次式で表される．

―6―CwEC A wE 2CwT Cw (1  3 )(14)

大分工業高等専門学校紀要

第 48 号((平成23 年 11 月)

Fig. 4 The effect of the wind frequencies f on the time
averaged power coefficients Cw of a rotor only turbine for U=5m/s, u~ =1.0m/s , i.e. u~ / U =20% with designed turbine

Fig. 7 The effect of the wind frequencies f on the time
averaged power coefficients Cw of a rotor only turbine for U=5m/s, u~ =1.0m/s , i.e. u~ / U =20% with designed turbine blade

bladeλD=5.0.

λD=3.7.

Fig. 5 The effect of the wind frequencies f on the time
averaged power coefficients Cw of a compact-typecompactturbine for
U =5m/s, u~ =1.0m/s , i.e. u~ / U =20% with designed turbine

Fig. 8 The effect of the wind frequencies f on the time
averaged power coefficients Cw of a compact-type turbine for
U =5m/s, u~ =1.0m/s , i.e. u~ / U =20% with designed turbine

bladeλD=5.0.

bladeλD =3.7.

Fig. 6 The effect of the wind frequencies f on the time
averaged power coefficients Cw of a long--type turbine for U=5m/s, u~ =1.0m/s , i.e. u~ / U =20% with designed turbine

Fig. 9 The effect of the wind frequencies f on the time
averaged power coefficients Cw of a long-type turbine for U=5m/s, u~ =1.0m/s , i.e. u~ / U =20% with designed turbine

bladeλD =5.0.

bladeλD =3.7.

―7―大分工業高等専門学校紀要

第 48 号((平成23 年 11 月)

６．定常風と変動風における最大出力係数
最大出力係数の比較
定常風の最大出力係数CwCと変動風のの最大出力係数CwWを比較するために，定常風の最大出力係数
最大出力係数に対する変動風出力の増加割合εを式(15)で定める．

C wW  C wC
100 [%]
CwC(15)

図11, 12に定常風と変動風での出力係数出力係数を比較した結果εの関係を示めす．
プロペラ単体では，設計周速比λD =3.7，，5.0ともに，周波数f=0.033，0.05Hzでは定常風の出力に比べべ変動風が大きく，周波数f=0.083，0.25Hzでは変動風の出力出力が小さい．プロペ
ラ単体において最も変動風の増加割合がが大きくなったのは，λD=3.7の翼車，周波数f=0.033Hzの場合場合であり，その増Fig. 10 The comparisons of the experimental andnumerical maximum power coefficients of the windturbines. The same color three bar graphs present thecases of rotor only, compact--type and long-type fromleft to light for each condition.

加割合は約6.7%であった．
コンパクトタイプでは，設計周速比λλD=3.7，5.0ともに，周波数f=0.033，0.05，0.083Hzでは定常風定常風の出力に比べ変動風が大きく，周波数f=0.25Hzでは変動風変動風が小さい．コンパ
クトタイプにおいて最も変動風の増加割合増加割合が大きくなったのは，周速比λD =3.7の翼車，f=0.033Hz=0.033Hzの場合で，その増加割合は約8.2%であった．
ロングタイプでは周速比 λD =3.7 において，すべての周波において数で定常風の出力に比べ変動風が上がっがった．周速比 λD=5.0 においては，周波数 f=0.033，0.050.05，0.083Hz では定常風の出力に比べ変動風の方が大きく，，周波数 f=0.25Hzでは変動風が小さい．ロングタイプにおいて最も変動風のロングタイプにおいて増加割合が大きくなったのは，周速比 λD =5.0 の翼車で周波数 f=0.05Hz の場合であり，その増加割合増加割合は約 6.1%であった．

Fig. 11 Comparisons of the maximum power
coefficients in steady and unsteady flows for designed
tip speed ratio λD=5.0.

また，周期の長い変動風から周期の短短い変動風になるにつれて，定常風に対する変動風の増加割合増加割合が小さくなる傾向にあることがわかる．これは，周期のの短い変動風に風車が追従できていないためと考えられる．
．
以上の結果より，最も定常風に対するする変動風の増加割合が大きかったのは，コンパクトタイプのコンパクトタイプの約8.2%であった．
しかし，ロングタイプはプロペラ単体やコンパクトタイプやコンパクトタイプに比べ，変動風の周期が短くなっても増加割合増加割合がそれほど小さくならないことがわかる．よって，，ロングタイプは最も，周期の短い変動風に対する風車の追従性追従性が高いと言える．
また，コンパクトタイプはλD =5.0の場合場合に比べλD =3.7が増加割合が大きいのに対し，ロングタイプはロングタイプは周速比λD =3.7よりもλD =5.0が大きい．λD =3.7の風車翼
風車翼はコンパクトタイプ用，λD =5.0の風車翼はロングタイプ用用に設計したものでFig. 12 Comparisons of the maximum power
coefficients in steady and unsteady flows for designed
tip speed ratio λD =3.7.

ある．このことより，それぞれの風車タイプタイプに合わせて設計した翼車を用いると，出力は若干小さいがさいが，定常風を基準とした変動風時の性能向上は期待できるできると考えられる．

―8―大分工業高等専門学校紀要

７．結

第 48 号(平成 23 年 11 月)

言

本研究では，プロペラ単体風車，コンパクトタイプおよびロングタイプ風レンズ風車に対して周期的変動風下で理論値と実験値の出力係数は，本実験風車，実験条
件において94%～102%の精度で一致する．



周波数f=0.25Hzのような高い周波数の変動風では翼車が変動風に追随することができず，理論値に対し翼車単体と風レンズ各風車の出力測定を行った．実験の条て実験値の出力係数が若干小さくなる．

件は時間平均風速を5m/s，振幅1m/s（20%）とし，変動風の周波数はf=0.033Hz，0.05Hz，0.083Hz，0.25Hzとした．
また，変動風下での出力係数の理論値と実験値を比較し理謝論予測精度を検証した．各種の実験から得られた結果をま辞本研究に対して，九州大学大学院応用力学研究所大屋裕

とめると以下の通りである．

二先生ならび烏谷隆先生に貴重なご助言を頂きました．ま【風車翼設計周速比の特性】

た，文部科学省科学研究費基盤研究（C）（課題番号：



出力は翼車の設計周速比3.7，5.0では，5.0の方が大き21560823）および九州大学応用力学研究所共同利用・共い．

同研究経費（課題番号 22ME-10）の援助を受けて行われ定常風での最大出力周速比に対応した翼車を使用すました．ここに記して謝意を表します．

ると，定常風の出力に対する変動風の出力の増加率が大きくなる．すなわち，変動風に対応しやすい風車となる．

参考文献(1) Ohya, Y., and Karasudani, T., A Shrouded Wind Turbine【定常風と変動風での性能】

Generating

定常風の最大出力係数に対する変動風の最大出力係数の増減比較の結果より以下のことが言える．
Output

PowerwithWind-lensTechnology, Journal of Energies, 2010, 3, pp.643-649.

(2) 烏谷，大屋，渡辺, 小型風車の風速および負荷の変動

低周波数変動風f=0.033Hz，0.05Hzでは，全タイプ実に対する応答について，日本風力エネルギー協会誌,

験風車において定常風に対して変動風の出力が増加

Vol. 31, No.2, 2007, pp.120-123.

する．特に，コンパクトタイプ風レンズ風車は定常(3) 烏谷，大屋，渡辺, 正弦変動する風速に対する固定ピ風に対する変動風の増加量が最も大きい．
High

ッチ翼風車の出力について，日本風力エネルギー協風レンズ風車は翼車単体より高周波数変動風に追随

会誌, Vol. 32, No.1, 2008, pp.128-132.

しやすい．特に，ロングタイプ風レンズ風車は，プ(4) 沖野, 伊庭, 烏谷, 大屋，渡辺, 大気乱流が風車の発ロペラ単体やコンパクトタイプに比べ高周波数変動電性能に与える影響, 第 29 回風力エネルギー利用シ風に対応できる．

ンポジウム講演論文集, 2007, pp.310-313.

【理論値と実験値との比較】
最高出力周速比における実験値と理論値との比較による結果から以下のことが言える．

―9―(2011.9.30 受付)