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

---
 599/CH2/EX2.12/example2_12.sce | 31 +++++++++++++++++++++++++++++++
 1 file changed, 31 insertions(+)
 create mode 100755 599/CH2/EX2.12/example2_12.sce

(limited to '599/CH2/EX2.12/example2_12.sce')

diff --git a/599/CH2/EX2.12/example2_12.sce b/599/CH2/EX2.12/example2_12.sce
new file mode 100755
index 000000000..7b641863c
--- /dev/null
+++ b/599/CH2/EX2.12/example2_12.sce
@@ -0,0 +1,31 @@
+
+clear;
+clc;
+printf("\t Example 2.12\n");
+// ammonia diffusing through inert air and air is non-diffusing
+ya1=0.1;
+ya2=0;
+T=(293);                             //temperature  in kelvin
+pt=1*1.013*10^5;                    //total pressure in pascal
+z=0.2*10^-2;                       //gas film thickness in m
+Dab=.185*10^-4;                   //diffusion coefficient in m^2/s
+R=8314;                          //universal gas constant
+    //part (i)when air is assumed to be stagnant and non-diffusing 
+
+Na=Dab*pt*log((1-ya2)/(1-ya1))/(z*R*T);    //diffusion flux in kmol/m^2*s
+mw=17;                                    //molecular weight of ammonia
+massflux=Na*mw;                          //mass flux of given NH3
+printf("\n  diffusion flux when total presssure is 1atm and air is non-diffusing :%f *10^-4 kg/m^2*s ",massflux/10^-4);
+    //part (ii) when pressure is increased to 10atm
+
+//Dab_1/Dab_2=pt_2/pt_1
+pt_2=10;                                         //final pressure  in atm
+pt_1=1;                                         //initially pressure was 1atm
+Dab_1=.185;                                    //initially diffusion coefficient was.185
+Dab_2=Dab_1*pt_1/pt_2;                        //for gases Dab is proportional to 1/pt
+Dab=Dab_2*10^-4;                             //new diffusion coefficient 
+pt=pt_2*1.013*10^5;                         //new total pressure
+Na=Dab*pt*log((1-ya2)/(1-ya1))/(z*R*T);    //diffusion flux in kmol/m^2*s
+printf("\n  diffusion  flux when pressure is increased to 10atm :%f *10^-5 kmol/m^2*s ",Na/10^-5);
+printf("\n \n so the rate of diffusion remains same on increasing the pressure");
+//end
\ No newline at end of file
-- 
cgit