§E15.34輸送方程式と特性曲線

最終更新

輸送方程式は、所与の速度場による物質や量の移動を記述する、一階偏微分方程式の基本的な例である。例えば、一定の速度で移動する濃度分布は、反応や外部からの供給がなければ、初期の形を保って平行移動する。速度が場所や時刻によって変わる場合には、初期分布の各点がどの経路をたどるかを求める必要がある。この経路を時空の曲線として表したものが特性曲線であり、特性曲線に沿う解の変化は常微分方程式に従う。特性曲線による解の構成は、線形・半線形輸送方程式を常微分方程式の理論と結び付け、初期値を与える曲面が流れを横切るという非特性条件の意味を明らかにする。本記事では、特性曲線による古典解の構成と初期値問題の基本的な性質を解説する。

1 速度場の流れと線形輸送方程式

定義 1.1 (特性曲線).n∈N≥1n\in\NN、U⊂R×RnU\subset\R\times\R^nを開集合、b∈C1(U;Rn)b\in C^1(U;\R^n)とする。(s,y)∈U(s,y)\in Uに対して、初期値問題

dXdt(t;s,y)=b(t,X(t;s,y)),X(s;s,y)=y\frac{dX}{dt}(t;s,y)=b(t,X(t;s,y)),\qquad X(s;s,y)=y

の最大解をX(⋅;s,y)X(\cdot;s,y)、その定義区間をIs,yI_{s,y}と書く。時空の曲線t↦(t,X(t;s,y))t\mapsto(t,X(t;s,y))を速度場bbの 特性曲線 (characteristic curve) という。

命題 1.2.n∈N≥1n\in\NN、U⊂R×RnU\subset\R\times\R^nを開集合、b∈C1(U;Rn)b\in C^1(U;\R^n)とし、X(⋅;s,y)X(\cdot;s,y)を速度場bbの最大解とする。s∈Rs\in\Rを固定し、A⊂{y∈Rn:(s,y)∈U}A\subset\{y\in\R^n:(s,y)\in U\}を開集合とする。

Ds,A={(t,y):y∈A, t∈Is,y},Ψs,A(t,y)=(t,X(t;s,y))D_{s,A}=\{(t,y):y\in A,\ t\in I_{s,y}\},\qquad \Psi_{s,A}(t,y)=(t,X(t;s,y))

とおく。このとき、Ds,AD_{s,A}は開集合であり、Ψs,A\Psi_{s,A}は開集合Ωs,A=Ψs,A(Ds,A)⊂U\Omega_{s,A}=\Psi_{s,A}(D_{s,A})\subset UへのC1C^1微分同相写像である。逆写像は

Ψs,A−1(t,x)=(t,X(s;t,x))\Psi_{s,A}^{-1}(t,x)=(t,X(s;t,x))

である。特に、時刻ttを固定すると、y↦X(t;s,y)y\mapsto X(t;s,y)は{y∈A:t∈Is,y}\{y\in A:t\in I_{s,y}\}からその像へのC1C^1微分同相写像である。

証明.UU上の自律場B(τ,x)=(1,b(τ,x))B(\tau,x)=(1,b(\tau,x))の最大流をφ\varphiとする。第一成分の微分は11であるから、

φr(s,y)=(s+r,X(s+r;s,y))\varphi_r(s,y)=(s+r,X(s+r;s,y))

である。最大流の定義域は開集合であり、§E13.15 系 2.2 (1)によりφ\varphiはC1C^1級である。したがって、Ds,AD_{s,A}は開集合であり、(t,y)↦X(t;s,y)(t,y)\mapsto X(t;s,y)はC1C^1級である。

§E13.15 定理 1.3により、t∈Is,yt\in I_{s,y}ならばs∈It,X(t;s,y)s\in I_{t,X(t;s,y)}であり、

X(s;t,X(t;s,y))=yX(s;t,X(t;s,y))=y

