Universidade Estadual de Campinas - fem.unicamp.brphoenics/EM974/PROJETOS/PROJETOS%201SEM-09/14/... · Estimativa do Coeficiente de Convecção Natural para o ar ... pelo sistema

Universidade Estadual de Campinas Faculdade de Engenharia Mecânica

EM 974 – MÉTODOS COMPUTACIONAIS EM ENGENHARIA TÉRMICA E

AMBIENTAL.

RELATÓRIO FINAL

Daniel Portioli Sampaio RA 032017 Thiago Hideki Sato RA 046721

2

Índice 1. Tema Proposto.............................................................................3

2. Objetivos.......................................................................................3

3. Transformadores Elétricos............................................................3

4. Bancada de Testes da Empresa...................................................5

5. Fotos Termográficas.....................................................................6

6. Estimativa do Coeficiente de Convecção Natural para o ar

Exterior..........................................................................................8

7. Implementação no Phoenics.........................................................8

8. Materiais e Métodos....................................................................10

8.1 Modelo da simulação inicial ........................................................10

8.2 Modelos das simulações finais.....................................................10

8.3 Obtenção dos valores de densidade e viscosidade..........................11

9. Calculo do hmed de Convecção Natural para o Ar.....................12

10. Resultados para a Simulação Inicial.........................................13

11. Resultados para a Alteração da Simulação Inicial....................16

12. Resultados para a Simulação de Três Tubos .........................22

13. Conclusão.................................................................................24

14. Bibliografia................................................................................25

15. Anexos......................................................................................25

15.1 Modelo 14 Tubos.....................................................................25

15.2 Arquivo Q1.............................................................................26

15.3 Arquivo Results.......................................................................30

3

1. Tema proposto:

Estudo de um radiador tubular usado em transformadores de potência.

2. Objetivos:

Simular no programa Phoencis através de volumes finitos um radiador

tubular de aço, com as simplificações de modelo necessária, sujeito a condição

de troca de calor por convecção natural com o ambiente. Visa-se também

comparar os resultados obtidos com fotos termográficas das chapas exteriores

dos elementos do radiador. Este radiador e usado no resfriamento de grandes

transformadores elétricos, onde suas espiras funcionam como uma fonte de

geração de calor. Os resultados esperados são a temperatura na parede do

trocador de calor e os diagramas de velocidade e pressão.

3. Transformadores Elétricos

Os transformadores elétricos têm o objetivo de elevar ou baixar os

valores de tensão das redes de transmissão através de enrolamentos de espiras

em um núcleo ferro magnético. Entretanto, a passagem de corrente elétrica

pelas espiras gera perda de potência na forma de calor que deve ser eliminada

pelo sistema de refrigeração.

A vida útil dos transformadores termina quando, após um estágio

gradativo de debilidade, falha a sua isolação. A vida do equipamento pode ser

definida como sendo o período necessário para que o isolamento do sistema

elétrico fique seriamente deteriorado, fazendo com que a probabilidade de falha

se torne excessivamente alta.

4

Figura 1 – Radiador

Um fator determinante no tempo de vida útil de um equipamento é o tipo

de material isolante empregado em sua construção. Este, por sua vez, tem seu

processo de envelhecimento variando em função da temperatura (o

envelhecimento acelera com o aumento da temperatura) no enrolamento.

Os fatores mais importantes a serem controlados são: A umidade, a

temperatura e a presença de partículas estranhas à composição da isolação,

isto é, prevenção de depósito ou mistura destas partículas durante a vida do

equipamento.

A combinação da isolação sólida e líquida nos transformadores esta

presente em diferentes partes de sua estrutura interna, entre espiras, núcleos,

bobinas e entre estes elementos e a carcaça.

Em alguns tipos de transformadores, o ar é utilizado como uma

vantagem adicional, no intuito de se reduzir a temperatura do óleo e do isolante

sólido em redor das fontes térmicas do equipamento. Isto é, além da troca de

calor natural efetuada pelo corpo do transformador com o ambiente, por meio

de aletas (canalização externa ao corpo do transformador para circulação do

óleo), aumenta-se a superfície de contato entre o óleo circulante e o meio

externo.

5

4. Bancada de Testes da Empresa

Para a realização do experimento foram pesquisados dados

experimentais para a comprovação via simulação numérica.

A bancada de teste na qual se pretende reproduzir esta apresentado na

figura 1 abaixo.

Figura 2 – Imagem da bancada.

Sendo assim, foi possível ter uma noção das dimensões do

equipamento.

É possível ver do equipamento uma região central onde se localiza onde

será aquecido o óleo por uma resistência elétrica para simular o aquecimento

do núcleo do transformador real. Nas laterais é possível identificar dois

radiadores de 20 elementos cada. A foto abaixo esta mostrando um detalhe da

resistência elétrica que vão ser utilizada para simular o núcleo do

transformador.

6

5. Fotos Termográficas

Para comparação com os valores experimentais foram utilizadas

algumas fotos termográficas de um e 20 elementos a fim de ver as diferenças

de temperaturas com o valor ambiente.

Abaixo estão inseridas na figura 2 exemplos de fotos termográficas do

experimento proposto para um elemento.

Figura 3 – exemplos de fotos termográficas do experimento

A seguir estão inseridas algumas imagens da distribuição de

temperatura na superfície de um radiador com 20 elementos.

7

Figura 4 – Valores para a variação da viscosidade e da densidade o óleo com a

temperatura pelo aquecimento

Foi necessária a pesquisa nos livros de transferência de calor dos dados

da variação da densidade do óleo com a temperatura pela dificuldade de

obtenção de dados experimentais, além da variação da viscosidade com os

aumentos de temperatura. A tabela a seguir mostra as propriedades para óleos

usados em transformadores elétricos.

8

Tabela 1 – Propriedades do óleo.

6. Estimativa do Coeficiente de Convecção Natural para o Ar Exterior

Como o problema de resfriamento do óleo de transformadores exige a

convecção natural tanto no escoamento externo como no escoamento interno.

Este fato resultaria de uma grande complexidade numérica e esforço

computacional, a troca térmica da convecção natural do ar externo foi

aproximada por um Hmed de convecção por um calculo do coeficiente de

convecção para uma placa plana vertical.

7. Implementação no Phoenics

Para a modelagem do radiador, primeiramente foram feitas algumas

simplificações do problema. A temperatura ambiente, essencial para a troca de

calor por convecção livre, foi considerada em torno de 23º Celsius, a força

gravitacional, importante para a circulação do óleo pela diferença de

densidade, foi considerada como 9,81 N/s2. Algumas propriedades do óleo

9

utilizado foram adquiridas de dados da empresa, como já demonstrado no

trabalho anterior, e outras de livros de transferência de calor.

Para a modelagem do domínio, foram utilizadas as dimensões aproximadas

as do produto fornecido pela empresa, no caso o tamanho de cada elemento

de aproximadamente 2 metros de altura por 0,5 metros largura, para a

simulação será usado um elemento bidimensional com 2m de altura com plates

em todas as áreas laterais para dissipação de calor, a espessura do elemento.

No lado direito do domínio foi colocado um plate como fonte de calor. Foi

adicionada também no canto inferior esquerdo uma pressão de referência. Já

no meio do domínio foi criado primeiramente um cubo com propriedades

adiabáticas para tentar visualizar somente a troca térmica do radiador sem a

interferência da troca térmica do tanque com ambiente e com o elemento, alem

disso este elemento servira como a divisão do escoamento para podermos

observar o movimento do óleo devido a convecção natural e diferença de

densidade.

Radiador B (1elemento) Potência (W) 4.411,00

Altura (m) 2

N.o. Elementos 1 Tmédia_entrada 74,35 Tmédia_saída 55,81

