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 --- 2126/CH3/EX3.13/13.sce | 30 ++++++++++++++++++++++++++++++ 1 file changed, 30 insertions(+) create mode 100755 2126/CH3/EX3.13/13.sce (limited to '2126/CH3/EX3.13') diff --git a/2126/CH3/EX3.13/13.sce b/2126/CH3/EX3.13/13.sce new file mode 100755 index 000000000..623846e6e --- /dev/null +++ b/2126/CH3/EX3.13/13.sce @@ -0,0 +1,30 @@ +clc +clear + +//input data +Po1=200 //Stagantion pressure at inlet in kPa +To1=303 //Stagnation temperature at inlet in K +M1=0.2 //Inlet Mach number from diagram +D=0.025 //inner tude diameter in m(missing data) +M2=0.8 //Outlet Mach number +f=0.005/4 //frictional factor +R=287 //Gas constant in J/kg-K + +//calculation +t1=0.992 //Static to Stagnation temperature ratio at entry from gas tables (M1,k=1.4,isentropic) +T1=To1*t1 //Static temperature in K +p1=0.973 //Static to Stagnation pressure ratio at entry from gas tables (M1,k=1.4,isentropic) +P1=Po1*p1 //Static pressure in kPa +p2=2.964 //Stagnation pressure ratio at inlet to critical state from gas tables (M1,k=1.4,fanno flow) +Pot=Po1/p2 //Stagnation pressure at critical state in kPa +X1=14.533 //frictional constant fanno parameter from gas tables,fanno flow tables @M1,k=1.4 +p3=1.038 //Stagnation pressure ratio at outlet to critical state from gas tables (M1,k=1.4,fanno flow) +Po2=Pot*p3 //Stagnation pressure at exit in kPa +X2=0.073 //frictional constant fanno parameter from gas tables,fanno flow tables @M2,k=1.4 +X3=X1-X2 //overall frictional constant fanno parameter +L1=(X3*D)/(4*f) //Length of the pipe in m +SPL=(1-(p3/p2))*100 //Percentage decrease in stagnation pressure in percent +ds=R*log(Po1/Po2) //Change of entropy in kJ/kg-K + +//output +printf('(A)Length of the pipe is %3.1f m\n (B)Percentage decrease in stagnation pressure is %3.2f percent\n (C)Change of entropy is %3.3f kJ/kg-K',L1,SPL,ds) -- cgit