重点課題⑥：革新的クリーンエネルギーシステムの実用化～高効率風力発電システム構築のための大規模数値解析～
東日本大震災後のエネルギー問題に対する国民の意識*1
・低炭素社会の実現を目標とした原子力に依存しない安全で持続的な社会の構築
・今後増やすべき発電用エネルギー源：自然エネルギー 70% （NHK世論調査）*1 独立行政法人 経済産業研究所：http://www.rieti.go.jp/jp/projects/research_activity/energy/→日本沿岸の潜在的エネルギー源を活用した大規模洋上ウィンドファームの開発が不可欠

風車の大型化（直径200m ：20MW)

洋上風車（2015年）

日本経済新聞電子版2013年7月18日

三菱重工プレスリリース 第 5427号

大規模洋上ウィンドファーム
K.S. Hansen, EERA DeepWind’2014, 2014

解決すべき課題：・ウィンドファームでは下流にある風車の発電量低下，翼の疲労破壊が問題
・気象変動の影響を考慮した発電量予測が困難
・直径（200m）の大規模風車の試作実験や発電量を実験だけで予測するのは困難→大規模数値解析によりウィンドファームの発電量の予測や寿命解析を行い，自然エネルギーを安定に利用する技術を確立

ポスト京が拓く「ウィンドファームの高効率化」
ポスト京を活用した1兆規模の流体構造連成解析によりウィンドファームの発電量を予測
・流れ場の大きさと時間規模：気象（10km ：年)，構造物(100m：分)，翼周り（0.1mm/0.001秒）
・考慮すべき物理現象：気流変動，翼の疲労破壊，流体制御，発電量制御

ウィンドファーム発電量：400万kw（2030年）実現への設計支援
洋上ウィンドファーム発電コストの改善（現状：陸上風車の2倍→ポスト京：36円/kWh以下）（１）最適配置による性能向上，稼働率改善，（２）先端流体制御による性能改善（風車単体10%)（３）ブレードの寿命評価によるメンテナンスコスト削減

京時代（現在）ポスト京時代（将来）

風車単体の大規模な非定常解析が可能となり、さらに空力・構造連成解析による風車の構造解析も可能となった。風車の後流の影響や地形の影響を考慮した非定常解析が可能となった。しかし、ウィンドファームのような風車群の流れ解析最適設計のための解析は実現できていない。

大気境界層、地形を考慮した大規
模ウィンドファームにおける風車間の流れの相互干渉を解析することにより、発電量の向上、ブレードの寿命改善、低コスト化が可能となる。
これらの技術によって環境に優しい安定した電力供給を実現すると共に我が国の風力発電産業の競争力向上に貢献する。

地形を考慮した風車解析

大規模解析（京）による最適化の事例
ポスト京では風車配置の最適化に適用予定

