空気 20°C: ν ≈ 1.5×10⁻⁵ m²/s, ρ ≈ 1.2 kg/m³。水 20°C: ν ≈ 1.0×10⁻⁶ m²/s, ρ ≈ 998 kg/m³。
平板上面に沿った流れ／緑実線＝Blasius プロファイル u/U／青破線＝δ_99／橙＝δ*／紫＝θ
ゼロ圧力勾配の平板上層流境界層について、Blasius 相似解は次の代表値を与えます。
局所 Reynolds 数。U は主流速度、x は前縁からの距離、ν は動粘度：
$$\mathrm{Re}_x = \frac{U\,x}{\nu}$$
99% 境界層厚さ・排除厚さ・運動量厚さ：
$$\delta_{99} \approx \frac{5.0\,x}{\sqrt{\mathrm{Re}_x}},\quad \delta^{*} \approx \frac{1.721\,x}{\sqrt{\mathrm{Re}_x}},\quad \theta \approx \frac{0.664\,x}{\sqrt{\mathrm{Re}_x}}$$
局所摩擦係数と壁面せん断応力。ρ は流体密度：
$$C_f = \frac{0.664}{\sqrt{\mathrm{Re}_x}},\qquad \tau_w = \tfrac{1}{2}\,C_f\,\rho\,U^2$$
形状係数 H = δ*/θ ≈ 2.59。Re_x ≳ 5×10⁵ で層流→乱流の遷移が始まり、以降は Blasius の式は適用できません。
Blasius 境界層シミュレーターとは
🙋
「境界層」って物の表面の流れが遅くなる薄い層のことですよね。Blasius 解って何がすごいんですか？
🎓
ざっくり言うと、Prandtl が立てた境界層方程式を「平板・層流・圧力勾配ゼロ」という一番簡単な状況で初めて解いてみせた、ターボ機械や航空力学の出発点になる古典解だ。上のシミュレーターで「評価位置 x」を動かすと δ_99 が √x で増えていくのがわかる。これが Blasius の核心だよ。
🙋
δ_99 って 99% の意味ですか？なんで 100% にしないんでしょう。
🎓
いいところに気づいたね。境界層は数学的には主流に「漸近」していくので、厳密に「ここまでが境界層」という境目はないんだ。だから慣習的に「流速が主流の 99% になる高さ」を境界層厚さと呼ぶ。Blasius では δ_99 ≈ 5x/√Re_x。シミュレーターのデフォルト（U=1 m/s, x=100 mm, ν=1.5×10⁻⁵ m²/s）だと約 6.1 mm になるはずだ。
🙋
「排除厚さ」「運動量厚さ」もカードに出てますね。δ_99 とどう違うんですか？
🎓
δ_99 は「見た目の厚さ」、δ* と θ は「効果の厚さ」と覚えるといい。δ* は境界層で減速した分だけ主流が押しのけられる相当の厚さ——例えばノズルの中で見ると「実効断面が δ* だけ狭くなる」。θ は運動量損失を厚さ次元に直したもので、壁の摩擦の積分そのものだ。Blasius では δ*/θ = 2.59、これは層流境界層を判別する特徴的な値だよ。
🙋
C_f が 1/√Re_x で減るのは、Re が大きくなると摩擦が減るっていうことですか？
🎓
「無次元化した摩擦係数」が減る、という意味だね。実際の壁面せん断応力 τ_w は τ_w = ½ C_f ρ U² だから、流速 U を上げれば U² の効きが C_f の 1/√Re_x の効きを大きく上回って、絶対値としての摩擦は増える。シミュレーターで U を 1→5 m/s に上げてみて。C_f は約 √5 倍小さくなるけど、τ_w は約 11 倍に増える。実機の抗力を見るときはこの両方を意識する必要がある。
🙋
Re_x が 5×10⁵ を超えると赤いタグになりますね。そこから先は使えないんですか？
🎓
そう、平板上の自然遷移はだいたい Re_x ≈ 5×10⁵ で始まる。それより先は乱流境界層になって、厚さは √x ではなく x⁴ᐟ⁵ で増え、C_f もずっと大きくなる（典型的には 0.0592/Re_x^(1/5)）。実機の翼や船体は前縁から少しすると乱流に切り替わるので、設計では Blasius と乱流相関式を組み合わせて使うのが基本だ。シミュレーターでは Re_x が 5×10⁵ を超えると警告を出して、Blasius の適用外であることを知らせている。
よくある質問
乱流境界層と比べると、層流境界層はどう違いますか？
同じ位置 x で比べると、層流境界層は乱流より薄く、摩擦係数も小さいですが、剥離（流れがはがれる現象）に弱いという特徴があります。乱流境界層は混合が活発な分だけ運動量が壁付近まで運ばれて剥離に強く、ゴルフボールのディンプルや航空機の乱流促進子（タービュレータ）は意図的に層流→乱流遷移を早めて剥離を防ぐ仕組みです。設計では「摩擦を減らしたい局所→層流」「剥離を防ぎたい局所→乱流」と使い分けます。
この結果から、平板全体の抗力（摩擦抵抗）はどう計算しますか？
局所摩擦係数 C_f を 0 から板長 L まで積分（または平均）して、平均摩擦係数 C_F = 1.328/√Re_L を得ます。これに動圧と濡れ面積を掛けると摩擦抗力 D = C_F · ½ρU² · (b·L) になります（b は板幅）。Re_L が 5×10⁵ を超えるなら、層流域と乱流域を分けて積分する「混合境界層モデル」を使うのが実務的です。シミュレーターは局所値のみを出すので、抗力評価には別途積分が必要です。
CFD で Blasius と比較するときは何に注意すべきですか？
第一に、壁近傍の格子解像度です。粘性下層を解像するには最初のセルの y+ を 1 以下にし、δ_99 の中に 10〜20 セルが入る厚さで揃えるのが目安です。第二に、入口条件と前縁形状。Blasius は薄い前縁を仮定するため、計算領域の入口を前縁に置くと圧力勾配が立って合いません。前縁を計算領域内に取り、滑り境界の助走区間を設けると Blasius と良く合います。第三に、層流計算を指定すること。標準で k-ε などをオンにすると合いません。
圧力勾配がある場合は使えますか？
使えません。Blasius は dP/dx = 0 を仮定した特解です。一般の楔形流れ（主流速度 U ∝ x^m）には Falkner–Skan 解が、もっと複雑な圧力勾配には積分形の境界層方程式（Thwaites の方法など）を使います。逆圧力勾配（dP/dx > 0）では層流境界層がすぐ剥離するため、実機ではこの領域がしばしば乱流遷移や流れ剥離の発端になります。Blasius 解は「圧力勾配ゼロの基準ケース」として、こうした拡張理論と比較する出発点になります。
実世界での応用
翼・船体・タービン翼の摩擦抵抗予測： 翼前縁から少しの区間は層流境界層が残るため、その範囲の摩擦抵抗を Blasius と平均化式 C_F = 1.328/√Re_L で推算します。乱流域は別の式（Schlichting 等）で計算し、合計して全体抗力を出すのが古典的な手順です。NACA の翼型設計データも、こうした層流／乱流の混合境界層モデルがベースになっています。
熱交換器の伝熱予測： 境界層と熱境界層は形が似ており（Pr = 1 でほぼ一致）、Blasius 解から局所 Nusselt 数 Nu_x ≈ 0.332 Re_x^(1/2) Pr^(1/3) が得られます。フィンや平板冷却の入口側で活躍する関係式で、シミュレーターで δ_99 が薄いほど熱境界層も薄く、熱伝達率が高いことに直結します。
CFD コードの検証ベンチマーク： Blasius は厳密な相似解を持つため、CFD ソフトの精度検証によく使われます。新しい乱流モデルや格子戦略を試すときの「層流ベースライン」として、δ_99 と C_f が解析解と一致するかをまず確認します。本シミュレーターの値を CFD ポスト処理の比較対象として使うことができます。
マイクロ流路・微小翼の設計： マイクロ流体デバイスや MAV（マイクロ航空機）の翼は Re が小さく、境界層全体が層流のままです。この場合 Blasius がほぼそのまま使え、流路の δ* で実効断面の補正、翼の C_F で抗力評価ができます。マクロな航空機より層流域の評価が支配的になるため、Blasius の理解が直結する分野です。
よくある誤解と注意点
最も多い誤解は、境界層厚さが「主流速度に依存しない」と考えてしまう ことです。式 δ_99 = 5x/√Re_x を見ると x が前面に出ていますが、Re_x = Ux/ν に U が含まれているため、結局 δ_99 = 5√(νx/U) で U の平方根に反比例 します。流速を 4 倍にすると境界層厚さは半分になります。シミュレーターで流速スライダーを 1 m/s から 4 m/s に動かしてみてください。δ_99 が約 6 mm から 3 mm に減るのが確認できます。
次に多いのが、排除厚さ $\delta^{*}$ と境界層厚さ $\delta_{99}$ を混同する ことです。Blasius では $\delta^{*}/\delta_{99} \approx 0.344$ で、$\delta^{*}$ は $\delta_{99}$ の約 1/3 しかありません。ノズル断面の実効補正をするときに $\delta_{99}$ を使うと過剰補正となり、流量や圧力損失の予測が大きくずれます。「内部流れの実効断面 → $\delta^{*}$」「外部流れの可視化的な厚さ → $\delta_{99}$」「抗力評価 → $\theta$」と用途で使い分けるのが基本です。シミュレーターでは 3 つを別色で同時表示しているので、その大小関係を視覚的に確認してください。
最後に、Blasius 解が「常に成り立つ」と思って遷移後の領域に適用してしまう ことです。Re_x が 5×10⁵ を超えると境界層は乱流に移行し、厚さは x⁴ᐟ⁵ で増え、C_f は約 0.0592/Re_x^(1/5) のオーダーになります。たとえば長さ 1 m の平板を 10 m/s の風が流れるとき、後端の Re_L ≈ 6.7×10⁵ なので、Blasius を全長に適用すると摩擦を大幅に過小評価します。シミュレーターでは Re_x ≥ 5×10⁵ で警告タグを表示するので、適用外領域に入っていないか必ず確認してから設計に使ってください。