1 速度場の流れと線形輸送方程式
定義 1.1 (特性曲線).n∈N≥1、U⊂R×Rnを開集合、b∈C1(U;Rn)とする。(s,y)∈Uに対して、初期値問題
dtdX(t;s,y)=b(t,X(t;s,y)),X(s;s,y)=yの最大解をX(⋅;s,y)、その定義区間をIs,yと書く。時空の曲線t↦(t,X(t;s,y))を速度場bの 特性曲線 (characteristic curve) という。
命題 1.2.n∈N≥1、U⊂R×Rnを開集合、b∈C1(U;Rn)とし、X(⋅;s,y)を速度場bの最大解とする。s∈Rを固定し、A⊂{y∈Rn:(s,y)∈U}を開集合とする。
Ds,A={(t,y):y∈A, t∈Is,y},Ψs,A(t,y)=(t,X(t;s,y))とおく。このとき、Ds,Aは開集合であり、Ψs,Aは開集合Ωs,A=Ψs,A(Ds,A)⊂UへのC1微分同相写像である。逆写像は
Ψs,A−1(t,x)=(t,X(s;t,x))である。特に、時刻tを固定すると、y↦X(t;s,y)は{y∈A:t∈Is,y}からその像へのC1微分同相写像である。
証明.U上の自律場B(τ,x)=(1,b(τ,x))の最大流をφとする。第一成分の微分は1であるから、
φr(s,y)=(s+r,X(s+r;s,y))である。最大流の定義域は開集合であり、§E13.15 系 2.2 (1)によりφはC1級である。したがって、Ds,Aは開集合であり、(t,y)↦X(t;s,y)はC1級である。
§E13.15 定理 1.3により、t∈Is,yならばs∈It,X(t;s,y)であり、
X(s;t,X(t;s,y))=yが成り立つ。よってΨs,Aは単射である。時刻sの断面から時刻tの断面への写像は、§E13.15 系 2.2 (3)の微分同相写像とその逆を、それぞれの断面へ制限したものである。両断面の接空間はともに{0}×Rnであるから、DyX(t;s,y)は可逆である。したがって、
DΨs,A(t,y)=(1b(t,X(t;s,y))0DyX(t;s,y))も可逆である。§E4.7 定理 2.1によりΨs,Aは各点の近傍でC1微分同相写像となる。その像Ωs,Aは開集合であり、単射性によって各局所逆写像は一つのC1逆写像を定める。合成則で得た等式から、この逆写像は(t,x)↦(t,X(s;t,x))である。▨
定理 1.3.n∈N≥1、U⊂R×Rnを開集合とし、b∈C1(U;Rn)、c,f∈C1(U;R)とする。s∈R、開集合A⊂{y:(s,y)∈U}、u0∈C1(A;R)に対し、速度場bの流れからDs,A、Ωs,Aを命題 1.2のように定める。このとき、
ut+b(t,x)⋅∇xu+c(t,x)u=f(t,x),u(s,y)=u0(y)を満たすu∈C1(Ωs,A;R)は一意に存在する。C(t,y)=∫stc(τ,X(τ;s,y))dτとおくと、その解は
u(t,X(t;s,y))=e−C(t,y)[u0(y)+∫steC(τ,y)f(τ,X(τ;s,y))dτ]((t,y)∈Ds,A)で与えられる。
証明.(t,y)∈Ds,Aに対して、主張の表示式の右辺をv(t,y)とする。Is,yはsを含む区間であるから、表示式の積分に現れる軌道はすべて定義されている。流れと係数はC1級である。各(t,y)の積分区間を含む軌道のコンパクト部分を覆う近傍をとると、近傍内の初期点に対しても同じ区間上で被積分関数とその初期点微分が連続である。したがって、積分を初期点について微分することによってC,v∈C1(Ds,A;R)を得る。
各y∈Aについて、§E13.3 定理 3.2によりt↦v(t,y)はIs,y上の方程式
∂tv(t,y)+c(t,X(t;s,y))v(t,y)=f(t,X(t;s,y)),v(s,y)=u0(y)の一意な古典解である。u=v∘Ψs,A−1と定めると、命題 1.2によりu∈C1(Ωs,A;R)である。§E4.3 定理 1.1により、
∂tv(t,y)=ut(t,X(t;s,y))+b(t,X(t;s,y))⋅∇xu(t,X(t;s,y))であるから、uは所与の偏微分方程式と初期条件を満たす。
u∈C1(Ωs,A;R)も同じ偏微分方程式と初期条件を満たすとする。連鎖律によりt↦u(t,X(t;s,y))は同じ一次線形方程式と初期条件を満たす。積分因子による一意性からu∘Ψs,A=v=u∘Ψs,Aであり、Ψs,Aの全射性からu=uである。▨
例 1.4.n∈N≥1、v∈Rn、γ∈R、s∈R、u0∈C1(Rn;R)とする。このとき、速度場b(t,x)=vの流れはX(t;s,y)=y+(t−s)vである。したがって、定理 1.3により、
ut+v⋅∇xu+γu=0,u(s,x)=u0(x)のR×Rn上の解は
u(t,x)=e−γ(t−s)u0(x−(t−s)v)である。γ=0のとき、解は各特性曲線上で初期値を保つ。
2 半線形輸送方程式と存在領域
定理 2.1.n∈N≥1、U⊂R×Rnを開集合、b∈C1(U;Rn)とする。G⊂U×Rを開集合、F∈C1(G;R)とする。s∈R、開集合A⊂{y:(s,y)∈U}、u0∈C1(A;R)は、すべてのy∈Aについて(s,y,u0(y))∈Gを満たすとし、X(⋅;s,y)を速度場bの最大解とする。各y∈Aに対して、
Z′(t;y)=F(t,X(t;s,y),Z(t;y)),Z(s;y)=u0(y)の最大解をZ(⋅;y)、その定義区間をJy⊂Is,yとする。このとき、
DF={(t,y):y∈A, t∈Jy},ΩF={(t,X(t;s,y)):(t,y)∈DF}はそれぞれ開集合であり、
u(t,X(t;s,y))=Z(t;y)は
ut+b(t,x)⋅∇xu=F(t,x,u),u(s,y)=u0(y)と(t,x,u(t,x))∈Gを満たすΩF上の一意なC1実数値解を定める。特に、各(s,y0)、y0∈Aの近傍で解が存在し、その点における初期値問題の解の芽は一意である。
証明.G上の自律場
B(τ,x,z)=(1,b(τ,x),F(τ,x,z))はC1級である。§E13.4 系 4.4により、各(s,y,u0(y))からの局所解は一意である。その第一成分はt、第二成分はX(t;s,y)であり、第三成分はZ(t;y)である。第二成分が第三成分に依存しないことから、拡大系の最大解の存在区間は、軌道X(⋅;s,y)とその上の反応方程式の解Z(⋅;y)が同時に存在する最大区間Jyに等しい。
拡大系の最大流の定義域の開性と§E13.15 系 2.2 (1)により、その流れとC1写像y↦(s,y,u0(y))との合成は開集合DF上でC1級である。特にZ∈C1(DF;R)である。命題 1.2の写像Ψs,AをDFへ制限すると、ΩFへのC1微分同相写像となる。したがってΩFは開集合であり、
u(t,x)=Z(t;X(s;t,x))はC1級である。連鎖律により、u(t,X(t;s,y))=Z(t;y)をtについて微分すると、主張の偏微分方程式を得る。Z(s;y)=u0(y)により初期条件も成り立つ。
同じΩF上の別のC1解uを各軌道へ制限すると、t↦u(t,X(t;s,y))はJy上でZと同じ反応方程式と初期条件を満たす。両解が一致する時刻の集合は、連続性によりJyの閉集合であり、§E13.4 系 4.4の局所一意性によりJyの開集合でもある。この集合はsを含み、Jyは区間であるから、両解はJy全体で一致する。よってu=uである。任意の局所解の定義域は、(s,y0)のまわりで短い特性曲線の族からなる近傍を含む。その近傍でも同じ反応方程式の一意性が成り立つから、解の芽は一意である。▨
系 2.2.n∈N≥1、I=(α,β)⊂Rを有界な開区間、s∈Iとする。b∈C1(I×Rn;Rn)、F∈C1(I×Rn×R;R)とし、あるC≥0がすべての(t,x,z)∈I×Rn×Rに対して
∣b(t,x)∣≤C(1+∣x∣),∣F(t,x,z)∣≤C(1+∣z∣)を満たすとする。このとき、任意のu0∈C1(Rn;R)に対して、
ut+b⋅∇xu=F(t,x,u),u(s,x)=u0(x)はI×Rn上で一意なC1実数値解をもつ。
証明.(x,z)∈Rn+1にノルム∥(x,z)∥=∣x∣+∣z∣を入れると、拡大した非自律場H(t,x,z)=(b(t,x),F(t,x,z))は
∥H(t,x,z)∥≤2C+C∥(x,z)∥を満たす。HはC1級であるから、各コンパクトな時間区間と閉球の積上で連続かつ状態変数について一様に Lipschitz 連続である。§E13.7 定理 4.1をp=2C、q=Cとして適用すると、すべての(y,u0(y))からの拡大系の解はI全体に存在する。得られた絶対連続解は、連続な右辺をもつ積分方程式を満たすためC1級である。速度場bについても同じ定理を正負両方向の時間に適用することによって、各(t,x)∈I×Rnから時刻sへ至る軌道が存在する。したがって定理 2.1の存在領域はΩF=I×Rnである。同定理により、I×Rn上のC1解は一意に存在する。▨
例 2.3.n∈N≥1、a>0とし、時空領域をR×Rnとする。b=0、F(t,x,z)=z2、u0(x)=aとする。このとき、各軌道はX(t;0,y)=yであり、反応方程式はZ′=Z2、Z(0)=aである。
Z(t)=1−ata(t<1/a)は初期条件を満たし、微分するとZ′=a2/(1−at)2=Z2である。よって定理 2.1により、解はu(t,x)=a/(1−at)である。t↑1/aでu(t,x)→+∞となるため、時刻1/aを含む近傍へ連続な解として延長することはできない。係数はC∞級であるが、反応z2は系 2.2の線形増大条件を満たさない。
同じ反応でn=1、u0(x)=exとすると、最大存在領域と解は
ΩF={(t,x):t<e−x},u(t,x)=1−texexである。任意のε>0に対してe−x<εとなるxが存在するため、[0,ε)×R全体に解を延長することはできない。したがって、各初期点の近傍での局所存在から、非コンパクトな初期断面全体に共通する存在時間は従わない。
例 2.4.n=1、b(t,x)=x2とする。このとき、y>0に対して、
X(t;0,y)=1−tyy(t<1/y)はX′=X2、X(0;0,y)=yを満たし、t↑1/yでX(t;0,y)→+∞となる。したがって、この軌道の最大存在区間は(−∞,1/y)である。y<0では同じ式の最大存在区間は(1/y,∞)、y=0ではX(t;0,0)=0が全時刻に存在する。
F=0ならば軌道上の反応方程式はZ′=0であるが、y>0から出る特性曲線の寿命は有限である。速度場x2は線形増大条件を満たさない。反応が全時刻に存在しても、軌道が有限時刻で無限遠へ離脱すると、その軌道による初期値の輸送はその時刻で終わる。
3 準線形方程式の特性曲線
例 3.1.A⊂Rを開区間、u0∈C1(A;R)とし、方程式
ut+uux=0,u(0,x)=u0(x)を考える。このとき、C1解uの特性曲線x=x(t)がx′=u(t,x(t))を満たすならば、連鎖律により
dtdu(t,x(t))=ut(t,x(t))+u(t,x(t))ux(t,x(t))=0である。したがって、x(0)=ξから出る曲線上ではu=u0(ξ)であり、x(t)=ξ+tu0(ξ)となる。
写像Q(t,ξ)=(t,ξ+tu0(ξ))の微分の行列式は1+tu0′(ξ)である。1+t0u0′(ξ0)=0ならば、§E4.7 定理 2.1により(t0,ξ0)の近傍で逆写像(t,x)↦(t,ξ(t,x))が存在する。
u(t,x)=u0(ξ(t,x))とおくと、x=ξ+tu0(ξ)の微分から
ξx=1+tu0′(ξ)1,ξt=−1+tu0′(ξ)u0(ξ)を得る。よって
ux=1+tu0′(ξ)u0′(ξ),ut=−1+tu0′(ξ)u0(ξ)u0′(ξ)であり、ut+uux=0が成り立つ。t0=0の場合、この局所解は初期条件u(0,x)=u0(x)も満たす。
A=R、u0(ξ)=−ξならば、x(t)=(1−t)ξであり、t<1で
u(t,x)=−1−tx,ux(t,x)=−1−t1となる。すべての曲線は時刻t=1にx=0へ集まり、t↑1でuxは非有界となる。したがって、初期値からのこの解を(1,0)の近傍へC1解として延長することはできない。各曲線上の値−ξも、交点における一つの連続な値と両立しない。
方程式ut+uux=0の速度は未知解uである。t<1で定まる速度場b(t,x)=−x/(1−t)はt=1へC1級で延長されないため、曲線の集積は注意 1.5の仮定に反しない。
4 非特性初期曲面
定義 4.1 (非特性条件).n∈N≥1、V⊂Rnを開集合、a∈C1(V;Rn)、Σ⊂VをC1超曲面とする。p∈Σにおいて
a(p)∈/TpΣが成り立つとき、Σはpで作用素a⋅∇に対して 非特性 (noncharacteristic) であるという。
定理 4.3.n∈N≥1、V⊂Rnを開集合、a∈C1(V;Rn)、F∈C1(V×R;R)とし、Σ⊂VをC1超曲面、g∈C1(Σ;R)とする。このとき、p∈Σにおいてa(p)∈/TpΣならば、pの開近傍W⊂Vとu∈C1(W;R)が存在し、
a(x)⋅∇u(x)=F(x,u(x))(x∈W),u∣W∩Σ=g∣W∩Σを満たす。この初期値問題の解の芽はpにおいて一意である。
証明.pの開近傍V0上でΣ∩V0={x∈V0:ρ(x)=0}とDρ(x)=0(x∈V0)を満たすC1定義関数ρをとる。非特性条件はDρ(p)a(p)=0と同値である。V0を小さくすると、Dρ(x)a(x)はV0上で一定の符号をもち、零にならない。pのまわりのΣのC1局所パラメータ表示をσ:O→Σ∩V0とし、O⊂Rn−1は開集合、σ(y0)=pとする。aの局所流をΦrとおき、
K(r,y)=Φr(σ(y))と定める。§E13.15 系 2.2 (1)によりKは(0,y0)の近傍でC1級である。その微分の列は
a(p), ∂y1σ(y0),…,∂yn−1σ(y0)である。後ろのn−1本はTpΣの基底であるから、非特性条件によりDK(0,y0)は可逆である。§E4.7 定理 2.1により、十分小さいε>0とy0の開近傍O0⊂Oを選ぶと、Kは(−ε,ε)×O0からその像へのC1微分同相写像となる。
拡大した自律場(x,z)↦(a(x),F(x,z))の、初期値(σ(y),g(σ(y)))から出る流れの第二成分をZ(r,y)とする。§E13.4 系 4.4と§E13.15 系 2.2 (1)により、εとO0を小さくするとZは(−ε,ε)×O0上でC1級であり、
∂rZ(r,y)=F(K(r,y),Z(r,y)),Z(0,y)=g(σ(y))を満たす。
εとO0を小さくしてW=K((−ε,ε)×O0)⊂V0とし、u=Z∘K−1と定める。ρ(K(0,y))=0であり、∂r(ρ(K(r,y)))=Dρ(K(r,y))a(K(r,y))は一定の符号をもつため、K(r,y)∈Σであることとr=0であることは同値である。連鎖律により、
∂rZ(r,y)=∇u(K(r,y))⋅∂rK(r,y)=a(K(r,y))⋅∇u(K(r,y))であるから、uは方程式を満たす。r=0ではu(σ(y))=Z(0,y)=g(σ(y))である。
同じ初期値をもつ別の局所解uに対し、pの近傍を小さくして各短い軌道r↦K(r,y)を両解の定義域に含める。u(K(r,y))はZ(r,y)と同じ反応方程式と初期条件を満たす。Fの第二変数に関する局所 Lipschitz 連続性と§E13.4 系 4.4により両者は一致するため、uとuはpの近傍で一致する。▨
例 4.4.R2の座標を(x,y)とし、方程式∂xu=0と初期曲線Σ={(x,0):x∈R}を考える。このとき、速度場はa=(1,0)であり、aは各点でTΣに属するため、初期曲線は非特性条件を満たさない。
C1解がu(x,0)=g(x)を満たすならば、初期曲線に沿って微分することによってg′(x)=∂xu(x,0)=0を得る。よって、非定数のg∈C1(R;R)に対して初期曲線全体の近傍で解は存在しない。特にg(x)=xならば、初期曲線のどの点の近傍にも解は存在しない。
一方、g=0の場合、任意のh∈C1(R;R)でh(0)=0を満たすものに対してu(x,y)=h(y)は方程式と初期条件を満たす。h(y)=λy、λ∈Rとすれば、任意の初期点で互いに異なる無数の解の芽を得る。したがって、非特性条件がない場合には、存在と一意性のいずれも一般には成り立たない。