2pi180 initial theta value m36 intital mach number

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: es outs1=obshk(1.4,M,theta,1); %%Store the physical quantities behind ...the first shock %%Compute physical properties behind the reflected oblique shock ...use Mach number M2 as an input and this is accessed as the 4th element in ...outs1 outs2=obshk(1.4,outs1(4),theta,1); while(~isnan(outs2(5))) %%Is the pressure ratio behind the reflected shock ... A real value? M=M-dM; %%Increment the Mach number %%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 value disp('Mach number') disp(M) Using this code we find to be PROBLEM 4 We can determine the Mach number by using both the Rayleigh Pitot formula and the oblique shock formulae. Guess the Mach Number M. Determine the free stream static pressure using the Rayleigh-Pitot formula and the measured pressure from the pitot probe. ( [ ) ( ) [ The pressure ratio across the oblique shock is then With the Mach number and half wedge angle as the input, we can determine the pressure ratio across the oblique shock using the obshk program. ( ) () Increment M till differs from by a specified tolerance. | | MATLAB CODE clc M=1.000; %%Input Mach number dM=0.01; %Mach Number...
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