【特別寄稿】 コーシーの積分公式の周辺-1 [複素解析]
【特別寄稿】 コーシーの積分公式の周辺-1
コーシーの積分定理と積分公式から、以下を導けます。複素平面をz=x+iyで表すとして、
(1) z=0を外点とする任意の領域Rで、
ただし∫dzはRの境界C上の複素積分。
(2) z=0を内点とする任意の領域Rで、
(3) z=0で形式的に、1/z=∞。
これを見た時「デルタ関数の積分だ!」と思ったものでした(係数2πi)。というのは、複素関数を実と虚部の成分関数でf(z)=u(x,y)+iv(x,y)と書けば、複素(線)積分はガウスの発散定理で、実質は実の領域積分に直せるからです。
数値計算を(間違って?(^^))仕事にし、物理数学と渋々直面した人は、必ず1回はグリーン関数法を経験します。グリーン関数法は、グリーンさんがグリーンの公式を使って線形偏微分方程式の便利解法として開発したものです。
領域Rで定義される線形偏微分方程式を表す線形作用素をLとし、Lf(x,y)=g(x,y)はRの境界Cで境界条件Bを満たすとします。ここでf(x,y)は未知関数,g(x,y)は既知関数です。C上でBを満たし、Rで、
を満たすG(x,y,ξ,η)がグリーン関数です。δ(ξ,η)は(ξ,η)に特異点を持つデルタ関数。Gが見つかれば、線形作用素Lには解の重ね合わせが効くので、
でf(x,y)を計算できます。∫dξdηはRで行います。
ところが境界条件Bを満たすグリーン関数をみつける計算の方が、Lf=gを解くよりよっぽど難しいなんて事態は、よく起こります。そこで出来るだけ簡単なグリーン関数を使用して、境界条件のずれは数値的に合わせようとグリーン関数法を現代っぽくアレンジしたのが、境界要素法です。
Lとして一番普及してるのは、非圧縮性渦無し完全流体と静電場の支配方程式に現れるラプラシアン
です。Δに対するグリーン関数は、
と書けます。とりあえずδは原点に特異点を持つとしますが、任意の位置に特異点を持つ場合は、上式の解を平行移動するだけです。境界要素法では出来るだけ簡単なグリーン関数を使いたいので、普通はφを等方的とし境界条件は付けません。この条件を付けたグリーン関数を、ラプラス方程式の基本解と呼びます。
等方的なφとすれば極座標に移った方が便利なので、Δをで書き変えます。
※等方的なのでの項は0になる。
δ(0)はr=0以外では0です。従って実質的には、r≠0で、
を解けば良い事になります。これは線形常微分方程式です。解の公式を適用し、
を得ます。AとBは積分定数。非圧縮性渦無し完全流体でも静電場でもφはポテンシャルですので、B=0に選べます。Aの決定にはデルタ関数の性質を使います。
∫dxdyは、原点を内点として含む任意の領域Rで行います。再びδの性質から、r≠0ではδ(0)=0なので、Rとして半径εの円が可能です。∫Δφdxdyはそのままでは積分できないのでガウスの発散定理を使い、円周上の線積分ε∫dθに直します。θ=0~2πです。そうするとφ(r)の具体的形から、
が得られて、2πA=1よりA=1/(2π)になります。
すなわちラプラス方程式の基本解は、
です。δが(ξ,η)に特異点を持つ場合は、とします。
(1)~(3)はデルタ関数の実用的な定義と言えますが、(1)~(3)のかわりに、基本解φがラプラス方程式、
を満たす事を、デルタ関数の定義にも出来るはずです。
(執筆 ddt³)
以上は、ネムネコ・ファミリーのddt²さんの「ねこ騙し数学」への投稿記事です。
見やすいように若干数式に手を入れました。
なお、δ関数とは
という特殊関数のこと。
物理学者のディラックが提案した関数。
ラプラシアンの極座標表示
になることは、
第11回 合成関数の微分法2
http://nekodamashi-math.blog.so-net.ne.jp/2017-06-30-3
の問2の(2)で、
φがのみの関数、すなわち、φ(r)のとき、ラプラシアンが
がなること、そして、微分方程式
この一般解が
であることは、
第12回 合成関数の(偏)微分の問題
http://nekodamashi-math.blog.so-net.ne.jp/2017-07-01-2
の問題1で示してあるケロよ。
一般解を求めるのは、
から簡単に出てくる。
渦なし条件とは、
のことで、φが全微分であること、すなわち
であることをあらわす。
ベクトルでかけば
ここで、
である。
コメント 0