From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 779/CH18/EX18.6/18_6.sce | 20 ++++++++++++++++++++ 1 file changed, 20 insertions(+) create mode 100755 779/CH18/EX18.6/18_6.sce (limited to '779/CH18/EX18.6') diff --git a/779/CH18/EX18.6/18_6.sce b/779/CH18/EX18.6/18_6.sce new file mode 100755 index 000000000..36a4b9ac1 --- /dev/null +++ b/779/CH18/EX18.6/18_6.sce @@ -0,0 +1,20 @@ +P1 = 1; P41 = 15; // P41 = P4/P1 +P21 = (P41)^(1/3); +P2 = P21*P1; n = 1.3; R = 0.287; +P3 = P21*P2; +P11 = P2; P12 = P1; +b = 0.45; s = 0.3; // Bore and stroke of cylinder +Vs = (%pi/4)*b^2*s; // Swept volume of the cylinder +V11 = 0.05*Vs; // Clearance volume +V1 = V11+Vs; +V12 = V11*(P11/P12)^(1/n); +disp("m3",V1-V12,"Effective swept volume of the LP cylinder is") +T1 = 291; T3 = T1; T5 = T1; +P43 = P21; // P4/P3 +T6 = T5*(P43)^((n-1)/n); +disp("K",T6,"Delivery temperature is") +P4 = 15; // Delivery pressure +V6_7 = (P1/P4)*(T6/T1)*(V1-V12); // V6-V7 +disp("m3",V6_7,"Volume of the air delivered") +W = ((3*n*R*T1)/(n-1))*((P21)^((n-1)/n)-1); +disp("kJ",W,"Work done per kg of the air is") -- cgit