Question

NOTE: IN THIS MATLAB PROGRAMMING EXERCISE, WRITE THE PROGRAM FIRST BY HAND, ON PAPER. THEN KEY YOUR PROGRAM INTO MATLAB, AND RUN IT TO PRODUCE OUTPUT. "WORK" FOR THIS EXERCISE WILL CONSIST OF YOUR HANDWRITTEN PROGRAM AND A SCREENSHOT OF YOUR PROGRAM IN MATLAB AS WELL AS THE OUTPUT YOUR PROGRAM PRODUCES MATLAB PROGRAMMING EXERCISE 8 (From physical chemistry; 4 pts.) Using the Redlich-Kwong equation of state for a real gas, we can calculate the work produced by that gas when the gas undergoes a reversible, isothermal expansion at temperature T, from initial molar volume Vi to final molar volume V2. That expression for reversible, isothermal work is: work-_RTIn where: R 0.082057; V1 = 0.5; V2 = 1.0; Calculate the work produced by each of the following gases, at T 693: GAS WORK PRODUCED Neon 1.4439 0.012049 Methane 31.784 0.029850 Carbon Dioxide 63.752 0.029677 Ammonia 86.660 0.026232

Answer 1

If you like the code, give it a thumbs up. Incase you need any clarifications let me know in comments.

Code starts below this line

%Initialize constants
R=0.082057;
V1=0.5;
V2=1.0;
T=693;

%Create vector for Gas, A and B values
Gas={'Neon','Methane','Carbon dioxide','Ammonia'};
A=[1.4439, 31.784, 63.752, 86.660];
B=[0.012049, 0.029850, 0.029677, 0.026232];

%Caculate work for each of the possible combinations
for i=1:length(A)
work(i)=-R*T*log((V2-B(i))/(V1-B(i)))-((A(i)/(B(i)*sqrt(T)))*log(((V2+B(i))*V1)/((V1+B(i))*V2)));
end

%Display the output where last bracket contains the work produced values
fprintf('Gas \t\t\t A \t\t B \t\t Work Produced\n');
for i=1:length(A)
fprintf('%s \t %f \t %f \t %f\n',cell2mat(Gas(i)),A(i),B(i),work(i));
end