徳城余市総合研究所トップページへ戻る
垂直衝撃波上流の全圧(Pt1,kPa)とマッハ数(M1)を与えて衝撃波通過後の全圧(Pt2,kPa)を計算
Octaveスクリプトファイル名：Pt1toPt2.m
-----------------------------------------------------------
# 垂直衝撃波上流の全圧(Pt1,kPa)とマッハ数(M1)を与えて衝撃波通過後の全圧(Pt2,kPa)を計算
# ーーー定数ーーー
k=1.4;
# ーーー入力ーーー
Pt1=input("Input:Pt1(kPa)=");
M1=input("Input:M1=");
# ーーー計算ーーー
Pt2=Pt1*((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("Pt2(kPa)=")
disp(Pt2)
-----------------------------------------------------------
(MM20120909)