1 files changed, 38 insertions, 0 deletions
diff --git a/581/CH11/EX11.4/Example11_4.sce b/581/CH11/EX11.4/Example11_4.sce
new file mode 100755
index 000000000..0d7f30ad9
--- /dev/null
+++ b/581/CH11/EX11.4/Example11_4.sce
@@ -0,0 +1,38 @@
+
+clear;
+clc;
+
+printf("\t Example 11.4\n");
+T1=373.15;          //temp.of tea,K
+XN2=0.7808;    //mole fraction of nitrogen
+XO2=0.2095;    //mole fraction of oxygen
+Xar=0.0093;    //mole fraction of
+a=[3.798 3.467 3.542];    //collisin diameter,m
+b=[71.4 106.7 93.3];    //lennard jones constant,K
+M=[28.02 32 39.95];    //molar masses,kg/kmol
+c=[0.9599 1.057 1.021];    //potential well function
+d=[1.8*10^-5 2.059*10^-5 2.281*10^-5];    //calculated viscosity,kg/(m*s)
+e=[1.8*10^-5 2.07*10^-5 2.29*10^-5 ];    // theoritical viscosity,kg/(m*s)
+f=[0.0260 0.02615 0.01787];    //theoritical thermal conducitvity,W/(m*K)
+i=1;
+while(i<4)
+u(i)=2.6693*10^-6*(M(i)*T1)^0.5/((a(i)^2*c(i)));    //viscosity,kg/(m*s)
+k(i)=0.083228/((a(i))^2*c(i))*(T1/M(i))^0.5    //thermal conductivity,W/(m*s)
+
+i=i+1;
+end
+umc=XN2*u(1)/0.9978+XO2*u(2)/1.008+Xar*u(3)/0.9435 ;   //calculated mixture viscosity,kg/(m*s)
+umc1=1.857*10^-5;
+printf("\t theoritical mixture viscosity is : % -5e kg/(m*s)\n",umc1);
+umd=XN2*e(1)/0.9978+XO2*e(2)/1.008+e(3)*Xar/0.9435;    //theoritical mixture viscosity,kg/(m*s)
+printf("\t calculated mixture viscosity is : % -5e kg/(m*s)\n",umd);
+
+kmc=XN2*k(1)/0.9978+XO2*k(2)/1.008+Xar*k(3)/0.9435;    //calculated thermal conducitvity,W/(m*K)
+kmc1=0.02623;
+printf("\t theoritical thermal conducitvity is : %f W/(m*K)\n",kmc1);
+kmd=XN2*f(1)/0.9978+XO2*f(2)/1.008+Xar*f(3)/0.9435;    //theoritical thermal conductivity, W/(m*K)
+printf("\t calculated thermal conducitvity is : %.5f W/(m*K)\n",kmd);
+Cp=1006    //mixture diffusivity,j/(kg*K)
+pr=umd*Cp/kmd;    //prandtl no.
+printf("\t prandtl no. is : %.3f\n",pr);
+//end