在绝大多数宏观流动模拟里，我们都默认“壁面无滑移”——流体微团一旦贴上壁面，就只能乖乖跟着壁面走，速度差为零。然而当流动尺度缩小到微米、纳米，或者气体稀薄到平均自由程与特征尺度可比拟时，壁面分子与流体分子之间的动量交换不再“完美”，于是流体在壁面处会出现一层可察觉的切向速度差，这种现象就叫壁面滑移（wall slip）。

它本质上是一条“边界处的速度不连续”，其大小通常用滑移长度 b 来量化：
u_slip = b · (∂u/∂n)|wall
其中 b 具有长度量纲，可视为“把速度剖面线性外推到零”时与壁面的距离。滑移长度越大，壁面给流体的“摩擦力”就越小，宏观上表现为减阻、流量增大或压降降低。

CFD 计算中如何把“滑移”写进代码？

1. 边界条件替换 在动量方程离散后的线性系统里，无滑移相当于把壁面节点速度直接钉死（Dirichlet）。滑移则把该节点速度作为未知量，把滑移长度或剪切应力写进矩阵的一行，变成 Robin 型条件： μ·(∂u/∂n) = C·(u_wall − u_fluid) 其中 C = μ/b，b 为滑移长度。
2. 近壁网格处理 若采用低 Re 湍流模型或层流解析，壁面第一层网格高度 Δy 应保证 y⁺≈1，否则滑移量会被数值黏性“吃掉”；若用大涡/雷诺平均，则需在壁面函数里把 Δu⁺ 的修正项显式写出，例如对数律叠加 Δu⁺(b⁺)。
3. 多尺度耦合 对微沟槽、多孔涂层等场景，可先用格子-玻尔兹曼或 MD 算出 b 与宏观参数（剪切率、压力、温度）的数据库，再用神经网络拟合成“滑移边界函数”，在宏观 RANS/LES 求解器里以自定义边界条件的方式调用。

COMSOL 里能用的“滑移”套餐
COMSOL 把滑移现象拆成三条主线，分别在“层流”、“滑移流”和“稀布/分子流”接口中开箱即用：

1. 通用“滑移”边界 位置：层流 → 壁 → 滑移 数学形式： τ_t = − (μ/b)·(u_t − u_wall) 用户只需填一个“滑移长度” b，可设为常数、表达式或随场变量函数。若 b→0，自动退化为无滑移；若 b→∞，则退化为完全滑移（剪切应力为零）。 适用：微通道液体、聚合物挤出、低渗透多孔壁等。
2. 增强 Navier 滑移（用户自定义） 当滑移长度与剪切率本身相关时（如粘弹性熔体），可在“弱贡献”里写入： b = b₀·|γ̇|^(n-1) 把 τ_t = − (μ/b)·u_t 作为附加弱形式积分到壁面单元上，实现剪切依赖滑移。
3. 滑移流接口（Kn 0.001–0.1） 位置：稀有fied flow → Slip flow 壁面自动加载 Maxwell 热-质滑移公式： u_s = σ_p·λ/μ·τ_t + σ_T·μ/(ρT)·∇T_t T_gas − T_wall = ζ_T·λ·(∂T/∂n) 其中 λ 为平均自由程，σ_p、σ_T、ζ_T 由“切向动量调节系数”α 计算而来，材料库已内置空气、氮气等参数。 该接口同时把连续性、动量、能量方程在克努森层外求解，省去解析纳米级边界层的网格负担。
4. 分子流/过渡流（Kn>0.1） 若通道尺寸继续缩小，COMSOL 提供“分子流”模块，采用线-of-sight 或 DSMC 方法，把壁面视为漫反射/镜面反射表面，滑移效应已隐含在散射核里，无需再额外指定滑移长度。

小结
壁面滑移不是“玄学”，而是把分子尺度动量交换的微观信息压缩进一个“可测长度 b”的边界参数。只要你给出合理的 b（实验、理论或跨尺度计算），在 COMSOL 中只需一行滑移长度，就能把微纳尺度的减阻、流量增强或温度跃变现象搬到宏观求解器里；当克努森数再往上走，就换用更底层的滑移流或分子流接口，让 Maxwell 或 DSMC 帮你自动处理那层“不听话”的边界。祝你在下一次仿真里，既能抓住边界层，又不被边界层“咬”住网格！

