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

---
 536/CH10/EX10.2/Example_10_2.sce | 57 ++++++++++++++++++++++++++++++++++++++++
 1 file changed, 57 insertions(+)
 create mode 100755 536/CH10/EX10.2/Example_10_2.sce

(limited to '536/CH10/EX10.2')

diff --git a/536/CH10/EX10.2/Example_10_2.sce b/536/CH10/EX10.2/Example_10_2.sce
new file mode 100755
index 000000000..5f223922a
--- /dev/null
+++ b/536/CH10/EX10.2/Example_10_2.sce
@@ -0,0 +1,57 @@
+clear;
+clc;
+
+printf("Example 10.2\n");
+
+th=[0 0 3 7 22 32 46 55 80 106];      //Time in hours
+tm=[0 26 5 36 16 38 50 25 22 25];     //Time in min
+
+//Conversion to kilo seconds
+for i=1:10
+    tm(i)=tm(i)*60;
+    th(i)=th(i)*3600;
+    tim(i)=(tm(i)+th(i))/1000;
+end
+
+L=[0 2.5 12.9 23.2 43.9 54.7 67.0 73.8 90.3 104.8];   //in mm
+
+Lo=L(1);
+
+//Calculations are done as indicated in the procedure
+//To obtain the values of x and y as below
+//For plotting x and t axis of graph
+x=L-Lo;
+
+y(1)=0;
+for j=2:10
+    y(j)=tim(j)/(L(j)-Lo);
+end
+
+plot2d(x,y);
+plot(x,y,'+');
+xtitle('t/(L-L0) vs (L-L0)','(L-L0) in mm','t/(L-L0) in ks/mm^2');
+
+//Calculation of slope
+s=(y(4)-y(3))/(x(4)-x(3))*10^3*10^6;
+printf("\nSlope is %.2e sec/m^2\n",s);
+
+Vl=22.4;                 //Molar volume in litres
+den=1540;                //Density in kg/m^3
+T0=273;
+T=321;
+vp=37.6;                 //vapour pressure in kPa
+P0=101.3;                //PRessue in kPa
+M=154;
+
+Ct=T0/(Vl*T);
+Ca=(vp*Ct)/P0;
+
+Cb1=Ct;
+Cb2=(P0-vp)*Ct/P0;
+Cbm=(Cb1-Cb2)/log(Cb1/Cb2);
+
+//Diffusivity calculation
+D=den*Cbm/(2*M*Ca*Ct*s);
+printf("\nDiffusivity is %.2e m^2/s\n",D);
+
+//End
-- 
cgit