program PROG2_CFD1

!     Escoamento laminar 1Dp, plenamente desenvolvido, num duto circular

!     Versão original 1.0 (12 Jun 03)

!     Versão atual    1.0 ( 6 Jul 06)

!     última alteração =    6 Jul 06

!     autor: Carlos Henrique Marchi (Curitiba, DEMEC/UFPR)

!     MODELO MATEMÁTICO (resumo)
!       Equação diferencial: (mi/r)*d(r*du/dr)/dr = C
!       Condição de contorno de Dirichlet em r = Rm: u(Rm) = 0
!       Condição de contorno de Neumann em r = 0: (du/dr)(0) = 0
!       r  = coordenada espacial (variável independente)
!       mi = viscosidade
!       u  = velocidade
!       C  = constante
!       Variáveis de interesse:
!          u(r)
!          umed = u média
!          umax = u máxima
!          F    = força do fluido sobre o duto
!       Solução analítica conhecida

!     MODELO NUMÉRICO (resumo)
!       Incógnita (variável primária, dependente): u
!       Método numérico: volumes finitos
!       Função de interpolação: CDS (variável primária u)
!       Propriedade constante: mi
!       Condição de contorno de Dirichlet em r = Rm
!       Condição de contorno de Neumann em r = 0
!       As condições de contorno são incorporadas aos volumes
!          com face coincidente com os contornos
!       umed: regra do retângulo
!       umax: ajuste de função quadrática
!       F: UDS
!       Malha uniforme
!       Solver: TDMA
!       Precisão: dupla
!       Linguagem FORTRAN 95
!       Aplicativo usado: Compaq Visual Fortran 6.6
!       Tipo de projeto: QuickWin Application
!       Expressão genérica do sistema de equações discretizado:
!       aP(i)*u(i) = aw(i)*u(i-1) + ae(i)*u(i+1) + bP(i) 
!            onde i = 1, 2, ... N

!     ARQUIVOS envolvidos no programa:
!       prog2_cfd1.f90 = programa principal
!       coef.f90       = calcula coeficientes e fontes do sistema linear
!       dados.f90      = lê e lista os dados do programa
!       result.f90     = resolve equações e gera listagens dos resultados
!       solvers.f90    = solver TDMA
!       varia.f90      = define todas as variáveis globais do programa
!       prog2_cfd1.ent = arquivo de dados do programa
!       *.txt          = listagem dos resultados
!       u.dat          = arquivo de dados para fazer gráfico
!       u.gnu          = arquivo de comandos para gerar gráfico
!       notepad.exe    = editor dos arquivos
!       Wgnuplot.exe   = visualizador do gráfico

! -----------------------------------------------

Seqüência para compilação dos módulos pela primeira vez após criar o projeto:

1) VARIA.f90

2) SOLVERS.f90

3) DADOS.f90

4) COEF.f90

5) RESULT.f90

6) PROG2_CFD1.f90