真夏日

ようやく部屋の掃除・布団干しができた。
さわやかな素敵なお休み。
でも夕方からは大学へ。。。

昨日途中で断念していたポリトロープ大気の半解析解を求める計算を再開。
数値計算の方は長時間の計算で収束解が得られていたので、これと比較するための半解析解。
Gさんによると１週間ぐらいかかる、かなり面倒な計算とのこと。
でも来週はイベントがいろいろあってあまり研究できなさそうだから、今日中に終わらせておきたい。

ということで、ごちゃごちゃプログラムを書いて、さっそく計算させてみることに。

・・・数値計算の結果と一致！

意外と簡単だった。
まあGさんからまとめノートをもらっていたから、比較的サクッとやれたわけだが。

さて、次のステップはエネルギー方程式を入れての断熱の場合の計算。
この前一度やっていて、一応それっぽい収束解は得られたのだが、γの値をポリトロープ指数と同じにしてもポリトロープ大気の場合の解と一致しない。
これって一致しなきゃいけない、、、よね？

ちょっともう少し考えてみよう。
プログラムのミスかもしれないし。
明日は完全休養だから、月曜日に（月惑をサボって？）やり直してみよう。

早く一層加熱大気に移りたい。