Post on 25-Nov-2018
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
Figura 11 – Diagrama de Temperatura.
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
Figura 15 – Diagrama de Temperatura.
O diagrama de temperatura se repetiu sem maiores alterações
Figura 16 – Diagrama de Velocidade.
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.
Figura 20 – Diagrama de Velocidade.
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
No PATCHes used for this Group
************************************************************
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
No PATCHes used for this Group
************************************************************
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
> OBJ, POSITION, 1.216216E-02,
4.444445E-02, 0.000000E+00
> OBJ, SIZE, 4.054054E-03, 1.111111E-
02, 1.000000E-01
> OBJ, GEOMETRY, cubet
> OBJ, ROTATION24, 1
> OBJ, TYPE, PRESSURE_RELIEF
> OBJ, PRES_RELIEF, 1.000000E+04,
1.000000E+01
> OBJ, NAME, ESQUERDA
> OBJ, POSITION, 6.000000E-01,
0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00,
2.000000E+00, 1.000000E-01
29
> OBJ, GEOMETRY, cube13
> OBJ, ROTATION24, 1
> 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, GEOMETRY, cube13
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, SURF_HEVT, 0.000000E+00,
4.400000E+03
> OBJ, NAME, CHãO
> OBJ, POSITION, 0.000000E+00,
0.000000E+00, 0.000000E+00
> OBJ, SIZE, 6.000000E-01,
0.000000E+00, 1.000000E-01
> OBJ, GEOMETRY, cube11
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, NAME, TETO
> OBJ, POSITION, 0.000000E+00,
2.000000E+00, 0.000000E+00
> OBJ, SIZE, 6.000000E-01,
0.000000E+00, 1.000000E-01
> OBJ, GEOMETRY, cube11
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, NAME, DIREITAB
> OBJ, POSITION, 0.000000E+00,
0.000000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 7.000000E-
01, 1.000000E-01
> OBJ, GEOMETRY, cube11
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, NAME, DIREITAA
> OBJ, POSITION, 0.000000E+00,
1.300000E+00, 0.000000E+00
> OBJ, SIZE, 0.000000E+00, 7.000000E-
01, 1.000000E-01
> OBJ, GEOMETRY, cube11
> OBJ, ROTATION24, 1
> OBJ, TYPE, PLATE
> OBJ, NAME, PARADE2
> OBJ, POSITION, 3.200000E-01,
1.000000E-01, 0.000000E+00
> OBJ, SIZE, 6.000000E-02,
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_W, 6.000000E+00,
2.300000E+01
> OBJ, LINR_HEVT_E, 6.000000E+00,
2.300000E+01
> OBJ, NAME, B10
> OBJ, POSITION, 4.500000E-01,
1.000000E-01, 0.000000E+00
> OBJ, SIZE, 6.000000E-02,
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_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: phoenics@cham.co.uk --------------------------------------------------------- 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 ************************************************************