summaryrefslogtreecommitdiff
path: root/629/CH9/EX9.6/example9_6.sce
diff options
context:
space:
mode:
authorprashantsinalkar2017-10-10 12:27:19 +0530
committerprashantsinalkar2017-10-10 12:27:19 +0530
commit7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch)
treedbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /629/CH9/EX9.6/example9_6.sce
parentb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff)
downloadScilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip
initial commit / add all books
Diffstat (limited to '629/CH9/EX9.6/example9_6.sce')
-rw-r--r--629/CH9/EX9.6/example9_6.sce30
1 files changed, 30 insertions, 0 deletions
diff --git a/629/CH9/EX9.6/example9_6.sce b/629/CH9/EX9.6/example9_6.sce
new file mode 100644
index 000000000..d3d30928f
--- /dev/null
+++ b/629/CH9/EX9.6/example9_6.sce
@@ -0,0 +1,30 @@
+clear
+clc
+//Example 9.6 LAMINAR/TURBULENT BOUNDARY LAYER ON FLAT PLATE
+Uo=30; //[m/s]
+L=3; //[m]
+v=1.51*10^-5; //viscosity [m^2/s]
+Re=Uo*L/v //Reynolds number
+
+//Average shear-stress coefficient
+Cf=(0.523/(log(0.06*Re))^2)-(1520/Re)
+printf("\nAverage shear-stress coefficient, Cf, for the plate = %.5f.\n\n",Cf)
+
+//Total shear force
+B=1; //[m]
+rho=1.2; //[kg/m^3]
+Fs=Cf*B*L*rho*(Uo^2/2) //[N]
+printf("\nTotal shear force on one side of plate = %.2f N.\n\n",Fs)
+
+Re_tr=5*10^5;
+xtr=Re_tr*v/Uo //transition point [m]
+//Laminar average shear-stress coefficient
+Cfl=1.33/Re_tr^(1/2)
+
+//Laminar shear force
+Fsl=(Cfl*B*xtr*rho*Uo^2)/2 //[N]
+printf("\nShear force due to laminar part = %.3f N.\n\n",Fsl)
+
+//Turbulent shear force
+Fst=Fs-Fsl //[N]
+printf("\nShear force due to turbulent part = %.2f N.\n",Fst) \ No newline at end of file