[论文解读] Green's functions for Neumann boundary conditions

Jerrold Franklin|arXiv (Cornell University)|Jan 29, 2012

Advanced Mathematical Modeling in Engineering参考文献 1被引用 3

一句话总结

本文推导了满足诺伊曼边界条件的偏微分方程的严格诺伊曼格林函数，纠正了教科书中普遍存在的遗漏。它引入了一个修正的格林函数，满足 $\mathcal{L}_{\mathbf{r}'}G_N(\mathbf{r},\mathbf{r}') = \delta(\mathbf{r}-\mathbf{r}') - 1/V$，确保与可解性约束 $\int \mathbf{dS} \cdot p(\mathbf{r}) \nabla\phi(\mathbf{r}) = \int \rho \, d\tau$ 一致，并提供了一个包含任意常数平均值 $\langle \phi \rangle$ 的解公式。该方法确保了有限体积问题的一致性，并解决了先前方法中的不一致性。

ABSTRACT

Green's functions for Neumann boundary conditions have been considered in Math Physics and Electromagnetism textbooks, but special constraints and other properties required for Neumann boundary conditions have generally not been noticed or treated correctly. In this paper, we derive an appropriate Neumann Green's function with these constraints and properties incorporated.

研究动机与目标

• 纠正教科书中对诺伊曼格林函数处理中的普遍错误，特别是忽略了关键的可解性约束。
• 推导一个满足非齐次微分方程 $\mathcal{L}_{\mathbf{r}'}G_N = \delta(\mathbf{r}-\mathbf{r}') - 1/V$ 和齐次诺伊曼边界条件的诺伊曼格林函数。
• 为有限体积中的诺伊曼问题提供一个一致的解框架，其中解包含一个任意常数 $\langle \phi \rangle$，其值等于 $\phi$ 的平均值。
• 证明标准方法中强制在 $G_N$ 上施加恒定法向导数的做法在非对称几何中会失败，而所提出的基于特征函数的方法则不会。

提出的方法

• 使用诺伊曼特征函数的完备集推导诺伊曼格林函数，确保展开中正交性和完备性。
• 在微分方程中引入 $-1/V$ 项，以满足可解性约束 $\int \mathbf{dS} \cdot p(\mathbf{r}) \nabla\phi(\mathbf{r}) = \int \rho \, d\tau$，该约束对解的存在性至关重要。
• 应用格林定理于修正的格林函数，推导出解公式：$\phi(\mathbf{r}) = \langle \phi \rangle + \int G_N \rho \, d\tau' - \int \mathbf{dS}' \cdot G_N p \nabla' \phi$。
• 确保格林函数满足齐次诺伊曼边界条件 $\mathbf{\hat{n}}' \cdot \nabla' G_N = 0$，在所有边界面上成立。
• 证明诺伊曼格林函数不一定是对称的，但可通过添加仅依赖于 $\mathbf{r}$ 的函数使其对称，而该操作不影响其在解公式中的作用。
• 将所提出的方法与存在缺陷的替代方法进行对比，即强制 $\mathbf{\hat{n}}' \cdot \nabla' G = 1/S$，表明该方法在一般几何（如长方体）中不成立。

实验结果

研究问题

• RQ1 为何教科书标准处理诺伊曼格林函数时未能考虑可解性约束 $\int \mathbf{dS} \cdot p \nabla\phi = \int \rho \, d\tau$？
• RQ2 能否为满足非齐次方程 $\mathcal{L}_{\mathbf{r}'}G_N = \delta(\mathbf{r}-\mathbf{r}') - 1/V$ 和齐次边界条件的有限域一致定义诺伊曼格林函数？
• RQ3 为何强制施加恒定法向导数 $\mathbf{\hat{n}} \cdot \nabla G = 1/S$ 的方法在非对称几何（如长方体）中会失败？
• RQ4 格林函数方程中 $-1/V$ 项的存在如何确保与诺伊曼问题约束的一致性？
• RQ5 解中的任意常数 $\langle \phi \rangle$ 的作用是什么，它与区域内 $\phi$ 的平均值有何关系？

主要发现

• 诺伊曼格林函数必须满足 $\mathcal{L}_{\mathbf{r}'}G_N(\mathbf{r},\mathbf{r}') = \delta(\mathbf{r}-\mathbf{r}') - 1/V$，以确保与诺伊曼问题可解性条件的一致性。
• 诺伊曼问题的解包含一个任意常数 $\langle \phi \rangle$，其值等于 $\phi$ 在体积内的空间平均值，可为简化而设为零。
• 所提出的格林函数在所有边界面上满足齐次诺伊曼边界条件 $\mathbf{\hat{n}}' \cdot \nabla' G_N = 0$，确保与格林函数零净通量的物理要求一致。
• 标准方法中强制施加恒定法向导数 $\mathbf{\hat{n}} \cdot \nabla G = 1/S$ 在一般几何中会失败，因为它对导数施加了不相容的要求，无法由光滑函数满足。
• 诺伊曼格林函数本身不具有内在对称性，但可通过添加仅依赖于 $\mathbf{r}$ 的函数使其对称，而该操作不会改变其在解公式中的作用。
• $-1/V$ 项自然地源于诺伊曼特征函数的展开，在内部问题（有限体积）中至关重要，仅在外部问题（无限域）中消失。