倉田奨励金研究報告書｜2018 年度（第 50 回）助成■受領№1363

高精度数値流体シミュレーションと最適化手法によるウィンドファームの風車配置の最適化

代表研究者

Goit Jay Prakash

近畿大学工学部機械工学科

1．研究⽬的助教

2．研究内容近年多数の⾵⾞を⼀ヶ所に設置し、集中的に発電する⼤型陸上・洋上ウィンドファームが増加してい上記の⽬的にしたがい、本助成期間に以下 3 点に着⽬し研究を実施した。

る。世界最⼤級の洋上⾵⼒発電所はイギリスのLondon Array (175基の⾵⾞)とWalney Extension(87基の⼤型⾵⾞)ウィンドファームで、それぞれ630MWと659MWの発電を⾏う。⽇本でも2019年に政府が洋上⾵⼒発電に適した11の「促進区域」を指定したことから⼤型洋上⾵⼒発電の開発が進んでいる。しかし、このような⼤型ウィンドファームでは、⾵⾞間の後流の相互作⽤により発電出⼒が低下し、単独の⾵⾞に⽐べて出⼒が40%まで低くなるケースもある(Barthelmieら, J. Phys: Conf. Ser. 2007)。ウィンドファーム内の⾵⾞の配置を決めるとき、⼀般的には簡易的なウェークモデル(Jensen,Ris-M-2411 1983) が使⽤される。しかし、このようなウェークモデルではウィンドファーム内の複雑な流れを考慮することが出来ないため、効果が⼤いに問われる。従って、ウェークの相互作⽤を考慮し、ウィンドファームの⾵⾞配置の改善により、ファーム全体の発電出⼒の向上を⽬的とし、以下3点に着⽬した。
(1) ドップラーライダーによるウィンドファーム開発サイトでの⾵況計測技術の開発(2) ウィンドファーム内の⾵⾞の再現と⾵の計算⽤の数値流体シミュレーションツールの開発(3) ウィンドファーム内の⾵⾞配置改善による発電出⼒の向上

2.1 ドップラーライダーによるウィンドファー
ム開発サイトでの⾵況計測技術の開発
研究に必要な実測データはドップラーライダーを⽤いて収集する予定であった。ドップラーライダー
は、⾼⾼度および幅広い範囲での⾵速の計測が可能であるという利点があるものの、従来の超⾳波⾵速計やカップ⾵速計と同程度計測精度が得られるかは明確ではなかった。本研究では、鉛直ライダー(Windcube V2)により計測を⾏い、隣接する観測タワーに設置した⾵速・⾵向計と検証を⾏った。ライ
ダーと観測タワーは図1 (a) と(b)に⽰す、産総研福島再⽣可能エネルギー研究所のものを使⽤した。ドップラーライダーの計測から⾵速を算出するためにDBS (Doppler Beam Swinging) 法を⽤いた。
図1(c)はドップラーライダーによる計測した10分間の平均⾵速と標準偏差を⾼度 57 m に設置された超⾳波⾵速計の計測値との⽐較を⽰す。ドップラーライダーの平均⾵速は超⾳波⾵速計の平均⾵速と精度良く⼀致しており、決定係数 (R2)は0.995と回帰線の傾きは1.0である。しかし、標準偏差の場合バラツキが⼤きく、決定係数は0.977である。図1(d)は乱流強度を⾵速の関数として⽰している。計測したサイトでは、ドップラーライダーから得られた乱流強度は超⾳波⾵速計の乱流強度に⽐べて約2%⾼倉田奨励金研究報告書｜2018 年度（第 50 回）助成

（a）観測タワー（b）鉛直ドップラーライダーと DBS法の概要（c）鉛直ドップラーライダーと超⾳波⾵速計の 10 分平均⾵速及び標
準偏差の⽐較
（d）乱流強度と平均⾵速の関係
図 1. 鉛直ドップラーライダーによる⾵況計測精度の評価
い結果が得られた。この差はあんまり⼤きくない全体をシミュレーションする必要があった。流体シ

ものの、本サイトに設置する⾵⾞の疲労荷重の予ミュレーションはラージエディシミュレーション

測に影響与えると思われる。さらに、計測した乱(LES)を⽤いて⾏った。その⽀配⽅程式はフィル

流強度分布は IEC (International Electrotechnical

ターされたナビエストークス⽅程式である：Commission) 規格に定義されている全ての乱流カ𝜕𝑢𝜕𝑥

テゴリーより⾼いため、サイトの乱流強度を考慮した ʻclass Sʼ ⾵⾞が必要となる。⽇本及びアジア諸国の場合、乱流強度が⾼い地域が多く存在する。ここでは⽰していないが、⾵⾞に作⽤する極値荷重評価に必要なピーク⾵速の⽐較、研究に良く⽤いられる
NREL-5MW⾵⾞を対象にドップラーライダーと超⾳波⾵速計を⽤いた場合の発電量分布と荷重を評価した。
2.2 数値流体シミュレーションツールの開発
流体シミュレーションではOpenFOAMというオープンソースツールを使⽤した。対象とした⾵⾞はロータ直径100mの⼤型⾵⾞であるため、⼤気境界層

𝜕𝑢
𝜕𝑡

𝜕
𝑢𝑢
𝜕𝑥
𝐹
𝑓1 𝜕𝑝
𝜌 𝜕𝑥

1
𝐶 𝜌𝑢 𝐴21
𝐶 𝜌𝑢2𝑑𝑥

ここで 𝑢0(1)
𝜈𝜕 𝑢
𝜕𝑥

𝜕𝜏
𝜕𝑥

𝑓(2)

