program PROG4_CFD

! Condução de calor unidimensional permanente em 
! parede plana com condutividade térmica variável

! Versão original 1.0 ( 9 Abr 09)
! Versão atual    1.0 ( 9 Abr 09)
! Última alteração =    6 Out 09

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

! MODELO MATEMÁTICO (equação diferencial):

!   d(k.dT/dx)/dx = 0

!   T  = temperatura (incógnita e variável primária) (C)
!   x  = coordenada espacial (m)
!   k  = condutividade térmica (W/m.K)
!   Condições de contorno (C.C.) do tipo Dirichlet (T prescrita)
!   Variáveis de interesse:
!     Temperatura no meio do domínio: Tc (C)
!     Média da temperatura: Tm (C)
!     Taxa de transferência de calor em x=0: q0 (W)
!     Taxa de transferência de calor em x=L: qL (W)
!     Média da norma l1 do erro numérico de T: Em (C/nó)
!   Solução analítica disponível para o seguinte problema:
!     k = exp(T)
!     TA = 0
!     TB = 1
!     A  = 1
!     L  = 1

! MODELO NUMÉRICO

!   Método numérico: volumes finitos
!   Malha unidimensional uniforme
!   Solver: TDMA
!   Precisão: dupla
!   Linguagem Fortran 90
!   Aplicativo usado: Microsoft Developer Studio, Fortran PowerStation 4.0
!   Tipo de projeto: Console Application
!   N = número de volumes de controle reais
!   Tipos de aproximação numérica:
!     Termo difusivo (dT/dx): CDS-2
!     k nas faces: método da interpolação linear dos k nodais
!     Tm: integração pela regra do retângulo
!     q0 e qL: CDS-2 com volume fictício
!   Aplicação das condições de contorno com volumes fictícios (P=0 e P=N+1):
!   Sistema de equações: ap(P)*T(P) = aw(P)*T(P-1) + ae(P)*T(P+1) + bp(P)
!   Chute inicial de T numérico = média das C.C.
!   A solução numérica pode ser obtida para quaisquer valores de TA, TB, A, L e k;
!     os valores de TA, TB, A e L devem ser fornecidos via arquivo de dados
!     os valores de k devem ser fornecidos dentro do programa-fonte

! ARQUIVOS envolvidos no programa (extensão a prog4_cfd_):
!   fonte.f90 = programa-fonte
!   dados.txt = arquivo de dados do programa
!   T.dat     = dados para os gráficos
!   i.dat     = dados para o gráfico das iterações
!   T.gnu     = arquivo de comandos para o Wgnuplot gerar gráfico de T(x)
!   E.gnu     = arquivo de comandos para o Wgnuplot gerar gráfico de erro(x)
!   I.gnu     = arquivo de comandos para o Wgnuplot gerar gráfico de iterações

! Outros arquivos
!   "caso"            = saída principal do programa
!   leia-me_prog4_cfd = dados básicos sobre o programa
!   Notepad.exe       = editor de arquivos de dados e resultados
!   Wgnuplot.exe      = programa gerador de gráfico

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