program Mach2D_5p6 ! Escoamento bidimensional permanente de fluido compressível viscoso laminar ! Versão atual 5.6 (17 Fev 09) ! Subversão atual 5.6.01 (17 Fev 09) ! Versão original 1.0 (Mai 89) ! última alteração = 15 Mar 09 ! Autor versão atual: Carlos Henrique Marchi ! DEMEC/UFPR, Curitiba, PR ----------------------------------------- MODELO FÍSICO ----------------------------------------- Fluido: Compressível Contínuo Invíscido ou viscoso Não-reativo Gás perfeito termicamente Sem radiação térmica Propriedade constante: constante do gás perfeito Propriedades variáveis: calor específico, razão entre os calores específicos, viscosidade absoluta e condutividade térmica Com ou sem troca de calor com as paredes Com ou sem dissipação viscosa Escoamento: Estado permanente Bidimensional Geometria plana ou axissimétrica Não-reativo Ideal ou laminar Tubeira: Condição inicial para tubeira circular Paredes adiabáticas ou não Sem radiação térmica Propriedade variável: temperatura da parede ----------------------------------------- MODELO MATEMÁTICO ----------------------------------------- 1) Equação de conservação da massa 2) Equação de conservação da quantidade de movimento linear em X 3) Equação de conservação da quantidade de movimento linear em Y 4) Equação de conservação da energia térmica 5) Equação de estado dos gases perfeitos Incógnitas: u = Componente do vetor velocidade na direção X v = Componente do vetor velocidade na direção Y p = Pressão estática T = Temperatura ro = Massa específica Sistema de coordenadas não-ortogonais ksi-eta -------------------------------------------------- MODELO NUMÉRICO -------------------------------------------------- Método numérico: volumes finitos Condições de contorno aplicadas com volumes fictícios Malha estruturada não-ortogonal de faces centradas Funções de interpolação: CDS para os termos difusivos CDS com correção adiada sobre o UDS para os termos advectivos Arranjo co-localizado de variáveis (NHT/1994) Formulação totalmente implícita no tempo para processo iterativo Solver MSI para 5 e 9 diagonais não-nulas Método de acoplamento pressão-velocidade: SIMPLEC Solução segregada (seqüencial) das equações --------------------------------------------------- MODELO COMPUTACIONAL --------------------------------------------------- Precisão: dupla Linguagem FORTRAN 95 Aplicativo usado: Compaq Visual Fortran 6.6 Tipo de projeto: Console Application ARQUIVOS envolvidos no programa: Executável: Mach2D_5p6p01.exe Dados: Mach2D_5p6.in Mach2D_5p6.ent Mach2D_5p6p01.txt Mach2D_5p6_dominio.gnu Mach2D_5p6_malha.gnu Mach2D_5p6_residuo.gnu Mach2D_5p6_velocidade_u.gnu Mach2D_5p6_velocidade_v.gnu Mach2D_5p6_temperatura.gnu Mach2D_5p6_pressao.gnu Mach2D_5p6_massa_especifica.gnu Mach2D_5p6_Mach.gnu Mach2D_5p6_empuxo.gnu Saídas: "caso".txt = saída principal do programa "caso".itn = variáveis de interesse ao longo das iterações "caso".bin = dados para reinício "caso".Richardson_3p0 = resultado das variáveis de interesse para o programa Richardson 3.0 "caso".contorno_sul.dat "caso".contorno_oeste.dat "caso".contorno_norte.dat "caso".contorno_leste.dat "caso".malha.dat "caso".residuo.gnu "caso".velocidade_u.dat "caso".velocidade_v.dat "caso".temperatura.dat "caso".pressao.dat "caso".massa_especifica.dat "caso".Mach.dat "caso".empuxo.dat Fontes: Mach2D_5p6_solvers.f90 Mach2D_5p6_dados.f90 Mach2D_5p6_malha.f90 Mach2D_5p6_coeficientes.f90 Mach2D_5p6_usuario.f90 Mach2D_5p6_resultados.f90 Mach2D_5p6_main.f90 Outros arquivos: Notepad.exe = editor de arquivos Wgnuplot.exe = programa gerador de gráfico (versão 3.5) ! ----------------------------------------------- Seqüência para compilação dos módulos pela primeira vez após criar o projeto: a mesma da lista de arquivos-fontes acima.