3 files changed, 102 insertions, 109 deletions
diff --git a/587/CH7/EX7.2/example7_2.sce b/587/CH7/EX7.2/example7_2.sce
index 780a50018..3e80d6250 100755
--- a/587/CH7/EX7.2/example7_2.sce
+++ b/587/CH7/EX7.2/example7_2.sce
@@ -1,44 +1,35 @@
-clear;
-clc;
-
-//Example7.2[Cooling of a Hot Block by Forced Air at High Elevation]
-//Given:-
-ReC=5*10^5;//critical Reynolds number
-v=8;//Velocity of air[m/s]
-T_air=20;//Initial Temperature of air[degree Celcius]
-T_plate=140;//Temperature of flat plate[degree Celcius]
-T_film=(T_air+T_plate)/2;//Film Temperature of air[degree Celcius]
-//Properties of air at film temperature[degree Celcius]
-k=0.02953;//[W/m.degree Celcius]
-Pr=0.7154;//Prandtl Number
-nu=2.097*10^(-5);//Kinematic Viscosity at 1 atm Pressure [m^2/s]
-nu_ac=nu/0.823;//Kinematic viscosity at pressure 0.823 atm[m^2/s]
-//Solution(a)
-L1=6;//Characteristic length of plate along the flow of air[m]
-w1=1.5;//width[m]
-ReL1=(v*L1)/nu_ac;//Reynolds number
-if(ReL1>ReC) then,
-    disp("Flow is not laminar")
-    //We have average Nusselt Number
-    Nu1=((0.037*(ReL1^(0.8)))-871)*(Pr^(1/3));
-    disp(ceil(Nu1),"Nusselt Number is")
-    h1=k*Nu1/L1;//[W/m^2.degree Celcius]
-    As1=w1*L1;//Flow Area of plate[m^2]
-    Q1=h1*As1*(T_plate-T_air);
-    disp("W",Q1,"Heat Flow Rate is")
-else,
-    disp("Flow is laminar")
-end
-//Solution(b)
-L2=1.5;//Characteristic length of plate along flow of air[m]
-ReL2=v*L2/nu_ac;//Reynolds Number
-if(ReL2<Rec) then,
-    disp("Flow is laminar")
-    Nu2=0.664*(ReL2^(0.5))*(Pr^(1/3));
-    disp(ceil(Nu2),"Nusselt Number is")
-    h2=k*Nu2/L2;//[W/m^2.degree Celcius]
-    Q2=h2*As1*(T_plate-T_air);
-    disp("W",ceil(Q2),"The heat transfer rate is")
-else,
-    disp("Flow is turbulent")
-end
+clear;
+clc;
+
+//Example7.2[Cooling of a Hot Block by Forced Air at High Elevation]
+//Given:-
+ReC=5*10^5;//critical Reynolds number
+v=8;//Velocity of air[m/s]
+T_air=20;//Initial Temperature of air[degree Celcius]
+T_plate=140;//Temperature of flat plate[degree Celcius]
+T_film=(T_air+T_plate)/2;//Film Temperature of air[degree Celcius]
+//Properties of air at film temperature[degree Celcius]
+k=0.02953;//[W/m.degree Celcius]
+Pr=0.7154;//Prandtl Number
+nu=2.097*10^(-5);//Kinematic Viscosity at 1 atm Pressure [m^2/s]
+nu_ac=nu/0.823;//Kinematic viscosity at pressure 0.823 atm[m^2/s]
+//Solution(a)
+L1=6;//Characteristic length of plate along the flow of air[m]
+w1=1.5;//width[m]
+ReL1=(v*L1)/nu_ac;//Reynolds number
+if(ReL1&gt;ReC) then,
+    disp("Flow is not laminar")
+    //We have average Nusselt Number
+    Nu1=((0.037*(ReL1^(0.8)))-871)*(Pr^(1/3));
+    disp(ceil(Nu1),"Nusselt Number is")
+    h1=k*Nu1/L1;//[W/m^2.degree Celcius]
+    As1=w1*L1;//Flow Area of plate[m^2]
+    Q1=h1*As1*(T_plate-T_air);
+    disp("W",Q1,"Heat Flow Rate is")
+else,
+    disp("Flow is laminar")
+end
+//Solution(b)
+L2=1.5;//Characteristic length of plate along flow of air[m]
+ReL2=v*L2/nu_ac;//Reynolds Number
+if(ReL2<ReC) then,="" disp("flow="" is="" laminar")="" nu2="0.664*(ReL2^(0.5))*(Pr^(1/3));" disp(ceil(nu2),"nusselt="" number="" is")="" h2="k*Nu2/L2;//[W/m^2.degree" celcius]="" q2="h2*As1*(T_plate-T_air);" disp("w",ceil(q2),"the="" heat="" transfer="" rate="" else,="" turbulent")="" end="" <="" div=""></rec)>
+\ No newline at end of file
diff --git a/587/CH7/EX7.6/example7_6.sce b/587/CH7/EX7.6/example7_6.sce
index eab4b82ae..9bb91d08f 100755
--- a/587/CH7/EX7.6/example7_6.sce
+++ b/587/CH7/EX7.6/example7_6.sce
@@ -1,31 +1,32 @@
-clear;
-clc;
-
-//Example7.6[Cooling of a Steel Ball by Forced Air]
-//Given:-
-rho=8055;//[kg/m^3]
-Cp=480;//[J/kg.degree Celcius]
-To=300;//Temp of oven[degree Celcius]
-Ta=25;//Temp of air[degree Celcius]
-va=3;//Velocity of air[m/s]
-Ts=200;//Dropped temp of surface of ball[degree Celcius]
-Ts_avg=(Ts+To)/2;//[degree Celcius]
-d=0.25;//[m]
-mu_s=2.76*10^(-5);//Dynamic Viscosity at average surface temperature[kg/m.s]
-//Properties of air at 25 degree Celcius
-k=0.02551;//[W/m.degree Celcius]
-nu=1.562*10^(-5);//kinematic viscosity[m^2/s]
-mu=1.849*10^(-5);//Dynamic viscosity of air at 25 degree C[kg/m.s]
-//Solution:-
-Re=va*d/nu;//[Reynolds Number]
-Nu=2+[(0.4*(Re^(1/2)))+(0.06*(Re^(2/3)))]*(Pr^(0.4))*((mu/mu_s)^(1/4));
-disp(ceil(Nu),"The Nusselt number is")
-h=k*Nu/d;//[W/m^2.degree Celcius]
-As=%pi*(d^2);//[m^2]
-Q_avg=h*As*(Ts_avg-Ta);//[W]
-disp("W",ceil(Q_avg),"The average rate of heat transfer from Newtons Law of cooling is")
-m=rho*%pi*(d^3)/6;//[kg]
-Q_total=m*Cp*(To-Ts);//[J]
-disp("J",Q_total,"The total heat transferred from the ball is")
-delta_t=Q_total/Q_avg;//[s]
+clear;
+clc;
+
+//Example7.6[Cooling of a Steel Ball by Forced Air]
+//Given:-
+rho=8055;//[kg/m^3]
+Pr = 0.7296;
+Cp=480;//[J/kg.degree Celcius]
+To=300;//Temp of oven[degree Celcius]
+Ta=25;//Temp of air[degree Celcius]
+va=3;//Velocity of air[m/s]
+Ts=200;//Dropped temp of surface of ball[degree Celcius]
+Ts_avg=(Ts+To)/2;//[degree Celcius]
+d=0.25;//[m]
+mu_s=2.76*10^(-5);//Dynamic Viscosity at average surface temperature[kg/m.s]
+//Properties of air at 25 degree Celcius
+k=0.02551;//[W/m.degree Celcius]
+nu=1.562*10^(-5);//kinematic viscosity[m^2/s]
+mu=1.849*10^(-5);//Dynamic viscosity of air at 25 degree C[kg/m.s]
+//Solution:-
+Re=va*d/nu;//[Reynolds Number]
+Nu=2+[(0.4*(Re^(1/2)))+(0.06*(Re^(2/3)))]*(Pr^(0.4))*((mu/mu_s)^(1/4));
+disp(ceil(Nu),"The Nusselt number is")
+h=k*Nu/d;//[W/m^2.degree Celcius]
+As=%pi*(d^2);//[m^2]
+Q_avg=h*As*(Ts_avg-Ta);//[W]
+disp("W",ceil(Q_avg),"The average rate of heat transfer from Newtons Law of cooling is")
+m=rho*%pi*(d^3)/6;//[kg]
+Q_total=m*Cp*(To-Ts);//[J]
+disp("J",Q_total,"The total heat transferred from the ball is")
+delta_t=Q_total/Q_avg;//[s]
 disp("hour",delta_t/3600,"The time of cooling is")
 \ No newline at end of file
