コアンダ効果シミュレーター 戻る
流体力学シミュレーター

コアンダ効果シミュレーター — 壁面付着噴流の付着角度と反力

円弧曲面に沿うコアンダ噴流の付着角度・剥離点・運動量反力を可視化。噴流厚さ・曲面半径から曲率パラメータ K=b/R を変えて、いつ噴流が剥離するかを学べます。

パラメータ設定
噴流速度 U_0
m/s
噴流厚さ b
mm
曲面半径 R
mm
曲面延長角 θ_geo
°

空気 ρ = 1.2 kg/m³、ν = 1.5×10⁻⁵ m²/s、K_crit = 0.25 を仮定した半経験モデルです。

計算結果
曲率パラメータ K=b/R
剥離限界角 θ_s_max
実際の付着角度 θ_e=min(θ_geo, θ_s)
鉛直方向反力 F_y/単位幅
曲面と噴流軌跡

灰色=円弧曲面/青実線=付着区間/青破線=剥離後の直線軌跡/赤×=剥離点/黄矢印=初期・最終速度ベクトル

理論・主要公式

曲面に沿う噴流は、遠心力 $\rho U^2/R$ と圧力差・表面摩擦の釣合いで付着が決まります。本ツールでは曲率パラメータと剥離限界角を線形にモデル化しています。

曲率パラメータ K。b は噴流厚さ、R は曲面半径:

$$K = \frac{b}{R}$$

剥離限界角 θ_s(簡易モデル、K_crit = 0.25、θ_max = 180°):

$$\theta_s = \theta_{\max}\left(1 - \frac{K}{K_{\text{crit}}}\right),\quad K < K_{\text{crit}}$$

実際の付着角度 θ_e と単位幅当たり運動量反力:

$$\theta_e = \min(\theta_{\text{geo}},\,\theta_s)$$ $$F_x = \rho U_0^2 b \sin\theta_e,\quad F_y = \rho U_0^2 b (1-\cos\theta_e)$$

本ツールは半経験的なシンプルモデルです。実機評価には CFD・風洞実験が必要です。

コアンダ効果シミュレーターとは

🙋
スプーンの背を蛇口の水に近づけると、水がスプーンに「吸い付く」みたいに曲がりますよね。あれって何ですか?
🎓
それがコアンダ効果だよ。ざっくり言うと、噴流が近くにある壁面に沿って曲がる現象だ。1930年にルーマニアの発明家アンリ・コアンダが報告したから、彼の名前が付いている。噴流が壁に近づくと、壁側だけ周りの空気を巻き込めなくて圧力が下がる。すると圧力差で噴流が壁にくっつく方向に引っ張られる。シミュレーターで「噴流厚さ b」を小さくしてみて。曲面に沿ってきれいに曲がるのが見えるよ。
🙋
じゃあ「曲率パラメータ K」のスライダーを動かすと、どうなるんですか?
🎓
K=b/R は、噴流厚さと曲面半径の比だ。直接動かすスライダーじゃないけど、b を大きくするか R を小さくすると K が増えて、噴流が曲がりきれなくなる。シミュレーターで「曲面半径 R」を 50mm くらいまで小さくしてごらん。剥離限界角 θ_s_max がガクンと下がって、赤い × 印(剥離点)が手前に現れるはず。経験的に K_crit ≈ 0.25 が限界で、それを超えると最初から剥離して付着しない。
🙋
「曲面延長角 θ_geo」と「剥離限界角 θ_s」は何が違うんですか?
🎓
θ_geo は曲面そのものの長さ(角度)で、θ_s は流体力学的に付着していられる限界角度。実際の付着角度 θ_e は両者の小さい方になる。例えば曲面が 90° しかなくて θ_s が 108° あれば、噴流は曲面の終わりまで付着して θ_e = 90°。逆に曲面が 180° あっても θ_s が 50° しかなければ、50° で剥がれて残りの 130° は空気中を直進する。
🙋
下のカードに「F_y = 4.80 N/m」って出てます。これは何の力ですか?
🎓
運動量変化に伴う反力だよ。噴流が θ_e だけ方向を変えると、その分の運動量変化が壁面に反力として伝わる。$F_y = \rho U_0^2 b(1-\cos\theta_e)$ で、付着角度が 90° なら $1-\cos 90°=1$ で最大の鉛直方向反力が出る。実は V-22 オスプレイのチルトローター制御や F1 のリアウィングの空力は、この運動量変化を意図的に使っているんだ。本ツールは半経験的なシンプルモデルだから、実機設計には CFD が必要だけど、現象の本質はつかめる。

よくある質問