が成り立つ。よってΨs,A\Psi_{s,A}は単射である。時刻ssの断面から時刻ttの断面への写像は、§E13.15 系 2.2 (3)の微分同相写像とその逆を、それぞれの断面へ制限したものである。両断面の接空間はともに{0}×Rn\{0\}\times\R^nであるから、DyX(t;s,y)D_yX(t;s,y)は可逆である。したがって、

DΨs,A(t,y)=(10b(t,X(t;s,y))DyX(t;s,y))D\Psi_{s,A}(t,y)= \begin{pmatrix} 1&0\\ b(t,X(t;s,y))&D_yX(t;s,y) \end{pmatrix}

も可逆である。§E4.7 定理 2.1によりΨs,A\Psi_{s,A}は各点の近傍でC1C^1微分同相写像となる。その像Ωs,A\Omega_{s,A}は開集合であり、単射性によって各局所逆写像は一つのC1C^1逆写像を定める。合成則で得た等式から、この逆写像は(t,x)↦(t,X(s;t,x))(t,x)\mapsto(t,X(s;t,x))である。▨

定理 1.3.n∈N≥1n\in\NN、U⊂R×RnU\subset\R\times\R^nを開集合とし、b∈C1(U;Rn)b\in C^1(U;\R^n)、c,f∈C1(U;R)c,f\in C^1(U;\R)とする。s∈Rs\in\R、開集合A⊂{y:(s,y)∈U}A\subset\{y:(s,y)\in U\}、u0∈C1(A;R)u_0\in C^1(A;\R)に対し、速度場bbの流れからDs,AD_{s,A}、Ωs,A\Omega_{s,A}を命題 1.2のように定める。このとき、

ut+b(t,x)⋅∇xu+c(t,x)u=f(t,x),u(s,y)=u0(y)u_t+b(t,x)\cdot\nabla_xu+c(t,x)u=f(t,x),\qquad u(s,y)=u_0(y)

を満たすu∈C1(Ωs,A;R)u\in C^1(\Omega_{s,A};\R)は一意に存在する。C(t,y)=∫stc(τ,X(τ;s,y)) dτC(t,y)=\int_s^t c(\tau,X(\tau;s,y))\,d\tauとおくと、その解は

u(t,X(t;s,y))=e−C(t,y)[u0(y)+∫steC(τ,y)f(τ,X(τ;s,y)) dτ]((t,y)∈Ds,A)u(t,X(t;s,y))= e^{-C(t,y)}\left[u_0(y)+\int_s^t e^{C(\tau,y)}f(\tau,X(\tau;s,y))\,d\tau\right] \qquad ((t,y)\in D_{s,A})

で与えられる。

証明.(t,y)∈Ds,A(t,y)\in D_{s,A}に対して、主張の表示式の右辺をv(t,y)v(t,y)とする。Is,yI_{s,y}はssを含む区間であるから、表示式の積分に現れる軌道はすべて定義されている。流れと係数はC1C^1級である。各(t,y)(t,y)の積分区間を含む軌道のコンパクト部分を覆う近傍をとると、近傍内の初期点に対しても同じ区間上で被積分関数とその初期点微分が連続である。したがって、積分を初期点について微分することによってC,v∈C1(Ds,A;R)C,v\in C^1(D_{s,A};\R)を得る。

各y∈Ay\in Aについて、§E13.3 定理 3.2によりt↦v(t,y)t\mapsto v(t,y)はIs,yI_{s,y}上の方程式

∂tv(t,y)+c(t,X(t;s,y))v(t,y)=f(t,X(t;s,y)),v(s,y)=u0(y)\partial_tv(t,y)+c(t,X(t;s,y))v(t,y)=f(t,X(t;s,y)),\qquad v(s,y)=u_0(y)

