model circuitoHidraulico_0 constant Real PI = 2*Modelica.Math.asin(1.0); constant Real g = 10.0; parameter Real rho = 1000.0; parameter Real p0 = 0.1e6; parameter Real A1 = 1.0; parameter Real A2 = 1.0; parameter Real H1 = 10.0; parameter Real H2 = 8.0; parameter Real R = 25e6; Real h1, h2; Real q; Real p1, p2; initial equation h1 = H1; h2 = H2; equation A1 * der(h1) = -q; A2 * der(h2) = q; p1 = p0 + rho * g * h1; p2 = p0 + rho * g * h2; p1 - p2 = R * q^2; annotation( Documentation(info = "Modelo com dois reservatórios hidráulicos e resistência de perda de carga constante.

Com incompabilidade no cálculo da vazão quando a perda de carga na tubulação fica numericamente negativa.
")); end circuitoHidraulico_0;