テスラノートテスラノート物理シミュレーションとAIとその応用。質量保存則を考えた時と同様に、流体の運動量保存則とはどのようなものになるかを考えていきます。しても上のような座標と微小体積を考えます。この微小体積における運動量保存則を考えていきます。ここで、微小体積中の流体の運動量が変化するとしたら、どのような量が変化するかを考えます。\\\\\end{align}$$となっているはずです。これは大丈夫でしょうか?。質量×速度=運動量の時間変化を考えているだけです。それは何によって変化するでしょうか?まず考えられるのは、「流速によって」です。そこで、前回と同じように\(x_1\)方向に流入と流出がある場合を考えてみます。微小体積の左の面から密度\(\rho\)で速度\(u_1\)で入ってくるとします。x_1}\)になるかと思います。u_1\)という運動量が流速によってどう変化する考えてみましょう。単純に\(x_1\)方向で加速、減速する場合。これは普通に、先ほどの質量と全く同じように考えてやれば大丈夫な気がします。\(x_1\)成分による運動量変化を\(M_{1x_1}\)とします。これらの単位はkg・m/s/sで、運動量の時間変化を考えています。\label{eq:scale-factor-2}\end{align}$$こんな風になるかと思います。今までなあなあにしてきましたがちゃんと微小体積の値について考えます。基本的に微小体積の値を代表しているのは一番左下の点、\(x_1,x_2,x_3\)が一番小さいところです。値が入っているのは左下の点ですが、左下の3つの面は左下の点の値で流入してくると考えます。x_3\)離れた点の一番右上の点に代表して値が入っていると考え、右上の3つの面はその値であると考えます。流入、流出という形にしていいのかと思うかもしれませんが、単に正方向を定義しているだけで、流入がマイナスになればすなわち流出なので、大丈夫です。他の方向からの成分で曲げられたりして変わる運動量はどのようになるのでしょうか?なかなか難しいところです。したら間違った理解をしているかもしれませんので、あればご指摘ください。現象としては、このような場合どのような現象なのか考えてみます。現象としては、別の方向からの速度を受けて、その方向に加速・減速するような現象になるはずです。元の流速との足し算になるような形になりそうではないでしょうか?今、下図のように\(x_1\)方向の流速に、\(x_3\)方向の流速があって、曲げられるような現象が起きるとします。流速はおそらく、その2つのベクトルを足したような流速になるでしょう。この場合、流速は速くなっています。運動量は変化しています。いやいや、\(x_1\)成分は変わってないから、\(x_1\)方向の運動量は変化してないじゃんと思ったかた、そうですよね。この辺がめんどくさいところです。x_3\)という運動量が時間的にどう変化するかです。u_3\)に運ばれて変化していて、変化した運動量との保存則を考えねばなりません。すごくややこしいですね。x_3\)の変化、M_{1x_2}はこうなります。\label{eq:scale-factor-3}\end{align}$$下のように、こうじゃないの?と思った方いるかもしれません。この辺がややこしいところです。いい説明あればコメントください。\label{eq:scale-factor-4}\end{align}$$だんだん見えてきました。一応、同じように、\(x_3\)成分もやっておきます。\label{eq:scale-factor-5}\end{align}$$これでだいたいわかりました。全体としてどうなるかと見ておきます。\label{eq:scale-factor-6}\end{align}$$よって、\(u_1\)のやつから算出してみましょう。\label{eq:scale-factor-7}\end{align}$$できました。添字にjを使ったのには、理由があります。いうか、i以外ならなんでもいいです。ここの添字はどうせ縮約を流のでなんでもいいですが、今は\(u_i\)等を考えているので、ここに抵触してはいけません。どこと縮約をとっていいかわからなくなります。他のも同時に書いてみます。\label{eq:scale-factor-8}\end{align}$$おお、なんかこれは、、、さらにテンソル表記でまとめられそうですね。\label{eq:scale-factor-10}\end{align}$$綺麗ですねー。これを移流による流体の運動量保存則です。運動量保存則はこれでいいのでしょうか?そんなわけありませんね。他にも、運動量に影響を与える要因があるはずです。考えられるものはなんでしょう?重力は?粘性とかの摩擦は?圧力は?次回からはこの辺を考えていきます。計算アルゴリズム-格子について前回、MAC法で方程式を導きました。運動量保存則が完成しました。プログラムは、以下の書籍の内容やプログラムを参考にして作ってます。コメントをキャンセルメールアドレスが公開されることはありません。