数値風況診断技術RIAM-COMPACTの開発と風車に対する地形性乱流の影響に関する研究-コンピュータシミュレーションが解決する風力発電の諸問題-

九州大学・応用力学研究所・風工学分野
内田 孝紀
takanori@riam.kyushu-u.ac.jp092-583-7776

非定常・非線形数値風況シミュレータ
RIAM-COMPACT®（リアムコンパクト）

九州大学応用力学研究所発の流体工学CFDモデル(LES乱流モデルを採用
LES乱流モデルを採用))
Research Institute for Applied Mechanics, Kyushu University,
COMputational Prediction of Airflow over Complex Terrain◆乱流モデル：次世代乱流モデルとして期待されているLES(Large-Eddy Simulation)の採用！
◆対象スケール：数(十)km以下の局所域スケールに的を絞る！
◆特長①：風に対する地形・建物の効果を高精度に予測！
◆特長②：安定成層，不安定成層など種々の大気安定度を考慮！
流れの局所的な増速
Flow

逆流域の形成(渦放出)

流れの衝突
流れの剥離流れの再付着【風に対する地形効果】

RIAM-COMPACT®の理論および解析手法■ 取り扱う流体 ： 非圧縮性・粘性流体■ 基礎方程式 ： （LESよりFilterを施した）連続の式，Navier-Stokes方程式■ 離散化手法 ： （有限）差分法■ 座標系 ： σ座標系（右図を参照）■ 変数配置 ： コロケート配置■ 空間項の離散化 ： 対流項を除いて2次精度中心差分■ 対流項の離散化 ： 3次精度風上差分（数値粘性項の係数0.5）■ 乱流モデル ： Large-Eddy Simulation■ SGSモデル ： 標準Smagorinskyモデル■ 圧力のPoisson方程式 ： SOR法
■ 地面境界条件 ： 非滑り条件

RIAM-COMPACT®(リアムコンパクト)の着眼点
流体(風況)のコンピュータシミュレーションにおける空間スケールによる分類＜(a)工学分野では数m以内の物体スケール＜
地上10mにおける水平風速
地上10mにおける水平風速と風速ベクトル，0418
台風と風速ベクトル，0418台風(b)数m～数十kmの局所域ｽｹｰﾙ
数時間程度の時間ｽｹｰﾙ(c)気象分野では数十km以上の広域空間スケール

工学分野では数m以内の物体スケール、気象分野では数十km以上の広域空間スケールが主な研究対象であった。
一方、その中間の数m～数十kmは私達の実生活に関連したごく身近な風環境にも関わらず、地物・地形による影響と時間変化の両者を忠実に再現可能な風況予測モデルの重要性は、未だ十分に認識されてなかった。
我々は、このスケールギャップに着目し、世界最先端の流体計算技術に基づいた斬新で汎用的な風況予測システム(RIAM-COMPACT®：リアムコンパクト)の開発に成功し、この分野の新たな道を拓いた。

RIAM-COMPACT®の応用分野
・山岳地形・市街地における自然エネルギーの有効利用(風力発電)
・巨大都市域(市街地)の大気環境予測と改善
・山間部の地形性局地強風の発生メカニズム解明
・山間部の送電鉄塔周辺の風害対策
・台風に伴う歴史的建造物の風害対策
・竜巻に伴う風害対策【風工学分野】

特に重点を置いているテーマ【鉄道分野】
・突風・強風時における鉄道の安全運行支援システム構築
・線路周辺の風害対策，風況マップ作成【航空・船舶分野】

竜巻シミュレータの開発

・離島空港建設のための風況アセスメント
・大型タンカー接岸ルート支援システムのための風況予測
・台風に伴う風害対策のためのハザードマップ作成(強風域特定)・山火事の延焼域の予測
・山火事の煙，火山ガス，大気汚染物質，花粉などの移流・拡散予測【森林分野】

風向変動を考慮した
ガス拡散予測【レジャー分野】
・ヨットレース，フィッシング，ゴルフ，バルーンなどを対象にした風情報配信サービス

福岡市の大規模風環境予測：福岡市からの委託PJ
福岡市役所 「都市計画基礎調査」のGISデータを利用

北北西の風330m

4,000m
3,800m◆計算領域：4,000m×3,800m×330m◆格子数：1,001×951×81点(約7,700万点)◆格子解像度：
水平方向 4m，鉛直方向 0.5m
＊上空ほど間隔を広げる◆代表長さ：50m

地上2mにおける主流方向速度成分の分布＋速度ベクトルの分布

北北西の風
天神地区を拡大

アクロス福岡最近の受賞

2015年 ISIT創立20周年記念
九州先端科学技術研究開発表彰

