initial commit / add all books

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)