Question

Code to Matlab

Make two plots on one page. mo ad tle veloclty as a function of time for 0sts4s. 33 The ideal gas equation states that PV RT 400 En van der Waals equation where P is the pressure, V is the volume, T is the temperature, R 0.08206 (L atm) (mol K) is the gas constant, and n is the number of moles. Real gases, especially at high pressures, deviate from this behavior. Their response can be modeled with the deal Gag oqualion 002 04 0608 Volume, (L)

Answer 1

Given data,

n=number of moles (Assume n=1)

P=Pressure (atm)

V=Volume (L) (0.065 <= V <=1)

T=Temperature (Assume T = 20C = 293K)

R=Gas constant = 0.08206 (L atm )/(mol K)

For CO₂ ,

a = L² atm / mol²

b = 0.04267 L/mol

Put all values in MATLAB command (define variables as given above)

>> n=1;R=0.08206;T=293;V=[0.065:0.005:1];

>> a=3.592;b=0.04267; %CO2 properties

>> P=(n*R*T)./V; % Ideal gas equation

>> Pw=((n*R*T)./(V-n*b))-((n^2*a)./(V.^2)); % Real gas equation

>> plot(V,P,'r',V,Pw,'k','LineWidth',2)
>> grid on
>> xlabel('Volume (L)');
>> ylabel('Pressure (atm)');
>> legend('Ideal Gas Eqaution','van der Waals Equation')

This will return a plot as shown in figure