Fluent多相流模型选型与设定避坑指南：VOF、Mixture、Eulerian实战

简介：这份PDF资料面向CFD仿真初学者与工程技术人员，系统梳理Fluent中多相流模型的选择逻辑与设定方法，帮助解决气-液、气-固、液-固等流动场景下模型选型困难的问题。内容从多相流分类入手，覆盖气泡流、液滴流、活塞流、分层自由面流、流化床、泥浆流、水力运输与沉降等典型模式，并逐一讲解VOF模型、混合模型与欧拉模型的原理、适用条件及局限性，同时对比欧拉-欧拉法与欧拉-拉格朗日法的差异，给出体积分率、相间耦合、计算精度与耗时等选择原则。资源包为1个PDF文件，大小约2.22MB，便于随时查阅与对照学习。目前已有653人学习，适合需要快速掌握多相流建模思路、提升仿真选型准确性的读者参考。

1. 多相流仿真翻车，八成栽在模型选错这一步

做流体仿真的人大多有过这种经历：几何画好了，网格也过了无关性验证，边界条件反复核对三遍，一跑起来残差曲线却像心电图，相体积分数场要么糊成一团，要么直接发散。回头查半天，问题往往不在网格也不在求解器，而是最开始那个下拉框——多相流模型选错了。Fluent 里多相流模型选择与设定这件事，看着只是点几下鼠标，实际上决定了你后面所有参数有没有意义。欧拉-欧拉和欧拉-拉格朗日两套框架，VOF、Mixture、Eulerian 三个模型，各自适用的物理场景完全不同，选错了再调松弛因子也是白费功夫。这篇笔记就按我自己的实操顺序，把模型怎么选、参数怎么设、哪里容易踩坑讲清楚，适合刚接手气液两相或颗粒流仿真的工程师，也适合做过单相流、第一次碰多相流的人。

2. 三个模型到底怎么选：从物理场景倒推，而不是从菜单顺序猜

选模型这件事，最忌讳的是打开 Fluent 看到下拉框里 VOF 排第一就顺手选它。正确的顺序是先想清楚你要回答什么问题，再倒推模型。我一般会先问三个问题：分散相体积分数大概多少？相界面是不是清晰可追踪？颗粒/气泡有没有明显的粒径分布和受力？这三个问题的答案基本就把模型锁死了。

2.1 VOF、Mixture、Eulerian 的适用边界

VOF 的核心是追踪自由液面，它求解的是各相体积分数的输运方程，界面通过几何重构保持锐利。适合的场景是分层流、溃坝、液滴撞击、气泡在静止液体中上升这类界面清晰、分散相体积分数很低的情况。它的代价是对界面附近的网格分辨率要求极高，界面越锐利，网格越要细，时间步也受 CFL 条件限制。

Mixture 模型是简化版的欧拉模型，它求解混合相的动量方程，用滑移速度描述相间相对运动。适合气泡/颗粒体积分数中等（大概 10% 到 20%）、相间耦合不算特别强、你又不想付出完整欧拉模型计算代价的场景。它的短板是滑移速度模型是经验性的，相间作用力描述不如 Eulerian 完整。

Eulerian 模型对每一相都求解独立的动量方程和连续性方程，相间通过曳力、升力、虚拟质量力等封闭模型耦合。适合鼓泡床、流化床、气力输送这类分散相体积分数高、相间作用强烈的场景。代价是方程数量随相数线性增长，收敛难度和计算量都明显上升。

这张表不是硬性分界线，体积分数只是参考。真正决定的是相间耦合强度和你要不要追踪界面。比如体积分数只有 5% 但你要看单个气泡变形，那还是 VOF；体积分数 15% 但相间作用很弱，Mixture 就够。

2.2 从物理问题到模型选择的三步判断

第一步，判断分散相是否连续。如果分散相在空间上是连通的（比如管道里的气塞、自由液面），选 VOF。如果分散相是离散的液滴、气泡、颗粒，进入第二步。

