EM 974 – Métodos Computacionais Em Engenharia Térmica …

1

EM 974 – Métodos Computacionais Em Engenharia Térmica e Ambiental

PROFESSOR: Eugênio Spano Rosa

ESCOAMENTO EM BIFURCAÇÕES TIPO “T” Branching Flows

Alunos: RA:

Fabio Arcos Nantes 043284

Marcelo Nogueira Grimberg 045027

Faculdade de Engenharia Mecânica - Universidade Estadual de Campinas

Campinas, SP, 2009.

2

3

ÍNDICE

ESCOAMENTO EM BIFURCAÇÕES TIPO “T” ............................................................................... 1

LISTA DE FIGURAS...................................................................................................................... 4

RESUMO .................................................................................................................................... 6

OBJETIVOS ................................................................................................................................. 6

INTRODUÇÃO ............................................................................................................................ 6

CONHECIMENTO E CURIOSIDADES ........................................................................................... 6

Arteriosclerose ...................................................................................................................... 7

Estenose Traqueal ................................................................................................................. 7

Escoamento petrolífero......................................................................................................... 9

REVISÃO DA LITERATURA ........................................................................................................ 10

ESCOAMENTO RAMIFICADO ................................................................................................... 10

RESULTADOS E MODELAGEM ................................................................................................. 11

Pré- simulação ......................................................................................................................... 11

Modelagem geométrica no PHOENICS ............................................................................... 11

Água ..................................................................................................................................... 12

Conclusão Parcial ................................................................................................................ 16

Simulação final ........................................................................................................................ 16

Modelagem geométrica no Pro-Engenieer ......................................................................... 16

Modelagem geométrica no PHOENICS ............................................................................... 17

Resultados da simulação final ................................................................................................. 22

CONCLUSÃO ............................................................................................................................ 34

REFERÊNCIAS ........................................................................................................................... 35

Web sites: ............................................................................................................................ 36

ANEXOS ................................................................................................................................... 37

Arquivo Result ..................................................................................................................... 37

Arquivo q1 ........................................................................................................................... 44

4

LISTA DE FIGURAS

Figura 1: Tubos usados na medicina.............................................................................................9

Figura 2: Vista da malha-junção tipo T .......................................................................................... 9

Figura 3: Tubo “T” inserido numa traquéia humana..................................................................... 9

Figura 4: Escoamento de petróleo em tubo transversal. .............................................................. 9

Figura 5: Domínio montado. ....................................................................................................... 12

Figura 6: Gráfico de convergência da simulação. ........................................................................ 13

Figura 7: Campo de pressão. ....................................................................................................... 14

Figura 8: Região susceptível à formação de vórtices. ................................................................. 14

Figura 9: Campo de velocidades. ................................................................................................ 15

Figura 10: Linhas de corrente. ..................................................................................................... 15

Figura 11: Domínio. ..................................................................................................................... 16

Figura 12: Tubulação "esculpida". ............................................................................................... 17

Figura 13: Geometria montada. .................................................................................................. 18

Figura 14: INLETS e OUTLET montados no domínio. ................................................................... 18

Figura 15: Perfil do escoamento encontrado por Rosa............................................................... 19

Figura 16: Perfil do escoamento encontrado no PHOENICS. ...................................................... 19

Figura 17: Vistas das malhas montadas no domínio. .................................................................. 20

Figura 18: Número total de células utilizadas. ............................................................................ 20

Figura 19: Gráfico de convergência da simulação....................................................................... 21

Figura 20: Diagrama esquemático dos equipamentos utilizados no trabalho experimental de

Rosa. .................................................................................................................................... 22

Figura 21: Campo de pressão. ..................................................................................................... 22

Figura 22: Campo de velocidades. .............................................................................................. 23

Figura 23: Linhas de corrente e “linha de divisão”. .................................................................... 23

Figura 24: Linhas de corrente e “linha de divisão” para Q=6,64 l/min. ...................................... 25

Figura 25: Linhas de corrente na região de sucção. .................................................................... 25

Figura 26: Vetores de velocidade na região de sucção. .............................................................. 26

5

Figura 27: Diagrama esquemático em três dimensões da superfície de “capturação”. ............. 27

Figura 28: Resultados da superfície de “capturação”. ................................................................ 27

Figura 29: Formação de vórtices contra-rotacionais .................................................................. 28

Figura 30: Região de “capturação” e vórtices contra-rotacionais que a atravessam. ................ 28

Figura 31: Vista de cima dos caminhos percorridos pelos vórtices contra-rotacionais. ............. 29

Figura 32: Vórtices contra-rotacionais ........................................................................................ 30

Figura 33: Velocidades XY no Plano Z25. .................................................................................... 31






Figura 39: Vista de vórtices contra-rotacionais no tubo principal. ............................................. 33

6

RESUMO

A bifurcação em T é presente em diversos sistemas de tubulação; de água, de

petróleo, de gases, sistemas hidráulicos e pneumáticos, equipamentos cirúrgicos e em

muitos outros.

Este trabalho focou o estudo do escoamento em ramificações em tubos tipo “T”.

Para isso, considerou-se o escoamento de água, em regime laminar, analisando sistemas

com diferentes parâmetros.

Além disso, comparou-se resultados de experimentos científicos com o software

de análise PHOENICS, estudando casos de escoamento com água e petróleo.

Entretanto, somente o trabalho de Rosa, o qual utilizou a água como fluido de trabalho,

mostrou-se uma bibliografia adequada e coerente para ser analisada.

OBJETIVOS

Este trabalho tem por objetivo explorar o escoamento ramificado em tubos tipo

“T”. Para isso, foram analisados diferentes fluidos (água e petróleo) escoando em

diferentes regimes de Reynolds. Assim, espera-se observar as diversas linhas de

correntes e a linha de divisão desse escoamento. Em paralelo, será observada e

analisada a variação do campo de pressão, velocidade e temperatura ao longo do tubo

por completo, e da ramificação “T”. Além disso, procurou-se observar os vórtices

contra rotacionais formados no tubo principal perto da ramificação.

INTRODUÇÃO

Escoamento multifásico de fluidos tem sido muito utilizado em várias aplicações

para diversas ocasiões. Indústrias de petróleo, cirurgias traqueais, redes de esgoto,

sistemas pneumáticos, entre outros, são exemplos de áreas em que o escoamento

bifásico de fluidos em conexões tipo T se mostra presente.

O estudo do escoamento neste trabalho focou dois assuntos principais.

Primeiramente, correlacionaram-se trabalhos científicos com os estudos e simulações no

software PHOENICS, e, em seguida, retrataram-se casos ligados à biomedicina e à

escoamento de petróleo.

CONHECIMENTO E CURIOSIDADES

Devido ao enorme interesse de desenvolvimentos de novas técnicas medicinais,

sempre se procura obter inovações que tragam algum benefício ao homem. Assim,

7

analisar um escoamento numa bifurcação tipo “T” é bastante importante quando se

pensa na distribuição arterial, venosa, ou até mesmo traqueal do corpo humano e de

animais (biomedicina).

Arteriosclerose Doenças como a arteriosclerose, termo genérico para espessamento e

endurecimento da parede arterial, é considerada a principal causa de morte no mundo

ocidental.

A aterosclerose, tipo de arteriosclerose, é uma doença que atinge artérias de

grande e médio calibre, como as artérias coronárias, as artérias carótidas e as artérias

dos membros inferiores. É caracterizada pelo depósito de gordura, cálcio e outros

elementos na parede das artérias, reduzindo seu calibre e trazendo um déficit sanguíneo

aos tecidos irrigados por elas. Manifesta-se clinicamente em 10% da população acima

de 50 anos, sendo isso apenas a ponta do iceberg, pois seu desenvolvimento é lento e

progressivo, e é necessário haver uma obstrução arterial significativa, de cerca de 75%

do calibre de uma artéria, para que surjam os primeiros sintomas isquêmicos (sintomas

derivados da falta de sangue).

A superfície interna irregular da artéria com arteriosclerose predispõe à

coagulação sanguínea neste local, com oclusão (entupimento) arterial aguda - trombose

- levando subitamente à falta de sangue para todos os tecidos nutridos por aquela artéria,

que podem entrar em isquemia (sofrimento) ou necrose (morte). Por esta razão, o

primeiro sinal de arteriosclerose pode ser a morte. Com risco tão elevado, é importante

diagnosticar precocemente a doença para detê-la e impedir suas manifestações.

Melhor que tratar é evitar o aparecimento da doença. Isso pode ser alcançado

com uma dieta alimentar equilibrada, não fumando e praticando regularmente exercícios

físicos [14]. Entretanto, caso o pior já esteja presente, a introdução de uma ramificação

em forma de tubo pode ser uma opção biomecânica adequada.

Assim, ajustando-se uma geometria adequada de um tubo ramificado às artérias

humanas, reativa-se a circulação sanguínea e previne-se uma tendência de um grave

quadro clínico.

Estenose Traqueal Estenose traqueal é um estreitamento da luz da traquéia, obstruindo

parcialmente a passagem do ar para os brônquios.

8

Causa irritação constante, tosse, inflamação, dificuldade respiratória e às vezes,

a mudanças nos pulmões e vias aéreas.

A doença ocorre em seres humanos, mas tem uma prevalência maior em animais

de pequeno porte, principalmente cães, é bem menos freqüente em gatos.

A traquéia é a estrutura que leva o ar da laringe até os brônquios. Assemelha-se

a uma mangueira de aspirador de pó – anéis rígidos interligados por um tecido flexível.

O estreitamento pode ter como causa fatores congênitos (fraqueza anormal nos

anéis da traquéia), trauma, necrose, edema, pressão externa devido a tumor ou abcesso,

desenvolvimento abaixo do normal.

A estenose tem como sinais clínicos: tosse forte, dificuldade respiratória,

cansaço mesmo com pouco esforço, em casos graves pneumonia e bronquite [15].

O tratamento é feito através de medicamentos ou cirurgia. Todavia, através do

desenvolvimento medicinal, a utilização de tubos traqueais de bifurcação tipo “T” vem

se tornando muito comum.

O tubo traqueal em T permite manter um calibre adequado na via aérea. Pode

funcionar como suporte da parede traqueal e no tratamento da estenose. Possui a

capacidade de manter esta função por períodos prolongados.

A rama externa dos tubos traqueais T evita o deslocamento e possibilita a

aspiração das secreções bronquiais.

Ambas ramas internas devem ser succionadas. Pode-se utilizar um cateter fino

conectado ao sistema de aspiração. Para guiar o cateter no sentido superior ou inferior

dentro do Stening T, pode-se inclinar a rama externa no sentido oposto ao que se vai

succionar. O mesmo procedimento pode ser efetuado através do canal do

fibrobroncoscópio

Finalizada a aspiração, a rama externa deve permanecer fechada de forma

permanente com a tapa fornecida para este fim. Dessa maneira, é possível a respiração

do ar com a umidade e temperatura resultante da sua passagem através da via área

superior.

O número inscrito indica o diâmetro externo (mm) da rama intra-traqueal [16].

9

Figura 1: Tubos usados na medicina Figura 2: Vista da malha-junção tipo T

Figura 3: Tubo “T” inserido numa traquéia humana.

Escoamento petrolífero Além disso, o escoamento monofásico de fluido em dutos tem sido usado em

várias aplicações industrial tais como transporte de petróleo, refrigeração ambiental,

irrigação, entre outros. Neste sentido, estudar o escoamento permanente de petróleo em

conexões tipo T com uma entrada de fluido e duas saídas torna-se muito relevante.

Figura 4: Escoamento de petróleo em tubo transversal. Pensando em todos esses casos e aplicando-os num software de análise,

consegue-se desenvolver um trabalho de cunho científico, aliando a área biológica à

mecânica industrial. Fora isso, tem-se a possibilidade de inovar mais uma área médico-

mecânica já que nenhuma referência informativa e bibliográfica vai tão a fundo quando

se fala em PHOENICS.

10

REVISÃO DA LITERATURA

Nesta etapa do projeto, estudaram-se artigos disponíveis na literatura.

Revisaram-se conceitos da Mecânica dos Fluidos e estudou-se escoamento ramificado.

Verificou-se a importância do escoamento em dutos com ramificações, pois

provoca alterações no campo de velocidade e na distribuição de pressões ao longo do

tubo e ramificação.

Assim, encontraram-se três artigos relacionados à branch flows:

• O trabalho de doutorado do Professor Eugênio Rosa, com o tema de

escoamento de água num tubo ramificado.

• O trabalho dos alunos da Universidade Federal de Campina Grande, os

quais analisaram o escoamento permanente de petróleo e óleo diesel em

conexões tipo T com uma entrada de fluido e duas saídas.

• O trabalho de Lemonnier e Hervieu (1991), que analisou a distribuição de

pressão em um escoamento monofásico nas vizinhanças da ramificação do

tubo T.

Entretanto, esta gama de literaturas não se mostrou adequada ao nosso projeto, já

que, os trabalhos desenvolvidos pelos alunos da Universidade Federal de Campina

