Problem 5.31, page 220, from Felder and Rousseau Elementary Principles of Chemical Processes book
https://www.academia.edu/30223589/Felder_and_Rousseau_Elementary_Principles_of_Chemical_Processes_3.pdf I do not know how to upload the images to this page. The link above has the book.
The inputs are stated in the problem. You only need to have y = []; to pass.
I worked the problem out. It was fun and challenging for me. I wanted to share. Work it out any way you want and learn. My answers for Ra are [5.6267 7.4260 5.5061]. Yours answer might be different. They might be right or wrong.
Here my code.
% Read in the values of Rf, Rt, Rp, the fuel gas component mole fraction xa, xb, xc, % xd, and xe, and the percent excess air PX, and to find the required value of Ra. % Methane = a , ethane = b, propane = c, n-butane = d, isobutane = e
% The linear proportionalities between the in and out signals and the % corresponding process variables may be determined by:
% Fuel Temperature (K)
Tf1 = 25+273; Tf2 = 35+273; Rt1 = 14; Rt2 = 27;
% Fuel Pressure (atm) assume Patm = 1 atm
Pf1 = 1; Pf2 = 1 + 20000/101325; Rp1 = 0; Rp2 = 6;
% Fuel Flow rate (L/h)
Vf1 = 0; Vf2 = 2000*1000; Rf1 = 0; Rf2 = 10;
% Air Flow Rate (L(STP/h)
Va1 = 0; Va2 = 100000*1000; Ra1 = 0; Ra2 = 25;
% data readings
Rf = [7.25,5.8,2.45]; Rt = [23.1, 7.5 46.5]; Rp = [7.5, 19.3, 15.8];
% data mole%
xa = [.81,.58,0]; xb = [.08, .31, 0]; xc = [.05, .06, .65]; xd = [.04, .05, .25]; xe = [.02, 0, .10];
% data percent excess air
PX = [.15,.23,.33];
% Linear equations R => x T(K)
slopeT = 1/((Rt2 - Rt1)/(Tf2 - Tf1)); Tb = Tf1 - slopeT*Rt1; %Rt = linspace(Rt1,Rt2,11);
T = Rt*slopeT + Tb;
% Linear equations R => x Pa(atm)
slopeP = 1/((Rp2 - Rp1)/(Pf2 - Pf1)); Pb = Pf2 - slopeP*Rp2; %Rp = linspace(Rp1,Rp2,11);
P = Rp*slopeP + Pb;
% Linear equations R => x V(L/h)
slopeRf = 1/((Rf2 - Rf1)/(Vf2 - Vf1)); Vb = Vf1 - slopeRf*Rf1; %Rf = linspace(Rf1,Rf2,11);
Vf = Rf*slopeRf + Vb;
% assume ideal gas behavior PV = nRT, R = constant find n(Kmol/h)
R = 0.08206; %(atm *L)/(gmol*k) nf = ((P.*Vf)./(R*T))/1000; % Kmol/h
% species' mole in the fuel mixer Kmol/h
nfa = nf.*xa; nfb = nf.*xb; nfc = nf.*xc; nfd = nf.*xd; nfe = nf.*xe;
% Theoretical O2 % combustion chemical equations:
% 1. CH4 + 2O2 > CO2 + 2H2O methane
% 2. C2H6 + 3.5O2 > 2CO2 + 3H2O ethane
% 3. C3H8 + 5O2 > 3CO2 + 4H2O propane
% 4. C4H10 + 6.502 > 4CO2 + 5H2O BUTANE
% 5. C4H10 + 6.502 > 4CO2 + 5H2O ISOBUTANE
% no2s, s = stoichiometric when its assumes the species are completely consumed.
nO2s = nfa*2 + nfb*3.5 + nfc*5 + nfd*6.5 + nfe*6.5; %(kmol O2/h)
% excess of O2
nO2_feed = PX.*nO2s + nO2s; %kmol O2/h
% assume 21 mole%O2 and 79 mole%N2
nair_feed = nO2_feed/.21; %Kmol air/h
% assume air behavior is ideally PV = nRT find V() using T(stp) =273K,P(stp) = 1atm % 1000gmol = 1 kmol
Vair = (nair_feed.*273*.08206)*(1000); %(L/h)
% Linear equations R => x V(L(STP)/h)
slopeVa = (Ra2 - Ra1)/(Va2 - Va1); Vab = Ra1 - slopeVa*Va1; Ra = Vair.*slopeVa + Vab
Please reformulate the statement to a real Cody problem. This is not a blogging site.
14321 Solvers
Return a list sorted by number of occurrences
1225 Solvers
231 Solvers
135 Solvers
167 Solvers