[論文レビュー] Green's functions for Neumann boundary conditions

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

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

ひとこと要約

本稿は、境界条件がノイマン型の偏微分方程式に対して、教科書で広く見られる省略を是正する厳密なノイマン・グリーン関数を導出する。修正されたグリーン関数は、$ olimits\mathcal{L}_{\mathbf{r}'}G_N(\mathbf{r},\mathbf{r}') = \delta(\mathbf{r}-\mathbf{r}') - 1/V$ を満たし、解の存在条件 $ olimits\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$ と同次境界条件を両方満たすノイマン・グリーン関数を導出すること。
• 有限体積領域におけるノイマン問題に対して一貫した解法フレームワークを提供すること。この場合、解には $\phi$ の空間平均に等しい任意の定数 $\langle \phi \rangle$ が含まれる。
• 標準的手法である $G_N$ に対して定常法線微分を $\mathbf{\hat{n}}' \cdot \nabla' G = 1/S$ と仮定する方法が、非対称な幾何形状では失敗することを示すこと。

提案手法

• 完全なノイマン固有関数の集合を用いてノイマン・グリーン関数を導出し、展開における直交性と完全性を保証する。
• 解の存在に不可欠な制約条件 $\int \mathbf{dS} \cdot p(\mathbf{r}) \nabla\phi(\mathbf{r}) = \int \rho \, d\tau$ を満たすために、微分方程式に $-1/V$ 項を組み込む。
• 修正されたグリーン関数に対してグリーンの定理を適用し、解の公式 $\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$ の体積領域における空間平均に等しく、簡略化のため 0 に設定できる。
• 提案されたグリーン関数は、すべての境界面で同次ノイマン境界条件 $\mathbf{\hat{n}}' \cdot \nabla' G_N = 0$ を満たし、グリーン関数自体にゼロのネットフラックスを要請する物理的要件と整合する。
• 標準的手法（$\mathbf{\hat{n}} \cdot \nabla G = 1/S$ を仮定）は、一般の幾何形状では失敗する。これは微分の要件が整合せず、滑らかな関数では満たせないからである。
• ノイマン・グリーン関数は本質的に対称的ではないが、その作用に影響を及げない $\mathbf{r}$ の関数を加えることで対称化可能である。
• $-1/V$ 項はノイマン固有関数展開から自然に導かれ、内部問題（有限体積）において不可欠であり、外部問題（無限領域）ではのみ消える。