Grande e por Lemonnier e Hervieu não apresentaram resultados coerentes e condizentes

ao estudo em questão.

A revisão da literatura considerada segue abaixo:

ESCOAMENTO RAMIFICADO

Analisando o experimento desenvolvido no trabalho científico do Professor

Eugênio Rosa [3], observou-se o escoamento de água desgaseificada num circuito

fechado. A tubulação tipo “T” analisada consistia num tubo central horizontal de 3

polegadas de diâmetro interno com 180 polegadas de comprimento antecedendo a seção

de teste. Fora isso, havia uma ramificação de diâmetro interno de 0,5 polegada,

perpendicular ao tubo central, onde era succionada uma fração do escoamento do duto

principal.

Assim, o objetivo deste experimento foi estudar a iteração entre os vórtices

rotacionais com os principais fluxos, estabelecer a relação da geometria e parâmetros do

escoamento, como o diâmetro do tubo e a taxa de escoamento. Além disso, estudou-se a

tensão de cisalhamento nas paredes, a distribuição circular ao longo do tubo, e, levou-se

a visualização de um par de vórtices contra-rotacionais.

11

Para a realização deste estudo, considerou-se uma variação do número de

Reynolds entre 500 a 1350, taxas de escoamento de 0.35 a 1.65 l/min, na ramificação e

de 1.90 a 6.65 l/min no tubo principal, e, taxa de variação do diâmetro da ramificação de

0.187 e 0.042.

Com isso, foram encontrados os seguintes resultados relevantes: a média de

valores de vorticidade axial atinge um máximo na ramificação, e, decai

exponencialmente. A distribuição da tensão de cisalhamento, na parede mais próxima

da ramificação (branch) chega de 2 a 4 vezes mais em relação ao aos valores de Poison.

Na parede oposta à ramificação, a distribuição da tensão de cisalhamento atinge um

mínimo na própria ramificação (branch), exibindo valores menores que os de Poison, e,

em alguns casos, indicando separação de escoamento.

RESULTADOS E MODELAGEM

Após um trabalho de análise, decidiu-se apenas comparar a solução

desenvolvida com o trabalho de Rosa. Isto porque, depois de feita a comparação dos

resultados obtidos nos outros trabalhos científicos listados acima, não se encontrou nada

relevante. Fora isso, considerou-se as outras referências duvidosas quanto aos resultados

exibidos e com pouca informação qualitativa.

Então, foi realizada uma análise comparativa com um modelo desenvolvido no

PHOENICS - software que permite ao usuário criar suas próprias rotinas e também

alterar as rotinas abertas do pacote - o qual está sendo estudado pelos autores deste

trabalho na disciplina EM 974 - Métodos Computacionais Em Engenharia Térmica e

Ambiental.

Pré- simulação

Num primeiro momento, focou-se em realizar uma pré-simulação dão problema

em questão. Para isso, considerou-se:

Modelagem geométrica no PHOENICS Através da análise do tutorial número 5, presente no site do PHOENICS (Polis),

montou-se um domínio preenchido por um sólido onde foi “esculpida” uma tubulação

principal e sua ramificação.

12

Primeiramente, para ilustrar o tubo “T”, montou-se uma composição de dois

tubos. Isto é, um tubo maior principal e um menor, na transversal, localizado na região

central.

Analisaram-se dois fluidos com diferentes viscosidades: água, correlacionando

os resultados com o trabalho experimental de doutorado de Rosa [3], e, petróleo,

correlacionando os resultados com o trabalho dos alunos da Unidade Acadêmica de

Engenharia Mecânica da Universidade Federal de Campina Grande [10]. Entretanto,

como comentado acima, só serão expressos os resultados utilizando água.

Água Nesta etapa do trabalho procurou-se repetir a mesmas condições de simulação

experimental feita por Rosa [3].

Para este caso, montou-se o domínio conforme a figura a seguir:

Figura 5: Domínio montado. A malha utilizada para simulação é ilustrada na figura 10. Uma malha

relativamente pouco refinada para que a convergência da simulação não demore muito a

ocorrer.

13

Nesta simulação, usou-se água a 20ºC como fluido de trabalho, numa pressão de

1 atm. Considerou-se um regime turbulento com modelo turbulento LVEL.

Com relação à geometria utilizada, construiu-se uma tubulação principal de

76.2mm de diâmetro e 300mm de comprimento, e, uma ramificação, na região central,

de 12.7mm de diâmetro.

Quanto à formulação do domínio, utilizou-se 1 INLET, na entrada do tubo

principal, e 2 OUTLETS (um na saída do tubo principal e outro na saída da

ramificação).

Com isso, montou-se um cenário com uma vazão volumétrica de 4l/min

escoando através da tubulação, seguindo a variação de 1.9 até 6.65 l/min do trabalho de

Rosa.

A figura 6 ilustra o gráfico de convergência do problema analisado. Neste

cenário, pode-se observar que a convergência se deu com 6% de erro na variável V1.

Figura 6: Gráfico de convergência da simulação.

A Figura 7 retrata o campo de pressão ao longo da tubulação. Nota-se que há

uma relativa variação de pressão na região da bifurcação, significando uma região

susceptível à formação de vórtices.

14

Figura 7: Campo de pressão.

A figura 8 mostra essa região em detalhe ampliado:

Figura 8: Região susceptível à formação de vórtices.

A figura abaixo retrata o campo de velocidades no plano central do tubo.

15

Figura 9: Campo de velocidades. Já a figura 10, mostra as linhas de corrente do escoamento de vazão volumétrica,

na entrada, de 4l/min.

Figura 10: Linhas de corrente.

16

Conclusão Parcial Na simulação com modelagem no PHOENICS, composta de dois tubos

separados, observou-se uma discrepância quanto ao resultado científico revisado.

Entretanto, sabe-se que, uma modelagem mais condizente com a realidade, utilizando,

por exemplo, o software Pro-Engenieer, resultaria em melhores resultados.

Utilizando a água como fluido de trabalho, a fim de retratar um possível

escoamento num caso medicinal e ilustrar o trabalho de Rosa [3], encontrou-se uma

região com diferença de pressão. Quer dizer, na região próxima à bifurcação, notou-se

uma área susceptível à formação de vórtices contra-rotacionais, entretanto, os mesmos

não foram visualizados. Fora isso, observou-se a pertinência do caminho das linhas de

corrente e do campo de velocidade, sendo que, a velocidade de escoamento é maior na

região de bifurcação e praticamente nula na saída, como o esperado.

Simulação final

Neste momento do trabalho, foram corrigidos os possíveis erros presenciados na

pré-simulação, com o intuito de melhorar a solução analisada.

Modelagem geométrica no Pro-Engenieer Para melhorar a simulação deste trabalho, pensou-se em modelar o tubo “T”

utilizando o programa computacional Pro-Engenieer. Assim, para ilustrar o tubo “T”,

montou-se uma única composição, isto é, uma única tubulação composta de um duto

com uma ramificação principal.

A figura a seguir ilustra o procedimento inicial utilizado para a modelagem do

tubo “T”.

Figura 11: Domínio.

A tubulação “esculpida” no domínio é ilustrada na figura abaixo.

17

Figura 12: Tubulação "esculpida".

Modelagem geométrica no PHOENICS Neste momento do trabalho foram determinados os melhores parâmetros a serem

utilizados e melhor forma de desenvolver a simulação para que a mesma parecesse ao

todo com a situação real. Assim, destacou-se a mudança dos seguintes pontos:

• Geometria

Com relação à geometria utilizada, manteve-se a tubulação principal de 76.2mm

de diâmetro e 300 mm de comprimento, e, uma ramificação, na região central, de

12.7mm de diâmetro.

Porém, realizou-se uma mudança significativa quanto à formulação do domínio.

Isto é, utilizam-se agora, dois INLETS e um OUTLET. Quer dizer, um INLET na

entrada do tubo principal e outro na entrada da ramificação, e, um OUTLET na saída do

tubo principal.

As figuras abaixo retratam o comentário acima:

18

Figura 13: Geometria montada.

Figura 14: INLETS e OUTLET montados no domínio.

Com isso, montou-se um cenário com uma vazão volumétrica de 6.64 l/min

escoando através da tubulação, seguindo o modelo de escoamento com perfil uniforme

analisado por Rosa. Este perfil pode ser visto na figura a seguir:

INLETS OUTLET

19

Figura 15: Perfil do escoamento encontrado por Rosa.

A seguir, retrata-se o perfil do escoamento encontrado no PHOENICS:

Figura 16: Perfil do escoamento encontrado no PHOENICS.

Fora isso, na simulação final utilizou-se um INLET com vazão de saída,

negativa, de intensidade iqual à 0,25*Q. Isto é, determinou-se que ¼ da vazão de

entrada do tubo principal se direcionava à saída da ramificação. Esta foi a razão de se

trocar o OUTLET inicialmente considerado (pré-simulação) por um INLET.

• Malha

A malha utilizada na simulação final foi muito mais discretizada e refinada. Isto

é, através da construção de uma malha 3D (três dimensões), aumentou-se o número de

células na região da unificação do branch com o tubo principal, e, em suas

proximidades. Assim, foi possível analisar com maior precisão as regiões críticas do

escoamento, de acordo com o trabalho de Rosa. Para tanto, fez-se uso de um objeto do

Perfil uniforme

20

tipo NULL, este tipo de objeto possibilita a criação de novas áreas para discretização do

modelo.

As figuras abaixo mostram as vistas das malhas e o número de elementos

presentes em cada região, a linha espessas em azul são referentes a geometria do tubo,

as em vermelho são referentes ao objeto do tipo NULL.

Figura 17: Vistas das malhas montadas no domínio.

Figura 18: Número total de células utilizadas.

21

• Convergência

Esse refinamento e detalhamento da malha levaram a uma solução “pesada”.

Quer dizer, a análise de uma região mais discretizada exigiu um maior custo

computacional, fazendo com que a solução convergisse em cerca de três horas com seis

mil iterações. Este processo gerou o seguinte gráfico de convergência, que é referente à

última etapa em que a análise foi rodada:

Figura 19: Gráfico de convergência da simulação.

Pela figura pode-se notar que os erros foram considerados razoavelmente baixos,

com apenas exceção da variável U1, correspondente à velocidade do fluxo no eixo x.

• Regime LVEL para LAMINAR

Na simulação final, usou-se água a 20ºC como fluido de trabalho, numa pressão

de 1 atm. Entretanto, foi corrigido um erro encontrado na pré-simulação, onde se havia

considerado um regime turbulento com modelo turbulento LVEL ao invés de regime

LAMINAR. Assim, com este novo modelo (regime laminar), encontraram-se resultados

bem mais condizentes com a realidade e com o experimento realizado por Rosa.

22

Resultados da simulação final

Para expressar os resultados obtidos no PHOENICS, fez-se uma comparação

visual com as fotos tiradas por Rosa em seus experimentos. Assim, compararam-se,

qualitativamente, os resultados de uma simulação com um software gráfico de alta

precisão, com medidas experimentais práticas.

O experimento de Rosa está esquematizado na figura 20:

Figura 20: Diagrama esquemático dos equipamentos utilizados no trabalho experimental de Rosa.

A Figura 21 retrata o campo de pressão ao longo da tubulação. Nota-se que há

uma relativa variação de pressão na região da bifurcação e a uma distância posterior à

mesma, significando então, uma região susceptível à formação de vórtices.

Figura 21: Campo de pressão.

23

A figura abaixo retrata o campo de velocidades no plano central do tubo.

Figura 22: Campo de velocidades.

Alguns conceitos importantes também foram observados no trabalho de Rosa, como o das linhas de corrente e da “linha de corrente de divisão”, conforme identificados na figura abaixo:

Figura 23: Linhas de corrente e “linha de divisão”.

24

Figura 28: Linha de corrente capturada

Linhas de corrente são linhas contínuas traçadas no líquido que, num mesmo

instante t considerado, mantém-se tangente em todos os pontos à velocidade do

escoamento. Podem também serem definidas como a família de curvas que, para cada

instante de tempo, são as envolventes do campo de velocidades num fluido.

As linhas de corrente são correspondentes diretamente à trajetória da partícula

no fluido.

Já a “linha de corrente de divisão”, correspondente à intersecção da superfície de

“capturação” com o plano simétrico, intercepta a parede do tubo em um ponto nodal. É

notável que esta linha termina a uma distância “n” da ramificação, que corresponde à

menos da metade do diâmetro do branch.

A figura a seguir mostra as linhas de corrente encontradas na solução

desenvolvida no PHOENICS:

25

Figura 24: Linhas de corrente e “linha de divisão” para Q=6,64 l/min. A próxima figura ilustra com mais detalhes a região de divisão das linhas de

correntes que são seccionadas pelo branch:

Figura 25: Linhas de corrente na região de sucção.

Abaixo, os vetores de velocidade na região da ramificaçã:

26

Figura 26: Vetores de velocidade na região de sucção.

Nota-se grande semelhança com a foto do trabalho experimental em questão,

