Academic Integrity: tutoring, explanations, and feedback — we don’t complete graded work or submit on a student’s behalf.

Can you please provide answers for all questions(1-5), including the modified fa

ID: 3210109 • Letter: C

Question

Can you please provide answers for all questions(1-5), including the modified falsi code and main program.

Find the temperature T CC) of methane (CH:) at a pressure p 101.3 kPa and specific volume u 0.54997 m/kg using the Benedict-Webb-Rubin equation of state. Solve the equation using your Modified Regula Falsi program. Note that: 1) Absolute temperature must be used (K C273.15). p is in kPa and 0 is in 2) The universal gas constant R = 8.31434 kJ/kmol-"K and the molecular weight of 3) With the above units, the eight constants in the Benedict-Webb-Rubin equation of m/kmol. methane is M = 16.043 kg/kmol. state for methane are as follows: a = 5.00 A = 187.91 b = 0.003380 B = 0.04260 c = 2.578 x 105 C = 2.286 x 106 = 1.244 x 10" = 0.0060 lease submit the following: A plot which shows the location all positive real roots T(K) of the equation for the given values of p and v. 1) 2) A flow chart of the Modified Regula Falsi algorithm. 3) A Tisting of your Matlab main program and all function programs used 4) Printed output from your program showing the appropriate root T CK) computed to an accuracy of at least 6 significant figures 5) The calculated value of the temperature T CCJ using the Benedict-Webb-Rubin equation of state, the ideal gas equation of state, and the relative error of each based on the value of temperature reported in CH, tables (-161.45 ). Present a least 5 significant figures in your results.

Explanation / Answer

%variables
a = 5.00;
A = 187.91;
b = 0.003380;
B = 0.04260;
c = 2.57*10^5;
C = 2.268*10^6;
alpha = 1.244*10^(-4);
gama = 0.0060;
M = 16.043;
R = 8.31434;
v = 3;

temperature = input('Enter temperature in celsius: ');
%converting to kelvin
T = temperature + 273.15;

% using the modified equation
p = (R*T/v) + (B*R*T-A-(C/T*T))*(1/(v*v)) + (b*R*T-a)/(v*v*v) + (alpha*alpha)/(v^6) + (c/((v^3)*(T^3)))*(1+(gama/v^2))*exp(-1*gama/(v*v))

%{
output:

Enter temperature in celsius: 100
p = -2.5097e+05

%}