§E15.23境界正則性

最終更新

楕円型方程式を弱形式で扱うと、境界条件を含む問題を一階の Sobolev 空間の上に定めることができる。しかし、弱形式から定まる作用素の定義域は、その定義だけでは二階微分の可積分性を明示しない。境界正則性は、領域の境界と境界条件を考慮して弱解の微分可能性を明らかにする理論であり、変分的な定式化を具体的な微分作用素へ結び付ける基本的な役割を担う。たとえば、零 Dirichlet 条件を課した Poisson 方程式では、領域の境界の形状が弱解の二階微分の可積分性に関わる。本記事では、純粋な境界条件のもとでの Sobolev 正則性と、Laplace 作用素およびその冪の定義域について解説する。

1 平坦な境界での差分評価

補題 1.1.n≥2n\geq2、0<r<R0<r<Rとし、QR=(−R,R)n−1×(0,R)Q_R=(-R,R)^{n-1}\times(0,R)、ΓR=(−R,R)n−1×{0}\Gamma_R=(-R,R)^{n-1}\times\{0\}とおく。A∈W1,∞(QR;Cn×n)A\in W^{1,\infty}(Q_R;\C^{n\times n})が、あるθ>0\theta>0に対して

ℜ(A(y)ξ⋅ξ‾)≥θ∣ξ∣2ほとんどすべての y∈QR,ξ∈Cn\Re(A(y)\xi\cdot\overline\xi)\geq\theta|\xi|^2 \quad\text{ほとんどすべての }y\in Q_R,\quad\xi\in\C^n

を満たすとする。v∈H1(QR)v\in H^1(Q_R)、F0∈L2(QR)F_0\in L^2(Q_R)、F∈H1(QR;Cn)F\in H^1(Q_R;\C^n)とする。さらに、次のいずれかが成り立つとする。

  1. Tv=0Tv=0がΓR\Gamma_R上で成り立ち、ΓR\Gamma_R上でトレースが零である試験関数ψ\psiに対して ∫QRA∇v⋅∇ψ‾=∫QR(F0ψ‾+F⋅∇ψ‾)\int_{Q_R}A\nabla v\cdot\overline{\nabla\psi} =\int_{Q_R}\left(F_0\overline\psi+F\cdot\overline{\nabla\psi}\right) が成り立つ。
  2. β∈W1,∞(ΓR)\beta\in W^{1,\infty}(\Gamma_R)、r0∈H1/2(ΓR)r_0\in H^{1/2}(\Gamma_R)が存在して、試験関数ψ\psiに対して ∫QRA∇v⋅∇ψ‾+∫ΓRβTv Tψ‾=∫QR(F0ψ‾+F⋅∇ψ‾)+∫ΓRr0Tψ‾\int_{Q_R}A\nabla v\cdot\overline{\nabla\psi} +\int_{\Gamma_R}\beta Tv\,\overline{T\psi} =\int_{Q_R}\left(F_0\overline\psi+F\cdot\overline{\nabla\psi}\right) +\int_{\Gamma_R}r_0\overline{T\psi} が成り立つ。

試験関数はH1(QR)H^1(Q_R)に属し、その台はQRQ_Rの上面と側面から正の距離をもつものとする。H1/2(ΓR)H^{1/2}(\Gamma_R)は平面上の差分ノルムで定める。このときv∈H2(Qr)v\in H^2(Q_r)であり、

∥v∥H2(Qr)≤C(∥F0∥L2(QR)+∥F∥H1(QR)+∥v∥H1(QR)+∥r0∥H1/2(ΓR))\|v\|_{H^2(Q_r)} \leq C\left(\|F_0\|_{L^2(Q_R)}+\|F\|_{H^1(Q_R)} +\|v\|_{H^1(Q_R)}+\|r_0\|_{H^{1/2}(\Gamma_R)}\right)

が成り立つ。Dirichlet 条件の場合にはr0r_0の項を除く。定数CCはn,r,R,θn,r,R,\thetaとA,βA,\betaの上記ノルムにのみ依存する。

証明.r<ρ<Rr<\rho<Rを取り、r0r_0をΓρ\Gamma_\rho上で変えずに、ΓR\Gamma_Rの縁の近傍で零となるように切断する。切断後の関数を∂QR\partial Q_Rの残りの部分で零とすると、境界のH1/2H^{1/2}ノルムはC∥r0∥H1/2(ΓR)C\|r_0\|_{H^{1/2}(\Gamma_R)}以下である。実際、切断関数の差を積の差に分けると、切断の微分を含む項の核は対角集合の近傍で可積分となり、異なる面の間の積分では台と縁の距離が正である。§E15.14 定理 2.1によって、Γρ\Gamma_\rho上でTW=−r0enTW=-r_0e_nとなるW∈H1(QR;Cn)W\in H^1(Q_R;\C^n)を

∥W∥H1(QR)≤C∥r0∥H1/2(ΓR)\|W\|_{H^1(Q_R)}\leq C\|r_0\|_{H^{1/2}(\Gamma_R)}

を満たすように取ることができる。H1H^1ベクトル場の弱い法線流束は(TW)⋅ν(TW)\cdot\nuである。これは滑らかな近似に§E15.14 系 4.3を適用し、H1H^1トレースの連続性で極限を取ることから従う。したがってΓρ\Gamma_\rhoに局在する試験関数に対して

∫Γρr0Tψ‾=∫QR(div⁡W ψ‾+W⋅∇ψ‾)\int_{\Gamma_\rho}r_0\overline{T\psi} =\int_{Q_R}\left(\operatorname{div}W\,\overline\psi+W\cdot\overline{\nabla\psi}\right)

である。F0F_0をF0+div⁡WF_0+\operatorname{div}W、FFをF+WF+Wに置き換えると、r0=0r_0=0の場合に帰着する。Dirichlet 条件の場合にはこの置換を行わない。

QrQ_r上で11、上面と側面の近傍で00となる実数値切断関数η\etaを取り、その台をQρ∪ΓρQ_\rho\cup\Gamma_\rhoに含める。1≤ℓ<n1\leq\ell<nと小さいh≠0h\ne0に対して

δhw(y)=w(y+heℓ)−w(y)h,ψ=−δ−h(η2δhv)\delta_h w(y)=\frac{w(y+he_\ell)-w(y)}h, \qquad \psi=-\delta_{-h}(\eta^2\delta_hv)

とおく。接方向の平行移動はΓR\Gamma_Rを保つ。Dirichlet 条件の場合にはTδhv=δhTv=0T\delta_hv=\delta_hTv=0であるから、ψ\psiは許された試験関数である。Robin 条件の場合には零トレースを課す必要がない。

平行移動による変数変換と積の差の公式により、次の二式が成り立つ。

∫p δ−hq‾=−∫δhp q‾,δh(A∇v)=A(⋅+heℓ)∇δhv+(δhA)∇v\int p\,\overline{\delta_{-h}q}=-\int\delta_hp\,\overline q, \qquad \delta_h(A\nabla v)=A(\cdot+he_\ell)\nabla\delta_hv+(\delta_hA)\nabla v

X=∥η∇δhv∥2X=\|\eta\nabla\delta_hv\|_2、M=∥v∥H1(QR)M=\|v\|_{H^1(Q_R)}とおく。線分上の積分表示

