力学におけるリュービルの定理の証明
📂動力学力学におけるリュービルの定理の証明
定理
ユークリッド空間 Rn と関数 f:Rn→Rn について次のようなベクトル場が微分方程式で与えられているとする。
x˙=f(x)
このシステムのフロー ϕt(⋅) と領域 D0⊂Rn について Dt:=ϕt(D0) をフローに従って時間 t が経過して移動した領域、そのボリュームを V(t)≡V(Dt) として示す。もし ∇⋅f=0 ならすべての D0⊂Rn と t∈R について次が成り立つ。
V(Dt)=V(D0)
数式の説明
∇⋅f はベクトル場の発散で、ベクトル場が広がったり集まったりする様子を示す。
V はベクトル場で与えられた領域のボリュームを意味する。リウヴィルの定理は言葉より数式の方が簡単だが、簡単に言えば、発散がどこでも 0 ならばフローが移動した領域のボリュームが変わらないということである。
証明
戦略: ベクトル微積分学を少し使えばよい。次のようにベクトル関数と各ベクトルの軸を表せば、解析学と線形代数の基本的なツールによる計算だけで直接演繹される。
f:=(f1,⋯,fn)x:=(x1,⋯,xn)
Part 1. t0=0 のとき dtdVt=t0=∫D0∇⋅fdx
ボリュームの定義: u∈Rn がベクトル関数 f:Rn→Rn によって f(u)=(f1(u),⋯,fn(u)) として変換されるとき、D のボリュームは次のようになる。
V(D)=∫D∂u∂f(u)du1du2⋯dun
ボリュームの定義に従って
V(t)=∫D0det∂x∂ϕt(x)dx
フロー ϕt(x) を t=0 でテイラー展開すると
ϕt(x)=x+f(x)t+O(t2)
x で偏微分すると n×n 単位行列 E に対して
∂x∂ϕt(x)=E+∂x∂ft+O(t2)
のように示すことができる。すると、その行列式は
det∂x∂ϕt(x)====det[E+∂x∂ft]+O(t2)1+tr[∂x∂f]t+O(t2)1+(dx1df1+⋯+dxndfn)t+O(t2)1+t∇⋅f+O(t2)
ここで tr はトレースTraceで、行列の対角成分の合計を意味する。両辺に ∫D0⋅dx を取ると
V(t)===∫D0det∂x∂ϕt(x)dx∫D01dx+∫D0t∇⋅fdx+O(t2)V(D0)+t∫D0∇⋅fdx+O(t2)
V(D0)=V(0) を左辺に移すと
V(t)−V(0)=t∫D0∇⋅fdx+O(t2)
両辺を (t−0) で割ると
t−0V(t)−V(0)=∫D0∇⋅fdx+O(t1)
テイラー近似を t=0 近傍で用いたので、t→0 のとき O(t1)→0 であるから、
dtdVt=0=∫D0∇⋅fdx
Part 2. t=t0 への拡張
y:=ϕt0(x) と仮定すると、Part 1.の議論を同じように進めて次を得る。
dtdVt=t0=∫D0∇⋅fdy
ある定数 c∈R に対して ∇⋅f=c なら、
dtdVt=t0=∫D0∇⋅fdy=∫D0cdy=c∫D0dy
y は x を時間 t0 だけフローに従って変換されたものであるから、∫D0dy は t0 時点のボリュームとなる。これは任意の t0∈R に対しても同様であるから、
V’=cV
上記の微分方程式は自明な解 V(t)=ectV(0) を持つ。
Part 3.
上記 Part 2. の最後の式で c=0 なら V(t)=ectV(0)=V(0) である。
■
参考文献