Simulação Numérica de Escoamento Reativo Quasi-Unidimensional em Motor-Foguete com Refrigeração Regenerativa Carlos Henrique Marchi (UFPR, gerente do projeto) Luciano Kiyoshi Araki (UFPR) Márcio Augusto Villela Pinto (UEPG) Cosmo Damião Santiago (UNIBRASIL) Fábio Alencar Schneider (UNICENP)
18
Embed
Simulação Numérica de Escoamento Reativo Quasi- Unidimensional em Motor- Foguete com Refrigeração Regenerativa Carlos Henrique Marchi (UFPR, gerente do.
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Simulação Numérica de Escoamento Reativo Quasi-
Unidimensional em Motor-Foguete com Refrigeração Regenerativa
Carlos Henrique Marchi (UFPR, gerente do projeto)
Luciano Kiyoshi Araki (UFPR)
Márcio Augusto Villela Pinto (UEPG)
Cosmo Damião Santiago (UNIBRASIL)
Fábio Alencar Schneider (UNICENP)
Problema
Divisão do problema em três subproblemas:• Câmara-Tubeira: Escoamento reativo, turbulento
de gases num motor-foguete.• Paredes: Condução de calor através das mesmas
entre os gases de combustão e o fluido refrigerante.
• Canais: Escoamento turbulento do fluido refrigerante nos canais ao redor da tubeira.
Figura 1: Motor-foguete bipropelente com refrigeração regenerativa.
Fonte: Marchi et al., Solução Numérica de Escoamentos em Motor-Foguete com Refrigeração Regenerativa, XXI Cilamce.
Figura 2: Detalhes nos canais de refrigeração.
Fonte: Marchi et al., Solução Numérica de Escoamentos em Motor-Foguete com Refrigeração Regenerativa, XXI Cilamce.
Fases
• Fase 1: Escoamento unidimensional reativo, sem transferência de calor.
• Fase 2: Escoamento unidimensional reativo, com escoamento regenerativo.
Metodologia
• Método dos Volumes Finitos.• Funções de interpolação de segunda ordem, co-
localizado.• Formulação apropriada a qualquer regime de
velocidades.• Malha uniforme.• Estimativa de erro GCI.
Modelos físicos
Monogás, com propriedades constantes; Monogás, com propriedades variáveis; Escoamento congelado; Escoamento em equilíbrio; Escoamento com taxa finita de reação.
Modelos químicos
Modelo Número de reações
Número de espécies
Espécies envolvidas
0 0 3 H2O, O2, H2
1 1 3 H2O, O2, H2
2 2 4 H2O, O2, H2, OH
3 4 6 H2O, O2, H2, OH, O, H
4 4 6 H2O, O2, H2, OH, O, H
5 8 6 H2O, O2, H2, OH, O, H
7 8 6 H2O, O2, H2, OH, O, H
10 6 8 H2O, O2, H2, OH, O, H, HO2, H2O2
9 18 8 H2O, O2, H2, OH, O, H, HO2, H2O2
Resultados – Fase 1
Estimativas de erros para temperatura (mod. monogás com propriedades constantes)
Analítico CEA (cong.) Mod. 3 (cong.) CEA (eq.) Mod. 3 (eq.) Mod. 31 (t. fin.) Mod. 32 (t. fin.) Mod. 5 (t. fin.) Mod. 7 (t. fin.) Mod. 10 (t. fin.) Perfil
Posição [m]
Tem
pera
tura
[K
]
Distribuição de temperaturas ao longo da tubeira
Resultados – Fase 1
0.0 0.1 0.2 0.3 0.4 0.5
0.75
0.80
0.85
0.90
0.95
CEA (cong.) Mod. 3 (cong.) CEA (eq.) Mod. 3 (eq.) Mod. 31 (t. fin.) Mod. 32 (t. fin.) Mod. 5 (t. fin.) Mod. 7 (t. fin.) Mod. 10 (t. fin.) Perfil
Posição [m]
Fra
ção
más
sica
Distribuição de frações mássicas de H2O ao longo da tubeira
Resultados – Fase 1
Tempo de CPU
Monogás, isentrópico Esc. com taxa finita
Malha Propriedades
constantes
Propriedades
variáveis
Esc.
congelado
(mod. 3)
Esc. em equilíbrio
(mod. 3) mod. 31 Mod. 32
80 volumes 3,08 s 0,750 s 0,969 s 2,34 min 17,0 min 27,3 min
2560 volumes 2,07 min 3,02 min 2,66 min 1,98 dia 1,66 h 2,61 h
10240 volumes 56,9 min 1,65 h 1,49 h --- --- ---
Resultados – Fase 2
Distribuição de 80 volumes
Resultados – Fase 2
Distribuição de temperaturas ao longo da tubeira
Resultados – Fase 2
Distribuição de temperaturas ao longo da parede da tubeira