Wall Distance (y+) Estimator
境界層を正しく捉える。解析解像度を科学的に定義する。
y+（ワイプラス）の物理：CFD解析の品質と境界条件を支配する壁面格子の科学
数値流体力学（CFD）のシミュレーション技術が飛躍的に発展した現代において、設計現場での活用は当たり前のものとなりました。しかし、その一方で「解析結果が現実と乖離している」「乱流モデルの選択が適切か分からない」といった技術的課題も散見されます。こうした問題の多くを解決する鍵を握るのが、y+（ワイプラス：無次元壁面距離）という物理量です。
流体が進む際、壁面（固体表面）との間には摩擦が生じます。この壁面近傍で発生する「境界層」という極めて薄い領域をいかに正確に捉えるかが、空力抵抗、熱伝達率、さらには剥離や失速といった複雑な現象の予測精度を決定します。y+ を正しく理解し、それに基づいたメッシュ生成（格子形成）を行うことは、CFD解析者が備えるべき最も基礎的でありながら最も高度な専門性の一つです。
図1：壁面からの距離と流速分布（粘性底層から対数領域まで）
1. y+ の物理的定義：なぜ「無次元」なのか
y+ は、壁面からの物理的な距離 y
を、流体の動粘性と「摩擦速度」を用いて無次元化した値です。なぜわざわざ無次元化する必要があるのでしょうか。それは、主流の速さや流体の種類が異なっても、境界層内の物理現象を共通のモノサシで測るためです。
- u* (摩擦速度)： 壁面におけるせん断応力 \( \tau_w \) を代表する速度スケール。 \( u* = \sqrt{\tau_w / \rho} \) で定義されます。
- ν (動粘性係数)： 流体の粘性を密度で割った値。分子レベルの拡散能力を示します。
- y： 壁面から最初のメッシュ（第一層）の重心までの物理距離。
この式から分かる通り、主流が速くなればなるほど（u* が増大する）、あるいは粘性が低くなるほど、同じ物理距離 \( y \) であっても y+ の値は大きくなります。つまり、高速流や低粘性流（高レイノルズ数流）を解析する際は、より薄いメッシュを用意しなければ、狙った y+ を維持できないことを意味します。
2. 境界層の三層構造と y+ の解釈
壁面近傍の流体の流れは、y+ の値に応じて劇的にその性質を変えます。一般的に、以下の3つの領域に分類されます。
① 粘性底層 (Viscous Sublayer)： y+ < 5
壁面に最も近い領域です。ここでは流体の粘性が支配的で、流速は壁からの距離に完全に比例します（\( u+ = y+ \)）。乱流の影響がほぼ無視できるほど穏やかな層です。翼の表面摩擦や、微細な熱伝達を厳密に解きたい場合は、この領域に第一層のメッシュを配置する（y+ ≈ 1）必要があります。これを「Low-Re（低レ）解析」や「境界解決型（Resolved）」と呼びます。
② 遷移領域 (Buffer Layer)： 5 < y+ < 30
粘性と乱流の効果が複雑に相殺し合う領域です。この領域は数学的な記述が難しく、CFDソルバー内での近似（壁面関数など）も誤差が出やすい「魔の領域」とされます。メッシュ作成時は、意図的にこの範囲に第一層が来ないように設計するのが熟練解析者のテクニックです。
③ 対数領域 (Log-law Region)： 30 < y+ < 300
壁面からある程度離れた領域で、流速分布が距離の対数関数に従う性質を持ちます。標準的な k-ε モデルなどに代表される「壁面関数（Wall Function）」を使用する場合、第一層はこの領域に配置します。境界層全体を細かく解く代わりに、この対数則を仮定して計算を簡略化することで、膨大なメッシュ数を節約し、計算コストを抑えることができます。
3. 解析目的に合わせたメッシュ戦略の選び方
全ての解析で y+ = 1 を目指すべきでしょうか？ 答えは「No」です。解析対象の規模と目的に応じて、以下の戦略を使い分けます。
戦略1：境界層解決 (Low Re Approach)
- 目標 y+： 1.0 以下（最大でも 5 未満）
- 適したモデル： k-ω SST, Spalart-Allmaras, 一部の k-ε 低レモデル
- 用途： 翼の剥離予測、エンジン内部の熱伝達、精密な空力解析
- 短所： 壁面付近のメッシュ数が膨大になり、計算時間が数倍〜数十倍に増える。
戦略2：壁面関数適用 (Wall Function / High Re Approach)
- 目標 y+： 30 〜 100（大きな解析対象なら 300 まで）
- 適したモデル： 標準 k-ε, 標準 k-ω
- 用途： 建築風環境、車全体の空力トレンド調査、空調解析
- 短所： 壁面付近の剥離や再付着といった局所的な現象の解像度は落ちる。
4. 推論ロジック：解析前に格子幅をどう決めるか
解析を始める前に第一層の厚み \( y \) を決めるために、本計算機は以下の工学的な経験式を利用しています。
- レイノルズ数の算出： 主流速 \( U \)、代表長さ \( L \)、動粘性 \( \nu \) から \( Re = UL/\nu \) を算出。
- 摩擦係数の推定： 平板の乱流境界層モデル（Schlichting の式など）を用いて、皮膚摩擦係数 \( C_f \) を推算。
- 壁面せん断応力の導出： \( \tau_w = 0.5 \cdot \rho \cdot U^2 \cdot C_f \)
- 摩擦速度の決定： \( u* = \sqrt{\tau_w / \rho} \)
- 物理距離の逆算： ターゲットの y+ を基に、 \( y = (y+ \cdot \nu) / u* \) を求める。
この予測値を用いてメッシュを作成し、解析を実行します。
5. CFD解析の落とし穴：y+ の「事後検証」が不可欠な理由
初心者が最も陥りやすい罠が、「予測計算した y+ だけで満足してしまうこと」です。本ツールを含め、解析前の予測はあくまで「平板の定常流」を仮定した近似に過ぎません。
実際の複雑な形状では、流速が急減する剥離域や、逆に加速する絞り部が存在します。解析が終了した後、必ず商用ソルバー（Ansys Fluent, Star-CCM+ 等）やオープンソース（OpenFOAM + ParaView）で「y-plus」の等値線図を確認してください。
もし、y+ = 1 を狙ったはずが、翼の後縁で y+ = 50 になっていたら、その領域の剥離予測は信用できません。より細かいメッシュに修正して再計算（リメッシュ）を行う勇気こそが、プロの解析者への第一歩です。
6. y+ に関するよくある質問 (FAQ)
Q. y+ が 10 になってしまいました。どうすればいいですか？
A. 最も避けるべき「遷移領域」に入っています。メッシュを大幅に細かくして 1 以下にするか、逆に少し粗くして（壁面関数を使う前提で）30 以上に振るのがセオリーです。
Q. 全ての壁で y+ を一定に保つべきですか？
A. 現実的には困難です。最も物理現象が重要な箇所（例：翼の前縁、熱交換器のフィン表面）でターゲット値を維持することを最優先し、影響の少ない箇所では多少の逸脱を許容するのが実務的な判断です。
Q. メッシュの「成長率」は y+ に関係ありますか？
A. 直接的な定義には含まれませんが、y+ = 1 を狙う場合、第一層からの拡大率（Growth Rate）を 1.1 〜 1.2 程度に抑えなければ、境界層全体の勾配を滑らかに捉えることができません。
まとめ
y+ は単なる無次元数ではありません。それは流体という気まぐれな存在と、物理モデルという数学的手段を繋ぐ「接点」です。適切な y+ 設定は、あなたの CFD 解析に客観的な根拠と高い説得力を与えてくれます。
本計算機を活用し、理論に基づいた高品質なシミュレーションを実現してください。Gojikara CFD は、エンジニアの皆様の正確な眼となるためのツールを提供し続けます。