Numerical simulations of transition and long-term response of a wind turbine airfoil

本論文は、低レイノルズ数におけるFFA-W3翼型の遷移および長期的な流れの応答を解析するためにNek5000およびEllipSysを用いた数値シミュレーションを提示し、ケルビン・ヘルムホルツ不安定性を捉えるEllipSysコードの能力を検証するとともに、低周波振動および潜在的なバブル崩壊に関連する揚力係数の緩やかな周期的変調を特定するものである。

要約

風力タービンのブレードを、巨大で回転する翼だと想像してみてください。飛行機の翼と同じように、効率よく機能するためには滑らかな気流が必要です。しかし、風がある特定の角度で当たると、空気が表面から「剥がれ」、混沌とした渦巻く塊を作ってしまうことがあります。この論文は、ハイテクな風洞実験のようなものです。ただし、物理的なモデルを使う代わりに、研究者たちはスーパーコンピューターの中に仮想的なモデルを構築し、空気がどのように振る舞うかを正確に観察しました。

彼らの発見の物語を、シンプルな概念ごとに解説します：

1. 仮想風洞

研究者たちは、巨大な風力タービン（10MW級）のブレードの特定の断面を研究したいと考えました。彼らは、NEK5000とELLIPSYSという2つの異なるコンピュータプログラムを使用して、ブレードの上を流れる空気をシミュレーションしました。

NEK5000を、あらゆる微細なディテールを捉える超高精度で高性能なカメラだと考えてください。非常に精密ですが、実行には時間がかかりコストも高いものです。一方、ELLIPSYSは、もう少し高速で効率的なカメラのようなものです。チームはまず、「より速い」カメラ（ELLIPSYS）が「高性能な」カメラと同じものを見ることができるかどうかを証明する必要がありました。その結果、ELLIPSYSは滑らかな空気の中にある極めて微細でかすかな波紋（滑らかにしすぎてしまうため）は見逃すものの、ブレードの性能に実際に影響を与える大きな混沌とした渦を捉えることには非常に優れていることがわかりました。

2. 風洞の幅はどのくらい必要か？

長いシミュレーションを実行する前に、彼らは仮想的な「風洞」をどの程度の幅にする必要があるかを判断しなければなりませんでした。風洞が狭すぎると、空気が押しつぶされて偽の結果が生じる可能性があります。逆に広すぎると、コンピュータの計算資源を無駄にしてしまいます。

彼らは、「狭い」風洞（翼幅の10%）を「広い」風洞（翼幅の20%）と比較テストしました。

• 例え： 川の流れを見ているところを想像してください。もし川の細い帯の部分だけを見ていたら、大きな波を見逃してしまうでしょうか？
• 結果： 彼らは、狭い風洞でも実際には十分であることを発見しました。大きな波や渦は、その狭い空間の中でも完璧に形成されました。これにより、より小さな、より狭いシミュレーションボックスを使用することで、コンピュータの計算時間を大幅に節約できることがわかりました。

3. 「バブル」と「フラップ」

この研究で最も興味深い部分は、翼の上部（吸い込み側）で起こりました。

• 剥離バブル（セパレーション・バブル）： 空気が翼の上を流れる際、一時的に表面から剥がれ、小さな再循環する空気のポケット、すなわち「層流剥離バブル（LSB）」を作り出します。これは、翼の表面にできる小さな一時的な渦巻きのようなものです。
• 不安定性： このバブルの中で、空気はただ静止しているわけではありません。空気は振動し、波のように丸まっていきます（池に広がる波紋のようなものです）。研究者たちは、これらの波がどのように成長するかを観察しました。彼らは、このバブルにおける主要な「ローラー」が、ケルビン・ヘルムホルツ不安定性と呼ばれる一種の不安定性であることを突き止めました。
• 発見： 彼らは、この「より速い」コンピュータプログラム（ELLIPSYS）が、これらの波がどのように成長し、バブルがどのように振る舞うかを正確に予測できることを確認しました。これは、超高精度なプログラムの結果と一致していました。