demonstrando que a solução analisada mostra-se dentro do esperado.

Outro ponto analisado foi a presença de vórtices contra-rotacionais na região

após a ramificação, quer dizer, vórtices formados a partir da superfície de “capturação”.

A superfície de “capturação” é uma superfície que divide o escoamento em dois

fluxos: um que continua escoando através do duto principal (com linhas de corrente

tangente à essa superfície) e outro que é sugado para a ramificação. Assim, os múltiplos

vórtices formados atravessam essa superfície em diferentes pontos criando uma

protuberância.

A figura a seguir ilustra o que foi levantado anteriormente:

27

Figura 27: Diagrama esquemático em três dimensões da superfície de “capturação”.

As figuras abaixo mostram a superfície de captura e a geração dos vórtices

contra-rotacionais que a atravessam. Também é possível visualizar o caminho das linhas

de corrente, distinguindo a região em que há sucção, em que há continuidade do fluxo

no tubo principal, em que há retrocesso de fluxo, isto é, formação de vórtices.

Figura 28: Resultados da superfície de “capturação”.

28

Figura 29: Formação de vórtices contra-rotacionais .

Figura 30: Região de “capturação” e vórtices contra-rotacionais que a atravessam.

Uma vista de cima da seção de teste (figura a seguir), com fluxo escoando da

esquerda para a direita, mostra os caminhos que os vórtices percorrem.

29

Figura 31: Vista de cima dos caminhos percorridos pelos vórtices contra-rotacionais.

Já a ilustração seguinte tirada por Rosa, consegue-se notar o semblante dos

vórtices formados através da adição de luz por um laser (He-He). A seqüência de

fotografias são expressas para 1d, 2d, 6d e 9d (d = diâmetro da ramificação) abaixo do

branch.

30

Figura 32: Vórtices contra-rotacionais

As figuras abaixo mostram os resultados alcançados no PHOENCIS utilizando o

Photon, pós-processador. Os vetores destas figuras representam as velocidades no plano

transversal ao tubo, ou seja, desconsideram a componente axial da velocidade. Através

dessas figuras foi possível observar com mais clareza a formação dos vórtices.

As figuras são referentes a região de sucção da ramificação, mas em diferentes

posições. A primeira figura ilustra o campo num plano logo a frente da sucção, no plano

de células Z=25.

31

Figura 33: Velocidades XY no Plano Z25.

A figura 34, referente ao plano Z=35, é a região de sucção propriamente dita.

Nota-se que todos os vetores apontam para o bocal.


Na região logo após a sucção, figura 35, as velocidades obtidas foram

semelhantes à região antes do tubo.

32


As figuras que seguem são das regiões posteriores à sucção, nestas figuras, 36,

37 e 38, é possível observar a formação dos vórtices.



33


Nota-se que os vórtices aumentam e se deslocam para baixo ao decorrer do

escoamento.

A figura 39 mostra os vórtices contra-rotacionais observados no trabalho de

Rosa:

Figura 39: Vista de vórtices contra-rotacionais no tubo principal.

34

CONCLUSÃO

Através da comparação do trabalho experimental de Rosa com os resultados

numéricos obtidos no PHOENICS, conseguiu-se entender inúmeros conceitos a respeito

do escoamento ramificado e retratar um caso prático para um modelo teórico.

Logo pela montagem da solução, pôde-se notar enorme correspondência com o

trabalho científico revisado. A importação de uma figura desenhada no software Pro-

Engenieer resultou num modelo teórico bem coerente com a realidade.

Quanto ao perfil de escoamento, pôde-se notar grande semelhança visual e

característica com o perfil visualizado no trabalho de Rosa.

Na figura das linhas de corrente, é bem notável a correspondência com a teoria.

Isto é, foi identificada a aparição de uma linha de corrente de divisão após a

ramificação, a qual corresponde à intersecção da superfície de “capturação” com o

plano simétrico.

Já as figuras que mostram os vórtices contra-rotacionais que são gerados a partir

da superfície de “capturação”, correlacionam-se com o diagrama teórico apresentado

(figura 27). Pode-se notar a presença dos vórtices após ramificação através de 150

linhas de correntes criadas ao longo de uma pequena região do INLET do tubo

principal.

Já as figuras geradas pelo comando Photon, possibilitaram ver mais claramente

os vórtices, assim como suas geometrias, posição e tamanhos em diferentes planos de

secção do tubo principal.

Assim, conseguiu-se retratar o trabalho de doutorado desenvolvido pelo

Professor Eugênio Rosa através de uma simulação de um problema de branch flow no

software PHOENICS. Os resultados mostraram-se bastante coerentes com a realidade

estudada por Rosa e resultaram em um enorme aprendizado para os alunos criadores

deste trabalho.

35

REFERÊNCIAS

[1] Cat. TM, ANVISA 1021290181

[2] Rodrigo Moura Da Silva*, Severino Rodrigues De Farias Netoº, Siderley Fernandes

Albuquerque*, Antonio Gilson Barbosa De Lima* , Simulação Numérica Do

Escoamento Bifásico De Petróleo E Água Em Conexões Tipo “T” , 8º CONGRESSO

IBEROAMERICANO DE ENGENHARIA MECANICA Cusco, 23 A 25 De Outubro

De 2007

[3] “Wall shear stress and circulation determinations from the measured streamline

pattern in a ninety degree pipe bifurcation”. Doutorado em Mechanical Engineering.

Case Western Reserve University, CWRU, Estados Unidos. Autor: Rosa,Eugênio. 1986.

[4] “Cardiovascular Fluid Mechanics” , editado por Gianni Pedrizzetti e Karl Perktold,

Springer Wien New York

[5] Journal of Applied Physiology Blood flow in branching vessels - Journal of Applied

Physiology - Pinchak and Ostrach.1976; 41: 646-658

[6] “Estudo do escoamento pistonado horizontal ar-água em tubulações com

ramificação "T" “, Tese de doutorado apresentada à comissão de Pós Graduação da

Faculdade de Engenharia Mecânica. Autor: Emerson dos Reis Orientador: Prof.

Leonardo Goldstein Júnior, Campinas, 2003

[7] R. W. Fox and A. T. McDonald, Introdução à Mecânica dos Fluidos, Ed.

Guanabara Dois S. A., Rio de Janeiro,Brasil. 1981.

[8] 2. F. A. Bastos, Problemas de Mecânica dos Fluidos, Editora Santuário, Aparecida.

1983.

[9] CFX 4.4, Solver Manual, CFX International, AEA Technology, UK. 1997.

36

[10] Rodrigo Moura da Silva*, Severino Rodrigues de Farias Netoº, Siderley Fernandes

Albuquerque*, Antonio Gilson Barbosa de Lima*, SIMULAÇÃO DO ESCOAMENTO

MONOFÁSICO DE PETRÓLEO EM JUNÇÕES TIPO “T” VIA MALHA NÃO-

ESTRUTURADA, 8º CONGRESSO IBEROAMERICANO DE ENGENHARIA

MECANICA Cusco, 23 a 25 de Outubro de 2007

[11] R. M. Silva, Simulação Numérica do Escoamento Bifásico Óleo-água em

Bifurcações Angulares, Monografia de Graduação, Programa de Recursos Humanos

para o Setor de Petróleo e Gás, ANP/UFCG-PRH-25, Universidade Federal de Campina

Grande, Campina Grande, Brasil. 77p. 2005.

[12] F. P. Incropera and D. P. DeWitt, Fundamentals of Heat and Mass Transfer, Ed.

John Wiley & Sons, New York, USA. 2002

[13] L. C. Burmeister, Convective Heat Transfer, John Wiley & Sons, New York,

USA. 1983.

Web sites: [14]

http://www.google.com.br/imgres?imgurl=http://www.stening.com.ar/images/productos

/DSC00656.jpg&imgrefurl=http://www.stening.com.ar/PTTML.htm&h=290&w=250&

sz=22&tbnid=Xfc9KQDIJKdFsM::&tbnh=115&tbnw=99&prev=/images%3Fq%3Dtub

o%2BT&hl=pt-BR&usg=__3aQ5wrA-

D8yMhAFbBKv8Uvkrcv8=&ei=rpfBSbfNEor2sAO-

w5DvBg&sa=X&oi=image_result&resnum=2&ct=image&cd=1

[15] http://www.becodosgatos.com.br/estenose.htm

[16] http://www.arteriosclerose.med.br/

37

ANEXOS

Arquivo Result ************************************************************

---------------------------------------------------------

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 24

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

also, other related settings are ...

USEGRX = T

USEGRD = T

>>> End of property-related data <<<

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

Number of F-array locations available is 17353664

Number used before BFC allowance is 14461565

Number used after BFC allowance is 14461565

biggest cell volume divided by average is 18.99432

at:

ix = 1 iy = 60 iz = 14

xg = 2.500000E-03 yg = 1.534798E-01

zg = 1.170000E-01

smallest cell volume divided by average is 4.4135153E-02

at:

ix = 39 iy = 33 iz = 37

xg = 4.733333E-02 yg = 7.852527E-02

zg = 1.516933E-01

ratio of smallest to biggest is 2.3235967E-03

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

Group 1. Run Title and Number

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

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

TEXT(BRANCH FLOW )

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

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

IRUNN = 1 ;LIBREF = 0

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

Group 2. Time dependence

STEADY = T

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

Group 3. X-Direction Grid Spacing

CARTES = T

NX = 67

XULAST = 8.620000E-02

XFRAC ( 1) = 5.800464E-02 ;XFRAC (14) = 2.988631E-01

XFRAC (27) = 4.361562E-01 ;XFRAC (40) = 5.638437E-01

XFRAC (53) = 6.995592E-01 ;XFRAC (66) = 9.419954E-01

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

Group 4. Y-Direction Grid Spacing

NY = 60

YVLAST = 1.562000E-01

YFRAC ( 1) = 3.201024E-02 ;YFRAC (13) = 3.440000E-01

YFRAC (25) = 4.660294E-01 ;YFRAC (37) = 5.250298E-01

YFRAC (49) = 6.168779E-01

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

Group 5. Z-Direction Grid Spacing

PARAB = F

NZ = 81

ZWLAST = 3.000490E-01

ZFRAC ( 1) = 2.888417E-02 ;ZFRAC (17) = 4.408447E-01

ZFRAC (33) = 4.956857E-01 ;ZFRAC (49) = 5.421267E-01

ZFRAC (65) = 5.902300E-01 ;ZFRAC (81) = 1.000000E+00

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

Group 6. Body-Fitted Coordinates

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

Group 7. Variables: STOREd,SOLVEd,NAMEd

ONEPHS = T

NAME( 1) =P1 ;NAME( 3) =U1

NAME( 5) =V1 ;NAME( 7) =W1

NAME(148) =DEN1 ;NAME(150) =PRPS

* 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,Y,N,N,Y)

SOLUTN(U1 ,Y,Y,N,N,N,Y)

SOLUTN(V1 ,Y,Y,N,N,N,Y)

SOLUTN(W1 ,Y,Y,N,N,N,Y)

SOLUTN(DEN1,Y,N,N,N,N,Y)

SOLUTN(PRPS,Y,N,N,N,N,N)

DEN1 = 148

PRPS = 150

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

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 (W1 ,Y,Y,Y,N,Y,N)

DIFCUT = 5.000000E-01 ;ZDIFAC = 1.000000E+00

GALA = F ;ADDDIF = F

ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1

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

Group 9. Properties used if PRPS is not

stored, and where PRPS = -1.0 if it is!

RHO1 = 9.982300E+02 ;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

EMISS = 0.000000E+00 ;SCATT = 0.000000E+00

RADIA = 0.000000E+00 ;RADIB = 0.000000E+00

ENUL = 1.006000E-06 ;ENUT = 0.000000E+00

PRNDTL(U1 ) = 1.000000E+00 ;PRNDTL(V1 ) = 1.000000E+00

PRNDTL(W1 ) = 1.000000E+00

PRT (U1 ) = 1.000000E+00 ;PRT (V1 ) = 1.000000E+00

PRT (W1 ) = 1.000000E+00

CP1 = 4.181800E+03 ;CP2 = 1.000000E+00

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

Group 10.Inter-Phase Transfer Processes

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

Group 11.Initial field variables (PHIs)

FIINIT(P1 ) = READFI ;FIINIT(U1 ) = READFI

FIINIT(V1 ) = READFI ;FIINIT(W1 ) = READFI

FIINIT(DEN1) = READFI ;FIINIT(PRPS) = -1.000000E+00

Parent VR object for this patch is: B1

PATCH(OB1 ,INIVAL, 1, 67, 1, 60, 1, 81, 1,

1)

INIT(OB1 ,PRPS, 0.000000E+00, 1.980000E+02)

INIADD = F

FSWEEP = 1

NAMFI =CHAM

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

Group 12. Patchwise adjustment of terms

38

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: INLET

PATCH(OB2 ,LOW , 1, 67, 1, 37, 1, 1, 1,

1)

