4. 物理量量纲与单位一致性检查
在多物理场耦合模型中,边界条件往往由不同来源拼接而成:结构场的热载荷来自热分析结果,流场的入口速度来自实验数据,电磁场的电流密度来自电路参数。这些数据在进入同一个模型之前,可能分别使用 W/m²、mW/mm²、kW/m²,或者 Pa、MPa、bar、atm,或者 m/s、km/h、ft/s。量纲与单位不一致不会让求解器报错,却会让边界值整体偏移几个数量级,最终表现为"温度高得离谱""压力几乎为零""流速完全不对"这类难以定位的现象。本节要解决的核心问题,就是建立一套可执行的单位制统一检查清单,并掌握单位不一致导致数量级错误时的排查方法。
一、量纲与单位的区别:先分清两个层次
量纲(dimension)描述物理量的本质属性,如温度的量纲是 Θ、压力的量纲是 M·L⁻¹·T⁻²、速度的量纲是 L·T⁻¹。单位(unit)则是量纲的具体度量约定,如 K 与 ℃ 同属温度量纲,Pa 与 bar 同属压力量纲。检查必须分两层:第一层确认量纲是否匹配,例如把热流密度(W/m²)误填到温度边界(K)上,这是量纲错误,属于物理意义错误;第二层确认单位换算系数是否正确,例如把 25 ℃ 直接写成 25 K,这是同量纲下的单位错误,属于数值错误。量纲错误通常源于边界类型选错,单位错误则更隐蔽,因为数值看起来"合理"。
需要特别注意的是偏移型单位。热力学温度与摄氏温度之间是 T(K) = t(℃) + 273.15,存在零点偏移;而压力单位之间、速度单位之间多为纯比例换算。偏移型单位不能只用乘系数处理,必须单独在清单中标注,否则 300 K 与 300 ℃ 的混淆会造成约 27 ℃ 与 573 K 的巨大差异。
二、单位制统一检查清单
建议在建模前建立一份覆盖本节常见物理量的检查清单,逐项确认"模型内部基准单位—数据来源单位—换算系数—是否含偏移"。下表给出可直接使用的示例清单,其中基准单位一列按国际单位制(SI)设定,实际项目可整体替换为自定义一致单位制,但必须全文统一。
| 物理量 | 量纲 | 模型基准单位 | 常见来源单位 | 换算关系 | 偏移 |
|---|---|---|---|---|---|
| 温度 | Θ | K | ℃、℉ | K = ℃ + 273.15;K = (℉−32)×5/9 + 273.15 | 有 |
| 压力 | M·L⁻¹·T⁻² | Pa | MPa、bar、atm、psi | 1 MPa = 1e6 Pa;1 bar = 1e5 Pa;1 atm = 101325 Pa | 无 |
| 速度 | L·T⁻¹ | m/s | km/h、ft/s | 1 km/h = 1/3.6 m/s;1 ft/s = 0.3048 m/s | 无 |
| 热流密度 | M·T⁻³ | W/m² | mW/mm²、kW/m² | 1 mW/mm² = 1000 W/m²;1 kW/m² = 1000 W/m² | 无 |
| 电流密度 | I·L⁻² | A/m² | A/mm²、mA/cm² | 1 A/mm² = 1e6 A/m²;1 mA/cm² = 10 A/m² | 无 |
清单的使用要点有三:其一,每个边界值在录入前必须能追溯到来源单位,不能凭经验假设;其二,换算系数写成"1 来源单位 = x 基准单位"的形式,避免乘除方向混淆;其三,对含偏移的量必须写成完整表达式而非单一系数。完成换算后,建议在模型中对边界值做一次量级预估,例如金属结构热流密度通常在 10³~10⁶ W/m² 量级,若换算后得到 10⁹ 或 10⁻³,就应立即回查单位。
三、案例:热流密度单位混用造成的数量级错误
设某电子散热模型以 SI 为基准单位,芯片底面热流密度由热测试报告给出为 0.8 mW/mm²。若建模者误将其按 0.8 W/m² 直接录入,则真实值应为 0.8 × 1000 = 800 W/m²,录入值仅为真实值的千分之一,仿真得到的芯片温度会明显偏低,甚至接近环境温度,与实测的 85 ℃ 严重不符。反过来,若把 800 W/m² 误当作 800 kW/m² 录入,则热流放大一千倍,温度会瞬间超出材料熔点,求解可能发散。两种错误方向相反,但根源相同:来源单位与基准单位之间缺少显式换算步骤。
再看压力边界。某流场入口给定 0.2 MPa,若模型基准为 Pa 而误录为 0.2 Pa,入口几乎无压差驱动,流速趋近于零;若把 0.2 MPa 误按 0.2 bar 处理,则实际只施加了 0.02 MPa,流量偏低约一个数量级。速度边界同理,入口 36 km/h 若直接写成 36 m/s,实际放大了 3.6 倍,雷诺数随之改变,湍流模型的选择都可能失效。
四、排查方法:从现象反推单位错误
当仿真结果出现异常时,可按"现象—量级—单位"的顺序反查。第一步看偏差倍数:若结果与预期相差约 10³、10⁶ 这类整数幂,几乎可以断定是单位前缀错误(m、k、M 的漏乘或误乘);若相差约 3.6 倍,优先怀疑 km/h 与 m/s;若温度相差约 273,优先怀疑 ℃ 与 K 的偏移未处理。第二步回到边界定义处,逐个核对来源单位与换算表达式,而不是只看最终数值。第三步用"手算一个点"验证:取边界上一点,用来源数据手算基准值,与模型输入值比对,二者不一致即定位到具体边界。
最后需要强调,单位一致性检查不是一次性动作,而应贯穿建模全过程:导入数据时检查来源单位,录入边界时执行换算,求解前做量级预估,结果异常时按倍数反查。把这份清单固化为建模前的固定步骤,才能让多物理场边界真正"对得上号",避免因单位问题浪费大量调试时间。