物理を解説 ♪
ツイッター用のシェアボタン フェイスブック用のシェアボタン はてなブックマーク用のシェアボタン ライン用のシェアボタン
ツイッター用のシェアボタン フェイスブック用のシェアボタン はてなブックマーク用のシェアボタン ライン用のシェアボタン

ラグランジュ微分

こんな簡単な話をこんな後の方に持ってきたわけ。
作成:2020/8/16

なぜ混乱しやすいのか

流体力学には 2 つの流儀が存在しているのだと以前に話したことがある.オイラーの方法とラグランジュの方法だ.オイラーの方法は流体の存在する空間全体を目の前に置いて流れの分布を眺めるようなイメージであり,ラグランジュの方法は流体の一部分だけに注目してその動きを追いかけていくようなイメージである.現在の主流はオイラーの方法であるから,あまり読者を混乱させないようにということで,ここまでオイラーの方法だけで話を進めてきたのだった.

今回説明する「ラグランジュ微分」というのは両方の流儀を混ぜたような概念である.

時々刻々と移動していく流体の「ある一部分」を追いかけながら,その「一部分」が持つ物理量Aの時間的な変化を考える.それはいかにもラグランジュ的な物の見方なのだが,それをオイラーの方法での数学設定を使って表したのがラグランジュ微分なのである.

「思想はラグランジュ,舞台設定はオイラー」という感じなのだ.だからこれはラグランジュの方法とは言えない.あくまでもオイラーの方法の枠組みの中での話になる.この辺りがまた初学者を混乱させてしまうのである.もちろん私も混乱していた.


考え方

話は簡単である.時刻t_{0}に位置(x_{0},y_{0},z_{0})にある流体の持つ何らかの物理量をA(x_{0},y_{0},z_{0},t_{0})と表すことにする.この地点にある流体が持つ速度は 数式 だが,面倒なので以後は変数を省いて単にv_x,v_y,v_zと表そう.Δt秒後にはこの流体はおよそ 数式 付近へと移動していることであろう.つまり,Δt秒が経過した後のこの移動後の地点での物理量Aは次のように表される. 数式 これを多変数関数のテイラー展開で近似的に表せば次のようになる. 数式 さて,この流体と一緒に移動している人にとって物理量AΔtの間にどれだけ変化したかといえば,今の結果を使って, 数式 のように書けるが,これを時間経過による変化の割合で表すと次のようになる. 数式 (1) 式の近似にはΔtの 2 乗以上に比例する程度の誤差が含まれていた.それをΔtで割った (2) 式にはΔtあるいはΔtの 2 乗以上に比例するような誤差が含まれているが,ここでΔt→0の極限を考えれば無視できるようになるので,それを等式を使って次のように表すことにする. 数式 これが「ラグランジュ微分」「物質微分」「実質微分」などと呼ばれるものの定義である.

先ほども言ったがもう一度言っておこう.これは流体と一緒に流れている人から見た,その人のいる地点での物理量\color{red}{A}の時間的な変化率である.

物理量Aを省いて,順序も次のように少し変えて表記した方が使いやすい. 数式 右辺の最初の 3 つの項は速度ベクトルvとベクトル演算子\nabla = ( ∂/∂x , ∂/∂y , ∂/∂z )との内積の形になっているので,次のように表記することもある. 数式 この物理量Aとして流体の速度を入れることもできるので,例えばv_xを入れたときには 数式 という計算をすればいいことになる.



別の導出法があるのでは?

ここまでの記事では徹底的にオイラーの方法に従うことを優先していたので,いかにもオイラーの方法らしい考え方で「連続の方程式」や「流体の運動方程式」を導出してきたのだった.しかしラグランジュ微分という,こんなに便利な概念があるのなら,それを使って導出することもできたのではないだろうか

例えば,「連続の方程式」は流体と一緒に流れている人にとっての「ごく周辺の質量」が変化しないことをラグランジュ微分で表せば導出できるのではないだろうか

「流体の運動方程式」の方だってラグランジュ微分で導出できそうである.前回やったように,流体の運動方程式の左辺はラグランジュ微分に書き換えて 1 つの項で表すことができる.流体と一緒に流れている人にとっての普通の運動方程式として解釈してやることができるだろう.逆に言えば,そのようなイメージで導出することもできたのではないだろうか

しかしそれをやってしまうとオイラーの方法で説明をしているのかラグランジュの方法で説明をしているのか,わけが分からなくなってしまうのではないかという心配があった.そのような導出法はラグランジュ微分と同じように,「思想はラグランジュ,舞台設定はオイラー」という感じの中途半端な理論構造になっているのである.

