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