企业官网建设流程全解析

COMSOL多物理场二次开发教程（11）：自定义物理场（上）——系数型 PDE 与通式 PDE 的 API 装配

版本与事实声明

• 版本锚点：COMSOL Multiphysics® 6.3。
• 已验证的官方事实：系数型 PDE 用于"把众多知名 PDE 写成系数形式"，其一般形式可含对时间与空间最高到二阶的导数（官方原文口径）；弱形式推导中出现的通量表达为Γ = −c∇u − αu + γ；官方说明"可向组件添加的任意形式PDE 接口都能继续添加Coefficient Form、General Form 或 Weak Form方程节点"，且 Coefficient/General/Weak Form PDE也提供边界、边、点版本；可用**辅助因变量 + Weak Contribution（PDEs and Physics）**追加额外弱方程。
• PDE 接口与系数节点的具体类型字符串未逐字确证，正文写作<COEFF_PDE_TYPE>/<GENERAL_PDE_TYPE>等占位符，并给出用Record Method在本机取得真名的动作（铁律 1、底账 U1/U3）。
• 示例方程与数值均为教学用途，不代表任何标准规定。

一句话结论：自定义物理场有三个能力阶梯，按"够用即止"的原则选：系数型 PDE（Coefficient Form PDE）用一组固定系数（c / α / β / a / γ / f等，具体属性名以官方文档为准）表达二阶线性/弱非线性方程，能覆盖扩散、对流、反应、波动、Helmholtz 等绝大多数工程方程；通式 PDE（General Form PDE）允许把方程写成任意守恒形式以表达强非线性与自定义通量；弱形式 PDE / 弱贡献才提供完全自由（第 12 篇）。能写系数型就不写通式，能写通式就不写弱形式。

〇、本篇要解决的认知问题

• Q1：什么时候必须自己写 PDE，而不能用官方物理场接口？
• Q2：系数型 PDE 的那组系数分别代表什么物理量？怎么把《传热学》《流体力学》里的方程搬进来？
• Q3：Coefficient / General / Weak 三者的关系是"三选一"还是"可叠加"？
• Q4：降维 PDE 接口（边界/边/点版本）解决什么特殊问题？
• Q5：怎么确认"我写的 PDE 真的等价于我想解的那个方程"？

一、机制解析

1.1 价值锚点：PDE 接口是"逃生门"，不是"起点"

先给一条能立刻省下大量时间的判断（经验法则）：

能找到的官方接口，绝不用 PDE 重写。
原因不是"难写"，而是官方接口附带了你没意识到的东西：变量命名、守恒结构、边界条件的自然形式、雅可比稀疏结构、与其它物理场的耦合变量、单位一致性与方程视图。

什么时候才必须用 PDE（三个真实场景）：

1. 方程超出官方接口的假设：官方接口假设了某种本构关系，而你的模型需要另一种（如非标准源项、非常规通量）；
2. 教学/研究用途：需要显式写出控制方程以验证格式或做数值实验；
3. 抽象耦合：你要把一个"只存在于论文里的方程"接到已有的温度场/浓度场上（这时 PDE 是最短的路径）。

1.2 系数型 PDE：一组系数读懂一张方程表

官方对系数型 PDE 的定位是：许多源自物理接口与其它领域的 PDE 都能写成"系数形式"——即用固定位置的一组系数去承载一个最多含时间与空间二阶导数的方程。

其结构（通量的Γ = −c∇u − αu + γ这一部分有官方原文佐证；完整系数项与界面/API 属性名以官方文档为准）：

ea·∂²u/∂t² + da·∂u/∂t + ∇·(−c∇u − αu + γ) + β·∇u + a·u = f └─ 二阶时间 ─┘ └─ 一阶时间 ─┘ └── 通量散度（空间二阶） ──┘ └对流┘ └吸收┘ └源项┘

把这组系数翻译成物理，是理解系数型 PDE 的全部关键：

示例：把稳态对流-扩散-反应方程搬进来

工程方程：

u_vel · ∇c − D ∇²c + k·c = 0 （稳态，1D 简化为）

对照表：β = u_vel、c = D、a = k、f = 0、da = 0（稳态）、ea = 0、α = 0、γ = 0。于是"写 PDE"就变成了"填一张系数表"——这就是系数型 PDE 的全部方法论。

