体育用品的空气动力学分析
理论与物理
概述
老师，体育用品的空气动力学分析主要针对哪些对象呢？
高尔夫球、足球、网球、自行车头盔、滑雪跳台服等，空气动力学性能直接影响其表现的用品都是分析对象。
尤其是球类，其表面纹理（凹坑、面板接缝）会对阻力产生巨大影响。与光滑球体的$C_D \approx 0.47$相比，高尔夫球的凹坑能将$C_D$降低至约0.25，几乎减半。
表面的凹凸竟然能如此改变阻力吗？
控制方程与阻力危机
球体的空气阻力可用下式表示。
其中$A = \pi d^2/4$是球体的正面投影面积。
球体的阻力存在一个称为“阻力危机”的剧烈现象。当雷诺数超过某个临界值时，边界层会从层流转变为湍流，分离点向下游移动，尾流区缩小。结果导致$C_D$急剧下降。
高尔夫球的凹坑促进了湍流转捩，从而在较低的雷诺数下就引发了阻力危机，对吧？
没错。高尔夫球的初速约为70m/s，$Re \approx 2 \times 10^5$。对于光滑球体，此时仍处于高阻力区，但得益于凹坑，高尔夫球进入了低阻力区。这使得飞行距离增加了近一倍。
马格努斯效应
旋转的球体会受到马格努斯力的作用。
旋转参数：
其中$\omega$是角速度，$d$是球体直径。$S$越大，偏转越明显。
| 运动项目 | 典型的$S$值 | 效果 |
|---|---|---|
| 高尔夫（后旋） | 0.1--0.3 | 升力延长飞行距离 |
| 足球（弧线球） | 0.1--0.5 | 横向偏转（弧线） |
| 网球（上旋球） | 0.2--0.6 | 向下的力改变弹跳 |
| 棒球（滑球） | 0.1--0.3 | 横向变化 |
足球中的无旋转任意球（电梯球）会飘忽不定，这也是空气动力学现象吗？
是的。无旋转（$S \approx 0$）时，球体后方的卡门涡街变得不稳定，横向力会随时间随机波动。这就是导致不规则轨迹变化的“Knuckle Effect”（无旋转飘球效应）的原因。在CFD中重现此现象必须使用LES。
体育用品特有的课题
棒球缝线创造魔球的流体力学
棒球的缝线可不只是装饰哦。缝线导致表面粗糙度不对称，使得只有一侧较早发生湍流转捩。这种左右不对称的转捩会产生横向力，从而形成“Gyroball”（螺旋球）或“Two-seam fastball”（二缝线快速球）的变化。在CFD中精确建模缝线时，只需将缝线角度改变几度，就能再现横向力方向反转的现象。棒球规则中“缝线的数量和配置是固定的”这一规定，实际上在空气动力学上保持着绝妙的平衡，这真是个有趣的话题。
各项的物理意义
- 时间项 $\partial(\rho\phi)/\partial t$：想象一下拧开水龙头的瞬间。最初水流会不稳定地喷溅，过一会儿才变成稳定的水流，对吧？描述这个“变化过程中”的就是时间项。心脏搏动导致血流脉动，发动机阀门每次开闭引起流动变化，这些都是非定常现象。那么定常分析是什么？就是“只观察经过足够长时间、流动稳定后的状态”——也就是将此项设为零。由于计算成本大幅降低，先用定常分析求解是CFD的基本策略。
- 对流项 $\nabla \cdot (\rho \mathbf{u} \phi)$：把落叶扔进河里会怎样？它会随水流被带到下游，对吧？这就是“对流”——流体运动搬运物质的效果。暖风的暖气能到达房间另一端，也是因为空气这个“搬运工”通过对流输送热量。这里有趣的是——此项包含“速度×速度”，因此是非线性的。也就是说，流速越快，此项会急剧增强，变得难以控制。这就是湍流的根本原因。常见的误解：“对流和传导差不多”→ 完全不一样！对流是流动搬运，传导是分子传递。两者效率有天壤之别。
- 扩散项 $\nabla \cdot (\Gamma \nabla \phi)$：有过在咖啡里倒入牛奶后放置不管的经历吗？即使不搅拌，过一会儿也会自然混合。那就是分子扩散。那么下一个问题——蜂蜜和水，哪个更容易流动？当然是水，对吧？因为蜂蜜的粘度（$\mu$）高，所以不易流动。粘度越大，扩散项越强，流体的运动就变得“粘稠”。雷诺数小的流动（缓慢、粘稠）中扩散占主导。相反，雷诺数大的流动中对流占压倒性优势，扩散则成为配角。
- 压力项 $-\nabla p$：推动注射器的活塞，液体就会从针头有力地射出，对吧？为什么呢？因为活塞侧压力高，针头侧压力低——这个压力差就是推动流体的力。水坝放水也是同样原理。天气图上等压线密集的地方会怎样？没错，会刮强风。“有压力差的地方就会产生流动”——这就是纳维-斯托克斯方程压力项的物理意义。这里容易误解的点是：CFD中的“压力”通常指表压而非绝对压力。切换到可压缩分析时结果突然出错，原因可能就是混淆了绝对压力/表压。
- 源项 $S_\phi$：受热的空气会上升——为什么呢？因为比周围空气轻（密度低），被浮力推上去了。这个浮力作为源项添加到方程中。此外，燃气灶火焰产生化学反应热、工厂电磁泵对金属熔液施加洛伦兹力……这些都是“从外部向流体注入能量或力”的作用，都用源项来表示。忘记源项会怎样？自然对流分析中如果忘记加入浮力，流体就完全不动——就像冬天房间里开了暖气，暖空气却不上浮一样，得到物理上不可能的结果。
假设条件与适用范围
- 连续介质假设：克努森数 Kn < 0.01（分子平均自由程 ≪ 特征长度）时成立
- 牛顿流体假设：剪切应力与应变速率呈线性关系（非牛顿流体需要粘度模型）
- 不可压缩假设（Ma < 0.3时）：将密度视为常数。马赫数0.3以上需考虑可压缩性效应
- Boussinesq近似（自然对流）：仅在浮力项中考虑密度变化，其他项使用恒定密度
- 不适用的情形：稀薄气体（Kn > 0.1）、超音速/高超音速流动（需要捕捉激波）、自由表面流动（需要VOF/Level Set等方法）
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 速度 $u$ | m/s | 入口条件中从体积流量换算时，注意截面面积的单位 |
| 压力 $p$ | Pa | 区分表压和绝对压力。可压缩分析中使用绝对压力 |
| 密度 $\rho$ | kg/m³ | 空气：约1.225 kg/m³@20°C，水：约998 kg/m³@20°C |
| 粘性系数 $\mu$ | Pa·s | 注意与运动粘性系数 $\nu = \mu/\rho$ [m²/s] 混淆 |
| 雷诺数 $Re$ | 无量纲 | $Re = \rho u L / \mu$。层流/湍流转捩的判断指标 |
| CFL数 | 无量纲 | $CFL = u \Delta t / \Delta x$。直接关系到时间步长的稳定性 |
数值解法与实现
数值方法
请讲解一下球体空气动力学分析中使用的数值方法。
旋转球体的非定常空气动力学分析计算成本很高。方法的区分使用很重要。
| 方法 | 网格数 | 用途 | 精度 |
|---|---|---|---|
| 定常RANS | 500万--2000万 | 平均$C_D$的概算 | 中 |
| URANS | 1000万--3000万 | 旋转球体的平均空气动力 | 中--高 |
| DDES | 3000万--1亿 | 非定常空气动力、涡结构 | 高 |
| LES | 5000万--3亿 | 无旋转飘球效应、阻力危机 | 最高 |
要解析高尔夫球的凹坑，网格数会变得非常庞大吧。
高尔夫球（直径42.7mm，凹坑300-500个）的LES分析需要1-3亿个网格。不过，也有研究采用等效粗糙度模型来替代凹坑效果的方法。即在壁面函数中设置等效砂粒粗糙度$K_s$。
旋转球体的网格策略
旋转球体的CFD网格有两种主要方法。
1. 滑移网格法
- 在球体周围设置旋转区域，使其物理旋转
- 精度高但计算成本高
- 时间步长: $\Delta t \cdot \omega \cdot d < 1°$（每步旋转小于1度）
2. MRF + 壁面旋转速度
- 定常近似。在球面应用旋转壁面条件
- 只能预测马格努斯力的定常分量
- 无法捕捉非定常涡的动力学
球面的边界层需要$y^+=1$吗？
预测阻力危机必须满足$y^+ < 1$。因为边界层的转捩位置决定了$C_D$，所以粘性底层的分辨率决定了精度。
足球分析示例
足球的面板接缝会影响其空气动力学性能。
2010年世界杯的Jabulani足球被说“飘忽不定”，正是因为这个空气动力学特性吧。
没错。由于面板数量少且接缝浅，导致阻力危机的转捩变得非常急剧，无旋转时的涡脱落不稳定，产生了巨大的横向力波动。
自行车空气动力学
自行车的空气动力学分析近年来也很盛行。骑手的身体占据了总阻力的70-80%。
- 头盔：通风孔布局不同，$C_D$会有5-10%的差异
- 骑手姿势：上半身角度不同，$C_DA$在0.20--0.35m^2之间大幅波动
- 跟骑：进入前车尾流区，阻力可降低30-40%
原来骑手的身体是最大的空气阻力源啊。看来姿势比器材更重要。
公路自行车空气动力学中骑手身体朝向最重要的原因
在公路自行车的CFD分析中，首先令人惊讶的事实是“骑手的身体占据了总阻力约70~80%”。无论把器材（车架、轮组）的空气动力学性能优化得多好，只要骑手稍微抬高头盔，这些优化效果就可能被抵消。职业车队的CFD分析中，会扫描骑手制作等身大的3D模型，以1度为单位改变骑行姿势角度，寻找使CdA（阻力面积）最小化的最优姿势。“姿势比器材更重要”是CFD定量证明的答案。
迎风格式（Upwind）
一阶迎风：数值扩散大但稳定。二阶迎风：精度提高但有振荡风险。高雷诺数流动中必不可少。
中心差分（Central Differencing）
二阶精度，但当Pe数 > 2时会产生数值振荡。适用于低雷诺数的扩散主导流动。
TVD格式（MUSCL、QUICK等）
通过限制器函数抑制数值振荡，同时保持高精度。对捕捉激波或陡峭梯度有效。
有限体积法 vs 有限元法
FVM：自然地满足守恒定律。CFD的主流方法。FEM：对复杂形状、多物理场问题有利。SPH等无网格法也在发展中。
相关主题
なった
詳しく
報告