引言：一个常见的困惑

当你在 COMSOL 中进行后处理，创建一维绘图组（1D Plot Group）并添加线图（Line Graph）时，X 轴数据（x-Axis Data）下拉菜单中会出现一个默认值：Arc length（弧长）。

很多用户会困惑：我画的是一条直线段，哪来的”弧”？这个弧长到底指什么？

如果你也曾盯着横坐标上的”弧长”标签感到疑惑，这篇文章将彻底澄清这个概念。

---

一、直观理解：弧长即”沿线行走的距离”

在 COMSOL 中，**弧长（Arc Length）**最朴素的定义是：

从你设定的起点开始，沿着所选几何路径累积的实际路径长度。

正如密歇根理工大学 CFD 课程教程中的经典解释：

“Arc Length means for the cut line you selected, move along the contour of the line, wherever it goes.”

（”弧长”意味着对于你选中的切割线，沿着这条线的轮廓移动，无论它通向何处。）

关键特征：

• 物理意义明确：横坐标的数值直接对应实际空间中的距离（单位通常是 mm 或 m）
• 路径跟随性：即使你的”切割线”（Cut Line）是弯曲的，弧长也会沿着曲线累积，而非简单的 x 或 y 投影
• 起点归零：弧长在路径的起始点始终为 0，向终点方向递增

---

二、数学本质：弧长参数化

从微分几何角度，COMSOL 中的弧长对应弧长参数化（Arc-length Parameterization）。

对于空间中任意曲线 $\vec{r}(t)$r(t)，其弧长 $s$s 定义为：$$s(t) = \int_{a}^{t} \|\vec{r}'(\tau)\| \, d\tau$$s(t)=∫at​∥r′(τ)∥dτ

其中：

• $\vec{r}(t)$r(t) 是曲线的参数方程
• $\|\vec{r}'(\tau)\|$∥r′(τ)∥ 表示速度向量的模长（即瞬时切向量的长度）
• 积分从起点 $a$a 累积到当前参数 $t$t

当使用弧长 $s$s 作为参数时，曲线满足单位速度条件：$$\|\frac{d\vec{r}}{ds}\| = 1$$∥dsdr​∥=1

这意味着：弧长坐标每增加 1 米，你就确切地在空间中移动了 1 米——没有畸变，没有缩放。

---

三、COMSOL 中的具体实现

1. 何时会出现”弧长”？

在以下数据集的线图中，X 轴可选择 Arc length：

• Cut Line 2D / 3D（切割线）
• Edge（几何边）
• Parametric Curve（参数化曲线）

2. 方向控制

COMSOL 提供了两个相关选项：

• Arc length：从起点到终点，沿曲线方向递增
• Reversed arc length：反向，从终点指向起点

这允许你自由控制绘图的方向，确保数据流符合物理直觉（例如从入口到出口）。

3. 与”表达式”（Expression）的区别

你可以将 X 轴从 Arc length 切换为 Expression（如 x, y, z 或 sqrt(x^2+y^2)），这实质上是投影到某个坐标轴或自定义度量上。而 Arc length 保持内禀（intrinsic）特性——它只关心你在曲线上走了多远，不关心曲线在空间中的朝向。

---

四、为什么使用弧长？三大优势

1. 几何无关性

无论切割线是直线、圆弧还是任意样条曲线，弧长坐标提供了统一的度量基准。你可以比较两条不同形状边界上的压力分布，只要它们的弧长范围相同。

2. 物理直观性

在流体力学中，沿壁面的距离；在传热中，沿散热器翅片的长度——这些过程量天然以路径长度衡量，而非笛卡尔坐标。

3. 后处理灵活性

当你使用附加平行线（Additional parallel lines）生成一系列切割线时，每条线的几何坐标系不同，但弧长坐标将它们统一到了相同的局部 1D 坐标系中，便于批量比较。

---

五、实战示例

示例 1：弯曲管道中的压力降

假设你建模了一个 S 形弯管，想查看沿管壁的压力分布：

• 使用 Cut Line 3D 沿管中心线创建切割线
• 选择 Arc length 作为 X 轴，Y 轴选择 p（压力）
• 所得曲线完美展示压力随流经距离的变化，而非随 x 坐标的变化

