The Three Reservoir Problem - Maple Help
For the best experience, we recommend viewing online help using Google Chrome or Mozilla Firefox.

Online Help

All Products    Maple    MapleSim


The Three Reservoir Problem

Introduction

Three reservoirs at different elevations are connected through a piping network at a single point, with an outflow from the common junction.  This application will calculate the flow rates, flow directions, and head at the common junction.

 

The Bernoulli equation (ignoring the losses associated with pipe fittings) for Reservoir 1 is

 

 z1=f1 L1D1V122 g+H

 

where z1 is the reservoir elevation, f1 is the friction factor, L1 is the pipe length, D1 is the pipe diameter, V1 is the liquid velocity, and H is the head at the common junction.  However, since flow can be either into or out of a reservoir (that is, a positive or negative velocity), we rewrite the Bernoulli equation thus

 

z1=f L1D1V1 V12 g+H

 

Similar equations can be defined for the other reservoirs.  These equations, along with the continuity equation,

 

Q1+Q2+Q3=Qoutflow

 

and a correlation for the friction factor, can be used to find the flow rates, flow directions, and the head at the junction.

 

> 

restart:withThermophysicalData:

Parameters

Liquid viscosity and density

> 

μ≔Propertyviscosity,water,pressure=1atm,temperature=298Kρ≔Propertydensity,water,pressure=1atm,temperature=298K

893.07×10-6⁢Pa⁢s

997.09⁢kgm3

(2.1)

Reservoir elevations

> 

z__1≔90 m:z__2≔85 m:z__3≔60 m:

 

Lengths, diameters, and cross-sectional areas of the pipes connecting the reservoirs to the junction

> 

L__1≔2000 m:L__2≔1500 m:L__3≔3000 m:

> 

D__1≔0.3 m:D__2≔0.25 m:D__3≔0.25 m:

> 

A__1≔π D__12/4: A__2≔π D__22/4: A__3≔π D__32/4:

Roughness of the pipe

> 

e≔0.0005⁢m:

 

Outflow from the junction

> 

Qoutflow≔0.03⁢m3 s−1:

 

Gravitational constant

> 

g≔9.81⁢m s−2:

Bernoulli Equations for the Three Reservoir System

Bernoulli analysis of the system

> 

sys≔ z__1&equals;fric__1⁢L__1 Q__1 Q__1&sol;2 g⁢D__1 A__12&plus;H&comma; z__2&equals;fric__2⁢L__2 Q__2 Q__2&sol;2 g⁢D__2 A__22&plus;H&comma; z__3&equals;fric__3 L__3 Q__3 Q__3&sol;2 g⁢D__3 A__32&plus;H&comma; Q__1&plus;Q__2&plus;Q__3&equals;Qoutflow&comma; Rey__1 &equals; 4 absQ__1 &rho;&sol;π D__1 &mu;&comma; Rey__2 &equals; 4 absQ__2 &rho;&sol;π D__2 &mu;&comma; Rey__3 &equals; 4 absQ__3 &rho;&sol;π D__3 &mu;&comma; fric__1 &equals; piecewiseRey__1<2500&comma; 64&sol;Rey__1&comma; 1&sol;1.8 log10e&sol;3.7 D__11.11&plus;6.9&sol;Rey__12&comma; fric__2 &equals; piecewiseRey__2<2500&comma; 64&sol;Rey__2&comma; 1&sol;1.8 log10e&sol;3.7 D__21.11&plus;6.9&sol;Rey__22&comma; fric__3 &equals; piecewiseRey__3<2500&comma; 64&sol;Rey__3&comma; 1&sol;1.8 log10e&sol;3.7 D__31.11&plus;6.9&sol;Rey__32&colon;

> 

estimates≔H&equals;1 m&comma;Q__1=0.1 m3 s−1&comma;Q__2=0.1 m3 s−1&comma;Q__3=0.1 m3 s−1&comma;Rey__1=1&comma;Rey__2&equals;1&comma;Rey__3&equals;1&comma;fric__1=1&comma;fric__2&equals;1&comma;fric__3&equals;1&colon;

Calculation of Flow Rates, Flow Directions, and Junction Head

Positive flow rates indicate flow out of a reservoir, while negative flow rates indicate flow into a reservoir

> 

Digits≔20&colon;

> 

results≔fsolvesys&comma;estimates

H=8.31×101⁢m&comma;Q__1=6.67×10-2⁢m3s&comma;Q__2=2.49×10-2⁢m3s&comma;Q__3=−6.16×10-2⁢m3s&comma;Rey__1=3.16×105&comma;Rey__2=1.42×105&comma;Rey__3=3.50×105&comma;fric__1=2.29×10-2&comma;fric__2=2.46×10-2&comma;fric__3=2.39×10-2

(4.1)
>