1.3 通式 PDE：当方程不能写成系数形式

通式 PDE（General Form PDE）允许你把方程写成守恒形式（通量 + 源项的自定义组合），从而表达：

• 通量非线性依赖于 u 或其梯度（如c = c(|∇u|)，非牛顿黏度、温度依赖导热的强形式）；
• 源项以任意非线性形式出现（如 Arrhenius 型源、带阈值/切换的速率）；
• 需要把方程按"守恒律书写习惯"呈现以便与文献逐项对照。

选择判据（最佳实践）：

如果方程的每个系数都能写成"参数或已求解量的简单表达式"，用系数型；如果通量本身需要自定义函数形态，用通式；如果连"通量/源项"的结构都不适用（如鞍点问题、要加约束的弱项），用弱形式（第 12 篇）。

1.4 三者的关系：可叠加，不是三选一

官方的关键说明（原文口径）：向组件添加任意形式的 PDE 接口后，还能继续添加 Coefficient Form、General Form 或 Weak Form 方程节点。

这意味着你可以构建混合系统，例如：

• 主方程用系数型（清晰、可读）；
• 某一个源项需要特殊结构时，为额外的因变量加一个通式方程节点；
• 需要补充一个约束时，用辅助因变量 + **Weak Contribution（PDEs and Physics）**加入额外弱方程（官方明示这条路径可用于处理薄层等情况）。

给自动化开发的启示：PDE 接口下的"方程节点"也是模型树节点，因此它们同样是可被 API 枚举、修改、扫描的对象——这决定了"自定义物理场"能不能被参数化扫描与优化（第 13/14 篇）。设计自定义物理场时，把可变部分做成参数或表达式，而不是写死在方程节点里，是能否自动化的分水岭。

1.5 降维 PDE：边界/边/点上的方程

官方说明：Coefficient Form PDE 还提供Lower Dimensions子菜单下的三个变体——Coefficient Form Boundary PDE (cb)、Coefficient Form Edge PDE (ce)、Coefficient Form Point PDE (cp)；官方另有一节专门讲"在边界、边、点上建模"。

它们的用途不是"简化"，而是在低维流形上求解独立方程：

• 边界 PDE (cb)：在 3D 体的表面上解一个 2D 方程（如表面吸附/表面反应、壳层扩散）；
• 边 PDE (ce)：在面上的曲线上解 1D 方程（如线反应器、管网）；
• 点 PDE (cp)：在边界的点上解 0D 方程（常微分/代数方程，如局部平衡）。

工程价值：避免为薄层做三维网格。一个 1 μm 的催化层如果用 3D 网格去解析，代价可能是几十倍的单元数；而"表面方程 + 体域方程"的两域耦合能以极小代价捕获同样物理。这也是官方提到"用辅助因变量 + 弱贡献处理薄层"的同一动机。

1.6 怎么验证"我写的 PDE 是对的"

三条验证手段（最佳实践，按成本递增）：

1. 退化到已知解：把系数设成能解析求解的常系数情形（1D 稳态扩散 = 线性分布），检查数值解与该解析解的相对误差；
2. 方程视图（Equation View）：官方说明它会显示弱表达式（Weak expression）与积分阶（Integration order），并且修改弱表达式会覆盖物理场接口给出的方程并出现警告。用它来"看见"你的系数最终被组装成什么弱式——这是把"我填的系数"和"实际求解的方程"对上的唯一可靠办法；
3. 守恒/对称性检查：无源时总通量应为零（守恒）；对称问题解应对称（对称性）。这两条能抓住"某个系数填错位置"的错误（因为填错位置常常仍能出图）。

二、完整代码与逐行剖析

代码 2-1：自定义对流-扩散-反应 PDE 的装配（Java，占位符需替换）