COVAL(OB2 ,P1 , FIXFLU , 1.104711E-01)

COVAL(OB2 ,U1 , 0.000000E+00, 0.000000E+00)

COVAL(OB2 ,V1 , 0.000000E+00, 0.000000E+00)

COVAL(OB2 ,W1 , 0.000000E+00, 1.106670E-04)

Parent VR object for this patch is: OUTLET

PATCH(OB3 ,HIGH , 1, 67, 1, 37, 81, 81, 1,

1)

COVAL(OB3 ,P1 , 1.000000E+03, 0.000000E+00)

COVAL(OB3 ,U1 , 0.000000E+00, 0.000000E+00)

COVAL(OB3 ,V1 , 0.000000E+00, 0.000000E+00)

COVAL(OB3 ,W1 , 0.000000E+00, 0.000000E+00)

Parent VR object for this patch is: OUT_RAMI

PATCH(OB4 ,NORTH , 26, 42, 60, 60, 27, 42, 1,

1)

COVAL(OB4 ,P1 ,In-Form:-2.767000E-

05*DEN1!IMAT<100!WHOL)

COVAL(OB4 ,U1 , 0.000000E+00, 0.000000E+00)

COVAL(OB4 ,V1 , 0.000000E+00, 0.000000E+00)

COVAL(OB4 ,W1 , 0.000000E+00, 0.000000E+00)

XCYCLE = F

EGWF = F

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

Group 14. Downstream Pressure For PARAB

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

Group 15. Terminate Sweeps

LSWEEP = 3000 ;ISWC1 = 1

LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1

= 1

SELREF = T

RESFAC = 1.000000E-02

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

Group 16. Terminate Iterations

LITER (P1 ) = 15 ;LITER (U1 ) = 5

LITER (V1 ) = 5 ;LITER (W1 ) = 5

ENDIT (P1 ) = 1.000000E-03 ;ENDIT (U1 ) = 1.000000E-03

ENDIT (V1 ) = 1.000000E-03 ;ENDIT (W1 ) = 1.000000E-03

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

Group 17. Relaxation

RELAX(P1 ,LINRLX, 1.000000E+00)

RELAX(U1 ,FALSDT, 1.000000E+00)

RELAX(V1 ,FALSDT, 1.000000E+00)

RELAX(W1 ,FALSDT, 1.000000E+00)

RELAX(DEN1,LINRLX, 1.000000E+00)