δhv(y)=∫01∂ℓv(y+theℓ) dt\delta_hv(y)=\int_0^1\partial_\ell v(y+the_\ell)\,dt

から、切断関数の台を含む固定した半直方体上で∥δhv∥2≤M\|\delta_hv\|_2\leq Mである。この表示はH1H^1関数について、ほとんどすべての接方向の直線上で成り立つ。

一様楕円性により、弱形式の主要項の実部はθX2\theta X^2以上である。∇η\nabla\etaとδhA\delta_hAを含む残りの項の絶対値はCM(X+M)C M(X+M)以下である。また、右辺は

∣∫F0ψ‾∣≤C∥F0∥2(X+M),∣∫F⋅∇ψ‾∣=∣∫δhF⋅∇(η2δhv)‾∣≤C∥F∥H1(X+M)\left|\int F_0\overline\psi\right| \leq C\|F_0\|_2(X+M), \qquad \left|\int F\cdot\overline{\nabla\psi}\right| =\left|\int\delta_hF\cdot\overline{\nabla(\eta^2\delta_hv)}\right| \leq C\|F\|_{H^1}(X+M)

を満たす。最初の不等式には、同じ線分上の積分表示をη2δhv\eta^2\delta_hvに適用した。

Robin の境界項は差分による部分積分の後に

∫ΓR(β(⋅+heℓ)∣ηδhTv∣2+(δhβ)Tv η2δhTv‾)\int_{\Gamma_R}\left(\beta(\cdot+he_\ell)|\eta\delta_hTv|^2 +(\delta_h\beta)Tv\,\eta^2\overline{\delta_hTv}\right)

となる。その絶対値は

C∥T(ηδhv)∥22+C∥Tv∥L2(ΓR)2≤εX2+CεM2C\|T(\eta\delta_hv)\|_2^2+C\|Tv\|_{L^2(\Gamma_R)}^2 \leq\varepsilon X^2+C_\varepsilon M^2

以下である。最後の評価では§E15.12 命題 2.4をηδhv\eta\delta_hvに適用し、vvにはトレースの有界性を用いた。ε\varepsilonと Young の不等式の定数を小さく取って主要項へ吸収すると、

∥∇δhv∥L2(Qr)≤C(∥F0∥2+∥F∥H1+M)\|\nabla\delta_hv\|_{L^2(Q_r)} \leq C\left(\|F_0\|_2+\|F\|_{H^1}+M\right)

を得る。

§E15.22 補題 1.2 (2)をQrQ_rの内部の各相対コンパクト開集合に適用する。評価の右辺が境界まで一様であるから、内部からの増大列で極限を取ると、∂ℓ∂jv∈L2(Qr)\partial_\ell\partial_jv\in L^2(Q_r)がすべてのℓ<n\ell<nとj≤nj\leq nについて従う。内部の試験関数に対する方程式を展開すると、超関数として

Ann∂n2v=−F0+div⁡F−∑(i,j)≠(n,n)Aij∂i∂jv−∑i,j(∂iAij)∂jvA^{nn}\partial_n^2v =-F_0+\operatorname{div}F -\sum_{(i,j)\ne(n,n)}A^{ij}\partial_i\partial_jv -\sum_{i,j}(\partial_iA^{ij})\partial_jv

である。右辺はL2(Qr)L^2(Q_r)に属する。ℜAnn≥θ\Re A^{nn}\geq\thetaであり、1/Ann∈W1,∞1/A^{nn}\in W^{1,\infty}であるから、∂n2v\partial_n^2vもL2(Qr)L^2(Q_r)に属し、同じ右辺による評価を満たす。▨

2 Dirichlet 条件と二階正則性

定理 2.1.n≥1n\geq1、Ω⊆Rn\Omega\subseteq\R^nを有界C1,1C^{1,1}開集合とする。A∈W1,∞(Ω;Cn×n)A\in W^{1,\infty}(\Omega;\C^{n\times n})、b∈L∞(Ω;Cn)b\in L^\infty(\Omega;\C^n)、c∈L∞(Ω)c\in L^\infty(\Omega)とし、あるθ>0\theta>0に対して

ℜ(A(x)ξ⋅ξ‾)≥θ∣ξ∣2ほとんどすべての x∈Ω,ξ∈Cn\Re(A(x)\xi\cdot\overline\xi)\geq\theta|\xi|^2 \quad\text{ほとんどすべての }x\in\Omega,\quad\xi\in\C^n

とする。

Lu=−div⁡(A∇u)+b⋅∇u+cuLu=-\operatorname{div}(A\nabla u)+b\cdot\nabla u+cu

に対して、u∈H01(Ω)u\in H^1_0(\Omega)がLu=f∈L2(Ω)Lu=f\in L^2(\Omega)を超関数の意味で満たすとする。このとき、u∈H2(Ω)u\in H^2(\Omega)であり、

∥u∥H2(Ω)≤C(∥f∥L2(Ω)+∥u∥L2(Ω))\|u\|_{H^2(\Omega)}\leq C\left(\|f\|_{L^2(\Omega)}+\|u\|_{L^2(\Omega)}\right)

が成り立つ。定数CCはn,Ω,θn,\Omega,\thetaと上記の係数ノルムにのみ依存する。

