› 潮流発電タービン Cp・TSR シミュレーター 戻る
海洋再生エネ・潮流発電

潮流発電タービン Cp・TSR シミュレーター — Lanchester-Betz

潮流発電タービンの心臓部「出力係数 Cp」と「翼端速度比 TSR」を、水平軸 HATT・垂直軸 Darrieus / Savonius・ダクト付の 4 形式で比較します。Lanchester-Betz 限界 16/27 ≈ 59.3% に対してどこまで出力を引き出せるか、抽出出力 kW・年間発電量 MWh・キャビテーション σ までリアルタイムに評価できます。

パラメータ設定
タービン形式
HATT/Darrieus/Savonius/ダクト付で Cp_max と最適 TSR が変わる
タービン径 D
m
潮流速度 V
m/s
大潮時のピーク流速。商用案件は 2〜4 m/s 級が主流
回転数 N
RPM
ブレード数 B
ピッチ角 θ
°
負値はフェザリング、+15°超は失速で出力低下
設置深度 h
m
海面下のローター中心深さ。σ 評価に使用
計算結果
—
翼端速度比 TSR
—
出力係数 Cp
—
Betz 限界対比 (%)
—
抽出出力 (kW)
—
年間発電量 (MWh/y)
—
キャビテーション σ
—
海中タービン断面図 — ブレード回転・潮流ベクトル

海面下に設置したタービンに潮流が流れ込み、ブレードが回転。右上の Cp ゲージは現在の出力係数を Lanchester-Betz 限界に対して表示します。

Cp – TSR 曲線（現在の形式）
形式別 Cp_max 比較
理論・主要公式

$$P = \tfrac{1}{2}\rho A V^{3} C_p,\qquad TSR = \frac{\Omega R}{V},\qquad C_{p,\max} = \frac{16}{27} \approx 0.593$$

抽出出力 P と翼端速度比 TSR。ρ=1025 kg/m³（海水）、A=πR² はローター掃過面積、Ω はローター角速度、V は自由流速。Lanchester-Betz 限界 16/27 が開放流れの理論上限。

$$C_p(\lambda) \approx C_{p,\max}\,\exp\!\left[-\left(\frac{\lambda-\lambda_{opt}}{3}\right)^{2}\right]$$

単純化した Cp(λ) ベル形曲線。λ=TSR、λ_opt はタービン形式依存（HATT 6・Darrieus 4・Savonius 1・ダクト付 5）。

$$\sigma = \frac{p_{atm}+\rho g h - p_{vap}}{\tfrac{1}{2}\rho U_{tip}^{2}}$$

キャビテーション数 σ。h は設置深度、U_tip=ΩR は翼端速度、p_vap≈2500 Pa（15°C）。σ<0.3 で高リスク。

潮流発電タービン 出力係数 Cp — Lanchester-Betz 限界

