4. 流动仿真的物理基础
底部填充胶（Underfill）在倒装芯片（Flip Chip）封装中的流动行为，本质上是微尺度下的粘性流体在毛细力驱动下的多相流问题。要准确仿真这一过程，必须建立基于连续介质力学的控制方程体系。本节将系统阐述Navier-Stokes方程的基本形式、毛细力驱动流动的物理模型，以及它们在底部填充仿真中的简化与适用条件。
4.1 流体力学基本方程：Navier-Stokes方程组
对于不可压缩的牛顿流体（如大多数环氧树脂型底部填充胶），其流动行为由质量守恒方程（连续性方程）和动量守恒方程（Navier-Stokes方程）共同描述。
4.1.1 连续性方程（质量守恒）
在直角坐标系下，不可压缩流体的连续性方程可写为：
∂u/∂x + ∂v/∂y + ∂w/∂z = 0其中 u、v、w 分别为 x、y、z 方向的速度分量。该方程表明，流体微元在单位时间内流入与流出的质量差为零，即体积流量守恒。
4.1.2 动量守恒方程（Navier-Stokes方程）
对于粘性不可压缩流体，Navier-Stokes方程的矢量形式为：
ρ (∂V/∂t + V·∇V) = -∇p + μ∇²V + ρg + F_s其中各项物理意义如下表所示：
| 符号 | 物理含义 | 单位 | 在底部填充中的角色 |
|---|---|---|---|
| ρ | 流体密度 | kg/m³ | 决定惯性力大小，微尺度下影响较小 |
| ∂V/∂t | 局部加速度（非定常项） | m/s² | 填充过程为瞬态，此项不可忽略 |
| V·∇V | 对流加速度（惯性项） | m/s² | 低雷诺数下可忽略（Stokes流） |
| -∇p | 压力梯度力 | Pa/m | 驱动流动的主要体积力之一 |
| μ∇²V | 粘性扩散力 | N/m³ | 微尺度下主导流动的阻力项 |
| ρg | 重力体积力 | N/m³ | 水平填充时可忽略，垂直填充时需考虑 |
| F_s | 表面张力源项（毛细力） | N/m³ | 底部填充的核心驱动力 |
0 = -∇p + μ∇²V + F_s这一简化大幅降低了计算复杂度，是大多数底部填充CFD仿真的理论基础。
4.2 毛细力驱动流动模型
底部填充胶的流动并非依靠外部泵压，而是由弯月面处的毛细压力差驱动。毛细力的大小由Young-Laplace方程描述，而流动速率则受Washburn方程控制。
4.2.1 Young-Laplace方程：毛细压力的来源
当液体与固体表面接触时，由于表面张力的存在，气-液界面会形成弯曲的弯月面。弯月面两侧的压力差（毛细压力）为：
ΔP_cap = σ (1/R₁ + 1/R₂) = σ · κ其中 σ 为液体表面张力系数，R₁、R₂ 为弯月面的两个主曲率半径，κ 为平均曲率。在平行平板间隙（芯片与基板之间）中，若接触角为 θ，间隙高度为 h，则毛细压力可近似为：
ΔP_cap ≈ (2σ cosθ) / h4.2.2 Washburn方程：流动前沿的推进
Washburn方程描述了毛细力驱动下液体在平行平板或毛细管中的流动前沿位置随时间的变化关系。对于底部填充的平行平板模型：
L(t) = √[ (h σ cosθ) / (3μ) · t ]其中 L(t) 为 t 时刻的流动前沿距离，h 为间隙高度，μ 为流体动力粘度。该方程揭示了几个关键规律：
- L ∝ √t：流动前沿位置与时间的平方根成正比，填充速度随时间逐渐减慢。
- L ∝ √h：间隙越大，相同时间内流动距离越远。
- L ∝ √(σ cosθ / μ)：表面张力越大、接触角越小（润湿性越好）、粘度越低，填充越快。
4.2.3 考虑焊球阵列的等效流动模型
实际Flip Chip封装中存在密集的焊球阵列（Bump Array），它们对流动产生显著的阻碍效应。通常采用多孔介质模型或等效阻力模型来近似：
μ_eff = μ · (1 + α · φ)其中 μ_eff 为等效粘度，φ 为焊球的面积占空比，α 为经验形状因子（通常取 2~5）。
4.3 多相流界面追踪方法
底部填充流动涉及气-液两相界面（空气与胶水）的演化。在CFD仿真中，常用的界面追踪方法包括：
| 方法 | 原理 | 优点 | 缺点 | 底部填充适用性 |
|---|---|---|---|---|
| VOF（流体体积法） | 通过体积分数 α 追踪界面 | 质量守恒好，适合复杂拓扑变化 | 界面锐利度不足，需高阶格式 | 最常用，推荐 |
| Level Set（水平集法） | 通过距离函数 φ 追踪界面 | 界面光滑，曲率计算准确 | 质量守恒较差，需重新初始化 | 适合毛细力精确计算 |
| Phase Field（相场法） | 基于扩散界面理论 | 物理基础严格，可处理接触角 | 计算量大，界面厚度非物理 | 学术研究常用 |
F_s = σ κ ∇α其中 κ 由界面曲率计算，∇α 为体积分数梯度，仅在界面处非零。
4.4 底部填充流动的物理模型总览
下图总结了底部填充流动仿真的物理模型层次与关键方程之间的逻辑关系：
4.5 关键无量纲参数
在底部填充流动仿真中，以下无量纲参数用于判断流动状态和主导物理机制：
| 无量纲数 | 定义 | 典型值范围 | 物理意义 |
|---|---|---|---|
| 雷诺数 Re | ρUL/μ | 10⁻⁴ ~ 10⁻² | 惯性力 / 粘性力 → 极低，粘性主导 |
| 毛细数 Ca | μU/σ | 10⁻³ ~ 10⁻¹ | 粘性力 / 表面张力 → 毛细力主导 |
| 邦德数 Bo | ρgL²/σ | < 0.1 | 重力 / 表面张力 → 重力可忽略 |
| 韦伯数 We | ρU²L/σ | < 10⁻⁴ | 惯性力 / 表面张力 → 惯性可忽略 |
(1) 可关闭重力项（Bo < 0.1）；
(2) 可忽略惯性项（Re << 1, We << 1）；
(3) 必须精确求解表面张力项（Ca 适中，毛细力是主要驱动力）；
(4) 时间步长需满足 Courant 数条件：Co = UΔt/Δx < 0.25，以保证界面追踪稳定。
4.6 本章小结
- 控制方程：不可压缩Navier-Stokes方程是底部填充流动仿真的基础，在低雷诺数下可简化为Stokes方程。
- 毛细驱动力：由Young-Laplace方程描述，ΔP_cap ∝ σcosθ/h，间隙越小毛细力越大。
- 流动速率：Washburn方程给出 L ∝ √(hσcosθ/μ · t)，是解析估算填充时间的经典模型。
- 界面追踪：VOF+CSF方法是工程仿真的主流选择，需注意Courant数控制和界面锐利度。
- 焊球效应：必须通过等效粘度或显式几何建模考虑焊球阵列的流动阻力。