斯托克斯流（低雷诺数）
斯托克斯流（低雷诺数）的理论基础
概述
老师，斯托克斯流和普通的NS方程有什么区别？
当雷诺数非常小（$Re \ll 1$）时，惯性力相对于粘性力可以忽略。从NS方程中去掉对流项（非线性项）后得到的线性方程称为斯托克斯方程。这在微小粒子沉降、微流体器件、生物流体力学等领域中很重要。
斯托克斯方程
稳态、不可压缩斯托克斯方程如下所示。
动量方程:
连续方程:
这是从NS方程中去掉 $\rho(\partial\mathbf{u}/\partial t + \mathbf{u}\cdot\nabla\mathbf{u})$ 后的形式。压力梯度和粘性力完全平衡。
在什么条件下才能忽略对流项？
用雷诺数判断。
| 对象 | 典型雷诺数 | 斯托克斯近似 |
|---|---|---|
| 细菌游动 | $10^{-5}$～$10^{-4}$ | 非常好 |
| 精子游动 | $10^{-2}$ | 好 |
| 粒子沉降（$d=10\mu$m, 水中） | $10^{-3}$ | 好 |
| 高粘度聚合物流 | $10^{-2}$～$10^{-1}$ | 大致有效 |
| 微流路（$L=100\mu$m） | $0.1$～$10$ | 视情况而定 |
斯托克斯阻力定律
当球体（半径 $R$）在粘性流体中以一致速度 $U$ 运动时，其阻力由斯托克斯在1851年推导。
这就是斯托克斯阻力定律。用阻力系数表示为，
其中 $Re = \rho U (2R)/\mu$ 定义。当 $Re < 1$ 时，该关系成立。
斯托克斯流的特性
斯托克斯流有独特的数学性质。
- 可逆性：时间反演后流动相同（倒放视频无法区分）
- 瞬时响应：因为没有惯性，对边界条件的变化立即响应
- 线性性：解的叠加原理成立。可以解析处理多个粒子的相互作用
- 唯一性：给定边界条件后，解唯一
可逆性很有趣。普通流动不可能这样。
在Taylor-Couette装置中，用粘性流体中的墨水条进行"搅拌"然后反向旋转就会恢复，这是一个经典的斯托克斯流可逆性演示。
细菌"反向游泳也能前进"的世界
细菌（如大肠杆菌）鞭毛旋转游动的机制只能在斯托克斯流（Re≪1）的世界中成立。有一个著名的结果叫"扇贝定理（Scallop theorem）"，它表明"可逆的运动（仅往复开关）在Re≪1的流中不产生推进力"。这就是为什么细菌不能用像人类那样的游泳方式的原因。细菌通过螺旋状旋转鞭毛实现非可逆运动，从而绕过这个矛盾。微机器人和医疗纳米机器的设计中，这个原理是重要的约束条件。
斯托克斯流（低雷诺数）的数值计算方法
斯托克斯方程的数值解法
请教一下斯托克斯方程的数值求解方法。
斯托克斯方程是线性的，所以比NS方程更容易求解。不过速度-压力的耦合（鞍点问题）仍然存在。
鞍点问题的定式化
用FEM对斯托克斯方程离散化后，得到如下鞍点型线性方程组。
其中 $A$ 是粘性项的刚度矩阵，$B$ 是散度算子的离散版本。
有零对角块看起来很难求解。
没错。需要速度-压力单元组合满足inf-sup条件（LBB条件）。
| 速度单元 | 压力单元 | LBB稳定性 | 名称 |
|---|---|---|---|
| P2（二次三角形） | P1（一次三角形） | 稳定 | Taylor-Hood单元 |
| P1+bubble | P1 | 稳定 | MINI单元 |
| Q2（二次四边形） | Q1（一次四边形） | 稳定 | 标准 |
| P1 | P1 | 不稳定 | 需要稳定化 |
| P1 | P0 | 不稳定 | 不可用 |
FVM求解
有限体积法中SIMPLE算法可以直接使用，但由于Re数极小，需要特殊考虑。
- 压力松弛系数：可以设置得比通常情况大（0.5～0.8）。因为没有对流项，所以更稳定
- 对流格式：不需要（对流项为零）
- 收敛：线性问题，通常数十次迭代就收敛
边界积分法（BEM）
利用斯托克斯流的线性性，有一个强大的方法是边界积分法。使用格林函数（斯托克斯点源），只在边界面上配置未知量，不需要内部。
这里 $G_{ij}$ 是Oseen-Burgers张量（自由空间的格林函数）。3D问题中，不需要体积网格，计算量大幅减少。
斯托克斯流中只用边界就能求解？NS方程中这不可能吧？
对。NS方程的非线性性导致不存在格林函数。斯托克斯流的线性性使这种方法成为可能。这是一个特权。
斯托克斯流的数值求解——边界元法(BEM)为什么比有限体积法更优
斯托克斯方程（低Re极限）是线性偏微分方程，因此边界元法（BEM）特别有效。BEM不需要对流体域内部进行网格划分，只对边界面进行离散化，所以3D问题中只需要2D网格。特别是对于外部流（开放区域），BEM自然处理无穷远的边界条件，而有限体积法必须人为设置"远场边界"。但BEM矩阵（N²的稠密矩阵）对网格数N的存储O(N²)和求解O(N³)开销大，解决方案是配合快速多极法（FMM: Fast Multipole Method）将其压缩到O(N log N)。这个方法对流体中的多粒子群、气泡群等多体问题特别有效。
斯托克斯流（低雷诺数）的实务应用
实践指南
请教一下在实务中进行斯托克斯流分析的情况。
主要应用领域如下。
粒子沉降速度的计算
球形粒子的终端沉降速度可以直接从斯托克斯阻力定律得出。当重力与阻力平衡时，
| 粒子 | 直径 | 流体 | 沉降速度 | Re |
|---|---|---|---|---|
| 花粉 | 30 $\mu$m | 空气 | 2.7 cm/s | 0.05 |
| 沙粒 | 100 $\mu$m | 水 | 0.8 cm/s | 0.08 |
| 血球 | 8 $\mu$m | 血浆 | 0.4 $\mu$m/s | $3\times10^{-6}$ |
| 雾滴 | 10 $\mu$m | 空气 | 0.3 cm/s | 0.02 |
当 $Re > 1$ 时需要修正吗？
没错。需要使用Oseen修正或Schiller-Naumann公式。
在CFD的粒子追踪（DPM模型）中，默认使用这个带修正的拖力模型。
微流体器件的设计
在微流路（宽度10～100 $\mu$m）中，$Re$ 为 $O(1)$～$O(10^{-2})$，斯托克斯近似有效。
| 设计参数 | 计算方法 |
|---|---|
| 压力损失 | $\Delta p = \frac{12\mu L Q}{wh^3}$（矩形截面近似） |
| 混合长 | $L_{\text{mix}} \sim Pe \cdot w$（Pe数 $= Uw/D$） |
| Dean流二次流 | 用 $De = Re\sqrt{d/R}$ 评估强度 |
我听说微流路中混合很困难？
这是因为斯托克斯流的可逆性。没有湍流混合，只能依靠扩散。所以使用锯齿形流路或鱼骨状结构来引入混沌对流的设计很常见。
CFD设置技巧
整理一下低Re数分析特有的注意点。
| 项目 | 设置 | 理由 |
|---|---|---|
| 乱流模型 | Laminar（层流） | Re < 1 不会出现乱流 |
| 对流格式 | 中心差分没问题 | Pe数很小 |
| 密度处理 | 不可压缩 | Ma数极小 |
| 浮力 | Boussinesq或密度差 | 沉降问题必须考虑重力 |
| 收敛 | 容易（数十次迭代） | 接近线性问题 |
低Re数分析收敛很容易呢。
是的。但物理现象变成尺度相关（表面张力、电渗流、范德瓦尔斯力等），需要正确识别支配力。
MEMS流道中的斯托克斯流——微阀设计中的CFD应用实例
MEMS器件的微细流道（宽1～100μm）中，雷诺数Re≪1的斯托克斯流（蠕动流）占主导，惯性力可以忽略。医疗药液输送泵的微阀设计中，阀体和阀座间隙（亚μm级）的斯托克斯流分析很关键。实际设计案例中，对硅微阀（狭缝宽0.5μm）的CFD分析（使用斯托克斯求解器），在阀前后压差0.1kPa条件下预测流量系数Cd，并评估制造公差±0.1μm对流量的灵敏度，这样就能把初期原型试制次数从3次减少到1次。相关文献发表在学术期刊Microfluidics and Nanofluidics上。
斯托克斯流（低雷诺数）的软件比较
低Re数流的商用工具支持
用什么求解器最适合斯托克斯流？
所有通用CFD求解器都能在层流模式下求解斯托克斯流。但对于特殊应用，FEM或BEM可能更有优势。
求解器比较
| 求解器 | 方法 | 低Re流的特点 | 应用领域 |
|---|---|---|---|
| Ansys Fluent | FVM | Laminar设置兼容。DPM粒子追踪 | 粒子沉降、过滤器 |
| STAR-CCM+ | FVM | Laminar + DEM耦合强大 | 粉体、浆液 |
| COMSOL | FEM | 有专用"Creeping Flow"物理模型 | 微流体、MEMS |
| OpenFOAM | FVM | icoFoam（层流专用） | 学术、微流体 |
| Stokesflow (BEM) | BEM | 边界元法。多粒子问题最优 | 生物流体、胶体 |
COMSOL有专门的"Creeping Flow"物理模型？
是的。COMSOL的"Creeping Flow"物理接口预设了斯托克斯方程（没有对流项），这样稳定性更好，即使在极低Re数下也能稳定计算。
粒子追踪模型比较
| 方法 | Fluent | STAR-CCM+ | COMSOL |
|---|---|---|---|
| DPM (Discrete Phase Model) | 支持 | 支持 | Particle Tracing |
| DEM (Discrete Element Method) | 有DEM耦合 | DEM标准配置 | 有限制 |
| 粒子间相互作用 | 简化模型 | 完全DEM对应 | 简化 |
| 布朗运动 | Brownian force选项 | 支持 | 支持 |
| 电泳力 | 通过UDF实现 | Field force | Electrophoresis |
适合微流体器件设计的工具
微流体设计中COMSOL特别强大。多物理耦合集成很容易，可以组合以下模块。
- Creeping Flow + Transport of Diluted Species（扩散混合）
- Creeping Flow + AC/DC Module（电渗流）
- Two-Phase Flow + Surface Tension（液滴操作）
多物理耦合时COMSOL的优势就显示出来了。
反过来，如果只是简单粒子沉降或管内压损计算，用Fluent或OpenFOAM的Laminar模式就足够了。工具选择取决于要解决的物理复杂度。
斯托克斯流对应CFD工具——COMSOL Multiphysics的蠕动流分析功能
在低雷诺数（斯托克斯流）分析领域，COMSOL Multiphysics得到许多研究者的青睐。COMSOL拥有"Creeping Flow（蠕动流）接口"，可以直接求解去掉惯性项的斯托克斯方程，用图形界面很容易设置。另外与电渗流（Electroosmotic Flow）等多物理耦合的设置很直观，在微流体研究中采用率很高。相比之下OpenFOAM没有StokesFoam，需要极端压缩simpleFoam的松弛系数来得到低Re解，这样效率很低。通用CFD（Fluent）虽然能得到低Re解，但收敛所需的迭代次数通常是斯托克斯专用求解器的10～100倍。
斯托克斯流（低雷诺数）的前沿研究
前沿主题
现在还在研究斯托克斯流吗？
生物流体力学、微机器人学、胶体科学等，低Re数世界是非常活跃的研究领域。
微型游动者的流体力学
从斯托克斯流的可逆性，可以导出"扇贝定理"（Scallop theorem, Purcell 1977）。时间可逆运动（仅1自由度开闭）不能产生推进。
微生物通过以下策略实现推进。
- 旋转鞭毛（大肠杆菌）：螺旋波传播，非互易运动
- 纤毛打动（草履虫）：有效打与回收打的不对称性
- 波动运动（精子）：正弦波传播
所以像扇贝那样开闭不能游动。
这是斯托克斯流可逆性的直接结果。微机器人设计中，实现非互易运动的方法是核心课题。
斯托克斯方程的高级解析解
球以外的形状对斯托克斯阻力的解析解也很重要。
| 形状 | 阻力 | 备注 |
|---|---|---|
| 球 | $F = 6\pi\mu RU$ | Stokes (1851) |
| 椭球体（长轴方向） | $F = \frac{16\pi\mu aU}{2\ln(2a/b) - 1}$ | Oberbeck (1876)，细长体近似 |
| 圆盘（法向） | $F = 16\mu RU$ | Lamb |
| 圆柱（无限长，2D） | 无解（斯托克斯矛盾） | 需要Oseen修正 |
2D圆柱没有解？
这是著名的"斯托克斯矛盾"。2维中粘性流的斯托克斯方程不能找到满足边界条件的解。需要用Oseen方程（对流项线性化版）来解决。3维中球形无此问题，解存在。
非牛顿流体的斯托克斯流
高粘度聚合物或血液不是牛顿流体。它们表现出剪切变稀性（shear-thinning）。
Cross模型和Carreau模型在CFD中常用。在CFD求解器中作为材料模型设置。
活性物质的群体运动
许多微型游动者的群体运动是"活性物质"，是物理学的最前沿课题。每个游动者产生的斯托克斯点源类的流场相互作用，导致有序图案或湍流行为的出现。
低Re数的世界看似简单，实际上物理很丰富。
从斯托克斯时代到现在已经170多年，新发现还在不断涌现。Re数小不代表物理简单。
斯托克斯流的前沿——微型游动者推进机制与CFD
细菌和精子等微生物（微型游动者）在 Re≪1 的斯托克斯流域中自我推进。此流域中成立"扇贝定理（Scallop Theorem）"，对称往复形状变化不能前进——这就是细菌旋转螺旋鞭毛前进的原因所在。最前沿研究利用斯托克斯流这一物理约束，设计纳米机器人（药物运输微型机器）的推进机制，CFD被广泛使用。边界元法（BEM: Boundary Element Method）因为流体-结构耦合计算成本小，是斯托克斯流分析的最优方法。MIT和哈佛大学的微型游动者CFD研究致力于医疗应用的实用化。
斯托克斯流（低雷诺数）的故障排除
故障排除
低Re数分析常见的故障有吗？
与高Re数不同的问题会出现。精度和物理模型问题是主要的。
1. 压力场不稳定（FEM的情况）
现象：压力出现棋盘状振荡，或压力为零附近的异常值。
原因：速度-压力单元的组合不满足inf-sup条件（LBB条件）。
对策：
- 使用Taylor-Hood单元（P2/P1）
- 采用稳定化技术（PSPG, GLS）
- COMSOL的Creeping Flow物理会自动选择稳定的定式化
2. 粒子的阻力与斯托克斯定律不符
CFD计算的球的拖力与 $6\pi\mu RU$ 偏差超过10%…
| 原因 | 对策 |
|---|---|
| 计算域太小 | 外部域至少球径的50倍 |
| 网格太粗 | 球表面至少30个单元（周向） |
| 壁面影响 | 远场边界处一致流条件 |
| Re > 0.1 | Oseen修正：$C_D = (24/Re)(1 + 3Re/16)$ |
斯托克斯解假设无限域，有限计算域中壁面效应不可避免。外部边界要设置足够远。
3. 微流路的流量与理论值不符
现象：与Hagen-Poiseuille理论流量差异大。
检查事项：
- 非圆形流道需要考虑截面形状的修正系数
- 入口助走区长度（$L_e \approx 0.06 Re D$，斯托克斯流中几乎为零）
- 壁面滑移条件：微观尺度下克努德森数 $Kn = \lambda/L$ 变得有限，no-slip可能失效
4. 非牛顿流体收敛困难
现象：幂律流体粘度局部变得非常大（或很小），无法收敛。
对策：
- 设置粘度上下限（$\mu_{\min} < \mu_{\text{eff}} < \mu_{\max}$）
- 使用Carreau模型（用零剪切粘度和无限剪切粘度定义上下限）
- 降低松弛系数
- 以牛顿流体（平均粘度）的解作为初值
5. 多粒子问题中相互作用被低估
粒子很多时，逐个计算不行吗？
孤立粒子假设仅在粒子间距足够大时有效。体积分率超过0.1时，粒子间的流体动力学相互作用不可忽视。
对策：
斯托克斯流也有这么多故障啊。
方程线性不代表分析简单。标度问题、壁面效应、非牛顿性等，低Re特有的课题要理解透彻。
斯托克斯流计算"不对"——怀疑雷诺数
斯托克斯流（Re≪1）分析最常见的失败是"低速就用斯托克斯近似"的思维误区。雷诺数不仅与流速有关，还取决于代表长度。1μm粒子在水中以相同流速运动可能Re≈0.001（完全斯托克斯域），但同样流速下1mm气泡可能Re≈1，接近慣性项生效范围。分析与实验不符时，首先检查问题的代表长度——这是斯托克斯流故障排除的铁律。
相关主题
更详细
错误