内田その他の受賞歴◆ NEDOの産業技術研究助成事業に採択．風車に対する地形効果に起因する風の乱れの影響を解明．この成果に対し，2010年 科学技術分野の「文部科学大臣表彰･若手科学者賞」を受賞◆ 2013年 内閣府主催の「第11回産学官連携功労者表彰・環境大臣賞」を受賞

本研究における計算機リソース～2015年1月

スーパーコンピュータNEC 「SX-ACE」 16ノード（最大理論性能4.1TFLOPS）
2015年3月より稼働開始

風車性能評価のためのCFD
Flow-レンズ風車の性能予測-

アクチュエータラインモデルを使用

スパン中央断面における主流方向速度成分の分布

計算条件など①一様流
3D (奥行方向も3D)
原点x=y=z=0
原点x=y=z=0

1.5D

D

発電量予測とウエイクの挙動の把握を目的とする！
10D■MX=741■MY=401■MZ=401
合計：約1億2千万点
Δxmin=Δymin=Δzmin=0.005D(注意：実際の計算では0.5Dが代表スケール)

計算機環境
・応用力学研究所計算機室所有のスパコン SXSX-ACE （NEC製）
NEC製）
・Memory 約26GBを使用
26GBを使用・ファイルサイズ 約7.5GBを使用
7.5GBを使用

計算条件など②一様流

Cii diffuser : つば10％

メモ

矩形格子近似法を用いるため，若干厚めにモデル化

ADM(アクチュエータディスクモデル)
アクチュエータディスクモデル)の範囲：■i=199～203, 厚さ0.02D■j=101～301■k=101～301
周速比：λ=4.0（最適値）

風車模型を使用した風洞実験
zyxy

MEL翼をベースに翼根を厚翼に改良

気流計測を実施x

ロータ直径D=1(m)配慮点◎翼のコード長に対するレイノルズ数の低下を避けるため，ロータ直径D=1(m)の風車模型を使用．
ロータ直径D=1(m)の風車模型を使用．
◎ブロッケージの影響を避けるため，側壁と上壁は外して実験を実施．

Flow

2D4D

Near wake

6D

Uref=8(m/s)

8D

10D

Far wake

ロータ直径に基づいたレイノルズ数Re(=UrefD/ν)≒5×105

◆Cw = P/0.5ρAU3
P is output power
A is swept area of turbine blades
U is approaching wind speed

計算結果①Power coefficients Cw◆λ = ωr/Uω is the angular frequency, 2πnr is the radius of a wind turbine rotor

Lens turbine 0.72（数値計算）
約1.9倍の性能向上

Bare turbine 0.38(数値計算)

Tip-speed ratio λ

広範囲のλに対するCwのプロット

計算結果②Flow
アクチュエータディスクモデルを使用

注目注目

スパン中央断面における主流方向速度成分の分布

計算結果③Flowアクチュエータディスクモデルを使用

注目

注目

スパン中央断面における主流方向速度成分の分布

Flow

計算結果④

Gap Flowの形成

計算格子図とスパン中央断面における主流方向速度成分の分布

航空レーザ測量を取り入れた
高解像度数値風況シミュレーション【克服すべき主な課題】
如何に精緻に（数値シミュレーションの1)地形の凹凸境界条件として）再現するかが極めて重要！
2)地表面粗度◆航空レーザ計測技術の有効利用
・地理情報システム(GIS)との連携を図り，独自の詳細地形データ構築技術を確立．
・最新の土地造成（盛土や切戸など）の状況や地表面粗度の再現

引用先 http://www.ajiko.co.jp/

・現況の敷地造成・搬入道路状況を設計CADデータに基づいて修正可能

70kWレンズ風車：伊都キャンパス
12.8m
15.44m

25.0m

No.1
No.2

【風車周辺における風環境の向上を目的とした数値風況予測の実施】

航空レーザ測量の実施：樹木の空間分布と樹木高さの把握

レンズ風車
樹木高さ：
平均10.8m

樹木高さ(m)=数値表層モデル（DSM; Digital Surface Model）－数値標高モデル（DEM; Digital Elevation Model)

伊都キャンパス周辺の風況特性(冬場：2015.10-12)

WNW

①②5km

③5km

④計算領域および地表面粗度のモデル化の範囲など10mの樹木高さとして再現

風向WNW(292.5度)

レンズ風車

モデル化のイメージ図

風向WNW(292.5度)

レンズ風車
No.1400m

No.2100m

計算結果①

鉛直断面内の速度ベクトル図と風速分布図（時間平均場）

地表面粗度無し

WNW

地表面粗度
有り
地表面粗度の影響

計算結果②
風向WNW(292.5度)

地表面粗度無し

Bird’ eye view

Rear view
流れは向こうからこちらへ

計算結果③

風車ハブ高さ(地上高25m)における主流方向風速の分布(Google Earthとの連携技術の開発)