Direct Numerical Simulation of Vertical-Axis Wind Turbine Near-Wake Dynamics

本研究は幾何学的に解像された直接数値シミュレーションを用いて、垂直軸風力発電機のブレード数を増やすことがブレード・渦間相互作用を通じて動的ストール渦の崩壊を加速し、近接後流をより迅速に鈍体力学へ遷移させることを明らかにし、この遷移および下流流入特性を支配する主要因はタービンの先端速度比ではなくブレード数であることを示している。

要約

垂直軸風力タービン（VAWT）を巨大な機械ではなく、空の川の中に直立して回転する風車として想像してみてください。この論文は、ブレードの周囲を巻き起こす目に見えない「気流」がくっきりと見えるほど、ブレードに極めて接近してズームインする、高解像度のスローモーション映画のようなものです。研究者たちはスーパーコンピュータを用いて、回転するブレードの直後にある「近接後流」と呼ばれる領域で空気がどのように振る舞うかを正確にシミュレーションしました。

彼らが発見したことを、シンプルに解説します。

全体像：なぜ「近接後流」を見るのか？

ほとんどの風力タービンは飛行機のプロペラのように水平軸で回転しますが、これらの垂直軸タービンはサラダスピナーのように回転します。研究者たちは知りたいと思いました：ブレードを通過した直後の空気には何が起こるのでしょうか？

彼らは、空気がタービンの後ろを単に滑らかに流れるわけではないことを発見しました。むしろ、空気は混沌とします。ブレードが空気を切り裂き、**動的ストール渦（DSV）**と呼ばれる巨大で渦巻く竜巻を作り出します。これらは、ブレードによって蹴り上げられ、タービンの後ろに尾を引く巨大で目に見えない渦のようなものです。

主な発見：「ブレード数」が速度よりも重要

チームは1 枚、2 枚、そして3 枚のブレードを持つタービンをテストしました。また、タービンの回転速度も変化させました。

ここには驚くべきひねりがあります：タービンの回転速度よりも、タービンが持つブレードの数の方がはるかに重要なのです。

• 1 枚ブレードタービン（ソロの奏者）：
  部屋で一人踊るダンサーを想像してください。彼らが回転すると、背後に一つ巨大で強力な空気の渦が生まれます。この巨大な渦（DSV）は、長く強さを保ち、その姿を維持します。それは、消散するのに長い時間を要する重く、ゆっくりとした移動をする雲のようです。この巨大な雲が長く留まるため、タービンの背後の空気は長い距離にわたって「乱雑」で混沌とした状態のままです。

• 3 枚ブレードタービン（トリオ）：
  次に、互いに近くで回転する三人のダンサーを想像してください。一人のダンサーが回転すると、彼らが直前に乱した空気（渦）は、次のダンサーによって即座に衝突されます。
  研究者たちはこれを**「渦の衝突（Vortex Impingement）」**と呼ぶ新しいメカニズムを発見しました。

  • アナロジー： 巨大な渦を石鹸の泡と想像してください。1 枚ブレードタービンでは、その泡は丸ごと遠くへ浮かび去ります。一方、3 枚ブレードタービンでは、次のブレードがピンとして働き、泡を早急に割ってしまいます。
  • 結果： 巨大で乱雑な渦は、遠くへ移動する前に、小さくて無害な泡（より小さな渦）に砕き散らされます。3 枚ブレードタービンの背後の空気は、1 枚ブレードのものよりもはるかに早く「静か」で秩序だった状態になります。

「交通渋滞」効果

この論文はまた、ブレード数が多いことが風にとってある種の「交通渋滞」を生み出すと説明しています。

• ブレードが多いと、タービンの中心を風が通るための開放空間が少なくなります。
• これにより、風は川の中の岩の周りを流れる水のように、タービンの周りを通らざるを得なくなります。
• これにより後流の振る舞いが変化します。低いブレード数で起こる巨大な回転渦に支配されるのではなく、後流は単純な固体の岩（「非流線型物体」）の後ろにできる後流のように振る舞い始めます。
• これが良い理由： 「岩のような」後流は、「渦」の後流よりもはるかに早く回復（再び滑らかな空気になる）します。