の一意な古典解である。u=v∘Ψs,A−1u=v\circ\Psi_{s,A}^{-1}と定めると、命題 1.2によりu∈C1(Ωs,A;R)u\in C^1(\Omega_{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))\partial_tv(t,y)=u_t(t,X(t;s,y))+b(t,X(t;s,y))\cdot\nabla_xu(t,X(t;s,y))

であるから、uuは所与の偏微分方程式と初期条件を満たす。

u~∈C1(Ωs,A;R)\widetilde u\in C^1(\Omega_{s,A};\R)も同じ偏微分方程式と初期条件を満たすとする。連鎖律によりt↦u~(t,X(t;s,y))t\mapsto\widetilde u(t,X(t;s,y))は同じ一次線形方程式と初期条件を満たす。積分因子による一意性からu~∘Ψs,A=v=u∘Ψs,A\widetilde u\circ\Psi_{s,A}=v=u\circ\Psi_{s,A}であり、Ψs,A\Psi_{s,A}の全射性からu~=u\widetilde u=uである。▨

例 1.4.n∈N≥1n\in\NN、v∈Rnv\in\R^n、γ∈R\gamma\in\R、s∈Rs\in\R、u0∈C1(Rn;R)u_0\in C^1(\R^n;\R)とする。このとき、速度場b(t,x)=vb(t,x)=vの流れはX(t;s,y)=y+(t−s)vX(t;s,y)=y+(t-s)vである。したがって、定理 1.3により、

ut+v⋅∇xu+γu=0,u(s,x)=u0(x)u_t+v\cdot\nabla_xu+\gamma u=0,\qquad u(s,x)=u_0(x)

のR×Rn\R\times\R^n上の解は

u(t,x)=e−γ(t−s)u0(x−(t−s)v)u(t,x)=e^{-\gamma(t-s)}u_0(x-(t-s)v)

である。γ=0\gamma=0のとき、解は各特性曲線上で初期値を保つ。

注意 1.5.n∈N≥1n\in\NN、U⊂R×RnU\subset\R\times\R^nを開集合、b∈C1(U;Rn)b\in C^1(U;\R^n)とする。同じ初期時刻ssの二つの初期点y1,y2y_1,y_2から出る特性曲線が、共通に存在する時刻ttでX(t;s,y1)=X(t;s,y2)X(t;s,y_1)=X(t;s,y_2)を満たすならば、命題 1.2の逆流を適用してy1=y2y_1=y_2を得る。したがって、異なる初期点から出る特性曲線は時空で交わらない。

2 半線形輸送方程式と存在領域

定理 2.1.n∈N≥1n\in\NN、U⊂R×RnU\subset\R\times\R^nを開集合、b∈C1(U;Rn)b\in C^1(U;\R^n)とする。G⊂U×RG\subset U\times\Rを開集合、F∈C1(G;R)F\in C^1(G;\R)とする。s∈Rs\in\R、開集合A⊂{y:(s,y)∈U}A\subset\{y:(s,y)\in U\}、u0∈C1(A;R)u_0\in C^1(A;\R)は、すべてのy∈Ay\in Aについて(s,y,u0(y))∈G(s,y,u_0(y))\in Gを満たすとし、X(⋅;s,y)X(\cdot;s,y)を速度場bbの最大解とする。各y∈Ay\in Aに対して、

Z′(t;y)=F(t,X(t;s,y),Z(t;y)),Z(s;y)=u0(y)Z'(t;y)=F(t,X(t;s,y),Z(t;y)),\qquad Z(s;y)=u_0(y)

の最大解をZ(⋅;y)Z(\cdot;y)、その定義区間をJy⊂Is,yJ_y\subset I_{s,y}とする。このとき、

DF={(t,y):y∈A, t∈Jy},ΩF={(t,X(t;s,y)):(t,y)∈DF}D_F=\{(t,y):y\in A,\ t\in J_y\},\qquad \Omega_F=\{(t,X(t;s,y)):(t,y)\in D_F\}

はそれぞれ開集合であり、

u(t,X(t;s,y))=Z(t;y)u(t,X(t;s,y))=Z(t;y)

