1 files changed, 28 insertions, 0 deletions
diff --git a/965/CH4/EX4.16/16.sci b/965/CH4/EX4.16/16.sci
new file mode 100644
index 000000000..0440c8aac
--- /dev/null
+++ b/965/CH4/EX4.16/16.sci
@@ -0,0 +1,28 @@
+clc;
+clear all;
+disp("temperation variation")
+L=60/2000;//m
+k=42.6;// W/(m.C)
+a=0.043;// m^2/h
+ti=440;// degree C
+h=235;// W/(m^ 2.C)
+ta=50;// degree C
+tau=4.3*60;//seconds
+Lc=L;// characteristic length
+Fo=a*tau/Lc^2;
+Bi=h*Lc/k;
+if Bi>1
+disp("Internal temperature gradients are not small and can not be neglected")
+end
+Bii=1/Bi;
+x=0;// midplane
+//(to-ta)/(ti-ta)=0.6
+to=0.6*(ti-ta)+ta;
+disp("Degree C",to,"Temperature at midplane  Tm =")
+ x=0.015/L;
+ Bii=6.06;
+ ti=to;
+ //(to-ta)/(ti-ta)=0.97
+ to=0.97*(ti-ta)+ta;
+ disp("Degree C",to,"Temperature inside the plate 15mm from the midplane To =")
+