Use mach number m2 as an input and this is accessed

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

View Full Document Right Arrow Icon
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

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.

Ask a homework question - tutors are online