HW4Problem11 - phi(i) = f(i)/fstoich; end figure hold on...

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

View Full Document Right Arrow Icon
AOE 4234 November 13, 2009 HW4 Sossi, Craig 904512344 1200 1300 1400 1500 1600 1700 1800 0 0.0025 0.005 0.0075 0.01 0.0125 0.015 0.0175 0.02 0.0225 0.025 0.0275 0.03 0.0325 0.035 Temperature Fuel-Air Ratio -0.8 -0.6 -0.4 -0.2 0 Equivalence Ratio % Craig Sossi % AOE 4234 - HW 4 - Problem 6.11 clear all close all clc CpO2 = 33.8; CpH2 = 30.2; CpN2 = 32.0; CpH2O = 40.0; MH2 = 2.02; MAir = 137.318; QR = 121000; fstoich = 2*MH2/MAir; T = linspace(1200,1800,100); f = zeros(length(T)); phi = zeros(length(T)); for i = 1:length(T) f(i) = (-fstoich/2)*((CpO2+3.76*CpN2)*(T(i)-600))/(QR+CpH2O*(T(i)-298)-CpO2/2*(T(i)-298)+CpH2*(221));
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: phi(i) = f(i)/fstoich; end figure hold on [AX,H1,H2] = plotyy(T,f,T,phi, 'plot' ); xlabel( 'Temperature' , 'fontsize' ,16) axis(AX(1)) ylim([0 .035]) set(gca, 'ytick' ,0:0.0025:0.035) set(get(AX(1), 'Ylabel' ), 'String' , 'Fuel-Air Ratio' , 'fontsize' ,16) set(get(AX(2), 'Ylabel' ), 'String' , 'Equivalence Ratio' , 'fontsize' ,16)...
View Full Document

This note was uploaded on 01/23/2012 for the course AOE 4234 taught by Professor Shinpaugh during the Fall '08 term at Virginia Tech.

Ask a homework question - tutors are online