RELAX(PRPS,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+06 ;VARMIN(U1 ) =-1.000000E+06

VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06

VARMAX(W1 ) = 1.000000E+06 ;VARMIN(W1 ) =-1.000000E+06

VARMAX(DEN1) = 1.000000E+10 ;VARMIN(DEN1) =-1.000000E+10

VARMAX(PRPS) = 1.000000E+10 ;VARMIN(PRPS) =-1.000000E+10

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

Group 19. Data transmitted to GROUND

USEGRD = T ;USEGRX = T

ASAP = T

PARSOL = T

SPEDAT(SET,DOMAIN,PHASE_1_MAT,I,67)

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

SPEDAT(SET,MAXINC,U1,R,1.00000E+01)

SPEDAT(SET,MAXINC,V1,R,1.00000E+01)

SPEDAT(SET,MAXINC,W1,R,1.00000E+01)

SPEDAT(SET,B1,DATFILE,C,tubostl)

SPEDAT(SET,OBJNAM,ÔB1,C,B1)

SPEDAT(SET,OBJTYP,ÔB1,C,BLOCKAGE)

SPEDAT(SET,B1,MATERIAL,R,1.98000E+02)

SPEDAT(SET,INLET,DATFILE,C,cylinder)

SPEDAT(SET,OBJNAM,ÔB2,C,INLET)

SPEDAT(SET,OBJTYP,ÔB2,C,INLET)

SPEDAT(SET,OUTLET,DATFILE,C,cylinder)

SPEDAT(SET,OBJNAM,ÔB3,C,OUTLET)

SPEDAT(SET,OBJTYP,ÔB3,C,OUTLET)

SPEDAT(SET,ARATIO,ÔB3,R,1.00000E+00)

SPEDAT(SET,OUT_RAMI,DATFILE,C,cylinder)

SPEDAT(SET,OBJNAM,ÔB4,C,OUT_RAMI)


SPEDAT(SET,SOURCE,R1!OB4,C,=-2.767000E-


SPEDAT(SET,OBJNAM,!OB5,C,B5)

SPEDAT(SET,OBJTYP,!OB5,C,NULL)

SPEDAT(SET,FACETDAT,NUMOBJ,I,5)

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:

* 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,Y,Y,Y,Y)

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

OUTPUT(W1 ,Y,N,Y,Y,Y,Y)

OUTPUT(DEN1,Y,N,Y,N,N,N)

OUTPUT(PRPS,Y,N,Y,N,N,N)

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

Group 22. Monitor Print-Out

IXMON = 67 ;IYMON = 54 ;IZMON = 74

NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = 10001

UWATCH = F ;USTEER = F

HIGHLO = F

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

Group 23.Field Print-Out & Plot Control

NPRINT = 3000 ;NUMCLS = 5

NXPRIN = 13 ;IXPRF = 1 ;IXPRL = 67

NYPRIN = 12 ;IYPRF = 1 ;IYPRL = 60

NZPRIN = 16 ;IZPRF = 1 ;IZPRL = 10000

XZPR = F ;YZPR = F

IPLTF = 1 ;IPLTL = 3000 ;NPLT = 150

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.500E-03 5.901E-03 7.702E-03 9.504E-03 1.131E-

02

1.311E-02 1.491E-02 1.671E-02 1.851E-02 2.031E-

02

2.211E-02 2.347E-02 2.439E-02 2.530E-02 2.622E-

02

2.714E-02 2.805E-02 2.897E-02 2.988E-02 3.080E-

02

3.171E-02 3.263E-02 3.355E-02 3.446E-02 3.538E-

02

3.629E-02 3.717E-02 3.802E-02 3.887E-02 3.971E-

02

4.056E-02 4.141E-02 4.225E-02 4.310E-02 4.395E-

02

4.479E-02 4.564E-02 4.649E-02 4.733E-02 4.818E-

02

4.903E-02 4.990E-02 5.081E-02 5.171E-02 5.262E-

02

5.352E-02 5.442E-02 5.533E-02 5.623E-02 5.714E-

02

5.804E-02 5.895E-02 5.985E-02 6.075E-02 6.166E-

02

6.256E-02 6.392E-02 6.574E-02 6.756E-02 6.938E-

02

7.120E-02 7.302E-02 7.484E-02 7.665E-02 7.847E-

02

8.029E-02 8.370E-02

Y-coordinates of the cell centres

2.500E-03 7.031E-03 1.109E-02 1.515E-02 1.921E-

02

2.327E-02 2.734E-02 3.140E-02 3.546E-02 3.952E-

02

4.358E-02 4.764E-02 5.170E-02 5.576E-02 5.982E-

02

6.389E-02 6.630E-02 6.706E-02 6.783E-02 6.859E-

02

6.935E-02 7.012E-02 7.088E-02 7.165E-02 7.241E-

02

7.318E-02 7.394E-02 7.470E-02 7.547E-02 7.623E-

02

39

7.700E-02 7.776E-02 7.853E-02 7.929E-02 8.005E-

02

8.082E-02 8.160E-02 8.241E-02 8.322E-02 8.403E-

02

8.484E-02 8.565E-02 8.646E-02 8.727E-02 8.808E-

02

8.889E-02 8.970E-02 9.051E-02 9.364E-02 9.908E-

02

1.045E-01 1.100E-01 1.154E-01 1.208E-01 1.263E-

01

1.317E-01 1.372E-01 1.426E-01 1.480E-01 1.535E-

01

Z-coordinates of the cell centres

4.333E-03 1.300E-02 2.167E-02 3.033E-02 3.900E-

02

4.767E-02 5.633E-02 6.500E-02 7.367E-02 8.233E-

02

9.100E-02 9.967E-02 1.083E-01 1.170E-01 1.257E-

01

1.306E-01 1.317E-01 1.328E-01 1.340E-01 1.351E-

01

1.363E-01 1.374E-01 1.385E-01 1.397E-01 1.408E-

01

1.419E-01 1.431E-01 1.441E-01 1.449E-01 1.458E-

01

1.466E-01 1.475E-01 1.483E-01 1.492E-01 1.500E-

01

1.508E-01 1.517E-01 1.525E-01 1.534E-01 1.542E-

01

1.551E-01 1.559E-01 1.568E-01 1.577E-01 1.586E-

01

1.595E-01 1.604E-01 1.613E-01 1.622E-01 1.631E-

01

1.640E-01 1.649E-01 1.658E-01 1.667E-01 1.676E-

01

1.685E-01 1.694E-01 1.703E-01 1.712E-01 1.721E-

01

1.730E-01 1.739E-01 1.748E-01 1.757E-01 1.766E-

01

1.775E-01 1.821E-01 1.902E-01 1.983E-01 2.065E-

01

2.146E-01 2.228E-01 2.309E-01 2.390E-01 2.472E-

01

2.553E-01 2.634E-01 2.716E-01 2.797E-01 2.878E-

01

2.960E-01

X-coordinates of the (higher) cell faces

5.000E-03 6.801E-03 8.603E-03 1.040E-02 1.221E-

02

1.401E-02 1.581E-02 1.761E-02 1.941E-02 2.121E-

02

2.301E-02 2.393E-02 2.485E-02 2.576E-02 2.668E-

02

2.759E-02 2.851E-02 2.942E-02 3.034E-02 3.126E-

02

3.217E-02 3.309E-02 3.400E-02 3.492E-02 3.583E-

02

3.675E-02 3.760E-02 3.844E-02 3.929E-02 4.014E-

02

4.098E-02 4.183E-02 4.268E-02 4.352E-02 4.437E-

02

4.522E-02 4.606E-02 4.691E-02 4.776E-02 4.860E-

02

4.945E-02 5.035E-02 5.126E-02 5.216E-02 5.307E-

02

5.397E-02 5.488E-02 5.578E-02 5.668E-02 5.759E-

02

5.849E-02 5.940E-02 6.030E-02 6.121E-02 6.211E-

02

6.301E-02 6.483E-02 6.665E-02 6.847E-02 7.029E-

02

7.211E-02 7.393E-02 7.574E-02 7.756E-02 7.938E-

02

8.120E-02 8.620E-02

Y-coordinates of the (higher) cell faces

5.000E-03 9.061E-03 1.312E-02 1.718E-02 2.124E-

02

2.531E-02 2.937E-02 3.343E-02 3.749E-02 4.155E-

02

4.561E-02 4.967E-02 5.373E-02 5.779E-02 6.185E-

02

6.592E-02 6.668E-02 6.744E-02 6.821E-02 6.897E-

02

6.974E-02 7.050E-02 7.127E-02 7.203E-02 7.279E-

02

7.356E-02 7.432E-02 7.509E-02 7.585E-02 7.661E-

02

7.738E-02 7.814E-02 7.891E-02 7.967E-02 8.044E-

02

8.120E-02 8.201E-02 8.282E-02 8.363E-02 8.444E-

02

8.525E-02 8.606E-02 8.687E-02 8.768E-02 8.849E-

02

8.930E-02 9.011E-02 9.092E-02 9.636E-02 1.018E-

01

1.072E-01 1.127E-01 1.181E-01 1.236E-01 1.290E-

01

1.344E-01 1.399E-01 1.453E-01 1.508E-01 1.562E-

01

Z-coordinates of the (higher) cell faces

8.667E-03 1.733E-02 2.600E-02 3.467E-02 4.333E-

02

5.200E-02 6.067E-02 6.933E-02 7.800E-02 8.667E-

02

9.533E-02 1.040E-01 1.127E-01 1.213E-01 1.300E-

01

1.311E-01 1.323E-01 1.334E-01 1.345E-01 1.357E-

01

1.368E-01 1.380E-01 1.391E-01 1.402E-01 1.414E-

01

1.425E-01 1.436E-01 1.445E-01 1.453E-01 1.462E-

01

1.470E-01 1.479E-01 1.487E-01 1.496E-01 1.504E-

01

1.513E-01 1.521E-01 1.530E-01 1.538E-01 1.547E-

01

1.555E-01 1.564E-01 1.573E-01 1.582E-01 1.591E-

01

1.600E-01 1.609E-01 1.618E-01 1.627E-01 1.636E-

01

1.645E-01 1.654E-01 1.663E-01 1.672E-01 1.681E-

01

1.690E-01 1.699E-01 1.708E-01 1.717E-01 1.726E-

01

1.735E-01 1.744E-01 1.753E-01 1.762E-01 1.771E-

01

1.780E-01 1.861E-01 1.943E-01 2.024E-01 2.105E-

01

2.187E-01 2.268E-01 2.350E-01 2.431E-01 2.512E-

01

2.594E-01 2.675E-01 2.756E-01 2.838E-01 2.919E-

01

3.000E-01

--- INTEGRATION OF EQUATIONS BEGINS ---

Flow field at ITHYD= 1, IZ= 16, ISWEEP= 3000, ISTEP=

1

Field Values of P1

IY= 60 none none none none

none


none

IY= 36 none none -4.647E-01 -4.813E-01

none

IY= 24 none -1.824E-01 -3.376E-01 -3.468E-01 -

1.940E-01

IY= 12 none -6.417E-02 -7.608E-02 -7.657E-02 -

6.551E-02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 -4.744E-02

IX= 66

Field Values of U1


none


none

IY= 36 none none 1.534E-03 -1.543E-03

none

IY= 24 none 3.709E-03 1.902E-03 -1.904E-03 -

3.685E-03

IY= 12 none 1.021E-03 4.086E-04 -4.055E-04 -

1.012E-03

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 none

IX= 66

Field Values of V1


none


none

IY= 35 none none 4.835E-04 5.349E-04

none

40

IY= 23 none 3.300E-03 4.490E-03 4.578E-03

3.285E-03

IY= 11 none 2.187E-03 2.701E-03 2.722E-03

2.249E-03

IX= 1 14 27 40

53

IY= 59 none

IY= 47 none

IY= 35 none

IY= 23 none

IY= 11 1.160E-03

IX= 66

Field Values of W1


none


none

IY= 36 none none 6.142E-03 9.050E-03

none

IY= 24 none 2.698E-02 3.651E-02 3.680E-02

2.883E-02

IY= 12 none 2.761E-02 2.775E-02 2.776E-02

2.762E-02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 1.133E-03

IX= 66

Field Values of DEN1


none


none

IY= 36 none none 9.982E+02 9.982E+02

none

IY= 24 none 9.982E+02 9.982E+02 9.982E+02

9.982E+02

IY= 12 none 9.982E+02 9.982E+02 9.982E+02

9.982E+02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 9.982E+02

IX= 66

Field Values of PRPS

IY= 60 blockage blockage blockage blockage

blockage


blockage

IY= 36 blockage blockage pil prop pil prop

blockage

IY= 24 blockage pil prop pil prop pil prop

pil prop


pil prop

IX= 1 14 27 40

53

IY= 60 blockage

IY= 48 blockage

IY= 36 blockage

IY= 24 blockage

IY= 12 pil prop

IX= 66


1

Field Values of P1

IY= 60 none none -4.296E+01 -4.310E+01

none

IY= 48 none none -3.715E+01 -3.714E+01

none

IY= 36 none none -3.431E+01 -3.105E+01

none

IY= 24 none -1.537E-01 -7.606E-01 -8.201E-01 -

1.764E-01

IY= 12 none -3.688E-02 -4.547E-02 -4.583E-02 -

3.782E-02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 -2.435E-02

IX= 66

Field Values of U1

IY= 60 none none 2.056E-03 -2.087E-03

none

IY= 48 none none 6.552E-03 -6.584E-03

none

IY= 36 none none 1.309E-01 -1.309E-01

none

IY= 24 none 1.128E-02 1.324E-02 -1.323E-02 -

1.134E-02

IY= 12 none 1.665E-03 7.086E-04 -7.046E-04 -

1.652E-03

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 none

IX= 66

Field Values of V1

IY= 59 none none 4.178E-02 1.704E-01

none

IY= 47 none none 4.925E-02 1.572E-01

none

IY= 35 none none 1.390E-01 1.759E-01

none

IY= 23 none 8.386E-03 2.743E-02 2.940E-02

8.951E-03

IY= 11 none 3.054E-03 4.094E-03 4.139E-03

3.172E-03

IX= 1 14 27 40

53

IY= 59 none

IY= 47 none

IY= 35 none

IY= 23 none

IY= 11 8.269E-04

IX= 66

Field Values of W1

IY= 60 none none -2.120E-02 -2.805E-02

none

IY= 48 none none -1.314E-02 -1.224E-02

none

IY= 36 none none 7.791E-02 5.796E-02

none

IY= 24 none 2.044E-02 3.220E-02 3.252E-02

2.262E-02

IY= 12 none 2.647E-02 2.636E-02 2.636E-02

2.645E-02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 6.485E-05

IX= 66


IY= 60 none none 9.982E+02 9.982E+02

none

IY= 48 none none 9.982E+02 9.982E+02

none

IY= 36 none none 9.982E+02 9.982E+02

none

IY= 24 none 9.982E+02 9.982E+02 9.982E+02

9.982E+02

IY= 12 none 9.982E+02 9.982E+02 9.982E+02

9.982E+02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 9.982E+02

IX= 66



blockage


blockage


blockage


pil prop


pil prop

IX= 1 14 27 40

53

IY= 60 blockage

41

IY= 48 blockage

IY= 36 blockage

IY= 24 blockage

IY= 12 pil prop

IX= 66


1

Field Values of P1


none


none

IY= 36 none none 1.673E-01 1.783E-01

none

IY= 24 none 2.423E-02 7.608E-02 8.046E-02

2.590E-02

IY= 12 none -4.465E-03 -1.498E-03 -1.347E-03 -

4.135E-03

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 -6.684E-03

IX= 66

Field Values of U1


none


none

IY= 36 none none -3.695E-04 3.953E-04

none

IY= 24 none 7.446E-03 6.450E-03 -6.452E-03 -

7.528E-03

IY= 12 none 1.661E-03 6.914E-04 -6.884E-04 -

1.648E-03

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 none

IX= 66

Field Values of V1


none


none

IY= 35 none none 6.939E-04 1.747E-03

none

IY= 23 none 5.050E-03 1.509E-02 1.581E-02

5.765E-03

IY= 11 none 2.946E-03 3.917E-03 3.958E-03

3.058E-03

IX= 1 14 27 40

53

IY= 59 none

IY= 47 none

IY= 35 none

IY= 23 none

IY= 11 2.247E-04

IX= 66

Field Values of W1


none


none

IY= 36 none none -4.031E-03 -6.329E-03

none

IY= 24 none 1.109E-02 1.455E-02 1.398E-02

1.286E-02

IY= 12 none 2.528E-02 2.471E-02 2.468E-02

2.521E-02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 -3.880E-04

IX= 66



none


none

IY= 36 none none 9.982E+02 9.982E+02

none

IY= 24 none 9.982E+02 9.982E+02 9.982E+02

9.982E+02

IY= 12 none 9.982E+02 9.982E+02 9.982E+02

9.982E+02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 9.982E+02

IX= 66



blockage


blockage


blockage


pil prop


pil prop

IX= 1 14 27 40

53

IY= 60 blockage

IY= 48 blockage

IY= 36 blockage

IY= 24 blockage

IY= 12 pil prop

IX= 66


1

Field Values of P1


none


none

IY= 36 none none 1.630E-01 1.667E-01

none

IY= 24 none 4.267E-02 1.074E-01 1.103E-01

4.580E-02

IY= 12 none 1.494E-02 2.129E-02 2.158E-02

1.566E-02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 6.721E-03

IX= 66

Field Values of U1


none


none

IY= 36 none none -8.879E-04 9.043E-04

none

IY= 24 none 1.679E-03 1.027E-03 -1.027E-03 -

1.816E-03

IY= 12 none 1.258E-03 4.842E-04 -4.836E-04 -

1.249E-03

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 none

IX= 66

Field Values of V1


none


none

IY= 35 none none -2.185E-04 -2.335E-04

none

IY= 23 none 1.697E-03 5.112E-03 5.200E-03

2.646E-03

IY= 11 none 2.328E-03 2.898E-03 2.920E-03

2.398E-03

IX= 1 14 27 40

53

IY= 59 none

IY= 47 none

IY= 35 none

IY= 23 none

42

IY= 11 -2.362E-04

IX= 66

Field Values of W1


none


none

IY= 36 none none 2.402E-03 3.522E-03

none

IY= 24 none 1.116E-02 2.078E-02 2.073E-02

1.342E-02

IY= 12 none 2.465E-02 2.400E-02 2.397E-02

2.458E-02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 -2.336E-04

IX= 66



none


none

IY= 36 none none 9.982E+02 9.982E+02

none

IY= 24 none 9.982E+02 9.982E+02 9.982E+02

9.982E+02

IY= 12 none 9.982E+02 9.982E+02 9.982E+02

9.982E+02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 9.982E+02

IX= 66



blockage


blockage


blockage


pil prop


pil prop

IX= 1 14 27 40

53

IY= 60 blockage

IY= 48 blockage

IY= 36 blockage

IY= 24 blockage

IY= 12 pil prop

IX= 66


1

Field Values of P1


none


none

IY= 36 none none 7.313E-03 7.380E-03

none

IY= 24 none 9.723E-03 1.392E-02 1.407E-02

1.015E-02

IY= 12 none 1.540E-02 1.860E-02 1.868E-02

1.590E-02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 5.167E-03

IX= 66

Field Values of U1


none


none

IY= 36 none none -4.441E-05 4.430E-05

none

IY= 24 none -1.785E-03 -3.013E-04 3.084E-04

1.796E-03

IY= 12 none -1.552E-04 -3.742E-05 3.863E-05

1.551E-04

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 none

IX= 66

Field Values of V1


none


none

IY= 35 none none 1.261E-05 2.466E-05

none

IY= 23 none -9.219E-05 8.975E-04 9.186E-04

1.082E-04

IY= 11 none 6.177E-04 5.731E-04 5.718E-04

6.182E-04

IX= 1 14 27 40

53

IY= 59 none

IY= 47 none

IY= 35 none

IY= 23 none

IY= 11 -3.868E-04

IX= 66

Field Values of W1


none


none

IY= 36 none none 1.705E-03 2.193E-03

none

IY= 24 none 2.125E-02 2.432E-02 2.432E-02

2.213E-02

IY= 12 none 2.416E-02 2.409E-02 2.407E-02

2.424E-02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 2.344E-03

IX= 66



none


none

IY= 36 none none 9.982E+02 9.982E+02

none

IY= 24 none 9.982E+02 9.982E+02 9.982E+02

9.982E+02

IY= 12 none 9.982E+02 9.982E+02 9.982E+02

9.982E+02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 9.982E+02

IX= 66



blockage


blockage


blockage


pil prop


pil prop

IX= 1 14 27 40

53

IY= 60 blockage

IY= 48 blockage

IY= 36 blockage

IY= 24 blockage

IY= 12 pil prop

IX= 66

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

Whole-field residuals before solution

with resref values determined by EARTH

43

& resfac= 1.000000E-02

variable resref (res sum)/resref (res sum)

P1 3.134E-07 8.825E-03 2.766E-09

U1 1.115E-09 1.936E+01 2.158E-08

V1 7.576E-09 2.817E+00 2.134E-08

W1 7.530E-09 3.655E+00 2.752E-08

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

Sources and sinks

Nett source of U1 at patch named: OB3 (OUTLET ) =-

3.980025E-08

Nett source of U1 at patch named: OB4 (OUT_RAMI) =

7.931286E-07

pos. sum= 7.931286E-07 neg. sum=-3.980025E-08

nett sum= 7.533284E-07

Nett source of V1 at patch named: OB3 (OUTLET ) =-

1.177284E-05

Nett source of V1 at patch named: OB4 (OUT_RAMI) =-

6.094207E-03

pos. sum= 0.000000E+00 neg. sum=-6.105980E-03

nett sum=-6.105980E-03

Nett source of W1 at patch named: OB2 (INLET ) =

2.685783E-03

Nett source of W1 at patch named: OB3 (OUTLET ) =-

1.687669E-03

Nett source of W1 at patch named: OB4 (OUT_RAMI) =

2.980590E-04

pos. sum= 2.983843E-03 neg. sum=-1.687669E-03

nett sum= 1.296173E-03

Nett source of R1 at patch named: OB2 (INLET ) =

1.104712E-01

Nett source of R1 at patch named: OB3 (OUTLET ) =-

8.285094E-02

Nett source of R1 at patch named: OB4 (OUT_RAMI) =-

2.762031E-02

pos. sum= 1.104712E-01 neg. sum=-1.104712E-01

nett sum= 0.000000E+00

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

spot values vs sweep or iteration number

IXMON= 67 IYMON= 54 IZMON= 74

TIMESTEP= 1

Tabulation of abscissa and ordinates...

ISWP P1 U1 V1 W1

1 0.000E+00 0.000E+00 0.000E+00 0.000E+00

151 0.000E+00 0.000E+00 0.000E+00 0.000E+00

301 0.000E+00 0.000E+00 0.000E+00 0.000E+00

451 0.000E+00 0.000E+00 0.000E+00 0.000E+00

601 0.000E+00 0.000E+00 0.000E+00 0.000E+00

751 0.000E+00 0.000E+00 0.000E+00 0.000E+00

901 0.000E+00 0.000E+00 0.000E+00 0.000E+00

1051 0.000E+00 0.000E+00 0.000E+00 0.000E+00

1201 0.000E+00 0.000E+00 0.000E+00 0.000E+00

1351 0.000E+00 0.000E+00 0.000E+00 0.000E+00

1501 0.000E+00 0.000E+00 0.000E+00 0.000E+00

1651 0.000E+00 0.000E+00 0.000E+00 0.000E+00

1801 0.000E+00 0.000E+00 0.000E+00 0.000E+00

1951 0.000E+00 0.000E+00 0.000E+00 0.000E+00

2101 0.000E+00 0.000E+00 0.000E+00 0.000E+00

2251 0.000E+00 0.000E+00 0.000E+00 0.000E+00

2401 0.000E+00 0.000E+00 0.000E+00 0.000E+00

2551 0.000E+00 0.000E+00 0.000E+00 0.000E+00

2701 0.000E+00 0.000E+00 0.000E+00 0.000E+00

2851 0.000E+00 0.000E+00 0.000E+00 0.000E+00

Variable 1 = P1 2 = U1 3 = V1 4 = W1

Minval= 0.000E+00 0.000E+00 0.000E+00 0.000E+00

Maxval= 0.000E+00 0.000E+00 0.000E+00 0.000E+00

Cellav= 0.000E+00 0.000E+00 0.000E+00 0.000E+00

1.00 +....+....+....+....+....+....+....+....+....+....+

. .

0.90 + +

. .

0.80 + +

. .

0.70 + +

. .

0.60 + +

. .

0.50 + +

. .

0.40 + +

. .

0.30 + +

. .

0.20 + +

. .

0.10 + +

. .

0.00 4..4.4..4.+4.4.+4.4.+4..4+4..4+.4.4+.4.4+.4..4.4..4

0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0

the abscissa is ISWP. min= 1.00E+00 max= 2.85E+03

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

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

residuals vs sweep or iteration number


ISWP P1 U1 V1 W1

1 1.241E+01 3.377E+01 6.483E+01 2.183E+04

151 8.122E-02 6.263E+01 1.213E+01 5.698E+01

301 1.025E-01 3.043E+01 6.270E+00 1.600E+01

451 5.695E-02 3.319E+01 6.041E+00 3.967E+01

601 7.341E-02 2.562E+01 4.795E+00 1.117E+01

751 4.150E-02 2.870E+01 4.937E+00 2.876E+01

901 5.337E-02 2.336E+01 4.020E+00 8.143E+00

1051 3.062E-02 2.505E+01 4.196E+00 2.107E+01

1201 3.854E-02 2.144E+01 3.438E+00 6.372E+00

1351 2.278E-02 2.265E+01 3.626E+00 1.548E+01

1501 2.876E-02 2.087E+01 3.183E+00 5.332E+00

1651 1.759E-02 2.160E+01 3.269E+00 1.176E+01

1801 2.148E-02 2.033E+01 2.956E+00 4.739E+00

1951 1.395E-02 2.061E+01 3.081E+00 8.923E+00

2101 1.675E-02 1.968E+01 2.941E+00 4.344E+00

2251 1.129E-02 2.010E+01 2.913E+00 7.039E+00

2401 1.292E-02 1.959E+01 2.856E+00 4.007E+00

2551 9.578E-03 2.013E+01 2.898E+00 5.645E+00

2701 1.074E-02 1.928E+01 2.772E+00 3.888E+00

2851 8.183E-03 1.953E+01 2.817E+00 4.837E+00

Variable 1 = P1 2 = U1 3 = V1 4 = W1

Minval= -4.806E+00 2.959E+00 1.020E+00 1.358E+00

Maxval= 2.518E+00 4.137E+00 4.172E+00 9.991E+00

1.00 4..2.+....+....+....+....+....+....+....+....+....+

. .

0.90 + +

. .

0.80 + +

. .

0.70 + +

. .

0.60 + +

. .

0.50 2 +

. 3 2 .

0.40 + 2 +

. 1 2 .

0.30 + 4 1 +

. 3 4 2 4 1 .

0.20 + 3 4 1 +

. 4 3 2 3 4 1 4 1 .

0.10 + 4 4 2 3 2 4 1 +

. 4 4 3 2 3 4 1 4 1 4

0.00 +....+....+....+....+....+....+.4..+.4.3+.4..3.4..3

0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0


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

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

SATLIT RUN NUMBER = 1 ; LIBRARY REF.= 0

Run completed at 19:20:08 on Tuesday, 16 June 2009

MACHINE-CLOCK TIME OF RUN = 5486 SECONDS.

TIME/(VARIABLES*CELLS*TSTEPS*SWEEPS*ITS) = 5.486E+03

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

44

Arquivo q1 *********************************************************

***

------------------------------------------------------

---

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 24

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

also, other related settings are ...

USEGRX = T

USEGRD = T

>>> End of property-related data <<<

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

***

Number of F-array locations available is 17353664

Number used before BFC allowance is 14461565

Number used after BFC allowance is 14461565

biggest cell volume divided by average is 18.99432

at:

ix = 1 iy = 60 iz =

14

xg = 2.500000E-03 yg = 1.534798E-01

zg = 1.170000E-01

smallest cell volume divided by average is 4.4135153E-

02 at:

ix = 39 iy = 33 iz =

37

xg = 4.733333E-02 yg = 7.852527E-02

zg = 1.516933E-01

ratio of smallest to biggest is 2.3235967E-03

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

***

Group 1. Run Title and Number

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

***

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

***

TEXT(BRANCH FLOW )

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

***

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

***

IRUNN = 1 ;LIBREF = 0

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

***

Group 2. Time dependence

STEADY = T

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

***

Group 3. X-Direction Grid Spacing

CARTES = T

NX = 67

XULAST = 8.620000E-02

XFRAC ( 1) = 5.800464E-02 ;XFRAC (14) = 2.988631E-01

XFRAC (27) = 4.361562E-01 ;XFRAC (40) = 5.638437E-01

XFRAC (53) = 6.995592E-01 ;XFRAC (66) = 9.419954E-01

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

***

Group 4. Y-Direction Grid Spacing

NY = 60

YVLAST = 1.562000E-01

YFRAC ( 1) = 3.201024E-02 ;YFRAC (13) = 3.440000E-01

YFRAC (25) = 4.660294E-01 ;YFRAC (37) = 5.250298E-01

YFRAC (49) = 6.168779E-01

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

***

Group 5. Z-Direction Grid Spacing

PARAB = F

NZ = 81

ZWLAST = 3.000490E-01

ZFRAC ( 1) = 2.888417E-02 ;ZFRAC (17) = 4.408447E-01

ZFRAC (33) = 4.956857E-01 ;ZFRAC (49) = 5.421267E-01

ZFRAC (65) = 5.902300E-01 ;ZFRAC (81) = 1.000000E+00

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

***

Group 6. Body-Fitted Coordinates

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

***

Group 7. Variables: STOREd,SOLVEd,NAMEd

ONEPHS = T

NAME( 1) =P1 ;NAME( 3) =U1

NAME( 5) =V1 ;NAME( 7) =W1

NAME(148) =DEN1 ;NAME(150) =PRPS

* 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,Y,N,N,Y)

