一种用于多物理场的原生权重边界元方法

《Engineering Analysis with Boundary Elements》：A born-weight boundary element method for multiphysics

　　

# Born权重边界元方法：一种富集的、非负的剖分单位

• 近谱精度，同时保持严格形状保真，无过冲现象
• 现有求解器的即插即用重参数化：未知量和矩阵不变
• 已通过Mie散射、精确隧穿和耦合界面进行验证
• 在非圆形边界上通过真实奇异积分再现了收敛性

## 引言

边界元方法长期以来因其降维特性而备受推崇：通过将控制偏微分方程重新表述为定义在区域边界上的积分方程，它仅对表面进行离散化，而非对体域进行离散化，并且能在不引入人工截断的情况下处理无限域和半无限域[3]、[4]、[5]。这些特性使BEM在声学、电磁学、位势理论和断裂力学中成为不可或缺的工具。然而，三个持续存在的困难限制了其向工业规模多物理场问题的扩展。

**1. 插值局限。** 多项式形状函数只能通过激进的网格细化来分辨高梯度特征——裂纹尖端奇异性、边界层、倏逝尾——这将使自由度数量急剧膨胀[6]。

**2. 稠密算子。** 影响矩阵完全填充，需要 $\mathcal{O}(N^2)$ 的存储和直接求解需要 $\mathcal{O}(N^3)$ 的计算量，这在几何复杂度尚未被充分利用之前就已经变得难以承受。

**3. 奇异积分。** 弱奇异、强奇异和超奇异积分需要精心设计的求积方案，其精确计算仍是实现错误的反复出现的来源[7]。

量子计算提供了一种对比鲜明的范式。$n$ 个量子比特的寄存器跨越 $2^n$ 维希尔伯特空间，利用叠加、纠缠和干涉的算法承诺为某些线性代数任务提供有利的标度特性[8]、[9]、[10]、[11]。尽管能够实现这些大规模收益的容错硬件尚不存在，但量子力学的数学结构可以移植到经典算法上，作为新近似方案的来源。这就是本文的前提。

**一句话概括：BWBEM是一种用于边界元的通用插值框架，基于非负的Born权重剖分单位，可以在边界积分方法的物理算子中保持不变地应用。本文其余所有内容都服务于这单一思想。** 该框架所替代的唯一经典对象是节点形状函数：取代固定的多项式帽函数，每个节点携带一个非负权重 $N_i(\mathbf{x}) = |c_i(\mathbf{x})|^2$，来源于物理自适应库[12]、[13]、[14]、[15]，从而将局部近似空间扩展以匹配预期的边界行为，并且由于它是节点值的凸组合，重构永远不会过冲。至关重要的是，这种富集既不改变边界积分方程，也不改变其核函数，也不改变未知量的数量或含义：它是现有求解器的即插即用修改。贯穿全文反复出现的三个特性——近谱精度、形状保真和无额外未知量——正是这种插值的特性，也正是图形摘要图1一览即所得出的结论。从相同的非负权重中还能引出另一项能力，它作为框架的结果呈现，而非第二个头条：算子级约束耦合通过将机械关联的自由度进行对称化影响项来关联它们，将界面连续性或裂纹面对称性等约束嵌入算子中，而非作为罚项或乘子方程附加。这在插值层面延续了双重边界元方法（DBEM）的单区域哲学[16]、[17]；其严格的断裂实现——Born权重双重边界元方法[2]——是同一框架的奇异极限（第5.4节）。

叠加和纠缠的术语仅作为一种诚实的标签保留，用于标记找到这些构造的启发式路径；这并非声称使用了量子硬件，也非声称所计算的期望值是决定性的加权和以外的任何东西。关于其地位做出了明确的决定。本文中的每一个构造都可用且确实在正文中以常规术语陈述和论证：Born权重基是一个非负的归一化（Shepard）剖分单位，其精度来源于此类方案的标准逼近理论（第6节）；耦合是算子级的对称化，其行为通过与精确耦合解直接对比来表征（第6.8节）；而"量子启发"的基本解仅仅是标准的修正Helmholtz格林函数（第3.4节）。因此，BWBEM以其自身的特点呈现——一种用于边界元的新数值近似方法、一种富集的剖分单位插值——并将量子术语视为一种启发性类比而非承重成分，正如张量网络文献[18]在经典硬件上借用量子态的线性代数结构的精神一样。读者读完本文后应将BWBEM视为边界元的一种新插值，而非现有插值的量子启发诠释：这种类比是设计方案的富有成效的启发，而本文的主张完全无需依赖它。

