乱流境界層シミュレーター 戻る
流体力学シミュレーター

乱流境界層シミュレーター — 平板上の1/7乗則とCf

平板上の乱流境界層を1/7乗則で可視化。流速・位置・動粘度・密度を変えて、境界層厚さ δ・摩擦係数 C_f・壁面せん断応力 τ_w を計算し、Blasius層流解との違いを学べます。

パラメータ設定
流速 U
m/s
位置 x
m
動粘度 ν
m²/s
密度 ρ
kg/m³

遷移目安 Re_x ≥ 5×10⁵ で乱流式が適用されます。それ未満では層流域として警告を表示。

層流域・乱流式不適
計算結果
Reynolds 数 Re_x
乱流境界層厚さ δ_99
局所摩擦係数 C_f
壁面せん断応力 τ_w
平板上の境界層と速度プロファイル

水色矢印=主流 U/灰板=平板/黄色破線=境界層厚さ δ_99/青曲線=乱流(1/7乗則)、橙曲線=層流(Blasius近似)

速度プロファイル u(y)/U

横軸=無次元速度 u/U/縦軸=壁からの無次元距離 y/δ(青=乱流 1/7乗則、橙=層流 放物線近似)

理論・主要公式

平板上の乱流境界層は、1/7乗則速度分布から導かれる経験式で記述されます。層流(Blasius)解との対比で乱流の性質が明確になります。

Reynolds 数。U は主流速度、x は前縁からの距離、ν は動粘度:

$$Re_x = \frac{U\,x}{\nu}, \qquad \text{遷移目安:}\; Re_x \ge 5\times 10^5$$

乱流境界層厚さ・排除厚さ・運動量厚さ(1/7乗則):

$$\frac{\delta_{99}}{x} = 0.37\,Re_x^{-1/5}, \quad \frac{\delta^{*}}{x} = 0.046\,Re_x^{-1/5}, \quad \frac{\theta}{x} = 0.036\,Re_x^{-1/5}$$

局所摩擦係数と壁面せん断応力:

$$C_f = 0.059\,Re_x^{-1/5}, \qquad \tau_w = \tfrac{1}{2}\,C_f\,\rho\,U^2$$

層流(Blasius)解との比較:

$$\frac{\delta_{99}}{x} = 5.0\,Re_x^{-1/2}, \qquad C_f = 0.664\,Re_x^{-1/2}$$

形状係数 H = δ*/θ ≈ 1.28(乱流)、≈ 2.59(層流 Blasius)。同じ Re_x でも乱流のほうが Cf が数倍大きく、これが乱流摩擦抵抗の本質です。

乱流境界層シミュレーターとは

🙋
流体力学の本で「境界層」って出てきたんですが、平板の上に薄い層があるイメージくらいしか分かりません。乱流境界層って何が違うんですか?
🎓
ざっくり言うと、平板の前縁から下流に行くにつれて、壁の摩擦に引きずられて遅くなる流体の領域が境界層だ。最初は層流で薄く穏やかに発達するんだけど、Re_x が 5×10⁵ くらいを超えると渦が立って乱流に遷移する。乱流になると渦が運動量を激しく混ぜるから、層流より境界層がぐんと厚くなって、しかも壁面の摩擦も増える。上のシミュレーターで流速 U を上げてみて。Re_x が一気に大きくなって乱流域に入るのが分かるよ。
🙋
境界層厚さの式 $\delta/x = 0.37\,Re_x^{-1/5}$ ってどこから来てるんですか?
🎓
これは「1/7乗則」と呼ばれる経験的な速度分布 $u/U = (y/\delta)^{1/7}$ を使って、運動量積分方程式を解いた結果なんだ。層流の Blasius解が $\delta/x = 5.0\,Re_x^{-1/2}$ なのと比較すると、指数が違うことに注目してほしい。乱流は Re_x の -1/5 でゆっくり減るから、Re_x が大きくなっても乱流境界層は層流ほど薄くならない。シミュレーターの理論カードに両方並べてあるよ。
🙋
摩擦係数 C_f が乱流のほうが大きいのは、なんで困るんですか?
🎓
航空機の翼や船の船体の摩擦抵抗が、そのまま燃費に直結するからだ。空気で U=20 m/s、x=1 m、ν=1.5×10⁻⁵ の標準条件だと Re_x ≈ 1.33×10⁶ で完全に乱流域。C_f は約 3.5×10⁻³ で、同じ Re_x の Blasius層流 (0.664·Re_x^(-1/2) ≈ 5.8×10⁻⁴) と比べると6倍以上だ。だから航空機の設計者は「いかに乱流遷移を遅らせるか」「層流翼を設計できるか」に必死になる。シミュレーターで流速を変えて Cf の値の動きを見てみて。
🙋
壁面せん断応力 τ_w はどう使うんですか?
🎓
τ_w を平板全体で積分すれば、平板にかかる総摩擦抵抗(粘性抗力)が出る。CFDで境界層を解くときも、τ_w を直接計算するか、壁関数 (wall function) で近似するかの選択が大きな分かれ目だ。実務では「平板で τ_w がいくらか」を頭の中に持っておくと、CFDの結果が桁オーダーで合っているかすぐ判断できる。本ツールの既定値(空気・U=20 m/s)なら τ_w ≈ 0.84 Pa = 844 mPa。覚えておくと便利な数字だよ。