「自己相似性」テスト

研究者たちは、タービンの背後にある乱れた空気が、いつ滑らかで予測可能な空気へと戻るのかを正確に知りたがりました。彼らは「自己相似性解析」と呼ばれる数学的なトリックを用いました。

• アナロジー： 煙の柱の形状を予測しようとしていると想像してください。最初は、煙は渦巻きの混沌とした塊です。しかし、最終的には予測可能なベルカーブの形状に落ち着きます。
• 彼らは、3 枚ブレードタービンの方が、1 枚ブレードタービンよりもはるかに早くこの予測可能な形状に落ち着くことを発見しました。「乱雑さ」はより早く消滅します。

未来への意味（論文によると）

この論文は特に、これらの発見がタービンを非常に接近させて配置する風力発電所にとって重要であると述べています。

• もし第一のタービンの真後ろに第二のタービンを置くと、第二のタービンが「吸い込む」空気は、第一のタービンが持つブレードの数に大きく依存します。
• 第一のタービンがブレード数少ない場合、第二のタービンは巨大で混沌とした渦に襲われます。
• 第一のタービンがブレード数多い場合、空気は第二のタービンに到達する頃にはすでに「回復」し、より滑らかになっています。

要約： この論文は、超高性能なコンピュータシミュレーションを用いて、垂直軸風力タービンにブレードを追加することが「渦の破壊者」として機能することを示しました。それは、巨大で乱雑な空気渦を小さな破片に砕き散らし、風がより早く滑らかな流れを取り戻すことを可能にします。これは、タービンを密集させて配置する風力発電所の設計において極めて重要です。

技術概要

技術サマリー：垂直軸風力発電機の近傍後流ダイナミクスに関する直接数値シミュレーション

問題提起
水平軸風力発電機（HAWT）が現在の研究を支配している一方で、垂直軸風力発電機（VAWT）は、重心が低いことや全方位への感度が低いことなど、洋上、都市部、遠隔地における応用に際して明確な利点を有する。VAWT の理論的な主要な利点の一つは、短い後流領域と強化された後流回復により、密に配置されたアレイにおいて高い電力密度を実現する可能性である。しかし、VAWT の近傍後流は、動的ストールに伴う複雑な流れの特徴、特に大規模な動的ストール渦（DSV）の形成と放出によって支配されている。遠方後流が鈍体後流のように振る舞うのに対し、近傍後流は、 tip-speed ratio（$\lambda$）や静的 solidity（$\sigma$）といった幾何学的および運転パラメータに強く依存する。既存の実験データは、幾何形状近傍または三次元的な解像度が不足していることが多く、アクチュエータライン法（ALM）と結合された高忠実度の大渦シミュレーション（LES）は、ブレード境界層を幾何学的に解像できていない。その結果、ALM は層流分離バブル（LSB）のダイナミクス、DSV の正確な形成、あるいはこれらの特徴を調節する重要なブレード - 渦相互作用（BVI）を捉えることができない。

手法
これらの限界に対処するため、著者らは、スペクトル/hp 要素法フレームワークである Nektar++ を用いた、直線ブレードを有するダリウス型 VAWT の幾何学的に解像された直接数値シミュレーション（DNS）を提示する。本研究は、1 枚、2 枚、3 枚のブレードを有する 3 つのタービン構成（NACA 0018 翼型）を、tip-speed ratio（$\lambda = 1.5, 2.5, 4.5$）の範囲で調査する。このパラメータ空間は、静的 solidity（$\sigma$）が 0.2、0.4、0.6 に対応し、動的 solidity（$\sigma_D$）は 0.47 から 0.94 の範囲にある。

