1 古典解の弱最大値原理
定理 1.1 (弱最大値原理).Ω⊆Rnを空でない有界開集合とする。実対称行列A=(aij)の成分と実関数b1,…,bn,cはΩ上で有界連続であり、あるθ>0に対して
ξ⋅A(x)ξ≥θ∣ξ∣2,c(x)≥0(x∈Ω, ξ∈Rn)を満たすとする。実関数u∈C2(Ω)∩C(Ω)が
Lu:=−i,j=1∑naij∂iju+b⋅∇u+cu≤0Ω 上を満たすとする。このとき、
Ωmaxu≤max{0,∂Ωmaxu}である。c=0ならばmaxΩu≤max∂Ωuである。
証明.Lu<0と仮定する。uが正の最大値を内部点xでとるならば、∇u(x)=0であり、D2u(x)は負半定値である。A(x)を正規直交基底で対角化すると−∑aij(x)∂iju(x)≥0となる。c(x)u(x)≥0も成り立つので、Lu(x)≥0となり、Lu<0に反する。従って狭義の不等式を満たす関数について主張の評価が成り立つ。
Lu≤0の場合には、k>0を十分大きくとると
L(ekx1)=ekx1(−a11k2+b1k+c)<0がΩ全体で成り立つ。実際、括弧内は−θk2+∥b1∥∞k+∥c∥∞以下である。任意のε>0に対してuε=u+εekx1はLuε<0を満たす。狭義の場合の評価をuεに適用し、Ω上の一様収束uε→uによってε↓0とすると、主張の評価を得る。
c=0の場合にはm=max∂Ωuとおく。L(u−m)=Lu≤0と境界上のu−m≤0に既に示した評価を適用すると、u−m≤0である。▨
系 1.2.Ω,A,b,c,Lは定理 1.1の仮定を満たすとする。実関数u,v∈C2(Ω)∩C(Ω)がLu≤LvをΩ上で、u≤vを∂Ω上で満たすとする。このとき、u≤vがΩ上で成り立つ。従って、与えられたfと境界値gに対するLu=f、u∣∂Ω=gの古典解は高々一つである。
証明.L(u−v)≤0と境界上のu−v≤0に定理 1.1を適用するとu−v≤0である。同じ方程式と境界値を満たす二解には、両方の順序でこの結論を適用することができる。▨
例 1.3. 区間(0,π)上でu(x)=sinxは−u′′−u=0と零境界値を満たすが、内部では正である。従って零階係数の非負性を外すと弱最大値原理は成り立たない。
区間(−1,1)上でv(x)=−coshxは−v′′+v=0を満たす。内部の最大値はv(0)=−1であり、境界値−cosh1より大きい。従ってc>0の場合には、弱最大値原理の右辺の0を除くことはできない。
2 Hopf の境界点補題と強最大値原理
補題 2.1.Ω,A,b,c,Lは定理 1.1の仮定を満たすとし、実関数u∈C2(Ω)∩C(Ω)はLu≤0を満たすとする。r>0、y∈Ω、x0∈∂Ωが
Br(y)⊆Ω,∣x0−y∣=r,ν=rx0−yを満たすとする。さらにM=u(x0)≥0かつu(x)<Mがすべてのx∈Br(y)について成り立つならば、
t↓0liminftu(x0)−u(x0−tν)>0である。uがx0の近傍へC1関数として延長される場合には、∂νu(x0)>0である。c=0の場合にはM≥0の仮定を除くことができる。
証明.D={x:r/2<∣x−y∣<r}とし、α>0に対して
h(x)=e−α∣x−y∣2−e−αr2とおく。z=x−yと書くと、
Lh=e−α∣z∣2(2αtrA−4α2z⋅Az−2αb⋅z)+ch.D上では∣z∣≥r/2、0≤h≤e−α∣z∣2である。従って
Lh≤e−α∣z∣2(2αni,jmax∥aij∥∞−α2θr2+2αr∥b∥∞+∥c∥∞).αを十分大きく固定するとLh<0となる。
内側の球面上ではu−Mの最大値が負なので、あるε>0に対してu−M+εh≤0が成り立つ。外側の球面上でも、連続性からu−M≤0でありh=0なので、同じ不等式が成り立つ。M≥0とc≥0から
L(u−M+εh)=Lu−cM+εLh≤0である。定理 1.1をD上で適用するとu−M+εh≤0を得る。0<t<r/2に対して
tM−u(x0−tν)≥εte−α(r−t)2−e−αr2⟶2εαre−αr2>0.C1延長が存在すれば左辺は∂νu(x0)に収束する。c=0の場合にはL(u−M)=Luであるため、Mの符号によらず同じ議論が成り立つ。▨
定理 2.2 (強最大値原理).Ω⊆Rnを空でない連結な開集合とする。実対称行列Aと実係数b,cは有界連続であり、あるθ>0に対して、すべてのx∈Ωとξ∈Rnについてξ⋅A(x)ξ≥θ∣ξ∣2、c(x)≥0を満たすとする。実関数u∈C2(Ω)が
−i,j∑aij∂iju+b⋅∇u+cu≤0を満たし、あるx∗∈Ωにおいてu(x∗)=supΩu≥0とする。このとき、uは定数である。c=0の場合には最大値の非負性を仮定する必要はない。
証明.M=u(x∗)、S={x∈Ω:u(x)=M}とおく。Sは空でない相対閉集合である。S=Ωと仮定すると、連結性によりSの相対境界に点zがある。B3ρ(z)⊆Ωとなるρ>0をとり、y∈Bρ(z)∖Sをとる。r=dist(y,S)は正であり、r≤∣y−z∣<ρである。距離を最小化する列はB2ρ(z)内にとることができるので、あるx0∈Sが∣y−x0∣=rを満たす。
Br(y)⊆Ω∖Sであり、Br(y)⊆Ωである。uをBr(y)に制限して補題 2.1を適用すると、ν=(x0−y)/rに対して∂νu(x0)>0を得る。一方、x0はΩ内の最大点なので∇u(x0)=0である。この二つの結論は両立しない。従ってS=Ωである。c=0の場合には、Hopf の補題もMの符号によらず適用することができる。▨
例 2.3. 球Br(0)上のu(x)=∣x∣2は−Δu=−2n<0を満たす。任意のx0∈∂Br(0)について、外向き単位法線はν=x0/rであり、∂νu(x0)=2r>0である。
正方形(0,1)2上のv(x,y)=−xyはΔv=0を満たし、内部で負であるが、頂点(0,0)ではv=0かつ∇v=0である。この頂点に内接球は存在しない。実際、中心(a,b)、半径rの内接球が頂点を通るならばr=a2+b2である一方、二つの座標軸から離れるためにはr≤min{a,b}が必要となり、矛盾する。
互いに素な二つの球の和集合上で、一方の球では1、他方の球では0と定めた関数は調和関数である。この関数は内部で非負の大域最大値をとるが定数ではないので、強最大値原理の連結性は必要である。
3 Sobolev 弱解の比較原理
定理 3.1 (弱解の最大値原理).Ω⊆Rnを有界 Lipschitz 開集合とし、A:Ω→Rn×nを有界可測行列とする。あるθ>0に対して、ほとんどすべてのx∈Ωとすべてのξ∈Rnについてξ⋅A(x)ξ≥θ∣ξ∣2とする。c∈L∞(Ω;R)は非負とし、実関数u∈H1(Ω)が
∫ΩA∇u⋅∇φ+cuφdx≤0(0≤φ∈H01(Ω))を満たすとする。実数M≥0に対してTu≤Mが境界の表面測度に関してほとんど至る所で成り立つとする。このとき、u≤MがΩ上ほとんど至る所で成り立つ。c=0の場合には任意のM∈Rに対してこの結論が成り立つ。
証明.w=(u−M)+とおく。§E15.8 命題 3.2と§E15.8 系 3.3によりw∈H1(Ω)であり、∇w=1{u>M}∇uである。§E15.12 命題 2.5によりTw=(Tu−M)+=0なので、§E15.12 定理 3.3からw∈H01(Ω)となる。試験関数φ=wを代入すると
θ∥∇w∥22≤∫ΩA∇w⋅∇w+cuwdx≤0.ここでw>0の点ではu>M≥0なのでcuw≥0である。従って∇w=0であり、§E15.8 定理 5.4によりw=0である。c=0の場合にはcuw=0であるため、Mの符号によらず同じ結論を得る。▨
系 3.2.Ω,A,cは定理 3.1の仮定を満たすとする。実関数u,v∈H1(Ω)が
∫ΩA∇(u−v)⋅∇φ+c(u−v)φdx≤0(0≤φ∈H01(Ω))と境界上のTu≤Tvを満たすとする。このとき、u≤vがほとんど至る所で成り立つ。従って同じ外力とトレースをもつ実弱解は一意である。
証明.u−vに定理 3.1をM=0として適用する。同じ外力とトレースをもつ二解については、差とその負の両方に適用する。▨
例 3.3.R>0とし、BR(0)⊆Rn上で
u(x)=2nR2−∣x∣2とおく。このとき、−Δu=1、Tu=0であり、u∈H01(BR(0))である。零関数とuの差は系 3.2の弱不等式を満たすのでu≥0が従う。実際、表示式からu(x)>0が内部の各点で成り立つ。
4 外力による有界性評価
定理 4.1 (Stampacchia の評価).Ω⊆Rnを空でない有界 Lipschitz 開集合とする。実有界可測行列Aが、あるθ>0に対してξ⋅A(x)ξ≥θ∣ξ∣2をほとんどすべてのx∈Ωとすべてのξ∈Rnについて満たすとする。0≤c∈L∞(Ω;R)、1≤q≤∞、q>n/2、f∈Lq(Ω;R)とする。実関数u∈H01(Ω)が
∫ΩA∇u⋅∇φ+cuφdx=∫Ωfφdx(φ∈H01(Ω))を満たすとする。このとき、
∥u∥L∞(Ω)≤C∥f∥Lq(Ω)である。定数CはΩ,n,q,θのみに依存する。
証明.n≥3の場合にはr=2n/(n−2)とおく。w∈H01(Ω)に対して、wj∈Cc∞(Ω)をwj→wがH1(Ω)で成り立つようにとる。このとき、零延長E0wjはE0wへH1(Rn)で収束する。§E15.10 定理 1.2を零延長に適用して
∥w∥r≤S∥∇w∥2を得る。n=2の場合にはr>2をr(1−1/q)>2となる有限指数に選ぶことができる。§E15.10 命題 2.1と§E15.8 定理 5.4により、同じ勾配による評価が成り立つ。どちらの場合も
γ=1−q1−r1>0,β=rγ=r(1−1/q)−1>1である。ただし1/∞=0とする。
k≥0に対してEk={x∈Ω:u(x)>k}、w=(u−k)+とおく。切断とトレースに関する§E15.8 系 3.3、§E15.12 命題 2.5、§E15.12 定理 3.3により、w∈H01(Ω)であり、∇w=1Ek∇uである。uw≥0なので、弱等式にwを代入し、§E10.22 定理 4.2を用いると
θ∥∇w∥22≤∫Ek∣f∣wdx≤∥f∥q∥w∥r∣Ek∣γ≤S∥f∥q∥∇w∥2∣Ek∣γ.中間の評価は1/q+1/r+γ=1に対して Hölder の不等式を繰り返して得られる。∥∇w∥2=0の場合も含めて
∥w∥r≤K∥f∥q∣Ek∣γ,K=S2/θとなる。h>kに対してw≥h−kがEh上で成り立つので、
∣Eh∣≤(h−kK∥f∥q)r∣Ek∣β.F=∥f∥q=0の場合には、弱等式にuを代入するとθ∥∇u∥22≤0となり、Poincaré の不等式からu=0である。F>0とし、a=r/(β−1)、
d=KF21+a/r∣Ω∣(β−1)/r,kj=d(1−2−j),Yj=∣Ekj∣とおく。水準集合の測度の不等式にh=kj+1、k=kjを代入すると
Yj+1≤(KF/d)r2r(j+1)Yjβである。Y0≤∣Ω∣を初期値として帰納するとYj≤∣Ω∣2−ajとなる。実際、この評価を右辺に代入したとき、jの指数はr−aβ=−aであり、定数因子はdの選び方から2−aになる。Ed⊆Ekjがすべてのjについて成り立つので∣Ed∣=0であり、u≤dを得る。−uは外力−fに対して同じ弱等式を満たすので、−u≤dも成り立つ。
n=1の場合には、Ωを含む有界区間Iをとり、uを零延長する。H01の近似により零延長はH01(I)に属する。§E13.29 定理 1.2 (2)と§E15.8 定理 5.4から
∥u∥∞≤CI∥u′∥2を得る。弱等式にuを代入し、Hölder の不等式を用いると
θ∥u′∥22≤∥f∥1∥u∥∞≤∣Ω∣1−1/q∥f∥qCI∥u′∥2.従ってq=1を含めて主張の評価が成り立つ。▨
例 4.2.n≥3、Ω=B1(0)⊆Rnとする。0≤χ≤1を満たし、B1/4(0)上で1となるχ∈Cc∞(B1/2(0))をとる。原点以外で
u(x)=χ(x)loglog(e/∣x∣)と定め、原点での値は任意に定める。このとき、uはH01(Ω)に属し、−Δu∈Ln/2(Ω)であるが、u∈/L∞(Ω)である。
§E15.10 例 2.2によりu∈W1,n(Ω)であり、uは本質的に非有界である。n≥3と∣Ω∣<∞から、Hölder の不等式によりu∈H1(Ω)である。0<r=∣x∣<1/4、t=log(e/r)、v(r)=logtとおくと
v′(r)=−rt1,v′′(r)=r2t1−r2t21,−Δv=r2tn−2+r2t21.uをBε(0)上で零、B2ε(0)の外でそのままにする滑らかな切断で近似する。切断関数の勾配がC/ε以下となるように選ぶと、追加される勾配の二乗積分は
O(εn−2(loglog(e/ε))2)⟶0である。uと∇uは二乗可積分なので、切断によって除かれる部分の関数と勾配のL2ノルムは零へ収束する。従って切断による近似はuへH1(Ω)で収束する。この近似はCc∞(Ω)に属するので、u∈H01(Ω)である。
B1/4(0)の外ではuは滑らかである。原点近傍のf=−Δuについては
∫01/4∣f(r)∣n/2rn−1dr≤C∫01/4rlog(e/r)n/2dr<∞となる。∫Ω∖Bε∇u⋅∇φを部分積分した際の内側境界項は、任意のφ∈Cc∞(Ω)に対してO(εn−2/log(e/ε))であり、零へ収束する。従って分布としても−Δu=fであり、原点に点質量は生じない。
さらに§E15.10 定理 1.2と有界集合上の Hölder の不等式により、φ↦∫fφはH01(Ω)上連続である。実際、n/(n−2)≤2n/(n−2)なので∥φ∥n/(n−2)≤C∥∇φ∥2である。密度によって弱等式はすべてのφ∈H01(Ω)に延長される。任意のK>0に対して、ある0<ρ<1/4が存在し、0<∣x∣<ρを満たすすべてのxについてu(x)>Kである。従って原点での値によらず、uは本質的に非有界である。従って定理 4.1の指数条件をq=n/2へ弱めることはできない。