Code Thermodynamics btech 3rd semester question to get results , plot graph if needed

A gas filled in a piston cylinder device, expands from 400 kPa and 0.05 m3

to a final volume of 0.2 m3

.

If the gas stisfies PVn=C, n and C are constants (take n=1.5), calculate the work done by the system.

Now, vary the value of n from 1.1 to 1.7, and plot the work versus polytropic exponent.

[Dyuman Joshi]

What does this question have to do with MATLAB? Specify it.

Answer 1

[VBBV]

p = 400*1e3;
v1 = 0.05;
v2 = 0.2;
n = 1.1:0.01:1.7;
W = (n./(n-1))*p*(v2-v1)

W = 1×61

1.0e+05 * 6.6000 6.0545 5.6000 5.2154 4.8857 4.6000 4.3500 4.1294 3.9333 3.7579 3.6000 3.4571 3.3273 3.2087 3.1000 3.0000 2.9077 2.8222 2.7429 2.6690 2.6000 2.5355 2.4750 2.4182 2.3647 2.3143 2.2667 2.2216 2.1789 2.1385

plot(n,W); xlabel('polytropic exponent [n]');ylabel('Work done')

W = @(n,p) (n./(n-1)).*p.*(v2-v1)

W = function_handle with value:
    @(n,p)(n./(n-1)).*p.*(v2-v1)

figure
fsurf(W,[1.1 1.7 100*1e3 400*1e3])
view(10,15)
axis square; zlabel('Work done')