diff options
Diffstat (limited to '3720/CH10')
| -rw-r--r-- | 3720/CH10/EX10.11/Ex10_11.sce | 15 | ||||
| -rw-r--r-- | 3720/CH10/EX10.12/Ex10_12.sce | 34 | ||||
| -rw-r--r-- | 3720/CH10/EX10.15/Ex10_15.sce | 15 | ||||
| -rw-r--r-- | 3720/CH10/EX10.2/Ex10_2.sce | 16 | ||||
| -rw-r--r-- | 3720/CH10/EX10.6/Ex10_6.sce | 37 | ||||
| -rw-r--r-- | 3720/CH10/EX10.8/Ex10_8.sce | 12 | ||||
| -rw-r--r-- | 3720/CH10/EX10.9/Ex10_9.sce | 16 |
7 files changed, 145 insertions, 0 deletions
diff --git a/3720/CH10/EX10.11/Ex10_11.sce b/3720/CH10/EX10.11/Ex10_11.sce new file mode 100644 index 000000000..7d56fcd9d --- /dev/null +++ b/3720/CH10/EX10.11/Ex10_11.sce @@ -0,0 +1,15 @@ +// Example 10_11
+clc;clear;funcprot(0);
+//Given data
+T=19;// °C
+D=30/100;// Diameter in m
+x=30/100;// Length of the tunnel in m
+V_b=4.0;// Velocity at beginning in m/s
+nu=1.507*10^-5;// m^2/s
+
+// Calculation
+Re_x=(V_b*x)/nu;// Reynolds number
+delta=((1.72*x)/(sqrt(Re_x)))*10^3;// The displacement thickness at the end of the test section in mm
+R=D/2;// Radius of the tunnel in m
+V_end=(V_b*(%pi*R^2))/(%pi*(R-(delta/1000))^2);// The average air speed at the end of the test section in m/s
+printf('\nThe average air speed at the end of the test section=%0.2f m/s',V_end);
diff --git a/3720/CH10/EX10.12/Ex10_12.sce b/3720/CH10/EX10.12/Ex10_12.sce new file mode 100644 index 000000000..93ce9eb47 --- /dev/null +++ b/3720/CH10/EX10.12/Ex10_12.sce @@ -0,0 +1,34 @@ +// Example 10_12
+clc;clear;funcprot(0);
+//Given data
+V=10.0;// m/s
+L=1.52;// m
+
+//Properties
+nu=1.516*10^-5;// m^2/s
+
+//Calculation
+//(a)
+x=L;// m
+Re_x=(V*x)/nu;// Reynolds number
+L=L*1000;// mm
+x=[0,L];// mm
+
+//For laminar case
+for(i=1:2)
+del_laminar(i)=(4.91*x(i))/sqrt(Re_x);// mm
+del_turbulenta(i)=(0.16*x(i))/(Re_x)^(1/7);// mm
+del_turbulentb(i)=(0.38*x(i))/(Re_x)^(1/5);// mm
+end
+xlabel('x,m');
+ylabel('delta,mm');
+x=x/1000;
+plot(x,del_laminar,'b',x,del_turbulenta,'r',x,del_turbulentb,'g');
+legend(['Laminar','Turbulent(a)','Turbulent(b)'],"in_upper_left");
+//(b)
+// For laminar boundary layer,
+C_fxl=0.664/sqrt(Re_x);
+// For turbulent boundary layer,
+C_fxt=0.027/(Re_x)^(1/7);
+printf('\nThe laminar boundary layer thickness at this same x-location=%0.2f mm \nThe turbulent boundary layer thickness at this same x-location=%0.1f mm \nThe local skin friction coefficient for the laminar boundary layer=%0.2e \nThe local skin friction coefficient for the turbulent boundary layer=%0.1e',del_laminar(2),del_turbulenta(2),C_fxl,C_fxt);
+// The answer vary due to round off error
diff --git a/3720/CH10/EX10.15/Ex10_15.sce b/3720/CH10/EX10.15/Ex10_15.sce new file mode 100644 index 000000000..bb4e74590 --- /dev/null +++ b/3720/CH10/EX10.15/Ex10_15.sce @@ -0,0 +1,15 @@ +// Example 10_15
+clc;clear;funcprot(0);
+//Given data
+T=20;// °C
+L=1.8;// Length in m
+w=0.50;// Width in m
+U=10;// Velocity of the flow in m/s
+delta_1=4.2/100;// Boundary layer thickness 1 in m
+delta_2=7.7/100;// Boundary layer thickness 2 in m
+nu=1.516*10^-5;// m^2/s
+rho=1.204;// kg/m3
+
+// Calculation
+F_d=(w*rho*U^2)*(4/45)*(delta_2-delta_1);// Drag force in N
+printf('\nThe total skin friction drag force=%0.2f N',F_d);
diff --git a/3720/CH10/EX10.2/Ex10_2.sce b/3720/CH10/EX10.2/Ex10_2.sce new file mode 100644 index 000000000..78f68fd39 --- /dev/null +++ b/3720/CH10/EX10.2/Ex10_2.sce @@ -0,0 +1,16 @@ +// Example 10_2
+clc;clear;funcprot(0);
+//Given data
+D=50*10^-6;// Diameter of spherical ash particle in m
+T=-50;// °C
+P=55;// kPa
+rho_p=1240;// The density of the particle in kg/m^3
+//Properties
+mu=1.474*10^-5;// kg/m.s
+rho_air=0.8588;// kg/m^3
+g=9.81;// The acceleration due to gravity in m/s^2
+
+//Calculation
+V=(D^2/(18*mu))*(rho_p-rho_air)*g;//The terminal velocity of this particle in m/s
+printf('\nThe terminal velocity of this particle,V=%0.3f m/s',V);
+Re=(rho_air*V*D)/mu;
diff --git a/3720/CH10/EX10.6/Ex10_6.sce b/3720/CH10/EX10.6/Ex10_6.sce new file mode 100644 index 000000000..7dc720be9 --- /dev/null +++ b/3720/CH10/EX10.6/Ex10_6.sce @@ -0,0 +1,37 @@ +// Example 10_6
+clc;clear;
+//Given data
+// Assume (vdot/L)_1=V1,(vdot/L)_2=V2;
+V1=2.00;// m^2/s
+V2=-1.00;// m^2/s
+gamma1=1.50;// m^2/s
+x_1=0;
+y_1=1;
+x_2=1;
+y_2=-1;
+x=1.0;
+y=0;// where all spatial coordinates are in meters.
+
+//Calculation
+//From fig.10-53,The vortex is located 1 m above the point (1, 0) and vortex velocity has positive i direction
+r_vortex=1.00;// m
+V_vortex=[gamma1/(2*%pi*r_vortex) 0];// m/s
+//Similarly, the first source induces a velocity at point (1, 0) at a 45° angle from the x-axis as shown in Fig. 10–53.
+r_source1=sqrt(2);// m
+V_source1=(V1)/(2*%pi*r_source1);// Resultant vector in m/s
+theta=45;// angle between two vectors
+// Function to find the velocity vector in i and j direction from resultant vector
+ function [X]=fric(f)
+ X(1)=f(1)^2 + f(2)^2-V_source1^2; // modulus(r)=sqrt(x^2+y^2)
+ X(2)=tand(theta)*f(1)-f(2);// theta=tan^-1(y/x)
+ endfunction
+
+ f=[0.01 0.01]; // Initial guess to solve X
+ V_source1_vec=fsolve(f,fric);// m/s (Calculating friction factor)
+
+//Finally, the second source (the sink) induces a velocity straight down i.e in the negative j direction
+r_source2=1.00;/// m
+V_source2=[0 (V2)/(2*%pi*r_source2)];// m/s
+V=V_vortex+V_source1_vec+V_source2;//The resultant velocity in m/s
+printf('\nThe resultant velocity, V = %0.3fi %1.0fj\n',V);
+
diff --git a/3720/CH10/EX10.8/Ex10_8.sce b/3720/CH10/EX10.8/Ex10_8.sce new file mode 100644 index 000000000..a379ca0ec --- /dev/null +++ b/3720/CH10/EX10.8/Ex10_8.sce @@ -0,0 +1,12 @@ +//Example 10_8
+clc;clear;funcprot(0)
+// Given values
+w=2.0;// Width in mm
+L=35.0;// Length in cm
+b=2.0;// Distance in cm
+v_dot=0.110;// The total volume flow rate in m^3/s
+u_starmax=0.159;// m/s
+// Calculation
+v_dotbyL=-(v_dot/(L/100));// Strength of line source in m^2/s
+u_max=-(u_starmax*(v_dotbyL/(b/100)));// Maximum speed along the floor
+printf('\nStrength of line source=%0.3f m^2/s \nMaximum speed along the floor,u_max=%0.2f m/s',v_dotbyL,u_max);
diff --git a/3720/CH10/EX10.9/Ex10_9.sce b/3720/CH10/EX10.9/Ex10_9.sce new file mode 100644 index 000000000..845941787 --- /dev/null +++ b/3720/CH10/EX10.9/Ex10_9.sce @@ -0,0 +1,16 @@ +// Example 10_9
+clc;clear;funcprot(0);
+//Given data
+V=5.0;// Uniform speed in mi/h
+x=16;// Length in ft
+T=50;// °F
+nu=1.407*10^-5;// The kinematic viscosity of water in ft^2/s
+
+// Calculation
+Re_x=(V*x)/nu;// The Reynolds number at the stern of the canoe
+Re_cr=1*10^5;// Critical Reynolds number
+if(Re_x>Re_cr)
+ printf('\nThe boundary layer is definitely turbulent by the back of the canoe.');
+else
+ printf('\nThe boundary layer is definitely laminar');
+end
|