第二步，估算分散相体积分数。在你要关注的区域里，分散相占的体积比例大概是多少。低于 10% 可以考虑 VOF 或 DPM，10% 到 20% 看 Mixture，超过 20% 基本就是 Eulerian。

第三步，判断相间作用力是否主导流动。如果曳力、升力、湍流扩散对分散相运动影响很大，Eulerian 的封闭模型更完整。如果只是想知道宏观混合效果，Mixture 的滑移速度够用。

提示：模型选择没有回头路。一旦选了 Eulerian 跑了几百步再想换 VOF，基本等于重做。所以前期花十分钟想清楚，比后面花十小时调参划算。

2.3 一个最小可复现的 VOF 算例配置

光说选型不够，得能跑起来。下面是一个二维溃坝的最小配置，用 VOF 追踪水气界面。这个算例我用来验证模型设置是否正确，跑通了再换自己的几何。

# Fluent Python 脚本：二维溃坝 VOF 最小配置 # 假设网格已导入，单位为米，重力方向 -y import ansys.fluent.core as pyfluent solver = pyfluent.launch_fluent(precision="double", processor_count=4) solver.tui.file.read_case("dam_break_2d.cas.h5") # 1. 求解器设置：瞬态、压力基、重力开启 solver.tui.define.models.unsteady_2d("yes") solver.tui.define.operating_conditions.gravity("yes", "0", "-9.81", "0") # 2. 多相流模型：VOF，两相，显式格式 solver.tui.define.models.multiphase.model("vof") solver.tui.define.models.multiphase.number_of_phases("2") solver.tui.define.models.multiphase.vof_parameters( "explicit", "yes", "0.25", "yes", "1e-5", "1" ) # 参数依次为：显式格式、开启界面锐化、界面压缩因子0.25、 # 开启体积分数截断、截断下限1e-5、截断上限1 # 3. 相设置：主相空气，次相水 solver.tui.define.phases.set_domain_properties.phases("air", "water") solver.tui.define.phases.set_domain_properties.phase_domain( "water", "material", "water-liquid" ) # 4. 表面张力：水气界面 0.072 N/m solver.tui.define.phases.set_domain_properties.surface_tension( "yes", "0.072", "constant" ) # 5. 初始化：patch 左侧为水 solver.tui.solve.initialize.compute_defaults.all_zones() solver.tui.solve.initialize.initialize_flow.all_zones() solver.tui.solve.patch( "yes", "water", "volume_fraction", "1", "0", "0.5", "0", "1", "0" ) # patch 区域：x 从 0 到 0.5，y 从 0 到 1，水体积分数设为 1 # 6. 时间步设置：固定步长 0.001 s，总步数 2000 solver.tui.solve.set.transient_controls.time_step_size("0.001") solver.tui.solve.set.transient_controls.number_of_time_steps("2000") solver.tui.solve.set.transient_controls.max_iter_per_time_step("40") # 7. 求解 solver.tui.solve.dual_time_iterate("2000", "40")

这段脚本的关键参数有三个。界面压缩因子 0.25 是控制界面锐利程度的，值越大界面越锐利但越容易发散，一般从 0.25 开始试。体积分数截断下限 1e-5 是为了避免出现负体积分数导致计算崩溃，这个值不能设太大，否则质量不守恒。时间步 0.001 秒对应 CFL 数大概在 1 附近，如果网格更细，时间步要相应减小。

跑通之后看水相体积分数云图，界面应该是一条相对清晰的斜线，而不是糊成一片的渐变。如果界面糊了，先检查网格分辨率，再检查界面压缩因子。

3. 相间作用力设定：Eulerian 模型里最容易被忽略的封闭项

选了 Eulerian 模型之后，很多人以为把相数设好、材料选好就能跑，结果发现气泡全部沉底或者颗粒全部吹飞。问题出在相间作用力上。Eulerian 模型对每一相求解动量方程，相与相之间怎么交换动量，靠的就是这些封闭模型。它们不是可选项，是必选项。

3.1 曳力模型的选择依据

