70 Tabela 3.10 – Or<strong>de</strong>ns verda<strong>de</strong>iras e assintótica das aproximações das variáveis Variável Or<strong>de</strong>ns Or<strong>de</strong>m Verda<strong>de</strong>iras Assintótica T c * * T m 2, 4, 6... 2 I 2, 3, 4... 2 L * * Importante notar que a variável T c t<strong>em</strong> seu valor dado pela solução no centro do volume <strong>de</strong> controle, ou seja, no valor nodal. Isto implica dizer, que não há relação que agregue maior or<strong>de</strong>m ao erro <strong>de</strong> truncamento, a não ser os inseridos pela aplicação das próprias FI’s. Para a variável L o mesmo se aplica, sendo suas or<strong>de</strong>ns <strong>de</strong> erros <strong>de</strong>finidas pela variável T c , conforme será <strong>de</strong>monstrado adiante. 3.4.3 Análise a Priori da Or<strong>de</strong>m do Erro <strong>de</strong> Discretização Final Será visto agora, qual a or<strong>de</strong>m do erro final <strong>de</strong> cada variável <strong>de</strong> interesse, consi<strong>de</strong>rando que sua obtenção está contaminada com erro <strong>de</strong> truncamento das funções <strong>de</strong> interpolação aplicadas (subit<strong>em</strong> 3.4.1), assim como da relação para obtenção da própria variável (subit<strong>em</strong> 3.4.2). As or<strong>de</strong>ns dos erros <strong>de</strong> truncamento final, conseguidas após consi<strong>de</strong>rar o erro <strong>de</strong> truncamento das FI`s e os erros inseridos pelas aproximações dos mo<strong>de</strong>los numéricos das próprias variáveis <strong>de</strong> interesse, são mostradas na Tab.3.11. L<strong>em</strong>brando que as aproximações para o termo difusivo foram todas feitas com o CDS-2 (on<strong>de</strong> necessário), exceto para o QUICK, que foi aproximado com CDS-2 e CDS-4. Para facilitar o entendimento, daqui <strong>em</strong> diante as referências passam a ser <strong>de</strong> QUICK2 quando o termo difusivo foi aproximado com CDS-2, e QUICK4 quando aproximado com CDS-4 Novamente, não preocupou-se no <strong>de</strong>correr do texto, <strong>de</strong> <strong>de</strong>monstrar como chega-se a este valor, mas caso o leitor tenha interesse, estas <strong>de</strong>monstrações faz<strong>em</strong> parte do apêndice B.
71 Tabela 3.11 – Or<strong>de</strong>ns verda<strong>de</strong>iras e assintótica das variáveis Função <strong>de</strong> Interpolação Termo Advectivo Variáveis <strong>de</strong> Interesse T c T m I L Or<strong>de</strong>ns Verda<strong>de</strong>iras Or<strong>de</strong>m Assintótica Or<strong>de</strong>ns Verda<strong>de</strong>iras Or<strong>de</strong>m Assintótica Or<strong>de</strong>ns Verda<strong>de</strong>iras Or<strong>de</strong>m Assintótica Or<strong>de</strong>ns Verda<strong>de</strong>iras Or<strong>de</strong>m Assintótica UDS 1, 2, 3... 1 1, 2, 3... 1 1, 2, 3... 1 1, 2, 3... 1 CDS-2 2, 4, 6... 2 2, 4, 6... 2 2, 3, 4... 2 2, 4, 6... 2 UDS-2 2, 3, 4... 2 2, 3, 4... 2 2, 3, 4... 2 2, 3, 4... 2 WUDS 2, 4, 6... 2 2, 4, 6... 2 2, 3, 4... 2 2, 4, 6... 2 PLDS 2, 4, 6... 2 2, 4, 6... 2 2, 3, 4... 2 2, 4, 6... 2 ADS 2, 4, 6... 2 2, 4, 6... 2 2, 3, 4... 2 2, 4, 6... 2 ALFA 1, 2, 3... 1 1, 2, 3... 1 1, 2, 3... 1 1, 2, 3... 1 TVD 2, 4, 6... 2 2, 4, 6... 2 2, 3, 4... 2 2, 4, 6... 2 QUICK2 2, 3, 4... 2 2, 3, 4... 2 2, 3, 4... 2 2, 3, 4... 2 QUICK4 3, 4, 5... 3 2, 3, 4... 2 2, 3, 4... 2 3, 4, 5... 3 3.5 SOLUÇÃO NUMÉRICA Serão <strong>de</strong>scritos a seguir, quais foram os meios e a estrutura computacional utilizada para obtenção das soluções numéricas da equação governante. Para obtenção da solução numérica, foram utilizados programas computacionais <strong>de</strong>senvolvidos <strong>em</strong> Fortran 2003 por meio do aplicativo Intel Visual Fortran 9.1, utilizando precisão quádrupla (REAL*16). O tipo <strong>de</strong> projeto para criação do programa foi o Fortran Console Application, versão release. As simulações foram feitas <strong>em</strong> computadores do Laboratório <strong>de</strong> Experimentação Numérica (LENA-1) da Universida<strong>de</strong> Fe<strong>de</strong>ral do Paraná. Foram utilizados dois computadores, o CFD-13 que possui processador Intel Core2Quad, com velocida<strong>de</strong> <strong>de</strong> barramento <strong>de</strong> 2.4 GHz, m<strong>em</strong>ória <strong>de</strong> 8GB <strong>de</strong> RAM e Windows XP 64 bits, e o CFD-8 com processador Intel Pentium IV, 3.0 GHz, 2.0 GB RAM.
- Page 1 and 2:
UNIVERSIDADE FEDERAL DO PARANÁ PRO
- Page 3 and 4:
TERMO DE APROVAÇÃO EDUARDO MATOS
- Page 5 and 6:
AGRADECIMENTOS Agradeço ao Governo
- Page 7 and 8:
ABSTRACT This work is based on Comp
- Page 9 and 10:
Figura A.2.2 - Volumes no contorno
- Page 11 and 12:
LISTA DE ABREVIATURAS ADS AIAA ANS
- Page 13 and 14:
LISTA DE SÍMBOLOS E Pe p L p V p E
- Page 15 and 16:
3.1.1 Fenômeno Físico............
- Page 17 and 18:
B.6 ANÁLISE A PRIORI DAS ORDENS VE
- Page 19 and 20: 19 Outra forma de análise possíve
- Page 21 and 22: 21 interesse nesta área é o Accel
- Page 23 and 24: 23 (PTC-61). Já a American Nuclear
- Page 25 and 26: 25 Verificar a posteriori as ord
- Page 27 and 28: 27 elementar, o menor volume de flu
- Page 29 and 30: 29 Com base no exposto acima, o tra
- Page 31 and 32: 31 q N N D f (2.7) g 1 ond
- Page 33 and 34: 33 caso de um domínio unidimension
- Page 35 and 36: 35 Considerando este aspecto, a pri
- Page 37 and 38: 37 numericamente difusiva. Este é
- Page 39 and 40: 39 2.8.5 PLDS (Power Law Difference
- Page 41 and 42: 41 A condição principal para um e
- Page 43 and 44: 43 2.9.1 CDS-2 O esquema CDS-2 é d
- Page 45 and 46: 45 Para utilização do método, de
- Page 47 and 48: 47 2.11.4 Erros de Iteração É a
- Page 49 and 50: 49 h h h h j-2 j-1 j j+1 j+2 Figura
- Page 51 and 52: 51 Voltando às relações desenvol
- Page 53 and 54: 53 as notações “1” e “2”
- Page 55 and 56: 55 2.12.1 Estimativas de Erro a Pri
- Page 57 and 58: 57 Esta ordem pode ser obtida com d
- Page 59 and 60: 59 3 PROCEDIMENTOS MATEMÁTICOS E N
- Page 61 and 62: 61 Fluido: newtoniano, incompress
- Page 63 and 64: 63 a aproximação UDS e UDS-2. Por
- Page 65 and 66: 65 Tabela 3.4 - Valor analítico da
- Page 67 and 68: 67 Como o intuito é verificar o im
- Page 69: 69 3.4.1 Análise a Priori da Ordem
- Page 73 and 74: 73 Por meio do gerenciador de tar
- Page 75 and 76: 75 primeira requer aproximadamente
- Page 77 and 78: 77 usam processo iterativo, pode-se
- Page 79 and 80: 79 Ve-se ainda que para a malha mai
- Page 81 and 82: 81 4.3.4 Variável L Os erros anali
- Page 83 and 84: 83 influência dos erros de arredon
- Page 85 and 86: 85 Tabela 4.6 - Ordens verdadeiras,
- Page 87 and 88: 87 Tabela 4.8 - Ordens verdadeiras,
- Page 89 and 90: 89 Para o primeiro caso, pode-se ve
- Page 91 and 92: 91 discrepância aumenta, sendo que
- Page 93 and 94: 93 4.5.4 Variável L Analisando a F
- Page 95 and 96: 95 conjuntos de FI`s, ao par formad
- Page 97 and 98: 97 Analisando agora os erros em ter
- Page 99 and 100: 99 Outra análise possível é a co
- Page 101 and 102: 101 Outro ponto a observar é que a
- Page 103 and 104: 103 Ainda com base na Fig. 4.21 o c
- Page 105 and 106: 105 5 CONCLUSÃO Neste último cap
- Page 107 and 108: 107 CDS-2 permite a obtenção de e
- Page 109 and 110: 109 computacionais mais reduzidos
- Page 111 and 112: 111 FERZIGER, J.H. ; PERIC, M.. Com
- Page 113 and 114: 113 ROACHE, P.J. Recent Contributio
- Page 115 and 116: 115 Δx Δx Δx w N-2 N-1 N N+1 Fig
- Page 117 and 118: 117 temos: Aplicando inicialmente p
- Page 119 and 120: 119 w * * * * W WW W w onde w (
- Page 121 and 122:
121 A.4.4.1 Nós P=1 a P=N Aplicand
- Page 123 and 124:
123 A.4.6 ADS / CDS-2 Consideraremo
- Page 125 and 126:
125 Considerando então a forma gen
- Page 127 and 128:
127 * E P 9. P 8. W E
- Page 129 and 130:
129 * * * * * 3. E 2. P W e P
- Page 131 and 132:
131 A.4.10.3 Nó P=2 Analisando ago
- Page 133 and 134:
133 70 60 Pe. x 15. 16 a * * *
- Page 135 and 136:
135 APÊNDICE B - Análise a priori
- Page 137 and 138:
137 E a equação geral do erro de
- Page 139 and 140:
139 que pode ser reescrita da forma
- Page 141 and 142:
141 Para o desenvolvimento do traba
- Page 143 and 144:
143 j 3. j1 j2 3 ii 2 1 iii 3 1
- Page 145 and 146:
145 j i 1 ii 2 j iii 3 1 iv 4 c
- Page 147 and 148:
147 Substituiremos agora as Eqs.(B.
- Page 149 and 150:
149 Consideramos na Eq.(B.6.75) o t
- Page 151 and 152:
151 contrário, se refinarmos a mal
- Page 153 and 154:
153 i j 27. j 1 27. j 1 j 2 j 2.
- Page 155 and 156:
155 '' ''' IV p 3 p 5 p 7 I .
- Page 157 and 158:
157 E o erro de truncamento dado po
- Page 159 and 160:
159 B.7.4 QUARTA VARIÁVEL DE INTER
- Page 161 and 162:
161 Conforme dito anteriormente, pa
- Page 163 and 164:
163 B.10 ANÁLISE A PRIORI DAS ORDE
- Page 165:
165 Como E( P ) é o erro de discr