核心贡献是Born权重插值本身：一种边界元的非负、剖分单位富集，在无需增加未知量且不改变影响矩阵的前提下，实现了近谱（超代数）精度以及形状保真——无节点间过冲——在Laplace、Helmholtz和屏蔽算子上均已验证，且在奇异极限下于姊妹研究[2]中进行了验证。三个组件为其提供支撑和补充，并额外提供一个提案。（i）一个提出的分类法将边界行为类别映射到候选振幅函数；其中只有高斯函数在此进行了定量验证，其他函数作为猜想和未来工作陈述，而非声称已验证的结果。（ii）一个组织性观察将"量子启发"的核归入修正Helmholtz族，将Laplace和Helmholtz核恢复为极限情况；这是一种组织便利，绝非新的格林函数，因为每种体制都用其自身的精确经典核求解。（iii）可选的算子级约束耦合将界面连续性和裂纹面对称性约束嵌入算子。（iv）针对精确参考的坦率验证毫不修饰地报告了富集方法获胜之处、仅与谱基匹配之处以及分级经典网格仍具竞争力之处——与定义其优势一样仔细地定义了方法的局限。全文中，复杂的穿孔板和二维井几何仅作为机制说明保留，决定性的奇异（断裂）情况则在专门的姊妹研究中发展。

用先验解空间知识富集边界元近似并非新事物：平面波和剖分单位富集已用于短波和断裂问题，最接近的相关方法是Simpson和Trevelyan的剖分单位富集双重边界元方法[15]，更广泛地说还有扩展有限元方法[14]。那些方法向近似中添加函数，$u = \sum_a N_a u_a + \sum_a N_a \sum_k \psi_k^a k_a$，引入富集系数 $k_a^a$ 作为额外未知量。BWBEM则重参数化现有的节点插值，保持未知量数量和影响矩阵 $\mathbf{H}$、$\mathbf{G}$ 的结构不变，代价是可达的局部空间维度等于单元维度。两者相关但不等价，本文不声称对富集方法的优越性；BWBEM所贡献的是一条系统的、保持非负性的通向非多项式形状函数的路径，它是现有求解器的即插即用修改。

独立于姊妹研究的奇异场结果，本文的独立贡献是一种保持正性、近谱的剖分单位富集边界元方法，它是现有求解器的即插即用修改——且每个限定词都有其分量。对于光滑边界数据，非负Born权重以超代数速率收敛，远超低阶单元；由于是节点值的凸组合，重构不会超出其范围，因此永远不会过冲；并且在不增加未知量或改变影响矩阵的情况下重参数化节点插值，现有代码可以几乎零成本采用。本文坦承，在光滑数据上，仅凭精度即可与纯谱基匹配——而这正是关键所在。高阶和谱基以振荡为代价换取精度（Runge和Gibbs过冲），这在重构量必须遵守物理边界时是致命的：非负浓度或绝对温度、有界压力、单调或非穿透界面，或传递给下游阈值或非线性的任何场，过冲都会使其不稳定。一种同时具有谱级精度和严格形状保真的基，在无需额外未知量且不改变求解器的情况下，触及了设计空间中低阶形状保真单元和高阶振荡单元均未占据的角落。这种组合——连同对其有效之处、无效之处以及决定性优势所在之处的坦率地图——正是本文独立于[2]所确立的。

姊妹论文[2]提供了全文引用的奇异场结果（其命题3的多项式再现障碍、Williams裂纹尖端验证以及真正奇异场的决定性优势），这是本文作者目前在本期刊待刊的论文。其结果在此作为姊妹结果引用，在那里建立而非在本文建立；上述独立主张不依赖于它们，读者可参阅[2]了解其证明和验证。如果姊妹论文对读者不可用，本文的每一项定量主张——光滑标量和粘性体制、声学散射、量子隧穿、耦合界面以及非圆形研究——均建立在第6节的自包含实验之上。

第2节回顾了经典边界积分方程及该方法借鉴的量子形式化要素。第3节将BWBEM公式发展为富集插值加约束耦合。第4节描述实现。第5节展示多物理场基准测试，第7节讨论机制、成本和局限。第8节总结。

## 章节摘要

### 经典边界积分方程