曳力是相间动量交换的主要来源。Fluent 内置了十几种曳力模型，常用的有 Schiller-Naumann、Morsi-Alexander、Grace、Tomiyama 等。选哪个取决于分散相的雷诺数和形状。

Schiller-Naumann 适合刚性球形颗粒，雷诺数在 1000 以内。Morsi-Alexander 适合球形颗粒但雷诺数范围更宽。Grace 适合气泡在液体中上升，考虑了气泡变形。Tomiyama 适合气泡在污染液体中的上升，考虑了表面活性剂的影响。

我一般先用 Schiller-Naumann 跑一个基线，看流场是否合理，再根据分散相的实际形态换更合适的模型。如果分散相是气泡且变形明显，Grace 或 Tomiyama 更合适。如果分散相是固体颗粒，Schiller-Naumann 或 Morsi-Alexander 就够。

# 设置气液两相曳力模型为 Grace solver.tui.define.phases.set_domain_properties.interaction_domain( "air", "water", "drag", "grace" ) # 设置气泡直径 2 mm，这是 Grace 模型需要的参数 solver.tui.define.phases.set_domain_properties.phase_domain( "air", "diameter", "0.002" )

气泡直径这个参数很关键。Grace 模型用它计算气泡的终端上升速度，设错了曳力会差一个量级。如果气泡粒径有分布，可以用群体平衡模型（PBM）耦合，但那是另一个话题了。

3.2 升力、壁面润滑与虚拟质量力什么时候要开

升力在剪切流中不可忽略。气泡或颗粒在速度梯度场中会受到垂直于相对速度方向的力，这个力会把分散相推向管壁或管中心，直接影响相分布。如果发现分散相在壁面附近异常聚集，先检查升力模型。

壁面润滑力在分散相靠近壁面时起作用，防止分散相穿透壁面。如果近壁面网格比较粗，这个力可以不开；如果近壁面网格很细且分散相有贴壁趋势，建议开启。

虚拟质量力在分散相密度远小于连续相时重要，比如气泡在水中。它的物理意义是气泡加速时周围液体也跟着加速，相当于气泡质量增加了。如果分散相是固体颗粒且密度比连续相大，这个力可以忽略。

# 开启升力、壁面润滑力、虚拟质量力 solver.tui.define.phases.set_domain_properties.interaction_domain( "air", "water", "lift", "yes" ) solver.tui.define.phases.set_domain_properties.interaction_domain( "air", "water", "wall_lubrication", "yes" ) solver.tui.define.phases.set_domain_properties.interaction_domain( "air", "water", "virtual_mass", "yes" )

注意：这些力不是开得越多越好。每开一个力就多一组经验系数要调，收敛难度也上升。我一般先只开曳力跑通，再逐个加，每加一个看残差和相分布有没有异常。

3.3 相间传质与传热的耦合设置

如果两相之间有蒸发、冷凝或化学反应，需要开启相间传质。Fluent 提供了蒸发-冷凝模型、沸腾模型、组分输运模型等。设置传质时最容易翻车的地方是传质系数和界面面积的估算。

传质速率等于传质系数乘以界面面积乘以浓度差。界面面积在 Eulerian 模型里是通过分散相体积分数和气泡/颗粒直径算出来的。如果气泡直径设得不准，界面面积就错了，传质速率跟着错。

# 开启蒸发-冷凝模型，水在 373 K 蒸发 solver.tui.define.phases.set_domain_properties.phase_domain( "water", "phase_change", "evaporation_condensation", "yes" ) solver.tui.define.phases.set_domain_properties.phase_domain( "water", "saturation_temperature", "373.15" ) solver.tui.define.phases.set_domain_properties.phase_domain( "water", "latent_heat", "2257000" )

饱和温度和潜热这两个参数必须和材料属性一致。如果材料库里的水是默认的，饱和温度是 373.15 K，潜热是 2257000 J/kg。如果用了自定义材料，这两个值要手动核对。

4. 求解器与松弛因子：多相流收敛难，多半是这里没配对