Tamb_estimada Tóleo_médio (C 65,08

CP_médio (J/kg.K) 2075,96 rmédio (kg/m³) 861,12

m (kg/s) 0,28 m (m³/s) 3,21E-04 m (m/s) 0,06

Tabela 2 - Valores experimentais reais que se pretende reproduzir.

Tabela 2 - Tabela de valores de variação das propriedades do óleo com a temperatura

10

8. Materiais e Métodos 8.1 Modelo da Simulação Inicial

Figura 5 - Esquema do transformador simplificado a ser simulado.

Na figura o radiador esta representado no lado esquerdo (em verde) e

devera ter a mesma área de um elemento do radiador real, e o núcleo do

transformador será representado na figura como uma fonte de calor (em

vermelho no esquema). Já o quadrado branco seria um bloco com paredes

externas adiabáticas. Será necessária também, a inserção de um ponto de

referencia de pressão alem de da gravidade. O domínio devera apresentar as

propriedades do óleo e responder à diferenças de densidade em função da

temperatura. As paredes do domínio NORTH e SOUTH também serão

consideradas adiabáticas. O sistema e uma cavidade, sendo assim, não

apresentarão entrada e nem saída de escoamento.

Ficam dependentes destes resultados iniciais os trabalhos futuros em

aproximações da situação física real.

8.2 Modelos das Simulações Finais

Para as simulações visou-se o aumento das características geométricas

do problema com a inserção de mais elementos para visualizar o

comportamento de termo sifão com este acréscimo de elementos. Para isso

Entrada de potência do núcleo do transformador

Pressão de Referência

Óle

o

Ambiente externo (temperatura

ambiente)

gravidade

11

usou-se basicamente as mesmas definições descritas acima mais distanciando

os elementos 6cm com largura de tubos 7cm

Os incrementos do numero de tubos foi de 2 e 13 tubos.

8.3 Obtenção dos Valores de Densidade e Viscosidade

Com os valores da tabela 1 foi possível a construção do gráfico abaixo

da variação da densidade do óleo pela variação da temperatura.

Densidade x temperaturay = -0,6152x + 891,42

860

870

880

890

900

910

920

930

-60 -40 -20 0 20 40 60

Temperatura (ºC)

Den

sid

ade

(Kg

/m³)

Gráfico 1 – Variação da densidade do óleo com a temperatura.

A seguir esta apresentada o gráfico da variação da viscosidade com a

temperatura. Foi adotada uma equação exponencial para inserção dos dados

no equacionamento do Phoenics. Como nosso problema está numa faixa entre

20 e 80ºC, foi utilizado um valor constante para a viscosidade no valor de

1.08e-6 m2/s.

12

Viscosidade x Temperatura

0

5000

10000

15000

20000

25000

30000

35000

-60 -40 -20 0 20 40 60

temperatura(ºC)

Vis

cosi

dad

e (m

²/s)

Gráfico 2 – Variação da viscosidade do óleo com a temperatura.

9. Calculo do hmed de Convecção Natural para o Ar

Para o calculo do valor médio do coeficiente de convecção natural foi

usado o seguinte programa de matlab;

clear all

close all

clc

L=2;

tsup=80; %temperatura da superfice média

tinf=23; %temperatura do ar

tfilme=(tsup+tinf)/2; %temperatura de filme =51.5C ou 324.5

g=9.81; %gravidade

alfa=26.2e-6;

v=18.4e-6;

beta=1/(tfilme+273);

Pr=0.707;

k=28.15e-3;

Ral=(g*beta*(tsup-tinf)*(L^3))/(v*alfa) %calculo de numero de Rayleigh=19.6292

Nul=0.68+(0.67*(Ral^0.25))/((1+(0.492/Pr)^.56)^.44) %=1.7639

H=Nul*k/L %=6

13

10. Resultados para a Simulação Inicial

Com os valores dos parâmetros das equações acima obtidas dos

gráficos foram adicionados no calculo da variação da densidade modelo pela

temperatura. Sendo assim, foi obtida a distribuição dos valores de densidade

para o radiador de um elemento.

Figura 6 – Variação da Densidade.

É possível verificar pela distribuição da densidade acima, que onde foi

adicionada a fonte de calor houve uma redução dos valores densidade como

esperado pela curva dos valores tabelados. Também foi possível verificar que

nos tubos coletores apresentam valores intermediários de densidade e na

chapa do radiador foi possível verificar os maiores valores de densidade. Um

problema identificado foi que pelas aproximações da geometria por um anel de

cantos vivos houve uma concentração de altos valores de densidade nestes

locais, fato que caracteriza a física do problema e pode ser uma fonte de erros.

A seguir esta apresentada a distribuição das velocidades do problema,

como se trata de convecção natural Acredita-se em baixos números de

Reynolds foi adotado o regime laminar de velocidades no modelo de

turbulência.

14

Figura 7 – Diagrama de Velocidade.

Foi possível verificar a recirculação do fluido de trabalho no domínio

proposto, este fato, deu uma maior confiança na veracidade da física do

problema esta sendo satisfeita pelo modelo. Também foi verificadas regiões de

velocidades próximas de 0m/s nos cantos vivos do domínio pelas

simplificações feitas não deve ser a física real do radiador

15

Figura 8 – Diagrama de Temperatura.

Já para o diagrama de temperaturas inserido acima, houve o

aquecimento das paredes onde foi colocada a fonte de calor e na parede onde

esta sendo adotada os tubos do radiador com convecção houve a diminuição

dos valores desta, como esperado pelo modelo. Entretanto, estávamos

esperando uma maior temperatura na parte superior direita do modelo, fato que

deve ser investigado.

Abaixo esta inserida o diagrama de pressões para o problema

Figura 9 –Diagrama de Pressão.

16

11. Resultados para a Alteração da Simulação Inicial

Na simulação inicial não foi possível inserir os valores de potencia

dissipada pelo radiador no problema, pois esta apresentou valores irreais

para temperatura.

Acredita-se apesar do óleo apresentar uma recirculação no termo sifão,

os valores de temperatura não condiziam com os valores apresentados

pelas fotos termográficas na revisão bibliográfica. Sendo assim retornou -

se na definição do problema para tentar resolver o problema da distribuição

da temperatura. Para isso foi alterada a malha computacional para 300 o

numero de elementos, conforme pode ser visto na figura abaixo.

Figura 10 – Malha simples.

Foi possível perceber do gráfico da variação da viscosidade com a

temperatura que a curva conduzia a um valor constante na faixa de operação

do problema.

Com esta alteração para valores constantes de viscosidade foi possível

obter as seguintes figuras de velocidade, temperatura , pressão

17


Este diagrama de temperaturas apesar de apresentar uma malha mais

grosseira apresentou valores mais condizentes com a resposta esperada, com

a região superior com maior temperatura e a região inferior com menores

temperaturas ate a região próxima a da resistência elétrica.

Figura 12 – Diagrama de Velocidade

O resultado para o gráfico de velocidades apresentou valores maiores

de velocidade que a simulação do primeiro teste

18

Figura 13 – Diagrama de Pressão.

Também houve uma melhora nos valores de pressão que passou a ter

valores positivos como o esperado.

Após este testes foi feito um primeiro refinamento na malha do domínio

que passou a ser de 6300 elementos que pode seu vista na figura abaixo

Figura 14 – Malha com primeiro refinamento.

Com esta malha foi obtido a seguintes figuras de velocidade, pressão e

temperatura.

19


O diagrama de temperatura se repetiu sem maiores alterações


Entretanto, o diagrama de velocidades apresentou uma grande alteração

principalmente na região próximas as paredes, apresentando valores máximos

da ordem de duas vezes os valores da simulação anterior.

Abaixo esta inserida a figura do diagrama de pressões para o primeiro

refinamento de malha.

20

Figura 17 – Diagrama de Pressão.

Como esperado os diagramas de pressão não apresentaram uma

alteração muito grande.

Como houve uma grande alteração nos valores de velocidade foi

proposto um maior refinamento de malha que chegou a valores 24000

elementos.

Figura 18 – Malha após segundo refinamento.

21

A seguir estão apresentados os diagramas de velocidades, pressões e

temperatura para um segundo refinamento.

Figura 19 –Diagrama de Temperatura.

Houve certa alteração dos valores de temperatura na parte inferior do

escoamento devido a recirculação do óleo no canto inferior direito.


22

Do diagrama de velocidades acima é possível comprovar a recirculação na

parte inferior direita do escoamento e também na parte superior central. Este

fenômeno pode não condizer com a situação real pois existe uma simplificação

nesta região onde caracterizaria o escoamento dentro de um tubo.

Figura 21 – Diagrama de Pressão

12. Resultados para a Simulação de Três Tubos

Para uma simulação com 3 tubos, foram feitos também três tipos de malha. O primeiro modelo conta com 380 elementos, o segundo com 3500 e o terceiro com 25500. Os resultados obtidos são mostrados a seguir.

23

Figura 22 – Malhas para o modelo de 3 tubos.

Figura 23 – Diagramas de Temperatura.

Figura 24 – Diagramas de Velocidade.

24

Figura 25 – Diagramas de Pressão.

13. Conclusão

Pode-se concluir que como se trata de um problema de convecção

natural, este demanda de muito tempo para sua solução. Quanto mais se

refinava a malha, maior era o erro final nos resultados de temperatura, porém

mais próximo da temperatura real se aproximava. Para uma simulação

adequada, seriam necessárias horas de simulação para uma convergência

aceitável.

Porém foi observado que um gradiente de temperatura é formado de

acordo com o esperado, somente seus valores que não estão totalmente de

acordo devido a convergência do problema. Para o diagrama de velocidade,

também foi observado que quanto maior o refinamento da malha, pior é o

resultado devido à baixa convergência do problema. Já para a pressão, não

foram obtidas grandes diferenças variando sua malha.

Foram testados também outros modelos, um com 14 elementos, mas

devido a não convergência dos resultados e o enorme tempo de simulação

requerido, foram feitas poucas análises em cima desse modelo. Alguns

resultados desse modelo podem ser visto no anexo.

25

14. Bibliografia

Frank K., Bohn, M. S., Principio de transferência de calor, Ed. Pioneira

Thomson Learning, São Paulo , SP 2003.

15. Anexos

15.1 Modelo 14 Tubos

Figura 26 – Malha Simples para modelo com vários tubos.

Figura 27 – Diagrama de pressão para modelo com vários tubos.

26

Figura 28 – Diagrama de temperatura para modelo com vários tubos.

Figura 29 – Diagrama de velocidade para modelo com vários tubos.

15.2 Arquivo Q1

<html><head><title>Q1</title>

<link rel="stylesheet" type="text/css"

href="/phoenics/d_polis/polstyle.css">

</head><body><pre><strong>

TALK=T;RUN( 1, 1)

************************************************************

Q1 created by VDI menu, Version 3.6, Date

11/05/06

CPVNAM=VDI;SPPNAM=Core

************************************************************

IRUNN = 1 ;LIBREF = 0

************************************************************

Group 1. Run Title

TEXT(Convecção Natural )

************************************************************

Group 2. Transience

STEADY = T

27

************************************************************

Groups 3, 4, 5 Grid Information

* Overall number of cells,

RSET(M,NX,NY,NZ,tolerance)

RSET(M,150,170,1)

************************************************************

Group 6. Body-Fitted coordinates

************************************************************

Group 7. Variables: STOREd,SOLVEd,NAMEd

ONEPHS = T

* Non-default variable names

NAME(143) =STAN ; NAME(144) =HTCO

NAME(145) =PRPS ; NAME(146) =TEM1

NAME(147) =ENUL ; NAME(148) =DEN1

NAME(149) =KOND ; NAME(150) =SPH1

* Solved variables list

SOLVE(P1 ,U1 ,V1 ,TEM1)

* Stored variables list

STORE(SPH1,KOND,DEN1,ENUL,PRPS,HTCO,STA

N)

* Additional solver options

SOLUTN(P1 ,Y,Y,Y,N,N,Y)

SOLUTN(TEM1,Y,Y,Y,N,N,Y)

************************************************************

Group 8. Terms & Devices

NEWRH1 = T

************************************************************

Group 9. Properties

PRESS0 = 1.000000E+05 ;TEMP0 =

2.730000E+02

RHO1 = GRND4

RHO1A = 8.620000E+02 ;RHO1B =-

6.152000E-01

RHO1C = 1.600000E+01

ENUL = 1.080000E-06

ENULA =-3.116500E+03 ;ENULB =-

1.215800E+02

ENULC = 7.420000E+00

CP1 = 2.075000E+03

ENUT = 0.000000E+00

DVO1DT = 1.180000E-04

PRNDTL(TEM1) = -1.100000E-01

************************************************************

Group 10.Inter-Phase Transfer Processes

************************************************************

Group 11.Initialise Var/Porosity Fields

FIINIT(PRPS) = -1.000000E+00 ;FIINIT(TEM1) =

1.000000E+02

No PATCHes used for this Group

INIADD = T

************************************************************

Group 12. Convection and diffusion adjustments


************************************************************

Group 13. Boundary & Special Sources

PATCH (BUOYANCY,PHASEM,0,0,0,0,0,0,1,1)

COVAL (BUOYANCY,U1 , FIXFLU , GRND2

)

COVAL (BUOYANCY,V1 , FIXFLU , GRND2

)

BUOYA = 0.000000E+00 ; BUOYB =-

9.810000E+00

BUOYC = 0.000000E+00

BUOYD = 8.620000E+02

BUOYE = 2.700000E+01

EGWF = T

************************************************************

Group 14. Downstream Pressure For PARAB

************************************************************

Group 15. Terminate Sweeps

LSWEEP = 10000

RESFAC = 1.000000E-03

************************************************************

Group 16. Terminate Iterations

LITER (P1 ) = 50 ;LITER (TEM1) = 50

************************************************************

Group 17. Relaxation

RELAX(P1 ,LINRLX, 1.000000E+00)

RELAX(TEM1,FALSDT, 1.000000E-01)

************************************************************

Group 18. Limits

28

VARMAX(U1 ) = 1.000000E+02 ;VARMIN(U1 )

=-1.000000E+02

VARMAX(V1 ) = 1.000000E+02 ;VARMIN(V1 )

=-1.000000E+02

VARMAX(TEM1) = 8.500000E+01

;VARMIN(TEM1) = 2.000000E+01

************************************************************

Group 19. EARTH Calls To GROUND Station

USEGRD = T ;USEGRX = T

ASAP = T

PARSOL = F

CONWIZ = T

SPEDAT(SET,OUTPUT,TECPLOT,C,YES)

************************************************************

Group 20. Preliminary Printout

ECHO = T

************************************************************

Group 21. Print-out of Variables

OUTPUT(U1 ,Y,N,N,Y,Y,Y)

OUTPUT(V1 ,Y,N,N,Y,Y,Y)

OUTPUT(ENUL,Y,N,N,N,N,N)

OUTPUT(KOND,Y,N,N,N,N,N)

OUTPUT(SPH1,Y,N,N,N,N,N)

************************************************************

Group 22. Monitor Print-Out

IXMON = 147 ;IYMON = 81 ;IZMON =

1

NPRMON = 100000

NPRMNT = 1

TSTSWP = -1

************************************************************

Group 23.Field Print-Out & Plot Control

NPRINT = 100000

ISWPRF = 1 ;ISWPRL = 100000


************************************************************

Group 24. Dumps For Restarts

GVIEW(P,0.000000E+00,0.000000E+00,1.000000E+0

0)

GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E

+00)

