RMHD

moppy 氏と YSN 氏から計算技術と物理との兼ね合い，意味合いの指導を受ける．なるどほなるほどと頷くこと頻り．眼を凝らせば物理が見えてくる．Cartoon 法ってそういうことか．
早速，円筒計算にとりかかる．こういう座標変換には測地線方程式を使って，，と思いつつ，このくらい別に相対論を持ち出すまでもないんだ．ヤコビアンが r だから，あとは div を間違えないように計算すればいいだけ．おお，上手いこと保存形式に落ち着く．しかし予想通り，余計な項も出てきた．これは右辺に置いておいて，スペクトル分解して，ヤコビアンを計算すればいいわけだ．今日は質量保存をまとめるところまで．
opening angle を動かしたいときは，球座標がベストだと思われるが，計算大変そう．今の系でどうやってエネルギーを deposit させるかが知恵の使いどころ．