多相流的收敛难度比单相流高一个量级，因为方程数量多了，相间耦合又强。很多人跑单相流习惯了默认松弛因子，换到多相流直接发散，然后开始怀疑网格。其实大部分时候是求解器设置和松弛因子的问题。

4.1 压力-速度耦合方式的选择

Fluent 提供了 SIMPLE、SIMPLEC、PISO、Coupled 几种压力-速度耦合方式。多相流里我一般用 SIMPLE 或 Coupled。SIMPLE 稳定但收敛慢，Coupled 收敛快但对内存要求高。如果计算资源够，Coupled 是更好的选择。

瞬态问题用 PISO 也可以，它对时间步的适应性更好，但每个时间步的计算量更大。如果时间步很小且流场变化剧烈，PISO 值得一试。

# 设置压力-速度耦合为 Coupled solver.tui.solve.set.p_v_coupling("coupled") # 如果收敛困难，改回 SIMPLE 并降低松弛因子 # solver.tui.solve.set.p_v_coupling("simple")

4.2 体积分数方程的求解策略

体积分数方程是多相流特有的，它的求解策略直接影响界面锐利度和收敛性。VOF 模型里，体积分数方程可以用显式或隐式格式。显式格式对时间步要求严格，但界面更锐利；隐式格式时间步可以大一些，但界面容易扩散。

Eulerian 模型里，体积分数方程的松弛因子一般设在 0.2 到 0.5 之间。如果相分布变化剧烈，先降到 0.2 跑几百步，稳定后再逐步提高。

# 设置体积分数方程的松弛因子 solver.tui.solve.set.relaxation("volume_fraction", "0.3") # 如果发散，降到 0.1 # solver.tui.solve.set.relaxation("volume_fraction", "0.1")

4.3 时间步长与库朗数的实操关系

瞬态多相流的时间步长不是随便设的。VOF 模型里，时间步长要满足 CFL 条件：库朗数等于速度乘以时间步长除以网格尺寸，一般要求小于 1，界面附近最好小于 0.5。Eulerian 模型里，时间步长可以大一些，但也要保证每个时间步内相间作用力不会导致体积分数突变。

我一般先用一个较大的时间步跑几十步看残差，如果残差震荡就减半，直到残差平稳下降。库朗数可以在 Fluent 里监控，如果超过 2 就说明时间步太大了。

# 监控库朗数 solver.tui.solve.monitors.residual.plot("yes") solver.tui.solve.monitors.courant_number.plot("yes") # 自适应时间步：根据库朗数自动调整 solver.tui.solve.set.transient_controls.adaptive_time_step("yes") solver.tui.solve.set.transient_controls.courant_number("1.5") solver.tui.solve.set.transient_controls.time_step_size_min("1e-6") solver.tui.solve.set.transient_controls.time_step_size_max("0.01")

自适应时间步是个好东西，但库朗数上限不要设太大，1.5 到 2 之间比较稳妥。最小时间步设 1e-6 是为了防止时间步被压到零导致计算卡死。

5. 多相流仿真的避坑清单：五条血泪经验

这一章不讲新东西，专门讲我踩过的坑。每一条都是现象、原因、解决三件套，照着排查能省不少时间。

5.1 残差震荡不收敛，先查相间作用力系数

现象：残差曲线在 1e-2 附近上下震荡，跑了几千步也不下降。

原因：相间作用力系数设得太大或太小，导致相间动量交换过强或过弱。比如曳力系数设大了，分散相被连续相拖得太紧，速度场剧烈震荡。

解决：先把所有相间作用力关掉，只留曳力，跑几百步看残差。如果稳定了，再逐个开启其他力，每开一个观察残差变化。找到导致震荡的那个力，把它的系数调小或换模型。

5.2 体积分数出现负值，检查截断和初始化

现象：计算过程中报错，提示某相的体积分数为负。

原因：体积分数方程的数值耗散导致负值，或者初始化时相分布不合理。

解决：开启体积分数截断，下限设 1e-5。检查初始化，确保每一相的体积分数在 0 到 1 之间。如果用了 patch，确认 patch 区域没有重叠。

