This preview shows page 1. Sign up to view the full content.
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
%%Display final theta value
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
theta=18.2*pi/180; %%Initial theta value
M=3.6; %%Intital Mach number
%%Initial calculations for theta = 18.2degre...
View Full Document
This homework help was uploaded on 03/31/2014 for the course AAE 334 taught by Professor Collicott during the Spring '09 term at Purdue University-West Lafayette.
- Spring '09