From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 650/CH3/EX3.10/10.sce | 13 +++++++++++++ 650/CH3/EX3.12/12.sce | 18 ++++++++++++++++++ 650/CH3/EX3.17/17.sce | 10 ++++++++++ 650/CH3/EX3.18/18.sce | 17 +++++++++++++++++ 650/CH3/EX3.19/19.sce | 22 ++++++++++++++++++++++ 650/CH3/EX3.2/2.sce | 8 ++++++++ 650/CH3/EX3.20/20.sce | 13 +++++++++++++ 650/CH3/EX3.4/4.sce | 17 +++++++++++++++++ 650/CH3/EX3.5/5.sce | 15 +++++++++++++++ 650/CH3/EX3.6/6.sce | 7 +++++++ 650/CH3/EX3.8/8.sce | 20 ++++++++++++++++++++ 650/CH3/EX3.9/9.sce | 10 ++++++++++ 12 files changed, 170 insertions(+) create mode 100755 650/CH3/EX3.10/10.sce create mode 100755 650/CH3/EX3.12/12.sce create mode 100755 650/CH3/EX3.17/17.sce create mode 100755 650/CH3/EX3.18/18.sce create mode 100755 650/CH3/EX3.19/19.sce create mode 100755 650/CH3/EX3.2/2.sce create mode 100755 650/CH3/EX3.20/20.sce create mode 100755 650/CH3/EX3.4/4.sce create mode 100755 650/CH3/EX3.5/5.sce create mode 100755 650/CH3/EX3.6/6.sce create mode 100755 650/CH3/EX3.8/8.sce create mode 100755 650/CH3/EX3.9/9.sce (limited to '650/CH3') diff --git a/650/CH3/EX3.10/10.sce b/650/CH3/EX3.10/10.sce new file mode 100755 index 000000000..a097d0781 --- /dev/null +++ b/650/CH3/EX3.10/10.sce @@ -0,0 +1,13 @@ +clc +rho=1000; //kg/m^3 +u=0.1; //Ns/m^2 +g=9.81; //m/s^2 +L=10; //m +H=2; //m +Q=14/3600; //m^3/s +d=0.05; //m + +dp=rho*g*(L+H) - (128*Q*u*L/%pi/0.05^4); +disp("Pressure drop across the valve =") +disp(dp) +disp("N/m^2") \ No newline at end of file diff --git a/650/CH3/EX3.12/12.sce b/650/CH3/EX3.12/12.sce new file mode 100755 index 000000000..4ea5e0207 --- /dev/null +++ b/650/CH3/EX3.12/12.sce @@ -0,0 +1,18 @@ +clc +Q=3*10^(-6); // m^3/s +u=0.001; // Ns/m^2 +W=1; +rho=1000; // kg/m^3 +g=9.81; // m/s^2 +d=1.016*10^(-4); // m + +theta=asind(3*Q*u/W/rho/g/d^3); +disp("Exact angle of inclination =") +disp(theta) + +d1=1.25*10^(-4); // m + +u1=W*rho*g*sind(theta)*(d1^3)/(3*Q); +disp("Viscosity of the second liquid =") +disp(u1) +disp("Ns/m^2") \ No newline at end of file diff --git a/650/CH3/EX3.17/17.sce b/650/CH3/EX3.17/17.sce new file mode 100755 index 000000000..a1e907537 --- /dev/null +++ b/650/CH3/EX3.17/17.sce @@ -0,0 +1,10 @@ +clc +u=1.5; // Ns/m^2 +v=0.5; // m/s +H=0.02/2; // m + +t=-u*3*v/H; +disp("The shear stress =") +disp(t) +disp("N/m^2") +disp("It acts in the opposite direction to the flow.") \ No newline at end of file diff --git a/650/CH3/EX3.18/18.sce b/650/CH3/EX3.18/18.sce new file mode 100755 index 000000000..4f25db397 --- /dev/null +++ b/650/CH3/EX3.18/18.sce @@ -0,0 +1,17 @@ +clc +N=600/60; // revolutions per sec +r=0.025; // m +t=400; // N/m^2 +l=0.002; // m + +w=2*%pi*N; + +u=t*l/w/r; +disp("Viscosity =") +disp(u) +disp("Ns/m^2") + +T=integrate('2*%pi*u*w/l*r^3', 'r', 0, r); +disp("Torque =") +disp(T) +disp("Nm") \ No newline at end of file diff --git a/650/CH3/EX3.19/19.sce b/650/CH3/EX3.19/19.sce new file mode 100755 index 000000000..d69cb1ea1 --- /dev/null +++ b/650/CH3/EX3.19/19.sce @@ -0,0 +1,22 @@ +clc; +u=0.153; //Ns/m^2 +r=0.05; // m +N=30; // rps +t=2/10^5; //s +L=0.2; // m + +tau=u*(2*%pi*N*r/t); + +F=tau*2*%pi*r*L; + +T=F*r; + +w=2*%pi*N; +P=T*w; + +disp("The torque on the bearing is found to be "); +disp(T); +disp("Nm"); +disp("and the power required to overcome the frictional resistance is "); +disp(P); +disp("W"); \ No newline at end of file diff --git a/650/CH3/EX3.2/2.sce b/650/CH3/EX3.2/2.sce new file mode 100755 index 000000000..03b3ccd87 --- /dev/null +++ b/650/CH3/EX3.2/2.sce @@ -0,0 +1,8 @@ +clc +Re=2000; +d=0.008; //m + +L1=0.058*Re*d; +disp("The furthest distance the fluid can flow into the 8 mm inside diameter pipe before fully developed laminar flow can exist is "); +disp(L1); +disp("m"); \ No newline at end of file diff --git a/650/CH3/EX3.20/20.sce b/650/CH3/EX3.20/20.sce new file mode 100755 index 000000000..87849c248 --- /dev/null +++ b/650/CH3/EX3.20/20.sce @@ -0,0 +1,13 @@ +clc; +t=0.0005; // s +P=22; // +r=300/60; // +R_1=0.1; // +R_2=0.0625; // + +w=2*%pi*r; + +u=2*t*P/(%pi*w^2*((R_1)^4-(R_2)^4)); +disp("The viscosity of the oil is found to be "); +disp(u); +disp("Nsm-2."); diff --git a/650/CH3/EX3.4/4.sce b/650/CH3/EX3.4/4.sce new file mode 100755 index 000000000..3518167d5 --- /dev/null +++ b/650/CH3/EX3.4/4.sce @@ -0,0 +1,17 @@ +clc +del_p=90*10^3; // N/m^2 +d=0.126; // m +R=0.126/2; // m +u=1.2; +L=60; // m +Rho=1260; + +Q=%pi * del_p * R^4 / (8*u*L); +disp("The glycerol delivery rate is "); +disp(Q); +disp("m^3/s"); + +Re=4*Rho*Q/(u*%pi*d); +disp("The Reynolds number is "); +disp(Re); +disp("As Re is below 2000, therefore confirming laminar flow."); diff --git a/650/CH3/EX3.5/5.sce b/650/CH3/EX3.5/5.sce new file mode 100755 index 000000000..4ea347f10 --- /dev/null +++ b/650/CH3/EX3.5/5.sce @@ -0,0 +1,15 @@ +clc +u=0.015; //Ns/m^2 +Q=0.004/60; //m^3/s +dp=100; +rho=1100; //kg/m^3 + +R=(8*u*Q/(%pi*dp))^(1/4); +Re=(4*rho*Q/(%pi*u*(2*R))); + +disp("Diameter of the pipe =") +disp(R) +disp("m") + +disp("Reynolds number =") +disp(Re) \ No newline at end of file diff --git a/650/CH3/EX3.6/6.sce b/650/CH3/EX3.6/6.sce new file mode 100755 index 000000000..a1a02600b --- /dev/null +++ b/650/CH3/EX3.6/6.sce @@ -0,0 +1,7 @@ +clc +u=0.03; //Ns/m^2 +Q=10^(-7); //m^3/s +dp=integrate('8*u*Q/%pi/0.005^4/(1-L)^4', 'L', 0, 0.5) +disp("Pressure difference =") +disp(dp) +disp("N/m^2") \ No newline at end of file diff --git a/650/CH3/EX3.8/8.sce b/650/CH3/EX3.8/8.sce new file mode 100755 index 000000000..18a04a849 --- /dev/null +++ b/650/CH3/EX3.8/8.sce @@ -0,0 +1,20 @@ +clc +u=0.1; // Ns/m^2 +d=0.1; //m +R=0.05; // m +Rho=900; //kg/m^3 + +v_max=2; // m/s +v=v_max/2; // m/s + +disp("At the pipe wall (r =R), therefore, the shear stress is"); +Tw=-2*u*v_max/R; +disp(Tw); +disp("N/m^2"); +disp("The negative sign indicates that the shear stress is in the opposite direction to flow."); + +disp("pressure drop per metre length of pipe is"); + +del_p=4*u*v_max/R^2; +disp(del_p); +disp("N/m"); \ No newline at end of file diff --git a/650/CH3/EX3.9/9.sce b/650/CH3/EX3.9/9.sce new file mode 100755 index 000000000..2d42012e8 --- /dev/null +++ b/650/CH3/EX3.9/9.sce @@ -0,0 +1,10 @@ +clc +u=0.032; // Ns/m^2 +Re=2000; // maximum value +Rho=854; +del_p=150; // N/m^2 + +d=(32*u^2*Re/(Rho*del_p))^(1/3); +disp("The maximum inside diameter is found to be ") +disp(d) +disp("m") \ No newline at end of file -- cgit