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

____________________________________________________
const
nst = 10;
ndx = 2000;
type
vet1 = array[1..nst] of double;
vet2 = array[0..ndx] of double;
enpc = object
l,d,a,f,ac,
alfa
: vet1;
n,s
: integer;
dx
: double;
dh,dq,h,q,
ep,fp,lp,mp,np : ^vet2;
hp,bp,cp,dp,gp : array[1..2] of ^vet2;
ck
: array[1..12] of vet1;
constructor alocmem;
destructor librmem;
procedure ctes;
procedure param;
procedure dupvarr;
end;
___________________________________________________
FIGURA 5.12 - Estrutura do objeto ENPC (KAVISKI e FRANCO,1997).
No objeto ENPC estão incluídos os parâmetros identificadores de cada ramo
da rede em vetores com os elementos, representando os trechos: comprimento (l),
diâmetro (d), área do conduto (a), resistência ao escoamento (f) celeridade (ac) e o
ângulo do eixo do conduto (alfa).
Estão incluídos em vetores com n elementos os coeficientes necessários
para a aplicação do método da dupla-varredura (ep, fp, lp, mp, np, hp, bp, cp, dp e
gp) e as variáveis dependentes (h, q, dh e dq). Os 12 coeficientes ck representam
algumas constantes auxiliares, para cada trecho do ramo do conduto.
No objeto ENPC estão incluídos um método construtor e um destrutor para
alocação e deslocação das variáveis dinâmicas. Estão incluídos 3 procedimentos:
ctes, para determinação dos coeficientes ck; param, para inicialização dos
coeficientes usados para aplicação do método de diferenças finitas; e dupvarr, que
soluciona o sistema de equações pelo método da dupla-varredura.
Nesta dissertação empregou-se o objeto ENPC para solucionar o estudo de
caso do protótipo, onde considerou-se a chaminé de equilíbrio como um contorno
103