証明.n≥2n\geq2の場合、境界近傍でΩ\Omegaをxn>φ(x′)x_n>\varphi(x')と表し、

Φ(y′,t)=(y′,t+φ(y′)),J=DΦ=(I0(∇φ)t1),v=u∘Φ\Phi(y',t)=(y',t+\varphi(y')),\qquad J=D\Phi=\begin{pmatrix}I&0\\(\nabla\varphi)^t&1\end{pmatrix},\qquad v=u\circ\Phi

とおく。det⁡J=1\det J=1であり、変換後の主要係数は

A~=J−1(A∘Φ)J−t\widetilde A=J^{-1}(A\circ\Phi)J^{-t}

である。φ∈C1,1\varphi\in C^{1,1}であるからA~∈W1,∞\widetilde A\in W^{1,\infty}であり、

ℜ(A~ξ⋅ξ‾)≥θ∣J−tξ∣2≥θ~∣ξ∣2\Re(\widetilde A\xi\cdot\overline\xi) \geq\theta|J^{-t}\xi|^2\geq\widetilde\theta|\xi|^2

が、座標近傍に依存するθ~>0\widetilde\theta>0に対して成り立つ。零トレースは座標変換で保たれる。変換後の方程式に補題 1.1をF=0F=0、

F0=(f−b⋅∇u−cu)∘ΦF_0=(f-b\cdot\nabla u-cu)\circ\Phi

として適用すると、各小さい境界近傍でvvのH2H^2ノルムが∥f∥2+∥u∥H1\|f\|_2+\|u\|_{H^1}によって抑えられる。§E15.9 補題 4.1により、同じ評価が元の座標のuuに対して成り立つ。境界の有限被覆と、残る内部のコンパクト集合に対する§E15.22 定理 2.2を合わせると、

∥u∥H2(Ω)≤C(∥f∥2+∥u∥H1(Ω))\|u\|_{H^2(\Omega)}\leq C(\|f\|_2+\|u\|_{H^1(\Omega)})

を得る。n=1n=1の場合には、有限個の区間上で

Au′′=−f+(b−A′)u′+cuA u''=-f+(b-A')u'+cu

と∣A∣≥θ|A|\geq\thetaを用いると同じ評価を得る。

H01H^1_0における弱形式でuu自身を試験関数に取る。§E15.16 定理 2.1により、あるα>0\alpha>0とκ≥0\kappa\geq0に対して

α∥u∥H12≤ℜ∫Ωfu‾+κ∥u∥22≤12∥f∥22+(κ+12)∥u∥22\alpha\|u\|_{H^1}^2 \leq\Re\int_\Omega f\overline u+\kappa\|u\|_2^2 \leq\tfrac12\|f\|_2^2+(\kappa+\tfrac12)\|u\|_2^2

である。得られたH1H^1評価をH2H^2評価に代入する。▨

定理 2.2.Ω⊆Rn\Omega\subseteq\R^nを有界C1,1C^{1,1}開集合とする。このとき、形式

q(u,v)=∫Ω∇u⋅∇v‾,u,v∈H01(Ω)q(u,v)=\int_\Omega\nabla u\cdot\overline{\nabla v},\qquad u,v\in H^1_0(\Omega)

に対応するL2(Ω)L^2(\Omega)上の Dirichlet Laplace 作用素ADA_Dは

D(AD)=H2(Ω)∩H01(Ω),ADu=−ΔuD(A_D)=H^2(\Omega)\cap H^1_0(\Omega),\qquad A_Du=-\Delta u

を満たす。

証明.u∈D(AD)u\in D(A_D)ならば、u∈H01(Ω)u\in H^1_0(\Omega)と−Δu=ADu∈L2-\Delta u=A_Du\in L^2が成り立つ。定理 2.1からu∈H2u\in H^2である。逆にu∈H2∩H01u\in H^2\cap H^1_0ならば、§E15.16 命題 1.7をA=IA=I、低階係数を零、ΓD=∂Ω\Gamma_D=\partial\Omegaとして適用することにより、

q(u,v)=∫Ω(−Δu)v‾(v∈H01(Ω))q(u,v)=\int_\Omega(-\Delta u)\overline v\quad(v\in H^1_0(\Omega))

を得る。これはu∈D(AD)u\in D(A_D)とADu=−ΔuA_Du=-\Delta uを意味する。▨

例 2.3.π<ω<2π\pi<\omega<2\piとし、平面の扇形

Ω={(rcos⁡ϑ,rsin⁡ϑ):0<r<1, 0<ϑ<ω}\Omega=\{(r\cos\vartheta,r\sin\vartheta):0<r<1,\ 0<\vartheta<\omega\}

を取る。α=π/ω∈(1/2,1)\alpha=\pi/\omega\in(1/2,1)とし、χ∈C∞([0,1])\chi\in C^\infty([0,1])を、r<1/4r<1/4で11、r>1/2r>1/2で00となるように取る。このとき

u(r,ϑ)=χ(r)rαsin⁡(αϑ)u(r,\vartheta)=\chi(r)r^\alpha\sin(\alpha\vartheta)

はH01(Ω)H^1_0(\Omega)に属し、−Δu∈L2(Ω)-\Delta u\in L^2(\Omega)であるが、u∉H2(Ω)u\notin H^2(\Omega)である。

実際、頂点の近傍では∣∇u∣=αrα−1|\nabla u|=\alpha r^{\alpha-1}であり、∫01/4r2α−1 dr<∞\int_0^{1/4}r^{2\alpha-1}\,dr<\inftyである。二本の直線辺ではu=0u=0、円弧の近傍でもu=0u=0であるから、零トレースの特徴づけ§E15.14 定理 2.1によりu∈H01u\in H^1_0である。極座標の Laplace 作用素から

−Δu=−(χ′′rα+(2α+1)χ′rα−1)sin⁡(αϑ)-\Delta u=-\left(\chi''r^\alpha+(2\alpha+1)\chi'r^{\alpha-1}\right)\sin(\alpha\vartheta)

であり、外力は頂点の近傍で零となる滑らかな関数である。一方、頂点の近傍では Hessian の二乗ノルムが

∣D2u∣2=2α2(α−1)2r2α−4|D^2u|^2=2\alpha^2(\alpha-1)^2r^{2\alpha-4}

となり、∫01/4r2α−3 dr=∞\int_0^{1/4}r^{2\alpha-3}\,dr=\inftyである。頂点の内角がπ\piより大きいため、この Lipschitz 開集合の境界はC1,1C^{1,1}ではない。

3 境界値と法微分の持ち上げ

定義 3.1.k≥0k\geq0を整数、Ω⊆Rn\Omega\subseteq\R^nを有界Ck+1,1C^{k+1,1}開集合とする。n≥2n\geq2のとき、有限個のグラフ座標γj:Bj→∂Ω\gamma_j:B_j\to\partial\Omegaと、境界上で和が11となる滑らかな切断関数χj\chi_jを取る。s=k+1/2s=k+1/2またはs=k+3/2s=k+3/2に対し、g∈Hs(∂Ω)g\in H^s(\partial\Omega)であるとは、各(χjg)∘γj(\chi_jg)\circ\gamma_jをBjB_jの外で零とした関数gjg_jがHs(Rn−1)H^s(\R^{n-1})に属することとする。ノルムを

∥g∥Hs(∂Ω)2=∑j∥gj∥Hs(Rn−1)2\|g\|_{H^s(\partial\Omega)}^2=\sum_j\|g_j\|_{H^s(\R^{n-1})}^2

で定める。n=1n=1の場合、Hs(∂Ω)H^s(\partial\Omega)は有限個の端点上の複素数値関数の空間とし、値の絶対値の二乗和の平方根をノルムとする。いずれの場合も、Hs(∂Ω)H^s(\partial\Omega)を ss次の境界 Sobolev 空間 (boundary Sobolev space of order s) という。

補題 3.2.定義 3.1の空間は座標と切断関数の選択に依存せず、それらを変えた二つのノルムは互いを定数倍で抑える。H1/2(∂Ω)H^{1/2}(\partial\Omega)は境界上の差分ノルムによる空間と一致する。また、整数m≥0m\geq0に対してa∈Wm+1,∞(∂Ω)a\in W^{m+1,\infty}(\partial\Omega)による乗算はHm+1/2(∂Ω)H^{m+1/2}(\partial\Omega)上有界である。ただし境界はCm,1C^{m,1}級以上であり、m=0m=0の場合には Lipschitz 境界とする。境界上のWm+1,∞W^{m+1,\infty}は各座標の局所表示にこの条件を課したものである。

証明. 双 Lipschitz な座標変換Ψ\Psiに対し、距離の比較と積分の変数変換により、次の不等式が成り立つ。

∬∣w(Ψ(x))−w(Ψ(y))∣2∣x−y∣d+1 dx dy≤C∬∣w(X)−w(Y)∣2∣X−Y∣d+1 dX dY,d=n−1,\iint\frac{|w(\Psi(x))-w(\Psi(y))|^2}{|x-y|^{d+1}}\,dx\,dy \leq C\iint\frac{|w(X)-w(Y)|^2}{|X-Y|^{d+1}}\,dX\,dY, \qquad d=n-1,

積awawの差をa(x)(w(x)−w(y))+(a(x)−a(y))w(y)a(x)(w(x)-w(y))+(a(x)-a(y))w(y)に分けると、aaが有界 Lipschitz 関数である場合のH1/2H^{1/2}乗算評価を得る。

m≥1m\geq1の場合、§E15.13 命題 3.2によってHm+1/2H^{m+1/2}ノルムをHmH^mノルムとmm階微分のH1/2H^{1/2}半ノルムで表す。Dm(w∘Ψ)D^m(w\circ\Psi)の連鎖律の各項は、(Djw)∘Ψ(D^jw)\circ\Psiと、Ψ\Psiのmm階までの微分の積である。Ψ∈Cm,1\Psi\in C^{m,1}ならば後者は有界 Lipschitz 関数であり、前者には上のH1/2H^{1/2}座標変換評価を適用することができる。逆座標にも同じ評価を適用する。積の Leibniz 公式では、Dα(aw)D^\alpha(aw)の各係数DγaD^\gamma aは∣γ∣≤m|\gamma|\leq mのとき Lipschitz 連続であるから、同じ乗算評価が各項を抑える。

座標の重なりごとに切断関数を挿入してこれらの評価を足すと、二つの有限座標系のノルムが互いを抑える。s=1/2s=1/2の境界差分ノルムとの同値は§E15.14 定理 2.1による。n=1n=1の場合は有限次元のノルムの比較と成分ごとの乗算である。▨

補題 3.3.Ω⊆Rn\Omega\subseteq\R^nを有界C1,1C^{1,1}開集合、A∈W1,∞(Ω;Cn×n)A\in W^{1,\infty}(\Omega;\C^{n\times n})とする。このとき、外向き単位法線ν\nu、行ベクトルν⋅A\nu\cdot A、関数aν=ν⋅Aνa_\nu=\nu\cdot A\nuは境界上で Lipschitz 連続な代表をもつ。さらに

ℜ(A(x)ξ⋅ξ‾)≥θ∣ξ∣2(ξ∈Cn)\Re(A(x)\xi\cdot\overline\xi)\geq\theta|\xi|^2 \quad(\xi\in\C^n)

がほとんどすべてのx∈Ωx\in\Omegaで成り立つならば、ℜaν≥θ\Re a_\nu\geq\thetaであり、aν−1a_\nu^{-1}も境界上で Lipschitz 連続である。これらの関数の各成分による乗算はH1/2(∂Ω)H^{1/2}(\partial\Omega)上有界である。

証明.Ω\Omegaがxn>φ(x′)x_n>\varphi(x')で表される座標では、

ν=(∇φ,−1)1+∣∇φ∣2\nu=\frac{(\nabla\varphi,-1)}{\sqrt{1+|\nabla\varphi|^2}}

である。∇φ\nabla\varphiは Lipschitz 連続であり、分母は11以上である。A∘ΦA\circ\Phiは半直方体上のW1,∞W^{1,\infty}関数であるから、その Lipschitz 連続な代表を底面へ制限することができる。積の差の公式によってν⋅A\nu\cdot Aとaνa_\nuの Lipschitz 連続性が従う。楕円性を境界へ連続に延長し、実単位ベクトルν\nuを代入するとℜaν≥θ\Re a_\nu\geq\thetaである。したがって

∣aν(x)−1−aν(y)−1∣≤θ−2∣aν(x)−aν(y)∣|a_\nu(x)^{-1}-a_\nu(y)^{-1}| \leq\theta^{-2}|a_\nu(x)-a_\nu(y)|

となる。有限被覆を用いると境界全体での Lipschitz 評価を得る。乗算の主張は§E15.14 命題 1.2による。n=1n=1の場合には各端点でν=±1\nu=\pm1として同じ結論を得る。▨

補題 3.4.k≥0k\geq0を整数、Ω⊆Rn\Omega\subseteq\R^nを有界Ck+1,1C^{k+1,1}開集合とする。このとき

Hk+2(Ω)⟶Hk+3/2(∂Ω)×Hk+1/2(∂Ω),u⟼(Tu,∂νu)H^{k+2}(\Omega)\longrightarrow H^{k+3/2}(\partial\Omega)\times H^{k+1/2}(\partial\Omega), \qquad u\longmapsto(Tu,\partial_\nu u)

は有界線形全射であり、有界線形右逆をもつ。ここで∂νu=ν⋅T∇u\partial_\nu u=\nu\cdot T\nabla uである。

さらにA∈Wk+1,∞(Ω;Cn×n)A\in W^{k+1,\infty}(\Omega;\C^{n\times n})が定数θ>0\theta>0に対する一様楕円性を満たすならば、任意のg∈Hk+1/2(∂Ω)g\in H^{k+1/2}(\partial\Omega)に対し、線形にw∈Hk+2(Ω)w\in H^{k+2}(\Omega)を選び、

Tw=0,ν⋅AT∇w=g,∥w∥Hk+2≤C∥g∥Hk+1/2Tw=0,\qquad \nu\cdot A T\nabla w=g,\qquad \|w\|_{H^{k+2}}\leq C\|g\|_{H^{k+1/2}}

とすることができる。

証明.n≥2n\geq2とする。半空間で接方向 Fourier 変換を

w^(ξ′,t)=∫Rn−1e−2πiy′⋅ξ′w(y′,t) dy′\widehat w(\xi',t)=\int_{\R^{n-1}}e^{-2\pi i y'\cdot\xi'}w(y',t)\,dy'

と定め、a=(1+4π2∣ξ′∣2)1/2a=(1+4\pi^2|\xi'|^2)^{1/2}とおく。境界値h∈Hk+3/2(Rn−1)h\in H^{k+3/2}(\R^{n-1})と垂直微分d∈Hk+1/2(Rn−1)d\in H^{k+1/2}(\R^{n-1})に対して

w^(ξ′,t)=(1+at)e−ath^(ξ′)+te−atd^(ξ′),t>0,\widehat w(\xi',t)=(1+at)e^{-at}\widehat h(\xi')+te^{-at}\widehat d(\xi'),\qquad t>0,

と定める。二つのプロファイルはt=0t=0でそれぞれ値と微分が(1,0)(1,0)、(0,1)(0,1)となる。接方向と垂直方向を合わせたjj階微分の二乗積分は、ttの積分でs=ats=atと置くと、

Cj∫Rn−1(a2j−1∣h^∣2+a2j−3∣d^∣2) dξ′C_j\int_{\R^{n-1}}\left(a^{2j-1}|\widehat h|^2+a^{2j-3}|\widehat d|^2\right)\,d\xi'

以下である。a≥1a\geq1であるから、0≤j≤k+20\leq j\leq k+2の和は

Ck(∥h∥Hk+3/22+∥d∥Hk+1/22)C_k\left(\|h\|_{H^{k+3/2}}^2+\|d\|_{H^{k+1/2}}^2\right)

以下となる。滑らかなデータでトレースを計算し、データの近似とトレースの連続性で一般のデータへ延長する。

局所座標Φ(y′,t)=(y′,t+φ(y′))\Phi(y',t)=(y',t+\varphi(y'))でs=(1+∣∇φ∣2)1/2s=(1+|\nabla\varphi|^2)^{1/2}とおく。Tu=hTu=hかつT∂t(u∘Φ)=dT\partial_t(u\circ\Phi)=dならば、連鎖律から

∂νu=∇φ⋅∇′hs−sd\partial_\nu u=\frac{\nabla\varphi\cdot\nabla'h}{s}-sd

である。所与の法微分qqに対して

d=∇φ⋅∇′hs2−qsd=\frac{\nabla\varphi\cdot\nabla'h}{s^2}-\frac q s

と選ぶ。∇φ,s−1∈Wk+1,∞\nabla\varphi,s^{-1}\in W^{k+1,\infty}であるから、補題 3.2により∥d∥Hk+1/2\|d\|_{H^{k+1/2}}は∥h∥Hk+3/2+∥q∥Hk+1/2\|h\|_{H^{k+3/2}}+\|q\|_{H^{k+1/2}}によって抑えられる。局所化したデータ(χjh,χjq)(\chi_jh,\chi_jq)に半空間の構成を適用し、データの台の近傍で11となる空間切断関数を掛けて元の座標へ戻す。§E15.9 補題 4.1によってHk+2H^{k+2}ノルムが保たれる。局所解を足すと、トレースの対が(h,q)(h,q)となる有界線形右逆を得る。

トレースの有界性は、§E15.9 定理 4.3による整数階の拡張、§E15.9 補題 4.1による平坦化、§E15.13 定理 6.1による値と一階微分のトレース評価から従う。法線の乗算には補題 3.2を用いる。滑らかな近似によって∇τTu=T∇τu\nabla_\tau Tu=T\nabla_\tau uも成り立つ。n=1n=1の場合には、各端点の近傍で値と一階微分を指定する一次多項式に、端点の近傍で11となる切断関数を掛け、互いに素な台をもつように足せばよい。トレースの評価は区間上の一次元の積分表示を各導関数に適用して得る。

零 Dirichlet トレースを課すとT∇w=ν∂νwT\nabla w=\nu\partial_\nu wである。したがって余法線条件はaν∂νw=ga_\nu\partial_\nu w=gとなる。aν−1∈Wk+1,∞(∂Ω)a_\nu^{-1}\in W^{k+1,\infty}(\partial\Omega)は、補題 3.3の逆数の計算をk+1k+1階まで微分することによって従う。境界データの対(0,aν−1g)(0,a_\nu^{-1}g)に既に構成した右逆を適用する。▨

系 3.5.定理 2.1と同じ領域、係数および一様楕円性の仮定の下で、u∈H1(Ω)u\in H^1(\Omega)がLu=f∈L2(Ω)Lu=f\in L^2(\Omega)とTu=h∈H3/2(∂Ω)Tu=h\in H^{3/2}(\partial\Omega)を満たすならば、u∈H2(Ω)u\in H^2(\Omega)であり、

∥u∥H2≤C(∥f∥2+∥h∥H3/2+∥u∥2)\|u\|_{H^2}\leq C\left(\|f\|_2+\|h\|_{H^{3/2}}+\|u\|_2\right)

が成り立つ。

証明.補題 3.4により、Tw=hTw=hと∥w∥H2≤C∥h∥H3/2\|w\|_{H^2}\leq C\|h\|_{H^{3/2}}を満たすwwを取る。§E15.14 定理 2.1によりu−w∈H01u-w\in H^1_0であり、L(u−w)=f−Lw∈L2L(u-w)=f-Lw\in L^2である。定理 2.1をu−wu-wに適用し、∥Lw∥2≤C∥w∥H2\|Lw\|_2\leq C\|w\|_{H^2}を用いる。▨

4 Neumann 条件と Robin 条件

定理 4.1.n≥1n\geq1、Ω⊆Rn\Omega\subseteq\R^nを有界C1,1C^{1,1}開集合とする。A∈W1,∞(Ω;Cn×n)A\in W^{1,\infty}(\Omega;\C^{n\times n})、b∈L∞(Ω;Cn)b\in L^\infty(\Omega;\C^n)、c∈L∞(Ω)c\in L^\infty(\Omega)とし、ほとんどすべてのx∈Ωx\in\Omegaとすべてのξ∈Cn\xi\in\C^nに対して

ℜ(A(x)ξ⋅ξ‾)≥θ∣ξ∣2,θ>0,\Re(A(x)\xi\cdot\overline\xi)\geq\theta|\xi|^2,\qquad\theta>0,

とする。β∈W1,∞(∂Ω)\beta\in W^{1,\infty}(\partial\Omega)、f∈L2(Ω)f\in L^2(\Omega)、g∈H1/2(∂Ω)g\in H^{1/2}(\partial\Omega)を固定する。u∈H1(Ω)u\in H^1(\Omega)が

Lu=−div⁡(A∇u)+b⋅∇u+cu=f,γν(A∇u)+βTu=gLu=-\operatorname{div}(A\nabla u)+b\cdot\nabla u+cu=f, \qquad \gamma_\nu(A\nabla u)+\beta Tu=g

を満たすとする。最初の等式は超関数の意味、第二の等式はH−1/2(∂Ω)H^{-1/2}(\partial\Omega)の意味であり、ggは⟨g,z⟩=∫∂Ωgz‾\langle g,z\rangle=\int_{\partial\Omega}g\overline zによって反双対の元とする。このときu∈H2(Ω)u\in H^2(\Omega)であり、

∥u∥H2(Ω)≤C(∥f∥L2(Ω)+∥g∥H1/2(∂Ω)+∥u∥L2(Ω))\|u\|_{H^2(\Omega)}\leq C\left(\|f\|_{L^2(\Omega)}+\|g\|_{H^{1/2}(\partial\Omega)}+\|u\|_{L^2(\Omega)}\right)

が成り立つ。CCは領域、楕円性の定数、係数の上記ノルムにのみ依存する。β=0\beta=0の場合が純粋な Neumann 条件である。

証明.div⁡(A∇u)=b⋅∇u+cu−f∈L2\operatorname{div}(A\nabla u)=b\cdot\nabla u+cu-f\in L^2であるから、弱い流束は§E15.14 定理 4.2によって定まる。仮定と§E15.14 系 4.3により、次の弱形式が成り立つ。

∫Ω(A∇u⋅∇v‾+(b⋅∇u+cu)v‾)+∫∂ΩβTuTv‾=∫Ωfv‾+∫∂ΩgTv‾(v∈H1(Ω))\int_\Omega\left(A\nabla u\cdot\overline{\nabla v} +(b\cdot\nabla u+cu)\overline v\right) +\int_{\partial\Omega}\beta Tu\overline{Tv} =\int_\Omega f\overline v+\int_{\partial\Omega}g\overline{Tv} \quad(v\in H^1(\Omega))

境界座標Φ(y′,t)=(y′,t+φ(y′))\Phi(y',t)=(y',t+\varphi(y'))に移ると、主要係数は定理 2.1の証明と同じA~\widetilde Aとなる。境界積分ではs=(1+∣∇φ∣2)1/2s=(1+|\nabla\varphi|^2)^{1/2}が掛かるので、Robin 係数とデータは

β~=s(β∘γ),g~=s(g∘γ)\widetilde\beta=s(\beta\circ\gamma),\qquad\widetilde g=s(g\circ\gamma)

となる。β~∈W1,∞\widetilde\beta\in W^{1,\infty}であり、局所化したg~\widetilde gは補題 3.2によってH1/2H^{1/2}に属する。補題 1.1を適用し、有限個の境界近傍と内部評価§E15.22 定理 2.2を合わせると、

∥u∥H2≤C(∥f∥2+∥g∥H1/2+∥u∥H1)\|u\|_{H^2}\leq C\left(\|f\|_2+\|g\|_{H^{1/2}}+\|u\|_{H^1}\right)

となる。n=1n=1の場合はAu′′=−f+(b−A′)u′+cuA u''=-f+(b-A')u'+cuによる評価を用いる。

v=uv=uとして§E15.16 定理 2.1を適用すると、

α∥u∥H12≤∥f∥2∥u∥2+C∥g∥H1/2∥u∥H1+κ∥u∥22\alpha\|u\|_{H^1}^2 \leq\|f\|_2\|u\|_2+C\|g\|_{H^{1/2}}\|u\|_{H^1}+\kappa\|u\|_2^2

を得る。Young の不等式で∥u∥H12\|u\|_{H^1}^2の小さい倍を左辺へ移すと、H1H^1ノルムも∥f∥2+∥g∥H1/2+∥u∥2\|f\|_2+\|g\|_{H^{1/2}}+\|u\|_2によって抑えられる。▨

系 4.2.Ω⊆Rn\Omega\subseteq\R^nを有界C1,1C^{1,1}開集合とする。このとき、形式

q(u,v)=∫Ω∇u⋅∇v‾,u,v∈H1(Ω)q(u,v)=\int_\Omega\nabla u\cdot\overline{\nabla v},\qquad u,v\in H^1(\Omega)

に対応するL2(Ω)L^2(\Omega)上の Neumann Laplace 作用素ANA_Nは

D(AN)={u∈H2(Ω):∂νu=0},ANu=−ΔuD(A_N)=\{u\in H^2(\Omega):\partial_\nu u=0\},\qquad A_Nu=-\Delta u

を満たす。

証明.u∈D(AN)u\in D(A_N)ならば、弱形式から−Δu=ANu∈L2-\Delta u=A_Nu\in L^2とγν(∇u)=0\gamma_\nu(\nabla u)=0を得る。定理 4.1によってu∈H2u\in H^2である。∇u∈H1\nabla u\in H^1に対する滑らかな近似と§E15.14 系 4.3により、γν(∇u)=ν⋅T∇u=∂νu\gamma_\nu(\nabla u)=\nu\cdot T\nabla u=\partial_\nu uである。逆にu∈H2u\in H^2と∂νu=0\partial_\nu u=0が成り立つならば、§E15.16 命題 1.7をA=IA=I、低階係数とβ\betaを零、ΓD=∅\Gamma_D=\emptysetとして適用すると、q(u,v)=∫(−Δu)v‾q(u,v)=\int(-\Delta u)\overline vがすべてのv∈H1v\in H^1に対して成り立つ。▨

例 4.3. 上半円Ω={(rcos⁡ϑ,rsin⁡ϑ):0<r<1, 0<ϑ<π}\Omega=\{(r\cos\vartheta,r\sin\vartheta):0<r<1,\ 0<\vartheta<\pi\}で、正の水平半径と円弧に Dirichlet 条件、負の水平半径に Neumann 条件を課す。原点の近傍で11、r>1/2r>1/2で00となる半径方向の滑らかな切断関数χ\chiに対して

u(r,ϑ)=χ(r)r1/2sin⁡(ϑ/2)u(r,\vartheta)=\chi(r)r^{1/2}\sin(\vartheta/2)

とおく。このとき、uuはH1(Ω)H^1(\Omega)に属し、指定した Dirichlet 部分で零となる。負の水平半径では

∂νu=r−1∂ϑu=12χ(r)r−1/2cos⁡(π/2)=0\partial_\nu u=r^{-1}\partial_\vartheta u =\tfrac12\chi(r)r^{-1/2}\cos(\pi/2)=0

である。f=−Δuf=-\Delta uは原点の近傍で零となり、切断関数が変化する環状領域に台をもつ滑らかな関数である。

このuuは混合条件の弱解でもある。原点から離れた試験関数には部分積分を適用することができる。原点で零となる半径方向の対数切断ηε\eta_\varepsilonを、r≤ε2r\leq\varepsilon^2で00、r≥εr\geq\varepsilonで11、その間でlog⁡(r/ε2)/log⁡(1/ε)\log(r/\varepsilon^2)/\log(1/\varepsilon)と定めると、

∫Ω∣∇ηε∣2≤πlog⁡(1/ε)⟶0\int_\Omega|\nabla\eta_\varepsilon|^2\leq\frac{\pi}{\log(1/\varepsilon)}\longrightarrow0

である。有界なH1H^1試験関数vvについてηεv→v\eta_\varepsilon v\to vがH1H^1で成り立つ。任意のH1H^1試験関数は実部と虚部をそれぞれ切り詰めることによって有界なものからH1H^1で近似され、Dirichlet 部分での零トレースも保たれる。したがって弱形式はすべての許された試験関数に延長される。

原点の近傍で∣D2u∣2=18r−3|D^2u|^2=\tfrac18r^{-3}であるから、その二乗積分は発散する。境界は条件の接合点で直線であるが、uuはその接合点のどの近傍でもH2H^2に属さない。純粋な境界条件を混合条件に置き換えると、境界の平坦さだけでは二階正則性は従わない。

5 高階の境界正則性

定理 5.1.k≥0k\geq0を整数、Ω⊆Rn\Omega\subseteq\R^nを有界Ck+1,1C^{k+1,1}開集合とする。A∈Wk+1,∞(Ω;Cn×n)A\in W^{k+1,\infty}(\Omega;\C^{n\times n})、b∈Wk,∞(Ω;Cn)b\in W^{k,\infty}(\Omega;\C^n)、c∈Wk,∞(Ω)c\in W^{k,\infty}(\Omega)とし、ほとんどすべてのx∈Ωx\in\Omegaとすべてのξ∈Cn\xi\in\C^nに対して

ℜ(A(x)ξ⋅ξ‾)≥θ∣ξ∣2,θ>0,\Re(A(x)\xi\cdot\overline\xi)\geq\theta|\xi|^2,\qquad\theta>0,

とする。L=−div⁡(A∇)+b⋅∇+cL=-\operatorname{div}(A\nabla)+b\cdot\nabla+c、f∈Hk(Ω)f\in H^k(\Omega)とし、u∈H1(Ω)u\in H^1(\Omega)がLu=fLu=fを超関数の意味で満たすとする。このとき、次が成り立つ。

  1. Tu=h∈Hk+3/2(∂Ω)Tu=h\in H^{k+3/2}(\partial\Omega)ならばu∈Hk+2(Ω)u\in H^{k+2}(\Omega)であり、 ∥u∥Hk+2≤Ck(∥f∥Hk+∥h∥Hk+3/2+∥u∥2)\|u\|_{H^{k+2}}\leq C_k\left(\|f\|_{H^k}+\|h\|_{H^{k+3/2}}+\|u\|_2\right) が成り立つ。
  2. β∈Wk+1,∞(∂Ω)\beta\in W^{k+1,\infty}(\partial\Omega)とし、γν(A∇u)+βTu=g∈Hk+1/2(∂Ω)\gamma_\nu(A\nabla u)+\beta Tu=g\in H^{k+1/2}(\partial\Omega)が反双対の意味で成り立つならばu∈Hk+2(Ω)u\in H^{k+2}(\Omega)であり、 ∥u∥Hk+2≤Ck(∥f∥Hk+∥g∥Hk+1/2+∥u∥2)\|u\|_{H^{k+2}}\leq C_k\left(\|f\|_{H^k}+\|g\|_{H^{k+1/2}}+\|u\|_2\right) が成り立つ。

定数はk,n,Ω,θk,n,\Omega,\thetaと指定した係数ノルムにのみ依存する。斉次 Dirichlet 条件ではh=0h=0とし、純粋な Neumann 条件ではβ=0\beta=0とする。

証明.補題 3.4によって、Dirichlet の場合にはTw=hTw=h、Robin の場合にはTw=0Tw=0とν⋅AT∇w=g\nu\cdot A T\nabla w=gを満たすw∈Hk+2w\in H^{k+2}を取る。∥w∥Hk+2\|w\|_{H^{k+2}}はそれぞれの境界データのノルムで抑えられる。Lw∈HkLw\in H^kであるから、u−wu-wの方程式に移ると境界条件は斉次になる。以下の斉次問題の評価をu−wu-wに適用すれば、二つの主張が従う。

k=0k=0の評価は定理 2.1と定理 4.1である。k≥1k\geq1とし、次数k−1k-1までの結論が成り立つと仮定する。帰納法の仮定によりu∈Hk+1u\in H^{k+1}であり、

∥u∥Hk+1≤C(∥f∥Hk−1+∥u∥2)\|u\|_{H^{k+1}}\leq C\left(\|f\|_{H^{k-1}}+\|u\|_2\right)

である。境界をCk+1,1C^{k+1,1}グラフで平坦化し、v=u∘Φv=u\circ\Phiとおく。変換後の係数は

A~=J−1(A∘Φ)J−t,b~=J−1(b∘Φ),c~=c∘Φ,f~=f∘Φ\widetilde A=J^{-1}(A\circ\Phi)J^{-t},\quad \widetilde b=J^{-1}(b\circ\Phi),\quad \widetilde c=c\circ\Phi,\quad \widetilde f=f\circ\Phi

であり、Robin 係数はβ~=s(β∘γ)\widetilde\beta=s(\beta\circ\gamma)である。主要係数と Robin 係数はWk+1,∞W^{k+1,\infty}、低階係数はWk,∞W^{k,\infty}、外力はHkH^kに属する。これらの事実は連鎖律と§E15.9 補題 4.1から従う。

接方向の多重指数α\alphaを∣α∣=k|\alpha|=kとして固定し、z=Dταvz=D_\tau^\alpha vとおく。v∈Hk+1v\in H^{k+1}であるからz∈H1z\in H^1である。局所弱形式を接方向にα\alpha回微分すると、

∫A~∇z⋅∇ψ‾+∫Γβ~TzTψ‾=∫(F0,αψ‾+Fα⋅∇ψ‾)+∫ΓrαTψ‾,\int\widetilde A\nabla z\cdot\overline{\nabla\psi} +\int_{\Gamma}\widetilde\beta Tz\overline{T\psi} =\int\left(F_{0,\alpha}\overline\psi+F_\alpha\cdot\overline{\nabla\psi}\right) +\int_{\Gamma}r_\alpha\overline{T\psi},F0,α=Dτα(f~−b~⋅∇v−c~v),Fα=−∑0<γ≤α(αγ)(DτγA~)∇Dτα−γv,rα=−∑0<γ≤α(αγ)(Dτγβ~)TDτα−γv\begin{aligned} F_{0,\alpha}&=D_\tau^\alpha(\widetilde f-\widetilde b\cdot\nabla v-\widetilde c v),\\ F_\alpha&=-\sum_{0<\gamma\leq\alpha}\binom\alpha\gamma (D_\tau^\gamma\widetilde A)\nabla D_\tau^{\alpha-\gamma}v,\\ r_\alpha&=-\sum_{0<\gamma\leq\alpha}\binom\alpha\gamma (D_\tau^\gamma\widetilde\beta)T D_\tau^{\alpha-\gamma}v \end{aligned}

となる。Dirichlet の場合には境界積分を除き、Tz=0Tz=0である。接方向の微分とトレースの交換は、整数階の滑らかな近似とトレースの連続性から従う。弱形式の微分も、滑らかな試験関数に接方向の部分積分を行った後に、H1H^1試験関数へ連続に延長することによって正当化される。

F0,αF_{0,\alpha}の各項に現れるvvの微分は高々k+1k+1階であるから、F0,α∈L2F_{0,\alpha}\in L^2である。FαF_\alphaをさらに一階微分しても、γ≠0\gamma\ne0によってvvの微分は高々k+1k+1階であり、係数の微分は高々k+1k+1階である。したがってFα∈H1F_\alpha\in H^1である。Dτα−γv∈H2D_\tau^{\alpha-\gamma}v\in H^2であり、Dτγβ~D_\tau^\gamma\widetilde\betaは Lipschitz 連続であるから、トレース評価と補題 3.2によってrα∈H1/2r_\alpha\in H^{1/2}である。各ノルムは

C(∥f~∥Hk+∥v∥Hk+1)C\left(\|\widetilde f\|_{H^k}+\|v\|_{H^{k+1}}\right)

以下となる。補題 1.1によってDταv∈H2D_\tau^\alpha v\in H^2であり、同じ右辺による局所評価を得る。従って全次数k+2k+2の微分のうち、法方向の微分を高々二個含むものはすべてL2L^2に属する。

方程式を

A~nn∂n2v=−f~−∑(i,j)≠(n,n)A~ij∂i∂jv+∑j(b~j−∑i∂iA~ij)∂jv+c~v\widetilde A^{nn}\partial_n^2v =-\widetilde f-\sum_{(i,j)\ne(n,n)}\widetilde A^{ij}\partial_i\partial_jv +\sum_j\left(\widetilde b^j-\sum_i\partial_i\widetilde A^{ij}\right)\partial_jv +\widetilde c v

と書く。3≤q≤k+23\leq q\leq k+2とし、全次数k+2k+2で法微分がq−1q-1個以下のものが既にL2L^2に属すると仮定する。この等式に∣δ∣=k|\delta|=k、δn=q−2\delta_n=q-2を満たすDδD^\deltaを施す。左辺で係数を微分しない項はA~nnDδ∂n2v\widetilde A^{nn}D^\delta\partial_n^2vである。右辺の全次数k+2k+2の各項は法微分が高々q−1q-1個であり、係数を一回以上微分する項に現れるvvの微分は高々k+1k+1階である。Dδf~∈L2D^\delta\widetilde f\in L^2であるから、右辺と左辺の交換子はすべて既に評価されたL2L^2関数である。∣A~nn∣≥θ~|\widetilde A^{nn}|\geq\widetilde\thetaによって割ると、法微分がqq個の項もL2L^2に属する。qqに関する帰納法により、v∈Hk+2v\in H^{k+2}とその局所評価を得る。

境界近傍では§E15.9 補題 4.1で元の座標へ戻し、内部には§E15.22 定理 3.1を適用する。有限被覆上の評価を足すと

∥u∥Hk+2≤Ck(∥f∥Hk+∥u∥Hk+1)\|u\|_{H^{k+2}}\leq C_k\left(\|f\|_{H^k}+\|u\|_{H^{k+1}}\right)

となる。帰納法のHk+1H^{k+1}評価を代入して斉次問題の結論を得る。n=1n=1の場合にはAu′′=−f+(b−A′)u′+cuA u''=-f+(b-A')u'+cuをkk回微分し、u∈Hk+1u\in H^{k+1}と∣A∣≥θ|A|\geq\thetaを使うと、同じ帰納的評価が成り立つ。▨

6 Laplace 作用素の冪と固有関数

系 6.1.Ω⊆Rn\Omega\subseteq\R^nを有界C∞C^\infty開集合、ADA_DをL2(Ω)L^2(\Omega)上の Dirichlet Laplace 作用素とする。このとき、任意の整数m≥1m\geq1に対して

D(ADm)={u∈H2m(Ω):TΔju=0(0≤j<m)},ADmu=(−Δ)muD(A_D^m)=\{u\in H^{2m}(\Omega):T\Delta^j u=0\quad(0\leq j<m)\}, \qquad A_D^m u=(-\Delta)^m u

が成り立つ。この空間上で∥u∥H2m\|u\|_{H^{2m}}と∥u∥2+∥ADmu∥2\|u\|_2+\|A_D^m u\|_2は同値なノルムである。

証明.m=1m=1は定理 2.2である。mmに対する主張を仮定し、

D(ADm+1)={u∈D(AD):ADu∈D(ADm)}D(A_D^{m+1})=\{u\in D(A_D):A_Du\in D(A_D^m)\}

を用いる。右辺のuuはH01H^1_0に属し、−Δu=ADu∈H2m-\Delta u=A_Du\in H^{2m}を満たす。定理 5.1によってu∈H2m+2u\in H^{2m+2}である。Tu=0Tu=0と、ADuA_Duに対する帰納法の境界条件からTΔju=0T\Delta^j u=0が0≤j≤m0\leq j\leq mで成り立つ。逆にu∈H2m+2u\in H^{2m+2}がこれらの境界条件を満たすならば、u∈D(AD)u\in D(A_D)であり、−Δu∈H2m-\Delta u\in H^{2m}は次数mmの帰納法の境界条件を満たす。従ってADu∈D(ADm)A_Du\in D(A_D^m)である。

高階評価をu,ADu,…,ADm−1uu,A_Du,\ldots,A_D^{m-1}uに適用すると、

∥u∥H2m≤Cm∑j=0m∥ADju∥2\|u\|_{H^{2m}}\leq C_m\sum_{j=0}^m\|A_D^j u\|_2

となる。§E15.19 定理 1.3の完全正規直交系を(eℓ)(e_\ell)、非負固有値を(λℓ)(\lambda_\ell)とすると、uℓ=(u,eℓ)2u_\ell=(u,e_\ell)_2に対し

∥ADju∥22=∑ℓλℓ2j∣uℓ∣2≤2∑ℓ(1+λℓ2m)∣uℓ∣2=2(∥u∥22+∥ADmu∥22)\|A_D^ju\|_2^2=\sum_\ell\lambda_\ell^{2j}|u_\ell|^2 \leq2\sum_\ell(1+\lambda_\ell^{2m})|u_\ell|^2 =2\left(\|u\|_2^2+\|A_D^mu\|_2^2\right)

である。逆向きの評価は、ADmu=(−Δ)muA_D^mu=(-\Delta)^muが2m2m階の定数係数微分作用素であることから従う。▨

系 6.2.Ω⊆Rn\Omega\subseteq\R^nを有界C∞C^\infty開集合とし、ADϕ=λϕA_D\phi=\lambda\phi、∥ϕ∥L2=1\|\phi\|_{L^2}=1とする。このときϕ∈C∞(Ω‾)\phi\in C^\infty(\overline\Omega)であり、任意の整数r≥0r\geq0に対して

∥ϕ∥H2r(Ω)≤Cr(1+λ)r\|\phi\|_{H^{2r}(\Omega)}\leq C_r(1+\lambda)^r

が成り立つ。また、任意の整数q≥0q\geq0に対し、2r>q+n/22r>q+n/2を満たす整数rrを取れば、

∥ϕ∥Cq(Ω‾)≤Cq,r(1+λ)r\|\phi\|_{C^q(\overline\Omega)}\leq C_{q,r}(1+\lambda)^r

が成り立つ。ここでC∞(Ω‾)C^\infty(\overline\Omega)は、すべての階数の導関数が境界まで連続に延長されることを意味する。

証明.ADϕ=λϕA_D\phi=\lambda\phiであるから、作用素の冪の定義によってϕ∈D(ADr)\phi\in D(A_D^r)とADrϕ=λrϕA_D^r\phi=\lambda^r\phiがすべてのr≥1r\geq1に対して成り立つ。系 6.1のノルム同値から∥ϕ∥H2r≤Cr(1+λr)≤2Cr(1+λ)r\|\phi\|_{H^{2r}}\leq C_r(1+\lambda^r)\leq2C_r(1+\lambda)^rを得る。r=0r=0は正規化の仮定による。§E15.10 系 4.2をp=2p=2、次数2r2rとして適用すると、2r>q+n/22r>q+n/2の下でH2rH^{2r}ノルムはCqC^qノルムを抑える。qqは任意であるから、すべての階数の導関数が境界まで連続である。▨

例 6.3.Ω=(0,ℓ)\Omega=(0,\ell)、ℓ>0\ell>0とする。このとき、各j≥1j\geq1に対して、次のϕj\phi_jは固有値λj\lambda_jに属する正規化 Dirichlet 固有関数である。

ϕj(x)=2ℓsin⁡jπxℓ,λj=(jπℓ)2,j≥1,\phi_j(x)=\sqrt{\frac2\ell}\sin\frac{j\pi x}\ell, \qquad \lambda_j=\left(\frac{j\pi}\ell\right)^2, \qquad j\geq1,

整数階の Sobolev ノルムを導関数の二乗和の平方根で取ると、任意の整数s,r≥0s,r\geq0について

∥ϕj(s)∥L2(0,ℓ)=λjs/2,∥ϕj∥H2r(0,ℓ)2=∑s=02rλjs\|\phi_j^{(s)}\|_{L^2(0,\ell)}=\lambda_j^{s/2},\qquad \|\phi_j\|_{H^{2r}(0,\ell)}^2=\sum_{s=0}^{2r}\lambda_j^s

である。偶数階微分は元の正弦関数の定数倍であるから、すべて端点で零となる。

一方、u(x)=x(ℓ−x)u(x)=x(\ell-x)は滑らかでu(0)=u(ℓ)=0u(0)=u(\ell)=0であるが、u′′=−2u''=-2は端点で零とならない。したがってu∈D(AD)u\in D(A_D)である一方、系 6.1によりu∉D(AD2)u\notin D(A_D^2)である。H4H^4正則性と Dirichlet 条件だけでは、二乗の定義域に属するための境界条件は満たされない。

前提記事