第6-4回 MUSCL法と流束制限関数[python]

今回は,第6-2回と第6-3回でそれぞれ説明したMUSCL法と流速制限関数を組み合わせます。 まず,なぜ流速制限関数を用いなければならないのか、考えます。 第6-2回で算出したMUSCL法は次の形でした。 $$u^L_{j+1/2}=u_j+\frac{1}{4}(1-k)(u_j-u_{j-1})+\frac{1}{4}(1+k)(u_{j+1}-u_j)$$ 上記の式に対して、次のように変形を行います。 $$u^L_{j+1/2}-u_{j+1}=u_j-u_{j+1}+\frac{1}{4}(1-k)(u_j-u_{j-1})+\frac{1}{4}(1+k)(u_{j+1}-u_j)$$ $$⇔u^L_{j+1/2}-u_{j+1}=\frac{1} … 続きを読む 第6-4回 MUSCL法と流束制限関数[python]