§E15.29Poisson 方程式と Newton ポテンシャル

最終更新

Poisson 方程式は、与えられた外力から未知関数を定める楕円型方程式の最も基本的な例の一つである。例えば、三次元の球内に一様に分布する外力がつくる Newton ポテンシャルは、球内では半径の二次式、球外では半径の逆数で表される。しかし、全空間の方程式を満たすポテンシャルが、領域の境界で指定された値をとるとは限らない。基本解に境界条件を反映した Green 関数は、外力から零境界値の解を定める積分核であり、弱形式による解の存在と具体的な積分表示を結びつける。本記事では、Newton ポテンシャルと Green 関数の基本的な性質および代表的な領域での積分表示について解説する。

1 Newton ポテンシャル

定義 1.1 (Newton ポテンシャル).n≥2n\geq2を整数とし、σn−1\sigma_{n-1}を単位球面Sn−1S^{n-1}の面積とする。−Δ-\Deltaの基本解

Φn(x)={∣x∣2−n(n−2)σn−1,n≥3,−12πlog⁡∣x∣,n=2\Phi_n(x)= \begin{cases} \displaystyle\frac{|x|^{2-n}}{(n-2)\sigma_{n-1}},&n\geq3,\\[4pt] \displaystyle-\frac{1}{2\pi}\log|x|,&n=2 \end{cases}

に対し、コンパクト台をもつf∈L2(Rn;C)f\in L^2(\R^n;\C)の Newton ポテンシャル (Newton potential) を

Nf(x)=∫RnΦn(x−y)f(y) dyNf(x)=\int_{\R^n}\Phi_n(x-y)f(y)\,dy

と定める。積分は絶対収束する点で定義し、それ以外では値を零とする。

定理 1.2.n≥2n\geq2を整数とする。このとき、コンパクト台をもつf∈L2(Rn;C)f\in L^2(\R^n;\C)に対して、Newton ポテンシャルの積分はほとんど至る所で絶対収束し、Nf∈Hloc2(Rn)Nf\in H^2_{\mathrm{loc}}(\R^n)および

−ΔNf=fin D′(Rn)-\Delta Nf=f\quad\text{in }\mathcal D'(\R^n)

が成り立つ。さらにf∈Cc∞(Rn;C)f\in C_c^\infty(\R^n;\C)ならばNf∈C∞(Rn)Nf\in C^\infty(\R^n)であり、∣x∣→∞|x|\to\inftyにおいて

Nf(x)={O(∣x∣2−n),n≥3,−12π(∫R2f)log⁡∣x∣+O(∣x∣−1),n=2Nf(x)= \begin{cases} O(|x|^{2-n}),&n\geq3,\\[2pt] \displaystyle-\frac{1}{2\pi}\left(\int_{\R^2}f\right)\log|x|+O(|x|^{-1}),&n=2 \end{cases}

が成り立つ。

証明. 任意の有界集合U⊂RnU\subset\R^nに対し、x∈Ux\in Uとy∈supp⁡fy\in\supp fの差をすべて含む球BMB_Mをとる。Φn1BM∈L1\Phi_n\mathbf1_{B_M}\in L^1なので、§E14.7 定理 2.2を∣Φn∣1BM|\Phi_n|\mathbf1_{B_M}と∣f∣|f|に適用すると、定義の積分はUUのほとんど至る所で絶対収束し、Nf∈L2(U)Nf\in L^2(U)となる。ffはコンパクト台の超関数でもあるから、§E15.6 定理 2.2と§E15.6 命題 1.2により−ΔNf=f-\Delta Nf=fである。§E15.21 定理 2.2を次数22、s=0s=0で用いるとNf∈Hloc2Nf\in H^2_{\mathrm{loc}}を得る。

f∈Cc∞f\in C_c^\inftyならば、畳み込みの微分をffに移すことにより∂αNf=Φn∗(∂αf)\partial^\alpha Nf=\Phi_n*(\partial^\alpha f)となり、§E15.6 命題 1.2によりNfNfは滑らかである。supp⁡f⊂BR\supp f\subset B_Rとし、∣x∣>2R|x|>2Rとする。n≥3n\geq3では∣x−y∣≥∣x∣/2|x-y|\geq|x|/2なので

∣Nf(x)∣≤2n−2(n−2)σn−1∣x∣2−n∥f∥1.|Nf(x)|\leq\frac{2^{n-2}}{(n-2)\sigma_{n-1}}|x|^{2-n}\|f\|_1.

n=2n=2では、∣y∣≤R|y|\leq Rに関して一様に

log⁡∣x−y∣−log⁡∣x∣=O(R/∣x∣)\log|x-y|-\log|x|=O(R/|x|)

である。両辺にf(y)f(y)を掛けて積分すると、二次元の展開を得る。▨

例 1.3.R3\R^3の球BR=B(0,R)B_R=B(0,R)に対してf=1BRf=\mathbf1_{B_R}とすると、r=∣x∣r=|x|の関数として

Nf(x)={R22−r26,0≤r≤R,R33r,r≥RNf(x)= \begin{cases} \displaystyle\frac{R^2}{2}-\frac{r^2}{6},&0\leq r\leq R,\\[4pt] \displaystyle\frac{R^3}{3r},&r\geq R \end{cases}

となる。実際、r,s>0r,s>0に対して球面上の極座標を用いると

14π∫S2dω∣x−sω∣=12∫−11dtr2+s2−2rst=(r+s)−∣r−s∣2rs=1max⁡{r,s}.\frac{1}{4\pi}\int_{S^2}\frac{d\omega}{|x-s\omega|} =\frac12\int_{-1}^1\frac{dt}{\sqrt{r^2+s^2-2rst}} =\frac{(r+s)-|r-s|}{2rs} =\frac1{\max\{r,s\}}.

したがってNf(x)=∫0Rs2/max⁡{r,s} dsNf(x)=\int_0^R s^2/\max\{r,s\}\,dsであり、s=rs=rで積分を分けると表示式を得る。r=0r=0では直接積分してNf(0)=∫0Rs ds=R2/2Nf(0)=\int_0^R s\,ds=R^2/2となる。

例 1.4.0≤f∈Cc∞(R2)0\leq f\in C_c^\infty(\R^2)、∫f=1\int f=1とする。このとき定理 1.2により

Nf(x)=−12πlog⁡∣x∣+O(∣x∣−1)⟶−∞.Nf(x)=-\frac{1}{2\pi}\log|x|+O(|x|^{-1})\longrightarrow-\infty.

したがって、三次元以上での無限遠の零条件は、二次元の非零総量をもつ外力に対しては成り立たない。

2 有界領域の Green 関数

定理 2.1.n≥2n\geq2を整数、Ω⊂Rn\Omega\subset\R^nを空でない有界 Lipschitz 開集合、y∈Ωy\in\Omegaとする。実関数χy∈Cc∞(Ω)\chi_y\in C_c^\infty(\Omega)をyyの近傍で11となるようにとり、

ρy=(Δχy)Φn(⋅−y)+2∇χy⋅∇Φn(⋅−y)\rho_y=(\Delta\chi_y)\Phi_n(\cdot-y) +2\nabla\chi_y\cdot\nabla\Phi_n(\cdot-y)

とおく。このとき、−Δvy=ρy-\Delta v_y=\rho_yを満たす実弱解vy∈H01(Ω)v_y\in H^1_0(\Omega)は一意に存在する。関数

K(x,y)=χy(x)Φn(x−y)+vy(x)(x≠y)K(x,y)=\chi_y(x)\Phi_n(x-y)+v_y(x)\qquad(x\ne y)

はχy\chi_yの選択に依存せず、

−ΔxK(⋅,y)=δy,K(⋅,y)−χyΦn(⋅−y)∈H01(Ω)-\Delta_xK(\cdot,y)=\delta_y,\qquad K(\cdot,y)-\chi_y\Phi_n(\cdot-y)\in H^1_0(\Omega)

を満たす。この二条件はK(⋅,y)K(\cdot,y)を超関数として一意に定める。KKをΩ\Omegaの Dirichlet Green 関数 (Dirichlet Green function) という。

関数hy=Φn(⋅−y)−K(⋅,y)h_y=\Phi_n(\cdot-y)-K(\cdot,y)はH1(Ω)H^1(\Omega)に属する調和関数であり、Ω\Omegaの内部で滑らかである。(x,y)↦hy(x)(x,y)\mapsto h_y(x)はΩ×Ω\Omega\times\Omega上で連続である。したがってKKは対角線の外で連続であり、対角線上の値を零とすれば可測である。

証明.∇χy\nabla\chi_yとΔχy\Delta\chi_yはyyの近傍で零なので、ρy∈Cc∞(Ω;R)\rho_y\in C_c^\infty(\Omega;\R)である。§E15.8 定理 5.4により、H01(Ω)H^1_0(\Omega)上の形式

a(u,v)=∫Ω∇u⋅∇v‾ dxa(u,v)=\int_\Omega\nabla u\cdot\overline{\nabla v}\,dx

は強圧的である。内積と形式は第一変数について線形とし、外力はv↦∫ρyv‾v\mapsto\int\rho_y\overline vという反線形汎関数として扱う。§E15.16 定理 2.2によりvyv_yが一意に存在する。係数と外力が実数値なので、共役をとった方程式の一意性によりvyv_yは実数値である。

基本解の等式§E15.6 定理 2.2に積の微分を適用すると

−Δ(χyΦn(⋅−y))=δy−ρy.-\Delta(\chi_y\Phi_n(\cdot-y))=\delta_y-\rho_y.

したがって定義したKKは二条件を満たす。同じ切断に対して二条件を満たす二つの超関数の差は、H01(Ω)H^1_0(\Omega)に属する調和関数である。その弱方程式に差自身を代入すると勾配が零となり、Poincaré 不等式により差は零である。二つの切断χy,χ~y\chi_y,\widetilde\chi_yの差とΦn(⋅−y)\Phi_n(\cdot-y)との積はCc∞(Ω)C_c^\infty(\Omega)に属するので、切断を替えても同じ一意性を適用することができる。

wy=(1−χy)Φn(⋅−y)w_y=(1-\chi_y)\Phi_n(\cdot-y)はΩ‾\overline\Omegaの近傍で滑らかであり、−Δwy=ρy-\Delta w_y=\rho_yである。したがってhy=wy−vy∈H1(Ω)h_y=w_y-v_y\in H^1(\Omega)は調和的であり、§E15.21 定理 2.2により内部で滑らかである。

y0∈Ωy_0\in\Omegaを固定し、y0y_0のある近傍V⋐ΩV\Subset\Omega全体で11となる切断χ\chiを一つ選ぶ。y∈Vy\in Vに対してこの切断を用いると、y↦ρyy\mapsto\rho_yは各整数k≥0k\geq0についてHk(Ω)H^k(\Omega)のノルムで連続である。実際、χ\chiの微分の台とVVを正の距離だけ離して選べば、ρy\rho_yの各微分は共通のコンパクト台上でyyに関して一様に連続となる。解写像の評価から

∥vy−vz∥H1(Ω)≤C∥ρy−ρz∥L2(Ω).\|v_y-v_z\|_{H^1(\Omega)}\leq C\|\rho_y-\rho_z\|_{L^2(\Omega)}.

内部の球U⋐U′⋐ΩU\Subset U'\Subset\Omegaにおいて、§E15.22 定理 3.1を差の方程式へ適用すると

∥vy−vz∥Hk+2(U)≤Ck(∥ρy−ρz∥Hk(U′)+∥vy−vz∥H1(U′)).\|v_y-v_z\|_{H^{k+2}(U)} \leq C_k\bigl(\|\rho_y-\rho_z\|_{H^k(U')}+\|v_y-v_z\|_{H^1(U')}\bigr).

k+2>n/2k+2>n/2とし、滑らかな球上の埋め込み§E15.10 系 4.2を用いると、y↦vyy\mapsto v_yはC(U‾)C(\overline U)のノルムで連続となる。wy(x)w_y(x)は(x,y)∈U×V(x,y)\in U\times Vの滑らかな関数なので、hy(x)=wy(x)−vy(x)h_y(x)=w_y(x)-v_y(x)は同時に連続である。xxとy0y_0は任意であったから、hhの連続性とKKの対角線外での連続性を得る。▨

命題 2.2.n≥2n\geq2、Ω⊂Rn\Omega\subset\R^nは空でない有界 Lipschitz 開集合とし、KKをその Dirichlet Green 関数とする。このとき、相異なるすべてのx,y∈Ωx,y\in\Omegaに対してK(x,y)≥0K(x,y)\geq0である。n≥3n\geq3ならば、さらに

0≤K(x,y)≤Φn(x−y)0\leq K(x,y)\leq\Phi_n(x-y)

が成り立つ。

証明.y∈Ωy\in\Omegaを固定する。hyh_yはyyの近傍で有界であり、Φn(x−y)→+∞\Phi_n(x-y)\to+\inftyであるから、十分小さいε>0\varepsilon>0についてK(x,y)>0K(x,y)>0が∣x−y∣=ε|x-y|=\varepsilon上で成り立つ。Ω∖B(y,ε)‾\Omega\setminus\overline{B(y,\varepsilon)}上のK(⋅,y)K(\cdot,y)はH1H^1に属する調和関数であり、外側の境界のトレースは零である。内側の境界のトレースは非負なので、§E15.18 定理 3.1を−K(⋅,y)-K(\cdot,y)に適用するとK(⋅,y)≥0K(\cdot,y)\geq0を得る。内部の連続性により、不等式はΩ∖B(y,ε)‾\Omega\setminus\overline{B(y,\varepsilon)}の各点で成り立つ。ε\varepsilonを任意に小さくすると、すべてのx∈Ω∖{y}x\in\Omega\setminus\{y\}についてK(x,y)≥0K(x,y)\geq0を得る。

n≥3n\geq3ならば、hyh_yの境界トレースはΦn(⋅−y)≥0\Phi_n(\cdot-y)\geq0である。§E15.18 定理 3.1を−hy-h_yに適用するとhy≥0h_y\geq0となり、K=Φn−hy≤ΦnK=\Phi_n-h_y\leq\Phi_nを得る。▨

3 Green 作用素の核表示

補題 3.1.n≥2n\geq2、Ω⊂Rn\Omega\subset\R^nは空でない有界 Lipschitz 開集合とする。形式a(u,v)=∫Ω∇u⋅∇v‾a(u,v)=\int_\Omega\nabla u\cdot\overline{\nabla v}、D(a)=H01(Ω)D(a)=H^1_0(\Omega)のL2L^2実現をADA_Dとし、G=AD−1G=A_D^{-1}とおく。KKをΩ\Omegaの Dirichlet Green 関数とする。このとき、任意の実関数f∈Cc∞(Ω)f\in C_c^\infty(\Omega)と任意のy∈Ωy\in\Omegaに対して

∫Ωf(x)K(x,y) dx=(Gf)(y)\int_\Omega f(x)K(x,y)\,dx=(Gf)(y)

が成り立つ。右辺にはGfGfの内部で滑らかな代表元を用いる。

証明.§E15.8 定理 5.4と§E15.16 定理 2.2により、各f∈L2f\in L^2に対してa(u,v)=(f,v)L2a(u,v)=(f,v)_{L^2}を満たすu∈H01u\in H^1_0が一意に存在する。§E15.16 定理 3.1によりADA_Dは自己共役であり、その逆GGは有界な自己共役作用素である。ker⁡AD={0}\ker A_D=\{0\}なので、GGは§E15.17 定理 3.1の Green 作用素に一致する。

実関数f∈Cc∞(Ω)f\in C_c^\infty(\Omega)に対してu=Gfu=Gfとおくと、§E15.21 定理 2.2によりuuは内部で滑らかである。qy=χyΦn(⋅−y)q_y=\chi_y\Phi_n(\cdot-y)とおき、supp⁡χy\supp\chi_yの近傍で11となるη∈Cc∞(Ω)\eta\in C_c^\infty(\Omega)をとる。−Δqy=δy−ρy-\Delta q_y=\delta_y-\rho_yを試験関数ηu\eta uに適用すると

∫Ωqyf=u(y)−∫Ωρyu.\int_\Omega q_y f=u(y)-\int_\Omega\rho_yu.

一方、u,vy∈H01(Ω)u,v_y\in H^1_0(\Omega)の弱方程式に互いを代入すると

∫Ωfvy=∫Ω∇u⋅∇vy=∫Ωρyu.\int_\Omega fv_y=\int_\Omega\nabla u\cdot\nabla v_y =\int_\Omega\rho_yu.

二式を加え、K(⋅,y)=qy+vyK(\cdot,y)=q_y+v_yを代入すると結論を得る。▨

補題 3.2.n≥2n\geq2、Ω⊂Rn\Omega\subset\R^nは空でない有界 Lipschitz 開集合とし、KKをその Dirichlet Green 関数とする。Ω‾⊂B(x0,R)\overline\Omega\subset B(x_0,R)とする。このとき、すべてのy∈Ωy\in\Omegaに対して

∫ΩK(x,y) dx≤R2−∣y−x0∣22n≤R22n\int_\Omega K(x,y)\,dx\leq\frac{R^2-|y-x_0|^2}{2n}\leq\frac{R^2}{2n}

が成り立つ。

証明.0≤ψj≤10\leq\psi_j\leq1、ψj∈Cc∞(Ω)\psi_j\in C_c^\infty(\Omega)を、各点で十分大きいjjに対してψj=1\psi_j=1となるように選び、fm=1−∏j=1m(1−ψj)f_m=1-\prod_{j=1}^m(1-\psi_j)とおく。このときfm∈Cc∞(Ω)f_m\in C_c^\infty(\Omega)、0≤fm↑10\leq f_m\uparrow1である。B(x)=(R2−∣x−x0∣2)/(2n)B(x)=(R^2-|x-x_0|^2)/(2n)とおくと−ΔB=1-\Delta B=1であり、B>0B>0がΩ‾\overline\Omega上で成り立つ。§E15.18 定理 3.1を−Gfm-Gf_mとGfm−BGf_m-Bに適用すると

0≤(Gfm)(y)≤B(y)(y∈Ω)0\leq(Gf_m)(y)\leq B(y)\qquad(y\in\Omega)

となる。各関数の内部での連続性により、不等式は各点で成立する。補題 3.1とK≥0K\geq0より、単調収束定理を適用して

∫ΩK(x,y) dx=lim⁡m→∞∫Ωfm(x)K(x,y) dx=lim⁡m→∞(Gfm)(y)≤B(y)\int_\Omega K(x,y)\,dx =\lim_{m\to\infty}\int_\Omega f_m(x)K(x,y)\,dx =\lim_{m\to\infty}(Gf_m)(y)\leq B(y)

を得る。▨

定理 3.3.n≥2n\geq2、Ω⊂Rn\Omega\subset\R^nは空でない有界 Lipschitz 開集合とし、KKをその Dirichlet Green 関数、G=AD−1G=A_D^{-1}を Dirichlet Laplace 作用素の Green 作用素とする。このとき、相異なるすべてのx,y∈Ωx,y\in\Omegaに対して

K(x,y)=K(y,x)K(x,y)=K(y,x)

が成り立つ。また、ある定数CΩ<∞C_\Omega<\inftyに対して

sup⁡y∈Ω∫ΩK(x,y) dx≤CΩ,sup⁡x∈Ω∫ΩK(x,y) dy≤CΩ\sup_{y\in\Omega}\int_\Omega K(x,y)\,dx\leq C_\Omega, \qquad \sup_{x\in\Omega}\int_\Omega K(x,y)\,dy\leq C_\Omega

である。任意のf∈L2(Ω;C)f\in L^2(\Omega;\C)に対して、積分はほとんど至る所で絶対収束し、

(Gf)(x)=∫ΩK(x,y)f(y) dy(Gf)(x)=\int_\Omega K(x,y)f(y)\,dy

がほとんど至る所で成り立つ。特に∫ΩK(x,y) dy=(G1)(x)\int_\Omega K(x,y)\,dy=(G1)(x)がほとんど至る所で成り立つ。

証明.補題 3.2とΩ\Omegaの有界性からK∈L1(Ω×Ω)K\in L^1(\Omega\times\Omega)である。任意の実関数f,g∈Cc∞(Ω)f,g\in C_c^\infty(\Omega)に対して、K(x,y)f(x)g(y)K(x,y)f(x)g(y)は絶対可積分なので、§E10.19 定理 3.2と補題 3.1により

∬Ω×ΩK(x,y)f(x)g(y) dx dy=∫Ω(Gf)g=∫Ωf(Gg)=∬Ω×ΩK(y,x)f(x)g(y) dx dy.\iint_{\Omega\times\Omega}K(x,y)f(x)g(y)\,dx\,dy =\int_\Omega(Gf)g =\int_\Omega f(Gg) =\iint_{\Omega\times\Omega}K(y,x)f(x)g(y)\,dx\,dy.

中央の等式はGGの自己共役性による。

L(x,y)=K(x,y)−K(y,x)L(x,y)=K(x,y)-K(y,x)とおく。任意の積試験関数とのLLの積分は零である。非負のρ∈Cc∞(B(0,1/2))\rho\in C_c^\infty(B(0,1/2))を∫ρ=1\int\rho=1となるようにとり、ρε(z)=ε−nρ(z/ε)\rho_\varepsilon(z)=\varepsilon^{-n}\rho(z/\varepsilon)とおく。R2n\R^{2n}上で積型の平滑化核ρε(x)ρε(y)\rho_\varepsilon(x)\rho_\varepsilon(y)を用いると、Ω×Ω\Omega\times\Omegaの各内部点で、十分小さいε\varepsilonに対するLLの平滑化は零となる。§E15.2 定理 4.5により平滑化は超関数としてLLに収束するので、L=0L=0が超関数として、したがってほとんど至る所で成り立つ。対角線外での連続性により、K(x,y)=K(y,x)K(x,y)=K(y,x)はx≠yx\ne yのすべてで成り立つ。

対称性と補題 3.2から両側の積分上界を得る。§E10.22 定理 8.1により

TKf(x)=∫ΩK(x,y)f(y) dyT_Kf(x)=\int_\Omega K(x,y)f(y)\,dy

はL2(Ω;C)L^2(\Omega;\C)上の有界作用素を定める。補題 3.1と対称性によりTKf=GfT_Kf=Gfが実数値のCc∞C_c^\infty上で成り立つ。複素線形性とCc∞(Ω)C_c^\infty(\Omega)のL2L^2における稠密性から、すべてのf∈L2f\in L^2に対して等式が成り立つ。Ω\Omegaは有界なので1∈L2(Ω)1\in L^2(\Omega)であり、f=1f=1を代入すると最後の等式を得る。▨

系 3.4.n≥2n\geq2、Ω⊂Rn\Omega\subset\R^nは空でない有界 Lipschitz 開集合とし、KKをその Dirichlet Green 関数とする。このとき、K∈L2(Ω×Ω)K\in L^2(\Omega\times\Omega)であるための必要十分条件はn≤3n\leq3である。n≤3n\leq3ならば、Dirichlet Laplace 作用素の完全正規直交固有系(ej,λj)(e_j,\lambda_j)に対して

K(x,y)=∑j=1∞λj−1ej(x)ej(y)‾in L2(Ω×Ω)K(x,y)=\sum_{j=1}^\infty\lambda_j^{-1}e_j(x)\overline{e_j(y)} \quad\text{in }L^2(\Omega\times\Omega)

が成り立つ。

証明.n≤3n\leq3ならば、§E15.19 系 4.2と§E15.19 定理 4.1により、表示された級数はL2L^2核K2K_2に収束し、GGを表す。K−K2∈L1(Ω×Ω)K-K_2\in L^1(\Omega\times\Omega)であり、二つの核が同じ作用素を表すことから、K−K2K-K_2と任意の積試験関数との積分は零である。定理 3.3の証明と同じ積型の平滑化によりK=K2K=K_2がほとんど至る所で成り立つ。

n≥4n\geq4でK∈L2K\in L^2と仮定すると、§E12.10 定理 5.1によりG=TKG=T_Kは Hilbert–Schmidt 作用素となる。これは§E15.19 系 4.2に反する。▨

4 球と半空間の Green 関数

命題 4.1.n≥2n\geq2、R>0R>0とし、BR=B(0,R)⊂RnB_R=B(0,R)\subset\R^nの Dirichlet Green 関数をGRG_Rと書く。0≠y∈BR0\ne y\in B_Rに対してy∗=R2y/∣y∣2y^*=R^2y/|y|^2とおく。このとき、

GR(x,y)={Φn(x−y)−(R∣y∣)n−2Φn(x−y∗),n≥3,Φ2(x−y)−Φ2(∣y∣R(x−y∗)),n=2G_R(x,y)= \begin{cases} \displaystyle\Phi_n(x-y)-\left(\frac R{|y|}\right)^{n-2}\Phi_n(x-y^*),&n\geq3,\\[4pt] \displaystyle\Phi_2(x-y)-\Phi_2\left(\frac{|y|}{R}(x-y^*)\right),&n=2 \end{cases}

である。y=0y=0ではGR(x,0)=Φn(x)−Φn(R)G_R(x,0)=\Phi_n(x)-\Phi_n(R)である。ここでΦn(R)\Phi_n(R)は半径RRでの基本解の値を表す。GRG_Rは対角線外で対称であり、−ΔxGR(⋅,y)=δy-\Delta_xG_R(\cdot,y)=\delta_yとGR(⋅,y)∣∂BR=0G_R(\cdot,y)|_{\partial B_R}=0を満たす。

証明.∣y∗∣=R2/∣y∣>R|y^*|=R^2/|y|>Rなので、表示式の補正項は球の閉包の近傍で滑らかである。動径関数の公式ΔF(∣x∣)=F′′(∣x∣)+(n−1)F′(∣x∣)/∣x∣\Delta F(|x|)=F''(|x|)+(n-1)F'(|x|)/|x|をF(r)=r2−nF(r)=r^{2-n}またはF(r)=log⁡rF(r)=\log rに用いると、補正項は調和的である。∣x∣=R|x|=Rのとき

∣x−y∗∣2=R2∣y∣2∣x−y∣2|x-y^*|^2=\frac{R^2}{|y|^2}|x-y|^2

なので、両次元の表示式は境界上で零となる。y=0y=0の補正項は定数であり、同じ境界条件を満たす。

任意の切断χy\chi_yに対して、表示式からχyΦn(⋅−y)\chi_y\Phi_n(\cdot-y)を引いた関数は閉球の近傍で滑らかであり、球面上で零となる。§E15.12 定理 3.3により、この差はH01(BR)H^1_0(B_R)に属する。定理 2.1の一意性により、表示式はBRB_Rの Green 関数である。対称性は定理 3.3から従う。▨

命題 4.2.n≥2n\geq2とし、H={x∈Rn:xn>0}H=\{x\in\R^n:x_n>0\}とする。y=(y′,yn)∈Hy=(y',y_n)\in Hとy∗=(y′,−yn)y^*=(y',-y_n)に対して

GH(x,y)=Φn(x−y)−Φn(x−y∗)G_H(x,y)=\Phi_n(x-y)-\Phi_n(x-y^*)

とおく。このときGH(x,y)=GH(y,x)G_H(x,y)=G_H(y,x)、GH(x,y)>0G_H(x,y)>0がx≠yx\ne yで成り立ち、

−ΔxGH(⋅,y)=δyin D′(H),GH(⋅,y)∣∂H=0-\Delta_xG_H(\cdot,y)=\delta_y\quad\text{in }\mathcal D'(H), \qquad G_H(\cdot,y)|_{\partial H}=0

を満たす。

証明.y∗∉H‾y^*\notin\overline Hなので、Φn(x−y∗)\Phi_n(x-y^*)はHHで滑らかかつ調和的である。基本解の等式から超関数の意味での方程式が従う。xn=0x_n=0ならば∣x−y∣=∣x−y∗∣|x-y|=|x-y^*|なので境界値は零である。x,y∈Hx,y\in Hでは

∣x−y∗∣2−∣x−y∣2=4xnyn>0.|x-y^*|^2-|x-y|^2=4x_ny_n>0.

Φn\Phi_nは半径に関して狭義単調減少であるからGH>0G_H>0となる。∣x−y∗∣2=∣x′−y′∣2+(xn+yn)2|x-y^*|^2=|x'-y'|^2+(x_n+y_n)^2はx,yx,yについて対称なので、GHG_Hも対称である。▨

命題 4.3.n≥2n\geq2、R>0R>0とし、GRG_RをBRB_Rの Dirichlet Green 関数とする。u∈C2(BR;C)u\in C^2(B_R;\C)の二階以下の偏導関数がすべて閉球まで連続に延長されるとする。このとき、任意のx∈BRx\in B_Rについて

u(x)=∫BRGR(x,y)(−Δu)(y) dy−∫∂BR∂νyGR(x,y)u(y) dσ(y)u(x)=\int_{B_R}G_R(x,y)(-\Delta u)(y)\,dy -\int_{\partial B_R}\partial_{\nu_y}G_R(x,y)u(y)\,d\sigma(y)

が成り立つ。ここでνy=y/R\nu_y=y/Rは球面の外向き単位法線である。

証明.x∈BRx\in B_Rを固定し、∣x∣+ε<s<R|x|+\varepsilon<s<Rに対してDs,ε=Bs∖B(x,ε)‾D_{s,\varepsilon}=B_s\setminus\overline{B(x,\varepsilon)}とおく。この領域においてF=u∇yGR(x,⋅)−GR(x,⋅)∇uF=u\nabla_yG_R(x,\cdot)-G_R(x,\cdot)\nabla uとおく。uuとGR(x,⋅)G_R(x,\cdot)はDs,εD_{s,\varepsilon}の閉包の近傍でそれぞれC2C^2級、滑らかであり、GR(x,⋅)G_R(x,\cdot)は調和的なので、§E15.12 定理 1.3により

∫Ds,εGR(x,y)(−Δu)(y) dy=∫∂Ds,ε(u∂νGR−GR∂νu) dσ.\int_{D_{s,\varepsilon}}G_R(x,y)(-\Delta u)(y)\,dy =\int_{\partial D_{s,\varepsilon}} \bigl(u\partial_\nu G_R-G_R\partial_\nu u\bigr)\,d\sigma.

s↑Rs\uparrow Rとすると、u,∇uu,\nabla uの境界までの連続性とGR∣∂BR=0G_R|_{\partial B_R}=0により、外側の境界積分は∫∂BRu∂νGR dσ\int_{\partial B_R}u\partial_\nu G_R\,d\sigmaに収束する。内側の球面の法線は−(y−x)/ε-(y-x)/\varepsilonであり、基本解の動径微分はΦn′(r)=−1/(σn−1rn−1)\Phi_n'(r)=-1/(\sigma_{n-1}r^{n-1})である。GR(x,y)−Φn(y−x)G_R(x,y)-\Phi_n(y-x)はxxの近傍で滑らかなので

∫∂B(x,ε)u∂νGR dσ=1σn−1εn−1∫∂B(x,ε)u dσ+o(1)⟶u(x).\int_{\partial B(x,\varepsilon)}u\partial_\nu G_R\,d\sigma =\frac1{\sigma_{n-1}\varepsilon^{n-1}} \int_{\partial B(x,\varepsilon)}u\,d\sigma+o(1)\longrightarrow u(x).

他方、内側の球面におけるGR∂νuG_R\partial_\nu uの積分の絶対値は、n≥3n\geq3ではO(ε)O(\varepsilon)、n=2n=2ではO(ε∣log⁡ε∣)O(\varepsilon|\log\varepsilon|)となり、零へ収束する。GRG_Rの特異性は局所可積分なので、体積積分にもε→0\varepsilon\to0の極限をとることができる。極限の等式をu(x)u(x)について解くと表示式を得る。▨

例 4.4.n≥2n\geq2、R>0R>0とする。このとき、BRB_Rの Green 関数の補正項hy(x)=Φn(x−y)−GR(x,y)h_y(x)=\Phi_n(x-y)-G_R(x,y)は、定理 2.1の連続性により、BR/2‾×BR/2‾\overline{B_{R/2}}\times\overline{B_{R/2}}上で有界である。n≥3n\geq3ならば、ある0<δ<R/40<\delta<R/4が存在して、y∈BR/4y\in B_{R/4}、0<∣x−y∣<δ0<|x-y|<\deltaに対して

GR(x,y)≥12Φn(x−y).G_R(x,y)\geq\frac12\Phi_n(x-y).

この領域での二乗積分は定数倍の

∫0δr2(2−n)rn−1 dr=∫0δr3−n dr\int_0^\delta r^{2(2-n)}r^{n-1}\,dr =\int_0^\delta r^{3-n}\,dr

を下回らないので、n≥4n\geq4では発散する。n=3n=3では0≤GR≤Φ30\leq G_R\leq\Phi_3と同じ動径積分から二乗可積分性が従う。n=2n=2の特異項は対数であり、∫0δr∣log⁡r∣2 dr<∞\int_0^\delta r|\log r|^2\,dr<\inftyである。

5 Poisson 積分による調和拡張

定理 5.1.n≥2n\geq2、R>0R>0とし、x∈BRx\in B_R、ζ∈∂BR\zeta\in\partial B_Rに対して

PR(x,ζ)=−∂νζGR(x,ζ)=R2−∣x∣2σn−1R∣x−ζ∣nP_R(x,\zeta)=-\partial_{\nu_\zeta}G_R(x,\zeta) =\frac{R^2-|x|^2}{\sigma_{n-1}R|x-\zeta|^n}

とおく。このとき、Poisson 核 (Poisson kernel)PRP_Rは正値であり、

∫∂BRPR(x,ζ) dσ(ζ)=1\int_{\partial B_R}P_R(x,\zeta)\,d\sigma(\zeta)=1

を満たす。任意のg∈C(∂BR;C)g\in C(\partial B_R;\C)に対して

u(x)=∫∂BRPR(x,ζ)g(ζ) dσ(ζ)u(x)=\int_{\partial B_R}P_R(x,\zeta)g(\zeta)\,d\sigma(\zeta)

はC∞(BR)∩C(BR‾)C^\infty(B_R)\cap C(\overline{B_R})に属する唯一の調和拡張である。

証明.命題 4.1の表示を、対称性を用いて第二変数について書く。x≠0x\ne0と∣ζ∣=R|\zeta|=Rに対して∣ζ−x∗∣=(R/∣x∣)∣ζ−x∣|\zeta-x^*|=(R/|x|)|\zeta-x|であり、全次元で

∇Φn(z)=−zσn−1∣z∣n.\nabla\Phi_n(z)=-\frac{z}{\sigma_{n-1}|z|^n}.

球面の法線ζ/R\zeta/Rに沿って微分すると

∂νζGR(x,ζ)=−(R2−ζ⋅x)−(∣x∣2/R2)(R2−ζ⋅x∗)σn−1R∣x−ζ∣n=−R2−∣x∣2σn−1R∣x−ζ∣n.\partial_{\nu_\zeta}G_R(x,\zeta) =-\frac{(R^2-\zeta\cdot x)-(|x|^2/R^2)(R^2-\zeta\cdot x^*)} {\sigma_{n-1}R|x-\zeta|^n} =-\frac{R^2-|x|^2}{\sigma_{n-1}R|x-\zeta|^n}.

二次元では対数の定数項の微分が零となるので、同じ式が成り立つ。x=0x=0ではGR(0,y)=Φn(y)−Φn(R)G_R(0,y)=\Phi_n(y)-\Phi_n(R)を微分すると同じ表示を得る。正値性は表示式から従い、質量が11であることは命題 4.3にu=1u=1を代入して従う。

xxがBRB_Rのコンパクト部分集合を動くとき、PRP_Rのxxに関する各微分はζ\zetaに関して一様に有界である。またGR(x,y)G_R(x,y)は、yyが球面の近傍にあるとき、内部のxxに関して調和的である。したがってΔxPR(x,ζ)=0\Delta_xP_R(x,\zeta)=0であり、積分の下で微分することによりu∈C∞(BR)u\in C^\infty(B_R)とΔu=0\Delta u=0を得る。

ζ0∈∂BR\zeta_0\in\partial B_Rとϵ>0\epsilon>0をとる。ggの連続性により、あるδ>0\delta>0に対して∣ζ−ζ0∣<δ|\zeta-\zeta_0|<\deltaならば∣g(ζ)−g(ζ0)∣<ϵ|g(\zeta)-g(\zeta_0)|<\epsilonである。核の正値性と質量の等式から

∣u(x)−g(ζ0)∣≤ϵ+2∥g∥∞∫∣ζ−ζ0∣≥δPR(x,ζ) dσ(ζ).|u(x)-g(\zeta_0)| \leq\epsilon+2\|g\|_\infty \int_{|\zeta-\zeta_0|\geq\delta}P_R(x,\zeta)\,d\sigma(\zeta).

x→ζ0x\to\zeta_0のとき、最後の積分では分母が正の定数で下から抑えられ、分子R2−∣x∣2R^2-|x|^2は零へ収束する。よってu(x)→g(ζ0)u(x)\to g(\zeta_0)である。一意性は§E15.18 系 1.2を二つの解の差の実部と虚部に適用して従う。▨

定理 5.2.n≥2n\geq2、H={xn>0}H=\{x_n>0\}とする。x=(x′,xn)∈Hx=(x',x_n)\in H、ζ∈Rn−1\zeta\in\R^{n-1}に対して

PH(x,ζ)=2xnσn−1(∣x′−ζ∣2+xn2)n/2P_H(x,\zeta) =\frac{2x_n}{\sigma_{n-1}(|x'-\zeta|^2+x_n^2)^{n/2}}

とおく。このとき、PHP_Hは正値で、∫Rn−1PH(x,ζ) dζ=1\int_{\R^{n-1}}P_H(x,\zeta)\,d\zeta=1を満たす。任意のg∈Cc(Rn−1;C)g\in C_c(\R^{n-1};\C)に対して

u(x)=∫Rn−1PH(x,ζ)g(ζ) dζu(x)=\int_{\R^{n-1}}P_H(x,\zeta)g(\zeta)\,d\zeta

はC∞(H)∩C0(H‾)C^\infty(H)\cap C_0(\overline H)に属する唯一の調和拡張である。ここでC0(H‾)C_0(\overline H)は、閉半空間上で連続であり、∣x∣→∞|x|\to\inftyで零へ収束する関数の空間である。

証明. 半空間の外向き法線は−en-e_nなので、命題 4.2の表示を第二変数で微分すると

−∂νyGH(x,y)∣y=(ζ,0)=∂ynGH(x,y)∣y=(ζ,0)=PH(x,ζ).-\partial_{\nu_y}G_H(x,y)\big|_{y=(\zeta,0)} =\partial_{y_n}G_H(x,y)\big|_{y=(\zeta,0)} =P_H(x,\zeta).

PHP_Hはxxに関して滑らかかつ調和的であり、表示式から正値である。t=xnt=x_n、η=(ζ−x′)/t\eta=(\zeta-x')/tと置換すると

∫Rn−1PH(x,ζ) dζ=2σn−1∫Rn−1(1+∣η∣2)−n/2 dη=1.\int_{\R^{n-1}}P_H(x,\zeta)\,d\zeta =\frac2{\sigma_{n-1}}\int_{\R^{n-1}}(1+|\eta|^2)^{-n/2}\,d\eta=1.

最後の積分は、極座標とr=tan⁡θr=\tan\thetaによって

σn−2∫0π/2sin⁡n−2θ dθ=σn−12\sigma_{n-2}\int_0^{\pi/2}\sin^{n-2}\theta\,d\theta =\frac{\sigma_{n-1}}2

となる。n=2n=2ではσ0=2\sigma_0=2とする。

ggは有界かつコンパクト台をもつので、内部のコンパクト集合上で積分を任意回微分することができ、uuは滑らかかつ調和的である。pt(z)=2t/(σn−1(∣z∣2+t2)n/2)p_t(z)=2t/(\sigma_{n-1}(|z|^2+t^2)^{n/2})とおく。任意のδ>0\delta>0について、極座標積分から

∫∣z∣≥δpt(z) dz≤Cntδ.\int_{|z|\geq\delta}p_t(z)\,dz\leq C_n\frac{t}{\delta}.

ggの一様連続性により、任意のϵ>0\epsilon>0に対してδ>0\delta>0を選び、∣z∣<δ|z|<\deltaならばすべてのx′x'で∣g(x′−z)−g(x′)∣<ϵ|g(x'-z)-g(x')|<\epsilonとすることができる。したがって

sup⁡x′∈Rn−1∣u(x′,t)−g(x′)∣≤ϵ+2Cn∥g∥∞tδ.\sup_{x'\in\R^{n-1}}|u(x',t)-g(x')| \leq\epsilon+2C_n\|g\|_\infty\frac{t}{\delta}.

これにより閉半空間への連続拡張を得る。

supp⁡g⊂{∣ζ∣≤M}\supp g\subset\{|\zeta|\leq M\}とし、∣x∣>2M|x|>2Mとする。∣x−(ζ,0)∣≥∣x∣/2|x-(\zeta,0)|\geq|x|/2かつxn≤∣x∣x_n\leq|x|なので

∣u(x)∣≤Cn∥g∥1∣x∣1−n⟶0.|u(x)|\leq C_n\|g\|_1|x|^{1-n}\longrightarrow0.

境界上でもggは無限遠で零であるからu∈C0(H‾)u\in C_0(\overline H)となる。

二つの拡張の差の実部をwwとする。任意のϵ>0\epsilon>0に対して、十分大きいRRでは∣w∣≤ϵ|w|\leq\epsilonがH‾∩{∣x∣≥R}\overline H\cap\{|x|\geq R\}上で成り立つ。半球H∩BRH\cap B_Rの平らな境界でw=0w=0、球面部分で∣w∣≤ϵ|w|\leq\epsilonなので、§E15.18 系 1.2により∣w∣≤ϵ|w|\leq\epsilonが半球全体で成り立つ。RRを任意の固定点を含むように選び、ϵ↓0\epsilon\downarrow0とするとw=0w=0を得る。虚部にも同じ議論を適用すると一意性が従う。▨

例 5.3.H={xn>0}H=\{x_n>0\}上でu(x)=xnu(x)=x_nはΔu=0\Delta u=0とu∣∂H=0u|_{\partial H}=0を満たす非零関数である。したがって、定理 5.2から無限遠で零となる条件を外すと、零境界値に対しても一意性が成り立たない。

例 5.4.n≥2n\geq2、R>0R>0、1≤i≤n1\leq i\leq nとする。このとき、球面上の関数g(ζ)=ζig(\zeta)=\zeta_iの調和拡張はu(x)=xiu(x)=x_iである。実際、xix_iは閉球上で連続かつ内部で調和的であり、境界値がggに等しいので、定理 5.1の一意性により

∫∂BRR2−∣x∣2σn−1R∣x−ζ∣nζi dσ(ζ)=xi\int_{\partial B_R}\frac{R^2-|x|^2}{\sigma_{n-1}R|x-\zeta|^n}\zeta_i\,d\sigma(\zeta)=x_i

が成り立つ。

前提記事

15 本の記事・単元を表示