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 --- 479/CH7/EX7.3/Example_7_3.sce | 47 +++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 47 insertions(+) create mode 100755 479/CH7/EX7.3/Example_7_3.sce (limited to '479/CH7/EX7.3') diff --git a/479/CH7/EX7.3/Example_7_3.sce b/479/CH7/EX7.3/Example_7_3.sce new file mode 100755 index 000000000..c40244580 --- /dev/null +++ b/479/CH7/EX7.3/Example_7_3.sce @@ -0,0 +1,47 @@ +//Chemical Engineering Thermodynamics +//Chapter 7 +//Ideal Gases + +//Example 7.3 +clear; +clc; + +//Given +P1 = 15;//initial pressure in Kgf/cm^2 +P2 = 1;//final pressure in Kgf/cm^2 +V1 = 0.012;//initial volume in m^3 +V2 = 0.06;//final volume in m^3 +T1 = 420;//initial temperature in K +M = 28;//molecular weight of the gas +Cp = 0.25;//specific heat at constant pressure in Kcal/Kg K +R = 1.98;//gas constant in Kcal/Kg mole K +R2 = 848;//gas constant in mKgf/Kgmole K +//Cv = a+0.0005*T1; Specific heat at constant volume + +//To Calculate the final temperature of the ideal gas, work done in an open and closed system,internal energy change for the process +//(a)Calculation of final temperature +//Using ideal gas law:(P*V)/(R*T) +T2 = (P2*V2*T1)/(P1*V1); +mprintf('(a)The final temperature is %d K',T2); + +//(b)Calculation of work in an open and closed system +//From equation 7.22(page no 147): P1*(V1^n)=P2*(V2^n) +n = (log(P2/P1))/(log(V1/V2)); +//From equation 7.25(page no 149) +W = ((P1*V1)-(P2*V2))/(n-1)*10^4;//work in mKgf +W1 = W/427;//Work in Kcal +mprintf('\n (b)The work in a closed system is %f Kcal',W1); +Ws = n*W1;//from equation 7.28(page no 149) +mprintf('\n The work in an open system is %f Kcal',Ws); + +//(c)Calculation of internal energy change +R1 = R/M;//gas constant in Kcal/Kg +Cv = Cp-R1;//specific heat at constant volume in Kcal/Kg K +a = Cv-(0.0005*T1); +m = (P1*10^4*V1*M)/(R2*T1);//mass of gas in Kg +function y = f(T) + y = m*(a+(0.0005*T)); +endfunction +del_E = intg(T1,T2,f);//internal energy change in Kcal/Kg +del_E1 = M*del_E;//internal energy change in Kcal/Kgmole +mprintf('\n (c)The internal energy change for the process is %f Kcal/Kgmole',del_E1); \ No newline at end of file -- cgit