・解析アプりケーション
FrontFlow/blue (大規模流体解析）
Adventure (大規模構造振動解析）
LANS3D（先進流体制御解析）
RIAM-COMPACT (数値風況解析)
REVOCAP （マルチ力学解析）

先進流体制御による発電量の改善

アプリの準備状況（京時代）
FFB

LANS3D• 大規模乱流解析による流体機械・自動車等の高精度予測• 先進流体制御による翼性能の改善，風車性能改善技術の提案- 自動車・船舶・流体機械等において100億メッシュ規模の大規模乱流解析を実現（産業利用における実績）- 気流の乱れを考慮した翼周りの流れ解析• 高精度流体解析による翼の最適制御方法の検討（数百ケース）• 流れの制御による風車の性能改善（10%)

翼周りの流れの最適化制御
流体機械解析

主流乱れを考慮した翼周りの流れ解析• FFB-FISTARによる風車ブレードの応力解析• FFB-ADVENTUREによる自動車車内音連成解
析• ・流れ場と構造振動・音までの連成解析システムを構築し，産業界での実証解析を実現

風車応力解析（三菱）

はく離制御による風車の性能改善

自動車車内解析（スズキ）

RIAM-Compact• 複雑地形を考慮した風車性能の予測• 複雑地形によって生じる流れ場を予測し，風力発電の設置に適切かどうかの判定が可能- 風況マップの作成による風車配置の最適化

大気流れが風車に及ぼす影響

風況調査結果と風車配置

ポスト京で開発するアプリの概要
WindFarm

LANS3D・FFB-ADVENTURE, RIAM-CPMPACT，
LANS3Dの特色を活かした大規模空力連成システ
ムによるマルチスケール解析ソフトの開発• プラズマアクチュエータ等の先進流体制御技術を組み込んだ大規模乱流解析• 高速ソルバーを利用した最適化解析
統合システム化

流れ場の制御技術の高度化

ウィンドファームの流れを再現

PAによる風車周りの流れの制御

統合システム化

FFB，FFB-ADVENTURE

統合システム化・FFBによる1000億規模大規模乱流解析
・FFB-ADVENTUREによる空力構造連成解析（寿命評価）

１０MW(２００ｍ）急須亜の解析（図は小型風車の解析事例）

RIAM-Compact (HPC)• 洋上流れ解析及び洋上風況マップの作成• 風車群解析のためのモデル化及び最適化

構想解析による
応力評価
湾内の流れ解析：風況マップ

複雑地形における風車周りの流れ場の評価

解析事例(NREL Offshore 5-MW Baseline Wind Turbine)

1億7000万要素 (富士通FX10 180時間 )5風車周りの流れ6

タンデム解析

Rotor Diameter
Tip Speed
Time period of one rotation
Distance between turbines
Uniform velocity
Elapsed time to pass between turbines
Ratio of Elapsed time to time period
Time step
Number of Time steps(one time unit)
Number of steps for one rotaion
Elapse time of one step
Nuber of Nodes
Number of Elements
Elapsed time of Computation(4 cycles of one time unit )D [m]V [m/s]T [sec]L [m]U [m/s]t [sec]t/T [-]
dt [sec]

L/D=0.0
L/D=0.5
1261126
179.16813
1 79.16813
15 3.1415927
5 3.1415927
0063
0.5
10.68770.135 10.6877
0.135
5.894628 3.7037037
1.178926 1.1789255
0.002441
0.001534 0.002441
0.001534

NTIME [Step]

2048

02415

2414.4394

NTIME*4[Step]
N [Step]
Te [sec]
MDOM
ME [million]

10000

10000
2048
30
96
80

10000

10000
2048
30
96
150

Time[hour]8484

L/D=1.0
126
179.16813
1
5 3.141593
126
110.6877
0.135
11.78926 7.407407
2.357851 2.357851
0.002441 0.001534

L/D=2.0
126
179.16813
1
5 3.141593
252
210.6877
0.135
23.57851 14.81481
4.715702 4.715702
0.002441 0.001534

L/D=3.0
126
179.16813
15 3.141593
378
310.6877
0.13535.36777 22.22222
7.073553 7.073553
0.002441 0.001534

4829 4828.879

9658 9657.758

14487 14486.64

24145 24144.39

20000

40000

60000

100000

187

20000
2048
33.6
96150

200

40000
2048
18
240
150

334

60000
2048
20480
300

L/D=51261
79.16813
1
5 3.141593
630
510.6877
0.135
58.94628 37.03704
11.78926 11.78926
0.002441 0.001534

100000
2048
20864
600

556

7

後流による干渉
1D

12 m

上流側の風車
風車間隔1Dにおける下流側の風車
風車間隔2Dにおける下流側の風車

1Dよりも2Dの方が流入速度の損失が少な
いyz1D と 2D の 間 で ，流入速度が最も小さくなると予想される

2D
2Dy

風車12m手前の平均風速分布z

風車10回転目における流れ場

下流側の風車に流入する前に平均風速が減少
8

～RIAM-COMPACTを用いた分析
開発目標・開発項目
・アクチュエータディスク，アクチュエータラインモデルを用いた風車ウエイクの再現
・大規模計算のための風車ウエイクモデルの開発と実装

研究成果と課題：・風車単体における風車ウエイクの再現に成功
・大型風車への適用と簡易ウエイクモデルの開発と実装
・実機風車による実スケールデータ(実測データ)による検証
2016/8/17

重点課題⑥-C 洋上風車

9

大気の乱れの影響評価
開発目標・開発項目
・ポスト「京」で実行可能な計算コードの改良・大規模ウィンドファームに適用可能な簡易風車ウエイクモデルの開発と実装
・流入変動気流の作成と検証
大u

小zzyy
時間t

時間t

LESのためお流入変動気流の作成例：
主流方向(x)の速度成分の分布
2016/8/17

重点課題⑥-C 洋上風車

10