avaliações numéricas de chaminés de equilÃbrio - ppgerha ...
avaliações numéricas de chaminés de equilÃbrio - ppgerha ...
avaliações numéricas de chaminés de equilÃbrio - ppgerha ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
zr : double = 1.0e1; {elevação do nível do<br />
reservatório (estático)}<br />
z0 : double = 9.35e0; {elevação da tomada <strong>de</strong> água -<br />
ponto "zero"}<br />
zt : double = 7.56e0; {elevação da base da chaminé}<br />
q0 : double = 9.64e-4; {vazão inicial}<br />
h1 : double = 1.91e0; {perda <strong>de</strong> carga inicial}<br />
g : double = 9.81e0; {aceleração da gravida<strong>de</strong>}<br />
ta : double = 1.2e2; {tempo <strong>de</strong> cálculo}<br />
tf : double = 0.2e0; {tempo <strong>de</strong> fechamento da válvula}<br />
kr : double = 1.0e0;<br />
kc : double = 1.0e0;<br />
procedure inic_dad;<br />
var<br />
aux : double;<br />
i,j : integer;<br />
begin<br />
nx[1] := 1000;<br />
{partição do tubo1}<br />
dx[1] := lt[1]/nx[1];<br />
nx[2] := trunc(lt[2]/dx[1]) + 1; {partição do tubo2}<br />
dx[2] := lt[2]/nx[2]; dt := 0.8e0*dx[1]/a[1];<br />
{<strong>de</strong>finição do <strong>de</strong>lta t}<br />
ndt := trunc(0.01e0/dt); {Definição do <strong>de</strong>lta t da valvula}<br />
(*<br />
if ndt < 200<br />
then<br />
ndt := 1<br />
else<br />
ndt := ndt div 100;<br />
*)<br />
writeln(dt,nx[1]:10,nx[2]:10,ndt:10); {plotagem dos<br />
valores}<br />
st[1] := (z0-zt)/lt[1];<br />
st[2] := 0.0e0;<br />
{seno do ângulo do tubo1}<br />
{seno do ângulo do tubo2}<br />
for i := 1 to 2 do<br />
begin<br />
a1[i] := ft[i]/(8.0e0*rt[i]);<br />
{celerida<strong>de</strong> do tubo}<br />
136