> DOM, SIZE, 6.000000E-01,

2.000000E+00, 1.000000E-01

> DOM, MONIT, 5.850000E-01, 9.500000E-

01, 5.000000E-02

> DOM, SCALE, 1.000000E+00,

1.000000E+00, 1.000000E+00

> DOM, SNAPSIZE, 1.000000E-02

> GRID, AUTO, F F T

> GRID, RSET_X_1, 37, 1.000000E+00

> GRID, RSET_X_2, 25, 1.000000E+00

> GRID, RSET_X_3, 18, 1.000000E+00

> GRID, RSET_X_4, 15, 1.000000E+00

> GRID, RSET_X_5, 18, 1.000000E+00

> GRID, RSET_X_6, 15, 1.000000E+00

> GRID, RSET_X_7, 22, 1.000000E+00

> GRID, RSET_Y_1, 9, 1.000000E+00

> GRID, RSET_Y_2, 51, 1.000000E+00

> GRID, RSET_Y_3, 50, 1.000000E+00

> GRID, RSET_Y_4, 51, 1.000000E+00

> GRID, RSET_Y_5, 9, 1.000000E+00

> GRID, RSET_Z_1, 1, 1.000000E+00

> OBJ, NAME, PAREDE1

> OBJ, POSITION, 1.500000E-01,

1.000000E-01, 0.000000E+00

> OBJ, SIZE, 1.000000E-01,

1.800000E+00, 1.000000E-01

> OBJ, GEOMETRY, cube

> OBJ, ROTATION24, 1

> OBJ, TYPE, BLOCKAGE

> OBJ, MATERIAL, 198,Solid with smooth-

wall friction

> OBJ, LINR_HEVT_E, 6.000000E+00,

2.300000E+01

> OBJ, NAME, B2


4.444445E-02, 0.000000E+00

> OBJ, SIZE, 4.054054E-03, 1.111111E-

02, 1.000000E-01

> OBJ, GEOMETRY, cubet


> OBJ, TYPE, PRESSURE_RELIEF

> OBJ, PRES_RELIEF, 1.000000E+04,

1.000000E+01

> OBJ, NAME, ESQUERDA


0.000000E+00, 0.000000E+00

> OBJ, SIZE, 0.000000E+00,

2.000000E+00, 1.000000E-01

29

> OBJ, GEOMETRY, cube13


> OBJ, TYPE, PLATE

> OBJ, LINR_HEVT, 6.000000E+00,

2.300000E+01

> OBJ, NAME, ENTRADA

> OBJ, POSITION, 0.000000E+00,

7.000000E-01, 0.000000E+00

> OBJ, SIZE, 0.000000E+00, 6.000000E-

01, 1.000000E-01



> OBJ, TYPE, PLATE

> OBJ, SURF_HEVT, 0.000000E+00,

4.400000E+03

> OBJ, NAME, CHãO


0.000000E+00, 0.000000E+00

> OBJ, SIZE, 6.000000E-01,

0.000000E+00, 1.000000E-01



> OBJ, TYPE, PLATE

> OBJ, NAME, TETO


2.000000E+00, 0.000000E+00

> OBJ, SIZE, 6.000000E-01,

0.000000E+00, 1.000000E-01



> OBJ, TYPE, PLATE

> OBJ, NAME, DIREITAB


0.000000E+00, 0.000000E+00

> OBJ, SIZE, 0.000000E+00, 7.000000E-

01, 1.000000E-01



> OBJ, TYPE, PLATE

> OBJ, NAME, DIREITAA


1.300000E+00, 0.000000E+00

> OBJ, SIZE, 0.000000E+00, 7.000000E-

01, 1.000000E-01



> OBJ, TYPE, PLATE

