第一次打开 COMSOL，左侧模型树一长串节点，几何、材料、物理场、网格、研究，每一步都点了，结果一运行就报错“未定义材料属性”或者“找不到解”。这几乎是每个新手都会经历的阶段。COMSOL 仿真入门真正卡人的地方，往往不是有限元理论，而是边界设置、网格划分和数据处理这三件事上的具体取舍：边界条件给多了算不动，给少了结果没意义；网格粗了看不出梯度，细了内存直接爆；好不容易跑出云图，却不知道怎么把一条曲线、一个积分值导出来写进报告。
这篇内容面向刚接触 COMSOL、已经能照着案例点完流程但独立建模就发怵的读者，也适合做了一段时间但结果总被质疑的同行。我会把边界设置、网格划分、数据处理这三块拆开，讲清楚每个选择背后的理由，再补上几个实际会遇到的卡壳现场。所有操作基于常见版本的功能布局，不同版本菜单文字可能略有差异，但逻辑一致。
1. 建模之前先把物理画面想明白，再动手画几何
1.1 从现象到方程：物理场不是越多越好
很多人上手就想着“多物理场耦合”，把固体力学、传热、流体全勾上，结果耦合项没配好，求解器直接不收敛。我的建议是先把要解释的现象用一句话说清楚：是结构在温度载荷下的变形，还是流体掠过表面带走热量，还是电磁场在介质里的分布。这句话决定了主物理场是谁，其余都是配角。
主物理场选错，后面边界条件和网格全都白搭。比如做电子器件散热，主物理场应该是传热，流体只是提供对流换热系数；如果你把流体也建成完整的三维湍流模型，计算量会翻好几倍，而真正影响芯片温度的可能只是那个换热系数。我见过某实验室的 A 同学，为了算一块散热板的温度分布，建了完整的风扇流场，网格量飙到几百万，最后结论和用经验换热系数算出来的差不到 3%。这不是说详细模型没价值，而是要先问自己：这次分析要回答的问题，需要多高的保真度。
判断依据可以按这个顺序走：先看是否有现成的经验关联式能替代某个物理场，再看该物理场是否只影响边界而不影响内部，最后才考虑全耦合。COMSOL 的多物理场节点里有很多预置耦合，比如热膨胀、共轭传热、焦耳热，选对了能省掉手动映射的麻烦，但前提是主物理场本身设置正确。
1.2 几何简化与维度选择：2D、2D 轴对称还是 3D
维度选择是新手最容易“过度建模”的地方。一个圆柱形加热棒，如果载荷和材料都沿周向均匀，用二维轴对称就能算，计算量可能只有三维的二十分之一。判断标准很简单：几何、材料、边界条件三者中，只要有一项沿周向变化，就不能用轴对称。
二维轴对称要求几何画在 r-z 平面内，且 r≥0，旋转轴是 r=0 那条边。设置时要在物理场里确认“轴对称”假设已经生效，否则软件会按普通二维平面处理，结果完全不对。我见过有人把轴对称模型算出来的“厚度方向”结果直接当成直径方向，单位没换算，结论差了十几倍。
三维建模时，对称面要用起来。一个几何和载荷都对称的支架，切一半加对称边界，单元数直接减半。但要注意，对称边界约束的是法向位移或法向通量，不是把所有方向都锁死。固体力学里的“对称”会约束法向位移，传热里的“对称”意味着法向热流为零，这两者在物理上是一致的，但误用成“固定约束”就会把自由收缩的边锁死，热应力虚高。
几何清理也不能省。COMSOL 的布尔操作和“虚拟操作”可以去掉小圆角、小孔、窄面，这些特征对网格质量杀伤力很大。一个 0.1 mm 的倒角放在 100 mm 的零件上，如果不处理，网格要么在倒角处畸形，要么整体加密到算不动。虚拟操作里的“忽略边”“形成复合域”能在不改变几何外观的前提下简化拓扑，值得花十分钟处理。
1.3 全局参数与单位：后面所有麻烦的源头
参数化建模不是高级功能，是入门就该养成的习惯。几何尺寸、材料属性、边界载荷都写成参数，放在“全局定义-参数”里，单位显式标注。COMSOL 支持带单位的表达式，比如L = 50[mm]、k = 45[W/(m*K)]，这样后续改尺寸不用回几何里一个个改。
单位混乱是报错的高发区。热通量单位是 W/m²，体热源是 W/m³，线热源是 W/m，三者差一个维度，填错量级结果能差几个数量级。我的做法是建一个参数表，把所有载荷按“物理量-符号-数值-单位-来源”列清楚，写进模型注释里。后面别人接手或者自己隔两周再看，能少很多困惑。
| 参数名 | 含义 | 示例值 | 单位 | 备注 |
|---|---|---|---|---|
| L | 特征长度 | 50 | mm | 几何主尺寸 |
| k | 导热系数 | 45 | W/(m*K) | 材料手册取值 |
| q0 | 面热源 | 1e4 | W/m² | 边界载荷 |
| h | 对流换热系数 | 8 | W/(m²*K) | 自然对流经验值 |
| Tinf | 环境温度 | 298.15 | K | 与初始温度区分 |
还有一点，初始条件和边界条件别打架。瞬态传热里初始温度设 293.15 K，边界又给一个 373.15 K 的固定温度，这没问题；但如果初始温度给 400 K，边界给 300 K，第一步就会出现剧烈梯度，时间步被迫压得很小，求解时间成倍增加。初始条件尽量贴近物理上的初始状态，这不是形式主义。
2. 边界条件写错，后面全白算：几类高频设置的取舍逻辑
2.1 边界条件到底在约束什么
边界条件不是“填个数”那么简单，它是在给偏微分方程定解。同一个物理场，不同的边界条件组合决定了方程有没有唯一解、解是否物理。数学上分三类：第一类给变量值（Dirichlet），第二类给通量（Neumann），第三类给变量与通量的线性关系（Robin）。对应到 COMSOL 里，固定温度是第一类，热通量是第二类，对流热通量是第三类。
理解这个分类的好处是能判断“过约束”。一个传热问题，如果所有外边界都给固定温度，内部热源产生的热量只能靠这些边界流出，物理上说得通；但如果同时在一个面上给固定温度又给热通量，软件会以其中一个为准或者直接冲突，结果不可信。流体问题更敏感，入口给速度、出口给压力是常规组合，入口给压力、出口也给压力，速度场就没有约束，容易发散。
我常用的自检方法是：把每个边界条件逐条读一遍，问自己“这条边在真实物理里发生了什么”。如果答不上来，说明这个条件是被迫加的，不是想清楚加的。
2.2 固体力学与传热里的常见组合
固体力学里，“固定约束”锁死全部位移，适合模拟焊接底座、螺栓连接面；“辊支承”只锁法向，适合模拟放在光滑平面上的零件；“边界载荷”给单位面积力，方向要注意是跟随法向还是全局坐标。热应力分析里，最容易犯的错是把所有外表面都固定，然后抱怨热应力大得离谱。真实结构总有一个自由收缩的方向，全固定等于给它加了无限刚度。
传热里，“热通量”的正方向是沿边界外法向流出，填正值表示散热，填负值表示加热。很多人在这里符号搞反，算出来的温度分布和预期完全镜像。“对流热通量”需要两个参数：换热系数 h 和环境温度 T_ext，h 的取值直接决定结果，自然对流空气通常 5~10 W/(m²*K)，强制风冷可能 20~100，水冷上千，取错一个量级结果就废了。
表面辐射在高温场景不能忽略。COMSOL 的“表面对环境辐射”需要发射率 ε 和环境温度，ε 对抛光金属可能只有 0.05，对黑体接近 1。我见过有人算 800 ℃ 的炉内件，只给对流不给辐射，结果温度低了上百度。辐射换热和温度四次方成正比，高温下占比会迅速上升，这个判断不能拍脑袋。
| 物理场 | 边界条件 | 物理含义 | 常见误用 |
|---|---|---|---|
| 固体力学 | 固定约束 | 约束全部位移分量 | 把应自由收缩的面全锁死，热应力虚高 |
| 固体力学 | 辊支承 | 只约束法向位移 | 用在斜面上忽略它按面法向约束 |
| 传热 | 热通量 | 单位面积热流，向外为正 | 符号填反，热流方向与预期相反 |
| 传热 | 对流热通量 | 换热系数与环境温度 | 自然对流场景取了强制对流的 h |
| 流体 | 入口/出口 | 速度或压力边界 | 两端都给压力，速度场无约束 |
| 电磁 | 散射边界 | 吸收外向波 | 离结构太近，反射污染近场 |
2.3 流体与电磁场的入口出口陷阱
流体入口给速度还是给压力，取决于你知不知道流量。知道流量给速度（或者质量流），知道压差给压力。入口给速度时，出口通常给压力，比如 0 Pa 表压；如果出口也给速度，整个域的质量守恒就被过度约束，求解器会用很奇怪的方式“消化”这个矛盾，结果看着收敛了但不可信。
壁面条件要区分无滑移和滑移。绝大多数固体壁面用无滑移，靠近壁面的速度梯度靠边界层网格捕捉；如果用了滑移，等于假设壁面没有摩擦，只适合对称面或者某些简化模型。我见过把管道壁设成滑移的模型，压降算出来只有实际的三分之一，原因就在这里。
电磁场里，“理想电导体”适合模拟金属表面，它强制切向电场为零；“散射边界条件”吸收外向波，但要求波垂直入射，斜入射时反射会变大，所以放置位置要离结构至少半个波长；“完美匹配层”吸收效果更好，但需要建一层额外的域并划分扫掠网格。端口条件要分清是波导端口还是集总端口，前者给模式场，后者给阻抗或功率，用错类型激励效率算不对。
2.4 检查边界条件是否“自洽”的三个动作
第一个动作是看模型树里有没有黄色警告。COMSOL 会在物理场节点上标出未定义的边界，比如“未指定材料的域”，这些警告不处理，求解一定出问题。第二个动作是画一遍边界选择。在“选择”列表里逐个点开边界条件，看高亮的边是不是你心里想的那条，尤其是导入几何后边界编号会变，凭记忆填编号十有八九错。
第三个动作是做量纲和量级估算。给了一个 1e6 W/m² 的热通量，先手算一下温升大概多少：qL/k，如果 L 是 10 mm、k 是 50，温升就是 200 K，这个量级合理不合理，心里要有数。算出来温升几万度，那肯定是单位或者数量级错了。这个习惯能拦住大部分低级错误。
3. 网格划分不是越细越好：把计算资源花在梯度大的地方
3.1 物理场控制与用户控制的适用边界
COMSOL 默认给的是“物理场控制网格”，九档精度从极粗到极细，软件根据物理场自动决定哪里加密。这个模式适合快速跑通流程、做趋势判断，但不适合最终出结果，因为你不知道它到底怎么分的。我一般先用“常规”档跑一遍，看收敛性和结果量级，再切到“用户控制网格”手动调。
用户控制网格的核心是“大小”节点，里面有最大单元尺寸、最小单元尺寸、最大单元增长率、曲率因子、狭窄区域分辨率。这五个参数不是随便填的：最大单元尺寸决定整体单元数，最小单元尺寸防止薄特征处单元过小导致时间步爆炸，增长率控制从细到粗的过渡（一般 1.2~1.5），曲率因子控制曲面处的单元密度（越小越密），狭窄区域分辨率控制薄缝里的层数。
一个实用起点：最大单元尺寸取几何特征尺寸的 1/10 到 1/20，最小单元尺寸取最大尺寸的 1/100 到 1/1000。这个范围能覆盖大多数场景，后面再根据梯度调整。填完之后看“统计”里的单元数和最小单元质量，质量低于 0.1 的单元要重点关注，它们往往是收敛困难的根源。
3.2 单元类型与几何拓扑的匹配
自由四面体是最通用的，适合任意复杂几何，但同样精度下单元数比六面体多。扫掠网格适合拉伸体、圆柱、薄板这类有一个方向拓扑一致的结构，生成六面体或棱柱，单元数少、精度高。映射网格适合规则的面或体，四边形或六面体分布均匀。边界层网格专门处理壁面法向的梯度，流体和传热都用得上。
选择逻辑可以按几何来：能扫掠就扫掠，不能扫掠再用自由四面体；薄壁结构优先边界层或映射；二维轴对称用自由三角形就够，规则区域可以映射成四边形。我见过有人对着一根圆管用自由四面体硬划，单元数比扫掠多五倍，结果还没扫掠准。
| 几何特征 | 推荐单元 | 理由 |
|---|---|---|
| 拉伸体、圆柱 | 扫掠+六面体 | 单元少、沿扫掠方向精度高 |
| 薄壁、涂层 | 边界层或映射 | 厚度方向需要分层捕捉梯度 |
| 复杂曲面铸件 | 自由四面体 | 适应任意拓扑，局部加密灵活 |
| 二维轴对称 | 自由三角形或映射四边形 | 计算量小，规则区可用结构化网格 |
| 流体近壁区 | 边界层+自由单元 | 壁面法向速度梯度需要多层单元 |
3.3 边界层网格与局部加密的实操参数
流体和传热问题里，壁面附近的梯度最陡，边界层网格的第一层厚度直接决定能不能捕捉到。层流里可以用经验公式估算边界层厚度，湍流里要看 y+ 目标。COMSOL 的边界层节点可以设层数、第一层厚度、厚度增长率，一般层数 5~10 层，增长率 1.2 左右。
第一层厚度怎么定？一个粗略办法是用特征长度除以雷诺数的平方根再乘个系数，但更稳妥的做法是先跑一遍，看壁面处的速度或温度梯度，再回来调。如果第一层太厚，壁面温度会算得偏高或偏低；如果太薄，单元长宽比过大，质量下降。我一般先给一个保守值，跑通之后做一次网格无关性检查，确认第一层厚度减半后结果变化小于 2%。
局部加密用“大小”节点配合选择，只对关键区域加密，比如热源附近、孔口、接触面。全局加密的代价是单元数按维度三次方涨，局部加密能控制在关键区域。加密之后记得检查过渡区的增长率，突变会导致单元畸形。
3.4 网格无关性验证怎么做才不算浪费时间
网格无关性不是把网格翻倍跑一遍看看，而是要有明确的监测量。先选一个对你结论最敏感的标量，比如最高温度、最大应力、流量、换热系数，然后做三套网格：粗、中、细，单元数大概按 1:2:4 递增。看这个标量随网格的变化，如果中网格到细网格的变化小于 2%，可以说基本无关；如果还在明显变化，继续加密或者检查是不是有奇点。
奇点是个坑。尖角、点载荷、理想接触边缘处的应力理论上无限大，网格越细应力越高，永远不收敛。这种位置的“网格无关性”没有意义，应该看离开奇点一定距离的路径上的应力，或者用圆角代替尖角。我见过有人为了一个尖角的应力值反复加密，算了三天，结论还是不能用。
还有一个省时间的做法：用二维模型或者轴对称模型做网格无关性，再把结论迁移到三维。梯度的相对分布是相似的，二维上花十分钟能定下来的网格密度，三维上直接照搬比例，能省不少机时。
4. 求解不收敛时，先别改网格：按顺序排查更省时间
4.1 残差曲线在说什么
求解不收敛，很多人的第一反应是“网格不够细”，于是加密重算，结果还是不收敛，时间全浪费了。正确的顺序是先看残差曲线。稳态求解的残差应该单调下降，如果震荡、平台、甚至上升，说明问题在方程本身或者边界条件，不在网格。
残差震荡常见于强非线性、接触、相变、湍流。这时候先看是哪个物理场的残差下不去，再到那个物理场里找原因。传热的残差下不去，可能是材料属性随温度变化太剧烈，或者边界条件里有突变；流体的残差下不去，可能是入口湍流参数给得不合理，或者网格长宽比太大。
COMSOL 的求解器日志里会显示每步的残差和阻尼。阻尼因子被压得很小，说明非线性太强，求解器在“小步试探”。这时候可以手动给初值，或者先算一个简化版本（比如去掉某个耦合）作为初值，再逐步加回来。
4.2 稳态与瞬态的初始化差异
稳态问题没有时间概念，求解器直接找平衡点。如果物理上存在多个稳态，或者从零场出发离解很远，就容易不收敛。解决办法是给一个好的初值：先用一个容易收敛的边界条件算一遍，把结果作为初值；或者用参数化扫描逐步加载载荷，从 10% 加到 100%。
瞬态问题相反，初始条件必须自洽。初始温度场和边界温度差太大，第一步就会产生剧烈梯度，时间步被压到极小。我一般先把初始场设成和边界条件接近的值，让系统“平滑启动”，再观察演化。如果关心的是启动瞬间的冲击，那就得接受小时间步，但要检查是不是有物理上不合理的突变。
瞬态的时间步长也不是越小越好。COMSOL 的“自由”时间步会根据局部误差自动调整，手动给最大步长限制是为了输出分辨率，不是为了收敛。步长太小会让求解时间成倍增加，而精度提升有限。
4.3 非线性、接触、相变场景的阻尼与步长
接触问题是最容易不收敛的场景之一。接触对的法向刚度、初始间隙、摩擦系数都会影响收敛。我的经验是先开“稳定接触”或者给一个小的接触弹簧，让接触在初始阶段就建立起来，再逐步调硬。摩擦系数从 0 开始加，比一上来给 0.3 容易收敛得多。
相变问题（比如融化、凝固）在相变区间材料属性剧烈变化，可以用表面对环境辐射里的分段函数，或者用 COMSOL 内置的相变材料节点，把潜热分摊到相变温度区间。区间给得太窄，属性变化太陡，求解器受不了；给得太宽，相变温度就不准。一般按实际相变温度范围的 5~10 K 来设，具体看材料。
非线性材料同理，塑性、超弹性、蠕变都需要逐步加载。参数化扫描配合“续解”功能，把上一步的解作为下一步的初值，比每次都从零开始稳得多。这个技巧在调试阶段特别有用，能快速定位是哪一段载荷区间出的问题。
// COMSOL Java API 示例：逐步加载边界载荷 for (int i = 1; i <= 5; i++) { double factor = i * 0.2; model.param().set("loadFactor", String.valueOf(factor)); model.study("std1").run(); }上面这段只是示意参数化加载的思路，实际用“参数化扫描”节点更省事，不需要写循环。写代码的好处是可以把多个工况串起来自动跑，适合做批量分析。
5. 后处理不是截图：把场数据变成能写进报告的结论
5.1 派生值、积分与平均值的正确用法
云图好看，但不能直接写进报告说“温度分布均匀”。报告需要的是数字：最高温度多少、平均换热系数多少、总热阻多少。这些要用“派生值”来算，而不是靠颜色估。
“平均值”要选对数据集和选择。选面平均还是体平均，结果差很多；选错边界，平均的是空气而不是固体。算平均温度时，如果关心的是出口截面，要用“截面平均值”，并且在数据集里明确是哪个截面。算总热功率用体积分，被积函数是体热源密度，单位是 W/m³，积分出来才是 W，单位不对结果差几个数量级。
“最大值”在奇点附近不可信，前面说过。如果一定要报最大值，注明位置和网格密度，或者用附近区域的平均值代替。我一般同时报最大值和 99 分位数，后者对奇点不敏感，更能反映整体水平。
| 需求 | 工具 | 注意点 |
|---|---|---|
| 面平均温度 | 派生值-平均值 | 选对数据集和边界，别选到对称面 |
| 体总热源功率 | 派生值-体积分 | 被积函数单位要是 W/m³ |
| 沿路径温度分布 | 截线+一维绘图组 | 截线坐标单位与几何一致 |
| 随时间变化的量 | 点图或全局图 | 时间步要与求解输出对应 |
| 表面换热系数 | 派生值-平均值 | 需要先定义热通量与温差 |
5.2 一维绘图组与截线：从二维云图到可对比曲线
云图适合看整体趋势，曲线适合做对比。COMSOL 的“截线 2D/3D”可以在任意位置切一条线，把沿线数据抽出来画成一维图。做散热分析时，我习惯沿散热路径切一条线，把不同方案的 temperature vs distance 画在同一张图里，比两张云图并排有说服力得多。
截线的坐标要写对。二维轴对称模型的横坐标是 r，不是 x；三维模型的截线要指定起点和终点坐标，单位与几何一致。切完之后在“一维绘图组”里选“线图”，把多个解或者多个参数工况的数据放进去，图例自动区分。导出时选“数据”而不是截图，CSV 里带坐标和数值，后面用别的工具重绘也方便。
还有一种常用的是“点图”，监控某个特定位置随时间或参数的变化。比如芯片中心点的温度随功率变化，用“全局图”或者“点图”直接出曲线，比每次截云图省事。
5.3 数据导出与外部处理
COMSOL 的导出功能支持 CSV、文本、图像。导出场数据时要注意坐标系和单位，二维轴对称导出的是 r-z 数据，不是 x-y；三维导出的网格节点顺序和外部工具可能不兼容，读进去之前先看一眼表头。
我一般把关键曲线导出成 CSV，用表格工具做拟合和对比。这样做的好处是版本可追溯，图可以重画，数据不会因为模型改动而丢失。导出图像时分辨率给到 300 dpi 以上，云图的颜色范围手动固定，不同方案之间才能公平对比。自动颜色范围会让差异被“拉伸”得很夸张，看着差别很大，实际数值只差一点点。
5.4 常见后处理误读
第一个误读是把“表面温度”当成“体内最高温度”。表面散热好，内部可能更热，尤其是厚壁或者内热源结构，必须看体最大值。第二个误读是把稳态结果当成瞬态峰值，瞬态过程中的过冲可能比稳态高不少，做安全评估要看峰值。
第三个误读是忽略网格对后处理的影响。派生值也是在离散场上算的，网格粗的地方积分误差大。做定量结论前，确认关键区域的网格已经通过无关性检查。第四个误读是单位。温度用 K 还是 ℃，热通量用 W/m² 还是 kW/m²，导出之后一定要核对量级，别把 300 K 写成 300 ℃。
6. 几个真实会遇到的“卡壳现场”复盘
6.1 案例一：温度场算出来比热源温度还高
某开发者做一个加热块的自然对流散热，边界给了对流热通量，热源给了体热源，算出来最高温度比热源表面温度还高，明显不对。排查顺序是这样的：先看热源总量，体热源密度乘体积，和输入功率对得上；再看对流换热系数，取的是 8 W/(m²*K)，量级没问题；最后看边界选择，发现“对流热通量”被加到了对称面上，而真正的散热面反而设成了“热绝缘”。
这个坑的根因是导入几何后边界编号变了，凭记忆选的选择集对不上。修复方法是用“选择”里的“相邻”功能重新选面，或者在几何里给关键面命名，后续按名称选。修复后最高温度回到合理范围，和手算的集总参数模型差 5% 以内。
这件事之后我养成一个习惯：每个边界条件建完之后，在图形窗口里单独显示一次，确认高亮的位置和预期一致。十秒钟的动作，能省掉几个小时的排查。
6.2 案例二：网格加密后结果反而“变了”
A 同学算一个薄板的热应力，粗网格算出来最大应力 120 MPa，加密之后变成 180 MPa，再加密到 200 MPa，开始怀疑人生。问题出在板的边缘有一个直角，热应力在尖角处理论上奇异，网格越细应力越高，永远不收敛。
处理办法不是继续加密，而是判断这个尖角在真实结构里是否存在。实际零件边缘有倒角或者圆角，建一个 0.5 mm 的圆角，应力集中系数立刻变得有界，网格加密后收敛到一个稳定值。如果几何上必须保留尖角，那就报告离开尖角 1 mm 处的应力，或者用截面平均，不要报尖角点的值。
这个案例说明，网格无关性验证之前，先确认模型里有没有奇点。有奇点的地方，加密是往错误的方向使劲。
6.3 案例三：参数化扫描中途失败
做参数化扫描时，前几个参数点算得很好，到某一个点突然不收敛。常见原因是那个参数下的物理状态发生了突变，比如材料进入塑性、流体从层流过渡到湍流、接触状态改变。解决办法是把扫描步长改小，在突变区间加密采样，或者用“续解”把上一个收敛解作为初值。
还有一种可能是参数本身越界了，比如几何尺寸改到负值，或者换热系数给到零导致方程退化。扫描之前给参数加个上下限检查，能提前拦住。COMSOL 的“参数化扫描”节点里可以设“重新求解”策略，选“上一步的解”作为初值，对连续变化的参数特别有效。
我自己做扫描的习惯是先用粗步长跑一遍，看结果曲线的形状，找到拐点区间，再在拐点附近细化。一次性给很密的扫描，算得慢，而且中间失败一个点后面的全停，反而浪费时间。
6.4 一个可以照着走的最小工作流
把上面的经验收拢成一个可复现的流程，适合新手拿来练手：先建二维轴对称模型，几何只保留主特征，全局参数写清楚；物理场只选主物理场，边界条件逐条确认高亮位置；网格先用物理场控制的“常规”档跑通，再切用户控制，最大单元尺寸取特征尺寸的 1/15；稳态求解看残差是否单调下降，不收敛先查边界和初值，别急着加密；后处理先算派生值，再画截线曲线，导出 CSV 留档；最后做一次三套网格的无关性检查，关键量变化小于 2% 再下结论。
这套流程跑通一个简单案例之后，再逐步加多物理场耦合、瞬态、非线性，比一上来就啃复杂模型有效得多。COMSOL 的学习曲线陡，但陡在概念和习惯上，不在操作上。把边界设置的物理含义想清楚，把网格花在梯度大的地方，把后处理做成可复现的数字而不是截图，这三件事做到了，大部分入门级的坑都能绕过去。
我个人在实际操作中的体会是，花在建模前思考上的时间，和后面调试的时间是反比关系。一个想清楚物理画面的模型，边界条件不会多也不会少，网格不会过密也不会过疏，后处理也知道该看哪个量。反过来，急着点下一步的模型，往往要在求解器报错和结果质疑之间来回折腾。
