4. 多物理场参数尺度差异的识别与量化
多物理场耦合问题的困难，往往不在于单个场的方程有多复杂，而在于不同物理场在空间尺度、时间尺度、参数量级上可能相差几个数量级。如果直接对原始方程做无量纲化或数值离散，很容易出现“一个场已经稳定、另一个场还在剧烈变化”的刚性现象。因此，本节的核心任务是：从控制方程和边界条件中提取特征尺度，计算各物理场参数的典型量级，并把这种差异用一张“参数尺度谱”直观地表达出来。
一、特征尺度从哪里来
特征尺度不是凭空假设的，它来自三类信息：几何尺寸、边界条件和材料物性。以微流控芯片散热为例，芯片流道宽度约为 W = 200 μm，基板厚度 H = 500 μm，流道长度 L = 2 cm；入口流速 U = 0.05 m/s；入口温度 T_in = 300 K，热源热流密度 q'' = 50 kW/m²；流体为水，ρ = 998 kg/m³，c_p = 4182 J/(kg·K)，k_f = 0.6 W/(m·K)，μ = 1.0×10⁻³ Pa·s；基板为硅，k_s = 150 W/(m·K)。
由此可以定义一组基础特征量：长度尺度取流道水力直径 D_h ≈ 200 μm（矩形截面近似），时间尺度取对流时间 t_c = L/U = 0.4 s，速度尺度取 U，温度尺度取温升参考值 ΔT = q''·H/k_s ≈ 0.17 K（仅作量级参考）。这些量构成了后续所有无量纲参数的“标尺”。
二、从方程中提取无量纲参数
对连续性方程和动量方程做无量纲化，可以得到雷诺数：
Re = ρ·U·D_h / μ
= 998 × 0.05 × 2.0e-4 / 1.0e-3
≈ 10Re ≈ 10 说明流动处于层流区，惯性项与黏性项同量级但黏性占优。对能量方程无量纲化，得到佩克莱数：
Pe = U·D_h / α_f, α_f = k_f/(ρ·c_p) ≈ 1.44e-7 m²/s
Pe ≈ 0.05 × 2.0e-4 / 1.44e-7 ≈ 69Pe ≈ 69 表明沿流动方向的对流输运远强于横向导热，温度场在流道横截面上接近均匀，而在流动方向形成明显温升梯度。再看固体与流体之间的共轭传热，用导热系数比和热扩散率比刻画：k_s/k_f = 250，α_s/α_f ≈ 6.4×10³（硅的热扩散率约 9.2×10⁻⁵ m²/s）。这两个比值说明：热量在硅基板中扩散极快，而在水中主要靠流动带走，固液界面的热流分配严重不对称。
三、时间尺度的分离与刚性
更值得警惕的是时间尺度的差异。流体对流时间 t_c = 0.4 s；而硅基板的热扩散时间 t_s = H²/α_s ≈ (5×10⁻⁴)²/9.2×10⁻⁵ ≈ 2.7×10⁻³ s；流道横向热扩散时间 t_f = D_h²/α_f ≈ (2×10⁻⁴)²/1.44×10⁻⁷ ≈ 0.28 s。可见固体热响应比流体对流快约两个数量级，若用显式时间推进，时间步长会被固体导热限制在毫秒量级，而总模拟时长却由对流决定，计算量被显著放大。这正是多物理场“刚性”的典型来源。
四、参数尺度谱的读法与用途
上图把四类尺度放在同一张“谱”上，横向位置表示量级，颜色区分物理场。读图时先看跨度：长度跨 2 个数量级、时间跨约 2.5 个数量级，说明几何与时间都存在明显多尺度特征；再看错位：固体热扩散时间落在谱的最左端，而对流时间在最右端，二者错位越大，数值刚性越强；最后看无量纲数：Re 与 Pe 决定了流动与传热的主导机制，k_s/k_f 与 α_s/α_f 决定了共轭界面的耦合强度。
这张谱的实用价值在于指导建模决策。例如，由于 Pe 较大且 α_s/α_f 极高，可以合理采用“流体沿流向一维简化 + 固体三维导热”的降阶策略；由于时间尺度错位明显，应优先选择隐式时间格式或稳态求解，而不是显式推进。反过来，如果某个参数落在谱的中间区域、与其它场同量级，就不能随意简化，必须保留完整耦合。
把上述过程固化为一张可填写的参数尺度表，便于在其它多物理场问题中复用：列出每个场的特征长度、特征时间、特征速度/温度、派生无量纲数及其量级，并标注“是否与其它场同量级”。同量级者保留完整耦合，跨 2 个数量级以上者考虑降阶或隐式处理。这样，尺度差异就从模糊的“感觉很难算”变成了可量化、可决策的工程判断。