6.1 欧拉-欧拉框架
文档摘要
6.1 欧拉-欧拉框架 多相流地图的第一块大陆：所有相都当连续介质、都在同一套网格上解各自的输运方程。本节讲清 VOF、Mixture、Eulerian 三种欧拉方法对"界面"的不同哲学、各自的算力代价与适用地盘，并用流化床案例演示 Eulerian 框架的完整配置思路。6.2 节的拉格朗日路线将以它为参照系登场。 界面哲学的三种答案 欧拉框架的出发点是把分散相也看成充满空间的连续流体，各相体积分数（相含率）加起来等于一。三种方法对"界面与相间作用"给出三种答案： VOF（流体体积法）。为界面分明的流动而生：它解一套动量方程、追踪体积分数场的输运，界面就是体积分数从一到零的等值面。它擅长自由面、晃动、射流破碎初期、气泡大到能被网格刻画的场景。
6.1 欧拉-欧拉框架
多相流地图的第一块大陆：所有相都当连续介质、都在同一套网格上解各自的输运方程。本节讲清 VOF、Mixture、Eulerian 三种欧拉方法对"界面"的不同哲学、各自的算力代价与适用地盘，并用流化床案例演示 Eulerian 框架的完整配置思路。6.2 节的拉格朗日路线将以它为参照系登场。
界面哲学的三种答案
欧拉框架的出发点是把分散相也看成充满空间的连续流体，各相体积分数（相含率）加起来等于一。三种方法对"界面与相间作用"给出三种答案：
VOF（流体体积法）。为界面分明的流动而生：它解一套动量方程、追踪体积分数场的输运，界面就是体积分数从一到零的等值面。它擅长自由面、晃动、射流破碎初期、气泡大到能被网格刻画的场景。它的边界条件也很明确：界面分辨精度直接受网格限制，界面厚度小于网格尺寸的物理（微小气泡、液滴）它看不清。
Mixture（混合模型）。不求界面只求平均：把各相当成互相穿透的连续体，解混合物的动量方程加滑移速度关系式，每个单元算个"平均相含率"。它比 Eulerian 便宜一个量级，适合颗粒或气泡较小且浓度适中的场景——旋流分离器的固液分级初算、沉降过程，都是它的地盘。
Eulerian（欧拉模型）。最贵也最全：每相各自一套动量方程与连续性方程，相间通过界面动量交换项（曳力、升力、湍流扩散力等）耦合。浓相颗粒流、流化床、鼓泡塔这类相间作用强烈、相含率高的场景只有它能正经描述。
图 6-1：三种欧拉方法的界面观
相间作用项：Eulerian 框架的配料表
Eulerian 模型的精度几乎全押在相间作用项上。曳力是最重的一项：Schiller-Naumann 适用于刚性球形颗粒低浓度，Ishii-Zuber 与 Gidaspow 各自覆盖气泡群与流化床的浓相区，Tomiyama 对界面变形气泡更准。升力在旋流气泡场景不可省，湍流扩散力决定颗粒会不会被湍流涡甩出主流区，壁面润滑力管气泡别贴壁扎堆。这些子模型各有经验系数与适用范围——工程上正确姿势是：按工况挑组合，用小型算例与文献或实验对标后再放大到整机。模型关的验收语言：曳力模型选型依据必须可查。
案例：鼓泡塔的气含率预测
背景。某化工厂鼓泡塔放大设计，需要预测不同表观气速下的气含率（气体占体积比）分布与循环流型。塔径两米，孔板布气，气泡尺寸毫米级——界面尺度远小于网格可分辨尺度，VOF 无望，定位到 Eulerian 框架。
操作。气液两相 Eulerian 配置：液相为主相；曳力选 Ishii-Zuber 加 Tomiyama 升力与湍流扩散力；气泡直径按孔板与气速关联式给定分布；湍流用 k-epsilon 加分散相湍流调制。先跑一维柱塞段对标文献气含率关联式，确认曳力组合合理，再上三维全塔。
结果。表观气速扫描四个工况，塔内平均气含率与工厂标定数据偏差在百分之十以内，预测出底部密相区与壁面下行循环流，为放大设计提供了取样口与内构件的位置依据。
解读。案例展示了"两级验证"工作法：先用一维或小算例验证相间作用模型组合，再上全尺寸。直接上全尺寸一旦偏差，你分不清是网格、湍流还是曳力的锅。分层锁定误差源，是多相仿真区别于单相的核心纪律。
变式。若颗粒是固相且浓度极高（流化床），曳力换 Gidaspow、附加摩擦应力模型，并评估是否升级到 DDPM 或 DEM 耦合以获得颗粒级信息；若气泡尺寸谱宽，考虑群体平衡模型（PBM）追踪尺寸分布而不是固定直径。
欧拉配置问答与选型对照
VOF 里的界面锐化要开吗。界面追踪精度与数值振荡是一对老冤家：几何重构类界面捕捉在库朗数受控时值得开启；开启后界面附近速度出现毛刺，先降库朗数再谈关掉锐化。
Mixture 什么时候会不够用。滑移关系式失真的场景：颗粒尺寸分布宽、相间速度差剧烈、需要颗粒浓度场的定量结论——这些信号出现就升级 Eulerian，别在 Mixture 里硬调参数。
Eulerian 的气泡直径怎么给。固定直径是最粗近似；关联式直径（按当地湍流与物性估算）是工程默认；尺寸谱重要时上群体平衡模型。三个档位的代价递增，先问结论对直径分布敏不敏感。
三种欧拉方法一分钟决策卡：
| 你的任务 | 首选 | 升级信号 |
|---|---|---|
| 自由面、晃动、充放液 | VOF | 界面毛刺先调库朗数 |
| 气泡小而密、只看场量 | Mixture | 需要浓度定量结论时升 Eulerian |
| 流化床、鼓泡塔、强相间作用 | Eulerian | 需要颗粒级细节升 DDPM 或 DEM |
| 大气泡可被网格刻画 | VOF 或 Eulerian | 按界面分辨率需求定 |
关于欧拉框架的算力与精度权衡，补一组量级感受：同一几何同一网格，Mixture 的每步耗时约为单相的一点五到两倍，Eulerian 双相约为两到三倍，VOF 因界面锐化与时间步约束常再翻倍。但更贵的框架不代表更准的结论——地盘错位时，贵的反而更糟。这组数字的用途是预算沟通：在项目计划里给多相任务留出正确的算力与工时余量，比开始后返工便宜得多。多相项目的常规预算经验：网格与前处理工时按单相的一点五倍计，计算时间按所选框架的倍率计，验证工时按两倍计（敏感性维度更多）。
最后补充湍流与多相的接口认知，因为它是欧拉框架里第二贵的误差源：分散相反过来影响连续相湍流（湍流调制），颗粒小而多会削弱湍流、大气泡会增强湍流生成。Fluent 提供分散相湍流模型开关，工程上的正确姿势是：先关掉调制算基准，再做开关对照评估其影响量级——影响超过结论容差就保留并写进模型声明。多相任务的湍流参数因此比单相多一层：既要有连续相湍流的口径说明，又要有调制开关的影响评估，两者都进验收单的模型栏。多相仿真的可信度就是这样一层层叠出来的，没有捷径，但有清单。
本节验收点
- 框架选择有相含率与界面尺度的量化依据，VOF 与 Eulerian 不混用错位。
- 相间作用项组合已列表，每个子模型有适用范围说明。
- 相间模型组合做过独立小算例对标，误差源分层锁定。
- 气泡或颗粒尺寸的取值与分布有出处（关联式或测量）。