summaryrefslogtreecommitdiff
path: root/536/CH10/EX10.4/Example_10_4.sce
blob: 1a1224c237b7cc9eba1b52156964ef765356ae05 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
clear;
clc;

printf("Example 10.4\n");

Q=3e-6; //Flow rate of water
Meu=1e-3; //Viscosity of water
D=1.5e-9; //diffusivity of carbon dioxide in water
rho=1e3; //Density of water

//the mean velocity of flow is governed by equation 3.87 in which sin(phi)is 
//put equal to unity for a vertical surface:
s=(Q*1e2*Meu*3/(rho*9.81))^(1/3);//Thickness of film

//A=Ux/Us=0.95;
A=0.95;
y=s*(1-A)^0.5;//The distance below the free surface
//using equation 10.108 and using tables of error fuctions
t=(1.305/1.822)^2
Us=rho*9.81*s^2/(2*Meu);//surface velocity
L=Us*t;//The maximum lend=gth of column
printf("\n The maximum length of column = %.2f m",L);