From f35ea80659b6a49d1bb2ce1d7d002583f3f40947 Mon Sep 17 00:00:00 2001 From: prashantsinalkar Date: Tue, 10 Oct 2017 12:38:01 +0530 Subject: updated the code --- 746/CH9/EX9.01/9_01.sce | 8 -------- 746/CH9/EX9.04/9_04.sce | 41 ----------------------------------------- 746/CH9/EX9.05/9_05.sce | 19 ------------------- 746/CH9/EX9.06/9_06.sce | 17 ----------------- 746/CH9/EX9.07/9_07.sce | 8 -------- 746/CH9/EX9.08/9_08.sce | 28 ---------------------------- 746/CH9/EX9.09/9_09.sce | 20 -------------------- 7 files changed, 141 deletions(-) delete mode 100755 746/CH9/EX9.01/9_01.sce delete mode 100755 746/CH9/EX9.04/9_04.sce delete mode 100755 746/CH9/EX9.05/9_05.sce delete mode 100755 746/CH9/EX9.06/9_06.sce delete mode 100755 746/CH9/EX9.07/9_07.sce delete mode 100755 746/CH9/EX9.08/9_08.sce delete mode 100755 746/CH9/EX9.09/9_09.sce (limited to '746/CH9') diff --git a/746/CH9/EX9.01/9_01.sce b/746/CH9/EX9.01/9_01.sce deleted file mode 100755 index 755aef7b6..000000000 --- a/746/CH9/EX9.01/9_01.sce +++ /dev/null @@ -1,8 +0,0 @@ -//static pressure// -pathname=get_absolute_file_path('9.01.sce') -filename=pathname+filesep()+'9.01-data.sci' -exec(filename) -//Change in static pressure between sections 1 and 2: -C=(((L-2*d1)/(L-2*d2))^4-1)*100; -printf("\n\nRESULTS\n\n") -printf("\n\nChange in static pressure between the sections 1 and 2: %.3f percent \n\n",C) diff --git a/746/CH9/EX9.04/9_04.sce b/746/CH9/EX9.04/9_04.sce deleted file mode 100755 index db1fd7c5a..000000000 --- a/746/CH9/EX9.04/9_04.sce +++ /dev/null @@ -1,41 +0,0 @@ -//Displacement thickness and stress// -pathname=get_absolute_file_path('9.04.sce') -filename=pathname+filesep()+'9.04-data.sci' -exec(filename) -//Reynolds number: -ReL=U*L/v -//FOR TURBULENT FLOW -//Disturbance thickness(in m): -dL1=0.382/ReL^0.2*L -//Displacement thickness(in m): -function y=f(n),y=dL1*(1-n^(1/7)) -endfunction -dl1=intg(0,1,f) -//Skin friction coefficient: -Cf1=0.0594/ReL^0.2 -//Wall shear stress(in N/m^2): -tw1=Cf1*0.5*d*U^2 -//For LAMINAR FLOW: -//Disturbance thickness(in m) -dL2=5/sqrt(ReL)*L -//Displacement thickness(in m): -dl2=0.344*dL2 -//Skin friction coefficient: -Cf2=0.664/sqrt(ReL) -//Wall shear stress(in N/m^2): -tw2=Cf2*0.5*d*U^2 -//COMPARISON OF VALUES WITH LAMINAR FLOW -//Disturbance thickness -D=dL1/dL2 -//Displacement thickness -DS=dl1/dl2 -//Wall shear stress -WSS=tw1/tw2 -printf("\n\nRESULTS\n\n") -printf("\n\nDisturbace thickness: %.3f m\n\n",dL1) -printf("\n\nDisplacement thickness: %.3f m\n\n",dl1) -printf("\n\nWall shear stress: %f N/m^2\n\n",tw1) -printf("\n\nCOMPARISON WIH LAMINAR FLOW\n\n\n") -printf("\n\n Disturbance thicknes: %.3f \n\n",D) -printf("\n\nDisplacement thickness: %.3f\n\n",DS) -printf("\n\nWall shear stress: %.3f \n\n",WSS) diff --git a/746/CH9/EX9.05/9_05.sce b/746/CH9/EX9.05/9_05.sce deleted file mode 100755 index 8e146f23f..000000000 --- a/746/CH9/EX9.05/9_05.sce +++ /dev/null @@ -1,19 +0,0 @@ -//force and power// -pathname=get_absolute_file_path('9.05.sce') -filename=pathname+filesep()+'9.05-data.sci' -exec(filename) -//Speed in m/s: -U=s*6076*0.305/3600 -//Reynolds number: -Re=U*L/v -//Drag coefficient: -Cd=0.455/log10(Re)^2.58-1610/Re -//Area(in m^2): -A=L*(W+D) -//Drag force(in N) -Fd=Cd*A*0.5*d*U^2 -//Power required to overcome skin friction drag(in W): -P=Fd*U -printf("\n\nRESULTS\n\n") -printf("\n\nDrag force: %f N\n\n",Fd) -printf("\n\nPower required to overcome skin friction drag: %.3f W\n\n",P) diff --git a/746/CH9/EX9.06/9_06.sce b/746/CH9/EX9.06/9_06.sce deleted file mode 100755 index 3786c2baa..000000000 --- a/746/CH9/EX9.06/9_06.sce +++ /dev/null @@ -1,17 +0,0 @@ -//Bending moment// -pathname=get_absolute_file_path('9.06.sce') -filename=pathname+filesep()+'9.06-data.sci' -exec(filename) -//Velocity in m/sec: -V=s*5/18 -//Reynolds number: -Re=d*V*D/u -//Value of Cd is obtained as: -Cd=0.35; -//Area(in m^2): -A=L^2; -//Moment about the chimney base(in N-m): -M0=Cd*A*D/4*d*V^2 -printf("\n\nRESULTS\n\n") -printf("\n\nBending moment at the bottom of the chimney: %.3f N-m\n\n",M0) - diff --git a/746/CH9/EX9.07/9_07.sce b/746/CH9/EX9.07/9_07.sce deleted file mode 100755 index 3f1d56c77..000000000 --- a/746/CH9/EX9.07/9_07.sce +++ /dev/null @@ -1,8 +0,0 @@ -//Time required// -pathname=get_absolute_file_path('9.07.sce') -filename=pathname+filesep()+'9.07-data.sci' -exec(filename) -//Time required to decelerate to 100 mph(in seconds): -t=(s1-s2)*2*w/(s1*s2)/Cd/d/A/g*3600/5280 -printf("\n\nRESULTS\n\n") -printf("\n\nTime required to decelerate to 100 mph: %.3f seconds\n\n",t) diff --git a/746/CH9/EX9.08/9_08.sce b/746/CH9/EX9.08/9_08.sce deleted file mode 100755 index 54e4c9eca..000000000 --- a/746/CH9/EX9.08/9_08.sce +++ /dev/null @@ -1,28 +0,0 @@ -//Optimum cruise speed// -pathname=get_absolute_file_path('9.08.sce') -filename=pathname+filesep()+'9.08-data.sci' -exec(filename) -//Plotting velocity with drag force -V=175:25:455; - -[m n]=size(V); -for i=1:n - CL(i)=2*W/p*(3600/V(i)/5280)^2/A; - Cd(i)=Cd0+CL(i)^2/%pi/ar; - Fd(i)=Cd(i)/CL(i)*W; - FD(i)=Fd(i)/1000; -end -plot(V,FD) -xtitle('Flight speed vs thrust','Flight Speed(in mph)','Drag Force(in 1000lbf)') -//Optimum cuise speed at speed level is obtained to be 320 mph from the graph. -Vosl=320; -//Ratio of speeds at 30000 ft and at sea level is given by: -r=sqrt(1/0.375); -//Stall speed at 30000ft is(in mph): -Vs3=Vssl*r; -//Optimum Cruise speed at 30000ft(in mph): -Vo3=Vosl*r; -printf("\n\nRESULTS\n\n") -printf("\n\nOptimum cruise speed at sea level: %.3f mph\n\n",Vosl) -printf("\n\nStall speed at 30000 ft: %.3f mph\n\n",Vs3) -printf("\n\nOptimum cruise speed at 30000 ft: %.3f\n\n",Vo3) diff --git a/746/CH9/EX9.09/9_09.sce b/746/CH9/EX9.09/9_09.sce deleted file mode 100755 index cef9919e3..000000000 --- a/746/CH9/EX9.09/9_09.sce +++ /dev/null @@ -1,20 +0,0 @@ -//Aerodynamic and Radius// -pathname=get_absolute_file_path('9.09.sce') -filename=pathname+filesep()+'9.09-data.sci' -exec(filename) -//Reynolds number: -//Value of wD/2V: -W=0.5*N*D/1000/V*2*%pi/60 -Red=V*D/v; -//For this value, CL is obtained as: -CL=0.3; -//Aerodynamic lift(in N): -FL=%pi/8*CL*(D/1000)^2*d*V^2; -//Radius of curvature of the path in the vertical plane(in m) with topspin: -Rts=V^2/(g+FL/(m/1000)); -//Radius of curvature without topspin(in m): -Rwts=V^2/g; -printf("\n\nRESULTS\n\n") -printf("\n\nAerodynamic lift acting on the ball:%.3f N\n\n",FL) -printf("\n\nRadius of curvature of the path when ball has topspin:%.3f m\n\n",Rts) -printf("\n\nRadius of curvature of the path when ball has topspin: %.3f m\n\n",Rwts) -- cgit