シミュレーションでは、動的な再メッシュ生成なしに回転幾何形状を処理するために、移動座標系（MRF）定式化を採用している。流れは、弦長に基づくレイノルズ数 $Re_c = 10^4$ における非圧縮性ナビエ - ストークス方程式としてモデル化されており、この値は都市環境や宇宙環境を代表し、かつ基礎物理学的な関連性を維持しつつ数値的な実現可能性を確保するために選択された。メッシュ解像度は、$\Delta y^+ \lesssim 1$ となるように DNS の要件を満たし、流線方向には 32 個のフーリエ平面を有している。数値的安定性のため、解像度が不足している領域にはスペクトル減衰粘性（SVV）が適用されている。この手法は Battisti らの実験データに対して検証されており、近傍後流の位置における時間平均された流線方向速度プロファイルにおいて強い一致を示している。

主要な貢献と結果

1. 動的ストールと BVI の解像：
   DNS は、境界層のロールアップ、DSV の形成、分離、および対流を含む、完全な動的ストールサイクルを成功裡に捉えている。新たな発見として、BVI に駆動される「渦衝突メカニズム」の同定がある。

   • 1 ブレード構成：DSV は形成され、分離し、一貫した大規模構造として近傍後流へ対流し、風下側で著しい流れの加速を引き起こす。
   • 2 ブレード構成：solidity の増加は DSV の対流を抑制し、近傍後流に到達する前に拡散してより小さな構造へと崩壊させる。
   • 3 ブレード構成：ブレード間隔の縮小は、風上通過の開始時に BVI を引き起こす。これらの相互作用は、動的ストールの初期段階（層流分離バブルの発達）において生じ、LSB を擾乱し、最終的な DSV のサイズを縮小させる。さらに、ブレードが風下領域に接近するにつれ、上流後流からの直接的な渦衝突により、DSV が複数の明確なコアへと「早期に崩壊」する。このメカニズムは、一貫した渦ダイナミクスから鈍体ダイナミクスへの遷移を加速させる。

2. 後流回復と鈍体ダイナミクスへの遷移：
   本研究は、後流の下流への進化を定量化する。高い solidity（より多くのブレード）は、より深刻な初期速度欠損をもたらすが、より早く、より迅速に回復を開始する。

   • 低 solidity（ブレード数=1）：回復は、大規模な DSV によって誘起される横方向混合によって支配され、回復が遅く、後流の非対称性が持続する。
   • 高 solidity（ブレード数=3）：BVI によって駆動される DSV の崩壊は、VAWT 後流を鈍体から区別する一貫した構造を除去する。その結果、後流は、本質的に速いせん断層関連の回復メカニズムへとより迅速に遷移する。

3. 自己相似性解析：
   自己相似性解析を近傍後流へ拡張し（以前は ALM 研究において遠方後流に限定されていた）、著者らは、後流がブレードによって生成された一貫構造への依存性を失い、自己相似的なガウス解へと収束する速度を定量化した。

   • 解析により、この遷移の速度を設定する上で、ブレード数は tip-speed ratio よりもより影響力があることが明らかになった。
   • ブレード数が多いほど、自己相似性の開始が早くなる。興味深いことに、遷移の速度は最高値の $\lambda$ ではなく、$\lambda=2.5$ でピークに達しており、特定の動的 solidity を超えるとプラトー効果が生じることを示唆している。

意義と主張
本論文は、幾何学的に解像された DNS アプローチが、ブレード数が VAWT の空力特性に影響を与える特定のメカニズムを同定するために不可欠であると主張している。主要な意義は、「BVI によって誘起される DSV の早期崩壊」というメカニズムの発見にあり、これはブロック効果とは独立して、鈍体後流ダイナミクスへの遷移を加速させる。

著者らは、これらの知見が、密に配置されたタービンアレイおよび結合ペア構成の設計に直接的な示唆を持つと述べている。下流ローターが経験する流入は、本質的に上流タービンのブレード数に依存するため、ブレード数の少ないペアは大きく一貫した DSV に遭遇するのに対し、ブレード数の多いペアはより拡散的で擾乱された場に出会うことになる。論文は、これらの異なる流入条件が下流性能にどのように影響するかを完全に理解するために、ペア構成の DNS が正当化されると結論付けており、高い solidity におけるより速い後流回復とより早期の自己相似性開始が、より短い下流距離においてそのようなアレイのより正確な解析モデル化を可能にするかもしれないと指摘している。