SOLUTN(U1 ,Y,Y,N,N,N,Y)

SOLUTN(V1 ,Y,Y,N,N,N,Y)

SOLUTN(W1 ,Y,Y,N,N,N,Y)

SOLUTN(DEN1,Y,N,N,N,N,Y)

SOLUTN(PRPS,Y,N,N,N,N,N)

DEN1 = 148

PRPS = 150

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

***

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 (W1 ,Y,Y,Y,N,Y,N)

DIFCUT = 5.000000E-01 ;ZDIFAC = 1.000000E+00

GALA = F ;ADDDIF = F

ISOLX = -1 ;ISOLY = -1 ;ISOLZ =

-1

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

***

Group 9. Properties used if PRPS is not

stored, and where PRPS = -1.0 if it is!

RHO1 = 9.982300E+02 ;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

EMISS = 0.000000E+00 ;SCATT = 0.000000E+00

RADIA = 0.000000E+00 ;RADIB = 0.000000E+00

ENUL = 1.006000E-06 ;ENUT = 0.000000E+00

PRNDTL(U1 ) = 1.000000E+00 ;PRNDTL(V1 ) =

1.000000E+00

PRNDTL(W1 ) = 1.000000E+00

45

PRT (U1 ) = 1.000000E+00 ;PRT (V1 ) =

1.000000E+00

PRT (W1 ) = 1.000000E+00

CP1 = 4.181800E+03 ;CP2 = 1.000000E+00

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

***

Group 10.Inter-Phase Transfer Processes

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

***

Group 11.Initial field variables (PHIs)

FIINIT(P1 ) = READFI ;FIINIT(U1 ) = READFI

FIINIT(V1 ) = READFI ;FIINIT(W1 ) = READFI

FIINIT(DEN1) = READFI ;FIINIT(PRPS) = -

1.000000E+00

Parent VR object for this patch is: B1

PATCH(OB1 ,INIVAL, 1, 67, 1, 60, 1, 81,

1, 1)

INIT(OB1 ,PRPS, 0.000000E+00, 1.980000E+02)

INIADD = F

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: INLET

PATCH(OB2 ,LOW , 1, 67, 1, 37, 1, 1,

1, 1)

COVAL(OB2 ,P1 , FIXFLU , 1.104711E-01)

COVAL(OB2 ,U1 , 0.000000E+00, 0.000000E+00)

COVAL(OB2 ,V1 , 0.000000E+00, 0.000000E+00)

COVAL(OB2 ,W1 , 0.000000E+00, 1.106670E-04)

Parent VR object for this patch is: OUTLET

PATCH(OB3 ,HIGH , 1, 67, 1, 37, 81, 81,

1, 1)

COVAL(OB3 ,P1 , 1.000000E+03, 0.000000E+00)

COVAL(OB3 ,U1 , 0.000000E+00, 0.000000E+00)

COVAL(OB3 ,V1 , 0.000000E+00, 0.000000E+00)

COVAL(OB3 ,W1 , 0.000000E+00, 0.000000E+00)

Parent VR object for this patch is: OUT_RAMI

PATCH(OB4 ,NORTH , 26, 42, 60, 60, 27, 42,

1, 1)

