徳城余市総合研究所トップページへ戻る
斜め衝撃波上流の全圧(Pt1,kPa)とマッハ数(M1)と偏角（Theta,deg）を与えて衝撃波通過後の全圧(Pt2,kPa)を計算
Octaveスクリプトファイル名：OPt1toPt2.m
-----------------------------------------------------------
# 斜め衝撃波上流の全圧(Pt1,kPa)とマッハ数(M1)と偏角（Theta,deg）を与えて衝撃波通過後の全圧(Pt2,kPa)を計算
# 適用範囲：M=2〜10程度
# ーーー定数ーーー
k=1.4;
# ーーー入力ーーー
Pt1=input("Input:Pt1(kPa)=");
M1=input("Input:M1=");
Theta=input("Input:Theta(deg)=");
# ーーー計算1ーーー
Ti=Theta*pi/180;
j=1.05/M1+0.5*Ti;
Theta=atan((2*cot(j)*((M1*sin(j))^2-1))/(M1^2*(k+cos(2*j))+2));
while abs(Ti-Theta)>0.000001
j=j+0.5*(Ti-Theta);
Theta=atan((2*cot(j)*((M1*sin(j))^2-1))/(M1^2*(k+cos(2*j))+2));
endwhile
Beta=j;
# ーーー計算2ーーー
Pt2=Pt1*((1+2*k/(k+1)*((M1*sin(Beta))^2-1))^(-1/(k-1)))*(((k+1)/2*(M1*sin(Beta))^2)/(1+(k-1)/2*(M1*sin(Beta))^2))^(k/(k-1));
# ーーー出力ーーー
disp("Pt2(kPa)=")
disp(Pt2)
-----------------------------------------------------------
(MM20121001)