blob: 318136c9fdfc2d936de2cbc887b96f93a977187b (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
|
clear;
clc;
printf("\n Example 7.2");
//The slurry contains 100kg whiting/m^3 of water
printf("\n Volume of 100 kg whiting = %f m^3",100/3000);
printf("\n Volume of cake = %f m^3",0.0333/0.6);
printf("\n Volume of liquid in cake = %f m^3",0.05556*0.4);
printf("\n Volume of filtrate = %.3f m^3",(1-0.0222));
printf("\n volume of cake/volume of filtrate v = %f",0.0556/0.978);
A = 10^(-4); //area in sq meters
deltaP = -1.65*10^(5); //P is in pascals
l = 0.01; //length is in meters
vol_flow_rate = 2*10^(-8); //Volume flow rate is in m^3/sec
u = 10^(-3); //vicosity is in Ns/m^2
r = poly([0],'r');
r1 = roots((10^4)*(2*10^(-8)*r)-1.65*10^(5)/(10^(-5)));
printf("\n r = %.2f*10^(13)/m^2",r1*10^(-13));
function[Lopt]=optimum()
Lopt = 1.161*10^(-3)*(900)^(0.5); //t = 900 secs
funcprot(0);
endfunction
printf("\n optimum frame thickness = %.1f mm",2*optimum()*1000);
//total cycle time = 1.015L^2 + 900
//rate of cake production R = L/(1.015L^2 + 900)
L = poly([0],'L');
L1 = roots(1.025*10^(6)*L^2 + 900 - 2.050*10^(6)*L^2);
printf("\n Frame thickness = %.2f mm",2*L1(1)*10^3);
|