よくある質問

教科書的な目安は 5×10⁵ ですが、実機では条件で大きくばらつきます。表面粗さ、上流の乱れ、圧力勾配、振動などがあると 1×10⁵ 程度まで前倒しになります。逆に滑らかな平板で乱れの少ない風洞では 3×10⁶ まで層流が維持されることもあります。本シミュレーターは便宜的に 5×10⁵ を遷移目安として警告表示しています。
Re_x が 5×10⁵〜10⁷ 程度の範囲で、実験と数%以内で一致する優れた近似式です。それ以上の高 Reynolds 数では Schlichting の対数則ベース式(Cf = (2·log10(Re_x) − 0.65)^(-2.3) 等)のほうが精度が高くなります。逆に Re_x < 5×10⁵ では層流域なので Blasius解 Cf = 0.664·Re_x^(-1/2) を使う必要があります。
流れの「状態」を一目で判断する指標になります。平板の層流 (Blasius) で H ≈ 2.59、乱流(1/7乗則)で H ≈ 1.28、剥離の直前で H ≈ 3〜4 まで急増します。CFDの可視化で H 分布をプロットすると、層流-乱流の遷移位置や剥離の前兆をきれいに捉えられるため、境界層解析で重宝されます。
本ツールの式は「圧力勾配ゼロ(dp/dx = 0)の平板」専用です。逆圧力勾配(dp/dx > 0、減速場)では境界層が厚くなりやすく剥離する可能性があり、順圧力勾配(dp/dx < 0、加速場)では境界層が薄く保たれます。翼型や曲面体では Thwaites法、Head法、より高度には積分境界層法やCFDが必要になります。あくまで「平板の基準値」として活用してください。

実世界での応用

航空機の翼と摩擦抵抗の予測:翼面上の境界層が層流のまま維持できるか、いつ乱流に遷移するかは、燃費と直結する最重要設計項目です。設計者は層流翼(NLF翼)の形状最適化や、人工的に乱流遷移を促す乱流遷移トリップの配置で、抵抗をコントロールします。1/7乗則とBlasius解の比較は、こうした「層流維持のメリット」を定量化する第一歩です。

船舶の船体摩擦抵抗:船舶の総抵抗の50〜80%は船体表面の摩擦抵抗で、その大半が乱流境界層によるものです。船級協会の標準摩擦線(ITTC 1957 など)は、平板の乱流摩擦係数 Cf の実験式に船型補正を加えた形をしています。設計流速・喫水長から Re を計算し、模型試験から実船性能を推定する「Froude則」の出発点も、まさに本ツールの式です。

CFD(数値流体解析)の壁関数と検証:RANSやLESで境界層を解く際、計算メッシュの y+ 分布や壁関数の選択が結果を大きく左右します。教科書的な平板境界層の Cf, δ, τ_w を本ツールで計算し、CFD結果と比較することで、メッシュ品質や乱流モデル設定の妥当性を素早く検証できます。設計現場での「sanity check」に欠かせません。

熱伝達率の予測:強制対流の熱伝達係数 h は、Reynolds-Colburn類似性により摩擦係数 Cf と密接に関係します(St ≈ Cf/2, Pr ≈ 1)。乱流境界層では Cf が大きいため、層流より対流熱伝達が桁違いに向上します。熱交換器や冷却フィンの設計では、敢えて乱流を促す乱れ生成リブを取り付けることがあります。

よくある誤解と注意点

最も多い誤解は、「乱流のほうが境界層が厚いから、抵抗が小さいはず」と直感的に思い込むことです。実際は逆で、乱流境界層は厚くなった上に摩擦も大きくなります。なぜなら、乱流の渦が運動量を激しく混合するため、壁近くの速度勾配 ∂u/∂y が層流より急峻になり、せん断応力 τ_w = μ·(∂u/∂y) が増えるからです。本ツールで Cf の値を Blasius解と比べると、同じ Re_x で乱流のほうが数倍大きいことが一目で分かります。「厚い=抵抗大」が乱流境界層の本質です。

次に多いのが、境界層厚さ δ_99 を「物理的に明確な境界」と勘違いすることです。実際の流体には不連続な境界は存在せず、δ_99 は単に「u が主流速度 U の 99% になる便宜的な位置」です。定義によっては δ_95 や δ_90 が使われることもあり、それぞれ値が変わります。設計に使うのは、より物理的に明確な排除厚さ δ* や運動量厚さ θ のほうが信頼できます。形状係数 H = δ*/θ も、δ_99 より遷移や剥離の判定に有効です。

最後に、1/7乗則の式が Re_x = 0 や非常に小さい Re_x でも成立すると思い込む注意点です。本式は遷移後の完全乱流域 (Re_x ≥ 5×10⁵ 程度) で導出されており、層流域で適用すると δ や Cf が物理的にあり得ない値になります。例えば Re_x = 10⁴ で乱流式を使うと δ/x = 0.37·(10⁴)^(-0.2) ≈ 0.059、つまり境界層厚さが位置の6%という非現実的な大きさになります。本ツールでは Re_x < 5×10⁵ で警告を表示するようにしているので、必ず確認してください。実務では、層流域では Blasius、乱流域では本式、と切り分ける必要があります。