8 files changed, 86 insertions, 0 deletions
diff --git a/3159/CH9/EX9.1/Ex9_1.sce b/3159/CH9/EX9.1/Ex9_1.sce
new file mode 100755
index 000000000..5c4bff11e
--- /dev/null
+++ b/3159/CH9/EX9.1/Ex9_1.sce
@@ -0,0 +1,25 @@
+// Calculate the critical free energy of nucleation of ice from water and critical radius
+clc
+del_t1 = 0// temperature difference in degree Celsius
+del_t2 = -5 // temperature difference in degree Celsius
+del_t3 = -40 // temperature difference in degree Celsius
+del_h = 6.02 // enthalpy of fusion in kJ/mol
+T_m = 273 // mean temperature
+Gamma = 0.076 // energy of ice water interface in J /m^2
+v = 19 // molar volume of ice
+printf("\n Example 9.1")
+printf("\n Part A")
+printf("\n At del_t = %d, there is no supercooling. So there is no critical radius",del_t1)
+printf("\n\n Part B")
+del_f = 16/3*%pi*(Gamma)^3*T_m^2/((del_h*1e3*1e6/v)^2*del_t2^2)
+r = 2*T_m*Gamma/(-del_h*1e3*1e6/v*del_t2)
+printf("\n Critical free energy of nucleation is %.1eJ",del_f)
+printf("\n Critical radius is %d angstrom",ceil(r*1e10))
+printf("\n\n Part C")
+temp_r = del_t3/del_t2
+del_f_ = del_f/temp_r^2
+r_ = r/temp_r
+
+ printf("\n Critical free energy of nucleation is %.1eJ",del_f_)
+printf("\n Critical radius is %d angstrom.", ceil(r_*1e10))
+
diff --git a/3159/CH9/EX9.1/Ex9_1.txt b/3159/CH9/EX9.1/Ex9_1.txt
new file mode 100755
index 000000000..e6ecccede
--- /dev/null
+++ b/3159/CH9/EX9.1/Ex9_1.txt
@@ -0,0 +1,12 @@
+
+ Example 9.1
+ Part A
+ At del_t = 0, there is no supercooling. So there is no critical radius
+
+ Part B
+ Critical free energy of nucleation is 2.2e-16J
+ Critical radius is 262 angstrom
+
+ Part C
+ Critical free energy of nucleation is 3.4e-18J
+ Critical radius is 33 angstrom. 
+\ No newline at end of file
diff --git a/3159/CH9/EX9.2/Ex9_2.sce b/3159/CH9/EX9.2/Ex9_2.sce
new file mode 100755
index 000000000..0b6e24c20
--- /dev/null
+++ b/3159/CH9/EX9.2/Ex9_2.sce
@@ -0,0 +1,15 @@
+// Calculate the change in del_f required to increase nucleation rate
+clc
+T= 300 // temperature in kelvin
+R = 8.314 // universal gas constant
+k = 2.303 // conversion factor
+a1 = 1e42
+a2 = 1e6 // nucleation rate
+a3 = 1e10 // nucleation rate
+printf("Example 9.2")
+I1 = (log10(a1)-log10(a2))*k*R*T // exponent factor
+I2 = (log10(a1)-log10(a3))*k*R*T// exponent factor
+del_f = I1-I2 // difference 
+a = 10^(log10(a3)-log10(a2))
+
+printf("\n  A change of %d KJ mol^-1 energy is required to increase nucleation factor from \n  %.0e m^-3 s^-1  to  %.0e m^-3 s^-1 ",ceil(del_f/1e3),a,a3)
diff --git a/3159/CH9/EX9.2/Ex9_2.txt b/3159/CH9/EX9.2/Ex9_2.txt
new file mode 100755
index 000000000..dbf59fd3d
--- /dev/null
+++ b/3159/CH9/EX9.2/Ex9_2.txt
@@ -0,0 +1,3 @@
+ Example 9.2
+  A change of 23 KJ mol^-1 energy is required to increase nucleation factor from 
+  1e+04 m^-3 s^-1  to  1e+10 m^-3 s^-1  
+\ No newline at end of file
diff --git a/3159/CH9/EX9.4/Ex9_4.sce b/3159/CH9/EX9.4/Ex9_4.sce
new file mode 100755
index 000000000..7dbc981a9
--- /dev/null
+++ b/3159/CH9/EX9.4/Ex9_4.sce
@@ -0,0 +1,13 @@
+// Calculate del_f_het as a fraction of del_f_homo
+clc
+Gamma_alpha_del = 0.5 // in J m^-2
+Gamma_alpha_beta = 0.5 // in J m^-2
+Gamma_beta_del = 0.01 // in J m^-2
+
+printf("\n Example 9.4") 
+theta = acos((Gamma_alpha_del -Gamma_beta_del)/Gamma_alpha_beta)
+del_f_ratio = 1/4*(2-3*cos(theta)+(cos(theta))^3)
+
+
+printf("\n del_f_het is %0.4fth fraction of  del_f_homo.",del_f_ratio)
+
diff --git a/3159/CH9/EX9.4/Ex9_4.txt b/3159/CH9/EX9.4/Ex9_4.txt
new file mode 100755
index 000000000..6f169cc76
--- /dev/null
+++ b/3159/CH9/EX9.4/Ex9_4.txt
@@ -0,0 +1,3 @@
+ 
+ Example 9.4
+ del_f_het is 0.0003th fraction of  del_f_homo. 
+\ No newline at end of file
diff --git a/3159/CH9/EX9.6/Ex9_6.sce b/3159/CH9/EX9.6/Ex9_6.sce
new file mode 100755
index 000000000..521afedd9
--- /dev/null
+++ b/3159/CH9/EX9.6/Ex9_6.sce
@@ -0,0 +1,12 @@
+// Calculate the free energy change during recrystallization
+clc
+mu = 45.5e9
+b = 2.55e-10
+n1 = 1e9 // initial dislocation density
+n2 = 1e13 // final dislocation density
+printf("\n Example 9.6") 
+E = 1/2*mu*b^2*n2
+del_g = E // as difference between initial and final dislocation energy is four order magnitude
+printf("\n Free energy change during recrystallization is %d J m^-3",-del_g)
+// Numerical value of answer in book is 14800
+
diff --git a/3159/CH9/EX9.6/Ex9_6.txt b/3159/CH9/EX9.6/Ex9_6.txt
new file mode 100755
index 000000000..028c84ef6
--- /dev/null
+++ b/3159/CH9/EX9.6/Ex9_6.txt
@@ -0,0 +1,3 @@
+ 
+ Example 9.6
+ Free energy change during recrystallization is -14793 J m^-3 
+\ No newline at end of file