diff --git a/587/CH7/EX7.8/example7_8.sce b/587/CH7/EX7.8/example7_8.sce
index 736b5f6c2..136c21cd3 100755
--- a/587/CH7/EX7.8/example7_8.sce
+++ b/587/CH7/EX7.8/example7_8.sce
@@ -1,35 +1,36 @@
-clear;
-clc;
-
-//Example7.8[Effect of insulation on Surface Temperature]
-//Given:-
-Ti=120;//Initial temp of hot water[degree Celcius]
-k_pipe=15;//W/m.degree Celcius 
-ri=0.008,ro=0.01;//Inner and outer radii[m]
-t=0.002;//Thickness of pipe[m]
-To=25;//Ambient temperature[degree Celcius]
-Ts=40;//Maximum Temp of outer surface of insulation[degree Celcius]
-hi=70,ho=20;//Heat transfer coefficients inside and outside of the pipe[W/m^2.degree Celcius]
-k_insu=0.038;//[W/m.degree Celcius]
-L=1;//section of pipe[m]
-//Solution:-
-//Areas of surfaces exposed to convection
-A1=2*%pi*ri*L;//[m^2]
-//Individual Thermal Resistances
-R_conv1=1/(hi*A1);//[degree Celcius/W]
-R_pipe=(log(ro/ri))/(2*%pi*k_pipe*L);//[degree Celcius/W]
-//R_insu=(log(r3/ri))/(2*%pi*k_insu*L)
-//R_conv2=1/(ho*2*%pi*r3*L)
-//R_total=R_conv1+R_conv2+R_pipe+R_insu
-//Q=(Ti-To)/R_total;
-//Q=(Ts-To)/R_conv2;
-//Equating both Q we get
-function[r]=radius(r3)
-    r(1)=1884*r3(1)*(0.284+0.0024+4.188*log((r3(1))/0.01)+(1/(125.6*r3(1))))-95;
-    deff('[r]=radius(r3)',['radius_3=1884*r3(1)*(0.284+0.0024+4.188*log((r3(1))/0.01)+(1/(125.6*r3(1))))-95'])
-endfunction
-    disp("m",xs,"The outer radius of the insulation is")
-    t=xs-ro;//[m]
-    disp("cm",100*t,"The minimum thickness of fibreglass insulation required is")
-    ///Correct output will be displayed after executing the codes once and then re-executin them
-    
-\ No newline at end of file
+clear;
+clc;
+
+//Example7.8[Effect of insulation on Surface Temperature]
+//Given:-
+Ti=120;//Initial temp of hot water[degree Celcius]
+k_pipe=15;//W/m.degree Celcius 
+ri=0.008,ro=0.01;//Inner and outer radii[m]
+t=0.002;//Thickness of pipe[m]
+To=25;//Ambient temperature[degree Celcius]
+Ts=40;//Maximum Temp of outer surface of insulation[degree Celcius]
+hi=70,ho=20;//Heat transfer coefficients inside and outside of the pipe[W/m^2.degree Celcius]
+k_insu=0.038;//[W/m.degree Celcius]
+L=1;//section of pipe[m]
+//Solution:-
+//Areas of surfaces exposed to convection
+A1=2*%pi*ri*L;//[m^2]
+//Individual Thermal Resistances
+R_conv1=1/(hi*A1);//[degree Celcius/W]
+R_pipe=(log(ro/ri))/(2*%pi*k_pipe*L);//[degree Celcius/W]
+//R_insu=(log(r3/ri))/(2*%pi*k_insu*L)
+//R_conv2=1/(ho*2*%pi*r3*L)
+//R_total=R_conv1+R_conv2+R_pipe+R_insu
+//Q=(Ti-To)/R_total;
+//Q=(Ts-To)/R_conv2;
+//Equating both Q we get
+function[r]=radius(r3)
+    r(1)=1884*r3(1)*(0.284+0.0024+4.188*log((r3(1))/0.01)+(1/(125.6*r3(1))))-95;
+    deff('[r]=radius(r3)',['radius_3=1884*r3(1)*(0.284+0.0024+4.188*log((r3(1))/0.01)+(1/(125.6*r3(1))))-95'])
+endfunction
+x0=[1] 
+[xs,fxs,m]=fsolve(x0',radius)
+    disp("m",xs,"The outer radius of the insulation is")
+    t=xs-ro;//[m]
+    disp("cm",100*t,"The minimum thickness of fibreglass insulation required is")
+    ///Correct output will be displayed after executing the codes once and then re-executin them
+\ No newline at end of file