示例 2：边界层速度剖面

在平板流动中，你创建垂直于壁面的多条切割线查看速度边界层：

• 若选择 y（垂直坐标），不同 x 位置的曲线因当地网格不同难以对齐
• 若选择 Arc length，所有曲线都从 0（壁面）开始，便于叠加比较无量纲化剖面

---

六、进阶技巧与避坑指南

技巧 1：归一化弧长

有时你想比较不同长度的边界（如短翅片 vs 长翅片），可在 X 轴表达式中使用：

其中 s 是弧长变量，length 是几何参数或积分值，实现归一化位置 。

技巧 2：与几何变量结合

COMSOL 内部使用变量 edgparal 表示弧长参数，edgparnal 表示归一化弧长参数（0 到 1）。在自定义表达式中可直接调用。

⚠️ 避坑：闭合环路的起点

对于闭合边界（如圆），COMSOL 会自动选取一个起点（通常是参数化起点）。如果你发现弧长图出现”断层”或突变，检查切割线是否跨越了几何接缝。

---

总结

COMSOL 一维绘图组中的弧长，本质上是沿指定路径的内禀距离坐标。它不是软件的随意设定，而是微分几何中弧长参数化的工程实现。

理解这一点，你就能：

• 自信地解释横坐标的物理意义
• 灵活选择 Arc length vs. 坐标投影
• 处理复杂几何路径的后处理结果

下次当你看到弧长从 0 延伸到 20 mm 时，你知道：这不是几何在”弯曲”，而是数据在沿着你关注的真实路径展开。

---

参考阅读：

• COMSOL Reference Manual: Line Graph Node
• MTU CM3110 CFD Lab Notes: Arc Length Explanation
• COMSOL Forum: “What does Arc Length mean” Discussion

01 先讲人话：什么叫“分散型”？

把咖啡倒进牛奶，一瞬间出现无数小液滴——它们不连续，却多到肉眼无法追踪单个界面。

COMSOL 把这种“离散相=无数小微粒”的场景统称为分散型两相流（Dispersed Multiphase Flow）。

核心思路：不画颗粒边界，只用“体积分数”α 描述空间占比，把两相当成“互相穿透的连续介质”来算，省时 90 %。

---

02 四兄弟亮相：一张选型图带走

接口	颗粒浓度	是否考虑滑移	颗粒-颗粒碰撞	计算量	典型场景
气泡流 Bubbly Flow	α < 10 %	连续场滑移	无	★	鼓泡塔、曝气池
混合物模型 Mixture	α < 10 %	连续场滑移	无	★	砂水泵、泥浆
欧拉-欧拉 Euler–Euler	α 任意	双速度场	可选	★★★	流化床、高浓浆
分散相模型 DPM	α < 10 %	单个粒子	可选	★★	喷雾、除尘、药片包衣

---

03 逐条拆解：它们到底差在哪？

① 气泡流——“给气泡集体发通行证”

• 只解一套动量方程，气泡当“第二相”拖在连续相里跑。
• 自带阻力+升力+壁面润滑，气泡直径可随传质/聚并实时变。
• 坑点：体积分数>10 % 时，连续相“被气泡反推”效应被低估，压降会飘。

② 混合物模型——“把两相绑成一根绳”

• 假设两相速度差=滑移速度，只求混合物动量方程，省内存 30 %。
• 可顺手打开“剪切诱导迁移”，模拟粗颗粒在弯管里的分离。
• 坑点：颗粒沉降速度必须提前给定，若粒径分布宽，误差会被放大。

③ 欧拉-欧拉——“给每相单独发驾照”

• 真正双速度场+双湍流，颗粒相也能有自己的 k-ε。
• 可开“颗粒相压力”+“摩擦-动能理论”，把堆积、流化、喷动一次打包。
• 坑点：方程数翻倍，网格>100 万时，64 GB 内存只是“起步价”。

④ DPM——“真的一个一个算”

• 把颗粒当质点，用牛顿第二定律拉格朗日追踪，百万颗粒也能画直方图。
• 可与任意连续场耦合，做喷雾干燥、静电除尘、药片包衣最直观。
• 坑点：体积分数>10 % 时，颗粒-颗粒碰撞、屏蔽效应全漏掉，结果会“过于乐观”。

---

