1 files changed, 38 insertions, 0 deletions
diff --git a/3557/CH6/EX6.13/Ex6_13.sce b/3557/CH6/EX6.13/Ex6_13.sce
new file mode 100644
index 000000000..e602d2229
--- /dev/null
+++ b/3557/CH6/EX6.13/Ex6_13.sce
@@ -0,0 +1,38 @@
+//Example 6.13//
+a=514;//K //Kelvin //Temperature
+b=273;//K //Kelvin //Temperature
+apt=a+b
+mprintf("apt = %i K for eta = 10^13.4P",apt)
+c=696;//K //Kelvin //Temperature
+spt=c+b//for eta=10^7.6P
+mprintf("\n spt = %i K",spt)
+i=(10^13.4); //P //Pascal //preexponential constant
+j=(10^7.6);//P // Pascal //preexponential constant
+f=8.314;//J/(mol K) //universal gas constant
+a1=log(i/j); //(Taking antilog of i and j to remove exponential term)
+//mprintf("\na1 = %f ",a1)
+b1=(1/apt)-(1/spt);//(subtracting the temperature)
+//mprintf("\nb1 = %e ",b1)
+Q=(a1/b1)*f
+mprintf("\nQ = %e J/mol  = 465kJ (As 1K = 10^3)",Q)
+eta0=i*%e^-((Q)/(f*apt))
+mprintf("\n eta0 = %e P",eta0)
+h=10^4;//given
+//for eta=10^4 P and eta=10^8 P
+//for eta = 10^4
+T=Q/((f)*log(h/eta0))
+mprintf("\n T = %i K = 858 degree C",T)
+//for eta=10^8P
+h1=10^8;//P //Pascal 
+T1=Q/((f)*log(h1/eta0))
+mprintf("\n T1 = %i K = 680 degree C",T1)
+//Therefore working range = 680 to 858 degree C
+
+//For melting range eta=50 to 500 P
+ eta=50;//P //Pascal
+T2=Q/((f)*log(eta/eta0))
+mprintf("\n T2 = %i K = 993 degree C",T2)
+ eta1 = 500;// P //Pascal
+T3=Q/((f)*log(eta1/eta0))
+mprintf("\n T3 = %i K = 931 degree C",T3)
+mprintf("\n melting range = 931 to 993 degree C")