4. ゆっくりとしたリズムの鼓動（大きな驚き）

ツールの検証を終えた後、彼らは非常に長い時間（空気が翼の横を通り過ぎる時間の50倍に相当する時間）シミュレーションを実行しました。ここで、彼らは特別なものを発見しました。

空気が高速で混沌とした動きをしている一方で、翼を押す力に非常にゆっくりとしたリズムの脈動があることに気づいたのです。

• 例え： ドラムのビートを想像してください。空気の高速な乱れは、素早い高音のドラムロールのようなものです。彼らが見つけたゆっくりとした脈動は、シミュレーション時間で48秒に一度発生する、深くゆっくりとした心拍のようなものです。
• 影響： このゆっくりとした心拍により、翼を押す力は上下に約**10.5%**も揺れ動きました。
• 実機タービンとの関連： これを実際の回転する風力タービンに当てはめると、このゆっくりとした脈動は、ブレードが7.7フル回転するごとに一度起こることに気づきました。

5. なぜこれが起こるのか？

研究者たちは、このゆっくりとした脈動は、空気が翼の上で「失速（ストール）」することと「失速解除（アンストール）」することのサイクルによって引き起こされていると考えています。

• サイクル： 空気が引っかかり、大きなバブルを作り出します。次に、何らかのきっかけでバブルが弾け、空気が再び滑らかに再付着します。その後、再び圧力が上昇し、バブルが形成されるというサイクルが繰り返されます。
• トリガー（引き金）： 彼らは、翼の上で空気が非常に強く後方へと渦巻いているため、「絶対不安定性」と呼ばれる状態が生じていると考えています。これは、空気が非常に乱れているため、自律的に振動せざるを得ない状態を指す専門的な表現です。

6. 結論

この研究は、コンピュータ・モデリングの成功例です。彼らは、まず「ゴールドスタンダード（最高基準）」と比較して検証を行えば、より速く効率的なコンピュータプログラム（ELLIPSYS）が複雑な風力タービンの物理現象を研究するために信頼できるものであることを証明しました。

彼らは、厚みのある風力タービンのブレードであっても、空気の流れにはおよそ8回転に一度起こる、非常にゆっくりとしたリズムの「呼吸」が存在することを発見しました。この呼吸は、揚力（タービンを回転させる力）を大きく増減させます。このゆっくりとしたリズムを理解することは、これが直ちにタービンを破壊することはありませんが、長年の運用において、これらのゆっくりとした大きな力の変動が材料の疲労につながる可能性があるため、非常に重要です。

要約すると： 彼らは仮想風洞を構築し、高速なコンピュータがその役割を果たせることを証明し、そして風力タービンのブレードには、その表面で空気の泡が形成されたり弾けたりすることによって生じる、ゆっくりとしたリズムの「心拍」が存在することを発見したのです。

技術概要

技術要約：風力タービン翼型の遷移および長期応答に関する数値シミュレーション

問題提起
高い流入角と厚い翼型形状に起因する風力タービンブレード上の流れの剥離は、性能および構造負荷に大きな影響を与える。先行研究では、剥離流における複数の固有振動数（後流の放出、ケルビン・ヘルムホルツ（KH）モード、および低周波振動（LFO）を含む）が特定されている。LFOは強い揚力変動を特徴とし、多くの場合、層流剥離泡（LSB）内における絶対不安定性に起因する、流れの周期的な失速および失速からの回復に関連していると考えられている。LFOは様々な翼型の研究で観察されてきたが、厚い風力タービン翼型の低い迎角において、それらがどのように発生するのか、またその条件下での具体的な駆動メカニズムは、さらなる調査が必要である。本研究では、翼弦レイノルズ数（$Re_c$）が $1 \times 10^5$ のDTU 10-MW参照風力タービン翼型（FFA-W3シリーズ）の一節に焦点を当てる。