04 三步选型口诀（背下来就能装大佬）

1. 看浓度： α < 10 % → ①②④ 随便挑；α > 10 % → 直接 ③ 欧拉-欧拉。
2. 看颗粒数： 颗粒数 < 10⁵ 且想追踪“单个”→ ④ DPM；否则 ①②③ 走连续场。
3. 看预算： 笔记本 16 GB 内存 → ①②；工作站 128 GB → ③ 随便玩；集群 → ④ 百万粒子也 OK。

---

05 30 秒案例对照

场景	推荐接口	关键设置
鼓泡塔（气含率 5 %）	气泡流	气泡直径 3 mm，升力系数 0.5
砂水输送（体积 8 %）	混合物	滑移速度模型选“Schiller-Naumann”
流化床（固含率 40 %）	欧拉-欧拉	打开颗粒相压力+Syamlal-O’Brien 曳力
农药喷雾（粒径 50 µm）	DPM	入口射流 1×10⁵ 颗粒/s，曳力选“Stokes-Cunningham”

---

06 写在最后

分散型两相流的核心只有一句话：“颗粒太多，界面不画，用体积分数代替。”

记住四兄弟的名字和浓度分界线，下次再看到“鼓泡塔+砂水+流化床”的混合需求，5 秒钟就能给出方案。

祝你颗粒不飞、气泡不炸、收敛条一路飘绿！

关键词：能斯特–普朗克、菲克定律、电迁移、扩散、COMSOL、电化学、传质

---

一、写在前面：为什么这两个名字总一起出现？

在电化学、膜分离、微流控乃至锂电池、燃料电池等多物理场耦合问题里，我们几乎总会看到两条“传质路线”：

1. 菲克定律（Fick’s Law）——“无外力”时的纯扩散；
2. 能斯特–普朗克方程（Nernst–Planck, NP）——“有电场”时的扩散+迁移。

它们既像“师生”又像“搭档”：菲克定律是 NP 方程在零电场下的特例；NP 方程则是菲克定律在带电离子世界里的“升级版”。理解这一点，是正确在 COMSOL 中选择接口、简化模型、解释结果的第一步。

---

二、物理图像：一条通量表达式的“进化史”

1. 菲克第一定律（1855）

Jdiff=−D∇c\mathbf{J}_\mathrm{diff} = -D\nabla cJdiff=−D∇c

• 仅由浓度梯度 ∇c 驱动
• 适用于电中性分子或“电场可忽略”的稀溶液区域
• 在 COMSOL 里对应 “稀物质传递” 接口，关闭电迁移项即可

2. 能斯特–普朗克方程（1890）

Ji=−Di∇ci−ziFRTDici∇ϕ+ciu\mathbf{J}_i = -D_i\nabla c_i – \frac{z_i F}{RT}D_i c_i\nabla\phi + c_i\mathbf{u}Ji=−Di∇ci−RTziFDici∇ϕ+ciu

• 三项依次为 扩散、电迁移、对流
• 电荷数 z_i 把电场 ∇ϕ 的作用量化进来
• 在 COMSOL 里对应 “电化学→能斯特–普朗克” 或 “稀物质传递+静电场” 耦合

3. 关系小结

场景	主导方程	COMSOL 实现
无电场、低浓度	菲克定律	稀物质传递（关闭电迁移）
有电场、稀溶液	NP 方程	稀物质传递 + 静电 或 腐蚀模块→NP
电场极强、电双层	修正 NP + Poisson	电渗流接口、EDL 边界条件

---

三、COMSOL 中的“三档”建模策略

1. 初级：纯扩散——“ Electroanalysis ”接口

• 默认关闭迁移项，仅求解菲克第二定律
• 适用于旋转圆盘电极、微电极计时安培法等“电流小到电场可忽略”的场景
• 边界直接给浓度或通量，求解速度快，易收敛

2. 中级：迁移+扩散——“Secondary Current Distribution + NP”

• 电解质电位 ϕ_l 与浓度 c_i 双向耦合
• 电导率 σ_l 随局部离子强度自动更新
• 适用于腐蚀、电镀、液流电池，其中电场显著但浓度变化不大

3. 高级：全耦合——“电渗流/电双层”模型

