2. 对标问题的数学物理描述

实验对标的第一步不是打开仿真软件，而是把实验台上真实发生的事，翻译成一组数学上可解、物理上可解释的方程。这个翻译过程就是数学物理描述：确定控制方程、边界条件、初始条件，并从实验现象中提取关键物理机制，据此给出合理的简化假设。翻译得好，仿真与实验的偏差主要来自模型误差；翻译得差，后面再精细的网格和再多的求解器设置都只是在解一个错误的问题。

一、控制方程：守恒律在具体物理场中的形式

控制方程本质上是守恒律在特定物理场中的具体写法。做对标时，先问三个问题：守恒量是什么（质量、动量、能量、电荷、组分浓度）？输运机制有哪些（对流、扩散、传导、辐射）？是否存在源项（化学反应、相变、焦耳热、体积力）？

以常见的散热器稳态对流换热为例，若流体不可压缩、物性近似为常数，控制方程可写为：

连续性：  ∇·u = 0
动量：    ρ(u·∇)u = -∇p + μ∇²u + ρg
能量(流体)：ρcp(u·∇T) = kf∇²T
能量(固体)：∇·(ks∇T) = 0

这里每个符号都对应实验中的一个可测量或可估算量：ρ、μ、cp、kf、ks 来自物性表或实测，u 来自流量计换算，T 来自热电偶。若实验中发现温差较大导致物性变化超过 10%，就应把 ρ(T)、μ(T) 写成温度函数，而不是继续用常物性——这正是从实验现象反推模型复杂度的典型判断。

二、边界条件：把实验台的约束翻译成数学约束

边界条件是实验与仿真之间最容易出错、也最值得逐条核对的部分。常见类型与实验对应关系如下：

实操中建议做一次“边界条件审计”：对每一条边界，写出实验中的物理载体、测量手段和不确定度。例如散热器底面接加热铜块，实测铜块温度 T_cu = 65 ℃、加热功率 Q = 20 W、接触面积 A = 1.6×10⁻³ m²，则名义热流 q = Q/A = 12.5 kW/m²。若接触热阻 R_c = 1.5×10⁻⁴ m²·K/W，界面温升 ΔT = q·R_c ≈ 1.9 K，说明底面温度应取 65 - 1.9 ≈ 63.1 ℃ 而非直接取 65 ℃。这类修正往往比加密网格更能改善对标一致性。

三、初始条件与瞬态对标

稳态问题不需要初始条件，但瞬态对标必须交代初始场。初始条件不是“随便给一个”，而应来自实验的起始状态：例如加热开始前试样已与环境平衡 30 min，则初始温度场可设为均匀 T₀ = 25 ℃；若试样刚从冰箱取出，则需给出非均匀初始分布或至少给出中心与表面的实测初值。初始条件设置不当，会导致早期温升曲线整体平移，容易被误判为材料参数错误。

四、从实验现象提取关键物理机制

实验数据不会直接告诉你“该用哪个方程”。需要从现象中识别主导机制。常用判断手段有三类：

• 量级比较：用无量纲数判断。例如自然对流与强制对流之比 Gr/Re²，若远小于 1 则强制对流主导，可忽略浮升力；若接近或大于 1，则必须考虑混合对流。
• 时间尺度比较：热扩散时间 τ_d = L²/α 与实验采样/加热时间比较。若 τ_d 远小于观测时间，可视为稳态；若相当，则必须做瞬态。
• 空间分布观察：若实测表面温差小于平均温升的 5%，可尝试用集中参数或一维简化；若局部热点明显，则必须保留三维几何与局部边界。

例如某散热器实验中，风速 2 m/s、特征长度 0.05 m，Re ≈ 6.7×10³，Gr ≈ 1.2×10⁵，Gr/Re² ≈ 2.7×10⁻³，说明强制对流占绝对主导，浮升力可忽略。这一判断直接决定了是否需要在动量方程中保留 ρg 项，也决定了网格是否需要延伸到远场。

五、简化假设的确定与记录

简化假设不是越少越好，而是“每条假设都有依据、都能被检验”。建议在建模前写一份假设清单，至少包含：假设内容、依据、可能影响、验证方式。例如：