ご返事遅くなり申し訳ございません。
ご教授頂きありがとうございました。

貼り付けたプログラムの事で質問です。 今長さ20ある中央でスロートになる管の中を空気が通る流れの数値解析を行っております。 フローチャートは 初期値→移流項→圧力項→境界条件→粘性項→境界条件→出力 という形です。 現状としては、移流項の空間微分の離散化を {U(I)+U(I+1)}　　　{U(I-1)+U(I)} ーーーーーー×F(I)-－－－－－－×F(I-1)　　　・・・（１） 　　　2　　　　　　　2 とするとスロート部でマッハ数１を超え、スロート部を右に少し行った所で衝撃波ができる結果となっていました。 そこで移流項の差分をもう少し精度をあげようとカワムラスキームやWENOスキームを...