• 引入 Poisson 方程求解空间电荷，获得 nm 级双电层
• 流体侧耦合 Navier–Stokes，壁面采用 Helmholtz–Smoluchowski 滑移速度
• 用于微流控芯片、纳滤膜、电渗泵

---

四、实战案例 1：纳米孔道中的离子选择性

问题：10 nm 单锥形纳米孔，两端电压 1 V，KCl 浓度 10 mM，如何计算 K⁺ 与 Cl⁻ 的通量差异？

建模要点

1. 几何：2D 轴对称，锥角 10°，尖端半径 5 nm
2. 物理场
   • “静电”求解 ∇²ϕ = –ρ/ε
   • “稀物质传递”添加 K⁺、Cl⁻，扩散系数 1.96×10⁻⁹ m² s⁻¹
   • 耦合方式：电迁移项自动调用 ∇ϕ
3. 边界
   • 左侧：ϕ = 1 V，c = 10 mM
   • 右侧：ϕ = 0 V，c = 10 mM
   • 壁面：无通量（硬墙）
4. 结果
   • 阳离子电流 > 阴离子电流，出现 电流整流
   • 在尖端出现 浓度极化——菲克定律单独无法预测此现象

---

五、实战案例 2：锂电正极 Li⁺ 互扩散+电迁移

问题：LiFePO₄ 颗粒半径 100 nm，表面嵌锂速率由 Butler–Volmer 控制，内部 Li⁺ 化学扩散系数 1×10⁻¹³ m² s⁻¹，如何获得倍率性能曲线？

建模要点

1. 几何：1D 球对称，r ∈ [0, 100 nm]
2. 物理场
   • “稀物质传递”接口仅保留菲克项（固相无电场）
   • 表面通量边界：–J_BV = i₀(c_s/c_max)^0.5(ηF/RT)
   • 颗粒群尺度用 “电池模块→单粒子模型” 封装
3. 结果
   • 5 C 放电容量保持 85 %，与实验吻合
   • 若错误地把液相 NP 方程套用到固相，会得到 错误扩散通量 一个量级

---

六、常见坑与调试技巧

症状	可能原因	解决
浓度出现负值	网格 Peclet 数 > 2	启用人工扩散 / 细化网格 / 使用 DG 方法
不收敛	电导率初始值零	给 σ_l 初值 1 S m⁻¹，或先跑 Secondary Current 分布
电流密度比实验小 10×	忘了 migration 项	检查“电迁移”复选框是否勾选
纳米孔模拟算不动	双电层太薄	使用 EDL 边界条件替代体网格解析

---

七、一张思维导图帮你选接口

浓度梯度？ 是
├─ 电场可忽略？
│  ├─ 是 → 稀物质传递（菲克）
│  └─ 否 → 能斯特–普朗克
│        ├─ 电中性 → 腐蚀/电池→NP
│        └─ 空间电荷 → 电渗流→EDL
└─ 固相扩散 → 固相扩散（菲克即可）

---

八、结语：从物理到代码，再回归物理

菲克定律像一把“瑞士军刀”，简单、直观，却只在“电场故事不精彩”时好用；能斯特–普朗克方程则把“电场”这位主角请回舞台，让传质与电荷真正握手。COMSOL 的价值，正是把这两段物理用同一套网格、同一支求解器串起来，让你在同一块屏幕上看见 扩散边界层 与 电迁移电流 如何此消彼长。

下次再打开模型向导时，不妨问自己三句话：

1. 我的离子感受到电场吗？
2. 浓度变化会反过来改变电场吗？
3. 我需要解析双电层吗？

答案一旦清晰，接口、边界、求解器设置就能“水到渠成”。祝你建模愉快，收敛更快！

---

参考文献与延伸阅读

COMSOL Corrosion Module User’s Guide, 2018.

Modeling Electroosmotic Flow and Electrical Double Layer, COMSOL Blog, 2023.

Lehigh University Thesis, 2023 — 电动力多孔介质传输.

CSDN 电化学传质笔记, 2021.

一篇让你秒懂“瞬态” vs “双向耦合”的爽文

---

01 先讲个故事：

凌晨 2 点，某车企工程师小王盯着屏幕抓狂：

“ 为啥我的排气管里 0.1 µm 的碳烟颗粒模拟得飞起，一到 10 µm 的 PM10 就全员‘撞墙’？！”

他隔壁做流化床的老李悠悠飘过：

