1. 多物理场耦合问题与量纲分析的角色
在真实工程系统中,热、流、固三类物理过程几乎从不孤立存在。以航空发动机涡轮叶片为例:高温燃气通过对流与辐射把热量传给叶片外表面,热量在叶片内部以导热方式向冷却通道传递,冷却空气在内部通道中流动并带走热量,而叶片本身在离心力、气动力与热应力共同作用下发生变形,变形又反过来改变冷却通道的几何形状与外部流场。这就是典型的热-流-固三场耦合问题。类似场景还出现在电子芯片散热、电池热失控、核反应堆安全分析、高超声速飞行器热防护等大量工程对象中。
多物理场耦合带来的第一个直接挑战,是参数尺度差异。仍以涡轮叶片为例,其几何尺度在 10⁻² m 量级,冷却通道水力直径约 10⁻³ m,边界层厚度约 10⁻⁵ m,而湍流耗散尺度可低至 10⁻⁶ m;时间上,燃气脉动周期约 10⁻⁴ s,叶片热响应时间常数约 10² s,蠕变变形则发生在 10⁴ h 量级。空间跨 4~5 个数量级、时间跨 9~10 个数量级,若不加处理直接离散,网格量与计算步数将不可承受。第二个挑战是量纲不一致:温度以 K 计、速度以 m/s 计、应力以 Pa 计、热流密度以 W/m² 计,不同物理量在耦合界面上的传递必须通过无量纲化才能比较与组合。第三个挑战是耦合强度差异:有些耦合是强耦合(如流固界面位移与压力互为边界条件),有些是弱耦合(如温度对材料参数的缓慢影响),若不加区分地全耦合求解,计算代价与收敛风险都会急剧上升。
量纲分析正是应对上述挑战的基础工具,其作用体现在三个层面。第一,模型简化:通过 Buckingham π 定理把多个有量纲参数组合成少数无量纲数,例如把流体密度 ρ、速度 U、特征长度 L、黏度 μ 组合为雷诺数 Re = ρUL/μ,把对流换热系数 h、导热系数 k、特征长度 L 组合为努塞尔数 Nu = hL/k。当某些无量纲数处于特定范围时,可判定某项物理效应可忽略,从而把三场耦合降为两场甚至单场问题。第二,实验设计:缩比实验必须保证关键无量纲数相等(相似准则),否则实验数据无法外推到原型。第三,数值模拟:无量纲化后的控制方程系数归一,可改善矩阵条件数、统一不同物理场的收敛判据,并让残差量级具有可比性。
上图把本节的核心逻辑串成一条链:先识别三个物理场及其耦合路径,再看到尺度与量纲差异,最后落到量纲分析的三大作用与课程任务链。需要强调的是,耦合路径的方向性必须明确——热通过改变物性影响流动,流动通过对流换热影响温度,流动通过压力载荷影响结构,结构通过几何变形影响流场,热通过热应力影响结构,结构通过接触导热影响温度。每一条路径都对应一组需要无量纲化的参数。
下面用一个可操作的示例说明量纲分析如何落地。假设我们要分析某冷却通道的对流换热强度,已知通道水力直径 D = 2 mm,冷却空气流速 U = 60 m/s,密度 ρ = 4.0 kg/m³,动力黏度 μ = 3.0×10⁻⁵ Pa·s,导热系数 k = 0.05 W/(m·K),比热容 cp = 1050 J/(kg·K)。先算雷诺数:Re = ρUD/μ = 4.0×60×0.002/3.0×10⁻⁵ ≈ 1.6×10⁴,处于湍流区。再算普朗特数:Pr = μcp/k = 3.0×10⁻⁵×1050/0.05 ≈ 0.63。选用 Dittus-Boelter 关联式 Nu = 0.023 Re0.8 Pr0.4,得 Nu ≈ 0.023×(1.6×10⁴)0.8×0.630.4 ≈ 43。由此 h = Nu·k/D ≈ 43×0.05/0.002 ≈ 1075 W/(m²·K)。这一组数值属于教学示例假设,用于演示流程,不代表任何具体型号发动机的实测结果。
得到 h 之后,就可以计算毕渥数 Bi = hL/ks(L 为固体特征尺寸,ks 为固体导热系数)。若叶片壁厚 L = 3 mm、镍基合金 ks ≈ 11 W/(m·K),则 Bi ≈ 1075×0.003/11 ≈ 0.29。Bi 远小于 1 时,固体内部温度梯度相对较小,可采用集总参数法把叶片壁面温度近似为均匀,从而把“热-固”强耦合降为“热场给固场施加均匀温度载荷”的单向耦合。这就是量纲分析直接指导模型简化的典型过程:不是凭经验拍脑袋,而是用无量纲数的大小给出可判定的依据。
在数值模拟层面,无量纲化的价值同样具体。把动量方程用特征速度 U、特征长度 L 归一后,黏性项系数变为 1/Re,对流项系数为 1;把能量方程归一后,对流项与扩散项之比由 Pe = Re·Pr 控制。这样不同物理场的方程系数都落在可比较的量级上,残差收敛判据可以统一设为 10⁻⁶ 量级,而不必为温度场设 10⁻³、为流场设 10⁻⁵ 这样彼此割裂的标准。同时,当 Re 很大时,可以判断某些低雷诺数修正项可忽略,从而减少耦合迭代次数。
本课程后续将沿着“识别物理场—列出参数—构造无量纲数—判断耦合强弱—指导简化与仿真”这条主线展开,覆盖热-流-固、电-热、流-化等多类耦合场景,并配合可复算的示例参数与判断依据。本节的任务是建立整体框架:多物理场耦合的困难来自尺度差异、量纲不一致与耦合强度差异,而量纲分析是把这些困难转化为可计算、可判断、可简化问题的核心工具。掌握这一框架,后续每一类耦合问题的处理都有了统一的出发点。