2. 流体力学与结构力学基础回顾

流固耦合仿真之所以比单一物理场仿真更难，根本原因在于它同时求解两套控制方程，并在交界面上反复交换信息。如果对这两套方程的形式、变量含义和耦合条件没有清晰认识，后续在 COMSOL 中设置接口、选择耦合类型、判断收敛困难时就会缺少依据。本节把流体侧、固体侧和交界面三部分的核心关系梳理清楚，并说明几个关键无量纲参数在工程判断中的作用。

一、流体侧：Navier-Stokes 方程与它的物理含义

不可压缩牛顿流体的运动由连续性方程和动量方程共同描述。连续性方程表达质量守恒：

∇·u = 0

动量方程表达牛顿第二定律在流体微元上的体现：

ρ(∂u/∂t + (u·∇)u) = -∇p + μ∇²u + f

其中 u 为速度矢量，p 为压力，ρ 为密度，μ 为动力黏度，f 为体积力。方程左侧第一项是当地加速度，第二项 (u·∇)u 是对流加速度，它使方程成为非线性；右侧依次为压力梯度力、黏性扩散力和体积力。理解这一点很关键：流固耦合中流体方程的非线性是迭代收敛困难的主要来源之一，而黏性项决定边界层和壁面剪切的解析精度。

在 COMSOL 中，层流接口直接求解上述方程；当雷诺数较高、惯性效应显著时，需要改用湍流模型（如 k-ε、k-ω SST），此时不再直接解析全部尺度，而是引入湍流黏度和附加输运方程。选择哪种模型，取决于流动特征而非个人偏好。

二、固体侧：固体力学基本方程

固体域通常在小变形假设下求解，其控制方程为运动方程：

ρ_s ∂²d/∂t² = ∇·σ + F_v

其中 d 为位移矢量，σ 为柯西应力张量，F_v 为体积力。配合几何方程 ε = ½(∇d + ∇dᵀ) 和本构方程 σ = C:ε（线弹性时 C 为弹性张量），构成完整的固体力学定解问题。对于大变形或柔性结构（如薄膜、软体、薄板大挠度），必须开启几何非线性，否则位移与应变关系失真，耦合结果会明显偏离实际。

结构侧的关键参数是弹性模量 E、泊松比 ν、密度 ρ_s 和阻尼。它们不仅决定结构自身响应，还通过耦合界面影响流体载荷的传递方式。例如 E 越小，结构越容易被流体压力推动，双向耦合的必要性就越高。

三、耦合界面条件：两套方程如何“对话”

流固交界面上必须同时满足两个条件。其一是运动学条件（位移/速度连续）：

u_fluid = ∂d_solid/∂t （界面法向速度一致）

其二是动力学条件（应力连续）：

σ_solid·n = -p·n + τ·n （界面法向与切向应力平衡）

前者保证流体不穿透固体、固体不脱离流体；后者保证界面两侧受力平衡。COMSOL 中“流固耦合”多物理场节点正是自动施加这两个条件，用户需要做的是正确指定耦合边界，并保证两侧网格在界面处协调。若使用分离式求解，界面上的数据传递顺序和松弛因子会直接影响稳定性。

四、无量纲参数：判断耦合行为与建模策略

在动手建模前，先用几个无量纲数估算物理特征，能避免选错模型。雷诺数 Re = ρUL/μ 表征惯性力与黏性力之比：Re 小则流动层流、稳定；Re 大则惯性主导，可能出现湍流、涡脱落，需要湍流模型和更细的网格。对于流固耦合，还要关注以下参数。

以水中柔性薄板为例：若板长 L=0.1 m，来流 U=0.5 m/s，水密度 ρ=1000 kg/m³，黏度 μ=1×10⁻³ Pa·s，则 Re≈5×10⁴，属于湍流范围；若板材料 E=2 MPa、ρ_s=1100 kg/m³，则质量比 m*≈1.1，属于低质量比问题。此时若采用分离式求解且松弛不足，界面迭代很容易发散，通常需要改用耦合式求解器或加入 Aitken 松弛。这些判断都来自无量纲参数的估算，而不是试错。

五、从方程到 COMSOL 设置的对应关系

理解方程之后，COMSOL 中的设置就有了明确指向：流体域选择“层流”或“湍流”接口，对应是否保留全部 N-S 项；固体域选择“固体力学”接口，并决定是否勾选几何非线性；交界面通过“流固耦合”多物理场节点自动施加两个界面条件；瞬态研究对应 ∂u/∂t 和 ∂²d/∂t² 项，稳态研究则忽略时间导数。求解器方面，双向耦合通常采用全耦合或分离式迭代，收敛判据应同时监控流体残差、固体残差和界面位移变化量。

本节建立的方程框架和参数判断，是后续所有流固耦合建模的基础。无论是阀门、管道、翼型还是生物软组织，其 COMSOL 实现都归结为：正确写出流体侧和固体侧的控制方程，在交界面上施加运动学与动力学条件，并用无量纲参数指导模型选择和求解策略。