“少年，你还没让颗粒‘还手’吧？单向耦合一时爽，双向耦合火葬场。”

——这就是今天的主角：瞬态粒子追踪（Transient）和双向耦合粒子追踪（Bidirectional Coupled）在 COMSOL 里的“江湖恩怨”。

---

02 一张图看懂“谁说了算”

维度	瞬态（单向）	双向耦合
颗粒对 fluid 的“话语权”	0	有，且实时
计算顺序	先算流场→再算颗粒	流场 & 颗粒同步迭代
典型颗粒体积分数	<1 %	1 %–10 %
电脑风扇噪音	笔记本轻哼	机房直升机
一句话人设	“流体是我爸，我只管被带飞”	“我命由我不由流”

---

03 模式 1：瞬态——“渣男”流体的单向奔赴

官方定义：先求解稳态/瞬态流场，再把速度场插值给颗粒，颗粒完全被动 。

COMSOL 操作：

1. 研究 1：稳态 → 层流
2. 研究 2：瞬态 → 粒子追踪

适合场景：

• 空气净化器里的花粉：体积分数 0.01 %，花粉再闹腾也掀不起风。
  • 药物气溶胶在肺部：0.5 µm 液滴， airway 气流说去哪儿就去哪儿。
    • 炫酷动画：把颗粒染成彩虹色，发个朋友圈“看，我的 CFD 求婚钻戒轨迹”

翻车现场： 小王把 5 % 质量负载的碳烟当成“瞬态”算，结果实验测得的压降比仿真高 30 %——颗粒其实集体刹车了气流，他却没让颗粒“还手”。

---

04 模式 2：双向耦合——“夫妻”吵架式实时互殴

官方定义：颗粒对流体施加体积力/动量源，流体再反过来影响颗粒，须迭代到收敛 。

COMSOL 操作：

研究 → 双向耦合粒子追踪（瞬态）

求解器自动循环：

1. 把当前颗粒力“甩”给流体方程
2. 流体更新速度场
3. 用新速度推颗粒
4. 重复 N 次直到“吵不动”

适合场景：

• 流化床锅炉：0.5 mm 石英砂 × thousands，床层膨胀 300 %，全靠颗粒“托”住气流。
• 泥浆泵送：10 % 体积分数的矿渣，颗粒集体“堵管”，压降翻倍。
• “沙尘暴”版显卡散热：沙漠基站进风口，0.2 mm 沙粒把风扇曲线硬生生压成“下垂狗”。

翻车现场： 老李第一次跑双向耦合，时间步长没改——0.01 s 的“大跨步”直接让颗粒穿墙而出，收敛曲线像心电图一样“蹦迪”。秘诀：先 segregated 后 fully coupled，再给牛顿法加点阻尼 。

---

05 实战锦囊：三步选对模式

1. 算体积分数 φv = (颗粒总体积)/(域体积) φv < 1 % → 瞬态，放心飞； 1 % < φv < 10 % → 双向耦合，准备加班； φv > 10 % → 考虑欧拉-欧拉颗粒流或 DEM，别硬刚。
2. 看“惯性指数” 斯托克斯数 St = τp / τf St << 1：颗粒乖乖跟着流； St ≈ 1：颗粒开始“叛逆”； St >> 1：颗粒当流场是空气，必须双向。
3. 问实验 压降、出口浓度、床层膨胀高度——只要实验比仿真“更倔强”，八成是你没让颗粒“还手”。

---

06 彩蛋：把“双向耦合”讲给老妈听

“妈，你知道咱俩为啥吵架不？因为我说一句你顶一句，这就是双向耦合。

要是我自言自语你不回，那叫瞬态单向——我一人说了算，你纯旁听。”

——老妈秒懂，还顺手给你发了 200 块红包买“更好的散热风扇”。

---

07 结语

下次再看到颗粒在 COMSOL 里“鬼畜”或“失踪”，先别拍键盘：

它们也许只是在抗议——“我们不想再做沉默的羔羊，我们想当改变风的方向的那粒沙。”

愿每一位 CFDer 都能在“瞬态”与“双向”之间，找到颗粒与流体最和谐的“婚姻模式”。

“颗粒若只如初见，何来流场空悲切。”
——共勉。

---

参考文献