COVAL(OB4 ,P1 ,In-Form:-2.767000E-


COVAL(OB4 ,U1 , 0.000000E+00, 0.000000E+00)

COVAL(OB4 ,V1 , 0.000000E+00, 0.000000E+00)

COVAL(OB4 ,W1 , 0.000000E+00, 0.000000E+00)

XCYCLE = F

EGWF = F

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

***

Group 14. Downstream Pressure For PARAB

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

***

Group 15. Terminate Sweeps

LSWEEP = 3000 ;ISWC1 = 1

LITHYD = 1 ;LITFLX = 1 ;LITC = 1

;ITHC1 = 1

SELREF = T

RESFAC = 1.000000E-02

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

***

Group 16. Terminate Iterations

LITER (P1 ) = 15 ;LITER (U1 ) = 5

LITER (V1 ) = 5 ;LITER (W1 ) = 5

ENDIT (P1 ) = 1.000000E-03 ;ENDIT (U1 ) = 1.000000E-

03

ENDIT (V1 ) = 1.000000E-03 ;ENDIT (W1 ) = 1.000000E-

03

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

***

Group 17. Relaxation

RELAX(P1 ,LINRLX, 1.000000E+00)

RELAX(U1 ,FALSDT, 1.000000E+00)

RELAX(V1 ,FALSDT, 1.000000E+00)

RELAX(W1 ,FALSDT, 1.000000E+00)

RELAX(DEN1,LINRLX, 1.000000E+00)

RELAX(PRPS,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+06 ;VARMIN(U1 ) =-1.000000E+06

VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06

VARMAX(W1 ) = 1.000000E+06 ;VARMIN(W1 ) =-1.000000E+06

VARMAX(DEN1) = 1.000000E+10 ;VARMIN(DEN1) =-1.000000E+10

VARMAX(PRPS) = 1.000000E+10 ;VARMIN(PRPS) =-1.000000E+10

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

***

Group 19. Data transmitted to GROUND

USEGRD = T ;USEGRX = T

ASAP = T

PARSOL = T

SPEDAT(SET,DOMAIN,PHASE_1_MAT,I,67)

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

SPEDAT(SET,MAXINC,U1,R,1.00000E+01)

SPEDAT(SET,MAXINC,V1,R,1.00000E+01)

SPEDAT(SET,MAXINC,W1,R,1.00000E+01)

SPEDAT(SET,B1,DATFILE,C,tubostl)

SPEDAT(SET,OBJNAM,ÔB1,C,B1)

SPEDAT(SET,OBJTYP,ÔB1,C,BLOCKAGE)

SPEDAT(SET,B1,MATERIAL,R,1.98000E+02)

SPEDAT(SET,INLET,DATFILE,C,cylinder)

SPEDAT(SET,OBJNAM,ÔB2,C,INLET)


SPEDAT(SET,OUTLET,DATFILE,C,cylinder)

SPEDAT(SET,OBJNAM,ÔB3,C,OUTLET)

SPEDAT(SET,OBJTYP,ÔB3,C,OUTLET)

SPEDAT(SET,ARATIO,ÔB3,R,1.00000E+00)

SPEDAT(SET,OUT_RAMI,DATFILE,C,cylinder)

SPEDAT(SET,OBJNAM,ÔB4,C,OUT_RAMI)


SPEDAT(SET,SOURCE,R1!OB4,C,=-2.767000E-


SPEDAT(SET,OBJNAM,!OB5,C,B5)

SPEDAT(SET,OBJTYP,!OB5,C,NULL)

SPEDAT(SET,FACETDAT,NUMOBJ,I,5)

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:

* 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,Y,Y,Y,Y)

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

OUTPUT(W1 ,Y,N,Y,Y,Y,Y)

OUTPUT(DEN1,Y,N,Y,N,N,N)

OUTPUT(PRPS,Y,N,Y,N,N,N)

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

***

Group 22. Monitor Print-Out

IXMON = 67 ;IYMON = 54 ;IZMON =

74

NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP =

10001

UWATCH = F ;USTEER = F

HIGHLO = F

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

***

Group 23.Field Print-Out & Plot Control

NPRINT = 3000 ;NUMCLS = 5

NXPRIN = 13 ;IXPRF = 1 ;IXPRL =

67

NYPRIN = 12 ;IYPRF = 1 ;IYPRL =

60

NZPRIN = 16 ;IZPRF = 1 ;IZPRL =

10000

XZPR = F ;YZPR = F

IPLTF = 1 ;IPLTL = 3000 ;NPLT =

150

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

46

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.500E-03 5.901E-03 7.702E-03 9.504E-03

1.131E-02

1.311E-02 1.491E-02 1.671E-02 1.851E-02

2.031E-02

2.211E-02 2.347E-02 2.439E-02 2.530E-02

2.622E-02

2.714E-02 2.805E-02 2.897E-02 2.988E-02

3.080E-02

3.171E-02 3.263E-02 3.355E-02 3.446E-02

3.538E-02

3.629E-02 3.717E-02 3.802E-02 3.887E-02

3.971E-02

4.056E-02 4.141E-02 4.225E-02 4.310E-02

4.395E-02

4.479E-02 4.564E-02 4.649E-02 4.733E-02

4.818E-02

4.903E-02 4.990E-02 5.081E-02 5.171E-02

5.262E-02

5.352E-02 5.442E-02 5.533E-02 5.623E-02

5.714E-02

5.804E-02 5.895E-02 5.985E-02 6.075E-02

6.166E-02

6.256E-02 6.392E-02 6.574E-02 6.756E-02

6.938E-02

7.120E-02 7.302E-02 7.484E-02 7.665E-02

7.847E-02

8.029E-02 8.370E-02

Y-coordinates of the cell centres

2.500E-03 7.031E-03 1.109E-02 1.515E-02

1.921E-02

2.327E-02 2.734E-02 3.140E-02 3.546E-02

3.952E-02

4.358E-02 4.764E-02 5.170E-02 5.576E-02

5.982E-02

6.389E-02 6.630E-02 6.706E-02 6.783E-02

6.859E-02

6.935E-02 7.012E-02 7.088E-02 7.165E-02

7.241E-02

7.318E-02 7.394E-02 7.470E-02 7.547E-02

7.623E-02

7.700E-02 7.776E-02 7.853E-02 7.929E-02

8.005E-02

8.082E-02 8.160E-02 8.241E-02 8.322E-02

8.403E-02

8.484E-02 8.565E-02 8.646E-02 8.727E-02

8.808E-02

8.889E-02 8.970E-02 9.051E-02 9.364E-02

9.908E-02

1.045E-01 1.100E-01 1.154E-01 1.208E-01

1.263E-01

1.317E-01 1.372E-01 1.426E-01 1.480E-01

1.535E-01

Z-coordinates of the cell centres

4.333E-03 1.300E-02 2.167E-02 3.033E-02

3.900E-02

4.767E-02 5.633E-02 6.500E-02 7.367E-02

8.233E-02

9.100E-02 9.967E-02 1.083E-01 1.170E-01

1.257E-01

1.306E-01 1.317E-01 1.328E-01 1.340E-01

1.351E-01

1.363E-01 1.374E-01 1.385E-01 1.397E-01

1.408E-01

1.419E-01 1.431E-01 1.441E-01 1.449E-01

1.458E-01

1.466E-01 1.475E-01 1.483E-01 1.492E-01

1.500E-01

1.508E-01 1.517E-01 1.525E-01 1.534E-01

1.542E-01

1.551E-01 1.559E-01 1.568E-01 1.577E-01

1.586E-01

1.595E-01 1.604E-01 1.613E-01 1.622E-01

1.631E-01

1.640E-01 1.649E-01 1.658E-01 1.667E-01

1.676E-01

1.685E-01 1.694E-01 1.703E-01 1.712E-01

1.721E-01

1.730E-01 1.739E-01 1.748E-01 1.757E-01

1.766E-01

1.775E-01 1.821E-01 1.902E-01 1.983E-01

2.065E-01

2.146E-01 2.228E-01 2.309E-01 2.390E-01

2.472E-01

2.553E-01 2.634E-01 2.716E-01 2.797E-01

2.878E-01

2.960E-01

X-coordinates of the (higher) cell faces

5.000E-03 6.801E-03 8.603E-03 1.040E-02

1.221E-02

1.401E-02 1.581E-02 1.761E-02 1.941E-02

2.121E-02

2.301E-02 2.393E-02 2.485E-02 2.576E-02

2.668E-02

2.759E-02 2.851E-02 2.942E-02 3.034E-02

3.126E-02

3.217E-02 3.309E-02 3.400E-02 3.492E-02

3.583E-02

3.675E-02 3.760E-02 3.844E-02 3.929E-02

4.014E-02

4.098E-02 4.183E-02 4.268E-02 4.352E-02

4.437E-02

4.522E-02 4.606E-02 4.691E-02 4.776E-02

4.860E-02

4.945E-02 5.035E-02 5.126E-02 5.216E-02

5.307E-02

5.397E-02 5.488E-02 5.578E-02 5.668E-02

5.759E-02

5.849E-02 5.940E-02 6.030E-02 6.121E-02

6.211E-02

6.301E-02 6.483E-02 6.665E-02 6.847E-02

7.029E-02

7.211E-02 7.393E-02 7.574E-02 7.756E-02

7.938E-02

8.120E-02 8.620E-02

Y-coordinates of the (higher) cell faces

5.000E-03 9.061E-03 1.312E-02 1.718E-02

2.124E-02

2.531E-02 2.937E-02 3.343E-02 3.749E-02

4.155E-02

4.561E-02 4.967E-02 5.373E-02 5.779E-02

6.185E-02

6.592E-02 6.668E-02 6.744E-02 6.821E-02

6.897E-02

6.974E-02 7.050E-02 7.127E-02 7.203E-02

7.279E-02

7.356E-02 7.432E-02 7.509E-02 7.585E-02

7.661E-02

7.738E-02 7.814E-02 7.891E-02 7.967E-02

8.044E-02

8.120E-02 8.201E-02 8.282E-02 8.363E-02

8.444E-02

8.525E-02 8.606E-02 8.687E-02 8.768E-02

8.849E-02

8.930E-02 9.011E-02 9.092E-02 9.636E-02

1.018E-01

1.072E-01 1.127E-01 1.181E-01 1.236E-01

1.290E-01

1.344E-01 1.399E-01 1.453E-01 1.508E-01

1.562E-01

Z-coordinates of the (higher) cell faces

8.667E-03 1.733E-02 2.600E-02 3.467E-02

4.333E-02

5.200E-02 6.067E-02 6.933E-02 7.800E-02

8.667E-02

9.533E-02 1.040E-01 1.127E-01 1.213E-01

1.300E-01

1.311E-01 1.323E-01 1.334E-01 1.345E-01

1.357E-01

1.368E-01 1.380E-01 1.391E-01 1.402E-01

1.414E-01

1.425E-01 1.436E-01 1.445E-01 1.453E-01

1.462E-01

1.470E-01 1.479E-01 1.487E-01 1.496E-01

1.504E-01

1.513E-01 1.521E-01 1.530E-01 1.538E-01

1.547E-01

1.555E-01 1.564E-01 1.573E-01 1.582E-01

1.591E-01

1.600E-01 1.609E-01 1.618E-01 1.627E-01

1.636E-01

1.645E-01 1.654E-01 1.663E-01 1.672E-01

1.681E-01

1.690E-01 1.699E-01 1.708E-01 1.717E-01

1.726E-01

1.735E-01 1.744E-01 1.753E-01 1.762E-01

1.771E-01

1.780E-01 1.861E-01 1.943E-01 2.024E-01

2.105E-01

2.187E-01 2.268E-01 2.350E-01 2.431E-01

2.512E-01

2.594E-01 2.675E-01 2.756E-01 2.838E-01

2.919E-01

3.000E-01

--- INTEGRATION OF EQUATIONS BEGINS ---


1

47

Field Values of P1


none


none

IY= 36 none none -4.647E-01 -4.813E-01

none

IY= 24 none -1.824E-01 -3.376E-01 -3.468E-01

-1.940E-01

IY= 12 none -6.417E-02 -7.608E-02 -7.657E-02

-6.551E-02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 -4.744E-02

IX= 66

Field Values of U1


none


none

IY= 36 none none 1.534E-03 -1.543E-03

none

IY= 24 none 3.709E-03 1.902E-03 -1.904E-03

-3.685E-03

IY= 12 none 1.021E-03 4.086E-04 -4.055E-04

-1.012E-03

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 none

IX= 66

Field Values of V1


none


none

IY= 35 none none 4.835E-04 5.349E-04

none

IY= 23 none 3.300E-03 4.490E-03 4.578E-03

3.285E-03

IY= 11 none 2.187E-03 2.701E-03 2.722E-03

2.249E-03

IX= 1 14 27 40

53

IY= 59 none

IY= 47 none

IY= 35 none

IY= 23 none

IY= 11 1.160E-03

IX= 66

Field Values of W1


none


none

IY= 36 none none 6.142E-03 9.050E-03

none

IY= 24 none 2.698E-02 3.651E-02 3.680E-02

2.883E-02

IY= 12 none 2.761E-02 2.775E-02 2.776E-02

2.762E-02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 1.133E-03

IX= 66



none


none

IY= 36 none none 9.982E+02 9.982E+02

none

IY= 24 none 9.982E+02 9.982E+02 9.982E+02

9.982E+02

IY= 12 none 9.982E+02 9.982E+02 9.982E+02

9.982E+02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 9.982E+02

IX= 66



blockage


blockage


blockage


pil prop


pil prop

IX= 1 14 27 40

53

IY= 60 blockage

IY= 48 blockage

IY= 36 blockage

IY= 24 blockage

IY= 12 pil prop

IX= 66


1

Field Values of P1

IY= 60 none none -4.296E+01 -4.310E+01

none

IY= 48 none none -3.715E+01 -3.714E+01

none

IY= 36 none none -3.431E+01 -3.105E+01

none

IY= 24 none -1.537E-01 -7.606E-01 -8.201E-01

-1.764E-01

IY= 12 none -3.688E-02 -4.547E-02 -4.583E-02

-3.782E-02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 -2.435E-02

IX= 66

Field Values of U1

IY= 60 none none 2.056E-03 -2.087E-03

none

IY= 48 none none 6.552E-03 -6.584E-03

none

IY= 36 none none 1.309E-01 -1.309E-01

none

IY= 24 none 1.128E-02 1.324E-02 -1.323E-02

-1.134E-02

IY= 12 none 1.665E-03 7.086E-04 -7.046E-04

-1.652E-03

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 none

IX= 66

Field Values of V1

IY= 59 none none 4.178E-02 1.704E-01

none

IY= 47 none none 4.925E-02 1.572E-01

none

IY= 35 none none 1.390E-01 1.759E-01

none

IY= 23 none 8.386E-03 2.743E-02 2.940E-02

8.951E-03

IY= 11 none 3.054E-03 4.094E-03 4.139E-03

3.172E-03

IX= 1 14 27 40

53

IY= 59 none

IY= 47 none

IY= 35 none

IY= 23 none

IY= 11 8.269E-04

IX= 66

Field Values of W1

IY= 60 none none -2.120E-02 -2.805E-02

none

IY= 48 none none -1.314E-02 -1.224E-02

none

48

IY= 36 none none 7.791E-02 5.796E-02

none

IY= 24 none 2.044E-02 3.220E-02 3.252E-02

2.262E-02

IY= 12 none 2.647E-02 2.636E-02 2.636E-02

2.645E-02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 6.485E-05

IX= 66


IY= 60 none none 9.982E+02 9.982E+02

none

IY= 48 none none 9.982E+02 9.982E+02

none

IY= 36 none none 9.982E+02 9.982E+02

none

IY= 24 none 9.982E+02 9.982E+02 9.982E+02

9.982E+02

IY= 12 none 9.982E+02 9.982E+02 9.982E+02

9.982E+02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 9.982E+02

IX= 66



blockage


blockage


blockage


pil prop


pil prop

IX= 1 14 27 40

53

IY= 60 blockage

IY= 48 blockage

IY= 36 blockage

IY= 24 blockage

IY= 12 pil prop

IX= 66


1

Field Values of P1


none


none

IY= 36 none none 1.673E-01 1.783E-01

none

IY= 24 none 2.423E-02 7.608E-02 8.046E-02

2.590E-02

IY= 12 none -4.465E-03 -1.498E-03 -1.347E-03

-4.135E-03

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 -6.684E-03

IX= 66

Field Values of U1


none


none

IY= 36 none none -3.695E-04 3.953E-04

none

IY= 24 none 7.446E-03 6.450E-03 -6.452E-03

-7.528E-03

IY= 12 none 1.661E-03 6.914E-04 -6.884E-04

-1.648E-03

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 none

IX= 66

Field Values of V1


none


none

IY= 35 none none 6.939E-04 1.747E-03

none

IY= 23 none 5.050E-03 1.509E-02 1.581E-02

5.765E-03

IY= 11 none 2.946E-03 3.917E-03 3.958E-03

3.058E-03

IX= 1 14 27 40

53

IY= 59 none

IY= 47 none

IY= 35 none

IY= 23 none

IY= 11 2.247E-04

IX= 66

Field Values of W1


none


none

IY= 36 none none -4.031E-03 -6.329E-03

none

IY= 24 none 1.109E-02 1.455E-02 1.398E-02

1.286E-02

IY= 12 none 2.528E-02 2.471E-02 2.468E-02

2.521E-02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 -3.880E-04

IX= 66



none


none

IY= 36 none none 9.982E+02 9.982E+02

none

IY= 24 none 9.982E+02 9.982E+02 9.982E+02

9.982E+02

IY= 12 none 9.982E+02 9.982E+02 9.982E+02

9.982E+02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 9.982E+02

IX= 66



blockage


blockage


blockage


pil prop


pil prop

IX= 1 14 27 40

53

IY= 60 blockage

IY= 48 blockage

IY= 36 blockage

IY= 24 blockage

IY= 12 pil prop

IX= 66


1

Field Values of P1


none


none

IY= 36 none none 1.630E-01 1.667E-01

none

IY= 24 none 4.267E-02 1.074E-01 1.103E-01

4.580E-02

49

IY= 12 none 1.494E-02 2.129E-02 2.158E-02

1.566E-02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 6.721E-03

IX= 66

Field Values of U1


none


none

IY= 36 none none -8.879E-04 9.043E-04

none

IY= 24 none 1.679E-03 1.027E-03 -1.027E-03

-1.816E-03

IY= 12 none 1.258E-03 4.842E-04 -4.836E-04

-1.249E-03

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 none

IX= 66

Field Values of V1


none


none

IY= 35 none none -2.185E-04 -2.335E-04

none

IY= 23 none 1.697E-03 5.112E-03 5.200E-03

2.646E-03

IY= 11 none 2.328E-03 2.898E-03 2.920E-03

2.398E-03

IX= 1 14 27 40

53

IY= 59 none

IY= 47 none

IY= 35 none

IY= 23 none

IY= 11 -2.362E-04

IX= 66

Field Values of W1


none


none

IY= 36 none none 2.402E-03 3.522E-03

none

IY= 24 none 1.116E-02 2.078E-02 2.073E-02

1.342E-02

IY= 12 none 2.465E-02 2.400E-02 2.397E-02

2.458E-02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 -2.336E-04

IX= 66



none


none

IY= 36 none none 9.982E+02 9.982E+02

none

IY= 24 none 9.982E+02 9.982E+02 9.982E+02

9.982E+02

IY= 12 none 9.982E+02 9.982E+02 9.982E+02

9.982E+02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 9.982E+02

IX= 66



blockage


blockage


blockage


pil prop


pil prop

IX= 1 14 27 40

53

IY= 60 blockage

IY= 48 blockage

IY= 36 blockage

IY= 24 blockage

IY= 12 pil prop

IX= 66


1

Field Values of P1


none


none

IY= 36 none none 7.313E-03 7.380E-03

none

IY= 24 none 9.723E-03 1.392E-02 1.407E-02

1.015E-02

IY= 12 none 1.540E-02 1.860E-02 1.868E-02

1.590E-02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 5.167E-03

IX= 66

Field Values of U1


none


none

IY= 36 none none -4.441E-05 4.430E-05

none

IY= 24 none -1.785E-03 -3.013E-04 3.084E-04

1.796E-03

IY= 12 none -1.552E-04 -3.742E-05 3.863E-05

1.551E-04

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 none

IX= 66

Field Values of V1


none


none

IY= 35 none none 1.261E-05 2.466E-05

none

IY= 23 none -9.219E-05 8.975E-04 9.186E-04

1.082E-04

IY= 11 none 6.177E-04 5.731E-04 5.718E-04

6.182E-04

IX= 1 14 27 40

53

IY= 59 none

IY= 47 none

IY= 35 none

IY= 23 none

IY= 11 -3.868E-04

IX= 66

Field Values of W1


none


none

IY= 36 none none 1.705E-03 2.193E-03

none

IY= 24 none 2.125E-02 2.432E-02 2.432E-02

2.213E-02

IY= 12 none 2.416E-02 2.409E-02 2.407E-02

2.424E-02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

50

IY= 36 none

IY= 24 none

IY= 12 2.344E-03

IX= 66



none


none

IY= 36 none none 9.982E+02 9.982E+02

none

IY= 24 none 9.982E+02 9.982E+02 9.982E+02

9.982E+02

IY= 12 none 9.982E+02 9.982E+02 9.982E+02

9.982E+02

IX= 1 14 27 40

53

IY= 60 none

IY= 48 none

IY= 36 none

IY= 24 none

IY= 12 9.982E+02

IX= 66



blockage


blockage


blockage


pil prop


pil prop

IX= 1 14 27 40

53

IY= 60 blockage

IY= 48 blockage

IY= 36 blockage

IY= 24 blockage

IY= 12 pil prop

IX= 66

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

***

Whole-field residuals before solution

with resref values determined by EARTH

& resfac= 1.000000E-02

variable resref (res sum)/resref (res sum)

P1 3.134E-07 8.825E-03 2.766E-09

U1 1.115E-09 1.936E+01 2.158E-08

V1 7.576E-09 2.817E+00 2.134E-08

W1 7.530E-09 3.655E+00 2.752E-08

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

***

Sources and sinks

Nett source of U1 at patch named: OB3 (OUTLET )

=-3.980025E-08

Nett source of U1 at patch named: OB4 (OUT_RAMI)

= 7.931286E-07

pos. sum= 7.931286E-07 neg. sum=-3.980025E-08

nett sum= 7.533284E-07

Nett source of V1 at patch named: OB3 (OUTLET )

=-1.177284E-05

Nett source of V1 at patch named: OB4 (OUT_RAMI)

=-6.094207E-03

pos. sum= 0.000000E+00 neg. sum=-6.105980E-03

nett sum=-6.105980E-03

Nett source of W1 at patch named: OB2 (INLET )

= 2.685783E-03

Nett source of W1 at patch named: OB3 (OUTLET )

=-1.687669E-03

Nett source of W1 at patch named: OB4 (OUT_RAMI)

= 2.980590E-04

pos. sum= 2.983843E-03 neg. sum=-1.687669E-03

nett sum= 1.296173E-03

Nett source of R1 at patch named: OB2 (INLET )

= 1.104712E-01

Nett source of R1 at patch named: OB3 (OUTLET )

=-8.285094E-02

Nett source of R1 at patch named: OB4 (OUT_RAMI)

=-2.762031E-02

pos. sum= 1.104712E-01 neg. sum=-1.104712E-01

nett sum= 0.000000E+00

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

***

spot values vs sweep or iteration number

IXMON= 67 IYMON= 54 IZMON= 74

TIMESTEP= 1


ISWP P1 U1 V1 W1

1 0.000E+00 0.000E+00 0.000E+00 0.000E+00

151 0.000E+00 0.000E+00 0.000E+00 0.000E+00

301 0.000E+00 0.000E+00 0.000E+00 0.000E+00

451 0.000E+00 0.000E+00 0.000E+00 0.000E+00

601 0.000E+00 0.000E+00 0.000E+00 0.000E+00

751 0.000E+00 0.000E+00 0.000E+00 0.000E+00

901 0.000E+00 0.000E+00 0.000E+00 0.000E+00

1051 0.000E+00 0.000E+00 0.000E+00 0.000E+00

1201 0.000E+00 0.000E+00 0.000E+00 0.000E+00

1351 0.000E+00 0.000E+00 0.000E+00 0.000E+00

1501 0.000E+00 0.000E+00 0.000E+00 0.000E+00

1651 0.000E+00 0.000E+00 0.000E+00 0.000E+00

1801 0.000E+00 0.000E+00 0.000E+00 0.000E+00

1951 0.000E+00 0.000E+00 0.000E+00 0.000E+00

2101 0.000E+00 0.000E+00 0.000E+00 0.000E+00

2251 0.000E+00 0.000E+00 0.000E+00 0.000E+00

2401 0.000E+00 0.000E+00 0.000E+00 0.000E+00

2551 0.000E+00 0.000E+00 0.000E+00 0.000E+00

2701 0.000E+00 0.000E+00 0.000E+00 0.000E+00

2851 0.000E+00 0.000E+00 0.000E+00 0.000E+00

Variable 1 = P1 2 = U1 3 = V1 4 = W1

Minval= 0.000E+00 0.000E+00 0.000E+00 0.000E+00

Maxval= 0.000E+00 0.000E+00 0.000E+00 0.000E+00

Cellav= 0.000E+00 0.000E+00 0.000E+00 0.000E+00

1.00 +....+....+....+....+....+....+....+....+....+....+

. .

0.90 + +

. .

0.80 + +

. .

0.70 + +

. .

0.60 + +

. .

0.50 + +

. .

0.40 + +

. .

0.30 + +

. .

0.20 + +

. .

0.10 + +

. .

0.00 4..4.4..4.+4.4.+4.4.+4..4+4..4+.4.4+.4.4+.4..4.4..4

0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0


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

***

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

***

residuals vs sweep or iteration number


ISWP P1 U1 V1 W1

1 1.241E+01 3.377E+01 6.483E+01 2.183E+04

151 8.122E-02 6.263E+01 1.213E+01 5.698E+01

301 1.025E-01 3.043E+01 6.270E+00 1.600E+01

451 5.695E-02 3.319E+01 6.041E+00 3.967E+01

601 7.341E-02 2.562E+01 4.795E+00 1.117E+01

751 4.150E-02 2.870E+01 4.937E+00 2.876E+01

901 5.337E-02 2.336E+01 4.020E+00 8.143E+00

1051 3.062E-02 2.505E+01 4.196E+00 2.107E+01

1201 3.854E-02 2.144E+01 3.438E+00 6.372E+00

1351 2.278E-02 2.265E+01 3.626E+00 1.548E+01

1501 2.876E-02 2.087E+01 3.183E+00 5.332E+00

1651 1.759E-02 2.160E+01 3.269E+00 1.176E+01

1801 2.148E-02 2.033E+01 2.956E+00 4.739E+00

1951 1.395E-02 2.061E+01 3.081E+00 8.923E+00

2101 1.675E-02 1.968E+01 2.941E+00 4.344E+00

2251 1.129E-02 2.010E+01 2.913E+00 7.039E+00

2401 1.292E-02 1.959E+01 2.856E+00 4.007E+00

2551 9.578E-03 2.013E+01 2.898E+00 5.645E+00

2701 1.074E-02 1.928E+01 2.772E+00 3.888E+00

2851 8.183E-03 1.953E+01 2.817E+00 4.837E+00

Variable 1 = P1 2 = U1 3 = V1 4 = W1

Minval= -4.806E+00 2.959E+00 1.020E+00 1.358E+00

Maxval= 2.518E+00 4.137E+00 4.172E+00 9.991E+00

1.00 4..2.+....+....+....+....+....+....+....+....+....+

. .

0.90 + +

. .

51

0.80 + +

. .

0.70 + +

. .

0.60 + +

. .

0.50 2 +

. 3 2 .

0.40 + 2 +

. 1 2 .

0.30 + 4 1 +

. 3 4 2 4 1 .

0.20 + 3 4 1 +

. 4 3 2 3 4 1 4 1 .

0.10 + 4 4 2 3 2 4 1 +

. 4 4 3 2 3 4 1 4 1 4

0.00 +....+....+....+....+....+....+.4..+.4.3+.4..3.4..3

0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0


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

***

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

***

SATLIT RUN NUMBER = 1 ; LIBRARY REF.= 0

Run completed at 19:20:08 on Tuesday, 16 June 2009

MACHINE-CLOCK TIME OF RUN = 5486 SECONDS.

TIME/(VARIABLES*CELLS*TSTEPS*SWEEPS*ITS) = 5.486E+03

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

***

Documents

EM 974 – Métodos Computacionais Em Engenharia Térmica …