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

zr : double = 1.0e1; {elevação do nível do
reservatório (estático)}
z0 : double = 9.35e0; {elevação da tomada de água -
ponto "zero"}
zt : double = 7.56e0; {elevação da base da chaminé}
q0 : double = 9.64e-4; {vazão final}
h1 : double = 1.91e0; {perda de carga inicial}
g : double = 9.81e0;
{aceleração da gravidade}
ta : double = 1.2e2; {tempo de cálculo}
tf : double = 0.2e0; {tempo p/abertura total da
válvula}
kr : double = 1.0e0;
kc : double = 1.0e0;
procedure inic_dad;
var
aux : double;
i,j : integer;
begin
nx[1] := 1000;
{partição do tubo1}
dx[1] := lt[1]/nx[1];
nx[2] := trunc(lt[2]/dx[1]) + 1; {partição do tubo2}
dx[2] := lt[2]/nx[2];
dt := 0.8e0*dx[1]/a[1]; {definição do delta t}
ndt := trunc(0.01e0/dt); {Definição do delta t da valvula}
(*
if ndt < 200
then
ndt := 1
else
ndt := ndt div 100;
*)
writeln(dt,nx[1]:10,nx[2]:10,ndt:10); {plotagem dos
valores}
st[1] := (z0-zt)/lt[1];
st[2] := 0.0e0;
{seno do ângulo do tubo1}
{seno do ângulo do tubo2}
142