私自身,そういう導出法は「ラグランジュの方法」に属するものだと勘違いしていた.

この他に「思想も舞台設定もラグランジュ」という導出方法もちゃんとある.そちらは見た目からして違ったものになるので「ラグランジュの連続方程式」や「ラグランジュの運動方程式」と呼ばれて区別される.そこから変数を変換してやればすでに導出した式と同じものが得られる.それらについては必要があれば別の機会に記事を書くかもしれない.

このように,流体力学には幾つかの入り口があって名前も紛らわしい.まだ良くわかっていない内に複数のやり方を次々と見せられると,一体どれがどういう考えで導かれたもので,本当に必要なのはどれで,それぞれの関係はどうなっているかというのが把握しきれず,整理されていない情報量の多さに圧倒されてしまうのである.


連続の方程式をラグランジュ微分で導出してみる

ラグランジュ微分を使った導出法というものがどういうものなのか,言ったからには見せておこう.これからやるのは連続の方程式の導出だけである.

流体と一緒に流れている観測者のごく周辺にある質量は密度ρと微小体積ΔVの積で表される.それが流れて行っても変化しないというのだから,この状況はラグランジュ微分を使って次のように表される. 数式 ラグランジュ微分の中身は偏微分であるから,普通の「関数の積の微分」のルールがそのまま使える. 数式 両辺をρΔVで割れば,第 1 項にρ,第 2 項にΔVを集めることが出来て,分かりやすくなるかもしれない. 数式 あまり分かりやすくもないか…….左辺第 2 項のΔVを何とか変形して,すでに知っている連続の方程式の形に近付けたい.微小体積をΔV=ΔxΔyΔzとして分解してやると左辺第 2 項は次のように変形できる.ここでも「積の微分」のルールを使うのである. 数式 形はあまり変わらないままだ.ちょっと物理的な解釈を挟んで何とかしないといけないだろう.例えばΔxというのは微小な幅を意味するのだから,ごく近くにある 2 つの点の座標値の差である.例えば,Δx=x_{2}-x_{1}と置いて考えたらいいだろう.しかしx座標のラグランジュ微分というのは実行可能なのだろうか座標xだって物理量なのだから,上の説明で使った物理量Aの代わりにxと置いてもいいのだろうが,ちゃんと意味を成すのか考えてみよう.とりあえずラグランジュ微分の定義に従って機械的に計算を実行してみると次のようになる. 数式 この計算を見て,途中に出てくる∂x/∂tは 0 ではなくv_xなのではないかと怪しむ人がいるかもしれない.基本に返って考えてみよう.物理量Aというのは変数(x,y,z,t)によって決まる量なのだった.Aの代わりに考えているxというのはこれと同様に変数(x,y,z,t)の組み合わせによって決まる量なのだが,今はたまたまその中のxのみによって決まる物理量であって,tとは無関係である.それで∂x/∂t =0となるのである.

(5) 式の意味は,位置xにあった流体から目を離さず追いかけていった時に,Δt秒後にはその流体のx座標は変化しているが,その変化の割合はどれくらいなのかというのを計算したことになる.そしてそれはv_xだという答えが出ているのだから,なるほど,確かにその通りだ.問題なく使えそうである.

確認は済んだので (4) 式に戻ってその第 1 項にあるラグランジュ微分を計算してみよう. 数式 (5) 式ではただのv_xという答えだったが,それは流体の初期位置xでの速度を意味しているのだった.今はx_{1}x_{2}は初期位置が違っているのでそれぞれの速度には違いがあることをちゃんと意識して答えを出したのである.そしてその速度差をΔv_xと書き表すことにした.

(4) 式の第 2 項,第 3 項も同じように計算できて,結局 (4) 式は次のように書き表される. 数式 これはΔx,Δy,Δzをそれぞれ 0 に近付ける極限を考えればdiv vと同じものになる.つまり,最初に仮定した微小体積ΔVとして無限小の領域を考えればそうなるということである.ということは (3) 式は次のように書き直せる. 数式 この第 1 項のラグランジュ微分を定義に戻って書き直すと 数式 となって,連続の方程式が再現されるのである.

色々な考察が必要で,なかなか大変である.私がいきなりこういうところから入るのを避けた理由が分かってもらえるだろうか.



趣味の物理学書店

趣味で量子力学2の広告バナー