> OBJ, NAME, PARADE2


1.000000E-01, 0.000000E+00

> OBJ, SIZE, 6.000000E-02,

1.800000E+00, 1.000000E-01





wall friction

> OBJ, LINR_HEVT_W, 6.000000E+00,

2.300000E+01

> OBJ, LINR_HEVT_E, 6.000000E+00,

2.300000E+01

> OBJ, NAME, B10


1.000000E-01, 0.000000E+00

> OBJ, SIZE, 6.000000E-02,

1.800000E+00, 1.000000E-01





wall friction

> OBJ, LINR_HEVT_W, 6.000000E+00,

2.300000E+01

> OBJ, LINR_HEVT_E, 6.000000E+00,

2.300000E+01

STOP

</strong></pre></body></html>

30

15.3 Arquivo Results ************************************************************ --------------------------------------------------------- CCCC HHH PHOENICS June 2006 - EARTH CCCCCCCC H (C) Copyright 2006 CCCCCCC See H Concentration Heat and Momentum Ltd CCCCCCC our new H All rights reserved. CCCCCC Web-site H Address: Bakery House, 40 High St CCCCCCC www.cham. H Wimbledon, London, SW19 5AU CCCCCCC co.uk H Tel: 020-8947-7651 CCCCCCCC H Fax : 020-8879-3497 CCCC HHH E-mail: [email protected] --------------------------------------------------------- This program forms part of the PHOENICS installation for: CHAM The code expiry date is the end of : may 2037 --------------------------------------------------------- ************************************************************ Information about material properties Total number of SPEDATs is 51 number of materials specified by SPEDATs is 1 solprp = 100 porprp = 198 vacprp = 199 !!!! The properties file is PROPS Properties being read from PROPS Properties have been read from PROPS Property-related data from gxprutil: PRPS is stored with initial value = =-1.000000E+00 Material properties used are... denst1 visclm tempr1 thrme1 speht1 also, other related settings are ... USEGRX = T USEGRD = T >>> End of property-related data <<< ************************************************************ Number of F-array locations available is 10000000 Number used before BFC allowance is 2430852 Number used after BFC allowance is 2430852 biggest cell volume divided by average is 1.043173 at: ix = 149 iy = 110 iz = 1 xg = 5.938637E-01 yg = 1.294000E+00 zg = 5.000000E-02 smallest cell volume divided by average is 0.9181433 at: ix = 79 iy = 166 iz = 1 xg = 3.141667E-01 yg = 1.950000E+00 zg = 5.000000E-02 ratio of smallest to biggest is 0.8801450 false-time relaxation for TEM1may be unwise. Try linrlx,0.5 instead.

************************************************************ -------- Recommended settings ------- CONWIZ = T activates settings based on refrho = 1.000000E+00 refvel = 1.000000E+01 reflen = 1.000000E+00 reftemp = 1.000000E+03 rlxdu1 = 5.000000E-01 rlxdv1 = 5.000000E-01 rlxdw1 = 5.000000E-01 Maximum change of U1 per sweep = 100.0000 Maximum change of V1 per sweep = 100.0000 Maximum change of TEM1 per sweep = 1000.000 relaxation and min/max values left at defaults may have been changed ************************************************************ ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(Convecção Natural ) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = T NX = 150 XULAST = 6.000000E-01 XFRAC ( 1) = 6.756757E-03 ;XFRAC ( 31) = 2.094595E-01 XFRAC ( 61) = 4.100000E-01 ;XFRAC ( 91) = 6.066666E-01 XFRAC (121) = 8.033333E-01 ************************************************************ Group 4. Y-Direction Grid Spacing NY = 170 YVLAST = 2.000000E+00 YFRAC ( 1) = 5.555556E-03 ;YFRAC ( 35) = 2.029412E-01 YFRAC ( 69) = 4.040000E-01 ;YFRAC (103) = 6.080000E-01 YFRAC (137) = 8.088235E-01 ************************************************************ Group 5. Z-Direction Grid Spacing PARAB = F NZ = 1 ZWLAST = 1.000000E-01 ************************************************************ Group 6. Body-Fitted Coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T NAME( 1) =P1 ;NAME( 3) =U1 NAME( 5) =V1 ;NAME(143) =STAN NAME(144) =HTCO ;NAME(145) =PRPS NAME(146) =TEM1 ;NAME(147) =ENUL NAME(148) =DEN1 ;NAME(149) =KOND NAME(150) =SPH1 * Y in SOLUTN argument list denotes: * 1-stored 2-solved 3-whole-field * 4-point-by-point 5-explicit 6-harmonic averaging SOLUTN(P1 ,Y,Y,N,N,N,Y) SOLUTN(U1 ,Y,Y,N,N,N,Y) SOLUTN(V1 ,Y,Y,N,N,N,Y) SOLUTN(STAN,Y,N,N,N,N,Y) SOLUTN(HTCO,Y,N,N,N,N,Y)

31

SOLUTN(PRPS,Y,N,N,N,N,Y) SOLUTN(TEM1,Y,Y,N,N,N,Y) SOLUTN(ENUL,Y,N,N,N,N,Y) SOLUTN(DEN1,Y,N,N,N,N,Y) SOLUTN(KOND,Y,N,N,N,N,Y) SOLUTN(SPH1,Y,N,N,N,N,Y) DEN1 = 148 VISL = 147 PRPS = 145 ************************************************************ Group 8. Terms & Devices * Y in TERMS argument list denotes: * 1-built-in source 2-convection 3-diffusion 4-transient * 5-first phase variable 6-interphase transport TERMS (P1 ,Y,Y,Y,N,Y,N) TERMS (U1 ,Y,Y,Y,N,Y,N) TERMS (V1 ,Y,Y,Y,N,Y,N) TERMS (TEM1,N,Y,Y,N,Y,N) DIFCUT = 5.000000E-01 ;ZDIFAC = 1.000000E+00 GALA = F ;ADDDIF = T NEWRH1 = T ISOLX = -1 ;ISOLY = -1 ;ISOLZ = 0 ************************************************************ Group 9. Properties used if PRPS is not stored, and where PRPS = -1.0 if it is! RHO1 = GRND4 ;TMP1 = 0.000000E+00 EL1 = 0.000000E+00 TSURR = 0.000000E+00 ;TEMP0 = 2.730000E+02 PRESS0 = 1.000000E+05 DVO1DT = 1.180000E-04 ;DRH1DP = 0.000000E+00 RHO1A = 8.620000E+02 ;RHO1B =-6.152000E-01 RHO1C = 1.600000E+01 EMISS = 0.000000E+00 ;SCATT = 0.000000E+00 RADIA = 0.000000E+00 ;RADIB = 0.000000E+00 ENUL = 1.080000E-06 ;ENUT = 0.000000E+00 PRNDTL(U1 ) = 1.000000E+00 ;PRNDTL(V1 ) = 1.000000E+00 PRNDTL(TEM1) = -1.100000E-01 PRT (U1 ) = 1.000000E+00 ;PRT (V1 ) = 1.000000E+00 PRT (TEM1) = 1.000000E+00 CP1 = 2.075000E+03 ;CP2 = 1.000000E+00 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initial field variables (PHIs) FIINIT(P1 ) = 1.000000E-10 ;FIINIT(U1 ) = 1.000000E-10 FIINIT(V1 ) = 1.000000E-10 ;FIINIT(STAN) = 1.000000E-10 FIINIT(HTCO) = 1.000000E-10 ;FIINIT(PRPS) = -1.000000E+00 FIINIT(TEM1) = 1.000000E+02 ;FIINIT(ENUL) = 1.080000E-06 FIINIT(DEN1) = 1.000000E-10 ;FIINIT(KOND) = 1.000000E-10 FIINIT(SPH1) = 1.000000E-10 Parent VR object for this patch is: PAREDE1 PATCH(OB1 ,INIVAL, 38, 62, 10, 161, 1, 1, 1, 1) INIT(OB1 ,PRPS, 0.000000E+00, 1.980000E+02) Parent VR object for this patch is: PARADE2 PATCH(OB9 ,INIVAL, 81, 95, 10, 161, 1, 1, 1, 1) INIT(OB9 ,PRPS, 0.000000E+00, 1.980000E+02) Parent VR object for this patch is: B10 PATCH(OB10 ,INIVAL, 114, 128, 10, 161, 1, 1, 1, 1) INIT(OB10 ,PRPS, 0.000000E+00, 1.980000E+02) INIADD = T FSWEEP = 1

