8 files changed, 145 insertions, 0 deletions
diff --git a/3557/CH5/EX5.1/Ex5_1.sce b/3557/CH5/EX5.1/Ex5_1.sce
new file mode 100644
index 000000000..9ca003014
--- /dev/null
+++ b/3557/CH5/EX5.1/Ex5_1.sce
@@ -0,0 +1,12 @@
+//Example 5.1//
+k=2.95*10^-4;// kg/(m^-4.s) //At 400 degree Celsius k rises
+k1=1.05*10^-8;//kg/(m^-4.s) // The value of k at 300 degree celsius
+R=8.314;//J/(mol.K) //universal gas constant
+T=673;//K //Kelvin //absolute temperature
+T1=573;//K //Kelvin //absolute temperature
+a=log(k/k1);// Taking antilog to remove exponential term
+//mprintf("a=%e ",a)
+c=(1/T)-(1/T1); //subtracting the term
+//mprintf("c = %e ",c)
+Q=(-(a/c))*R //cross multiplication of the term
+mprintf("Q = %e J/mol = 328 kJ/mol",Q)
diff --git a/3557/CH5/EX5.2/Ex5_2.sce b/3557/CH5/EX5.2/Ex5_2.sce
new file mode 100644
index 000000000..8ff2c18c7
--- /dev/null
+++ b/3557/CH5/EX5.2/Ex5_2.sce
@@ -0,0 +1,17 @@
+//Example 5.2//
+nv=2.29*10^-5;//the fraction of aluminium lattice sites vacant 
+Ev=0.76;//eV //elevtrom volts
+k=86.2*10^-6;//eV //Boltzmann's constant
+T=673;//K //Kelvin //absolute temperature
+T1=933;// K //Kelvin //absolute temperature
+//At 400degree C(=673K)
+a=Ev/(k*T)// solving the exponential raise to equation
+//mprintf("a = %f ",a)
+C=nv*%e^a
+mprintf("C = %f",C)
+//At 660 degree C (=993K)
+b=Ev/(k*T1)//solving the exponential raise to equation
+//mprintf("b = %f ",b)
+N=C*%e^-b
+mprintf("\nN = %e ",N)
+mprintf("\nor roughly nine vacancies occur for every 10,000 lattice sites ")
diff --git a/3557/CH5/EX5.3/Ex5_3.sce b/3557/CH5/EX5.3/Ex5_3.sce
new file mode 100644
index 000000000..407456640
--- /dev/null
+++ b/3557/CH5/EX5.3/Ex5_3.sce
@@ -0,0 +1,28 @@
+//Example 5.3//
+c1=5;//at % //drop in carbon concentration 5 to 4 at %
+c2=4;// at % //drop in carbon concentration 5 to 4 at %
+x1=1;//mm //millimetre
+x2=2;//mm //millimetre
+d=(c1-c2)/(x1-x2)
+mprintf("d = %i at percent /mm",d)
+a=7.63;//g/cm^3 //gram per cubic centimeter
+b=0.6023*10^24;//atoms //Avgardo's number
+c=55.85;//g //atomic mass of iron (from appendix 1)
+p=a*(b/c)
+mprintf("\np = %e atoms/cm^3",p)
+a1=0.01;//given
+c1=1;//mm //millimetre
+d1=10^6;//cm^3/m^3
+e1=10^3;//mm/m
+d2=-((a1*p)*c1)*(d1)*(e1)
+mprintf("\nd2 = %e atoms/m^4",d2)
+D0=20*10^-6;//m^2/s //preexponential constant
+Q=142000;//J/mol //activation energy for defect motion
+R=8.314;//J/mol/K //universal gas constant
+T=1273;//K //Kelvin // absolute temperature
+Dc=D0*(%e^-(Q/(R*T)))
+mprintf("\nDc = %e m^2/s",Dc)
+c2=(-8.23*10^29);//atoms/m^4
+J=-Dc*c2
+mprintf("\nJ =%e atoms/(m^2.s)",J )
+
diff --git a/3557/CH5/EX5.4/Ex5_4.sce b/3557/CH5/EX5.4/Ex5_4.sce
new file mode 100644
index 000000000..8919c8d19
--- /dev/null
+++ b/3557/CH5/EX5.4/Ex5_4.sce
@@ -0,0 +1,19 @@
+//Example 5.4//
+cx=0.5;//carbon content
+b=1;//given
+e=b-cx
+mprintf("e = %f ",e)
+c=0.4755;//As z= 0.45 therefore erf (z) is obtained //Interpolating table 5.1 gives
+d=0.5205;//As z=0.50 therefore erf(z) is obtained //Interpolating table 5.1 gives
+g=0.45;//given
+z=(((e-c)/(d-c))*(e-g))+g
+mprintf("\nz = %f",z)
+x=1*10^-3;//Using the diffusivity from sample problem 5.3
+D=2.98*10^-11;//m^2/s //Arrhenius equation
+m=(x^2)/(4*(z^2)*D)
+//mprintf("\nm = %e ",m)
+i=1;//h //hour
+j=3.6*10^3;//s //second
+t=m*(i/j)
+mprintf("\nt = %f h",t)
+
diff --git a/3557/CH5/EX5.5/Ex5_5.sce b/3557/CH5/EX5.5/Ex5_5.sce
new file mode 100644
index 000000000..36daf22c6
--- /dev/null
+++ b/3557/CH5/EX5.5/Ex5_5.sce
@@ -0,0 +1,10 @@
+//Example 5.5//
+x=1*10^-3;//m// Using the diffusivity from sample problem 5.3
+D=2.98*10^-11;//m^2/s //arrhenius equations
+a=0.95;//from the figure 5.11
+d=(x^2)/((a^2)*(D))// calculating the value of d
+mprintf("d = %e h",d)
+b=1;//h //hour
+c=3.6*10^3;//s //second
+t=d*(b/c)
+mprintf("\nt = %f h",t)
diff --git a/3557/CH5/EX5.6/Ex5_6.sce b/3557/CH5/EX5.6/Ex5_6.sce
new file mode 100644
index 000000000..bc065f965
--- /dev/null
+++ b/3557/CH5/EX5.6/Ex5_6.sce
@@ -0,0 +1,14 @@
+//Example 5.6//
+x=0.75*10^-3;//m //meter //given
+t=3.6*10^4;//s //seconds //time
+a=0.95;//given
+D=(x^2)/((a^2)*(t))
+mprintf("D = %e m^2/s",D)
+b=20*10^-6;//m^2/s //preexponential constant
+c=142000;//J/mol //activation energy for defect motion
+d=8.314;//J/(mol.K)//universal gas constant 
+e=c/d
+//mprintf("\ne = %e",e)
+y=(-log(D/b))
+T1=inv(y/e)
+mprintf("\nT1 = %i K  = 952 degree C",T1)
diff --git a/3557/CH5/EX5.7/Ex5_7.sce b/3557/CH5/EX5.7/Ex5_7.sce
new file mode 100644
index 000000000..06e98d276
--- /dev/null
+++ b/3557/CH5/EX5.7/Ex5_7.sce
@@ -0,0 +1,14 @@
+//Example 5.7//
+D=(-1.0*10^-8);//m^2/s //constant diffusion coefficient
+ch=1.5;//kg/m^3 //constant surface concentrationsof the diffusing species
+ct=0.3;//kg/m^3 //constant surface concentrationsof the diffusing species
+x=5*10^-3;//m //meter //solid of thickness
+y=(-D)*(((ch-ct)/(x)))
+//mprintf("y = %e kg/m^2 h",y)
+t=3.6*10^3;//s/h //time
+J=y*t
+mprintf("J = %e kg/m^2.h",J)
+//The total mass of hydrogen being purified will then be this flux times the membrane area
+A=0.2;//m^2 //membrane area
+m=J*A
+mprintf("\nm = %e kg/h",m)
diff --git a/3557/CH5/EX5.8/Ex5_8.sce b/3557/CH5/EX5.8/Ex5_8.sce
new file mode 100644
index 000000000..b5dc7aff6
--- /dev/null
+++ b/3557/CH5/EX5.8/Ex5_8.sce
@@ -0,0 +1,31 @@
+//Example 5.8//
+
+cx=0.01;// distance of x
+c0=0;////for initially pure A
+c=cx-c0
+mprintf("c = %f ",c)
+a=1;//given
+e=a-c
+mprintf("\ne = %f ",e)
+b=0.9928;//As z= 1.90 erf(z)=0.9928 //Interpolating table 5.1 gives
+d=0.99;//Interpolating table 5.1 gives
+f=0.9891;//As z=1.80 erf(z)=0.9891 //Interpolating table 5.1 gives
+h=1.90;//given
+i=1.80;//given
+z=-((((b-d)/(b-f))*(h-i))-h)
+mprintf("\nz = %f ",z)
+D=1*10^-10;//m^2/s// grain boundary
+D1=1*10^-14;//m^2/s // volume of bulk grain
+t=1;//h //hour //time
+t1=3.6*10^3;//s/h //time
+x=2*z*sqrt(D*t*t1)
+mprintf("\nx = %e m ",x)
+a1=10^3;//(As 1milli = 10^-3)
+a2=a1*x
+mprintf(" = %f mm",a2)
+//(b) For comparison
+x1=2*z*sqrt(D1*t*t1)
+mprintf("\nx1 = %e m ",x1)
+b1=10^6;//(As mew = 10^-6)
+b2=b1*x1
+mprintf(" = %f mew m",b2)