は

ut+b(t,x)⋅∇xu=F(t,x,u),u(s,y)=u0(y)u_t+b(t,x)\cdot\nabla_xu=F(t,x,u),\qquad u(s,y)=u_0(y)

と(t,x,u(t,x))∈G(t,x,u(t,x))\in Gを満たすΩF\Omega_F上の一意なC1C^1実数値解を定める。特に、各(s,y0)(s,y_0)、y0∈Ay_0\in Aの近傍で解が存在し、その点における初期値問題の解の芽は一意である。

証明.GG上の自律場

B(τ,x,z)=(1,b(τ,x),F(τ,x,z))\mathcal B(\tau,x,z)=(1,b(\tau,x),F(\tau,x,z))

はC1C^1級である。§E13.4 系 4.4により、各(s,y,u0(y))(s,y,u_0(y))からの局所解は一意である。その第一成分はtt、第二成分はX(t;s,y)X(t;s,y)であり、第三成分はZ(t;y)Z(t;y)である。第二成分が第三成分に依存しないことから、拡大系の最大解の存在区間は、軌道X(⋅;s,y)X(\cdot;s,y)とその上の反応方程式の解Z(⋅;y)Z(\cdot;y)が同時に存在する最大区間JyJ_yに等しい。

拡大系の最大流の定義域の開性と§E13.15 系 2.2 (1)により、その流れとC1C^1写像y↦(s,y,u0(y))y\mapsto(s,y,u_0(y))との合成は開集合DFD_F上でC1C^1級である。特にZ∈C1(DF;R)Z\in C^1(D_F;\R)である。命題 1.2の写像Ψs,A\Psi_{s,A}をDFD_Fへ制限すると、ΩF\Omega_FへのC1C^1微分同相写像となる。したがってΩF\Omega_Fは開集合であり、

u(t,x)=Z(t;X(s;t,x))u(t,x)=Z(t;X(s;t,x))

はC1C^1級である。連鎖律により、u(t,X(t;s,y))=Z(t;y)u(t,X(t;s,y))=Z(t;y)をttについて微分すると、主張の偏微分方程式を得る。Z(s;y)=u0(y)Z(s;y)=u_0(y)により初期条件も成り立つ。

同じΩF\Omega_F上の別のC1C^1解u~\widetilde uを各軌道へ制限すると、t↦u~(t,X(t;s,y))t\mapsto\widetilde u(t,X(t;s,y))はJyJ_y上でZZと同じ反応方程式と初期条件を満たす。両解が一致する時刻の集合は、連続性によりJyJ_yの閉集合であり、§E13.4 系 4.4の局所一意性によりJyJ_yの開集合でもある。この集合はssを含み、JyJ_yは区間であるから、両解はJyJ_y全体で一致する。よってu~=u\widetilde u=uである。任意の局所解の定義域は、(s,y0)(s,y_0)のまわりで短い特性曲線の族からなる近傍を含む。その近傍でも同じ反応方程式の一意性が成り立つから、解の芽は一意である。▨

系 2.2.n∈N≥1n\in\NN、I=(α,β)⊂RI=(\alpha,\beta)\subset\Rを有界な開区間、s∈Is\in Iとする。b∈C1(I×Rn;Rn)b\in C^1(I\times\R^n;\R^n)、F∈C1(I×Rn×R;R)F\in C^1(I\times\R^n\times\R;\R)とし、あるC≥0C\geq0がすべての(t,x,z)∈I×Rn×R(t,x,z)\in I\times\R^n\times\Rに対して

∣b(t,x)∣≤C(1+∣x∣),∣F(t,x,z)∣≤C(1+∣z∣)|b(t,x)|\leq C(1+|x|),\qquad |F(t,x,z)|\leq C(1+|z|)

を満たすとする。このとき、任意のu0∈C1(Rn;R)u_0\in C^1(\R^n;\R)に対して、