NAMFI =CHAM ************************************************************ Group 12. Patchwise adjustment of terms Patches for this group are printed with those for Group 13. Their names begin either with GP12 or & ************************************************************ Group 13. Boundary & Special Sources Parent VR object for this patch is: B2 PATCH(P-REL1 ,CELL , 3, 5, 4, 6, 1, 1, 1, 1) COVAL(P-REL1 ,P1 , 1.000000E+04, 1.000000E+01) COVAL(P-REL1 ,U1 , 0.000000E+00, 0.000000E+00) COVAL(P-REL1 ,V1 , 0.000000E+00, 0.000000E+00) COVAL(P-REL1 ,TEM1, 0.000000E+00, SAME ) Parent VR object for this patch is: ESQUERDA PATCH(OB3 ,EWALL , 150, 150, 1, 170, 1, 1, 1, 1) COVAL(OB3 ,V1 , 1.000000E+00, 0.000000E+00) Parent VR object for this patch is: ENTRADA PATCH(OB4 ,WWALL , 1, 1, 61, 110, 1, 1, 1, 1) COVAL(OB4 ,V1 , 1.000000E+00, 0.000000E+00) COVAL(OB4 ,TEM1, FIXFLU , 4.400000E+03) Parent VR object for this patch is: CHãO PATCH(OB5 ,SWALL , 1, 150, 1, 1, 1, 1, 1, 1) COVAL(OB5 ,U1 , 1.000000E+00, 0.000000E+00) Parent VR object for this patch is: TETO PATCH(OB6 ,NWALL , 1, 150, 170, 170, 1, 1, 1, 1) COVAL(OB6 ,U1 , 1.000000E+00, 0.000000E+00) Parent VR object for this patch is: DIREITAB PATCH(OB7 ,WWALL , 1, 1, 1, 60, 1, 1, 1, 1) COVAL(OB7 ,V1 , 1.000000E+00, 0.000000E+00) Parent VR object for this patch is: DIREITAA PATCH(OB8 ,WWALL , 1, 1, 111, 170, 1, 1, 1, 1) COVAL(OB8 ,V1 , 1.000000E+00, 0.000000E+00) PATCH(BUOYANCY,PHASEM, 1, 150, 1, 170, 1, 1, 1, 1) COVAL(BUOYANCY,U1 , FIXFLU , GRND2 ) COVAL(BUOYANCY,V1 , FIXFLU , GRND2 ) Parent VR object for this patch is: PAREDE1 PATCH(OC1_E ,WEST , 63, 63, 10, 161, 1, 1, 1, 1) COVAL(OC1_E ,TEM1, 6.000000E+00, 2.300000E+01) Parent VR object for this patch is: ESQUERDA PATCH(OC3 ,EAST , 150, 150, 1, 170, 1, 1, 1, 1) COVAL(OC3 ,TEM1, 6.000000E+00, 2.300000E+01) Parent VR object for this patch is: PARADE2 PATCH(OC9_W ,EAST , 80, 80, 10, 161, 1, 1, 1, 1) COVAL(OC9_W ,TEM1, 6.000000E+00, 2.300000E+01) Parent VR object for this patch is: PARADE2

32

PATCH(OC9_E ,WEST , 96, 96, 10, 161, 1, 1, 1, 1) COVAL(OC9_E ,TEM1, 6.000000E+00, 2.300000E+01) Parent VR object for this patch is: B10 PATCH(OC10_W ,EAST , 113, 113, 10, 161, 1, 1, 1, 1) COVAL(OC10_W ,TEM1, 6.000000E+00, 2.300000E+01) Parent VR object for this patch is: B10 PATCH(OC10_E ,WEST , 129, 129, 10, 161, 1, 1, 1, 1) COVAL(OC10_E ,TEM1, 6.000000E+00, 2.300000E+01) XCYCLE = F EGWF = F BUOYA = 0.000000E+00 ; BUOYB =-9.810000E+00 BUOYC = 0.000000E+00 BUOYD = 8.620000E+02 BUOYE = 2.700000E+01 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 10000 ;ISWC1 = 1 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1 SELREF = T RESFAC = 1.000000E-03 ************************************************************ Group 16. Terminate Iterations LITER (P1 ) = 50 ;LITER (U1 ) = 10 LITER (V1 ) = 10 ;LITER (TEM1) = 50 ENDIT (P1 ) = 1.000000E-03 ;ENDIT (U1 ) = 1.000000E-03 ENDIT (V1 ) = 1.000000E-03 ;ENDIT (TEM1) = 1.000000E-03 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX, 5.000000E-01) RELAX(U1 ,LINRLX, 5.000000E-01) RELAX(V1 ,LINRLX, 5.000000E-01) RELAX(STAN,LINRLX, 1.000000E+00) RELAX(HTCO,LINRLX, 1.000000E+00) RELAX(PRPS,LINRLX, 1.000000E+00) RELAX(TEM1,FALSDT, 1.000000E-01) RELAX(ENUL,LINRLX, 1.000000E+00) RELAX(DEN1,LINRLX, 5.000000E-01) RELAX(KOND,LINRLX, 1.000000E+00) RELAX(SPH1,LINRLX, 1.000000E+00) OVRRLX = 0.000000E+00 EXPERT = F ;NNORSL = F ************************************************************ Group 18. Limits VARMAX(P1 ) = 1.000000E+10 ;VARMIN(P1 ) =-1.000000E+10 VARMAX(U1 ) = 1.000000E+02 ;VARMIN(U1 ) =-1.000000E+02 VARMAX(V1 ) = 1.000000E+02 ;VARMIN(V1 ) =-1.000000E+02 VARMAX(STAN) = 1.000000E+10 ;VARMIN(STAN) =-1.000000E+10 VARMAX(HTCO) = 1.000000E+10 ;VARMIN(HTCO) =-1.000000E+10 VARMAX(PRPS) = 1.000000E+10 ;VARMIN(PRPS) =-1.000000E+10 VARMAX(TEM1) = 8.500000E+01 ;VARMIN(TEM1) = 2.000000E+01 VARMAX(ENUL) = 1.000000E+10 ;VARMIN(ENUL) = 1.000000E-10 VARMAX(DEN1) = 1.000000E+10 ;VARMIN(DEN1) = 1.000000E-06 VARMAX(KOND) = 1.000000E+10 ;VARMIN(KOND) =-1.000000E+10

VARMAX(SPH1) = 1.000000E+10 ;VARMIN(SPH1) = 1.000000E-10 ************************************************************ Group 19. Data transmitted to GROUND USEGRD = T ;USEGRX = T ASAP = T PARSOL = F CONWIZ = T SPEDAT(SET,OUTPUT,TECPLOT,C,YES) SPEDAT(SET,OBJNAM,!OB1,C,PAREDE1) SPEDAT(SET,OBJTYP,!OB1,C,BLOCKAGE) SPEDAT(SET,PAREDE1,MATERIAL,R,1.98000E+02) SPEDAT(SET,OBJNAM,!OC1_E,C,PAREDE1) SPEDAT(SET,OBJTYP,!OC1_E,C,USER_DEFINED) SPEDAT(SET,OBJNAM,!OC1_E,C,PAREDE1) SPEDAT(SET,OBJTYP,!OC1_E,C,USER_DEFINED) SPEDAT(SET,B2,DATFILE,C,cubet) SPEDAT(SET,OBJNAM,^OB2,C,B2) SPEDAT(SET,OBJTYP,^OB2,C,PRESSURE_RELIEF) SPEDAT(SET,OBJNAM,^P-REL1,C,B2) SPEDAT(SET,OBJTYP,^P-REL1,C,USER_DEFINED) SPEDAT(SET,OBJNAM,!OB3,C,ESQUERDA) SPEDAT(SET,OBJTYP,!OB3,C,PLATE) SPEDAT(SET,OBJNAM,!OC3,C,ESQUERDA) SPEDAT(SET,OBJTYP,!OC3,C,USER_DEFINED) SPEDAT(SET,OBJNAM,!OB4,C,ENTRADA) SPEDAT(SET,OBJTYP,!OB4,C,PLATE) SPEDAT(SET,OBJNAM,!OB5,C,CHãO) SPEDAT(SET,OBJTYP,!OB5,C,PLATE) SPEDAT(SET,OBJNAM,!OB6,C,TETO) SPEDAT(SET,OBJTYP,!OB6,C,PLATE) SPEDAT(SET,OBJNAM,!OB7,C,DIREITAB) SPEDAT(SET,OBJTYP,!OB7,C,PLATE) SPEDAT(SET,OBJNAM,!OB8,C,DIREITAA) SPEDAT(SET,OBJTYP,!OB8,C,PLATE) SPEDAT(SET,OBJNAM,!OB9,C,PARADE2) SPEDAT(SET,OBJTYP,!OB9,C,BLOCKAGE) SPEDAT(SET,PARADE2,MATERIAL,R,1.98000E+02) SPEDAT(SET,OBJNAM,!OC9_W,C,PARADE2) SPEDAT(SET,OBJTYP,!OC9_W,C,USER_DEFINED) SPEDAT(SET,OBJNAM,!OC9_W,C,PARADE2) SPEDAT(SET,OBJTYP,!OC9_W,C,USER_DEFINED) SPEDAT(SET,OBJNAM,!OC9_E,C,PARADE2) SPEDAT(SET,OBJTYP,!OC9_E,C,USER_DEFINED) SPEDAT(SET,OBJNAM,!OC9_E,C,PARADE2) SPEDAT(SET,OBJTYP,!OC9_E,C,USER_DEFINED) SPEDAT(SET,OBJNAM,!OB10,C,B10) SPEDAT(SET,OBJTYP,!OB10,C,BLOCKAGE) SPEDAT(SET,B10,MATERIAL,R,1.98000E+02) SPEDAT(SET,OBJNAM,!OC10_W,C,B10) SPEDAT(SET,OBJTYP,!OC10_W,C,USER_DEFINED) SPEDAT(SET,OBJNAM,!OC10_W,C,B10) SPEDAT(SET,OBJTYP,!OC10_W,C,USER_DEFINED) SPEDAT(SET,OBJNAM,!OC10_E,C,B10) SPEDAT(SET,OBJTYP,!OC10_E,C,USER_DEFINED) SPEDAT(SET,OBJNAM,!OC10_E,C,B10) SPEDAT(SET,OBJTYP,!OC10_E,C,USER_DEFINED) SPEDAT(SET,FACETDAT,NUMOBJ,I,10) SPEDAT(SET,MATERIAL,198,L,T) ************************************************************ Group 20. Preliminary Printout ECHO = T ************************************************************ Group 21. Print-out of Variables INIFLD = F ;SUBWGR = F * Y in OUTPUT argument list denotes:

