Stokes流 — CAE术语解说

好问题。一般来说，指的是雷诺数 Re < 1 的领域。特别是在微生物游动、血液中血细胞流动、微小气泡上升等情况下，Re经常小于0.01。例如，直径10μm的细胞以1mm/s的速度在水中运动，其雷诺数约为0.01。

不是感觉，而是通过无量纲化的尺度分析严格证明的。设代表速度为U，长度为L，则惯性项的量级为

完全正确。控制方程组如下所示。

改变是根本的。由于没有非线性项，不需要迭代（如牛顿-拉夫逊法）。离散化后得到**线性方程组**。具体地，速度和压力作为未知数的联立一次方程

主要有两个问题。第一，系数矩阵右下角的零矩阵导致整个矩阵不定。第二，如果速度和压力的离散化空间组合不合适，压力场会出现数值振荡（棋盘式压力）。要防止这种现象，必须使用满足inf-sup条件的单元对，如Taylor-Hood单元（P2-P1）或MINI单元。

取决于问题规模。未知数在10万左右的中小规模问题，直接法求解器（MUMPS、PARDISO）既确定又快速。但数百万以上的大规模问题，直接法会消耗巨大内存，必须采用迭代法。此时可以使用Uzawa法或对块矩阵应用预处理共轭梯度法（PCG）（例如对压力Schur补元使用AMG预处理），这在COMSOL和OpenFOAM的`smoothSolver`中已采用。

首先是**确认雷诺数**。从进口速度、代表长度和流体动粘度计算Re，确保它确实小于1。这是Stokes流分析的前提条件。其次，在物理模型选择的「非压缩流动」中明确选择「Stokes流」或「蠕动流」。对于Ansys Fluent，虽然没有专门的Stokes流模式，但通过将「粘性模型」设为Laminar并设置极低的流速来实现事实上的Stokes流，但拥有专门模式的软件会更方便。

完全正确。Stokes流不会形成速度边界层，因为方程不是抛物型而是椭圆型。因此，壁面不需要极细的边界层网格，这只是浪费计算资源。相反，应该确保整个域内具有品质良好的均匀网格来精确表达几何形状（特别是角和狭缝）。单元次数应选择满足inf-sup条件的组合（例如二次速度/一次压力）。

用标准基准问题验证最可靠。例如，**圆柱周围的Stokes流**很适合。当直径D的无限长圆柱被垂直一致流U冲击时，单位长度的阻力F可以与理论公式

Fluent没有「Stokes流」专用物理模型。因此，选择Laminar模型，设置极低的流速（如0.001mm/s）来实质上将惯性项影响压低至接近零的「变通办法」。求解器设置选择「Pressure-Based」，压力-速度耦合求解法推荐「Coupled」。但需要注意，这种方法不是严格的Stokes求解器，取决于网格和求解器设置，可能还会残留微小惯性影响。

**COMSOL Multiphysics**明显领先。在「CFD模块」或「微粒流模块」中，有明确的「Stokes流」界面。使用它会自动将之前提到的线性方程组用优化的求解器（直接法MUMPS或PARDISO）高效求解。另外，**Abaqus/CFD**（目前正逐步废止）和开源的**OpenFOAM**的`icoFoam`求解器（非压缩层流）也可以用极低的雷诺数来实现。OpenFOAM还内置`smoothSolver`这种适合Stokes类问题的预处理求解器。

这里体现了不同软件的核心差异。Ansys Fluent和STAR-CCM+主要擅长结构网格（六面体），但复杂形状也能用非结构网格。COMSOL和Abaqus基于有限元法，从一开始就以四面体等非结构网格为标准。特别是COMSOL可以关闭「边界层网格」的添加，对Stokes流来说就很合理。无论哪种，都要注意单元质量（纵横比、偏斜角），通过COMSOL的「网格品质检查」功能或Ansys Meshing的「Element Quality」指标来监控生成过程。

这是经典的**压力振荡问题**。原因是速度和压力的离散化组合不满足inf-sup条件。具体来说，一次单元速度加一次单元压力（P1-P1）等组合就会出现。对策是：1) 把单元类型改为安定化的组合如Taylor-Hood单元（P2-P1），或2) 引入压力稳定化方案（PSPG）。COMSOL在选择物理界面时会自动选择稳定的单元，但如果用户修改了设置就可能遇到这个问题。