{[ promptMessage ]}

Bookmark it

{[ promptMessage ]}

PS5_solution

Use mach number m2 as an input and this is accessed

This preview shows page 1. Sign up to view the full content.

This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: theta,1); while(~isnan(outs2(5))) %%Is the pressure ratio behind the reflected shock ... A real value? theta=theta+dtheta; %%Increment the angle of the compression corner %%Perform the oblique shock calculations once more outs1=obshk(1.4,M,theta,1); outs2=obshk(1.4,outs1(4),theta,1); end %%Display final theta value disp('theta') disp(theta*180/pi) Using this code we find to be PROBLEM 3c The MATLAB program decreases in steps of 0.001 till the obshk program puts out a value of NaN for the reflected oblique shock Pressure ratio. The final value for is smaller than the Mach number of 3.6 in part (a). This trend is expected since decreasing the Mach number decreases the value of M2, the Mach number behind the first shock. In other words, the upstream Mach number for the reflected shock is reduced but the deflection angle is kept the same - this scenario favors Mach reflection. MATLAB CODE clear clc theta=18.2*pi/180; %%Initial theta value M=3.6; %%Intital Mach number dM=0.001; %%Initial calculations for theta = 18.2degre...
View Full Document

{[ snackBarMessage ]}

Ask a homework question - tutors are online