5.3 界面糊成一片，网格和时间步都要查

现象：VOF 模型跑出来的界面不是锐利的，而是很宽的渐变带。

原因：网格分辨率不够，或者界面压缩因子太小，或者时间步太大。

解决：先加密界面附近的网格，至少保证界面过渡带内有 3 到 5 层网格。然后提高界面压缩因子到 0.5 试试。如果还不行，减小时间步，让库朗数降到 0.5 以下。

5.4 气泡全部沉底或颗粒全部吹飞，曳力模型选错了

现象：Eulerian 模型里，气泡不上升反而沉底，或者颗粒不沉降反而被吹走。

原因：曳力模型和分散相的实际形态不匹配，导致曳力方向或大小错误。

解决：检查曳力模型是否适合当前的分散相。气泡用 Grace 或 Tomiyama，颗粒用 Schiller-Naumann 或 Morsi-Alexander。检查气泡直径或颗粒直径是否设对，这个参数直接影响曳力大小。

5.5 计算几分钟就发散，松弛因子先降一半

现象：计算刚开始几步就发散，残差直接飙到 1e+3 以上。

原因：松弛因子太大，或者初始场太差。

解决：把所有松弛因子降到默认值的一半，跑几百步稳定后再逐步提高。如果还是发散，检查边界条件是否合理，特别是压力出口和速度入口的搭配。

提示：发散的时候不要急着改网格。先降松弛因子，再查时间步，最后才动网格。这个顺序能解决八成问题。

6. 从跑通到跑准：多相流模型的验证与进阶技巧

跑通一个多相流算例只是第一步，跑准才是目的。验证这件事，我一般分三层：网格无关性、时间步无关性、实验数据对比。前两层是自证，第三层是实证。

网格无关性不是把网格加密一倍看结果变不变，而是要看关键物理量。比如气泡上升的终端速度、溃坝前沿的推进距离、流化床的压降。这些量在网格加密后变化小于 5%，才能说网格无关。时间步无关性同理，看的是关键物理量的时间历程，不是残差曲线。

实验数据对比是最难的，因为多相流的实验数据本身就不好找，而且边界条件很难完全复现。我的做法是找经典算例，比如溃坝的 Martin 和 Moyce 实验数据，或者气泡上升的 Grace 经验关联式。这些数据在文献里能查到，用来验证模型设置是否合理。

进阶技巧方面，群体平衡模型（PBM）是绕不开的。如果分散相的粒径分布很宽，单一粒径假设就不成立了。PBM 把粒径分成若干组，每组求解一个体积分数方程，相间作用力按组计算。代价是计算量成倍增加，但相分布预测会准很多。

另一个技巧是自适应网格细化（AMR）。VOF 模型里，界面附近网格自动加密，远离界面的地方网格自动粗化。这样既保证了界面分辨率，又控制了计算量。AMR 的触发条件一般用体积分数梯度，梯度大的地方加密。

# 开启自适应网格细化，基于体积分数梯度 solver.tui.solve.adaptation.control.enable("yes") solver.tui.solve.adaptation.control.criteria("volume_fraction") solver.tui.solve.adaptation.control.refine_threshold("0.1") solver.tui.solve.adaptation.control.coarsen_threshold("0.01") solver.tui.solve.adaptation.control.max_level("3") solver.tui.solve.adaptation.control.interval("10")

这段设置的意思是：每 10 个时间步做一次自适应，体积分数梯度大于 0.1 的网格加密，小于 0.01 的网格粗化，最多加密 3 层。加密阈值和粗化阈值之间要有足够的间隔，否则网格会反复加密粗化，浪费计算资源。

最后说一个我自己的习惯：每次开新算例，先用最粗的网格、最大的时间步、最简的模型跑一个基线，确认流程能走通，再逐步加复杂度。这个习惯让我避免了很多“一上来就调精细参数，结果发现是边界条件写错了”的尴尬。多相流仿真没有后悔药，每一步都留个能回退的版本，比什么都重要。

希望帮到你。

本文还有配套的精品资源，点击获取