33

* 1-field 2-correction-eq. monitor 3-selective dumping * 4-whole-field residual 5-spot-value table 6-residual table OUTPUT(P1 ,Y,N,Y,Y,Y,Y) OUTPUT(U1 ,Y,N,N,Y,Y,Y) OUTPUT(V1 ,Y,N,N,Y,Y,Y) OUTPUT(STAN,Y,N,Y,N,N,N) OUTPUT(HTCO,Y,N,Y,N,N,N) OUTPUT(PRPS,Y,N,Y,N,N,N) OUTPUT(TEM1,Y,N,Y,Y,Y,Y) OUTPUT(ENUL,Y,N,N,N,N,N) OUTPUT(DEN1,Y,N,Y,N,N,N) OUTPUT(KOND,Y,N,N,N,N,N) OUTPUT(SPH1,Y,N,N,N,N,N) WALPRN = T ************************************************************ Group 22. Monitor Print-Out IXMON = 147 ;IYMON = 81 ;IZMON = 1 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = 10001 UWATCH = F ;USTEER = F HIGHLO = F ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 10000 ;NUMCLS = 5 NXPRIN = 30 ;IXPRF = 1 ;IXPRL = 150 NYPRIN = 34 ;IYPRF = 1 ;IYPRL = 170 IPLTF = 1 ;IPLTL = 10000 ;NPLT = 500 ISWPRF = 1 ;ISWPRL = 100000 ITABL = 3 ;IPROF = 1 ABSIZ = 5.000000E-01 ;ORSIZ = 4.000000E-01 NTZPRF = 1 ;NCOLPF = 50 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20 No PATCHes yet used for this Group ************************************************************ Group 24. Dumps For Restarts SAVE = T ;NOWIPE = F NSAVE =CHAM *** grid-geometry information *** X-coordinates of the cell centres 2.027E-03 6.081E-03 1.014E-02 1.419E-02 1.824E-02 2.230E-02 2.635E-02 3.041E-02 3.446E-02 3.851E-02 4.257E-02 4.662E-02 5.068E-02 5.473E-02 5.878E-02 6.284E-02 6.689E-02 7.095E-02 7.500E-02 7.905E-02 8.311E-02 8.716E-02 9.122E-02 9.527E-02 9.932E-02 1.034E-01 1.074E-01 1.115E-01 1.155E-01 1.196E-01 1.236E-01 1.277E-01 1.318E-01 1.358E-01 1.399E-01 1.439E-01 1.480E-01 1.520E-01 1.560E-01 1.600E-01 1.640E-01 1.680E-01 1.720E-01 1.760E-01 1.800E-01 1.840E-01 1.880E-01 1.920E-01 1.960E-01 2.000E-01 2.040E-01 2.080E-01 2.120E-01 2.160E-01 2.200E-01 2.240E-01 2.280E-01 2.320E-01 2.360E-01 2.400E-01 2.440E-01 2.480E-01 2.519E-01 2.558E-01 2.597E-01 2.636E-01 2.675E-01 2.714E-01 2.753E-01 2.792E-01 2.831E-01 2.869E-01 2.908E-01 2.947E-01 2.986E-01

3.025E-01 3.064E-01 3.103E-01 3.142E-01 3.181E-01 3.220E-01 3.260E-01 3.300E-01 3.340E-01 3.380E-01 3.420E-01 3.460E-01 3.500E-01 3.540E-01 3.580E-01 3.620E-01 3.660E-01 3.700E-01 3.740E-01 3.780E-01 3.819E-01 3.858E-01 3.897E-01 3.936E-01 3.975E-01 4.014E-01 4.053E-01 4.092E-01 4.131E-01 4.169E-01 4.208E-01 4.247E-01 4.286E-01 4.325E-01 4.364E-01 4.403E-01 4.442E-01 4.481E-01 4.520E-01 4.560E-01 4.600E-01 4.640E-01 4.680E-01 4.720E-01 4.760E-01 4.800E-01 4.840E-01 4.880E-01 4.920E-01 4.960E-01 5.000E-01 5.040E-01 5.080E-01 5.120E-01 5.161E-01 5.202E-01 5.243E-01 5.284E-01 5.325E-01 5.366E-01 5.407E-01 5.448E-01 5.489E-01 5.530E-01 5.570E-01 5.611E-01 5.652E-01 5.693E-01 5.734E-01 5.775E-01 5.816E-01 5.857E-01 5.898E-01 5.939E-01 5.980E-01 Y-coordinates of the cell centres 5.556E-03 1.667E-02 2.778E-02 3.889E-02 5.000E-02 6.111E-02 7.222E-02 8.333E-02 9.444E-02 1.059E-01 1.176E-01 1.294E-01 1.412E-01 1.529E-01 1.647E-01 1.765E-01 1.882E-01 2.000E-01 2.118E-01 2.235E-01 2.353E-01 2.471E-01 2.588E-01 2.706E-01 2.824E-01 2.941E-01 3.059E-01 3.176E-01 3.294E-01 3.412E-01 3.529E-01 3.647E-01 3.765E-01 3.882E-01 4.000E-01 4.118E-01 4.235E-01 4.353E-01 4.471E-01 4.588E-01 4.706E-01 4.824E-01 4.941E-01 5.059E-01 5.176E-01 5.294E-01 5.412E-01 5.529E-01 5.647E-01 5.765E-01 5.882E-01 6.000E-01 6.118E-01 6.235E-01 6.353E-01 6.471E-01 6.588E-01 6.706E-01 6.824E-01 6.941E-01 7.060E-01 7.180E-01 7.300E-01 7.420E-01 7.540E-01 7.660E-01 7.780E-01 7.900E-01 8.020E-01 8.140E-01 8.260E-01 8.380E-01 8.500E-01 8.620E-01 8.740E-01 8.860E-01 8.980E-01 9.100E-01 9.220E-01 9.340E-01 9.460E-01 9.580E-01 9.700E-01 9.820E-01 9.940E-01 1.006E+00 1.018E+00 1.030E+00 1.042E+00 1.054E+00 1.066E+00 1.078E+00 1.090E+00 1.102E+00 1.114E+00 1.126E+00 1.138E+00 1.150E+00 1.162E+00 1.174E+00 1.186E+00 1.198E+00 1.210E+00 1.222E+00 1.234E+00 1.246E+00 1.258E+00 1.270E+00 1.282E+00 1.294E+00

34

