２次元数値解析による浮体型波力発電装置特性の妥当性と最適形状の検討 ＊鈴 木 正 己 *1 ， 窪 木 利 有 *2
永 田 修 一 *3 ， 瀬戸口俊明 *3
Investigation of Validity of Characteristics and Optimal Profile
Predicted by 2-D Numerical Analysis of Floating Type Wave
Power Generating System with OWC
Masami SUZUKI *4 , Toshiari KUBOKI,
Shuichi NAGATA and Toshiaki SETOGUCHI*4Department of Mechanical Engineering, Graduate School of Engineering, The University of To-

kyo, 7−3−1 Hongo, Bunkyo−ku, Tokyo, 113−8656 Japan

A floating type wave power generating system has an Oscillating Water Column (OWC). The device
captures the wave energy using the heaving, pitching and surging motion of the device and the heaving
motion of OWC. The two-dimensional numerical method for analyzing a floating body with OWC type wave
energy conversion device is introduced where the eigenfunction expansion method is described under thecondition that the linear water wave theory is applicable. It is confirmed that these solutions give good
agreement with the several experimental results in the previous paper. The two-dimensional optimal profiles are eliminated according to the different conditions from the high efficiency, i.e. the minimum size per aunit power and so on.

Key Words :1. 緒

Oscillating Water Column, Wave Power Generating System, Numerical Analysis, Fluid machinery, Ocean Engineering

言

ー変換率と高い空気室内圧力が得られる形状の解明が理論と実験により進められている．

振動水柱(Oscillating Water Column : OWC)型の波力発電装置は空気室と空気タービンにより構成されている．空気室は波浪エネルギーを往復空気流のエネル
ギーに変換する装置であり，空気タービンはこの往復流により駆動される．空気タービンは構造が簡単で，往復流中でも常に一方向に回転できるウェルズタービン(1) が微小振幅波の線形理論を基礎とした解析結果(3) ～ (8)は実験データのバラツキなどを考慮すると実験結果との一致は比較的良好である.ただし，個々の空気室形状に対する特性は模型実験による確認も重要であり，実験と数値解析のそれぞれの利点を活かしつつ研究を進めることが求められている．

主に用いられている．このタービンは一定回転数の下では流量とタービン差圧が比例関係にあり(2) ，空気室からは線形負荷とみなすことができる．空気室特性として，
波浪から空気流へのエネルギー変換率（２次元に限り
効率とも一致）や空気室内の圧力振幅，室内水面の上下動振幅が重要であり，より小さな浮体で高いエネルギ＊

原稿受付 2007 年 02 月 28 日．
*1
東京大学大学院工学系研究科（〒113-8656 東京都文京区本郷7-3-1）.
*2
太平洋学会（〒108-0073 東京都港区三田 4-15-29 三田東急アパート 3F）.
*3
佐賀大学海洋エネルギー研究センター（〒840-8502 佐賀県佐賀市
本庄町一番地）.
E-mail : wells@cfdl.t.u-tokyo.ac.jp

本研究では数値解析により浮体型波浪エネルギー
変換装置の特性を把握し，最適浮体形状について検討を行う．まず，数値解析結果の妥当性を確認するた
め，解析対象とする浮体の実験結果と比較し検証を行う．
これらの検証を基に，エネルギー吸収特性に加え，建造コスト，すなわち，浮体の小型化も視野に最適浮体形
状について検討を行う．
2. 数値解析による推定精度の検証
2・1 検証対象浮体

海洋科学技術センター（現，

