Voltage and current density in a rectangle with point in- and ejection via Green's functions and $q$-analysis

技术摘要：通过格林函数与 q-分析求解矩形区域内的电压与电流密度

问题陈述
本文探讨了一个经典问题，即计算在二维矩形区域内，单位电流从源点 $s$ 注入并从目标点 $d$ 喷出的电势与电流密度。虽然三维矩形束中的电势可以通过标准的变量分离法得到双傅里叶级数解，但将其降维至二维极限（当束厚度 $c \to 0$ 时）会带来显著的计算挑战。具体而言，电流密度的双傅里叶级数表示收敛极慢，导致数值评估变得不切实际。此外，论文指出，在二维连续体中，两个点电荷之间的“有效电阻”是发散的，这是由点源的零维特性引起的建模伪影。

方法论
作者采用了多阶段分析方法：

1. 三维公式化： 首先在具有诺依曼（Neumann）边界条件（即侧壁无电流流出）的三维矩形平行六面体中求解该问题。利用变量分离法，将电势表示为包含 $x$ 和 $y$ 的余弦项以及 $z$ 的双曲余弦项的双傅里叶级数。
2. 维度缩减： 通过取厚度 $c \to 0$ 的极限来推导二维解。这需要一个重整化论证（$\rho I / c = 1$）来处理点注入的奇异性，从而有效地将问题转化为一个带有相反 $\delta$ 函数的泊松方程。
3. 拟格林函数： 意识到对于具有单个 $\delta$ 源的诺依曼边界条件，标准的格林函数并不存在（因为需要满足零净电荷要求），作者构建了一个“拟格林函数”，代表两个相反 $\delta$ 源之差。该函数最初被表述为单傅里叶级数，其收敛速度比双级数快，但在源/汇线处仍表现出收敛缓慢及导数不连续的问题。
4. q-分析与椭圆函数： 核心的方法论创新在于将单傅里叶级数转化为 $q$-级数。通过利用高斯多项式和雅可比 $\theta$ 函数的性质，作者将无穷级数重新表示为一个涉及参数 $q = e^{-2\pi b/a}$ 的无穷乘积。这种变换利用了 $q$-乘积的快速收敛特性（其中 $q$ 通常非常小，例如对于正方形，$q = e^{-2\pi} \approx 1/535$）。

主要贡献与结果

• 精确解析公式： 论文推导出了矩形内电势和电流密度大小的精确、闭式解析公式。电势通过一个拟格林函数 $K(x, y; \xi, \eta)$ 表示，该函数被分解为一个线性项和一个 $q$-函数 $T_q$。
• 计算效率： 主要结果证明了使用 $q$-级数可以极高效率地计算电流密度。论文表明，对于正方形区域，该级数的收敛速度极快，以至于仅需单项（$k=0$）即可提供在视觉和数值上与包含多项的级数无异的精度。这与传统双傅里叶级数的缓慢收敛形成了鲜明对比。
• 解决发散问题： 论文严谨地证明了在二维连续体中，两个点源之间的有效电阻是无穷大的（发散的）。它澄清了这是点源模型（$\delta$ 函数）的结果，而非物理现实，并指出通过对有限面积进行建模可以解决这种发散。
• 与特殊函数的联系： 该工作明确地将矩形中的拉普拉斯方程解与椭圆函数理论联系起来，特别展示了该解如何与雅可比 $\theta$ 函数及高斯多项式相关联。

意义与主张
论文将自身定位为从经典傅里叶分析向更高级的 $q$-分析领域的探索之旅。作者声称，所推导的 $q$-级数公式是解决该问题“精度最高且计算量最小”的方法。其意义在于：

• 为一个标准方法（双傅里叶级数）因收敛缓慢而失效的问题，提供了计算上可行的解决方案。
• 在有界区域内的静电学与椭圆函数理论之间搭建了一座数学桥梁，呼应了昂萨格（Onsager）在二维伊辛模型上的工作。
• 强调了二维连续体网络中点源模型的局限性（特别是关于有效电阻），以及为了物理一致性采用有限面积建模的必要性。

作者并未提出新的实验装置或超出随机几何图（HRGG）分析背景之外的广泛应用；在该背景下，电流密度被推测与节点位于最短路径上的概率成正比。这项工作始终保持在解析数学与理论物理的范畴内。