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

1 view (last 30 days)
vikram kumar on 2 Dec 2022
Edited: VBBV on 2 Jan 2023
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 on 2 Dec 2022
What does this question have to do with MATLAB? Specify it.

VBBV on 2 Jan 2023
Edited: VBBV on 2 Jan 2023
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') 