🙋
潮流発電って、海の中に「水中の風車」を置くやつですよね？風力発電とどこが違うんですか？
🎓
基本原理は同じで、流体の運動エネルギーをローターで受けて発電する。ただ流れる媒体が「空気 → 海水」に変わるだけで、密度が約 800 倍違うんだ。だから同じローター径・同じ流速なら理論出力は 800 倍。実際の潮流流速は風速より遅いから、現場の出力比はだいたい 10〜20 倍くらいになる。スコットランドの MeyGen プロジェクトみたいに、4 基で 6 MW を出している商用機がもう動いてるよ。
🙋
なるほど！じゃあ Cp っていう係数が大事って聞いたんですが、これは何ですか？さっきの 16/27 って何の数字なんでしょう。
🎓
Cp は「出力係数 (Power Coefficient)」で、流れ込んでくるエネルギーのうち何 % をローターで取り出せたかの効率指標。式で書くと Cp = P / (½ρAV³) で、A はローターの掃過面積、V は流速。そして 16/27 ≈ 0.593 が有名な Lanchester-Betz 限界。Lanchester が 1919 年、Betz が 1920 年に独立に「開放流れ中の単一ローターはこれ以上は理論的に取れない」と証明したんだ。風力も潮流も同じ上限で、商用 HATT で 0.40〜0.45、ダクト付なら見かけ上 0.55+ までいける。
🙋
えっ、ダクト付って Betz 超えてもいいんですか？それは「タービンの限界突破」みたいなものですか？
🎓
良い質問。実は Betz の式は「ローター面積基準」で書かれているから、ダクト（ディフューザー）で流れを集めるとローター面を通る流量自体が増える。だから「ローター面積基準の Cp」は 0.6〜0.7 まで行く。でもダクト出口面積基準で考え直せば 16/27 を破ってない。理論を破ってるんじゃなくて、定義の取り方の問題なんだ。Sabella D10（フランス）や OpenHydro が実用化していて、潮流が弱い海域でも実装しやすいメリットがある。
🙋
TSR っていうスライダーを動かすと Cp がベル型に変化しますね。これは何ですか？
🎓
TSR は「翼端速度比 (Tip Speed Ratio)」で、TSR = ΩR/V＝翼端の周速÷流速。低すぎると後流の旋回損失が大きく、高すぎると抗力と乱流剥離が増える。だから最適 TSR が存在する。HATT は 5〜7、Darrieus は 4、Savonius（抗力型）は 1 付近が最適点。実機では流速が満ち引きで変わるので、ピッチ角と回転数を可変制御して常に TSR_opt に張り付くようにする。これが Cp_max を維持するコツだよ。
🙋
最後にキャビテーション σ ってありますね。これは何のチェックですか？
🎓
翼端は最速で動くから、そこで局所圧力が水の飽和蒸気圧（15°C で約 2500 Pa）を下回ると気泡が発生する＝キャビテーション。気泡が崩壊するとマイクロジェットで翼面がピッティング損傷する。σ がその余裕度。σ<0.3 で高リスク、1.0 以上で安全。対策は深度を深くする、回転数を下げて翼端速度を抑える、専用断面（NACA 6 系等）を使う、の 3 つ。プロペラ船の世界では昔からの問題で、潮流タービンでもまったく同じ判断が必要になる。

よくある質問

理想的なディスクアクチュエータ理論（運動量理論）から導かれる、開放流れ中の単一タービンが流体から抽出できる運動エネルギーの理論最大割合です。値は 16/27 ≈ 0.593（59.3%）で、1919 年に Lanchester、1920 年に Betz が独立に導出しました。風車・潮流タービン共通の上限ですが、ダクト（diffuser augmented）で流れを集めるとローター面を通過する流量が増え、ローター面積基準の Cp は 0.6〜0.7 に達することがあります。
TSR = ΩR/V（翼端の周速÷自由流速）が低すぎるとブレード後流の旋回損失が大きく、高すぎるとブレード間の抗力損失と乱流剥離が支配的になります。水平軸潮流タービン（HATT）では TSR=5〜7 付近で Cp が最大になり、Darrieus 垂直軸は TSR=4 付近、Savonius（抗力型）は TSR=1 付近が最適です。ローター回転数を流速に合わせて変えるピッチ・速度可変制御で、潮の満ち引きを通じて Cp を最大点近くに維持するのが商用機の主流です。
出力は流体密度 ρ に比例し、海水の密度は 1025 kg/m³ と空気 1.225 kg/m³ の約 837 倍です。同じ流速・同じローター径なら理論出力は約 800 倍。実際の潮流流速は風速より小さい（典型 2〜3 m/s vs 風 8〜12 m/s）ものの、流速の 3 乗で効くため、それでも 10〜20 倍規模の出力密度になります。さらに潮汐は天体力学で決まるので 12.4 時間周期で正確に予測でき、電力系統に組み込みやすいのも風力にない利点です。
σ = (p_atm + ρgh − p_vap) / (½ρU²) で、翼端付近で局所圧力が飽和蒸気圧を下回ると気泡が生じ（キャビテーション）、崩壊時の衝撃で翼面がピッティング損傷します。σ < 0.3 で高リスク、0.3〜1.0 で要注意、1.0 以上で概ね安全とされます。対策は (1) 設置深度を深くして p を上げる、(2) ローター径を大きくし回転数を下げて翼端速度を抑える、(3) 翼断面をキャビテーション特化型（NACA 6-series 等）に変える、の 3 つが主流です。

実世界での応用

MeyGen（スコットランド・ペントランド海峡）：世界最大級の潮流発電商用プロジェクト。SIMEC Atlantis Energy が運営する 4 基の AR1500（径 18m、定格 1.5 MW）で合計 6 MW を系統連系。本ツールのデフォルトパラメータ（D=18m, V=2.5m/s）に近い設計で、実機の Cp は 0.40〜0.45 のレンジに収まります。フェーズ拡張で 86 MW、最終的に 398 MW を目指す計画です。

