HW4Problem11

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

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

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));
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