1. 热流固耦合仿真概述
热流固耦合（Thermo-Hydro-Mechanical，简称 THM）研究的是温度场、流体场与固体力学场之间的双向或多向相互作用。它并不是把三个物理场简单叠加，而是强调“一个场的变化会改变另一个场的控制方程参数或边界条件”。例如温度升高会改变流体黏度与固体热膨胀，流体压力会改变固体有效应力与孔隙率，固体变形又会反过来改变流道几何与渗透率。理解这种“互为因果”的闭环关系，是使用 COMSOL Multiphysics 进行多物理场建模的起点。
一、多物理场耦合的基本概念
在 COMSOL 中，多物理场耦合通过“物理场接口 + 多物理场节点”实现。物理场接口负责描述单个场的控制方程，多物理场节点负责把不同场的变量传递到对方的方程中。常见的耦合方式有三类：
- 单向耦合：A 场影响 B 场，但 B 场不反馈 A 场。例如先算稳态温度场，再把温度作为热载荷施加到结构上做热应力分析。计算量小，适合温度变化对变形影响微弱、或变形对流动影响可忽略的场景。
- 双向耦合：两个场互相影响，需要迭代求解。例如热对流中流速影响温度分布，温度又通过浮力改变流速，形成自然对流。
- 全耦合与分离耦合：全耦合把所有场变量放在同一方程组中一次性求解，收敛快但内存需求高；分离耦合按场依次求解并迭代，内存友好但可能需要更多迭代次数。COMSOL 默认对多数多物理场问题采用全耦合，用户可在求解器配置中切换。
耦合的“强度”决定了求解策略。若耦合项在方程中表现为弱源项，分离求解通常足够；若耦合项出现在主对角或强非线性项中（如流固耦合中的动网格与压力耦合），则全耦合更稳健。
二、COMSOL Multiphysics 界面与建模流程
COMSOL 的建模流程可以概括为“几何—材料—物理场—边界—网格—求解—后处理”七步，但实际工作中这些步骤往往需要反复回溯。新建模型时，首先在“模型向导”中选择空间维度（三维、二维轴对称、二维等）与物理场接口。对于 THM 问题，典型组合是：
- 固体力学接口：描述位移、应力、应变，可加入热膨胀与孔隙弹性。
- 传热接口：描述温度场，可包含固体传热、流体传热与多孔介质传热。
- 流体流动接口：层流、湍流或达西流，取决于介质是自由流体还是多孔介质。
- 多物理场节点：如“热膨胀”“温度耦合”“流固耦合”“多孔弹性”等，用于建立场间变量传递。
几何建模时，建议先建立简化二维模型验证耦合逻辑，再扩展到三维。材料属性中，热导率、比热容、密度、弹性模量、泊松比、渗透率、孔隙率、流体黏度等参数都可能随温度或压力变化，COMSOL 允许把它们写成温度或压力的函数表达式。边界条件要特别注意：温度场常用热绝缘、温度、热通量；流场常用入口流速/压力、出口压力、壁面；固体场常用固定约束、指定位移、载荷。网格方面，流固耦合区域需要动网格或 ALE 方法，边界层网格对传热与流动精度影响显著。
三、热-流-固三场耦合的典型工程场景
THM 耦合在工程中非常普遍。地热开采中，低温水注入高温岩体，岩体冷却收缩、渗透率变化，同时流体被加热，形成典型的 THM 闭环；核废料地质处置中，衰变热使缓冲材料膨胀，孔隙水压力升高，可能影响围岩稳定性；锂电池热管理中，电芯发热引起结构膨胀，冷却液流动又带走热量；管道输送高温高压流体时，管壁热应力与流体压力共同作用，可能引发屈曲或疲劳。这些场景的共同点是：温度、压力、变形三者相互制约，忽略任一耦合都可能得到错误结论。
上图给出了三场之间的主要耦合路径。可以看到，温度场通过黏度、密度、热膨胀影响流体与固体；流体场通过对流换热和压力作用影响温度与固体；固体场通过变形改变流道几何、渗透率与传热路径。实际建模时，不必一开始就打开所有耦合项，而应根据工程问题的物理主导机制，先判断哪些耦合是“必须保留”的，哪些可以暂时忽略。例如地热开采中，温度对流体黏度和岩石渗透率的影响通常不可忽略，而固体变形对温度场的反馈在浅层低应力条件下可能较弱，可先做单向或部分双向耦合验证，再逐步增加耦合项。
四、建模前的关键判断
在动手建模之前，建议先回答三个问题：第一，三个场的时间尺度是否接近？如果温度变化以年计、流体流动以天计、固体变形以秒计，可能需要分步或准静态处理。第二，耦合是单向还是双向？可以通过量级估算判断，例如热膨胀应变与弹性应变的比值、浮力与惯性力的比值。第三，几何是否会发生大变形？若流固界面位移超过特征长度的 5%～10%，通常需要动网格或 ALE 方法，否则可忽略几何非线性。这三个判断直接决定物理场接口选择、求解器配置与计算成本，是 THM 仿真从“能算”走向“算得对”的关键一步。