代用電荷法

代用電荷法（だいようでんかほう、テンプレート:Lang-en-short）または電荷重畳法は数値解析手法の一つである。基本解近似解法（テンプレート:Lang-en-short）ともいう。構造力学や電界計算の分野で広く使われている。偏微分方程式に対するメッシュフリー法であり、有限個の基本解を境界条件が満たされるように重ね合わせて近似解を構成する。誤差が境界で最大値に達する性質があり、誤差評価を容易にしている。原理が簡単で、プログラムが容易、高速、高精度であるが、非線型問題には適用できない。テンプレート:要検証範囲、その後日本で大きく研究が進んだ。宅間董により種々の電界計算に応用され、村島定行により汎用の解析法として確立された。

導入

例として、ディリクレ境界条件の2次元ラプラス方程式 $$\{\begin{array}{cc}\mathrm{\Delta }u(𝒙)=0& (𝒙\in \mathrm{\Omega }),\\ u(𝒙)=f(𝒙)& (𝒙\in \partial \mathrm{\Omega })\end{array}$$ に代用電荷法を適用する。ただし、テンプレート:Mathはテンプレート:Math上の有界な単連結領域とする。

代用電荷法は解を基本解の重ね合わせ（線型結合）で近似する。2次元ラプラス方程式であれば対数ポテンシャル $$G(𝒙,𝒔)=-\frac{1}{2\pi }\log |𝒙-𝒔|$$ を重ね合わせて $$u(𝒙)\approx u_N(𝒙)=\sum _{i=1}^N{Q}_{i}G(𝒙,{𝒔}_i)$$ とすることが多い。式中のテンプレート:Mvarは電荷点（charge point）と呼ばれ、テンプレート:Mathの外部から適当に選ばれる。

テンプレート:Mvarの値は選点法により決定される。すなわち、境界テンプレート:Mathから拘束点（collocation point）と呼ばれる点テンプレート:Mvar（テンプレート:Math）を適当に選び、これらの点において近似解が拘束条件 $$u_N({𝒙}_k)=f({𝒙}_k)(k=1,2,\dots ,N)$$ を満たすようにする。このとき、テンプレート:Mvarは線型方程式系 $$\left(\begin{array}{cccc}G({𝒙}_1,{𝒔}_1)& G({𝒙}_1,{𝒔}_2)& \cdots & G({𝒙}_1,{𝒔}_N)\\ G({𝒙}_2,{𝒔}_1)& G({𝒙}_2,{𝒔}_2)& \cdots & G({𝒙}_2,{𝒔}_N)\\ ⋮& ⋮& \ddots & ⋮\\ G({𝒙}_N,{𝒔}_1)& G({𝒙}_N,{𝒔}_2)& \cdots & G({𝒙}_N,{𝒔}_N)\end{array}\right)\left(\begin{array}{c}{Q}_1\\ Q_2\\ ⋮\\ Q_N\end{array}\right)=\left(\begin{array}{c}f({𝒙}_1)\\ f({𝒙}_2)\\ ⋮\\ f({𝒙}_N)\end{array}\right)$$ の解となる。

代用電荷法の特徴として、境界における誤差で内部における誤差を評価できることが挙げられる。問題の解テンプレート:Mathが存在し、さらにテンプレート:Mathが閉包テンプレート:Mathにおいて連続ならば、調和関数の最大値原理から誤差に関する等式 $$\underset{𝒙\in \mathrm{cl}(\mathrm{\Omega })}{\max }|u(𝒙)-u_N(𝒙)|=\underset{𝒙\in \partial \mathrm{\Omega }}{\max }|f(𝒙)-u_N(𝒙)|$$ が成立するテンプレート:Sfnテンプレート:Sfn。

変種

室田の不変スキーム

室田の不変スキームは、対数ポテンシャルの線型結合に定数項テンプレート:Mathを加えた近似解 $$u_N(𝒙)=Q_0-\frac{1}{2\pi }\sum _{i=1}^N{Q}_i\log |𝒙-{𝒔}_i|$$ を、拘束条件に $$\sum _{i=1}^N{Q}_i=0$$ という制約を追加して構成する手法であるテンプレート:Sfn。

ラプラス方程式の解は座標系のスケール変換テンプレート:Mathと境界条件の平行移動テンプレート:Mathに対して不変だが、通常の代用電荷法で構成される近似解は不変にならない。室田の不変スキームで構成される近似解は、この不変性を満たす点に特徴がある。

参考文献

• テンプレート:Cite journal
• テンプレート:Cite book
• テンプレート:Cite journal
• テンプレート:Cite journal

関連論文

• テンプレート:Cite journal
• テンプレート:Cite journal
• テンプレート:Cite journal
• テンプレート:Cite journal
• テンプレート:Cite journal
• テンプレート:Cite journal
• テンプレート:Cite book

関連項目

• 境界要素法
• 数値解析
• 有限要素法
• 有限体積法

外部リンク

• 代用電荷法（基本解近似解法）Charge Simulation Method (Method of Fundamental Solutions)
• 代用電荷法と双極子法の理論的・実験的研究
• 代用電荷法と数値等角写像に関する研究

テンプレート:偏微分方程式の数値解法