球周围的流动
球周围流动的理论基础
概述
老师，球周围的流动与圆柱有什么区别？
由于球是三维物体，其尾流结构本质上是三维的。不像圆柱那样存在可以用二维计算的区域。另外，称为阻力危机（drag crisis）的、依赖于Re数的 $C_D$ 急速变化在工程实用中非常重要。高尔夫球的凹陷就利用了这一现象。
Stokes 流（极低Re数）
请从最基础的地方开始讲解。
当Re $\ll 1$ 时，惯性项可以忽略，Stokes方程存在精确解。
这就是Stokes阻力定律。$a$ 是球的半径，$U_\infty$ 是均匀来流速度。写成阻力系数的形式：
$C_D = 24/Re$，当Re增大时阻力系数下降了？
但是阻力的绝对值 $F_D$ 与 $U_\infty$ 成正比增加。$C_D$ 下降是因为惯性力的参考量 $\frac{1}{2}\rho U_\infty^2$ 增长更快。
Re 数的流动过渡
让我整理球周围流动的分类。
| Re 范围 | 流动状态 | 特征 |
|---|---|---|
| Re < 20 | 稳定轴对称 | 无分离~微小回流 |
| 20 < Re < 210 | 稳定轴对称尾流 | 环状回流区增长 |
| 210 < Re < 270 | 稳定非轴对称 | 平面对称性丧失（规则分岔） |
| 270 < Re < 800 | 周期性涡脱落 | 发卡形涡的规则脱落 |
| 800 < Re < $3 \times 10^5$ | 亚临界区 | 湍流尾流，$C_D \approx 0.44$ |
| $3 \times 10^5$ < Re < $4 \times 10^5$ | 阻力危机 | $C_D$ 从 $0.44$ 急降到 $0.1$ |
| Re > $4 \times 10^5$ | 超临界区 | 湍流边界层分离 |
与圆柱相比，涡脱落开始的Re数较高呢。
正是如此。由于球是三维的，流体能够绕过物体，因此后流不稳定性延迟了。此外，与圆柱的涡街不同，球的涡脱落是发卡形涡的脱落，这是一种更加复杂的结构。
经验阻力相关式
中等Re数下 $C_D$ 怎样估算？
通常采用Schiller-Naumann相关式。
更宽范围的方程式是Morrison方程。
该公式从 $Re = 10^{-1}$ 到 $10^6$ 范围内连续拟合，包含阻力危机。
Morrison方程相当复杂，但能表达阻力危机，看起来很方便。
高尔夫球凹陷使飞行距离增加一倍的悖论
把表面光滑的球和布满凹陷的球以相同速度抛出，哪个飞得更远——人们直觉上会认为"光滑球空气阻力小"，但实际结果相反。凹陷通过使边界层湍流化，将分离点向后移动，从而显著缩小尾流宽度。结果是凹陷球的阻力系数约为光滑球的一半。一般的高尔夫球飞行距离据说是没有凹陷情况下的约2倍。这是球周围流动中的"阻力危机"现象，通过人工设计在较低Re数下提前触发，是流体力学上非常优雅的解决方案。在用CFD求解球的阻力系数时，表面粗糙度模型的影响直接关系到这一点。
球周围流动的数值计算方法
数值方法的选择
计算球周围流动适合用什么方法？
由于球本质上是三维问题，与圆柱相比计算成本高得多。
| Re 范围 | 推荐方法 | 网格规模参考 |
|---|---|---|
| Re < 300 | DNS | 50万~200万单元 |
| 300 < Re < $10^4$ | LES | 500万~5000万单元 |
| $10^4$ < Re | RANS (SST) / DES | 200万~2000万单元 |
网格策略
球的网格怎样生成？
理想做法是在球面周围使用O型（球壳型）网格。
- 壁面第一层：$y^+ < 1$（壁面解析）。Re = 1000时 $\Delta r / D \approx 5 \times 10^{-3}$
- 球面方向分割：赤道附近细化（追踪分离点移动），极点附近可粗糙
- 尾流区域：球心下游延伸$30D$以上。追踪发卡形涡的发展
- 外边界距离：球心到外边界$20D$以上（阻塞比<0.25%）
极点特异性（网格汇聚到一点的问题）怎样处理？
这是个好问题。标准的球面坐标系网格在极点处单元退化。对策有：
- 立方球面：将立方体各面映射到球面。不存在极点特异性
- 非结构网格：球面附近棱柱层，外部四面体/多面体
- 重叠网格：用另一个网格块覆盖极点附近区域
工程实践中，STAR-CCM+的多面体网格或Fluent的poly-hexcore很方便。自动在球面生成棱柱层，外部用多面体填充。
轴对称计算的应用
Re较低时用轴对称计算可以吗？
但Re > 210后轴对称性被破坏，必须进行完整三维计算。即使"轴对称计算效果很好"，也不能外推到高Re。这是很危险的。
粒子追踪耦合
球周围流动也与粒子沉降有关吧。
对。沉降球的终端速度由Stokes定律给出：
"球周围直接用球形网格会失败"的原因
不少工程师在模拟球周围流动时会想到"用球形外边界域"，看似整洁但往往失败。原因是超椭球或球形外边界在远场的流向不清楚，边界条件设定困难。标准做法是"将球放在立方体域中心，仅球面高分辨率网格"。这样入口、出口、侧面边界条件明确，结构化或多面体网格都容易处理。流向网格延伸到球直径的30~40倍，横向15~20倍，不然会产生壁面干扰。看似整洁的设计往往不如明确的边界条件来得实用。这是工程经验的一个教训。
球周围流动的工程应用
分析方法
用CFD分析球周围流动，步骤如何？
分步骤讲解。
1. 确认Re数并选择方法：计算 $Re = \rho U D / \mu$，选择DNS/LES/RANS
2. 利用对称性：Re < 210可用轴对称计算。210 < Re < 270时某些情况可利用平面对称性将计算域减半
3. 网格生成：球面棱柱层加外部非结构网格
4. 边界条件：入口均匀流，出口压力固定，球面无滑移，远场滑移
5. 初始化和计算：先低Re稳定计算，再转非定常
6. 统计积累：舍去初始瞬态后开始统计平均
验证用基准
有什么数据可以验证结果的正确性？
可以与以下数据对比。
| Re | $C_D$ | 尾流长度 $L_w/D$ | 分离角 | 出处 |
|---|---|---|---|---|
| 50 | 1.57 | 1.06 | 130° | Johnson & Patel (1999) |
| 100 | 1.09 | 1.70 | 127° | Johnson & Patel (1999) |
| 200 | 0.77 | 2.2 | 116° | Tomboulides & Orszag (2000) |
| 300 | 0.66 | — | — | DNS各种 |
| $10^4$ | 0.41 | — | $\sim 82°$ | 实验数据 |
$C_D = 24/Re$ 的Stokes解与Re=100时的 $C_D = 1.09$ 相差很大。Stokes的话是0.24。
Stokes定律仅在Re $\ll 1$ 时有效。即Re = 1时也需要Oseen修正（$C_D = (24/Re)(1 + 3Re/16)$）。
壁面压力分布的确认
除了 $C_D$，还要检查什么量？
将球面压力系数 $C_p(\theta)$ 作为角度函数绘图。在驻点（$\theta = 0$）处 $C_p = 1$，分离点附近急速变化。这个分布与基准相符时，说明不只积分量（$C_D$）正确，局部流场结构也对。
工程应用技巧
现场需要注意什么？
几个重要点。
- 湍流强度的影响：入口湍流强度高时阻力危机在较低Re出现。需根据风洞自由流湍流强度（TI = 0.1%~5%）调整
- 壁面粗糙度：如高尔夫球凹陷，表面粗糙度改变过渡位置，对 $C_D$ 影响很大。RANS壁函数中可设置粗糙度参数 $k_s$
- 非定常力评价：Re > 270时涡脱落产生横向力。需评估不仅是时均 $C_D$，还有 $C_L$ 的RMS值
- 支撑结构影响：实验中用支杆支撑球。与CFD比较时要注意支杆有无
雨滴与球的形状——"下落时变椭圆形"是都市传说
"雨粒是泪滴形"这是广为人知的误解。实际上小雨滴（直径1mm以下）几乎是球形，较大的变成汉堡包形（扁平球体）。这种形状由球周围流体动力和表面张力的平衡决定。雨滴的终端速度分析以Stokes流（低Re球阻力）为出发点，大雨滴时Re数增高，非线性效应开始起作用。用CFD模拟"下落变形液滴"需与VOF方法结合，但首先必须准确重现"刚体球的阻力"作为验证步骤。身边的雨粒物理与球流动理论直接关联。
球周围流动的软件比较
工具特性对比
适合球周围流动的CFD工具有哪些？
从球分析专用功能来比较。
| 工具 | 球面网格生成 | 粒子追踪 | 阻力危机处理 |
|---|---|---|---|
| Ansys Fluent | Inflation Layer自动生成、Poly-hexcore | DPM模块 | Transition SST支持 |
| STAR-CCM+ | 多面体+棱柱层 | 拉格朗日多相 | k-kl-omega遷移模型 |
| OpenFOAM | snappyHexMesh | DPMFoam, MPPICFoam | 手工实现（Langtry-Menter等） |
| COMSOL | Physics-controlled meshing | Particle Tracing | SST + γ-Re_θ |
Ansys Fluent 设置示例
用Fluent求解球周围流动的推荐设置是什么？
以亚临界区（Re = $10^4$）的RANS计算为例：
- 网格：Poly-hexcore。球面Inflation Layer（15层、增长率1.2、初层 $y^+ \approx 1$）
- 求解器：Pressure-Based, Transient, Coupled
- 湍流模型：SST $k$-$\omega$。要捕捉过渡可用Transition SST ($\gamma$-$Re_\theta$)
- 空间离散：Least Squares Cell Based (gradient), Second Order Upwind (momentum), Second Order (pressure)
- 时间步长：$\Delta t = 0.01 D / U_\infty$ 左右
- 监测：球表面 $C_D$, $C_L$ 每步输出
OpenFOAM 设置示例
OpenFOAM呢？
Re = 100稳定流解析用simpleFoam足够。
```
constant/transportProperties: nu 0.01; (U=1, D=1 时Re=100)
system/fvSchemes:
divSchemes: div(phi,U) Gauss linearUpwind grad(U);
system/fvSolution:
SIMPLE: nNonOrthogonalCorrectors 1;
relaxationFactors: U 0.7; p 0.3;
```
非定常（Re = 300）则改用pimpleFoam，ddt改为backward（二阶隐式）。不需要湍流模型，直接计算（laminar）就行。
粒子沉降分析的工具选择
粒子沉降这类应用，用哪个工具好？
视用途而定。
- 单个粒子详细分析：Resolved CFD（球体网格化）。OpenFOAM的
pimpleFoam足够 - 多粒子分散相分析：Unresolved CFD-DEM。FluentDPM或OpenFOAM的
DPMFoam - 浓厚浆液：Euler-Euler二流体模型。Fluent的Mixture/Eulerian或STAR-CCM+的Eulerian Multiphase
- 流化床：CFD-DEM耦合。CFDEM coupling（OpenFOAM + LIGGGHTS）很有名
CFDEM可以免费用吗？
开源的。OpenFOAM和粒子模拟器LIGGGHTS耦合。到数十万粒子规模计算时间还能接受。
药物粉末吸入器——球形微粒的空气动力学是生死攸关
哮喘治疗用的粉末吸入器（DPI）中，药剂需以球形微粒（直径1~5μm）的形态深入肺部。太大会卡在气道，太小会被呼出不起作用。粒子在气道内的运动由球周围Stokes阻力与重力、惯性竞争决定，用粒子追踪的CFD（DPM分析）是装置设计的核心技术。各大药企（葛兰素史克、诺华等）在吸入器的CFD模拟上投入巨资，球的阻力系数精度直接关系到"药能否有效"。球周围流动变成了医疗器械性能基准，这是少有人讲的重要话题。
球周围流动的前沿研究
尾流的分岔结构
球的尾流过渡比圆柱更复杂吧。
Johnson & Patel (1999)之后，球的尾流过渡被深入研究。从轴对称到非轴对称的过渡是规则分岔（叉形），Re $\approx 210$时发生。此时尾流有一个对称面，呈"单翼"结构。
继续增大Re到$\approx 270$，Hopf分岔出现，周期性发卡形涡脱落开始。涡脱落方位被Re = 210时选定的对称面束缚。
与圆柱的卡门涡街交替脱落不同啊。
球的情况是单侧连续脱落发卡涡。这由于后流对称性不同——圆柱是平面对称，球是轴对称——导致的。
阻力危机的机制
阻力危机的CFD模拟可能吗？
阻力危机由边界层层流→湍流过渡将分离点后移引起。CFD要捉住这一点必须用过渡模型。
- $\gamma$-$Re_\theta$过渡模型（Langtry & Menter）：用RANS预测过渡。Fluent、CFX、STAR-CCM+都支持
- LES：格子足够细则自然捕捉过渡，但成本巨大
- DNS：Re ~ $10^4$为止已有实现（Rodriguez等2011）
Achenbach (1972)的实验表明光滑球的阻力危机在Re $\approx 3.7 \times 10^5$发生。这是自由流湍流强度为$0.45\%$的情况。
旋转球与 Magnus 效应
球旋转时怎样？
旋转球受Magnus力。用旋转参数 $\alpha = \omega a / U_\infty$ 无量纲化，横向力系数 $C_L$ 为：
足球的无转飘球（butterfly effect）是无旋转（$\alpha \approx 0$）下球的剥离点位置随机变动，横向力不稳定振荡的现象。这在阻力危机Re域最明显。
浸没边界方法
涉及球运动的问题可用浸没边界法吧。
IBM在固定直交格子上用罚函数或内插表示物体。球的沉降、碰撞、多体问题都适用。
- Direct Forcing IBM：Uhlmann (2005)方法。粒子系大规模DNS用的标准
- 体积罚函数：物体内速度强制为零。实现简单
- 切割单元法：物体与格子交界处离散化保守性处理
多个球相互作用的问题每步都重建网格不现实吧。
太空垃圾（碎片）再入与Stokes阻力
坠落地球的卫星碎片在大气层中能烧尽到什么程度的计算中，各种碎片形状的阻力系数很关键。高空时大气密度薄，Kn数（Knudsen数）很大，连续体流体力学不适用。但高度80km以下，连续体球周围流动理论可用。JAXA及各国航天机构在计算"碎片撞击地面的概率"时，用球、圆柱、平板等基本形状的阻力系数组合模拟碎片群的轨迹。"准确知道球的阻力系数"背后是航天安全保障的一角，这是鲜为人知的事实。
球周围流动的问题排查
常见问题
球周围流动计算常见失败原因有哪些？
1. 轴对称计算的后流不非对称
原因：即使Re > 210，如果网格对称性太强或初始条件完全轴对称，轴对称解仍可能输出。
对策：改为完整三维计算，在初始条件加入微弱的非对称扰动。或者如果用的是楔形二维网格，改为360度完整网格。
2. Stokes解的不匹配（低Re数）
Re = 0.1算出来与 $C_D = 24/Re$ 不符。
检查要点：
- 计算域大小：低Re数Stokes流的影响范围很广。外边界离球不足$100D$时，会受边界影响（Stokes流的影响以$1/r$衰减，极其深远）
- 入口、侧面条件：入口太近会影响球上游
- Oseen修正：即Re = 0.1也有几%的偏差。$C_D = (24/Re)(1 + 3Re/16) = 240.45$，比$C_D = 240$相差0.2%
3. 阻力危机没被捕捉
纯RANS SST $k$-$\omega$中，湍流模型的假设使边界层从一开始就是湍流。不会再现层流剥离→过渡→湍流再附着的过程，故不出现阻力危机。
对策：启用$\gamma$-$Re_\theta$过渡模型。Fluent中Models → Viscous → Transition SST (4-equation) 选择。
4. 尾流涡结构崩坏
发卡涡出不来。
检查要点：
| 症状 | 原因 | 对策 |
|---|---|---|
| 涡在下游消失 | 数值扩散 | 使用二阶及以上精度格式。提高中心差分混合比例 |
| 涡形状不自然 | 网格各向异性 | 尾流域宽高比控制在3以下。等向网格最优 |
| 涡脱落不规律 | 统计时间短 | 至少统计50个周期以上 |
| 3D模式不出现 | 网格方位角分辨率不足 | 方位角最少64分割（5.6度间隔）以上 |
5. 计算发散（高Re数LES）
高ReLES中用中心差分在格子Re数大的区域容易震荡、发散。
对策：
- 混合格式（中心差分+少量风上）。Fluent的Bounded Central Differencing就是这种
- 利用SGS模型扩散。Dynamic Smagorinsky按局部自适应调整
- 细化网格降低格子Re数（根本但成本高）
问题出现时先从网格和边界条件怀疑是对的。
对。最稳妥的方法是从低Re验证过的情况开始，逐步增加Re。Re=100时验证$C_D = 1.09$、$L_w/D = 1.7$，再往上才有把握。
"CD与实验差3成"——过渡模型前先确认这些
CFD算出的球阻力系数比实验高30%，这种事时常发生。常见两个原因是"计算域太窄的阻塞效应"和"表面网格粗导致球面解析度不足"。域阻塞率（球截面/域截面）超5%时壁面加速流动导致阻力增大。这点实验也一样，风洞试验标准做阻塞补正。CFD应让阻塞率在1%以下。这样还不符再怀疑湍流或过渡模型精度的调查顺序很重要。外观整洁不如边界条件明确。
相关主题
更多细节
错误