%% Vaba välja- ja ülevooluga puhverpaagi mudel (vt. nt)
%
% Mahuühikutes:
% Q - püsivoog läbi paagi, qi - sissevoolu muutus, qo - väljavoolu muutus [m3/s]

H = 3 %püsinivoo (head) [m]
A = 0.46 % nivoo pinna pindala [m2]
Ao = 0.00185 % väljalaske toru ristlõike pindala [m2] algväärtus
g = 9.81; % raskuskiirendus [m/s2]
roo = 1000; % vedeliku tihedus [kg/m3]

%
% Massiühikutes \frac{\textrm{dM}}{\textrm{dt}}=\left(w_i -w_o \right), vool on võrdeline rõhkude vahega w_o =\frac{1}{R}{\left(p_1 -p_2 \right)}^{\frac{1}{\alpha }} , p_1 =\rho \textrm{gH}\left(t\right),\;p_2 \approx P_a
% ehk  \rho A\frac{\textrm{dH}\left(t\right)}{\textrm{dt}}=w_i -\rho A_o \sqrt{2\textrm{gH}\left(t\right)} [kg/s]

Wi = 14.3 % sissevool püsireziimis [kg/s] (=Wo)
Vo = sqrt(2*g*H) % väljavoolu kiirus püsireziimis [m/s]
Wo = roo*Ao*Vo % väljavool püsireziimis [kg/s]

%
% lineariseeritud mudelina C\frac{\textrm{dh}\left(t\right)}{\textrm{dt}}=-\frac{1}{R}h\left(t\right)+u\left(t\right)

R = sqrt(2*g*H)/(roo*g*Ao) % väljalaske hüdrauliline takistus [m.s/kg]
C = roo*A % mahtuvus [kg/m]
T = R*C % ajakonstant [s]