境界層 (Boundary Layer) — CAE用語解説

ざっくり言うと、壁面のすぐ近くで流体の速度がゼロ（壁に張り付いた状態）から主流の速度まで急激に変化する薄い層のことだ。この層の中では粘性の影響が非常に大きくて、壁面の抗力や熱伝達のほとんどがここで決まる。例えば飛行機の翼の表面摩擦抗力は、まさに境界層内の速度勾配で生じるんだよ。

それがすべりなし条件（no-slip condition）で、粘性流体の基本原理だ。分子レベルで壁面に張り付いた流体が隣の層を引っ張る——その連鎖で壁からの距離とともに速度が回復していく。その回復が完了するまでの領域が境界層。プラントル（Prandtl）が1904年に提唱した概念で、流体力学の教科書には必ず出てくる最重要テーマの一つだよ。

大きく違うよ。層流境界層は流体が整然と層状に流れていて、速度分布はなめらかな放物線に近い形になる。厚みの成長が遅いから壁面摩擦も小さい。一方、乱流境界層は壁面のごく近く（粘性底層）を除いて乱流混合が激しく、速度分布がずっと「ふっくら」している——壁に近いところで急激に速度が立ち上がるんだ。

そのとおり、乱流の方が摩擦抗力は大きい。でも工業的な流れはReynolds数が高いから、ほとんどの場合で乱流境界層を扱うことになる。例えば自動車の表面だとRe=10⁶〜10⁷のオーダーで、前縁のごく短い距離を除けばほぼ全域が乱流境界層だ。逆にマイクロ流路やMEMSデバイスだとRe=1〜100程度で層流のまま、ということもある。

y⁺は壁面からの距離を摩擦速度で無次元化したものだ。定義はこうなる：

$$y^+ = \frac{y\,u_\tau}{\nu}, \qquad u_\tau = \sqrt{\frac{\tau_w}{\rho}}$$

$y$ は壁面からの距離、$u_\tau$ は摩擦速度、$\tau_w$ は壁面せん断応力、$\rho$ は密度、$\nu$ は動粘性係数だ。CFDでは「第一層セルの中心が $y^+$ いくつの位置にあるか」が、乱流モデルの選択とメッシュ設計を直接左右する最重要パラメータだよ。

それは使う乱流モデルによって変わるんだ。整理するとこんな感じ：

• 壁関数を使う場合（Standard k-ε等）: $y^+ = 30\text{--}300$。粘性底層を飛ばすから、逆に $y^+$ が小さすぎると壁関数の仮定が崩れてエラーの原因になる。
• 壁面分解する場合（SST k-ω、Spalart-Allmaras等）: $y^+ \approx 1$。粘性底層の中にセルを配置して直接解く。

実務で一番よく使われるSST k-ωモデルは $y^+ \approx 1$ が推奨だけど、壁関数との自動切替機能を持つソルバーも多いから、$y^+ = 1\text{--}5$ あたりに入っていれば大抵は大丈夫だよ。

半分正解で半分危険な理解だね。壁関数は対数則を使って粘性底層を「スキップ」する手法だから、第一層セルを対数則領域（$y^+ = 30\text{--}300$）に配置することで計算コストを大幅に下げられる。でも「粗くていい」のは壁法線方向の第一層だけで、壁に沿った方向（流れ方向や周方向）はちゃんと解像しないと意味がない。

あるよ、結構多い。壁関数の前提は「壁面近くの流れが対数則に従う」ことだけど、次のような場合はこの仮定が崩れるから壁面分解が必要になる：

• 流れの剥離・再付着：逆圧力勾配が強い領域では対数則が成り立たない。自動車のリアウィンドウや翼の失速領域が典型例。
• 強い加速・減速：ノズル内の加速流やディフューザの減速流。
• 熱伝達を精度よく求めたい場合：電子機器冷却やガスタービン翼の温度予測など、壁面近傍の温度勾配が重要な問題。
• 遷移を含む流れ：層流→乱流の遷移を捉えたい場合。

こういう場面では素直に $y^+ \approx 1$ のメッシュを作るべきだ。計算コストは上がるけど、壁関数で誤った結果を出すよりはずっとマシだよ。

まず第一層厚さ $\Delta s$ は、目標の $y^+$ から逆算する。流れの条件からまず壁面せん断応力を推定する必要があるんだけど、平板の近似式を使うと実用的に計算できる：

$$C_f \approx \frac{0.058}{Re_L^{0.2}}$$ $$\tau_w = \frac{1}{2}\,C_f\,\rho\,U_\infty^2$$ $$\Delta s = \frac{y^+\,\mu}{\sqrt{\rho\,\tau_w}}$$

例えば空気（$\rho=1.225$ kg/m³、$\mu=1.8\times10^{-5}$ Pa·s）で $U_\infty=30$ m/s、代表長さ $L=1$ m、$y^+=1$ を狙うと、$\Delta s \approx 0.015$ mm くらいになる。

適当にやってるエンジニア、実は多いけど（笑）、ちゃんと指針がある。成長率（Growth Rate）は 1.1〜1.2 が推奨。1.3以上にするとセル間の品質が急激に悪化して、数値拡散が増えたり収束が悪くなる。層数は、インフレーション層の最外層が境界層厚さ $\delta$ を超えるまで積むのが理想で、だいたい 15〜25層 が目安だ。

ほぼすべてのCFDソルバーが壁面の $y^+$ をポスト処理で出力できるよ。OpenFOAMなら yPlus ユーティリティ、FluentならContour Plotで Wall Y+ を選べばいい。重要なのは、注目する壁面全体で $y^+$ が目標範囲に入っていることを確認すること。局所的に $y^+$ が大きくはずれている箇所があれば、そこだけメッシュを細かくする必要がある。解析結果を信頼する前に、必ず $y^+$ の分布を確認するのが鉄則だよ。