initial commit / add all books

3 files changed, 82 insertions, 0 deletions

diff --git a/1309/CH4/EX4.3/Figure4_3.jpeg b/1309/CH4/EX4.3/Figure4_3.jpeg
new file mode 100755
index 000000000..ff1b509de
--- /dev/null
+++ b/1309/CH4/EX4.3/Figure4_3.jpeg
diff --git a/1309/CH4/EX4.3/Result4_3.pdf b/1309/CH4/EX4.3/Result4_3.pdf
new file mode 100755
index 000000000..2f20eac4b
--- /dev/null
+++ b/1309/CH4/EX4.3/Result4_3.pdf
diff --git a/1309/CH4/EX4.3/ch4_3.sce b/1309/CH4/EX4.3/ch4_3.sce
new file mode 100755
index 000000000..b78cb8e36
--- /dev/null
+++ b/1309/CH4/EX4.3/ch4_3.sce
@@ -0,0 +1,82 @@
+clc;

+clear;

+printf("\t\t\tChapter4_example3\n\n\n");

+hc=30;

+L=0.24;

+k=1.25;

+c=890;

+rou=550;

+Bi=hc*L/k;

+alpha=k/(rou*c);

+printf("The value of diffusivity is %.2e sq.m/s",alpha);

+Tc=150;

+T_inf=600;

+T_i=50;

+printf("\nThe Biot number is %.2f,",Bi);

+if Bi<0.1 then

+    n=0;

+else if  Bi>0.1 then

+        n=1;

+    end

+end

+select n

+case 0 then 

+    disp('The Lumped capacity approach is applicable');

+case 1 then

+    disp('Since value of Biot number is greater than 0.1, Lumped capacity approach would not give accurate results, so figure 4.6 is to be used');

+    reciprocal_Bi=1/Bi;

+    dimensionless_temp=(Tc-T_inf)/(T_i-T_inf);

+    Fo=0.4; //the value of Fourier Number from figure 4.6(a)

+    t=L^2*Fo/alpha;

+    printf("The required time is %d s = %.1f hr",t,t/3600);

+end

+// reading values of dimensionless temperature from figure 4.6(b) using reciprocal of Biot number

+x_per_L=[0 0.2 0.4 0.6 0.8 0.9 1.0];

+[n,m]=size(x_per_L);

+printf("\nThe choosen values of x/L are: \n");

+disp(x_per_L);

+printf("\n Values for dimensionless temperature for corresponding values of x/L:")

+dim_T=[1.0 .97 .86 .68 .48 .36 .24]; // value for dimensionless temperature for corresponding value of x/L

+disp(dim_T);

+printf("the temperature profile with distance is\n");

+printf("\tx/L\t\t");

+for j=1:m

+    printf("%.2f\t",x_per_L(1,j));

+    

+end

+printf("\n");

+printf("(T-T_inf)/T_i-T_inf)\t");

+for i=1:m

+    printf("%.2f\t",dim_T(i));

+end

+T=zeros(1,m);

+x=zeros(1,m);

+for i=1:m

+    T(1,i)=dim_T(1,i)*dimensionless_temp*(T_i-T_inf)+T_inf;

+    x(1,i)=x_per_L(1,i)*L;

+end

+printf("\n\tx,cm\t\t");

+for i=1:m

+    X(1,i)=x(1,i)*100;

+    printf("%.1f\t",X(1,i));

+end

+printf("\nT, degree celsius\t");

+for i=1:m

+    printf("%d\t",T(1,i));

+end

+plot2d(X,T,rect=[0,0,24,600]);

+a=gca();

+newticks=a.x_ticks;

+newticks(2)=[0;4;8;12;16;20;24];

+newticks(3)=['0';'4';'8';'12';'16';'20';'24'];

+a.x_ticks=newticks;

+newticks1=a.y_ticks;

+newticks1(2)=[0;100;200;300;400;500;600];

+newticks1(3)=['0';'100';'200';'300';'400';'500';'600'];

+a.y_ticks=newticks1;

+xlabel('x,cm');

+ylabel('t,degree celsius');

+title('Temperature profile in the 24-cm slab after 2.5 hr.');

+filename='Temperature profile in the 24-cm slab after 2.5 hr.';

+xgrid(1);

+xs2jpg(0,filename);