(3)
(4)

𝑢 , 𝑢 , 𝑢 速度、𝑝 は圧⼒、𝜌は空気密度、𝜈は動粘性係数、 𝜏 はサブグリッドスケール(SGS) モデルである。𝑓 は⾵⾞が流れに与える⼒を⽰す。後で説明するように、本研究ではこの⼒は
actuator-disk model (ADM)によりモデル化した。

倉田奨励金研究報告書｜2018 年度（第 50 回）助成
⼤気境界を数値計算で正確に再現する際、精度の⾼い流⼊境界条件が課題である。本研究では流⼊⾵を⽣成するために⾵⾞無しの領域で主流⽅向に周期境界条件を設けて計算を⾏った。そして、スパン̶鉛直⽅向の断⾯の各時間ステップの⾵速データを保存した。実際の⾵⾞やウィンドファームの流体計算

2.3 ウィンドファーム内の⾵⾞配置改善による発電出⼒の向上
開発した数値流体シミュレーションツールを⽤いてウィンドファーム内の流体解析を⾏った。本研究では、2 種類配置（格⼦型とスタッガード型）のウィンドファームを対象とした。両ウィンドフ
ァームには直径 D=100m の 16 台の⾵⾞がある。

時に保存していた断⾯のデータを流⼊⾵の境界条件

⾵⾞間距離は流れ⽅向にとして使⽤した。に 𝑺𝒛

𝑺𝒙𝟕𝑫、スパン⽅向𝟒. 𝟓𝑫とした。スタッガード配置の場合、各⾵⾞は上流と下流の⾵⾞よりスパン⽅向に2.25D ずらして設置した。計算領域は流れ⽅向、鉛直⽅向、スパン⽅向にそれぞれ 5.1 km, 1 km, 2km をとり、格⼦数は 850, 200 と 400 にした。合計格⼦点数は 68,000,000 と⼤型計算になったため、東京⼤学が運営する Oakforest-PACS スーパーコ
ンピューターを使⽤した。最⼤で 544 のプロセッサを使⽤し、並列計算を⾏った。
流⼊⾵に関しては，洋上ウィンドファームの建設が予定されているサイトの年間平均⾵速(7.5m/s)になるように流⼊⾵を⽣成した。図 3 はウィ

図 2. Actuator disk モデルの概要

⾵⾞周りの流れのシミュレーションでは、⾵⾞の適切なモデリングが⼤きな課題となる。しかし、⽻根上に細かい格⼦を⽣成し解析するのは実⽤的ではない。そこで、本研究では⾵⾞のロータのモデリン

ンドファーム内の⾵⾞ハブ⾼さでの⽔平⾯の平均⾵速を⽰す。ここで⿊い縦線が⾵⾞の位置を⽰す。
格⼦型ウィンドファームでは、列毎に⾵速の低い後流が出来ていることが観察される。その結果下流にある⾵⾞位置での⾵速がウィンドファーム上グをactuator diskモデル (図2を参照)を⽤いて⾏っ流の⾵速に⽐べ⼤幅に低くなる。⼀⽅、スタッガた。Actuator disk モデルの⼒は式 (3) のように定ードでは流れ⽅向の⾵⾞間距離が 2 倍になり、上義し、式 (2) の単位質量当たりの⼒𝑓 は式 (4)から流の⾵⾞の後流がある程度回復し、下流の⾵⾞の求めた。ここで、𝐶 は推⼒係数，𝑢 は⾵⾞上流側に位置での⾵速が⽐較的に⾼い。

おけるロータと同じ⾯積での平均⾵速、𝐴はロータ⾯積である。

（a）格子型配置(Rectangular)（b）スタッガード配置（Staggered）

図 3. ウィンドファーム内の平均⾵速場．ハブ⾼さでの⽔平⾯を⽰している倉田奨励金研究報告書｜2018 年度（第 50 回）助成
from LiDAR-measured wind speeds”, 41st
Wind Energy Symposium of JWEA, Tokyo,December 2019.

図 4. 時間平均及び横列平均発電出⼒の⽐較
この 2 種類の配置の発電出⼒の時間平均及び横列平均を⾵⾞列の関数として図 4 に⽰す。結果は横⼀列⽬の平均発電出⼒で正規化している。格⼦型配置の場合、下流の列では発電出⼒が急に低下している。⼀⽅、格⼦型配置に⽐べスタッガード配置の場合、下流の列の発電出⼒は⾼い。今回の計算のウィンドファームの⼤きさでは、配置改善による発電出⼒向上率は 51%である。したがって、ウィンドファームの配置の最適化は発電出⼒の向上につながると思われる。しかし、ウィンドファームが⼤きくなればなるほど、最適化の効果が低下する。そして、最適な配置を算出することも難しくなると思われる。

３．発表（研究成果の発表）(1) Goit, J. P., Shimada, S., Kogaki, T., “CanLiDARs Replace Meteorological Masts inWind Energy?”, Energies, Vol. 12 (3680), pp.
1-24, September 2019.
(2) Jay Prakash Goit，嶋⽥進，⼩垣哲也,“ドップラーライダーによる⻑期間の⾵速及び乱流計測の評価”, ⽇本機械学会 第97期 流体⼯学部
⾨ 講演会 豊橋

2019年11⽉(3) Goit, J. P., Shimada, S., Kogaki, T., “Reliabilityof long-term Lidar-based wind measurementsfor various wind energy applications”, 72nd
Annual Meeting of the APS DFD, Seattle,WA, USA, November 2019(4) Goit, J. P., Shimada, S., Kogaki, T., “Accuracyof wind turbine power and loads estimation