// ===== 用系数型 PDE 写稳态对流-扩散-反应：u_vel*dc/dx - D*d2c/dx2 + k*c = 0 =====// 占位符说明：<COEFF_PDE_TYPE> 等必须用 Record Method 在本机取得真名（铁律 1）clearModel(model);model.component().create("comp1",true);model.component("comp1").geom().create("geom1",1);model.component("comp1").mesh().create("mesh1");model.component("comp1").geom("geom1").create("i1","Interval");model.component("comp1").geom("geom1").feature("i1").setIndex("coord",0.20,1);// 20 cmmodel.component("comp1").geom("geom1").run();// --- 参数表：把"可扫描的部分"全部参数化（决定能否自动化） ---model.param("default").set("D_AB","1e-9[m^2/s]");// 扩散系数model.param("default").set("u_vel","0.01[m/s]");// 对流速度model.param("default").set("k_rxn","1e-3[1/s]");// 一阶反应速率model.param("default").set("c_in","10[mol/m^3]");// 入口浓度model.param("default").set("L_dom","0.20[m]");// 域长（用于核对）// --- 添加 PDE 接口（类型字符串用录制取得） ---model.component("comp1").physics().create("cdeq","<COEFF_PDE_TYPE>","geom1");// [1]// --- 设置系数（属性名以官方文档/录制结果为准） ---// 目标方程：β·∇c + ∇·(−c∇u) + a·u = f （稳态，da=ea=0，α=γ=0）// 对应：c 系数 <- D_AB ；β 系数 <- u_vel ；a 系数 <- k_rxn ；f <- 0// 以下为"按录制结果替换"的骨架（属性名字符串必须核对，不可猜）：// model.component("comp1").physics("cdeq").feature("cfeq1").set("c", "D_AB"); // [2]// model.component("comp1").physics("cdeq").feature("cfeq1").set("be", new String[]{"u_vel"}); // [3]// model.component("comp1").physics("cdeq").feature("cfeq1").set("a", "k_rxn");// model.component("comp1").physics("cdeq").feature("cfeq1").set("f", "0");// model.component("comp1").physics("cdeq").feature("cfeq1").set("da", "0"); // 稳态 -> 0// --- 因变量名与初值（决定后处理里这个场叫什么） ---// model.component("comp1").physics("cdeq").prop("ShapeProperty").set("c", "cA"); // [4] 以官方文档为准debugLog("请在界面上确认因变量名（后处理表达式要用它）");// --- 边界条件：入口浓度 + 出口（用录制取得边界特征类型） ---// model.component("comp1").physics("cdeq").create("in1", "<DIRICHLET_TYPE>", 0); // [5] 点// model.component("comp1").physics("cdeq").feature("in1").selection().set(1);// model.component("comp1").physics("cdeq").feature("in1").set("r", "c_in - cA"); // [6] Dirichlet 形式// --- 研究 ---model.study().create("std1");model.study("std1").create("stat","Stationary");model.study("std1").feature("stat").setSolveFor("/physics/cdeq",true);// [7]model.study("std1").createAutoSequences("all");model.sol("sol1").runAll();// --- 验证：退化检查（把 u_vel 与 k_rxn 置零 -> 纯扩散应得线性分布） ---model.param("default").set("u_vel","0[m/s]");model.param("default").set("k_rxn","0[1/s]");model.sol("sol1").runAll();// [8]debugLog("纯扩散时应为线性浓度分布；非线性即装配有误");

逐行剖析

• [1]<COEFF_PDE_TYPE>：PDE 接口类型字符串——本文刻意留白。用三步录制法取得（第 07 篇 1.3 节）。界面上"数学 > 系数型 PDE"添加后录制一次即可。
• [2]~[3] 系数赋值：属性名字符串（如"c"、"be"、"a"、"f"、"da"）必须核对。注意α、β这类希腊字母系数在不同版本/界面里的内部属性名不是简单等于希腊字母，猜测必然失败——这就是为什么本段用注释占位而不是直接给"看起来对"的字符串。
• [3] 注意β（对流）是向量：在 1D 里也只传数组new String[]{"u_vel"}；3D 里要传三个分量。维度不同、数组长度不同，这是最容易被忽略的一致性问题。
• [4] 设定因变量名：这决定了后处理里写cA还是c。自定义物理场最容易被忽略的一步——名字没定好，后面所有导出/绘图表达式都会写错。先在界面上确认，再写进脚本。
• [5][6] 边界条件：Dirichlet 类边界在 PDE 接口里常以"残差形式"给出（如r = c_in - cA，即令残差为零）。这个写法与官方物理场接口里"直接给温度值"的写法语义不同，是 PDE 使用者最常见的困惑点之一；具体属性名以录制结果为准。
• [7]setSolveFor("/physics/cdeq", true)：自定义物理场同样是物理场，同样必须显式声明被求解（沿用第 07 篇的教训）。
• [8] 退化检查：把对流与反应置零后重跑，验证是否退化为线性分布。这是"我写的 PDE 对不对"的最短验证路径，成本极低，收益极高——建议每个自定义 PDE 都先跑这一步再进业务。

