summaryrefslogtreecommitdiff
path: root/539/CH5/EX5.5/Example_5_5.sce
diff options
context:
space:
mode:
Diffstat (limited to '539/CH5/EX5.5/Example_5_5.sce')
-rwxr-xr-x539/CH5/EX5.5/Example_5_5.sce23
1 files changed, 23 insertions, 0 deletions
diff --git a/539/CH5/EX5.5/Example_5_5.sce b/539/CH5/EX5.5/Example_5_5.sce
new file mode 100755
index 000000000..f16339a69
--- /dev/null
+++ b/539/CH5/EX5.5/Example_5_5.sce
@@ -0,0 +1,23 @@
+//Diffusion Coefficient Activation Energy and Preexponential Calculations
+
+clear;
+clc;
+
+printf("\tExample 5.5\n");
+
+//From graph log D ad 1/T are deducted
+inv_T1=0.8*10^-3; //Reciprocal of temp. in K^-1
+inv_T2=1.1*10^-3; //Reciprocal of temp. in K^-1
+logD1=-12.4;
+logD2=-15.45;
+
+R=8.31; //Gas law Constant in J/mol-K
+
+Qd=-2.3*R*(logD1-logD2)/(inv_T1-inv_T2);
+printf("\nActivation energy is %d kJ/mol",Qd/1000);
+
+//For calculating Peexponential factor
+D0=10^(logD2+(Qd*inv_T2/(2.3*R)));
+printf("\nPreexponential factor D0 is %.1f * 10^-5 m^2/s\n",D0/10^-5);
+
+//End \ No newline at end of file