avaliações numéricas de chaminés de equilíbrio - ppgerha ...

a2[i] := a[i]*a[i]/g;
celeridade}
a3[i] := 0.5e0/dx[i];
end;
{coeficiente que inclui a
{coeficiente denominador}
b1 := 0.5e0*kr/g; {coeficiente perda de carga da tomada}
b2 := 0.5e0*dt/ac; {coeficiente de perda de carga de
saída}
b3 := 0.5e0*kc/(ac*ac*g); {coeficiente fator para o
cálculo}
for i := 1 to 2 do
begin
aux := q0/at[i];
for j := 0 to nx[i] do
v[i,2]^[j] := aux;
end;
{velocidade inicial}
h[1,2]^[0] := zr - b1*aux*aux;; {perda de carga no tubo1}
aux := h1/nx[1];
{distribuição da perda de
carga no tubo1}
for j := 1 to nx[1] do
{perda de carga em cada
ponto}
h[1,2]^[j] := zr - j*aux;
hch[2] := h[1,2]^[nx[1]];
qch[2] := 0.0e0;
qvs[2] := q0;
{nível inicial na chaminé}
{vazão inicial na chaminé}
{vazão inicial da valvula}
h[2,2]^[0] := hch[2];
{distribuição da perda de
carga no tubo2}
aux := (h1/lt[1])*lt[2]/nx[2];
for j := 1 to nx[2] do
h[2,2]^[j] := hch[2] - j*aux;
end;
procedure lax;
var
vm,hm,
dvx,dhx : double;
i,j : integer;
137