代码 2-2：从"方程视图"反查组装结果（操作清单）

目标：确认系数型 PDE 最终被组装成什么弱表达式 步骤（界面操作，不可替代）： 1) 打开模型，展开 组件 > 系数型 PDE > 因变量/方程节点 2) 右键该节点 -> 显示 -> 方程视图（Equation View） 3) 查看两列： - Weak expression : 该节点的弱表达式（你的系数最终如何进入弱式） - Integration order : 该表达式的积分阶（多项式最高可被无系统误差积分） 4) 修改 Weak expression 会【覆盖】物理场接口给出的方程，并出现警告标记 -> 这与"改系数"是两条完全不同的路径，不可混用

剖析

• 方程视图是"可读回"的唯一官方窗口：官方明确它显示弱表达式与积分阶，并且修改它会覆盖原方程并给出警告。这意味着：它既能用来核对，也能用来接管。
• 接管要慎重：一旦在方程视图里改了弱表达式，该节点就与"系数设置"脱钩了——后续改系数不再生效。在自动化脚本里尤其危险，因为脚本可能反复改系数却看不到效果。建议：要么全程用系数，要么明确用弱形式节点（第 12 篇），不要混用。
• 积分阶（Integration order）是很多人不知道的实用信息：它告诉你"该表达式的多项式最高能被无系统误差地积分"。当你加入高梯度/强非线性源项时，若积分阶不足会产生系统性误差（表现为网格加密后解持续漂移）。它是诊断"数值收敛但解不对"的一把钥匙。

三、常见报错与排查

报错 3-1：create("cdeq","<猜的类型>","geom1")报未知类型。
现象：PDE 接口创建失败。根因：类型字符串猜测错误（这是底账 U1/U3明确列为不可断言的一类信息）。解法：三步录制法取真名。不要继续试第二个猜测值。

报错 3-2：系数设了但解完全没变。
现象：改系数后结果不变。根因：三类——① 属性名不对（set了一个不存在的属性，可能静默无效果）；② 该节点的方程被"方程视图"改写过，系数已不生效；③ 改的是另一个因变量的方程节点。解法：用方程视图查看该节点的弱表达式是否反映了你的系数；确认setSolveFor里的物理场 tag 与你要改的接口一致。

报错 3-3：自定义场的后处理表达式报"未定义变量"。
现象：导出/绘图写cA找不到。根因：因变量名没设对（默认名可能是u、u2、c之类）。解法：在界面上确认因变量名（PDE 接口的因变量设置），或用 API 读取（属性名以官方文档为准）；也可先手工导一次表看表头里叫什么。

报错 3-4：网格加密后解持续漂移。
现象：解不收敛于一个值，越加密越变。根因：积分阶不足（高梯度/强非线性源项被欠积分），或方程本身在离散意义下不守恒（系数填错位置）。解法：在方程视图里查积分阶并按需提高；同时用退化算例与守恒检查排除"系数填错位置"。

报错 3-5：薄层模型网格量爆炸。
现象：3D 网格单元数远超预期，求解极慢。根因：对薄层（催化剂层、涂层、膜）做三维解析。解法：改用降维 PDE——用 Coefficient Form Boundary PDE（cb）/Edge PDE（ce）/Point PDE（cp）在低维流形上解方程，或用辅助因变量 + Weak Contribution 处理薄层（官方明示的两条路径）。

四、动手练习

• 练习 1（取得真名）：在界面上添加一个系数型 PDE 接口与一个方程节点，用录制取得它们与系数字符串的真名。判定：录制得到的create(...)与set(...)语句能原样在 Java Shell 跑通，且界面显示的系数值随之改变。
• 练习 2（系数翻译）：为稳态对流-扩散-反应方程u_vel·∇c − D∇²c + k·c = 0填写完整的系数对照表（c/α/β/a/γ/f/da/ea八行）。判定：表格八行齐全，其中c=D、β=u_vel、a=k、f=0、da=ea=α=γ=0五项正确。
• 练习 3（退化验证）：按代码 2-1 的 [8] 步把对流与反应系数置零，验证退化为 1D 纯扩散的线性浓度分布。判定：浓度剖面与线性函数的最大相对偏差 < 1%（可用导出数据在 Python 里算），并解释为什么这能证明"系数装配正确"。
• 练习 4（能力阶梯，思考题）：说明为什么"能写系数型就不写通式、能写通式就不写弱形式"。验证要点（至少 3 点）：① 系数型的所有系数都是可枚举、可扫描的属性，自动化成本最低；② 通式提供自定义通量/源的自由度但牺牲了系数的可枚举性，参数化更麻烦；③ 弱形式完全自由但需要自己保证守恒与雅可比结构，且改弱表达式会覆盖接口方程（官方警告），维护成本最高。

