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 --- 172/CH14/EX14.6/ex6.sce | 29 +++++++++++++++++++++++++++++ 1 file changed, 29 insertions(+) create mode 100755 172/CH14/EX14.6/ex6.sce (limited to '172/CH14/EX14.6/ex6.sce') diff --git a/172/CH14/EX14.6/ex6.sce b/172/CH14/EX14.6/ex6.sce new file mode 100755 index 000000000..ec041b50a --- /dev/null +++ b/172/CH14/EX14.6/ex6.sce @@ -0,0 +1,29 @@ +//ques6 +//isothermal steady state processes +clear +clc +//from table A.2 +P1=8;//pressure at state 1 in MPa +P2=0.5;//pressure at state 2 in MPa +T1=150;//Temperature at state 1 in K +Pr1=P1/3.39;//Reduced pressure at state 1 +Pr2=P2/3.39;//Reduced pressure at state 2 +Tr1=T1/126.2;//Reduced temperature +T2=125;//temperature at state 2 +//from compressibility chart h1*-h1=2.1*R*Tc +//from zero pressure specific heat data h1*-h2*=Cp*(T1a-T2a) +//h2*-h2=0.5*R*Tc +//this gives dh=h1-h2=-2.1*R*Tc+Cp*(T1a-T2a)+0.15*R*Tc +R=0.2968;//gas constant for given substance +Tc=126.2;//K, Constant temperature +Cp=1.0416;//heat enthalpy at constant pressure in kJ/kg +dh=(2.35)*R*Tc+Cp*(T2-T1);// +printf('Enthalpy change = %.2f kJ/kg \n',dh); +//change in entropy +//ds= -(s2*-s2)+(s2*-s1*)+(s1*-s1) +//s1*-s1=1.6*R +//s2*-s2=0.1*R +//s2*-s1*=Cp*log(T2/T1)-R*log(P2/P1) +//so +ds=1.6*R-0.1*R+Cp*log(T2/T1)-R*log(P2/P1); +printf(' Entropy Change = %.4f kJ/kg.K ',ds); \ No newline at end of file -- cgit