This preview shows page 1. Sign up to view the full content.
Unformatted text preview: rogram increments in steps of 0.001 degrees till the obshk program puts out a value of
NaN for the reflected oblique shock pressure ratio. The final value for is greater than the compression
corner angle in part (a). This trend is expected since increasing the compression corner angle decreases
the value of M2, the Mach number behind the first shock while simultaneously increasing the value of
the deflection angle at the upper wall. In other words, the upstream Mach number for the reflected
shock is lower and the deflection angle is higher  this scenario favors Mach reflection. MATLAB CODE
clear
clc
theta=0*pi/180; %%Initial theta value
dtheta=0.001*pi/180; %%Incremental angle
M=3.6; %%Intital Mach number
%%Initial calculations for theta = 0 degrees
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),...
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 UniversityWest Lafayette.
 Spring '09
 COLLICOTT

Click to edit the document details