Sabella D10（フランス・ブルターニュ沖）：水平軸ダクト無しタイプで径 10m、定格 1 MW。ダクト無しの代わりにブレードを 6 枚にし、低流速海域（2 m/s 前後）でも安定発電できる設計。ヴェルダン島の系統に直結して離島電力を賄った実績があり、本ツールの形式選択「水平軸 HATT」+「径 10m」で再現できます。

OpenHydro / ダクト付タービン：ナバル・エナジー社が開発したダクト型タービンで、ローターをリング状ハブの内側に収めた特徴的なフォルム。本ツールで「ダクト付」を選ぶと Cp_max が 0.55 に上がるのを確認できます。カナダ Bay of Fundy（潮位差 16m）で実証試験が行われましたが、極端な乱流環境で機械的損傷が頻発し、技術課題の難しさも浮き彫りになりました。

BEM 解析・CFD の事前検討：詳細な Blade Element Momentum 法（BEM）や CFD（OpenFOAM の interFoam / k-ω SST）を回す前に、本ツールのような Cp(TSR) 概算で「形式選定・ローター径・回転数の当たり」を付けます。CFD で出した Cp が概算と桁違いなら、メッシュ品質や乱流モデルの設定ミスを疑うサニティチェックにも使えます。

よくある誤解と注意点

最大の落とし穴が、「ダクト付タービンが Betz 限界を破る」という誤解です。実際にはダクトが集流効果で流量を増やしているだけで、ダクト入口の自由流断面積で Cp を再定義すれば 16/27 は破られていません。論文を読むときは「Cp の基準面積が ローター面積 / ダクト入口面積 / 仮想自由流面積 のどれか」を必ず確認してください。比較対象を間違えると、ダクト付が水平軸の 50% 増しに見えたり、逆に等価に見えたりします。本ツールはローター面積基準で書いているので、ダクトの Cp_max=0.55 は便宜上の値です。

次に、「年間発電量＝定格出力 × 8760 時間」という単純計算。潮流は 12.4 時間周期で流速が 0 ⇔ ピークを正弦的に変動するため、平均流速は√(2)/2 倍に近く、平均出力はピークの 1/8〜1/4 になります。さらに保守・潮流方向反転時のシャットダウンも入って、設備利用率（capacity factor, CF）は 30〜40% が現実値。本ツールでは CF=0.35 を仮定していますが、サイト個別の海象データで補正してください。CF=0.5 を超える数字を提示する事業計画は楽観的すぎる可能性が高いです。

最後に、キャビテーション σ の評価で「平均的な静圧」だけ見て安全と判断する罠。実際には翼端の前縁付近で局所圧力が静圧より 2〜5 倍低く落ち込むため、平均 σ が 1.0 を超えていても局所的にキャビテーションが発生することがあります。BEM や CFD で翼面の Cp 分布を確認し、最小圧係数 Cp_min < −(σ − 0.05) になっていないかを必ず別途チェックしてください。設置深度を深くするのが最も効きますが、潮位差が大きい海域では干潮時の海面下クリアランスも考慮する必要があります。

使い方ガイド

  1. ロータ直径（0.5～12 m）、潮流速度（0.5～3.0 m/s）、回転数（5～200 rpm）、翼枚数を入力します
  2. シミュレータが翼端速度比 TSR と出力係数 Cp を自動計算し、Lanchester-Betz 限界（16/27 = 0.593）との比較を表示します
  3. 抽出出力（kW）、年間発電量（MWh/y）、キャビテーション数 σ を確認し、設計パラメータの最適化に活用します

具体的な計算例

直径 3.5 m の水平軸タービン（HATT）が秋田県沖の潮流 1.8 m/s、120 rpm で運転される場合：翼端速度比 TSR ≈ 4.2、最適 Cp ≈ 0.48（Betz 限界の 81%）、抽出出力 ≈ 24.3 kW。年間潮流稼働率 38% を仮定すると年間発電量 ≈ 80.4 MWh/y。キャビテーション数 σ ≈ 1.1（閾値 0.8 以下で気泡発生）で、翼面圧力低下による侵食リスク評価が可能です。

実務での注意点