徳城余市総合研究所トップページへ戻る
垂直衝撃波上流の静温(Ts1,K)とマッハ数(M1)を与えて衝撃波通過後の静温(Ts2,K)を計算
Octaveスクリプトファイル名：Ts1toTs2.m
-----------------------------------------------------------
# 垂直衝撃波上流の静温(Ts1,K)とマッハ数(M1)を与えて衝撃波通過後の静温(Ts2,K)を計算
# ーーー定数ーーー
k=1.4;
# ーーー入力ーーー
Ts1=input("Input:Ts1(K)=");
M1=input("Input:M1=");
# ーーー計算ーーー
Ts2=Ts1*(1+2*(k-1)/(k+1)^2*(k*M1^2+1)/M1^2*(M1^2-1));
# ーーー出力ーーー
disp("Ts2(K)=")
disp(Ts2)
-----------------------------------------------------------
(MM20120830)