「対流項を含む偏微分方程式の数値解法」で使用したスプレッドシートを公開するにゃ。


CTFSのエクセル版は


https://docs.google.com/spreadsheets/d/e/2PACX-1vQZgPy2aBZNO5s-q5Vs-3sL0qwV_uF8Oqsj9ZIGkYVZZl24ZBvRxqm4t5SuudMwWWzN41uAqga9dobb/pub?output=xlsx


odsファイルは


https://docs.google.com/spreadsheets/d/e/2PACX-1vQZgPy2aBZNO5s-q5Vs-3sL0qwV_uF8Oqsj9ZIGkYVZZl24ZBvRxqm4t5SuudMwWWzN41uAqga9dobb/pub?output=ods


だにゃ。


 


計算結果だけを知りたいヒトは、


https://docs.google.com/spreadsheets/d/e/2PACX-1vQZgPy2aBZNO5s-q5Vs-3sL0qwV_uF8Oqsj9ZIGkYVZZl24ZBvRxqm4t5SuudMwWWzN41uAqga9dobb/pubhtml


 


CTFSを用いた計算では、計算領域の右端x=10（L列のセル）の値を0として計算してある。


このことの是非については、１５時公開予定の「ネムネコの脳内仮想会議」という記事に書くにゃ。


計算結果を見るとわかると思うけれど、計算法上、このセルに０以外の値が入るのは、１０ステップ計算をしたあと――上流から１ステップ計算すると、その値は隣接する右のセルに伝わるため――。しかもその値は、少なくとも２０ステップくらいまでは無視できるくらいの小さな値。だから、少なくとも２０ステップくらいまでは、この計算は全体的に信用できると思うにゃ。記事中に示してあるのは、x=1〜4に相当するC〜F列で、オーダー的に考えて、５０ステップくらいまでは、まぁ信用できると思う。


 


風上（上流）差分を用いたものはコチラ↓


エクセル版


https://docs.google.com/spreadsheets/d/e/2PACX-1vTxVAt5H4Yfkm8Zh2DlDedU6bkYh8GZ9yODCdKwgOt52FfxSqBzyvnEKwta_sh2tLPJEpWsl6HlmU8x/pub?output=xlsx


OpenDocument版


https://docs.google.com/spreadsheets/d/e/2PACX-1vTxVAt5H4Yfkm8Zh2DlDedU6bkYh8GZ9yODCdKwgOt52FfxSqBzyvnEKwta_sh2tLPJEpWsl6HlmU8x/pub?output=ods


 


計算結果だけを知りたいヒトは、ここを見るといいと思うにゃ↓


https://docs.google.com/spreadsheets/d/e/2PACX-1vTxVAt5H4Yfkm8Zh2DlDedU6bkYh8GZ9yODCdKwgOt52FfxSqBzyvnEKwta_sh2tLPJEpWsl6HlmU8x/pubhtml


 


 


「おい、⑨ネコ。お前のことだから、このスプレッドシート以外にFortranプログラムも作って計算したはずだ。でなければ、あの記事の最後に出ている２つのグラフは書けなかったに違いないのだから。それも公開しろ！！」


「さて、何のことやら(^^ゞ」