V-22 の主翼はローターからの強い下降流(ダウンウォッシュ)を直接受けるため、ホバリング時の効率低下や振動の原因になります。コアンダ効果を利用したフラップ位置と翼面形状の最適化で、下降流を翼面に沿わせて整流したり、推力方向の微調整に使ったりします。一部の研究機(例:HondaJet 風洞研究、X-32 系の研究)ではより積極的に推力偏向制御に応用されています。
2010年代の F1 では「コアンダ排気」が流行しました。エキゾーストガスをサイドポッドの曲面に沿って流し、リアディフューザー入り口まで導いてダウンフォースを稼ぐ手法です。リアウィングのフラップ形状も、上面で噴流を付着させて剥離を遅らせ、ストール直前まで揚力(ダウンフォース)を引き出すために最適化されます。曲率パラメータ K の管理が設計の核心です。
天井埋め込み型エアコンの吹き出しノズルは、コアンダ効果で天井面に沿わせて気流を遠くまで送り、室内全体への混合を促します。ダイソンのブレードレス扇風機は、リング内側の薄いスリットから噴出した高速気流が、リング外側の曲面に沿って流れることで周囲空気を巻き込み、入力空気量の 15〜20 倍の大風量を作り出しています。いずれも壁面付着と周囲流体の巻き込みを使う点で同じ原理です。
本ツールは曲率パラメータ K と剥離限界角 θ_s を線形でモデル化した半経験的なシンプルモデルです。実際の剥離挙動は Re 数、壁面粗さ、噴流の初期乱れ、外乱(音響・振動)に強く依存し、定量予測には CFD(RANS や LES)や風洞試験が必要です。本ツールは「曲率を大きくすると剥離が早まる」「θ_e から運動量反力が決まる」といった現象の定性的理解を得るための教育用ツールです。実機設計の数値を直接使うことはできません。

実世界での応用

航空機の高揚力装置:STOL(短距離離着陸)機や艦載機では、フラップ前縁から噴流を吹き出し、コアンダ効果で翼面に付着させて剥離を遅らせる「ブローイング・フラップ」が古くから使われています。NASA の QSRA(Quiet Short-Haul Research Aircraft)など、上面ブロー・ジェット・フラップ機が低速で巨大な揚力を発生できるのはこの原理です。

ティルトローター機・推力偏向:V-22 オスプレイをはじめとするティルトローター機では、ローター下流の翼面形状とフラップ位置をコアンダ効果を踏まえて設計し、ダウンウォッシュの整流や推力ベクター制御に活用します。Bell-Boeing による研究では、コアンダ効果を利用した翼形状で操縦応答性が大幅に改善されることが報告されています。

F1・モータースポーツの空力:2010年代の F1 で流行した「コアンダ排気」、リアウィングやディフューザーの曲面設計、サイドポッドの形状は、いずれもコアンダ効果で剥離を遅らせて空気を意図した方向に流すための工学的最適化の結果です。曲率パラメータの管理がレギュレーション内で許される設計余地の鍵となります。

HVAC・産業用ノズル:天井埋め込み型エアコン、レンジフード、産業炉のガスバーナーノズル、製紙工場のエアナイフなど、噴流を壁面に沿わせて遠くまで運ぶ用途で広く使われます。ダイソンのブレードレス扇風機は、薄いスリットから出た噴流が曲面に付着して周囲空気を引き込む「インダクション」効果で大風量を作る代表例です。

よくある誤解と注意点

最も多い誤解は、「コアンダ効果はベルヌーイの定理だけで説明できる」と考えてしまうことです。一般向け解説では「壁側で流れが速くなって圧力が下がる」と説明されがちですが、これは現象の半分だけ。実際には噴流が壁面に近づくと周囲流体の巻き込み(エントレインメント)が壁側で阻害され、その結果生じる低圧領域が噴流を壁面に引き寄せます。粘性と乱流による運動量交換の効果が本質で、ベルヌーイの定理(非粘性・非回転)だけでは正確に説明できません。シミュレーターで b と R を変えると K=b/R で挙動が決まるのは、この相互作用の幾何的な「曲げの厳しさ」を表しているからです。

次に多いのが、「K が小さければ何度でも付着し続ける」と思い込むことです。本ツールの簡易モデルでは θ_s_max が最大 180°(半円分)ですが、実際には壁面摩擦と周囲流体の巻き込みで噴流のエネルギーが消費され、十分長い曲面では運動量が減衰して剥離します。さらに Re 数が低いと層流状態で曲げに弱く、Re 数が高いと乱流の運動量混合で曲げに強くなるなど、Re 依存性も無視できません。本ツールはこれらを線形 K モデルで簡略化しているため、絶対値よりも「曲率を増やすと剥離が早まる」という傾向の理解に使ってください。

最後に、このシミュレーターは「壁面付着の限界角度と運動量反力」を計算するもので、噴流の速度プロファイルや圧力分布は扱っていない点に注意してください。実際の剥離点予測には壁面圧力勾配・粘性応力・乱流モデルを含む CFD 計算が必要で、F1 や航空機開発では数百万メッシュの RANS/LES 計算と風洞試験を組み合わせます。本ツールは現象の本質をつかむ第一歩であり、F_y = 4.8 N/m といった数値は単位幅当たりの理論最大値であって、実機の設計値ではないと理解してください。