1 files changed, 31 insertions, 0 deletions
diff --git a/1073/CH3/EX3.1/3_1.sce b/1073/CH3/EX3.1/3_1.sce
new file mode 100755
index 000000000..12e938423
--- /dev/null
+++ b/1073/CH3/EX3.1/3_1.sce
@@ -0,0 +1,31 @@
+clc;
+clear;
+//Example 3.1
+mu=10^-3        //N.s/m^2
+//At distance y from surface
+//ux=a+by+cy^2+dy^3
+//At y=0,ux=0 therefore a=0
+//i.e tao=0
+//At edge of boundary layer,ie y=del
+//ux=u_inf
+//At y=o,c=0
+//At y=del,ux=b*del+d*del^3
+
+//Therefore, b=-3*d*del^3
+//d=-u_inf/(2*del^2)
+//b=3*u_inf/(2*del)
+
+//For velocity profile,we have:
+//del/x=4.64*(Nre_x)^(-1/2)
+
+//Evaluate N re_x
+
+x=75;        //[mm]
+x=x/1000;        //[m]
+u_inf=3;            //[m/s]
+rho=1000            //[kg/m^3] for air
+Nre_x=u_inf*rho*x/mu        //Reynold number
+//Substituting the value,we get
+del=x*4.64*(Nre_x^(-1/2))       //[m]
+printf("\nBoundary layer thickness is del=%f m or %f mm",del,del*1000);
+printf("\nWrong units in answer of book,m and mm are wrongly interchanged");