7.1 自定义 PDE 与物理场开发器
7.1 自定义 PDE 与物理场开发器
本节摘要：当软件自带的物理场接口没有你要的方程，数学接口家族给你三级台阶：系数型偏微分方程、广义形式、弱形式节点；物理场开发器再把调通的方程封装成可分发的自制接口。本节用系数型接口复算一维导热并完成对账，走通从裸方程到正式接口的全路径。
什么时候需要自己写方程
判断标准就一条：你要的物理在物理场接口列表里找不到对应物，或现有接口的表达力不够。典型场景有几类：非标准的输运问题（比如带温度漂移项的粒子输运）、企业自有的等效本构（专利材料模型）、简化代理模型（把三维问题降维成标量场方程）。反过来的告诫同样重要：接口列表里已有的东西别自造——自造接口的全部维护成本、边界条件、文档都由你承担，而已有接口经过千锤百炼。先翻接口列表，再进数学接口，这个顺序别倒。
三级台阶：系数型、广义形式、弱形式
数学接口家族里最常用的是系数型偏微分方程接口，它把一大类二阶线性与拟线性方程统一写成：ea·utt + da·ut + ∇·(−c∇u − αu + γ) + β·∇u + a·u = f。名字唬人，对应关系极简单：c 是扩散系数，β 是对流项，a 是反应（吸收）项，f 是源，da 与 ea 管瞬态与二阶瞬态。第 4.4 节的扩散-反应方程里，c = D、a = k、f = 0——一个标量输运问题用这组系数就能完整表达。广义形式把对流项写成通量的一部分，对流主导（佩克莱数大）时数值性质更好；弱形式节点则直通 5.1 讲的底层，直接输入弱表达式，试验函数 test(u) 显式登场，表达力最强、门槛也最高。
动手复算一维稳态导热对账，把台阶踩实。新建一维组件，几何一条 0.1 m 的区间；数学接口→系数型偏微分方程（c），研究选稳态。设置里逐项填：c 填 k_solid（先在参数里建 238[W/(m*K)]），a 填 0，f 填 Q_vol（建参数 1e6[W/m^3] 均匀内热源）。边界：左端"狄利克雷边界条件"填 293.15 K，右端默认（系数型接口的默认零通量与传热接口一致，出处正是 5.1 的自然边界条件）。网格两百个单元足够，计算后对照解析解：均匀源、一端定温一端绝热的稳态解是抛物线 T(x) = T0 + f·(L·x − x²/2)/c，最高温在绝热端，T_max = T0 + f·L²/(2c) = 293.15 + 1e6 × 0.01/476 ≈ 294.15 K。派生值一读，仿真与解析解在小数点后第三位握手——你亲手用"裸方程"复现了传热接口的全部行为，从此这个接口对你不再是黑盒。
物理场开发器：把方程变成接口
方程在某个模型里调通了，下一个问题是怎么让同事不用懂细节也能用。物理场开发器（随核心版提供）就是干这个的：它用一套定义界面让你创建自制物理场接口——声明场变量、把系数与源做成"特性节点"、配上图标与默认设置，最终导出为一个可用文件，同事双击安装后在物理场列表里看到与你手写接口完全同款的东西。定义语言是基于 XML 的声明式描述，不用写代码，工作量集中在把方程拆成特性与默认值两件事上。适用场景有明确的边界：团队层面反复使用的专用模型（企业本构、标准化的等效方法）值得封装；一次性的研究模型，写在数学接口里带注释就够了，封装纯属过度工程。
⚠️ 自定义方程最容易翻车的地方不是表达式写错，而是边界条件与量纲：自造接口没有现成的"温度""接地"这类语义化边界，一切边界都要你自己在方程层面定义，量纲自查要当成固定动作。
本节要点回顾
- 先翻接口列表再进数学接口：已有接口千锤百炼，自造接口自己养。
- 三级台阶按表达力排：系数型覆盖标量输运主力，广义形式利于对流主导，弱形式直通底层。
- 复算对账是脱黑盒的最短路：系数型接口复现一维导热，解析解 294.15 K 与仿真三位小数握手。
- 物理场开发器面向团队复用：XML 声明式封装、免代码，一次性模型不值得封装。
- 自定义方程的翻车点在边界语义与量纲，自查清单化才能稳定交付。
自定义方程的收敛自查
自定义方程没有内置接口那套成熟的初始化与容差策略，收敛问题要自己兜底。上手顺序建议固定：先解一个去掉源项或非线性项的退化版本，确认系数矩阵本身健康，再全量打开观察残差在哪一步开始失控——第一步就发散，多半是初值离物理解太远；中途发散，则常见于边界条件互相矛盾或量纲不统一。所有系数尽量写成带物理量纲的表达式，让单位检查替你把关，这一条能拦住大半的翻车。排查无果时，把方程两端的量纲逐项手推一遍，往往比反复调容差更快找到病根。
自定义 PDE 的调试技巧：初值和边界条件的合理性是自定义方程能否收敛的关键——系数取值偏离物理量级（比如扩散系数给了 10⁻¹² 而不是 10⁻⁵）会让求解器在第一步就卡住或产生非物理的解。排查方法：先用简化方程（比如把非线性项删掉只留线性扩散）确认求解器本身工作正常，再逐项加回复杂度。另外，COMSOL 的物理场开发器（Physics Builder）可以把一个自定义 PDE 封装成一个可复用的用户自定义物理场接口，让组态体验与内置物理场一致——如果你需要在多个项目中反复使用同一个自定义方程，花时间做一次 Physics Builder 封装的回报是长期性的。