summaryrefslogtreecommitdiff
path: root/896/CH13/EX13.4/4.sce
blob: 4c2bef7b335454ec70f2fea746d0e76cac3f4d26 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
clc
//Example 13.4
//Calculate the pressure gradient
D=0.152//m
V_avg=3.04//m/s
rho=1000//kg/m^3
n=0.41//dimentionless
K=0.66//dimentionless
R_pl=8*rho*V_avg^(2-n)*D^n/(K*(2*(3*n+1)/n)^n)//dimentionless
printf("The reynolds number is %f\n",R_pl);
f=0.004//dimentionless (fanning friction factor)
//Let DP denote the pressure gradient
DP=4*f*(rho/D)*(V_avg^2/2)/1000//KPa/m
printf("The pressure gradient is %f KPa/m",DP);