球の抗力（Stokes/Oseen解） — トラブルシューティングガイド
より充実した内容を sphere-drag.html でご覧いただけます。
検証ベンチマークとしての球の抗力
一様流中の球の抗力は、CFDの精度検証に最もよく使われる解析解ベンチマークです。基準となる理論は3層あります。
| 理論 | 抗力係数 | 適用範囲 |
|---|---|---|
| Stokes解（慣性完全無視） | \( C_D = \dfrac{24}{Re} \) | \( Re \lesssim 0.1 \) |
| Oseen補正（慣性の一次補正） | \( C_D = \dfrac{24}{Re}\left(1 + \dfrac{3}{16}Re\right) \) | \( Re \lesssim 1 \) |
| 標準抗力曲線（実験相関：Schiller-Naumann等） | \( C_D = \dfrac{24}{Re}(1 + 0.15\,Re^{0.687}) \) | \( Re \lesssim 800 \) |
ここで \( Re = \rho U d / \mu \)（直径基準）、\( C_D = F_D / (\frac{1}{2}\rho U^2 \cdot \frac{\pi d^2}{4}) \)（投影面積基準）です。このページは「CFDで球を解いたら理論と合わない」ときの原因を、頻度順に潰していきます。
Re=0.05でCFDを回したのに、CDがStokes解より15%も大きいんです。Stokes解って厳密解ですよね？ CFDが間違ってるんでしょうか…。
その症状、犯人はほぼ確実に計算領域の狭さだよ。Stokes流の撹乱は \( 1/r \) でしか減衰しない——粘性流れの影響が驚くほど遠くまで届くんだ。だから領域境界（外側の壁や遠方境界）が球の50倍程度離れていても、壁の閉塞効果で抗力が数%上乗せされる。15%増なら境界はおそらく20D前後じゃないかな。低Reの球は「理論が厳密なのにCFDの設定要求が一番厳しい」ベンチマーク——そこがこの問題の教育的価値でもあるんだよ。
症状0: 定義の取り違え（最初に3分で確認）
- Reの長さスケール — 直径か半径か。半径でRe計算すると理論値と系統的に2倍ずれる。文献・ツール・自分の後処理で統一されているか
- 基準面積 — \( C_D \) は投影面積 \( \pi d^2/4 \) 基準。表面積 \( \pi d^2 \) で割ると4倍ずれる
- 抗力の成分 — 圧力抗力＋摩擦抗力の合計か。低Reでは摩擦が2/3を占める（Stokes解で圧力1/3・摩擦2/3）ので、片方だけ拾うと大きく外れる
- 速度の基準 — 一様流速か、閉塞補正後の速度か（風洞比較時）
「倍・半分・3割」のきれいな比率でずれているときは、物理でも数値でもなく、ほぼ定義ミスです。
症状1: 低ReでCDが理論より大きい——領域サイズの閉塞効果
Stokes流の長距離減衰（\( u' \sim 1/r \)）のため、円筒・角柱領域の側面境界は「見えない壁」として抗力を増やします。壁補正の古典（円管内の球のBohlin/Habermanの補正）では、球径と管径の比 \( \lambda = d/D \) に対し補正が \( 1/(1 - 2.104\lambda + \dots) \) の形で効き、\( \lambda = 0.02 \)（領域が球の50倍）でも約4%の増加になります。実務の選択肢は次の3つです。
- 領域を十分大きく — 目標精度1%なら半径方向に球径の100倍以上が目安。低Reベンチマークの領域は「大きすぎるくらいが正しい」
- 壁補正式で理論側を補正 — 有限領域の理論値（壁補正付きStokes抗力）と比較する。領域を広げる計算コストを回避できる賢い方法
- 領域サイズのスタディ — 25D/50D/100Dで抗力を計算し、\( 1/r \) 系の外挿で無限遠値を推定。閉塞が原因かどうかの診断を兼ねる
症状2: Re=1〜10で「Stokes解と合わない」
それは合わなくて正常です。Stokes解の誤差はRe=1で既に10%超、Re=10では数倍——比較すべき基準の選択ミスであり、CFDの欠陥ではありません。Re>0.1ではOseen補正、Re>1では標準抗力曲線（実験相関）を基準にします。逆に言うと、CFD検証としての使い分けは：①Re≪1でStokes解と比較（拡散項・境界処理の精密検証）、②Re=10〜100で標準曲線と比較（対流項・後流の検証）、の2点セットで行うと、コードの別の側面を検査できます。相関式自体にも±数%のばらつきがあるため、Re>1の合否判定は±5〜10%を目安にします。
症状3: Reを上げたら収束しない・答えが揺れる——流れ構造の変化
| Re域 | 流れの状態 | 解析上の帰結 |
|---|---|---|
| Re < 20程度 | 付着した定常軸対称流 | 定常・軸対称（2D）解析で可 |
| 20 < Re < 210程度 | 定常軸対称の剥離渦輪 | 後流の解像（球後方の細分）が精度を左右 |
| 210 < Re < 270程度 | 軸対称性の破れ（定常非軸対称） | 軸対称モデルは原理的に不可。3D化必須 |
| Re > 270程度 | 非定常渦放出 | 非定常3D解析。CDは時間平均で評価 |
「Re=300で定常解析が収束しない」は、ソルバーの不調ではなく定常解が存在しない物理の答えです（収束診断の一般則）。また軸対称解析をRe=250で使って「収束したが実験と合わない」場合、対称性の強制が非物理な解を作っています。Re域と解析タイプの対応表を最初に確認してください。
症状4: メッシュを細かくしても収束しない・値が漂う
- 球面近傍の解像 — 表面の圧力・せん断の積分が抗力。境界層（低Reなら球径オーダーの粘性層）に十分な層数と、球面の滑らかな離散化（粗い多面体近似は表面積・曲率誤差が直撃）
- 後流の解像 — Re>20では剥離域の解像不足が圧力抗力を歪める。後流方向の細分を忘れない
- 収束判定 — 低Re・大領域は収束が遅い。抗力の反復履歴が完全に平坦になるまで回す（残差だけで打ち切らない）
- 境界条件 — 側面は滑り壁でなく遠方流入/開放が望ましい（滑り壁でも閉塞は残る）。出口は球から十分下流（低Reでも上流側の影響が大きい点に注意——Stokes流は前後対称に近い）
推奨検証プロトコル（この順で1日）
- 定義チェック — Re・面積・力成分の定義を明文化（症状0のチェック項目）
- Re=0.05・領域100D — Stokes解±1〜2%を目標。外れたら領域スタディ（25/50/100D）で閉塞を単離
- メッシュ収束 — 3水準でGCI（手順）。抗力は積分量なので比較的よく収束する
- Re=100 — 標準抗力曲線±5%。後流細分の効果を確認
- 記録 — 領域サイズ・メッシュ・収束基準・比較基準（式）を検証台帳に残し、以後のプロジェクトの標準設定にする
球の抗力って地味な問題だと思っていましたが、ここまで奥があるんですね…。
この問題の価値は「厳密解があるのに、まともに合わせるには実務スキルの総動員が要る」ところなんだ。定義の規律・領域設計・メッシュ収束・流れ構造の知識・比較基準の選択——実プロジェクトで必要な検証スキルが全部この1問に入っている。新しいソルバーや新人の腕試しに球抗力を課すのは、そういう理由だよ。ここで±2%を出せる人は、実問題でも信頼できる検証ができる。
関連記事：解析解ベンチマークの一覧、GCIのトラブルシュート、ソルバー収束の確認と診断。
関連トピック
なった
詳しく
報告