1.306E+00 1.318E+00 1.329E+00 1.341E+00 1.353E+00 1.365E+00 1.376E+00 1.388E+00 1.400E+00 1.412E+00 1.424E+00 1.435E+00 1.447E+00 1.459E+00 1.471E+00 1.482E+00 1.494E+00 1.506E+00 1.518E+00 1.529E+00 1.541E+00 1.553E+00 1.565E+00 1.576E+00 1.588E+00 1.600E+00 1.612E+00 1.624E+00 1.635E+00 1.647E+00 1.659E+00 1.671E+00 1.682E+00 1.694E+00 1.706E+00 1.718E+00 1.729E+00 1.741E+00 1.753E+00 1.765E+00 1.776E+00 1.788E+00 1.800E+00 1.812E+00 1.824E+00 1.835E+00 1.847E+00 1.859E+00 1.871E+00 1.882E+00 1.894E+00 1.906E+00 1.917E+00 1.928E+00 1.939E+00 1.950E+00 1.961E+00 1.972E+00 1.983E+00 1.994E+00 Z-coordinates of the cell centres 5.000E-02 X-coordinates of the (higher) cell faces 4.054E-03 8.108E-03 1.216E-02 1.622E-02 2.027E-02 2.432E-02 2.838E-02 3.243E-02 3.649E-02 4.054E-02 4.459E-02 4.865E-02 5.270E-02 5.676E-02 6.081E-02 6.486E-02 6.892E-02 7.297E-02 7.703E-02 8.108E-02 8.514E-02 8.919E-02 9.324E-02 9.730E-02 1.014E-01 1.054E-01 1.095E-01 1.135E-01 1.176E-01 1.216E-01 1.257E-01 1.297E-01 1.338E-01 1.378E-01 1.419E-01 1.459E-01 1.500E-01 1.540E-01 1.580E-01 1.620E-01 1.660E-01 1.700E-01 1.740E-01 1.780E-01 1.820E-01 1.860E-01 1.900E-01 1.940E-01 1.980E-01 2.020E-01 2.060E-01 2.100E-01 2.140E-01 2.180E-01 2.220E-01 2.260E-01 2.300E-01 2.340E-01 2.380E-01 2.420E-01 2.460E-01 2.500E-01 2.539E-01 2.578E-01 2.617E-01 2.656E-01 2.694E-01 2.733E-01 2.772E-01 2.811E-01 2.850E-01 2.889E-01 2.928E-01 2.967E-01 3.006E-01 3.044E-01 3.083E-01 3.122E-01 3.161E-01 3.200E-01 3.240E-01 3.280E-01 3.320E-01 3.360E-01 3.400E-01 3.440E-01 3.480E-01 3.520E-01 3.560E-01 3.600E-01 3.640E-01 3.680E-01 3.720E-01 3.760E-01 3.800E-01 3.839E-01 3.878E-01 3.917E-01 3.956E-01 3.994E-01 4.033E-01 4.072E-01 4.111E-01 4.150E-01 4.189E-01 4.228E-01 4.267E-01 4.306E-01 4.344E-01 4.383E-01 4.422E-01 4.461E-01 4.500E-01 4.540E-01 4.580E-01 4.620E-01 4.660E-01 4.700E-01 4.740E-01 4.780E-01

4.820E-01 4.860E-01 4.900E-01 4.940E-01 4.980E-01 5.020E-01 5.060E-01 5.100E-01 5.141E-01 5.182E-01 5.223E-01 5.264E-01 5.305E-01 5.345E-01 5.386E-01 5.427E-01 5.468E-01 5.509E-01 5.550E-01 5.591E-01 5.632E-01 5.673E-01 5.714E-01 5.755E-01 5.795E-01 5.836E-01 5.877E-01 5.918E-01 5.959E-01 6.000E-01 Y-coordinates of the (higher) cell faces 1.111E-02 2.222E-02 3.333E-02 4.444E-02 5.556E-02 6.667E-02 7.778E-02 8.889E-02 1.000E-01 1.118E-01 1.235E-01 1.353E-01 1.471E-01 1.588E-01 1.706E-01 1.824E-01 1.941E-01 2.059E-01 2.176E-01 2.294E-01 2.412E-01 2.529E-01 2.647E-01 2.765E-01 2.882E-01 3.000E-01 3.118E-01 3.235E-01 3.353E-01 3.471E-01 3.588E-01 3.706E-01 3.824E-01 3.941E-01 4.059E-01 4.176E-01 4.294E-01 4.412E-01 4.529E-01 4.647E-01 4.765E-01 4.882E-01 5.000E-01 5.118E-01 5.235E-01 5.353E-01 5.471E-01 5.588E-01 5.706E-01 5.824E-01 5.941E-01 6.059E-01 6.176E-01 6.294E-01 6.412E-01 6.529E-01 6.647E-01 6.765E-01 6.882E-01 7.000E-01 7.120E-01 7.240E-01 7.360E-01 7.480E-01 7.600E-01 7.720E-01 7.840E-01 7.960E-01 8.080E-01 8.200E-01 8.320E-01 8.440E-01 8.560E-01 8.680E-01 8.800E-01 8.920E-01 9.040E-01 9.160E-01 9.280E-01 9.400E-01 9.520E-01 9.640E-01 9.760E-01 9.880E-01 1.000E+00 1.012E+00 1.024E+00 1.036E+00 1.048E+00 1.060E+00 1.072E+00 1.084E+00 1.096E+00 1.108E+00 1.120E+00 1.132E+00 1.144E+00 1.156E+00 1.168E+00 1.180E+00 1.192E+00 1.204E+00 1.216E+00 1.228E+00 1.240E+00 1.252E+00 1.264E+00 1.276E+00 1.288E+00 1.300E+00 1.312E+00 1.324E+00 1.335E+00 1.347E+00 1.359E+00 1.371E+00 1.382E+00 1.394E+00 1.406E+00 1.418E+00 1.429E+00 1.441E+00 1.453E+00 1.465E+00 1.476E+00 1.488E+00 1.500E+00 1.512E+00 1.524E+00 1.535E+00 1.547E+00 1.559E+00 1.571E+00 1.582E+00 1.594E+00 1.606E+00 1.618E+00 1.629E+00 1.641E+00 1.653E+00 1.665E+00 1.676E+00 1.688E+00 1.700E+00 1.712E+00 1.724E+00 1.735E+00 1.747E+00 1.759E+00 1.771E+00 1.782E+00 1.794E+00 1.806E+00 1.818E+00 1.829E+00

35

1.841E+00 1.853E+00 1.865E+00 1.876E+00 1.888E+00 1.900E+00 1.911E+00 1.922E+00 1.933E+00 1.944E+00 1.956E+00 1.967E+00 1.978E+00 1.989E+00 2.000E+00 Z-coordinates of the (higher) cell faces 1.000E-01 --- INTEGRATION OF EQUATIONS BEGINS --- Flow field at ITHYD= 1, IZ= 1, ISWEEP= 7532, ISTEP= 1 Field Values of P1 IY= 170 4.184E+03 4.184E+03 4.184E+03 4.184E+03 4.184E+03 IY= 136 3.332E+03 3.332E+03 none none none IY= 102 2.465E+03 2.465E+03 none none none IY= 68 1.590E+03 1.590E+03 none none none IY= 34 7.317E+02 7.317E+02 none none none IX= 1 31 61 91 121 Field Values of U1 IY= 170 2.858E-03 1.597E-02 1.064E-02 9.705E-04 -3.988E-05 IY= 136 3.732E-04 2.297E-05 none none none IY= 102 2.272E-05 -1.083E-04 none none none IY= 68 -9.211E-04 2.145E-04 none none none IY= 34 3.803E-05 6.359E-04 none none none IX= 1 31 61 91 121 Field Values of V1 IY= 169 7.833E-03 2.352E-04 -1.028E-03 -6.054E-05 1.240E-04 IY= 135 5.462E-02 1.653E-03 none none none IY= 101 7.734E-02 1.262E-04 none none none IY= 67 3.336E-02 1.431E-03 none none none IY= 33 2.812E-03 1.319E-03 none none none IX= 1 31 61 91 121 Field Values of STAN IY= 170 1.000E-10 1.000E-10 1.000E-10 1.000E-10 1.000E-10 IY= 136 1.000E-10 1.000E-10 none none none IY= 102 1.000E-10 1.000E-10 none none none IY= 68 1.000E-10 1.000E-10 none none none IY= 34 1.000E-10 1.000E-10 none none none IX= 1 31 61 91 121 Field Values of HTCO IY= 170 1.000E-10 1.000E-10 1.000E-10 1.000E-10 1.000E-10 IY= 136 1.000E-10 1.000E-10 none none none IY= 102 1.000E-10 1.000E-10 none none none IY= 68 1.000E-10 1.000E-10 none none none

IY= 34 1.000E-10 1.000E-10 none none none IX= 1 31 61 91 121 Field Values of PRPS IY= 170 pil prop pil prop pil prop pil prop pil prop IY= 136 pil prop pil prop 1.970E+02 1.970E+02 1.970E+02 IY= 102 pil prop pil prop 1.970E+02 1.970E+02 1.970E+02 IY= 68 pil prop pil prop 1.970E+02 1.970E+02 1.970E+02 IY= 34 pil prop pil prop 1.970E+02 1.970E+02 1.970E+02 IX= 1 31 61 91 121 Field Values of TEM1 IY= 170 8.500E+01 8.490E+01 8.489E+01 8.499E+01 8.497E+01 IY= 136 8.500E+01 8.464E+01 8.500E+01 8.500E+01 8.500E+01 IY= 102 8.500E+01 8.344E+01 8.500E+01 8.500E+01 8.500E+01 IY= 68 8.354E+01 8.150E+01 8.500E+01 8.500E+01 8.500E+01 IY= 34 8.082E+01 8.076E+01 8.500E+01 8.500E+01 8.500E+01 IX= 1 31 61 91 121 Field Values of ENUL IY= 170 1.080E-06 1.080E-06 1.080E-06 1.080E-06 1.080E-06 IY= 136 1.080E-06 1.080E-06 none none none IY= 102 1.080E-06 1.080E-06 none none none IY= 68 1.080E-06 1.080E-06 none none none IY= 34 1.080E-06 1.080E-06 none none none IX= 1 31 61 91 121 Field Values of DEN1 IY= 170 6.418E+02 6.418E+02 6.418E+02 6.418E+02 6.418E+02 IY= 136 6.418E+02 6.420E+02 1.000E-06 1.000E-06 1.000E-06 IY= 102 6.418E+02 6.427E+02 1.000E-06 1.000E-06 1.000E-06 IY= 68 6.427E+02 6.439E+02 1.000E-06 1.000E-06 1.000E-06 IY= 34 6.443E+02 6.444E+02 1.000E-06 1.000E-06 1.000E-06 IX= 1 31 61 91 121 Field Values of KOND IY= 170 1.100E-01 1.100E-01 1.100E-01 1.100E-01 1.100E-01 IY= 136 1.100E-01 1.100E-01 0.0 0.0 0.0 IY= 102 1.100E-01 1.100E-01 0.0 0.0 0.0 IY= 68 1.100E-01 1.100E-01 0.0 0.0 0.0 IY= 34 1.100E-01 1.100E-01 0.0 0.0 0.0 IX= 1 31 61 91 121 Field Values of SPH1 IY= 170 2.075E+03 2.075E+03 2.075E+03 2.075E+03 2.075E+03 IY= 136 2.075E+03 2.075E+03 1.000E+20 1.000E+20 1.000E+20 IY= 102 2.075E+03 2.075E+03 1.000E+20 1.000E+20 1.000E+20

