1 files changed, 32 insertions, 0 deletions
diff --git a/914/CH14/EX14.5/ex14_5.sce b/914/CH14/EX14.5/ex14_5.sce
new file mode 100755
index 000000000..b0f45eb05
--- /dev/null
+++ b/914/CH14/EX14.5/ex14_5.sce
@@ -0,0 +1,32 @@
+clc;
+warning("off");
+printf("\n\n example14.5 - pg730");
+// given
+T=323;  //[K] - temperature
+P=1;  //[atm] - pressure
+Dab_experimental=7.7*10^-6;  //[m^2/sec]
+DPM_A=1.9;  // dipole moment of methyl chloride
+DPM_B=1.6;  // dipole moment of sulphur dioxide
+Vb_A=5.06*10^-2;  // liquid molar volume of methyl chloride
+Vb_B=4.38*10^-2
+Tb_A=249;  // normal boiling point of methyl chloride
+Tb_B=263;  // normal boiling point of sulphur dioxide
+del_A=((1.94)*(DPM_A)^2)/(Vb_A*Tb_A);
+del_B=((1.94)*(DPM_B)^2)/(Vb_B*Tb_B);
+del_AB=(del_A*del_B)^(1/2);
+sigma_A=(1.166*10^-9)*(((Vb_A)/(1+1.3*(del_A)^2))^(1/3));
+sigma_B=(1.166*10^-9)*(((Vb_B)/(1+1.3*(del_B)^2))^(1/3));
+etaabykb=(1.18)*(1+1.3*(del_A^2))*(Tb_A);
+etabbykb=(1.18)*(1+1.3*(del_B^2))*(Tb_B);
+sigma_AB=(1/2)*(sigma_A+sigma_B);
+etaabbykb=(etaabykb*etabbykb)^(1/2);
+Tstar=T/(etaabbykb);
+sigmaDnonpolar=1.602;
+sigmaDpolar=sigmaDnonpolar+(0.19*(del_AB^2))/Tstar;
+patm=1;
+Ma=50.488;  //[kg/mole] - molecular weight of methyl chloride
+Mb=64.063;  //[kg/mole] - molecular weight of sulphur dioxide 
+D_AB=(1.8583*(10^-(27))*(((T^3)*((1/Ma)+(1/Mb)))^(1/2)))/(patm*(sigma_AB^(2))*sigmaDpolar);
+printf("\n\n Dab=%em^2/sec",D_AB);
+printf("\n\n The Chapman Enskog prediction is about 8 percent higher");
+