五、小结与下一篇预告

本篇打开了"自定义物理场"这扇门：PDE 是逃生门而非起点（能找到的官方接口绝不用 PDE 重写）；系数型 PDE 的核心是"把方程翻译成一张系数表"（c扩散、β对流、a吸收、f源、da/ea时间项，属性名以官方文档为准）；系数→通式→弱形式是能力阶梯，按需上升；三者可叠加（官方明示可继续添加 Coefficient/General/Weak 方程节点），也能用辅助因变量 +Weak Contribution追加额外弱方程；降维 PDE（cb/ce/cp）用于薄层建模以省下巨量网格；验证靠退化算例 + 方程视图 + 守恒/对称性三件套，其中方程视图会显示弱表达式与积分阶，修改弱表达式即覆盖接口方程并触发警告。

第 12 篇《自定义物理场（下）》将把阶梯的最后一级讲透：弱形式 PDE 与弱贡献。我们会看清弱形式推导为什么能"降低导数阶数、让自然边界条件显式化"（官方给出的两条优势），Γ = −c∇u − αu + γ这个通量在变分里如何变成边界项，以及怎样用"弱贡献"给已有物理场追加一个方程——并给出泊松与对流扩散两个方程的弱形式写法与验算流程。

本篇认知问题回显（FAQ）

Q1：什么时候必须自己写 PDE，而不能用官方物理场接口？
A：只有三种场景：方程的假设超出官方接口（需要非标准本构、通量或源项）；教研用途需显式写出控制方程以验证格式；把论文里的抽象方程接到已有温度场/浓度场上。找到的官方接口绝不用 PDE 重写，因为官方接口附带变量命名、守恒结构、自然边界条件、雅可比稀疏结构、单位一致性与方程视图等隐含资产。

Q2：系数型 PDE 的系数分别代表什么，怎么把工程方程搬进来？
A：结构为"二阶时间项 ea + 一阶时间项 da + 通量散度 ∇·(−c∇u − αu + γ) + 对流 β·∇u + 吸收 a·u = f"。其中 c 是扩散/传导系数（导热 k、扩散 D）、β 是对流速度、a 是吸收或一阶反应速率、f 是体积源项、da 是时间一阶项（传热中为 ρ·cp）、ea 是时间二阶项；完整系数项与 API 属性名以官方文档为准。把方程写成该结构后逐项对应填表即可。

Q3：Coefficient / General / Weak 三者的关系是三选一还是可叠加？
A：可叠加。官方说明向组件添加任意形式的 PDE 接口后，还能继续添加 Coefficient Form、General Form 或 Weak Form 方程节点；Coefficient/General/Weak Form PDE 也提供边界、边、点版本。因此可以主方程用系数型、特殊源项用通式、额外约束用辅助因变量 + Weak Contribution。

Q4：降维 PDE 接口解决什么特殊问题？
A：解决薄层建模的网格代价问题。Coefficient Form Boundary PDE (cb)、Edge PDE (ce)、Point PDE (cp) 让你在体表面、曲线、点上求解独立方程（表面反应、线反应器、局部平衡等），从而避免为 1 μm 级薄层做三维网格解析，代价可相差数十倍单元数。

Q5：怎么确认我写的 PDE 真的等价于想解的方程？
A：三件套：① 退化验证——把系数设为能解析求解的常系数情形（如纯扩散应得线性分布），检查数值解相对误差；② 方程视图——它显示 Weak expression 与 Integration order，可核对系数最终组装成的弱表达式，且修改弱表达式会覆盖接口方程并出现警告；③ 守恒与对称性检查——无源时总通量应为零、对称问题解应对称，用于捕捉"系数填错位置"却仍能出图的错误。