2. 流体-固体-热耦合物理基础

在 Ansys Fluent 的流固耦合与热变形联合仿真中，真正决定结果是否可信的，不是网格有多密、迭代有多快，而是三个物理场是否被正确地“接”在一起。流体域遵循 Navier–Stokes 方程，固体域遵循热弹性力学方程，两者在交界面上必须同时满足运动学、动力学和热学三类守恒条件，热变形则通过界面位移把温度场的影响反馈给流场。本节围绕这条主线，把耦合的物理基础讲清楚。

一、流体域：Navier–Stokes 方程与能量方程

流体域求解的核心是不可压缩（或可压缩）Navier–Stokes 方程。以不可压缩牛顿流体为例，连续性方程与动量方程为：

∇·u = 0
ρ(∂u/∂t + u·∇u) = -∇p + ∇·(μ(∇u + ∇uᵀ)) + ρg

其中 u 为速度矢量，p 为压力，ρ 为密度，μ 为动力黏度。当涉及热变形时，还必须求解能量方程：

ρc_p(∂T/∂t + u·∇T) = ∇·(k∇T) + Φ

这里 T 为温度，c_p 为定压比热，k 为导热系数，Φ 为黏性耗散项。对大多数低速热流固耦合问题，黏性耗散可忽略，但对高剪切流动（如窄间隙润滑、高速密封）则不能随意丢弃。Fluent 中这些方程以有限体积法离散，压力-速度耦合常用 SIMPLE 或 Coupled 格式，动量与能量对流项建议采用二阶迎风格式，以保证温度梯度与速度梯度在交界面附近的分辨率。

二、固体域：热弹性力学方程

固体域不求解流动，而是求解热弹性平衡方程。在准静态假设下，忽略惯性项，控制方程为：

∇·σ + f = 0
σ = C : (ε - ε_th)
ε_th = α(T - T_ref) I

其中 σ 为应力张量，ε 为总应变，ε_th 为热应变，α 为线膨胀系数，T_ref 为无应力参考温度，C 为弹性刚度张量。这里有一个极易出错的约定：热应变的正负号取决于参考温度的定义。若 T > T_ref，则 ε_th 为正，材料膨胀；若 T < T_ref，则收缩。全文及图示均采用“以 T_ref 为零应变状态、升温为正膨胀”的约定，切换参考温度时必须重新说明。

在 Fluent 的流固耦合工作流中，固体域通常由 Ansys Mechanical 求解，Fluent 通过 System Coupling 传递界面温度与热流，Mechanical 返回界面位移与温度。若固体变形量级很小（例如相对特征尺寸小于 1%），可只做单向热-结构映射；若变形会显著改变流道，则必须做双向耦合。

三、流固交界面守恒条件

交界面是耦合的“合同”，必须同时满足三类条件：

• 运动学条件：流体速度等于固体界面速度，无滑移时 u_f = ∂d_s/∂t；
• 动力学条件：流体应力与固体应力在界面法向平衡，σ_f·n = σ_s·n；
• 热学条件：温度连续 T_f = T_s，热流连续 q_f = q_s，即 k_f ∂T/∂n = k_s ∂T/∂n。

实际仿真中，温度连续通常由 System Coupling 自动保证，但热流连续依赖两侧网格与时间步的匹配。若流体侧边界层网格过粗，界面热流会被低估，导致固体温度偏低、热变形偏小。因此交界面附近流体侧建议保证 y+ 在 1 附近或使用增强壁面处理，固体侧在界面处至少布置 3～5 层细化网格。

四、热变形驱动机制与耦合分类

热变形的驱动链可以概括为：流场温度分布 → 界面热流 → 固体温度梯度 → 热应变 → 位移与应力 → 流道几何改变 → 流场再分布。这条链是否闭合，决定了耦合类型的选择。

按耦合强度可分为：单向耦合（只传递温度，不反馈位移）、双向弱耦合（每步交换温度与位移，但固体惯性可忽略）、双向强耦合（考虑固体惯性、大变形或接触）。按时间推进方式又可分为稳态耦合与瞬态耦合。工程中常见的换热器管板、涡轮叶片、电子散热器，多数属于双向弱耦合；而阀门启闭、密封条压紧等涉及接触与大位移的问题，则需要强耦合或至少在每个时间步内多次迭代。

五、耦合参数与判断依据

选择耦合类型时，可用一个简单判据：若固体最大热变形 δ 与流道特征尺寸 L 之比 δ/L < 1%，且温度变化不引起流动分离或间隙显著变化，单向耦合通常足够；若 δ/L 在 1%～5% 之间，建议双向弱耦合；若 δ/L > 5% 或存在接触、屈曲、材料非线性，则必须强耦合。此外，流固热耦合的时间尺度差异也需注意：流体热响应时间 τ_f ≈ L²/α_f，固体热响应时间 τ_s ≈ L²/α_s，若两者相差两个数量级以上，可采用不同时间步的准稳态策略。