From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001 From: prashantsinalkar Date: Tue, 10 Oct 2017 12:27:19 +0530 Subject: initial commit / add all books --- 3720/CH10/EX10.12/Ex10_12.sce | 34 ++++++++++++++++++++++++++++++++++ 1 file changed, 34 insertions(+) create mode 100644 3720/CH10/EX10.12/Ex10_12.sce (limited to '3720/CH10/EX10.12/Ex10_12.sce') 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 -- cgit