ut+b⋅∇xu=F(t,x,u),u(s,x)=u0(x)u_t+b\cdot\nabla_xu=F(t,x,u),\qquad u(s,x)=u_0(x)

はI×RnI\times\R^n上で一意なC1C^1実数値解をもつ。

証明.(x,z)∈Rn+1(x,z)\in\R^{n+1}にノルム∥(x,z)∥=∣x∣+∣z∣\|(x,z)\|=|x|+|z|を入れると、拡大した非自律場H(t,x,z)=(b(t,x),F(t,x,z))H(t,x,z)=(b(t,x),F(t,x,z))は

∥H(t,x,z)∥≤2C+C∥(x,z)∥\|H(t,x,z)\|\leq 2C+C\|(x,z)\|

を満たす。HHはC1C^1級であるから、各コンパクトな時間区間と閉球の積上で連続かつ状態変数について一様に Lipschitz 連続である。§E13.7 定理 4.1をp=2Cp=2C、q=Cq=Cとして適用すると、すべての(y,u0(y))(y,u_0(y))からの拡大系の解はII全体に存在する。得られた絶対連続解は、連続な右辺をもつ積分方程式を満たすためC1C^1級である。速度場bbについても同じ定理を正負両方向の時間に適用することによって、各(t,x)∈I×Rn(t,x)\in I\times\R^nから時刻ssへ至る軌道が存在する。したがって定理 2.1の存在領域はΩF=I×Rn\Omega_F=I\times\R^nである。同定理により、I×RnI\times\R^n上のC1C^1解は一意に存在する。▨

例 2.3.n∈N≥1n\in\NN、a>0a>0とし、時空領域をR×Rn\R\times\R^nとする。b=0b=0、F(t,x,z)=z2F(t,x,z)=z^2、u0(x)=au_0(x)=aとする。このとき、各軌道はX(t;0,y)=yX(t;0,y)=yであり、反応方程式はZ′=Z2Z'=Z^2、Z(0)=aZ(0)=aである。

Z(t)=a1−at(t<1/a)Z(t)=\frac{a}{1-at}\qquad(t<1/a)

は初期条件を満たし、微分するとZ′=a2/(1−at)2=Z2Z'=a^2/(1-at)^2=Z^2である。よって定理 2.1により、解はu(t,x)=a/(1−at)u(t,x)=a/(1-at)である。t↑1/at\uparrow1/aでu(t,x)→+∞u(t,x)\to+\inftyとなるため、時刻1/a1/aを含む近傍へ連続な解として延長することはできない。係数はC∞C^\infty級であるが、反応z2z^2は系 2.2の線形増大条件を満たさない。

同じ反応でn=1n=1、u0(x)=exu_0(x)=e^xとすると、最大存在領域と解は

ΩF={(t,x):t<e−x},u(t,x)=ex1−tex\Omega_F=\{(t,x):t<e^{-x}\},\qquad u(t,x)=\frac{e^x}{1-te^x}

である。任意のε>0\varepsilon>0に対してe−x<εe^{-x}<\varepsilonとなるxxが存在するため、[0,ε)×R[0,\varepsilon)\times\R全体に解を延長することはできない。したがって、各初期点の近傍での局所存在から、非コンパクトな初期断面全体に共通する存在時間は従わない。

例 2.4.n=1n=1、b(t,x)=x2b(t,x)=x^2とする。このとき、y>0y>0に対して、

X(t;0,y)=y1−ty(t<1/y)X(t;0,y)=\frac{y}{1-ty}\qquad(t<1/y)

はX′=X2X'=X^2、X(0;0,y)=yX(0;0,y)=yを満たし、t↑1/yt\uparrow1/yでX(t;0,y)→+∞X(t;0,y)\to+\inftyとなる。したがって、この軌道の最大存在区間は(−∞,1/y)(-\infty,1/y)である。y<0y<0では同じ式の最大存在区間は(1/y,∞)(1/y,\infty)、y=0y=0ではX(t;0,0)=0X(t;0,0)=0が全時刻に存在する。

