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/CH14/EX14.10/Example_14_10.sce | 49 ++++++++++++++++++++++++++++++++++++++
 1 file changed, 49 insertions(+)
 create mode 100755 479/CH14/EX14.10/Example_14_10.sce

(limited to '479/CH14/EX14.10')

diff --git a/479/CH14/EX14.10/Example_14_10.sce b/479/CH14/EX14.10/Example_14_10.sce
new file mode 100755
index 000000000..f2d9716eb
--- /dev/null
+++ b/479/CH14/EX14.10/Example_14_10.sce
@@ -0,0 +1,49 @@
+//Chemical Engineering Thermodynamics
+//Chapter 14
+//Thermodynamics of Chemical Reactions
+
+//Example 14.10
+clear;
+clc;
+
+//Given
+//SO2(A) + (1/2)O2 (B) - SO3(C)
+//Basis: 1 Kgmole SO2
+n_A = 1;//Kgmole of SO2 fed
+n_B = n_A;//Kgmole of O2 fed
+T1 = 273+400;//Temperature in K at which reactants enter
+To = 298;//room temperature in K
+del_H = -23490;//Standar heat of reaction at 25 deg cel from example 14.6 in Kcal/Kgmole
+
+//At T1,
+C_A_T1 = 11.0; C_B_T1 = 7.4;//in Kcal/Kgmole
+//Assume the various temperature be
+T = [913 1073 1373 1573];//in K
+//Mean specific heats of the components A,B & C at various temperature are given below in Kcal/Kgmole K
+C_A = [11.6 11.8 12.3 12.5];
+C_B = [7.7 7.8 8.0 8.2];
+C_C = [  16.6 17.2 18.2 18.6];
+
+//To Calculate the final temperature for various conversions and the maximum conversion that can be attained in a single reactor operating adiabatically
+//In equation 14.18 (page no 307), H2-Hr = K  & (Hp-H3)= L(say)
+K = ((n_A*C_A_T1)+(n_B*C_B_T1))*(To-T1);//in Kcal/Kgmole
+mprintf('Adiabatic reaction temp in K             pecentage conversion of SO2');
+for i = 1:4
+    n_C(i) = (-K-(C_A(i)*(T(i)-To))-(C_B(i)*(T(i)-To)))/((T(i)-To)*(C_C(i)-C_A(i)-0.5*C_B(i))+del_H);
+mprintf('\n   %d',T(i));
+mprintf('                                              %f',n_C(i)*100);
+end
+
+clf;
+//plot(T,n_C*100);
+plot2d(T,n_C*100,style =2);
+//Now equilibrium conversion at various temperature taken from figure 14.7 (page no 325) are given as
+Ta = [850 900 1000 1100 1200 1300 1400];
+n_C1 = [93.5 88.2 69.0 49.0 37.0 21.5 6.25];
+plot2d(Ta,n_C1,style=3);
+xtitle("Temperature vs Percentage Conversion","Temperature in K","% Conversion");
+//From the graph,it can be seen that the curve cut each other approximately at the temp
+T1 = 1140;//in Kelvin
+C = interpln([Ta;n_C1],T1);
+mprintf('\n\n The maximum conversion that can be attained is %d percent',C);
+//end
\ No newline at end of file
-- 
cgit