海洋研究開発機構）が実証実験（9），(10） に先駆けて実施した基礎実験の２次元浮体形状について，領域分割法－1－

２次元数値解析による浮体型波力発電装置特性の妥当性と最適形状の検討

による数値解析(6) を行い，実験結果との比較により，解出されており，定式化の詳細は前報(3),(6) に記述されて析精度の検討を行う．図１は海洋科学技術センターで

いる．

(7)

実施された実験装置の概略図を示している ．実験は本論の数値解析および実験は規則波を対象にして水槽に浮体を固定した状態と弛緩係留した状態の２種いるため，計算における入力条件は入射波の波高（ま類で実施されている．係留時の浮体は弱い２本のバネ

たは振幅），係留索のバネ定数，空気室（OWC）の負荷で弛緩係留され，無拘束に近い状態になっている．また，ダンピング係数である．ここでは線形問題として取扱っ浮体の重心は浮心の付近となるように調整が行われてているため，負荷ダンピング係数は定数で，その定義式

いる．波高計は船首前方 6m に入射波用，船尾から後は方 1m に透過波用が設置され，反射波は船首前方 1m

D = AW p / zに設置した波高計を上流側へゆっくりと移動し，波の腹で表され，AW：空気室内の水線面積，p：空気室内のと節を計測することにより算出している．波浪エネルギー

圧力， z ：空気室内の水面を空間平均した上下速度では浮体前面の空気室で吸収し，往復空気流に変換され

ある．

(1)

実験では空気タービンの代わりに，空気室の負荷とる構造となっている．空気室内の圧力は浮体前面の空して取付けの容易なオリフィスが従来から用いられてい気室上部に設けた圧力孔で計測されている．
模型が小型であるため，空気タービンの代わりに空る．オリフィスの圧力損失は流量の２乗に比例するため，気室の負荷としてオリフィスが用いられ，実験が行われ非線形負荷となる．このため，線形負荷と２乗負荷の平ている．空気室の負荷特性はタービン直径や回転数に均エネルギーが一致するように等価線形化する．オリフ

より変化する．これを評価するため，空気室上部にある

ィスの圧力損失 p はオリフィスの開口比（オリフィスの絞り面積比）を変化させ

p=て，負荷特性の影響を検討している．実験は３種類の開1
C p ρ a z 22(2)

口比 1/120（オリフィス径 φ 14mm），1/150（ φ 12.6mm），で表すことができ，これを式(２)に代入し線形化すると負1/200（φ 10.7mm）について実施されている．入射波の目

荷ダンピング係数 D は

標波高は 20mm と 30mm の２種類で，波周期は 0.5～D=

1.5 秒まで 0.1 秒刻みで計測されている．

8 1
C p ρ a z AW
3π 2

(3)

船尾にある３つの空気室は１室につきφ 25mm の穴がと表される．ここで，Cp：オリフィスの圧力損失係数，ρa：

８個開いており，エネルギー吸収用ではないため各室

空気の密度， z ：空気室内水面の上下速度の振幅，の開口比はそれぞれ 1/20 と大きい．

8/3π は等価線形化の補正係数である．

(6)

図１および図２に示す浮体

タービンを負荷とするとき種類（ウェルズタービンや衝型波浪エネルギー変換装置の特性は前後揺れ，上下

動タービンなど）や運転方法によりダンピング係数が異

揺れ，縦揺れ，４空気室の室内水面の上下変動に対すなってくる．定回転運転するときウェルズタービンはターる合計７つの運動方程式を連立して解かれる．これによ

ビン差圧と流量が比例するため，線形負荷となりダンピ

り，前後揺れ ξ ，上下揺れ ζ ，縦揺れ θ の変位，エネル

ング係数は波高の影響を受けないが，衝動タービンはギー吸収を行う船首空気室内の水面の上下変位 z と後オリフィスに近い特性になるため，波高の影響を受ける．

部３室の水面の上下変位が算出される．運動方程式中また，ウェルズタービンも回転数が変化するときには非には入射波により浮体や空気室に加えられる波浪強制

線形負荷となる．このように非線形負荷は特性を評価し力，浮体や空気室内の水面が動揺することにより浮体

辛いように思えるが，オリフィスのように等価線形化してや空気室に作用する流体力が現れる．この流体力は微負荷ダンピング係数を表すことにより統一的な評価が可

小振幅波の線形理論に基づいた領域分割法により算能である．実際の洋上では不規則に波高や周期が変2・2 数値解析化するため，波浪エネルギー変換装置はこの不規則波
60.5Wave Generator

Wave Probe

中で駆動される．そのため，著者ら (11)は実証実験の解Spring1

1

Wave Probe

Wave Absober
Model

Wave Probeに近くなることに着目し，タービンのダンピング係数を等価線形化し，不規則波特性を解析する方法も提案してWave Tank
11.5
20

析結果からタービン流量の確率密度分布がガウス分布

Unit : mいる．このような等価線形化による方法は設計などにおいても実用的で有用な方法である．

Fig. 1 Experimental apparatus－2－２次元数値解析による浮体型波力発電装置特性の妥当性と最適形状の検討

丸印，計算結果は破線（拘束状態）と実線（係留状態）

時間平均した単位幅当たりの空気室出力WOWC は，で表示している．また，数値解析は実験と同じ波周期に

空気室奥行きを B とし，

WOWC =D ωz

2

2

ついて計算し，折線により表示を行っている．入射波の(4)B

波高は設定した目標波高から多少ずれるため，計測さで表され，単位峰幅当たりの入射波の波パワーWi は，Wi =f ( k0 h ) =れた入射波の振幅 a i / B も図中に添えている．入射波以外の計測データは空気室内の圧力振幅 p / ρW gai ，ρW g 2 ai 2
4ω f (k0 h)

反射波の振幅 ar / ai ，透過波の振幅 at / ai である．図(5)2 cosh 2 (k0 h)
2k0 h + sinh(2k0 h)

中にはこれらに加えて数値解析から得た空気室内の水面変位の振幅 z / ai も示している．浮体の動揺特性に関しては前後揺れ(surging) ξ ，上下揺れ(heaving) ζ ，であることから，空気室効率（エネルギー吸収効率）

縦揺れ(pitching) θ を示し，重心まわりの振幅( ξ , ζ , θ )ηOWC は，と位相( ε ξ,ε ζ ,ε θ ) により表している．位相は空気室前面

W
ηOWC = OWC
Wi（船首）位置における入射波の上下変位を基準に表示(6)している．縦揺れは時計回りの角度を正とし，波形勾配となる．ここに，D ：タービンやオリフィスなどの負荷ダン（波高と波長の比）で無次元化している．実験による波

ピング係数，ω ：入射波の角周波数，z：船首空気室内水

高および圧力の読取り誤差は 1mm，1mmAq（10Pa）以面の上下変位，a i ：入射波の振幅，k0 ：入射波の波数，下である．しかし，これを超える実験結果のバラツキが

h ：水深である．

生じており，より大きな実験系の総合的誤差が内在する

2・3 数値解析結果の検証(7)

図２は数値解析結ことがわかる． このため，数値解析で推定される数値果を検証するため示した１例であり，入射波の設定波高

誤差 (3),(6) （0.5% 以下）とともに実験誤差を踏まえた特性

20mm(a i / B = 0.063 )，開口比 R=1/150 について供試浮把握を行っている．
実験範囲にある波長 λ / B=2～18 の間では数値解析

体の実験結果と比較されている．実験結果は白丸と黒

0.1

40

4

6

0

8 10 12 14 16 18 20

0

2

4

6

8 10 12 14 16 18 20

1.4
1.2
1.0
0.8
0.6
0.4
0.2
0.0
360

θ / (2ai/λ)

2

180
90

180
90

0

2

4

6

8 10 12 14 16 18 20
λ/B

0

1.4
1.2
1.0
0.8
0.6
0.4
0.2
0.0
1.4
1.2
1.0
0.8
0.6
0.4
0.2
0.0
0

B

2

4

6

8 10 12 14 16 18 20

2

4

6

8 10 12 14 16 18 20
λ/B

0.7
0.6
0.5
0.4
0.3
0.2
0.1
0.0
180
90

εθ(°)

270

εζ(°)

270

εξ(°)

60
20

1.4
1.2
1.0
0.8
0.6
0.4
0.2
0.0
360

0

at / ai

80

η (%)

0.0
1.4
1.2
1.0
0.8
0.6
0.4
0.2
0.0
0

ζ / aiξ / aip / ρwg ai

ai / B

0.2

1.4
1.2
1.0
0.8
0.6
0.4
0.2
0.0
100

a r / ai

Condition
Fixed
Mooredai /B≅ 0.063, R≅1/150

z / aiExp. Cal.

0
-90

0

2

4

6

8 10 12 14 16 18 20
λ/B

-180

0

Fig.2 Pressure in the air chamber, elevation of the water surface in the air chamber, efficiency, reflected waves,transmitted waves, surging, heaving, and pitching of the floating body against the wavelength, λ , at the opening
ratio of the air chamber, R = 1 /150 .

－3－

２次元数値解析による浮体型波力発電装置特性の妥当性と最適形状の検討の結果（実線および破線）は実験データのバラツキ程度動揺したときの効率は急激な低下を起こしている．ここに収まっており，空気室内の圧力 p，反射波 ar，透過波で重要となる運動は上下揺れと縦揺れであり，上下揺at，効率ηOWC は波の振幅 ai や空気室の開口比 R にかれの振幅 ζ は最大効率点で極大値となっている．このとかわらず，実験結果に匹敵する計算精度が得られていきの位相は上下揺れε ζ ，縦揺れε θ ともに−90°であり，合る．ただし，短波長（ λ / B < 8 ）側では実験と計算結果の成された浮体前面の上下変位の位相は−90°となるため，差が増大する傾向にある．これは短波長や高波高にな

浮体前面では入射波に対して位相が 90°遅れて運動をる程，波形勾配が大きくなり微小振幅波の線形理論にしていることがわかる．図中には示していないが，このとよる推定精度が低下するとともに，実験では波の周波数

き空気室内水面の上下変位は固定と動揺状態ともに位が高くなるため波形が変形し易くなるなど，生成波形や相が−80°前後になっている．

計測など実験系全体の誤差も増大し数値解析と実験結

基本的に浮体が固定された状態では，波長（波周果との差が増す要因になっていると考えられる．また，期）の影響が少なく，空気室奥行き B に対する波長λが浮体の角部からのはく離による損失の増加も影響して３倍から 12 倍の間は 50％以上の効率ηOWC が得られていると予測される．

いる．浮体が動揺すると最大効率点よりも長波長（長周空気室内の圧力 p は波長 λ / B = 8 付近に最大値が期）側では効率が急激に低下する．これは浮体の波乗

現れ，負荷ダンピング係数の増加とともに圧力が増加す

りが顕著に現れ，縦揺れの位相が入射波と同位相に近る傾向を示すものの，図２以外の実験結果や計算結果

づき，波が透過しやすく，透過波の波高が急激に増大を考察すると波の振幅 a i / B や開口比 R による負荷ダするためである．最大効率点よりも短波長（短周期）側でンピング係数の顕著な影響は受けていない．

は，入射波と空気室内の振動水柱との連成運動が支配

最大効率点は波長が空気室奥行きの 6～8 倍，的になっており，動揺による透過波の増大はあるものの

λ /B=6 ～ 8 付近に現れ，図２に示す波の振幅 a i /B ≅反射波に顕著な違いはなく，空気室内の圧力 p や効率

0.063 ，開口比 R=1/150 では 80%と高い効率ηOWC が得

ηOWC は固定時と同様の傾向を示している．

られている．この波長では実験が行われたすべての波

空気室効率について考えると，動揺時においても高，開口比の条件において，反射波および透過波の振80％程度の高い効率が達成されている．そのため，今幅が極値となり，一旦小さな値になって効率が向上して後重要な点は，さらに長い波長で高い効率が得られる

いる．λ / B=8 より大きな波長では透過波の振幅は急激

浮体形状の開発である．すなわち，建造コストを低く抑に増大することになる．これに対して，反射波は波長がえるためには，より短い浮体で高い効率を得る必要があ

λ / B=8 より短くなると，その波高が大きくなる．このためる．そして，狭い同調点付近のみでなく，広い範囲で高λ / B=8 より短い波長で効率は低下する．また，開口比

い効率が得られることも必要である．これを踏まえて，次が小さくなるにつれて，負荷ダンピング係数が過大にな

章で種々の形状について検討し，より適切な浮体形状るため，効率は低下するようになる．線形負荷であればの提案を行う．

波高により負荷ダンピング係数が変化しないため，微3. 浮体形状による特性比較

小振幅波の線形理論が成立する範囲では波高が変化
例えば，無限水深で波周期が８秒の波は波長が 100mしても効率は変化せず一定となる．しかし，本研究は，になり，図２の浮体形状では最大エネルギー吸収効率

非線形負荷であるオリフィスを用いているため，入射波

点λ /B=8 で運転すると，約 75m もの浮体長が必要となの振幅が増大すると負荷ダンピング係数が大きくなり，る．有義波高 2m，有義周期 8sec の波浪エネルギーは効率の低下につながっている．

単位幅あたり約 16kW/m であることを考慮すると，エネ

図２の下部には浮体の動揺特性を示している．浮体

ルギー吸収効率を優先した設計は，加えられる波浪エ

動揺は重心まわりの振幅(ξ , ζ , θ ) と位相( ε ξ ,ε ζ ,ε θ ) によ

ネルギーに比較して非常に大きな浮体を生み出すことり表している．位相は空気室前面（船首）位置におけるになる．このため，高いエネルギー吸収効率を優先する

入射波の上下変位を基準に表示している．また，εθ はのではなく，浮体サイズを小さくし，建造コストの低減と±180°の範囲で表示しているため，表示上±180°付近でエネルギー吸収効率を加味した総合的な評価方法のは 360°の跳びが生じる．縦揺れθ は時計回りの角度を確立が求められる．このため，２つの評価関数

正とし，波形勾配（波高と波長の比）で無次元化している．最大効率点は浮体の固定時，動揺時にかかわらず，
波長が空気室奥行きの 6～8 倍付近に現れているが，8倍を越えると浮体が固定された状態に比較し，自由に－4－

2
Ld / ( λMAXηOWC )(7)

λMAX / λMIN(8)

２次元数値解析による浮体型波力発電装置特性の妥当性と最適形状の検討を導入する．浮体サイズを表す重要な因子として，浮体る．浮体内に空気室が存在する形状で，空気室を浮体

断面積 Ld，波長λ M A X ，エネルギー吸収効率ηOWC を挙前方に設置した形状と後方に設置した形状について検げることができる．同一浮体形状でのエネルギー吸収効討を行っている．ここで，空気室奥行き B，浮体長 L，喫率は浮体長と波長の比および浮体形状の影響を受ける．

水 d，水深 h としている．また，x1 の位置を前方，後方に浮体の喫水と波長の比は浮体のサイズと浮体形状を表移動することで，空気室が船首，船尾に設置された浮す因子でもある．そして，浮体形状はエネルギー吸収効

体形状を表現できる．
図４は船尾に空気室を有する浮体特性を示している．

率に影響をおよぼし，高いエネルギー吸収効率では同じエネルギーを得るための浮体幅を狭くでき，結果とし

水深 h が空気室奥行き B の 40 倍，喫水 d を３倍とするて浮体を小型化できることになる．これらを考慮し，浮体浮体形状について，浮体長の影響を検討している．浮

サイズに関する評価関数を導出すると式(7)が得られる．

体長 L は空気室奥行き B の 3～12 倍まで変化させ，波式(8)は運転範囲を表し，大きい値ほど広い運転範囲が長ごとに負荷ダンピング係数を変化させて最大エネル得られることを示している．本研究における運転範囲の

ギー吸収効率を求め，その負荷ダンピング係数 D，浮選定は解析結果から波長の影響を考慮し，エネルギー

体と空気室内水面との相対変位振幅 zOWC，空気室内

吸収効率が 50%以上としている．なお，運転範囲中で最圧力振幅 p，エネルギー吸収効率ηOWC，浮体重心位置も短い波長λ M I N ，最も長い波長λ M A X としている．

の上下変位振幅 ζ ，縦揺れ変位振幅 θ を示している．
最大エネルギー吸収効率は全長 L/B=4 の浮体で生じ，

図３は検討対象とした基本的な浮体形状を示してい波長 と空気室奥行きの比はλ /B=22 付近となっている．

z
x0

ae

なお，作動範囲は浮体長が L/B=8，12 と長くなるほど拡x3

Lx1 B x2

大する．
ここで，エネルギー吸収効率が 50%以上を作動範囲

iω tとして，式(7), (8)の評価関数により最適浮体形状を考察
xdhTable 1. Typical effects of the length for the device with OWC atstern side, where is over 50% in the efficiency.

-d1
-d2

L/B
3
4
8
12

Fig.3 Schematic view of floating device with OWC.

λ/B

L/λMAX

2
Ld / λMAX
19.5 ∼ 28.3
23.0 ∼ 44.5
29.0 ∼ 47.0


0.141
0.180
0.255


0.01490.01210.0163λMAX/λMIN

1.45
1.92
1.62

L

L/B=4

zOWC / ai

L/B=3

3

4

2

8

12

1

1.0

8

12

3

0.5
0.0
5

ζ / ai

L/B=3

3
2

4
8

0

0

10

20

8

30

40

50
λ/B

60

40

3

0
80
60

L/B=3

4

40

0

8
12

0

Fig.4 Effect of the length of floating devices which have OWC at stern side.

－5－

12

60

20
12

3
12

L/B=4

80

20

1

L/B=12

3

0
100
L/B=4

8

2

1

4

L/B=8

L/B=4

4

0
1.5

ηOWC (%)

L/B=3

5

4

θ / (2ai /λ)

40

0.1

p / ρw g ai

d=3

D/ρ√ gB3

5

B=1
0.5

10

20

30

40

50
λ/B

60

２次元数値解析による浮体型波力発電装置特性の妥当性と最適形状の検討

する．表１はこの結果を表し，L/B=8 の形状が最も小さく，考察する．表２はこの結果を表しており，エネルギー吸作動範囲も広く，最適な形状となっている．L/B=8 の作

収効率が 50%以上を作動範囲としている．この結果，評動範囲が広くなる理由は，動揺によりλ /B=45 付近の波

2は L/B=1.5 の形状が最も小さく，作動
価関数 Ld / λMAX長にも共振が生じ，効率が向上するためである．

範囲λ M A X /λ M I N は L/B=6 が最も広くなる．L/B=1.2 は

図５は船首に空気室を有する浮体の特性を示してい

全域にわたってエネルギー吸収効率が 50%以下と低くる．水深 h が空気室奥行き B の 40 倍，喫水 d を１倍となり，短すぎる浮体は適さないことが分かる．エネルギーする浮体形状について，浮体長の影響を検討している．

吸収効率および作動範囲的には浮体長が長い L/B=6

浮体の全長 L/B は 1.2 ～6.0 倍までの 5 種類とし，波長が優れているが，浮体サイズが大きくなるため，経済的ごとに，負荷ダンピング係数を変化させて最大エネルギには浮体サイズが小さく，作動範囲も 2.63 倍と十分にー吸収効率を求め，その負荷ダンピング係数 D，浮体

Table 2. Typical effects of the length for the device with OWC at
front side, where is over 50% in the efficiencyと空気室内水面との相対変位振幅 zOWC，空気室内圧
力振幅 p，エネルギー吸収効率ηOWC，浮体重心位置の

L/B
1.2
1.5
2.0
4.0
6.0

上下変位振幅 ζ ，縦揺れ変位振幅 θ を示している．最大エネルギー吸収効率は全長 L/B=6 の浮体に生じ，波長と空気室奥行きの比はλ /B=6.2 付近に現れる．作
動範囲は浮体長が長くなるほど拡大するがこと分かる．
ここで，式(7), (8)の評価関数により最適浮体形状を

λ/B

L/λMAX

2
Ld / λMAX
4.0 ∼ 10.6
3.5 ∼ 9.3
3.1 ∼ 10.9
3.0 ∼ 13.0


0.142
0.215
0.367
0.462


0.01330.02310.03370.0355λMAX/λMIN

2.63
2.63
3.57
4.35

L

1.2

3

4

2
1

p / ρw g ai

2.0

6

0
1.5

L/B=1.2

L/B=1.5

2
4

1.0

1.5
1.0
0.5

L/B=1.56

0.5

L/B=1.54

ζ / aiL/B=2.0

L/B=4.0

2

3
2
1.2

1
0

L/B=6.0

0

2

4

6

4

θ / (2ai /λ)

0.0
5

L/B=6
4

2

80

L/B=6

4

2

60

1.5

40

1.2

20

1.2

L/B=1.5

1.2

1.5

0.0
100

ηOWC (%)

0.1

d=1

40

D/ρ√ gB3

0.2

2

4

zOWC / ai

5

B=1

6

8 10 12 14 16 18 20
λ/B

0
60
50
40
30
20
10
0

L/B=1.5

1.2

2
4

0

4

2

6

6

8 10 12 14 16 18 20
λ/B

Fig.5 Effect of the length of floating devices which have OWC at front side.
1.5

2
1

0.5

0.5
0

2

4

6

8 10 12 14 16 18 20
λ/B

0.6
0.4
0.2
0.0
5

θ / (2ai /λ)

1.0

0.0

0.8

1.0

0.0
100

ηOWC (%)

p / ρ w g ai

0
1.5

1.0

ζ / aizOWC / ai

D/ρ√ gB3

3

80
60
40

B

4
3
2

20

1

0

0

0

2

4

6

8 10 12 14 16 18 20
λ/B

Fig.6 Characteristics of Mighty whale type device.

－6－

0

2

4

6

8 10 12 14 16 18 20
λ/B２次元数値解析による浮体型波力発電装置特性の妥当性と最適形状の検討

広い L/B=1.5 が最適な形状である．浮体長の短い

動揺させてエネルギーを吸収する構造になるため，大きL/B=1.5 で作動範囲が広くなった理由は，動揺によりな動揺により設置機器への強度面での負担が大きくな

λ /B=10.2 付近の波長で共振が生じ，効率が向上するる点が挙げられる．

ためである．図２に示すように一般的には動揺によりエ

このように，より現実的な浮体形状の選定には，機器

ネルギー吸収効率が低下するが，動揺による共振が長の実装などを含めた構造，建造コストや入射波の出現

波長側で効果的働くと作動範囲が拡大し浮体サイズの率を考慮することが必要である．ただし，初期設計の段低減に寄与できることが分かる．

階では大まかな形状を決定する必要がある点や小型化

図６は代表的な浮体形状として，海洋科学技術センの可能性を示すことができるため，ここに示した簡易的ターにより実証実験がなされたマイティホエール( 9 ) を挙な手法は有効である．

げ，この浮体形状について２次元数値解析を行った結

5. 結

言

果である．この浮体形状は高いエネルギー吸収効率98%と広い作動範囲λ M A X /λ M I N =2.59 を得ることができる．
ただし，波長λ /B=8.6 付近に現れる共振によるエネル

本研究では空気室を有する波浪変換エネルギー装
置特性を数値解析により推定し，最適浮体形状について検討を行った．数値解析結果の妥当性を検証するた

ギー吸収効率の明確な寄与は確認できない．

め，水面上に浮体が固定された状態と無拘束に近い状
4. 最適浮体形状の検討

態で動揺するときの基礎特性を実験と比較し，良好な

マイティホエール型浮体形状との比較を通して最適推定精度が得られることを確認するとともに，浮体の動浮体形状について考察を行う．図７は 3 種類の浮体形

揺特性と空気室特性との関連を把握しつつ，最適浮体状，マイティホエール，空気室が船首および船尾にある形状の提案を行い，以下の結論を得た．

形状についてサイズ比較を行っている．浮体サイズに影( 1 ) 浮体の動揺特性と付随する空気室の特性は領響する最も厳しい条件は運転範囲の中で最も長波長の域分割法により取扱っている．解析は供試体の形状に波であり，この波長に対応した浮体長が最も長くなる．こ合わせ，計算領域を８つに分割して実施した．

のため，50%以上のエネルギー吸収効率となる波長を作( 2 ) 実験では空気室内の圧力，反射波，透過波を動範囲とし，運転すべき最大波長でのサイズを比較して計測し，空気室内の圧力から空気室効率（エネルギー

いる．最適浮体形状は(b)の空気室が船首に設置された吸収率）が算出されている．数値解析は実験と同じ項目

型式で，浮体長と空気室奥行きの比が L/B=1.5 の浮体を算出するとともに，浮体の動揺特性の結果も算出してである．同一波長の波に対して，浮体断面積で比較すいる．これらの数値解析結果は実験結果と比較され，実ると最適浮体形状はマイティホエール形状の浮体の験が行われた波長の範囲では計算結果と実験結果は1/6 程度のサイズに抑えられている．ただし，課題も考全ての実験点で良く一致する結果が得られた．実験デえられ，小型化とともにタービン発電機などの実装スペータのバラツキや誤差を考慮すると，本解析方法によるースや配置などの面で自由度が少なくなる点や浮体を結果は実験結果に匹敵する精度を期待できる．
( 3 ) 実験結果および計算結果から，空気室内の圧L / λ M A X = 0.602

力は波長 λ /B=8 付近が最大値となり，最大効率点は波L/B=5
d/B=0.8
L d / λ 2M A X = 0 . 0 5 8
λ MAX/ λ MIN= 2 . 5 9

d / λ M A X = 0.096

長が空気室奥行きの 6 ～ 8 倍，λ /B=6 ～ 8 付近に現れ，波の振幅 a i / B ≅ 0.063，開口比 R=1/150 では 80%と高い効率が得られている．この波長では反射波および透過(a) Mighty Whale type device
L / λ M A X = 0.11

d / λ M A X = 0.086

波の振幅が極小値をとり，一旦，小さな値となり効率の
L/B=1.5
d/B=1.2
L d / λ 2M A X = 0 . 0 0 9 2
λ MAX/ λ MIN= 3 . 5向上に寄与した後，これより大きな波長では反射波および透過波の振幅が急激に増大するため効率は急激に低下する．これに対して，反射波の波高は波長が短くな(b) Optimal device with OWC installed at the front side

るつれ増大し，効率の低下に強く影響を与えている．
L / λ M A X = 0.19
d / λ M A X = 0.024

L/B=8
d/B=2.5
L d / λ 2M A X = 0 . 0 1 1 3
λ MAX/ λ MIN= 1 . 9 5(c) Optimal device with OWC installed at the stern side( 4 ) 浮体の動揺による影響は最大効率点より長波長側で急激に現れ，透過波および反射波の波高が急増して効率が急激に低下する．これに対して浮体が固定状態のときは効率の変化が非常に緩やかである．最大効

Fig.7 Size Comparison with Mighty whale type device and optimal
profile devices

率点の位相は上下揺れ，縦揺れともに入射波に対して－7－２次元数値解析による浮体型波力発電装置特性の妥当性と最適形状の検討( 4 ) Suzuki, M. and Arakawa, C., Numerical Methods to
Predict Characteristics of Air Chamber for Fixed
Terminator Type of Wave Energy Converter, Transaction of the Japan Society of Mechanical Engineers ,Series B, Vol.62, No.604 (1996), pp.4135-4141.

90°遅れて，反射波高，透過波高ともに極小値を取り，特に反射波高はゼロに近い値となっている．最大効率
点よりも短波長側は動揺の影響は少なく，固定状態の特性と同じ傾向を波長に対して示している．これは入射( 5 ) Suzuki, M. and Arakawa, C., Numerical methods topredict characteristics of oscillating water column for
terminator type of wave energy, Proc 13th Int Offshore
and Polar Eng Conf, Honolulu, Hawaii, USA, ISOPE,
CD−ROM (2003), pp.333−340.

波と空気室内の水柱の動揺の関係が浮体動揺に比較して支配的であることを示している．
( 5 ) 最適浮体形状を評価するため，浮体サイズと作( 6 ) Suzuki, M., Washio, Y., and Kuboki, T., Analysing
Method of Floating Type Wave Energy Conversion
Device with Air Chamber, Transaction of the Japan
Society of Mechanical Engineers , Series B, Vol.72,
No.718 (2006), pp.1529-1536.

2
ηOWC とλ M A X /λ M I N な る 評
動範囲について， Ld / λMAX

価関数を導出した．
( 6 ) マイティホエール形状の浮体は高い空気室効率と広い作動範囲を有しており，最大効率は λ /B=6 付近( 7 ) Suzuki, M. Washio, Y., and Kuboki, T., Characteristics of Floating Type Wave Energy Conversion Device
with Air Chamber, Transaction of the Japan Society of
Mechanical Engineers, Series B, Vol.72, No.722(2006), pp.2488-2495.

に現れ，99%にも達している．しかし，ここに提案した最適浮体形状に比較して大きな浮体となる．
( 7 ) 提案した最適浮体形状は船首に空気室を有し，( 8 ) Maeda, H., Kinoshita, T., Masuda, K. and Kato, W.,
Fundamental research on oscillating water column
wave power absorbers, Journal of Energy Resource
Technology, ASME, Vol.107 (1985), pp.81−86.

L/B=1.5 ， d/B=1.2 で あ り ， サ イ ズ を 表 す 評 価 関 数
2Ld / λMAXがマイティホエール型，船尾空気室型に比して，最も小さな値，0.0092 となっている．

文( 9 ) Washio, Y., Osawa, H., Nagata, Y., Fujiii, F., Furuyama, H., and Fujita, T., The offshore floating typewave power device “Mighty Whale”: open sea tests,

献Proc 10th Int Offshore and Polar Eng Conf, Seattle,
ISOPE, Vol.1 (2000), pp.373−380.

( 1 ) Raghunathan, S., The wells air turbine for wave energy conversion, Prog. Aerospace Sci., Vol. 31 (1995),
pp.335 −386.

(10)( 2 ) Suzuki, M. et al., Analysis and Design on Guide
Vanes Characteristics of Wells Turbine for Wave
Power Generator, Transaction of the Japan Society of
Mechanical Engineers, Series B, Vol.65, No.636(1999), pp.2718-2725.

Osawa, H., Nagata, Y., Miyajima, S., and Maeda, H.,
A Design of Mooring System for Floating Wave Energy Converter in Shallow Water，Journal of the Society of Naval Architects of Japan ， Vol.182 (1997) ，pp.341−348.

(11)Suzuki, M. and Arakawa, C., Method of Selection on
Wells Turbine Profile for Wave Power Generating
System, Transaction of the Japan Society of Mechanical Engineers, Series B, Vol.69, No.683 (2003),
pp.1621-1627.

( 3 ) Suzuki, M., Washio, Y., and Kuboki, A., Characteristics of Floating Type Wave Power Generating Systemwith Oscillating Water Column, Proc 15th Int Offshore
and Polar Eng Conf, Seoul, Korea, ISOPE, CD-ROM(2005), pp.529−536.

－8－