36

IY= 68 2.075E+03 2.075E+03 1.000E+20 1.000E+20 1.000E+20 IY= 34 2.075E+03 2.075E+03 1.000E+20 1.000E+20 1.000E+20 IX= 1 31 61 91 121 ************************************************************ Whole-field residuals before solution with resref values determined by EARTH & resfac= 1.000000E-03 variable resref (res sum)/resref (res sum) P1 9.506E-07 3.652E+01 3.471E-05 U1 3.102E-09 1.123E+06 3.484E-03 V1 1.345E-05 8.067E+02 1.085E-02 TEM1 1.624E-01 2.550E+03 4.142E+02 ************************************************************ Sources and sinks Nett source of U1 at patch named: OB5 (CHãO ) = 7.504486E-05 Nett source of U1 at patch named: OB6 (TETO ) =-5.048280E-05 pos. sum= 7.504486E-05 neg. sum=-5.048280E-05 nett sum= 2.456206E-05 Nett source of V1 at patch named: OB3 (ESQUERDA) = 1.548831E-03 Nett source of V1 at patch named: OB4 (ENTRADA ) =-1.244211E-03 Nett source of V1 at patch named: OB7 (DIREITAB) =-7.087897E-05 Nett source of V1 at patch named: OB8 (DIREITAA) =-1.136035E-03 Nett source of V1 at patch named: BUOYANCY = 1.720275E+02 pos. sum= 1.720291E+02 neg. sum=-2.451125E-03 nett sum= 1.720266E+02 pos. sum= 0.000000E+00 neg. sum= 0.000000E+00 nett sum= 0.000000E+00 Nett source of TEM1 at patch named: OB4 (ENTRADA ) = 2.640000E+02 Nett source of TEM1 at patch named: OC1_E (PAREDE1 ) =-6.387757E+01 Nett source of TEM1 at patch named: OC3 (ESQUERDA) =-7.100232E+01 Nett source of TEM1 at patch named: OC9_W (PARADE2 ) =-6.388308E+01 Nett source of TEM1 at patch named: OC9_E (PARADE2 ) =-6.389434E+01 Nett source of TEM1 at patch named: OC10_W (B10 ) =-6.389427E+01 Nett source of TEM1 at patch named: OC10_E (B10 ) =-6.391757E+01 pos. sum= 2.640000E+02 neg. sum=-3.904691E+02 nett sum=-1.264691E+02 ************************************************************ spot values vs sweep or iteration number IXMON= 147 IYMON= 81 IZMON= 1 TIMESTEP= 1 Tabulation of abscissa and ordinates... ISWP P1 U1 V1 TEM1 1 9.641E+02 -1.962E-05 4.685E-05 8.500E+01 501 1.945E+03 1.479E-04 3.787E-03 8.500E+01 1001 1.943E+03 2.024E-04 1.242E-03 8.500E+01 1501 1.940E+03 -1.059E-04 1.090E-02 8.492E+01 2001 1.939E+03 -2.286E-04 5.709E-03 8.473E+01 2501 1.937E+03 1.892E-04 6.473E-03 8.431E+01 3001 1.936E+03 -1.616E-04 2.060E-03 8.397E+01 3501 1.935E+03 -6.858E-05 3.432E-03 8.358E+01 4001 1.934E+03 -2.624E-04 5.717E-03 8.361E+01 4501 1.932E+03 2.339E-04 -1.838E-03 8.346E+01

5001 1.931E+03 -9.814E-05 4.628E-03 8.332E+01 5501 1.930E+03 1.290E-04 3.099E-03 8.324E+01 6001 1.929E+03 7.539E-04 2.472E-04 8.302E+01 6501 1.928E+03 -3.532E-04 3.822E-03 8.273E+01 7001 1.927E+03 -9.755E-06 2.119E-03 8.259E+01 7501 1.927E+03 1.430E-04 1.376E-03 8.244E+01 Variable 1 = P1 2 = U1 3 = V1 4 = TEM1 Minval= 9.641E+02 -3.532E-04 -1.838E-03 8.244E+01 Maxval= 1.945E+03 7.539E-04 1.090E-02 8.500E+01 Cellav= 1.874E+03 3.071E-05 3.301E-03 8.381E+01 1.00 4..4.+.4..3..1.+.1..1..1.+.1..1..1.+.1..2..1.+.1..1 . 4 . 0.90 + 4 + . . 0.80 + + . 4 . 0.70 + + . 3 . 0.60 + 3 4 3 + . 2 . 0.50 + 2 2 3 + . 3 4 4 2 3 2 0.40 + 3 4 3 + . 4 . 0.30 2 3 4 3 + . 3 2 2 4 3 0.20 + 2 + 3 2 3 . 0.10 + 2 2 4 + . 4 . 0.00 1....+....+....+....+....+....3....+....+..2.+....4 0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 the abscissa is ISWP. min= 1.00E+00 max= 7.50E+03 ************************************************************ ************************************************************ residuals vs sweep or iteration number Tabulation of abscissa and ordinates... ISWP P1 U1 V1 TEM1 1 3.974E+04 1.700E+07 2.529E+07 2.311E-01 501 1.018E+01 7.676E+05 2.218E+03 2.437E+03 1001 1.006E+01 7.058E+05 1.805E+03 1.890E+03 1501 1.477E+01 1.135E+06 1.418E+03 2.644E+03 2001 1.535E+01 7.930E+05 1.260E+03 2.346E+03 2501 2.061E+01 1.498E+06 1.137E+03 3.096E+03 3001 2.053E+01 1.136E+06 1.176E+03 2.739E+03 3501 2.434E+01 1.749E+06 1.215E+03 3.045E+03 4001 2.779E+01 1.920E+06 1.047E+03 3.347E+03 4501 2.447E+01 9.768E+05 6.012E+02 3.131E+03 5001 3.174E+01 1.889E+06 8.116E+02 3.504E+03 5501 3.231E+01 1.821E+06 1.006E+03 2.981E+03 6001 3.069E+01 9.919E+05 5.982E+02 2.989E+03 6501 3.678E+01 1.690E+06 8.521E+02 3.219E+03 7001 3.639E+01 1.460E+06 7.394E+02 3.168E+03 7501 3.493E+01 1.001E+06 6.859E+02 2.621E+03 Variable 1 = P1 2 = U1 3 = V1 4 = TEM1

37

Minval= 2.309E+00 1.347E+01 6.394E+00 -1.465E+00 Maxval= 1.059E+01 1.665E+01 1.705E+01 8.162E+00 1.00 3....+....+....+.4..+..4.+.4..4..4.+.4..4..4.+.4..+ . 4 4 4 4 4 4 0.90 + + . . 0.80 + + . . 0.70 + + . . 0.60 + + . . 0.50 + + . . 0.40 + + . . 0.30 + 2 2 2 2 + . 2 2 2 .

0.20 + + . 2 2 1 1 1 1 1 1 0.10 + 3 3 3 1 1 1 1 2 2 2 . 2 1 3 3 3 3 3 3 3 3 . 0.00 4..1.+.2..+....+....+....+....3....+....3....+.3..3 0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 the abscissa is ISWP. min= 1.00E+00 max= 7.50E+03 ************************************************************ ************************************************************ SATLIT RUN NUMBER = 1 ; LIBRARY REF.= 0 Run completed at 14:19:18 on Thursday, 25 June 2009 MACHINE-CLOCK TIME OF RUN = 1408 SECONDS. TIME/(VARIABLES*CELLS*TSTEPS*SWEEPS*ITS) = 1.833E-06 ************************************************************

Documents

Universidade Estadual de Campinas - fem.unicamp.brphoenics/EM974/PROJETOS/PROJETOS%201SEM-09/14/... · Estimativa do Coeficiente de Convecção Natural para o ar ... pelo sistema