summaryrefslogtreecommitdiff
path: root/629/CH6/EX6.3/example6_3.sce
blob: 1d9d2d66312e61eaad7899349ecc05cd303a0f37 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
clear
clc
//Example 6.3 FORCE ON A NOZZLE
d1=0.06; //[m]
A1=%pi*d1^2/4 //area[m^2]
d2=0.01; //[m]
rho=1.22; //density[kg/m^3]
p1=3.7*1000; //[Pa]
//Bernoulli equation, p1+rho*v1^2/2=rho*v2^2/2
v2=sqrt(2*p1/(rho*(1-(d2/d1)^4))) //Exit velocity[m/s]
printf("\n The air speed at the exit of the nozzle = %.1f m/s.\n\n",v2)
v1=v2*(d2/d1)^2 //Inlet velocity[m/s]
m=rho*A1*v1 //mass flow rate[kg/s]
//Force on flange
F=m*(v2-v1)-p1*A1 //[N]
printf("\n The force on the flange = %.2f N.\n",-F)