考虑一个线性椭圆算子 $\mathcal{L}$ 控制域 $\Omega \subset \mathbb{R}^d$ 上的场 $u$，边界为 $\Gamma = \partial\Omega$，

$$\mathcal{L}\, u(\mathbf{x}) = 0, \quad \mathbf{x} \in \Omega,$$

其中 $\mathcal{L}$ 可以是Laplace算子 $\nabla^2$、Helmholtz算子 $\nabla^2 + k^2$、修正Helmholtz算子 $\nabla^2 - \kappa^2$，或弹性静力学Navier算子。Green第二恒等式给出边界积分方程

$$c(\boldsymbol{\xi})\, u(\boldsymbol{\xi}) + \int_\Gamma H(\mathbf{x}, \boldsymbol{\xi})\, u(\mathbf{x})\, d\Gamma(\mathbf{x}) = \int_\Gamma G(\mathbf{x}, \boldsymbol{\xi})\, q(\mathbf{x})\, d\Gamma(\mathbf{x}),$$

其中 $G(\mathbf{x}, \boldsymbol{\xi})$ 是 $\mathcal{L}$ 的基本解，$H = \partial G / \partial n$ 是其法向导数，$q = \partial u / \partial n$ 是边界通量，$c(\boldsymbol{\xi})$ 是

### 叠加插值作为富集剖分单位

在具有节点 $\mathbf{x}_1$、$\mathbf{x}_2$ 的双节点边界元上，将公式(4)中的二态与每个点 $\mathbf{x}$ 关联，将振幅 $c_1(\mathbf{x})$、$c_2(\mathbf{x})$ 视为位置的函数。引入对角可观测量 $U = u_1 |1\rangle\langle 1| + u_2 |2\rangle\langle 2|$，其特征值为节点值。定义插值场为状态 $|\psi(\mathbf{x})\rangle$ 中 $U$ 的期望值，并在正交基 $\{|1\rangle, |2\rangle\}$ 上展开，得

$$u_q(\mathbf{x}) = \langle \psi(\mathbf{x}) | U | \psi(\mathbf{x}) \rangle = (c_1^* \langle 1| + c_2^* \langle 2|)\, U\, (c_1 |1\rangle + c_2 |2\rangle) = |c_1(\mathbf{x})|^2\, u_1 + |c_2(\mathbf{x})|^2\, u_2,$$

交叉项因正交性消失，

### 实现

BWBEM流水线保留了BEM求解器的熟悉阶段——状态初始化、矩阵组装、耦合、奇异性处理和求解——并为每个阶段增加一个富集或约束步骤。该公式适用于在计算机代数系统中对核导数进行符号验证；相比之下，第6节的定量结果由专门编写的求解器产生，其奇异处理采用谱精确的Kress乘积求积

### 多物理场应用

BWBEM现在在五个基准测试上进行了检验，这些测试共同涵盖了扩散、传播、倏逝、弹性静力和粘性体制。从一开始就必须明确每项的证据地位，因为它们权重并不相同。扩散（Laplace）、传播（Helmholtz）和倏逝（屏蔽）标量体制在第6节中与精确的人工制造参考进行了定量验证，粘性（Stokes）体制也通过人工Stokeslet

### 在规范问题上的数值验证

前面的基准测试展示了公式的适用范围，但本着姊妹断裂研究[2]的精神，本文的主张必须通过与精确参考的对照实验来确定，而非通过断言。本节报告了这样的研究。其核心（第6.1–6.4节）刻意简洁——单一几何、精确的人工解、富集被简化为其本质——使得唯一变量是插值基——之后紧跟四个

### 讨论：机制、成本与局限

表13按每项测试的机制汇总了基准，第6节的实测研究确立了公式预测的两种模式。首先，与自然的第一直觉相反，富集剖分单位的最强之处不在尖锐特征，而在光滑边界场，其非负Born权重提供了保持正性的近谱精度，这是经典低阶单元在相同自由度下无法匹敌的。

其次，作为结论与未来展望，本文提出了BWBEM作为一个核心理念：一种面向边界元的、增强的非负单位分解插值——即Born权重基——它将经典的节点形函数推广为一套物理适配的振幅平方权重库，并可在边界积分法的各个算子中不加修改地直接适用。本文的其他组成部分是对这一理念的支撑，而非与之竞争。约束耦合在算子内部表达了界面连续性条件。利益冲突声明：作者声明不存在可能影响本文所报告工作的已知竞争性经济利益或个人关系。A. Portela