プログラミングの記事を書きます。熱力学的な内部状態の正体を古典論を元に理解することを目的とします。議論するのは、混合気体ではなく、単一の気体のみで話を進めます。混合気体の場合は、>>ドルトンの分圧の法則を使えば簡単に「混合気体での理想気体の状態方程式」は導出できるのですが、本質を外れそうなので、本記事は単一気体のみとします。(場合によっては質量も)」など・・・・それは実態は何なのでしょう?体積と質量などは感覚的にわかるでしょう。体積イメージ通り閉じた空間のことです。※開いていても良いが、開いていると今考えようとしている何か実態が逃げていきそうなので・・・それは視覚的にもとてもわかりやすいでしょう。空間が膨張すれば気体は膨張しますし、空間が縮小すれば気体は縮小するでしょう。質量これもイメージができるでしょう。理想気体の気体は、その名の通り気体分子なので分子の質量のことを言っています。温度や圧力は?これらがわかりづらいな~って思ったので、この記事を書き始めました。「内部状態の実態は何か?」でもあります。それを気体分子という分子を考え、その運動によって内部状態を記述できるという話をしていきます。ここでは、人が観測しうる圧力や温度とは何か?ということに関心があります。ミクロな視点での興味ではなく、マクロな視点(粒子集団が圧力や温度を生み出す)で見つつ、「マクロな物理量を、気体分子の運動によるもの」と考えて、力学的に考えた時、どのように解釈できるかということに関心ごとがあります。問題設定ですが、下記の絵のように、体積\(L^{3}\)の閉じた空間に\(N\)個の質量\(m\)の気体分子(粒子)が存在しているとしましょう。問題を簡単にするために、1次元(\(x\)方向のみ)で、かつひとつの粒子の運動について考えます。このように、\(x\)方向に\(v_{x}\)の速度をもっていた粒子が壁に当たると、衝撃力のようなものを与えるというものでイメージできるでしょう。それは、粒子の運動量変化として、壁に力を与えます。これが、圧力というものであると解釈できます。具体的には、この衝撃力は次のように考えることができます。(1)\end{align*}と書けます。t\)秒間に何回、この灰色の壁に力を与えるでしょうか?算数での考え方です。粒子が1辺が\(L\)の空間に対して、1往復したときに1回衝突すると考えれば良いです。(2)\end{align*}ですね。わかりにくければ、下記の絵のように箱をつなげればわかりやすいですかね。(3)\end{align*}の衝撃力を与えることになります。ここで、\(x\)向だけを考えましたが、\(x\)方向は何か特別な方向でしたでしょうか?\(x\)方向は、問題を簡単に考えるために選んだ方向であって、\(y\)でも\(z\)でも良かったのです。この場合は方向に特別な意味はなく、等方的だと考えましょう。(4)\end{align*}と書けます。\(x,y,z\)のどの方向かは問題にならないので、\(v_{x}=v_{y}=v_{z}\)としても良いでしょう。(5)\end{align*}とすることができます。(6)\end{align*}となります。先ほどまではひとつの粒子に対して考えていましたが、N個の粒子へ考え方を拡張しましょう。(6)式を、\(N\)倍するだけでしょうか?いや、もう少しちゃんと考えよう。(A)粒子同士の衝突は?ひとつ目が粒子同士の衝突の話です。粒子同士衝突するとは、Aという粒子とBという粒子が同じ位置(粒子の直径分離れているが)に来るときに衝突すると考えるべきでしょうが、そもそも粒子同士に分子間力が働いているので、接近しただけで粒子同士の運動は乱されます。(6)のよりも幾分か小さい値であると予想されます。本記事の冒頭に書きましたが、今は理想気体を考えているですから、(1)分子間力を考慮しない(2)希薄であるという元での理論なので、この(A)による影響は初めから考えないものとしているのです。「(A)粒子同士の衝突(分子間力なども)」の影響は除外できます。(B)すべての粒子が同じ\(v\)で運動するのでしょうか?イメージしてみましょう。社会に出れば活発的な人もいれば、冷めた人もいるでしょう。粒子だって同じです(笑)「あ~運動したくないよ~」ってやつもいれば、「目をぎらぎらさせて日々走り回っている気違いなほどアクティブな人」もいるでしょう。それらを含めて、平均化された集団としての評価で特徴が決められます。速度\(v\)は\(N\)個の粒子を考えるときは、集団としての平均速度\(v\)として考えるのが正しいのです。(7)\end{align*}壁に与える力ここで粒子の運動量変化はわかりましたが、その分壁へ力を与えているはずです。その力を\(F\)とします。t\)秒間に与える力、すなわり力積です。(8)\end{align*}ここで再び協調しておきますが、この\(F\)は一瞬での力を表しているものではありません。t\)秒間での平均化された力を意味しています。t\)秒間で、強い力で壁に激突するときもあれば、弱い力で激突しているときもあるでしょう。t\)秒間平均すると\(F\)であるという意味。(9)\end{align*}と書けます。圧力ようやく(8)式で、壁に与える力というものを数式で表現することができました。初めに知ることができるのは、圧力でしょう。(10)\end{align*}です。(10)\end{align*}と書くことができました。ここで、「おや、\(L^3\)は箱の体積ではないか」ということに気付きます。(11)\end{align*}となります。ここで、「気体粒子の運動」を考えると、圧力の意味が見えてきました。「多数の粒子が動き回って、ポカポカ壁に激突したときの単位面積当たりの力なのか」と。温度ここで、(11)式が何かに見えてこないでしょうか??(\(\diamondsuit\))式の理想気体の状態方程式です。(12)\end{align*}とすればなおさら、(\(\diamondsuit\))と形が似ています。温度というものは、(12)式の右辺で何か表現できそうであることがわかります。「お、運動エネルギー\(\frac{1}{2}m\bar{v}^2\)というのを作ってみよう。(15)\end{align*}と書けるでしょう。(17)\end{align*}と置きましょう。これを[ボルツマン定数]と呼びます。(18)\end{align*}とすることができました。粒子の平均の運動エネルギーと関連付けられる量であることがわかります。温度というものの実態は、「気体分子の運動を考える」ならば運動エネルギーと等価なものだったのです。粒子の平均の運動エネルギーと関連付けられる量このように解釈すれば良いのかと。(17)\end{align*}なのですよね。お気づきでしょうが、この3は(5)式からずっと引きづって来ていますよね。3次元空間だから3なのです。正確には、3自由度だから3なのです。回転運動が入る2原子分子であれば5自由度(並進運動3自由度+回転の2自由度)だったら?5自由度だから5なのです。回転運動のエネルギーが2乗であらわされる場合だけですよ。x^2\)のような場合は?その場合も2乗なので、1自由度とカウントできます。ここで例をひとつ出しましょう。x^2\end{align*}とかけて、この2乗の形が重要なのです。※ここでの\(p\)は圧力ではなく、運動量\(p=mv\)です。座標と運動量で表現されたエネルギー量のことを、特に「ハミルトニアン」と呼びます。解析力学の範疇では、もはや座標\(x\)も運動量\(p\)もそれ自体そのような名前を名付けること自体に意味を成しません。力学的状態を位相空間での運動で考えるとき、座標と運動量という変数に対して正準変換すると数学的には同じ扱いにされてしまうのです。(18)\end{align*}となります。お分かりいただけましたでしょうか?ポテンシャル(運動エネルギー、調和振動子ポテンシャル)が2次関数で書けるとき、その数だけ\(\frac{1}{2}k_{B}T\)が足されるのです。(19)\end{align*}と表現されていきます。これをエネルギー等分配則と言います。温度や圧力は知っているが実態は?温度や圧力の実態が何かを考えるということはしないのに、なぜ温度や圧力という言葉を使って物事を考えることができるのでしょう。別に不思議なことではありません。何も理論的な背景を知らない人が、生まれ持って理論的に考えるなんてことはしません。無意識に観測したデータや数値を、客観的に大きいか小さいかしか理解していないはずです。普通に過ごしていると、「暑くてすごしにくいな~」「寒くて凍えそうだ」と感じるのは、その名の通り感覚器官が感じているだけですよね。それを定量的に気温として観測し、数値化するから「暑い・寒い」の程度が数値的に定義付られます。その温度を、ず~っと身近に感じているとその意味を考えなくなります。「暑いから・・・・こうなる」「冷たいと・・・・・こうなる」と・・・・それは、理論的に理解しているというわけではなく、単に相関関係が頭にしみついているだけです(個人的にそう思っています)。圧力の場合は、力という考え方で説明することができるでしょう。ペンの先で自分の顔に押し付けるととても痛いですが、手のひらで顔を押してもまあ知れているでしょう。両者の押している合計の力は同じですが、感覚的に「ペンの先」で押した方が痛いと感じるのは、単面積当たりにかかっている圧力がとても大きからですよね。実際は、単位面積当たりの力を両者で同じにしたら、手のひらで押そうがそりゃー痛いはずです。(ペンの先より面積が広い分、広範囲で顔面が痛いことになるが)このように圧力が、どのように作用するかは力という考え方を借りて理解することができます。気体が存在している状態で「圧力」とは何を意味しているのでしょう?さらに、温度とはどのように理解すればよろしいでしょう?そこに気体分子があると”考えられる”気体の圧力考えるときに、そこに何があるかを考えるべきでしょう。気体分子がそこにある以上は、気体分子の運動が「圧力」を生みだし、「圧力」や「温度」というものを我々に観測させていると言えるでしょう。正確な表現をしておくことにします。「そこに何があると考えれば様々な問題にロジックの飛躍無しに答えることがことができるでしょう」とした方がいいかも知れません。気体分子が存在しているということを前提にしているけれども、「誰か気体分子が存在していることを目で確認したのでしょうか?」と反論がありそうです。前段で強調しているのです。低温にすると、「気体も液体や固体になるので理想気体ではなくなりますね」などという意味で強調しているわけではなく、理論はその現象を理解し予測さえできれば良いので、気体分子があると考えて物事の話を進めるという程度の意味しか持っていないということを強調しておきたいです。2.7Kにしてもようやく液体になるヘリウムなどはほとんどの温度域で気体なのです。言って、「結構な温度域で気体粒子として扱えますね」とはならないのです。※あえて気体粒子と、”粒子”という言葉を強調しました。粒子と言っているのは、そこに分子という粒々があると考えると1つ・2つ・・・・と数えていけば良いですですが、「では粒々同士の隙間には何も存在しない空間があるのでしょうか?」という点からも、よく考えてみましょう。することで、古典論にも限界があることが見えてくるのではないでしょうか?それは、「それだけの低温域で量子力学的効果が出てこないのか?」という問いにつながるからです。※詳細は省きますが、物理学では「波」と「粒子」は定義が分かれています。上記の低温域(とてもミクロな現象など)では、古典的な力学の考えを放棄し、量子力学の考え方に以降しなければなりません。下記のように強調したいのです。「そこに気体分子があるからこのように考えられる」ではなく、「気体分子があると考えてこのように考える」というのが正しい表現であると言えるでしょう。「気体分子がそこにある」と考えて、圧力と温度が何かを理解していこうという姿勢に立って始めて納得いく議論へと展開できるのです。ミクロな視点では限界があるなぜ上記のようなごちゃごちゃ述べたかと言いますと、これが理論の全てではないからです。環境下が違えば、通用しない理論なんていくらでもあります。それが低温領域であり、ミクロな現象をミクロな視点で理解しようとしたときに古典論では限界があるということを言いたいのです。いって悲観的になる必要もないでしょう。物事の何を見ようとしているかによって、考え方を変えるだけで、ずいぶんと見通しの良い理論展開ができるのです。シンプルな仮定でも現象をできれば、それで十分なのです。何もてんこ盛りの理論でごっちゃがえすこともないでしょう。現象と違うのであれば、理論を少し追加すればいいだけなのです。※ただ、量子力学はそんな古典力学のちょっとした補正程度ではないですが(-_-)こちらの記事もどうぞ定積比熱と定圧比熱の温度依存性について!多項式で近似する。気体が膨張すると温度は下がるのか。それとも温度は上がるのか。-【CAE解析の目的とは?】CAEは実機の代替えなのか?精度を無視してパラメータスタディすることは意味がないのか?【Udemyがお勧め理由5つ】プログラミングを勉強するなら動画学習おすすめ【まとめ】比熱の温度依存性について。古典論と量子論から理解する。【理想気体の比熱の温度依存性(7)】シュレディンガー方程式の極座標表記(2原子分子)から動径方向だけを抽出しようコメントを残すコメントをキャンセルプロフィールカマキリ(^^)大学以上の物理を勉強してブログに書いていきます。大学の専攻は物性理論で、Fortranを使って数値計算をしていました。CAE解析(流体解析)・毎日note更新しています。内容は、「現在勉強中の内容」「日々思ったこと」を書いています。1DCAEの操作方法などを動画にしています。youtube・Qiitaではプログラミング言語の基本的な内容をまとめています。(10)ブログをメールで購読メールアドレスを記入して購読すれば、更新をメールで受信できます。