4. 参数扫描与响应面初步分析

参数扫描是连接“单次仿真”与“优化算法”的桥梁。它的核心任务不是求最优解，而是用一批可控的仿真样本，把设计变量与目标量、约束量之间的映射关系显式地记录下来，形成一张可供后续拟合、插值或代理模型训练的数据表。在 COMSOL 中，这一过程通常由“研究 → 参数化扫描”或“优化 → 参数扫描”节点完成，前者适合快速批量计算，后者便于与优化模块共享变量定义。

一、先明确三类量：设计变量、目标量、约束量

在动手扫描之前，必须先把物理模型中的量按角色分类，否则扫描结果会变成一堆无法解释的曲线。以常见的散热器稳态热-流耦合分析为例：

• 设计变量：翅片厚度 t_fin（0.8–2.0 mm）、翅片间距 p_fin（3–8 mm）、入口流速 v_in（0.5–2.5 m/s）。它们是扫描的自变量，必须给出明确的上下限和步长。
• 目标量：热源表面最高温度 T_max、压降 Δp。前者希望最小化，后者希望不超过允许值。
• 约束量：材料许用温度 T_limit、泵功率上限 P_pump_max、几何可制造性下限。约束量在扫描阶段只记录、不惩罚，留到优化阶段处理。

分类完成后，建议在 COMSOL 的“全局定义 → 参数”中集中声明这些量，并给每个量加上单位。单位缺失是后续提取数据时最常见的错误来源之一。

二、参数扫描的三种典型组织方式

COMSOL 的参数化扫描支持“全部组合”“指定组合”和“手动输入”三种模式。选择哪一种，取决于变量数量和计算预算。

对于本节示例，三个变量各取 4 个水平，全部组合共 64 次求解。若单次稳态求解约 8 秒，总时间约 8–10 分钟，属于可接受范围。若变量增加到 5 个、每变量 5 水平，全组合将膨胀到 3125 次，此时应改用拉丁超立方或正交试验设计。

三、扫描设置的关键细节

在“研究 1 → 步骤 1：稳态 → 参数化扫描”中，逐行添加变量。每个变量的“单位”列必须与参数定义一致；步长建议用“步长”而非“点数”，便于核对。求解前勾选“重用解”，可让相邻参数点以上一次解为初值，显著缩短收敛时间，但要注意：当参数变化跨越流动状态转变（如层流到湍流）时，重用解可能收敛到非物理解，此时应关闭该选项并改用零初值。

另一个容易被忽略的设置是“同时求解”与“依次求解”。对于强耦合多物理场，依次求解更稳定；对于弱耦合问题，同时求解更快。扫描阶段建议先用依次求解跑通全流程，再根据耗时决定是否切换。

四、从结果中提取目标量与约束量

扫描完成后，COMSOL 会在“结果 → 派生值”中生成随参数变化的表格。推荐使用“全局计算”或“体最大值/积分”算子，把 T_max、Δp、P_pump 等量定义为派生值，再通过“表格”节点导出为 CSV。导出时务必勾选“包含参数列”，否则数据表会丢失自变量信息。

下面是一段典型的导出数据表结构（示例数据，非实测）：

t_fin(mm), p_fin(mm), v_in(m/s), T_max(K), delta_p(Pa), P_pump(W)
0.8, 3.0, 0.5, 342.1, 18.4, 0.0092
0.8, 3.0, 1.0, 331.7, 62.3, 0.0623
0.8, 3.0, 1.5, 326.4, 128.9, 0.1934
...
2.0, 8.0, 2.5, 318.9, 41.2, 0.1030

拿到这张表后，不要急于做优化。先做三件事：检查是否有求解失败或明显异常的点；绘制单变量趋势图，确认物理上合理；计算目标量之间的相关性，判断是否存在强耦合。

五、响应面初步分析：从数据表到趋势判断

响应面初步分析的目标是回答两个问题：哪些变量对目标量影响最大？目标量与变量之间大致是线性、二次还是单调关系？在 COMSOL 中，可以用“结果 → 二维绘图组 → 表面图”直接绘制 T_max 对 t_fin 和 v_in 的响应面，固定 p_fin 为中间值。更系统的做法是把 CSV 导入 MATLAB 或 Python，用最小二乘拟合二次响应面：

% 示例：二次响应面拟合（MATLAB）
X = [t_fin, p_fin, v_in];
y = T_max;
mdl = fitlm(X, y, 'quadratic');
disp(mdl);

拟合后重点看三项：各变量的 p 值（显著性）、R²（拟合优度）、以及交互项系数。若 t_fin×v_in 的系数显著，说明两个变量对温度存在耦合影响，后续优化不能简单分开处理。

六、扫描阶段的常见问题与处理

参数扫描最容易出现三类问题。第一类是不收敛：通常发生在参数极端组合处，处理方法是缩小该变量的扫描范围，或改用更宽松的求解器容差。第二类是数据错位：导出时参数列与结果列顺序不一致，务必在导出后核对前几行。第三类是物理异常：例如压降随流速增大反而减小，这往往意味着网格在某个参数点失效，需要检查网格质量或改用变形几何。

完成扫描和初步响应面分析后，你应当得到一张结构清晰的数据表、一组趋势图，以及关于变量显著性和耦合关系的初步结论。这张表和这些结论，就是下一阶段构建代理模型和运行优化算法的输入。

这张图概括了本节的核心链路：设计变量经参数化扫描生成样本，样本以数据表形式承载目标量与约束量，再通过趋势图和二次拟合完成响应面初步分析。实际工程中，扫描范围、水平数和求解设置需要根据模型收敛性和计算预算反复调整，但“先扫描、再分析、后优化”的顺序不应颠倒。