1 files changed, 62 insertions, 0 deletions
diff --git a/1040/CH3/EX3.2/Chapter3_Ex2.sce b/1040/CH3/EX3.2/Chapter3_Ex2.sce
new file mode 100644
index 000000000..4686ddaea
--- /dev/null
+++ b/1040/CH3/EX3.2/Chapter3_Ex2.sce
@@ -0,0 +1,62 @@
+//Harriot P.,2003,Chemical Reactor Design (I-Edition) Marcel Dekker,Inc.,USA,pp 436.
+//Chapter-3 Ex3.2 Pg No. 96
+//Title:Residence time and heat generation for four STR's in series
+//========================================================================================
+clear
+clc
+// COMMON INPUT
+X_A=0.95;//Given conversion
+t_batch=6;//Batch time to reach the desired conversion
+N=4//No.of reactors in series
+X_final=X_A;
+
+//CALCULATION (Ex3.2.a)
+k=log((1/(1-X_A)))/t_batch;//Refer equation 3.29 Pg No. 90
+t_1=((1/(1-X_A))^(1/N)-1)/k;//Refer equation 3.40 Pg No. 94
+t_Tot=N*t_1;
+
+//OUTPUT (Ex3.2.a)
+mprintf('\n OUTPUT Ex3.2.a');
+mprintf('\n==================================================================');
+mprintf('\nThe total residence time of the four reactors in series= %f hr',t_Tot);
+
+//=======================================================================================
+
+//Title:Heat generation in CSTR in Series
+//=============================================================================================================
+//CALCULATION (Ex3.2.b)
+t_1=((1/(1-X_final))^(1/N)-1)/k;//Refer equation 3.40 Pg No. 94
+for i=1:N
+    X(i)=1-(1/(1+k*t_1)^(i));
+end
+
+delQ_by_Q(1)=(X(1))/X_final; //  Ratio of heat generated in 1st reactor 
+for i=1:N-1
+    delQ_by_Q(i+1)=(X(i+1)-X(i))/X_final;  // Ratio of heat generated in 2nd, 3rd and 4th reactors 
+end
+
+//OUTPUT (Ex3.2.b)
+mprintf('\n========================================================================================\n')
+mprintf('\n OUTPUT Ex3.2.b');
+mprintf('\n==================================================================');
+mprintf('\nReactor vessel \t Conversion \t Fraction of total heat released \n')
+mprintf('\n==================================================================')
+for i=1:N
+    mprintf('\n %d \t \t %0.3f \t \t \t %0.3f \n',i,X(i),delQ_by_Q(i))
+end
+
+//FILE OUTPUT
+fid=mopen('.\Chapter3-Ex2-Output.txt','w');
+mfprintf(fid,'\n OUTPUT Ex3.2.a');
+mfprintf(fid,'\n==================================================================');
+mfprintf(fid,'\nThe total residence time of the four reactors in series= %f hr',t_Tot);
+    mfprintf(fid,'\n==================================================================')
+    mfprintf(fid,'\nReactor vessel \t Conversion \t Fraction of total heat released \n')
+    mfprintf(fid,'\n==================================================================')
+for i=1:N
+    mfprintf(fid,'\n %d \t \t %0.3f \t \t \t %0.3f \n',i,X(i),delQ_by_Q(i))
+end
+mclose(fid);
+
+
+//=============================================================END OF PROGRAM================================