From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 587/CH7/EX7.2/example7_2.sce | 44 ++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 44 insertions(+) create mode 100755 587/CH7/EX7.2/example7_2.sce (limited to '587/CH7/EX7.2') diff --git a/587/CH7/EX7.2/example7_2.sce b/587/CH7/EX7.2/example7_2.sce new file mode 100755 index 000000000..780a50018 --- /dev/null +++ b/587/CH7/EX7.2/example7_2.sce @@ -0,0 +1,44 @@ +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