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 --- 914/CH2/EX2.6/ex2_6.sce | 20 ++++++++++++++++++++ 1 file changed, 20 insertions(+) create mode 100755 914/CH2/EX2.6/ex2_6.sce (limited to '914/CH2/EX2.6') diff --git a/914/CH2/EX2.6/ex2_6.sce b/914/CH2/EX2.6/ex2_6.sce new file mode 100755 index 000000000..3b859df23 --- /dev/null +++ b/914/CH2/EX2.6/ex2_6.sce @@ -0,0 +1,20 @@ +clc; +warning('off'); +printf("\n\n example2.6 - pg36"); +// given +d=1.2*10^-2; //[m] - diameter of the hole +Ca1=0.083; //[kmol/m^3] +Ca2=0; //[kmol/m^3] +L=0.04; //[m] - thickness of the iron piece +Dab=1.56*10^-3; //[m^2/sec] - diffusion coefficient of CO2 +A=(%pi*d^2)/4; //area +// (a)using the formula (Na/)A=(Ja/A)=-Dab(delCa/delx) +intdCa=integrate('1','Ca',Ca2,Ca1); +intdx=integrate('1','x',0,0.04); +g=(intdCa/intdx)*Dab; +printf("\n\n (a) The molar flux with respect to stationary coordinates is\n (Na/A)=%fkmol/m^2*sec",g); +// using the formula na/A=(Na/A)*Ma +Ma=44.01; //[kg/mol] - molcular weight of co2 +na=(intdCa/intdx)*Dab*Ma*A*(3600/0.4539); +printf("\n\n The mass flow rate is %flb/hr",na); + -- cgit