球阻力有哪些解析解？

黏性主导的慢速流动中，Stokes解 Cd = 24/Re 是精确解。稍微计入惯性的Oseen解为 Cd = (24/Re)(1 + 3Re/16)。更高雷诺数没有解析解，用Schiller-Naumann公式等经验式作参考值。

Stokes解能用到什么范围？

与经验式相比，雷诺数0.1时偏小3%，0.5时8.5%，1时13%。Oseen解到0.5附近吻合良好，但2时偏大11%，10时偏大66%。CFD验证一般雷诺数0.1以下用Stokes解，约0.1～1用Oseen解。

低雷诺数CFD验证要注意什么？

计算域大小。Stokes流中球的影响只与距离成反比缓慢衰减，外边界太近阻力就偏大。圆筒中的球，计算域为直径10倍时+26%，50倍时+4.4%，200倍时仍有+1.1%。要把计算域取得足够大，或加入壁面效应修正。

球的阻力（Stokes解、Oseen解）

分类：V&V·解析解 | 2026-10-01 改写

CAE visualization for sphere drag theory - technical simulation diagram

球的抗力（Stokes/Oseen解）

理论：慢速流中的球

概述

🙋

我要用CFD模拟微粒沉降，想先用球的阻力验证。和Stokes公式比较就行吗？

🎓

选得好，但有两点要注意。第一，Stokes解只在雷诺数很小时精确，接近1时惯性影响使偏差超过一成。不选择与验证雷诺数相符的参考值，就会把正确的CFD误判为“不吻合”。第二是计算域大小。慢速流中球扰动的范围很广，外边界近了阻力就会偏大。把握这两点，球阻力就是很好的验证问题。

阻力系数公式

$$ C_D^{Stokes} = \frac{24}{Re}, \qquad C_D^{Oseen} = \frac{24}{Re}\left(1 + \frac{3}{16}Re\right), \qquad C_D^{SN} = \frac{24}{Re}\left(1 + 0.15\,Re^{0.687}\right) $$

$Re = \rho U d/\mu$ 以球直径为基准。Stokes解是忽略惯性的精确解，Oseen解近似计入远场惯性，Schiller-Naumann公式（SN）是雷诺数约800以下的经验式。

$$ v_t = \frac{(\rho_p - \rho_f)\,g\,d^2}{18\mu} \ \ (\text{Stokes}), \qquad v_t = \sqrt{\frac{4(\rho_p - \rho_f)\,g\,d}{3\rho_f C_D}} $$

沉降颗粒的终端速度。Stokes式只在雷诺数小时正确，大颗粒要把阻力系数和雷诺数反复迭代求解。

Coffee Break 随笔闲谈

Stokes、Oseen与佯谬

1851年斯托克斯在研究摆的空气阻力时，推导出慢速流中球的阻力 3πμdU。但用同样近似求解圆柱时却不存在解（Stokes佯谬），也有人指出即使对球，远处惯性也不能忽略。1910年奥森提出近似计入远场惯性的方法，解决了这一矛盾。后来密立根用Stokes定律处理油滴下落速度测量了电子电荷，并获得诺贝尔奖。

算例

例1：解析解的适用范围

雷诺数	Stokes	Oseen	经验式（SN）	Stokes偏差	Oseen偏差
0.1	240	244	247	−3.0%	−1.2%
0.5	48.0	52.5	52.5	−8.5%	+0.1%
1	24.0	28.5	27.6	−13.0%	+3.3%
2	12.0	16.5	14.9	−19.5%	+10.8%
10	2.40	6.90	4.15	−42.2%	+66.2%
100	0.24	4.74	1.09	−78.0%	+334%

例2：计算域大小

慢速流中，位于圆筒中心的球，其阻力随圆筒直径 D 与球直径 d 之比增加多少（Haberman-Sayre壁面修正）：

计算域 D/d	5	10	20	50	100	200
阻力增加	+68%	+26%	+11.7%	+4.4%	+2.2%	+1.1%

🙋

计算域为直径10倍，阻力也偏大26%？

🎓

慢速流中，球周围的速度扰动只与距离成反比缓慢衰减。所以外边界即使很远，在那里固定速度也会约束流动，使阻力变大。按高雷诺数计算中“直径10～20倍”的习惯建计算域，低雷诺数时会偏差很大；要做到1%精度的验证，需要约直径200倍的计算域。也可以不扩大计算域，而用壁面修正公式修正参考值再比较。

例3：砂粒沉降速度

粒径	Stokes	经验式迭代	雷诺数	Stokes偏大
50 µm	2.25 mm/s	2.18 mm/s	0.11	3%
100 µm	9.00 mm/s	7.98 mm/s	0.80	13%
200 µm	36.0 mm/s	24.8 mm/s	4.96	45%
500 µm	225 mm/s	78.6 mm/s	39.2	187%

验证步骤

确定验证的雷诺数，选择相应参考值（Stokes、Oseen或经验式）。
把计算域取得足够大（低雷诺数时直径100倍以上），或加入壁面修正。
加密球表面附近网格，用三级以上确认网格依赖性。
把阻力分为压差阻力和摩擦阻力，与Stokes解的比例（1:2）比较。
改变雷诺数，确认与经验式吻合的趋势。

Coffee Break 随笔闲谈

“怎么算都大两成”

为验证微粒CFD，在雷诺数0.05下计算球阻力，网格再细也比Stokes解大约两成。问题不在网格，而是计算域只有球直径的12倍。慢速流中即使远处边界也有影响，把计算域扩大到直径200倍后差异约为1%。用高雷诺数流动的感觉建计算域，低雷诺数时会出现意想不到的误差。

常见错误

错误与对策

错误	后果	对策
Re 1时与Stokes解比较	把正确CFD判为错误	用Oseen或经验式
计算域太小	阻力偏大	直径100～200倍
大颗粒用Stokes沉降速度	严重高估速度	与阻力系数迭代
表面附近网格粗	摩擦阻力错误	解析边界层
只比较总阻力	漏掉误差相互抵消	分为压差与摩擦

🙋

我想试试相关计算。

🎓

可以用球阻力计算、沉降速度计算和阻力系数计算试算。其他流动验证问题参考圆管层流、库埃特流和顶盖驱动方腔流。