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 --- 2063/CH1/EX1.1/1_1.sce | 26 ++++++++++++++++++++++++++ 1 file changed, 26 insertions(+) create mode 100755 2063/CH1/EX1.1/1_1.sce (limited to '2063/CH1/EX1.1/1_1.sce') diff --git a/2063/CH1/EX1.1/1_1.sce b/2063/CH1/EX1.1/1_1.sce new file mode 100755 index 000000000..19c744842 --- /dev/null +++ b/2063/CH1/EX1.1/1_1.sce @@ -0,0 +1,26 @@ +clc +clear +//Input data +V1=0.5;//Initial Volume before the commencement of compression in m^3 +P1=1;//Initial pressure before the commencement of compression in bar +T1=300;//Initial temperature in K +P2=12;//Final pressure at the end of compression stroke in bar +Q=220;//Heat added during the constant volume process in kJ +r=1.4;//Isentropic constant for air +R=0.287;//Characteristic Gas constant in kJ/kg K +Cv=0.718;//Specific heat of mixture in kJ/kg K + +//Calculations +r1=(P2/P1)^(1/r);//Compression ratio +T2=T1*(r1)^(r-1);//Final temperature after the end of compression stroke in K +V2=(P1*T2*V1)/(P2*T1);//Final volume after the end of compression stroke in m^3 +m=(P1*10^5*V1)/(R*T1*1000);//Mass of air flowing in kg +T3=(Q/(m*Cv))+T2;//Temperature after constant volume heat addition in K +P3=(P2*T3)/T2;//Pressure after constant volume heat addition in K +V3=V2;//Volume at 3 +P4=P3*(1/r1)^(r);//Pressure after isentropic expansion in bar +V4=V1;//Volume after isentropic expansion in m^3 +T4=T3*(1/r1)^(r-1);//Temperature at the end of isentropic expansion in K + +//Output +printf('(a)The pressures at 1 is %3.0fbar\n (b)Pressure at 2 is %3.0fbar\n (c)Pressure at 3 is %3.2fbar\n (d)Pressure at 4 is %3.2fbar\n (e)Temperature at 1 is %3.1fK\n (f)Temperature at 2 is %3.1fK\n (g)Temperature at 3 is %3.0fK\n (h)Temperature at 4 is %3.0fK\n (i)Volume at 1 is %3.0fm^3\n (j)Volume at 2 is %3.5fm^3\n (k)Volume at 3 is %3.5fm^3\n (l)Volume at 4 is %3.0fm^3',P1,P2,P3,P4,T1,T2,T3,T4,V1,V2,V3,V4) -- cgit