手法
本研究では、非圧縮性ナビエ・ストークス方程式を解くために、サブグリッドスケールモデルを用いない壁解像型ラージ・エディ・シミュレーション（LES）を採用している。以下の2つのソルバーを使用する：

1. NEK5000: 数値散逸が極めて少なく、高精度なスペクトル要素法（SEM）コード。
2. ELLIPSYS: DTU/Risøで開発された有限体積法コード。従来はRANSに使用されてきたが、ここではインプリシットLESとして適用される。

研究には、ELLIPSYSをNEK5000に対して系統的に検証するプロセスが含まれる。主な手法的なステップは以下の通りである：

• ドメイン感度: 関連する流れ構造の解像度を確保するため、$L_z = 0.1c$ と $L_z = 0.2c$ のシミュレーションを比較することで、スパン方向のドメイン幅（$L_z$）を評価する。
• ソルバーの検証: 平均流場、摂動の進化、およびスペクトル特性について、両ソルバー間の比較を行う。
• 安定性解析: パラボリック安定方程式（PSE）およびスペクトル固有直交分解（SPOD）を利用して、不安定性（トールミーン・シュリヒトヒンゲン（TS）モードおよびKHモード）の成長と、コヒーレントモードの空間構造を解析する。
• 長期積分: 長期的な現象を捉えるため、計算効率の高いELLIPSYSソルバーを用いて、50通過時間（$T = 50c/U_\infty$）にわたる流れの進化をシミュレートする。

主な結果

• ドメイン幅: スパン方向の幅は、主要な摂動の進化を捉え、時間平均の流れの統計量を再現するために、翼弦の10%（$L_z = 0.1c$）で十分であることが判明した。これは、より広いドメイン（$L_z = 0.2c$）の結果と一致している。
• ソルバーの検証: ELLIPSYSは、平均流場および最も増幅された摂動の進化においてNEK5000と密接な一致を示した。しかし、ELLIPSYSは数値散逸が高いため、付着した境界層内のTS波の振幅を過小評価した。逆に、ELLIPSYSは層流剥離泡（LSB）内におけるKHモードの進化を正確に予測し、PSEの予測とよく一致した。
• 遷移メカニズム: 吸込面の遷移は、剥離したせん断層上で発達するKH不安定性によって駆動される。最も増幅されたKHモードの周波数は $f \approx 15$ であり、その下流への進化は両方のソルバーおよびPSEによって良好に捉えられている。
• 低周波振動（LFO）: 長期シミュレーションにより、垂直抗力係数（$C_N$）の緩やかな変調が明らかになった。その振幅は10.5%であり、周期は48通過時間であった。これは周波数 $f = 0.021$ およびストラハル数 $St = 0.0012$ に相当する。
  • この周波数は、主要なKH不安定性の周波数よりも約720倍低い。
  • DTU 10-MWタービンの文脈では、この周期は約7.7回転分に相当する。
  • 観測された $St$は、他の翼型について報告されている値よりも低く、比較的小さな実効迎角（$3.1^\circ - 3.3^\circ$）で発生している。

意義および主張
本論文は、ELLIPSYSが遷移流の壁解像型LESのための実行可能なツールであり、長時間の積分においてスペクトル要素法コードに代わる計算効率の高い選択肢であることを確立している。本研究は、高い逆流（吸込側で最大-17%）による絶対不安定性に誘発される可能性がある、流れの周期的な失速および失速からの回復が、観察されたLFOの原因であることを確認した。

重要な知見は、厚い翼型（厚さ24%）において、低い迎角でLFOが特定されたことである。この条件は、文献に見られる薄型翼型の失速とは異なる。著者らは、翼型の両側における逆流が絶対不安定性を許容するのに十分であり、それがLSBの完全な失速状態への周期的なバーストを引き起こしている可能性があると述べている。本研究は、厚い翼型形状が低周波振動の発生と特性にどのように影響するかについての理解に寄与しており、これらの現象が標準的な翼型で記録されているものよりも低い入射角およびストラハル数で発生し得ることを強調している。