COMSOL Particle Tracing Module User’s Guide

COMSOL 中国《粒子追踪模块简介》

COMSOL 5.6 粒子追踪模块文档

COMSOL 5.3a Introduction to Particle Tracing

CSDN 博客《Comsol 多孔介质粒子流动模拟与轨迹追踪实战案例》

一、先分清：什么叫“分离”两相流

在 COMSOL 的语境里，只要两相各自占据连续且清晰的空间，中间只出现一条或有限条可辨识的界面，就叫分离（separated）两相流。典型场景：

• 微通道里一段水柱被油柱推着走；
• 液膜沿壁面下滑，气体在中心 core 区域；
• 毛细管出口尚未破裂的液射流。

一旦液滴开始破裂、合并，拓扑结构瞬间爆炸，就滑入“分散”范畴，必须改用水平集/相场。今天的主角是只能留在“分离”区间的三套方法：

1. 水平集（Level Set）
2. 相场（Phase Field）
3. 动网格（Moving Mesh）

---

二、一张表看差异

维度	水平集	相场	动网格
界面厚度	人为设定 ε（通常 2–3 单元）	人为设定 ε + 迁移率 M	零厚度，网格即界面
拓扑变化	✅ 自动破裂/合并	✅ 自动破裂/合并	❌ 不允许（网格撕裂直接报错）
表面张力	连续表面力 CSF	同左，能量泛函更守恒	直接在边界加力，精度最高
守恒性	质量略漂移	能量+质量守恒最好	质量绝对守恒
计算量	中	中	最小（二维可省 50 %）
参数门槛	γ 重初始化参数	M 迁移率要调	几乎零参数，但需手动给接触角
模块归属	CFD / 微流体	CFD / 微流体	仅微流体
典型应用	液滴破裂、射流、T 型乳液	微流控、表面张力主导、需守恒	无破裂：液膜、弯月面、毛细上升

---

三、实战选型 3 句口诀

1. “会断”→ 直接排除动网格 只要预判液柱会断成滴，动网格就会因为“边界无法分裂”而崩溃 。
2. “要守恒”→ 先相场再水平集 相场把界面嵌入自由能泛函，能量守恒最好，适合表面张力 < 0.1 N m⁻¹ 的微尺度 。
3. “只求快”→ 先动网格再转场 二维液膜、无破裂毛细上升，先用动网格跑粗网格，5 min 出结果，再决定是否转场方法 。

---

四、调参黑名单 & 灰名单

方法	灰名单（要调但好办）	黑名单（一错算到崩溃）
水平集	ε = hmax/2，γ = umax	入口边界放在 φ=0.5 区 → 不收敛
相场	ε 同上，M≈2umaxhmax/(3√2σ)	M 过大 → 界面“糊”；M 过小 → 时间步爆炸
动网格	光顺系数 0.3–0.5	接触角 ≠ 几何初始角 → 第 1 步就飞射

---

五、案例 30 秒对比

T 型微通道生成 50 μm 液滴（Ca = 0.01，二维轴对称）

接口	自由度	生成周期误差	CPU 时间	备注
动网格	42 k	6 %	0.4 ×	液滴一旦颈缩→ 报错终止
水平集	87 k	3 %	1 ×	顺利断裂，质量漂移 < 1 %
相场	89 k	2 %	1.1 ×	断裂点与实验吻合最好

---

六、一张思维导图（保存即可）

会破裂？  
├─ 是 → 动网格❌  
│   ├─ 要守恒？  
│   │  ├─ 是 → 相场  
│   │  └─ 否 → 水平集  
└─ 否 → 动网格✅（最快）

---

七、结语

• 先问“会不会断”，再问“要不要守恒”，最后看“算力够不够”，三句话就能选出 COMSOL 分离两相流的正确接口。
• 入口边界永远放在“单相域”，别让初始界面穿过入口，这是 90 % 不收敛的元凶。
• 把这篇导图存成手机壁纸，下次建模直接对照，再也不用翻 Help 文档！

---

参考

LeMoine J., Comparison of Two-Phase Numerical Modelling Techniques in Applications with Electrohydrodynamics, 2023.

jishulink 博客，《COMSOL 两相流之水平集法》，2020.

知乎专栏，《COMSOL 中的多相流建模与仿真：第 1 部分》，2021.