WOLFRAM|DEMONSTRATIONS PROJECT

Dynamic Simulation of a Gravity-Flow Tank

​
intial flow rate/design flow rate
0.6
The dynamic behavior of a tank and pipe system is described by the following coupled ordinary differential equations:
dv
dt
=
g
L
h-
K
F
g
c
ρ
A
p
2
v
,
A
T
dh
dt
=
F
0
-F
,
where
L
is the pipe length,
g
is the acceleration of gravity,
F
0
is the inlet flow rate to the tank,
F
is the outlet flow rate from the tank,
h
is the tank's height,
ρ
is the fluid's density,
v
is the fluid's velocity in the pipe,
g
c
=32.1740
is a unit conversion factor,
A
p
and
A
T
are the cross-sectional area of the pipe and the tank, respectively, and
K
F
is a proportionality coefficient that appears in the frictional force
K
F
L
2
v
.
The first equation represents a force balance on the outlet line and the second equation is the continuity equation on the liquid in the tank. We have the following relation:
F=
A
p
v
. If one neglects the dynamic behavior of the pipe, then we get the classic equation:
F∝
h
.
This Demonstration shows that it is important to take into account the dynamic behavior in the pipe (see blue curve) since in some cases it allows the prediction of tank overflow, a phenomenon that will not be observed otherwise (see red curve). All numerical values used in the present Demonstration can be found in Luyben's book.