F=0F=0ならば軌道上の反応方程式はZ′=0Z'=0であるが、y>0y>0から出る特性曲線の寿命は有限である。速度場x2x^2は線形増大条件を満たさない。反応が全時刻に存在しても、軌道が有限時刻で無限遠へ離脱すると、その軌道による初期値の輸送はその時刻で終わる。

3 準線形方程式の特性曲線

例 3.1.A⊂RA\subset\Rを開区間、u0∈C1(A;R)u_0\in C^1(A;\R)とし、方程式

ut+uux=0,u(0,x)=u0(x)u_t+u u_x=0,\qquad u(0,x)=u_0(x)

を考える。このとき、C1C^1解uuの特性曲線x=x(t)x=x(t)がx′=u(t,x(t))x'=u(t,x(t))を満たすならば、連鎖律により

ddtu(t,x(t))=ut(t,x(t))+u(t,x(t))ux(t,x(t))=0\frac{d}{dt}u(t,x(t))=u_t(t,x(t))+u(t,x(t))u_x(t,x(t))=0

である。したがって、x(0)=ξx(0)=\xiから出る曲線上ではu=u0(ξ)u=u_0(\xi)であり、x(t)=ξ+tu0(ξ)x(t)=\xi+t u_0(\xi)となる。

写像Q(t,ξ)=(t,ξ+tu0(ξ))Q(t,\xi)=(t,\xi+t u_0(\xi))の微分の行列式は1+tu0′(ξ)1+t u_0'(\xi)である。1+t0u0′(ξ0)≠01+t_0u_0'(\xi_0)\ne0ならば、§E4.7 定理 2.1により(t0,ξ0)(t_0,\xi_0)の近傍で逆写像(t,x)↦(t,ξ(t,x))(t,x)\mapsto(t,\xi(t,x))が存在する。

u(t,x)=u0(ξ(t,x))u(t,x)=u_0(\xi(t,x))

とおくと、x=ξ+tu0(ξ)x=\xi+t u_0(\xi)の微分から

ξx=11+tu0′(ξ),ξt=−u0(ξ)1+tu0′(ξ)\xi_x=\frac{1}{1+t u_0'(\xi)},\qquad \xi_t=-\frac{u_0(\xi)}{1+t u_0'(\xi)}

を得る。よって

ux=u0′(ξ)1+tu0′(ξ),ut=−u0(ξ)u0′(ξ)1+tu0′(ξ)u_x=\frac{u_0'(\xi)}{1+t u_0'(\xi)},\qquad u_t=-\frac{u_0(\xi)u_0'(\xi)}{1+t u_0'(\xi)}

であり、ut+uux=0u_t+u u_x=0が成り立つ。t0=0t_0=0の場合、この局所解は初期条件u(0,x)=u0(x)u(0,x)=u_0(x)も満たす。

A=RA=\R、u0(ξ)=−ξu_0(\xi)=-\xiならば、x(t)=(1−t)ξx(t)=(1-t)\xiであり、t<1t<1で

u(t,x)=−x1−t,ux(t,x)=−11−tu(t,x)=-\frac{x}{1-t},\qquad u_x(t,x)=-\frac{1}{1-t}

となる。すべての曲線は時刻t=1t=1にx=0x=0へ集まり、t↑1t\uparrow1でuxu_xは非有界となる。したがって、初期値からのこの解を(1,0)(1,0)の近傍へC1C^1解として延長することはできない。各曲線上の値−ξ-\xiも、交点における一つの連続な値と両立しない。

方程式ut+uux=0u_t+u u_x=0の速度は未知解uuである。t<1t<1で定まる速度場b(t,x)=−x/(1−t)b(t,x)=-x/(1-t)はt=1t=1へC1C^1級で延長されないため、曲線の集積は注意 1.5の仮定に反しない。

