球体阻力与Stokes/Oseen解析解的对比——故障排查指南

球体阻力（Stokes/Oseen解） — 故障排除指南

作为验证基准的球体阻力

均匀来流中球体所受的阻力，是CFD精度验证中使用最多的解析解基准。作为参照的理论有三个层次。

其中 \( Re = \rho U d / \mu \)（以直径为特征长度），\( C_D = F_D / (\frac{1}{2}\rho U^2 \cdot \frac{\pi d^2}{4}) \)（以投影面积为参考面积）。本页按出现频率从高到低，逐条排查"用CFD算球体却对不上理论"的原因。

症状0：定义搞错（先花3分钟确认）

• Re的特征长度——是直径还是半径。用半径算Re会与理论值系统性地差2倍。文献、软件、自己的后处理三者是否统一
• 参考面积——\( C_D \) 以投影面积 \( \pi d^2/4 \) 为准。若除以表面积 \( \pi d^2 \) 就差4倍
• 阻力的分量——是否为压差阻力＋摩擦阻力之和。低Re下摩擦占三分之二（Stokes解中压差占1/3、摩擦占2/3），只取其中一项会偏离很大
• 速度的基准——是来流速度，还是阻塞修正后的速度（与风洞对比时）

凡是偏差呈"2倍、一半、三成"这类整齐比例的，既不是物理问题也不是数值问题，基本都是定义搞错了。

症状1：低Re下阻力系数大于理论值——计算域尺寸的阻塞效应

由于Stokes流的扰动衰减很慢（\( u' \sim 1/r \)），圆柱形或方形计算域的侧边界会充当"看不见的墙"而抬高阻力。经典的壁面修正（管内球体的Bohlin/Haberman修正）中，修正量以球径与管径之比 \( \lambda = d/D \) 按 \( 1/(1 - 2.104\lambda + \dots) \) 的形式起作用，即使 \( \lambda = 0.02 \)（计算域为球径的50倍）也仍有约4%的增量。实务上有三条路可走。

1. 把计算域取得足够大——目标精度1%的话，径向至少取球径的100倍为宜。低Re基准的计算域"大到觉得过分"才是正确的
2. 用壁面修正式修正理论侧——与有限计算域的理论值（带壁面修正的Stokes阻力）比较。这是回避扩大计算域计算量的聪明做法
3. 做计算域尺寸研究——在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下与标准曲线比较（对对流项与尾流的验证）——这样可以考察代码的不同侧面。关联式本身也有±百分之几的离散，因此Re>1时的合格判定以±5～10%为宜。

症状3：Re一调高就不收敛、答案来回晃——流动结构变了

"Re=300下定常分析不收敛"并非求解器出了毛病，而是定常解本就不存在这一物理事实给出的答案（收敛诊断的通则）。另外，在Re=250下用轴对称分析而出现"收敛了却对不上实验"时，是对称性的强制造出了非物理的解。请先确认Re区间与分析类型的对应表。

症状4：网格加密了也不收敛、数值漂移

• 球面附近的分辨率——阻力是表面压力与剪切的积分。边界层（低Re时是球径量级的黏性层）内要有足够的层数，球面本身也要离散得光滑（粗糙的多面体近似会把表面积与曲率误差直接带进结果）
• 尾流的分辨率——Re>20时分离区分辨率不足会扭曲压差阻力。别忘了在尾流方向加密
• 收敛判据——低Re加上大计算域，收敛很慢。要一直算到阻力的迭代历程完全平坦为止（不要只看残差就收工）
• 边界条件——侧边界宜用远场入流/开放边界而非滑移壁（滑移壁同样残留阻塞）。出口要离球足够远的下游，同时注意低Re下上游侧的影响同样大——Stokes流前后近乎对称

推荐验证流程（按此顺序，一天完成）

1. 定义检查——把Re、参考面积、力的分量的定义写成文字（症状0的检查项）
2. Re=0.05、计算域100D——以Stokes解±1～2%为目标。若偏离，用计算域研究（25/50/100D）把阻塞单独分离出来
3. 网格收敛——三套网格做GCI（步骤）。阻力是积分量，收敛得相对较好
4. Re=100——标准阻力曲线±5%。确认尾流加密的效果
5. 留档——把计算域尺寸、网格、收敛判据、比较基准（公式）记入验证台账，作为后续项目的标准设置

相关文章：解析解基准一览、GCI故障排查、求解器收敛的确认与诊断。