計算スキーム探り

ここ最近熱中して安定で精度がよくて計算コストが少ない計算スキームはないかなあといろいろと探っている。

最近流行のHLLDを実装したいのだけれど、その準備として
Lax-Friedrich+MUSCL+TVDというスキームと
GMUSTE-1というスキームを実装してみた。

衝撃波管の問題(SOD)の問題がきれいに解けるととても気持ちいい

来週は地球惑星連合大会だがしばらく熱中し続けよう