4 非特性初期曲面

定義 4.1 (非特性条件).n∈N≥1n\in\NN、V⊂RnV\subset\R^nを開集合、a∈C1(V;Rn)a\in C^1(V;\R^n)、Σ⊂V\Sigma\subset VをC1C^1超曲面とする。p∈Σp\in\Sigmaにおいて

a(p)∉TpΣa(p)\notin T_p\Sigma

が成り立つとき、Σ\Sigmaはppで作用素a⋅∇a\cdot\nablaに対して 非特性 (noncharacteristic) であるという。

注意 4.2.a∈C∞(V;Rn)a\in C^\infty(V;\R^n)のとき、一階作用素a⋅∇a\cdot\nablaの主表象は2πi a(x)⋅ξ2\pi i\,a(x)\cdot\xiであり、Σ\Sigmaがppで非特性であることは、Σ\Sigmaのppにおける余法線ν∈Rn∖{0}\nu\in\R^n\setminus\{0\}(TpΣT_p\Sigmaと直交する非零ベクトル)が§E15.15 例 4.1の意味の特性方向でないこと、すなわち2πi a(p)⋅ν≠02\pi i\,a(p)\cdot\nu\ne0が成り立つことと同値である。

定理 4.3.n∈N≥1n\in\NN、V⊂RnV\subset\R^nを開集合、a∈C1(V;Rn)a\in C^1(V;\R^n)、F∈C1(V×R;R)F\in C^1(V\times\R;\R)とし、Σ⊂V\Sigma\subset VをC1C^1超曲面、g∈C1(Σ;R)g\in C^1(\Sigma;\R)とする。このとき、p∈Σp\in\Sigmaにおいてa(p)∉TpΣa(p)\notin T_p\Sigmaならば、ppの開近傍W⊂VW\subset Vとu∈C1(W;R)u\in C^1(W;\R)が存在し、

a(x)⋅∇u(x)=F(x,u(x))(x∈W),u∣W∩Σ=g∣W∩Σa(x)\cdot\nabla u(x)=F(x,u(x))\quad(x\in W),\qquad u|_{W\cap\Sigma}=g|_{W\cap\Sigma}

を満たす。この初期値問題の解の芽はppにおいて一意である。

証明.ppの開近傍V0V_0上でΣ∩V0={x∈V0:ρ(x)=0}\Sigma\cap V_0=\{x\in V_0:\rho(x)=0\}とDρ(x)≠0D\rho(x)\ne0(x∈V0x\in V_0)を満たすC1C^1定義関数ρ\rhoをとる。非特性条件はDρ(p)a(p)≠0D\rho(p)a(p)\ne0と同値である。V0V_0を小さくすると、Dρ(x)a(x)D\rho(x)a(x)はV0V_0上で一定の符号をもち、零にならない。ppのまわりのΣ\SigmaのC1C^1局所パラメータ表示をσ:O→Σ∩V0\sigma:O\to\Sigma\cap V_0とし、O⊂Rn−1O\subset\R^{n-1}は開集合、σ(y0)=p\sigma(y_0)=pとする。aaの局所流をΦr\Phi_rとおき、

K(r,y)=Φr(σ(y))K(r,y)=\Phi_r(\sigma(y))

と定める。§E13.15 系 2.2 (1)によりKKは(0,y0)(0,y_0)の近傍でC1C^1級である。その微分の列は

a(p), ∂y1σ(y0),…,∂yn−1σ(y0)a(p),\ \partial_{y_1}\sigma(y_0),\ldots,\partial_{y_{n-1}}\sigma(y_0)

である。後ろのn−1n-1本はTpΣT_p\Sigmaの基底であるから、非特性条件によりDK(0,y0)DK(0,y_0)は可逆である。§E4.7 定理 2.1により、十分小さいε>0\varepsilon>0とy0y_0の開近傍O0⊂OO_0\subset Oを選ぶと、KKは(−ε,ε)×O0(-\varepsilon,\varepsilon)\times O_0からその像へのC1C^1微分同相写像となる。

