Fluent压力速度耦合方案的选取

Fluent压力基求解器里“Pressure-Velocity Coupling”那一栏的下拉菜单，在 Solution Methods 面板里，稳态和瞬态都在同一位置。路径 /solve/set/p-v-coupling 对应五种方案：20-SIMPLE、21-SIMPLEC、22-PISO、24-Coupled、25-Fractional Step，TUI 直接敲编号就行。D:\Fluent_work\solver_setup\coupling_schemes.txt 里存着五种方案的默认松弛因子对照表。

SIMPLE 是默认选项，几乎所有不可压缩流动都能跑，内存需求最低，但收敛慢，每步迭代压力修正方程之后速度场不能完全满足动量方程，得反复迭代。松弛因子默认压力0.3动量0.7，收敛困难的时候往下压，压力到0.2动量到0.5，压得太低收敛更慢，压得太高直接发散。稳态计算起步用SIMPLE最稳，算到残差降了一个数量级之后再看要不要换。

SIMPLEC 是 SIMPLE-Consistent，做了速度增量近似，压力修正方程里少了一项，收敛速度比 SIMPLE 快，亚松弛因子可以给得更大，压力修正松弛因子一般设1.0。简单层流问题用 SIMPLEC 通常比 SIMPLE 快不少。但网格歪曲度大的时候，压力修正松弛因子给到1.0会导致不稳定，得退到0.7，或者干脆换回 SIMPLE。Skewness Correction 那个勾选框在高歪曲度网格上打开，但打开之后求解方案的鲁棒性会变差，得权衡。

PISO 多了邻域校正和网格歪曲校正，压力修正方程解完之后还要再做一步动量修正，压力场更准，收敛加速。瞬态计算强烈推荐用 PISO，尤其是时间步长给得大的时候。亚松弛因子设1.0或者接近1.0，PISO 在大时间步下也能稳住。大涡模拟不建议用 PISO，LES 本身时间步很小，PISO 每一步多出来的计算量反而拖慢整体效率，用 SIMPLE 或者 SIMPLEC 更划算。网格质量差的稳态问题也可以用 PISO，只开 Skewness Correction 的时候压力和动量松弛因子加起来等于1，压力0.3动量0.7；只开 Neighbor Correction 的时候所有方程松弛因子设1.0；两个都开就按 Neighbor Correction 的规则来。

Coupled 把动量方程和压力基连续方程耦合起来一起解，收敛速度比分离式算法快很多，文献里说最多快5倍，但内存需求是分离式算法的1.5到2倍，一百万网格大概吃6到7个G。稳态单向流用 Coupled 更高效，鲁棒性也强。瞬态计算时间步大或者网格质量差的时候，Coupled 基本是必须的。VOF 多相流、多模型耦合的工况也优先用 Coupled。有个坑是 Coupled 跟 NITA 求解器不兼容，选了 Coupled 就不能用 Non-Iterative Time Advancement。多孔介质问题有时候 Coupled 不收敛，退回去用分离式算法就行。用 Coupled 的时候控制求解过程有两个手段，Courant 数默认200，多相流和燃烧问题降到10到50，Courant 数小计算稳，大算得快但可能发散。Pseudo Transient 是另一个控制手段，比普通 Coupled 更稳，网格质量差的时候打开，内流用 Automatic 加 Conservative，外流用 Automatic 加 User-Specified 并把特征长度填进去，CHT 问题里固体时间尺度用 Automatic，2020R2 之前默认的拟时间步太小导致固体热计算极慢，得手动调。

Fractional Step 只在 NITA 模式下可用，对不可压缩流保持了形式上的二阶时间精度，不需要每个时间步做昂贵的外迭代，LES 这种计算量大的瞬态模拟用它效率高。但它跟 VOF 多相流不兼容，做 VOF 的时候别选这个，老老实实用 PISO 或者 Coupled。

实际选的时候按工况分。稳态不可压缩，SIMPLE 或者 SIMPLEC 起步，SIMPLEC 快一点，网格歪曲度大或者收敛困难就退回 SIMPLE。稳态高速可压缩，Coupled 配 Pseudo Transient。瞬态不可压缩，PISO 为主，时间步小的话也可以试试 SIMPLE。瞬态 VOF，PISO 或者 Coupled，Courant 数锁死在0.25以下。瞬态 LES，SIMPLE 或者 SIMPLEC，别用 PISO 浪费计算资源。多孔介质，分离式算法优先。D:\Fluent_work\solver_setup\下面那个弯管稳态 case 用的 SIMPLEC，残差降到1e-5用了八百步，换 Coupled 之后四百步就到了，内存占用从3.2G涨到6.1G，跑大模型之前先看机器内存够不够。

不收敛的时候先别急着换方案，把松弛因子往下降一档再跑五十步看看残差趋势。残差震荡不降，降松弛因子。残差发散，换方案。残差先降后升，初始化问题。残差降不动了但监测的物理量还在漂，收敛判据太松，残差标准从1e-3收紧到1e-5。SIMPLE 收敛困难换 SIMPLEC 试试，SIMPLEC 不稳定换 Coupled 加 Pseudo Transient，Coupled 内存不够换回 SIMPLE 加多重网格。D:\Fluent_work\solver_setup\下面那个多孔介质 case 就是 Coupled 不收敛换回 SIMPLE 跑通的，残差从震荡变成单调下降，代价是迭代步数多了三倍，那个内存占用是……

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删