拡大した自律場(x,z)↦(a(x),F(x,z))(x,z)\mapsto(a(x),F(x,z))の、初期値(σ(y),g(σ(y)))(\sigma(y),g(\sigma(y)))から出る流れの第二成分をZ(r,y)Z(r,y)とする。§E13.4 系 4.4と§E13.15 系 2.2 (1)により、ε\varepsilonとO0O_0を小さくするとZZは(−ε,ε)×O0(-\varepsilon,\varepsilon)\times O_0上でC1C^1級であり、

∂rZ(r,y)=F(K(r,y),Z(r,y)),Z(0,y)=g(σ(y))\partial_rZ(r,y)=F(K(r,y),Z(r,y)),\qquad Z(0,y)=g(\sigma(y))

を満たす。

ε\varepsilonとO0O_0を小さくしてW=K((−ε,ε)×O0)⊂V0W=K((-\varepsilon,\varepsilon)\times O_0)\subset V_0とし、u=Z∘K−1u=Z\circ K^{-1}と定める。ρ(K(0,y))=0\rho(K(0,y))=0であり、∂r(ρ(K(r,y)))=Dρ(K(r,y))a(K(r,y))\partial_r(\rho(K(r,y)))=D\rho(K(r,y))a(K(r,y))は一定の符号をもつため、K(r,y)∈ΣK(r,y)\in\Sigmaであることとr=0r=0であることは同値である。連鎖律により、

∂rZ(r,y)=∇u(K(r,y))⋅∂rK(r,y)=a(K(r,y))⋅∇u(K(r,y))\partial_rZ(r,y)=\nabla u(K(r,y))\cdot\partial_rK(r,y) =a(K(r,y))\cdot\nabla u(K(r,y))

であるから、uuは方程式を満たす。r=0r=0ではu(σ(y))=Z(0,y)=g(σ(y))u(\sigma(y))=Z(0,y)=g(\sigma(y))である。

同じ初期値をもつ別の局所解u~\widetilde uに対し、ppの近傍を小さくして各短い軌道r↦K(r,y)r\mapsto K(r,y)を両解の定義域に含める。u~(K(r,y))\widetilde u(K(r,y))はZ(r,y)Z(r,y)と同じ反応方程式と初期条件を満たす。FFの第二変数に関する局所 Lipschitz 連続性と§E13.4 系 4.4により両者は一致するため、uuとu~\widetilde uはppの近傍で一致する。▨

例 4.4.R2\R^2の座標を(x,y)(x,y)とし、方程式∂xu=0\partial_xu=0と初期曲線Σ={(x,0):x∈R}\Sigma=\{(x,0):x\in\R\}を考える。このとき、速度場はa=(1,0)a=(1,0)であり、aaは各点でTΣT\Sigmaに属するため、初期曲線は非特性条件を満たさない。

C1C^1解がu(x,0)=g(x)u(x,0)=g(x)を満たすならば、初期曲線に沿って微分することによってg′(x)=∂xu(x,0)=0g'(x)=\partial_xu(x,0)=0を得る。よって、非定数のg∈C1(R;R)g\in C^1(\R;\R)に対して初期曲線全体の近傍で解は存在しない。特にg(x)=xg(x)=xならば、初期曲線のどの点の近傍にも解は存在しない。

一方、g=0g=0の場合、任意のh∈C1(R;R)h\in C^1(\R;\R)でh(0)=0h(0)=0を満たすものに対してu(x,y)=h(y)u(x,y)=h(y)は方程式と初期条件を満たす。h(y)=λyh(y)=\lambda y、λ∈R\lambda\in\Rとすれば、任意の初期点で互いに異なる無数の解の芽を得る。したがって、非特性条件がない場合には、存在と一意性のいずれも一般には成り立たない。

前提記事