UNIVERSIDADE ESTADUAL DE CAMPINAS

FACULDADE DE ENGENHARIA MECÂNICA

EM 974 Métodos Computacionais em Engenharia Térmica e Ambiental

Prof. Responsável: Eugênio Spanó Rosa

IDENTIFICAÇÃO

NOME RA Bruno Rugno Delati 085784

Andre Stuart Nogueira 085756

TURMA: B GRUPO: 6

TÍTULO DO TRABALHO

Análise Aerodinâmica da Carenagem do Formula Elétrico 2012.

AVALIAÇÃO ETAPA IV

1.

(20%)

Apresentação e Organização: o texto é claro e objetivo, a formatação

do trabalho apresenta o trabalho de forma organizada e de fácil leitura,

as tabelas e gráficos complementam as informações, os gráficos são

claros e objetivos, as variáveis utilizadas foram definidas

propriamente, as variáveis possuem definição das dimensões.

Bom Médio Fraco

2.

(10%)

Introdução: apresentar a motivação que levou a desenvolver o

trabalho, em que área ele se aplica e o objetivo do trabalho, isto é, o

que o grupo pretende alcançar.

Bom Médio Fraco

3.

(10%)

Revisão da Literatura: tomar conhecimento se há trabalhos similares

na literatura, se há dados experimentais disponíveis.

Bom Médio Fraco

4.

(20%)

Implementação no Phoenics: anexar o arquivo Q1 e destacar em texto,

os grupos do Q1 que contêm as maiores contribuições do

desenvolvimento do projeto. Deixar claro o domínio computacional,

as condições de contorno empregadas e as propriedades dos materiais.

Bom Médio Fraco

4.

(20%)

Resultados numéricos: apresentar teste de malha e os resíduos

numéricos. Apresentar os resultados numéricos em termos de gráficos

do problema juntamente com um texto explicando o significado dos

gráficos.

Bom Médio Fraco

5.

(20%)

Análise: nesta seção o grupo vai interpretar os resultados obtidos para:

fundamentar como se comporta o fenômeno estudado e tirar

conclusões de projeto. Por último é apresentado uma conclusão geral

do trabalho.

Bom Médio Fraco

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

Introdução A SAE (Society of Automotive Egineers) é uma organização fundada nos Estados Unidos por engenheiros do setor da mobilidade, terrestre, aéreo e naval. Na década de 80 a SAE International começou a promover a atividade estudantil denominada de Formula SAE, para incentivar jovens estudantes de engenharia a ingressar no mercado da engenharia de mobilidade e transportes. O que começou nos Estados Unidos, tomou enormes proporções sendo que hoje em dia esta aividade é denominada Formula Series, com competições nos Estados Unidos, Australia, Italia, Reino Unido, Alemanha e Brasil. Esta ultima organizada pela SAE Brasil. Aqui, o evento começou a ser promovido em 2004 a fim de aprimorar a formação profissional dos participantes, incentivar o trabalho em equipe, levantar fundos, planejar e analisar os custos através de simulação de desenvolimento de um novo produto. É um desafio de engenharia, não se limita ao desenvolvimento dos protótipos em pista, mas sim a todo o desenvolvimento, analisando também o projeto e a viabilidade de produção. Os dois melhores coloados na competição brasileira ganham as inscrições para disputar as etapas norte americanas. A competição Formula SAE Eletrico é a mais nova competião estudantil de engenharia do país. O principal objetivos é que as equipes desenvolvam e consruam protótipos puramente elétricos para participar de provas teoricas e práticas presentes na competição. Além disso a competição integra os estudates às empresas , já que a parceria entre equipes e emprsas é imprescindível para a arrecadação de fundos, know-how, maquinário e outrs elementos indispensáveis para a construção do carro. A equipe Unicamp E-Racing foi fudada em Novemro de 2011 com o intuito de ser a melhor equipe brasileira na categoria Formula SAE Eletrico. A equipe competirá na 1ª Competição de Formula SAE Eletrico do Brasil. A competição é divida entre provas estáticas e provas dinâmicas. Nas estáticas o carro não está na pista enquanto que as dinâmicas são todas provas contra o relógio. Provas Estáticas: Inspeção de segurança Determina a possibilidade ou não de o veículo e os pilotos competirem nas provas dinâmicas. Os juízes avaliam se os participantes estão de acordo com as regras. Nesta prova não há pontuação. Prova de Custos É feita uma análise de todos os custos do projeto e seu processo de manufatura, privilegiando os projetos mais baratos e que se encaixam melhor no conceito do lean manufacturing. Prova de Design Apresentação do projeto aos juízes, onde são discutidas as soluções, o projeto, os ajustes e as validações. Apresentação Os participantes devem desenvolver um business case, no qual se defende que o veículo atende melhor às demandas de mercado e é economicamente viável para produção em grande escala. Após passar pela inspeção de segurança, o veículo está autorizado a participar dos eventos dinâmicos, onde a performance será testada na pista. As provas dinâmicas são constituídas por: Aceleração Os veículos devem percorrer uma pista de 75m partindo do zero. O melhor tempo ganha.

Skid Pad Um circuito em forma de "8", no qual é testada a dinâmica em curva, também premiando o menor tempo. Autocross Um circuito completo de aproximadamente 1/2 milha, onde se avalia a manobrabilidade e dinâmica dos veículos. Enduro Uma prova de 22 km, onde os carros são levados ao seu limite, testando a performance, confiabilidade e resistência dos protótipos. Aerodinâmica na Formula SAE Na ultima competição de Formula SAE em Michigam, considerada o “mundial” dentre todas as etapas, 120 times de inumeras universidades do planeta se registraram para o evento. Como esse número de equipes vem crescendo com o passar dos anos - no Brasil também é interessante notar que em 2008 tinhamos 17 equipes enquanto que em 2011, 22 equipes disputaram – a performance dos prototipos tambem vem crecendo junto. Muita tenologia tem sido desenvolida para este tipo de competição, como aerofólios, difusores, apêndices. Melhorias aerodinamica tem se tornado comum nas equipes de Formula SAE. Essa melhorias permitem com que o veíulos adiquiram maior tração ao mesmo tempo que reduzem o peso do sistema, porque com estudos aerodinamicos é possivel tem um aumento de downfoce e redução no drag, ganhando estabilidade e performance. Como qualquer outra competição internacional, o Formula SAE Brasil está se tornando cada vez mais competitivo e de alto nivel e assim tambem espera-se que aconteça com o Formula SAE Eletrico.

O principal objetivo deste trabalho é de levantar o valor do Coeficiente de arrasto do protótipo de Formula SAE Elétrico da equipe Unicamp E-Racing. A princípio a ideia era de ser realizada uma simulação 3D do veículo, porém tivemos que mudar os planos, pois tal simulação demandaria muito processamento e tempo. Com isso optamos por rodar uma simulação 2D do veículo, na qual pudermos calcular o valor aproximado do coeficiente de arrasto e analisar a esteira de ar no plano XY. O domínio então foi estabelecido com X representando o comprimento, Y a altura e Z a profundidade (unitária porque a simulação foi 2D). Com X=20m, Y=5m e Z=0,1m. Foi adotado esses valores grandes em X e Y para simular o formula elétrico andando em espaço aberto.

Figura 1 - Domínio 2D utilizado.

Além disso, temos um INLET no plano YZ referenciado em x=0 e um OUTLET no plano YZ referenciado em x=20. Também foi estabelecido o chão como um PLATE. As características destes são:

INLET: velocidade em x=16.67m/s;

OUTLET: manteve a pressão de referencia constante;

PLATE: velocidade em x=16.67m/s. Geometria O assembly completo do Formula Elétrico foi desenhado e montado no Software Creo da PTC. Para a simulação 2D produzimos uma fatia na porção central do assembly de modo a ter a maior área de seção longitudinal. Este arquivo foi exportado para .STL para que fosse reconhecido pelo PHOENICS e foi alocado como um BLOCKAGE.

Figura 2 - Seção do Formula Elétrico com simulação do piloto com capacete.

Escolha do modelo de turbulência Para se definir o modelo de turbulência antes é necessário o calculo do numero de Reynolds do escoamento. Tomando as medidas do domínio como 0,1m x 5m, pode-se calcular o diâmetro hidráulico como:

At – Área da seção transversal At = 0,1 x 5 = 0,5m². P – Perímetro molhado P = 2*(0,1+5) = 10,2m. Isso nos dá Dh =0,1961 m. E o número de Reynolds, por sua vez, é calculado por:

V=17m/s (60km/h,adotada); =0.1961m; =1.554. m²/s. Temos então Re=210359,524 . Para valores acima de 200.000, o escoamento é considerado turbulento. Portanto escolhemos o modelo de turbulência KMODEL.

Critérios Numéricos Para um primeiro teste colocamos 1500 iterações e uma malha com 1040 elementos distribuidos em Nx=40 e Ny=26.

Figura 3 - Malha com 1040 elementos.

Assim,

Figura 4 - Mapa de velocidades para malha de 1040 elementos.

Figura 5 - Mapa de pressões para malha de 1040 elementos.

Os erros associados a cada variável estão representados na figura a baixo:

Figura 6 - Erros para malha de 1040 elementos.

As forças sobre o corpo são retiradas do arquivo result: Integrated force on object: FORMULA ----------------------------------- Fx = 1.099308E+01 (Pressure= 1.088846E+01, Friction= 1.046255E-01) Fy = 1.852440E+01 (Pressure= 1.851572E+01, Friction= 8.681241E-03) Force unit vector: 5.103403E-01 8.599725E-01 1.499462E-04 Além disso os valores de resíduo para esta simulação são mostrados: Whole-field residuals before solution with resref values determined by EARTH & resfac=1.0E-03 variable resref (res sum)/resref (res sum) P1 4.179E-04 1.140E-01 4.765E-05 U1 7.299E-03 1.718E+00 1.254E-02 V1 4.115E-04 5.504E+00 2.265E-03 KE 7.739E-04 4.027E+00 3.116E-03 EP 1.059E-02 1.555E+01 1.647E-01 Temos então que refinar a malha para avaliar se não as forças se estabilizam. Quando isso acontecer quer dizer que não é necessário mais refiná-la. Com isso montamos uma malha com 3680 elementos distribuidos em Nx=80 e Ny=46. Aqui tivemos que aumentar para 5000 iterações para que a simulação pudesse convergir.

ERROS

Figura 7 - Malha com 3680 elementos.

Assim,

Figura 8 - mapa de velocidades para malha de 3680 elementos.

Figura 9 - Mapa de pressões para malha de 3680 elementos.

Novamente a simulação convergiu, como pode ser observado a baixo:

Figura 10 - Erros para malha de 3680 elementos.

As forças sobre o corpo deram: Integrated force on object: FORMULA ----------------------------------- Fx = 1.690814E+01 (Pressure= 1.679388E+01, Friction= 1.142634E-01) Fy = 8.725848E+00 (Pressure= 8.714083E+00, Friction= 1.176525E-02) Force unit vector: 8.886406E-01 4.586040E-01 3.794711E-04 Temos aqui uma variação em relação à malha anterior, o que nos sugere refinarmos ainda mais a malha. Álém disso os resíduos foram: Whole-field residuals before solution with resref values determined by EARTH & resfac=1.0E-03 variable resref (res sum)/resref (res sum) P1 2.418E-04 1.803E-01 4.361E-05 U1 4.368E-03 3.006E+00 1.313E-02 V1 3.327E-04 4.552E+00 1.515E-03 KE 6.167E-04 1.731E+00 1.067E-03 EP 3.818E-03 1.298E+00 4.954E-03 Com isso temos agora uma malha com 5320 elementos distribuidos em Nx=95 e Ny=56. Aqui também tivemos que usar 5000 iterações.

ERROS

Figura 11 - Malha com 5320 elementos.

Assim,

Figura 12 Mapa de velocidades para malha com 5320 elementos.

Figura 13 - Mapa de pressões para malha com 5320 elementos.

E os erros:

Figura 14 - Erros para malha com 5320 elementos.

As forças sobre o corpo: Integrated force on object: FORMULA ----------------------------------- Fx = 1.659012E+01 (Pressure= 1.647946E+01, Friction= 1.106615E-01) Fy = 9.580073E+00 (Pressure= 9.566530E+00, Friction= 1.354339E-02) Force unit vector: 8.659855E-01 5.000690E-01 2.178543E-04 Os resíduos para esta malha foram: Whole-field residuals before solution with resref values determined by EARTH & resfac=1.0E-03 variable resref (res sum)/resref (res sum) P1 1.994E-04 2.309E-01 4.602E-05 U1 3.594E-03 2.379E+00 8.550E-03 V1 2.818E-04 3.428E+00 9.662E-04 KE 5.164E-04 2.985E+00 1.542E-03 EP 2.747E-03 1.846E+00 5.070E-03 Foi montada mais uma configuração esta com 6160 elementos em Nx=110 e Ny=56.

ERROS

Figura 15 - Malha com 6160 elementos.

Assim,

Figura 16 - Mapa de velocidades para malha de 6160 elementos.

Figura 17 - Mapa de pressões para malha de 6160 elementos.

E os erros:

Figura 18 - erros para malha com 6160 elementos.

As forças sobre o corpo: Integrated force on object: FORMULA ----------------------------------- Fx = 1.572264E+01 (Pressure= 1.561001E+01, Friction= 1.126203E-01) Fy = 9.973325E+00 (Pressure= 9.961365E+00, Friction= 1.196007E-02) Force unit vector: 8.444388E-01 5.356520E-01 1.380849E-04 E os resíduos: Whole-field residuals before solution with resref values determined by EARTH & resfac=1.0E-03 variable resref (res sum)/resref (res sum) P1 1.953E-04 2.315E-01 4.520E-05 U1 3.523E-03 8.907E-01 3.138E-03 V1 2.690E-04 1.979E+00 5.323E-04 KE 4.675E-04 2.594E+00 1.212E-03 EP 2.273E-03 1.605E+00 3.647E-03 Podemos ver que as simulações convergiram, pois os erros estão bem pequenos e além disso pode-se ver que os resíduos estão extremamente pequenos também o que confere à simulação bons resultados. Mesmo sem a estabilização nos valores das forças, estas estão com a mesma ordem de grandeza. Isso nos leva a parar com o refino da malha. Muito provavelmente se continuarmos com o refinamento, não teremos mudanças nos resultados obtidos. Somado a isso se pode notar que a velocidade média na região do probe se manteve constante nas ultimas 3 simulações, com mudanças apenas nos centésimos, o que reafirma que a malha não precisa ser mais refinada.

ERROS

Pode-se ver pela figura 17 que a maior parte do downforce está sendo gerado na parte de trás do formula. Isso é representado pelas pressões negativas. Na frente ocorre o contrário o que nos leva a pensar que o uso de uma asa dianteira seria uma boa ideia para induzir o downforce na dianteira. Resultados Os resultados das forças e velocidade média no probe para cada malha podem ser observados na seguinte tabela:

MALHA (elementos) Fx (N) Fy (N) Vx_média (m/s)

1040 10,9931 18,5244 16,83987

3680 16,9081 8,7258 16,93920

5320 16,5901 9,5801 16,93925

6160 15,7226 9,9733 16,93620

Lembrando que o Probe foi alocado na esteira de ar do escoamento.

O coeficiente de arrasto, Cx é calculado pela expressão:

Onde Fr é a força de arrasto, é a densidade do fluido, U é a velocidade do fluido e A é a área da seção transversal à direção do escoamento. Calculando a área da maior seção transversal da fatia do Formula Eletrico como:

( )

Sendo 0,1m a largura do domínio e 1,1m a altura máxima na maior seção transversal. Assumindo que ρ=1,225 kg/m³. Temos:

Ux Fx Cd

16,93987 10,9931 0,5754

16,9392 169081 0,8746

16,93925 16,5901 0,8582

16,9362 15,7226 0,8136

A seguir temos o gráfico que representa a tabela. Podemos ver que o C cresce linearmente com o refinar da malha. Lembrando-se que estes valores são apenas comparativos, pois a simulação é 2D. Podemos comparar esses valores mas não podemos dizer que estes valores são reais para o veículo completo.

Seguindo o processo descrito por McBeath [1] e também por Wordley [9], a velocidade máxima teorica de um veículo sem aparatos aerodinâmicos pode ser calculado pela equação do BHP. Desconsiderando a resistencia a rolamento o BHP no S.I pode ser calculado por:

Sendo A=0,11m². Levando-se em conta que a potencia do Formula Eletrico é de 45kW, rearranjando a equação a cima paa a velocidade, temos:

√

Essa velocidade é um absurdo e está muito acima do que é visto nas competçõe de Formula SAE e Formula SAE Eletrico, então assumindo uma nova velocidade máxima em torno de 120km/h será possíve calcular o BHP que será sacrificado para vencer a resistencia do ar:

Com isso para o Cd=0,8136 que foi o valor obtido através da malha mais refinada, temos que o BHP é 2,029kW, ou seja, para vencer a resistencia do ar o Formula Eetrico irá sacrificar 2,029kW de potencia para isso.

0

0,1

0,2

0,3

0,4

0,5

0,6

0,7

0,8

0,9

1

0 5 10 15 20

Cd

Força de arrasto em "x"

Cd X Força

Series1

Conclusões finais Conseguimos obter valores aceitaveis de Cd paa um veíulo que possui uma careagem que privilegia o escoamento mas que não possui aerofolios, difusor ou outros apendices aerodinamicos. Estes valores de Cd são para fins comparativos entre as malhas vito que a simulação é 2D. Lógicamente devem estar perto dos valores reais do veículo, porem uma analise eficiente 3D demandaria muito tempo e processamento, não sendo compatível com o intuito da disciplina. A simulação 3D é um bom tema para um trabalho futuro. Além disso pela analise do BHP sacrificado, podemos ter noção do quanto de potencia perderemos para vencer a resistencia do ar. Um carro que tenha menor Cd, com certeza sacificará menor potencia do motor paravencer essa resitencia, e é isso um dos faores que leva ao estudo aerodinamico de veiculos. Isso é muito importante no setor omercial de linha das industias automobilisticas, porque scificando menos potencia para vencer a resistencia do ar, pode-se ter um carro mais econômico, e economia de combustivel “pesa muito” no bolso do consumidor nos dias de hoje. A grande parte da população que compra carros para locomoção, ir e vi do trabalho e etc necessita de um carro economico, que gaste menos gasolina, e isto pode ser alcançado dentre tantos fatores pela redução do BHP sacrificado. E para essa redução é necessário a redução do Cd do veículo. Mesmo sem melhorias aerodinâmicas (tais como aerofolios, difusor e etc), no presente momento o projeto do Formula elétrico conta com uma carenagem com Cd da ordem de grandeza de 0,85 o que é muito bom e que sacrifica um pouco mas de 2kW de potencia do veículo. Visto que a potencia do mesmo é de 45kW, temos 4% da potencia total sendo sacrificada para vencer a resistencia do ar, o que podemos dizer que não é muito. Por fim, a ferramenta computacional de CFD PHOENICS mostrou-se eficiente pois nos permitiu analisar o Cd de nosso veículo, mesmo que sendo numa simulação 2D, com resltados satisfatórios. O arquivo results em anexo é da simulação com a última malha, a mais refinada com 6160 elementos em Nx=110 e Ny=56. Bibliografia:

1. McBeth S., Competition Car Downforce, 1998. 2. Çengel Yunus A, Mecânica dos Fluidos, McGraw Hill, 2007. 3. Aono A./Souza Henrique de., Estudo do Escoamento em uma Carenagem de

Formula SAE pelo uso do Phoenics, 2009. Disciplina EM974/FEM. 4. Wordley Scott/Saunders Jeff, Aerodynamics for Formula SAE: Initial design and

performance prediction, SAE paper number: 2006-010806. 5. Liang Li, Flavio/Molina de Almeida, Ederson, Analise comparativa de motocicleta

esportiva com e sem o “efeito bolha”. 6. Gillespie, T.D “Fudamentals of Vehicle Dynamics”, Society of Automotive Engineers

Inc. 1992 7. http://www.cham.co.uk/ChmSupport/polis.php [POLIS] (acesso em 20/05/2012) 8. Ranieri, V. “Arrasto Aerodinâmico no Mini baja”. Disciplina EM974/FEM. 9. Wordley, S. “Aerodynamics for Formula SAE: Initial Design and Performance

Prediction”, Monash University – SAE International 2006.

ARQUIVO RESULT ************************************************************ --------------------------------------------------------- CCCC HHH PHOENICS October 2010 - EARTH CCCCCCCC H (C) Copyright 2010 CCCCCCC See H Concentration Heat and Momentum Ltd CCCCCCC our new H All rights reserved. CCCCCC Web-site H Address: Bakery House, 40 High St CCCCCCC www.cham. H Wimbledon, London, SW19 5AU CCCCCCC co.uk H Tel: 020-8947-7651 CCCCCCCC H Fax : 020-8879-3497 CCCC HHH E-mail: phoenics@cham.co.uk --------------------------------------------------------- This program forms part of the PHOENICS installation for: CHAM The code expiry date is the end of : Sep 2002 --------------------------------------------------------- Running with 32-bit executable ************************************************************ Initial estimated storage requirement is 10000000 Information about material properties Total number of SPEDATs is 17 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 PRPS is stored with initial value= =-1. Material properties used for phase 1 are density laminar viscosity turbulent viscosity mixing length >>> End of property-related data <<< ************************************************************ Number of F-array locations available is 10000000 Number used before BFC allowance is 748710 Number used after BFC allowance is 748710 >> Current turbulence model constants << They may be changed by inserting in Q1 SPEDAT(KECONST,name of constant,R,value) CMU =0.5478 CD =0.1643 CMUCD =0.090004 C1E =1.44 C2E =1.92 AK =0.41 EWAL =8.6

biggest cell volume divided by average is 12.78154 at: ix = 83 iy = 56 iz = 1 xg =12.066252 yg =4.640397 zg =0.05 smallest cell volume divided by average is 1.5713096E-02 at: ix = 32 iy = 2 iz = 1 xg =4.07875 yg =0.032429 zg =0.05 ratio of smallest to biggest is 1.2293590E-03 Number used after PARSOL allowance is 748710 ************************************************************ -------- Recommended settings ------- CONWIZ=T activates settings based on refrho =1. refvel =10. reflen =1. reftemp =1000. rlxdu1 =0.5 rlxdv1 =0.5 rlxdw1 =0.5 Maximum change of U1 per sweep= 100.0000 Maximum change of V1 per sweep= 100.0000 Maximum change of KE per sweep= 0.1000000 Maximum change of EP per sweep= 1000.000 relaxation and min/max values left at defaults may have been changed ************************************************************ ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(No title has been set for this run. ) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = T NX = 110 XULAST =20. XFRAC(1)=6.666667E-03 ;XFRAC(23)=0.153333 XFRAC(45)=0.239375 ;XFRAC(67)=0.379725 XFRAC(89)=0.697075 ************************************************************ Group 4. Y-Direction Grid Spacing NY = 56 YVLAST =5. YFRAC(1)=6.000001E-03 ;YFRAC(12)=0.037242 YFRAC(23)=0.129967 ;YFRAC(34)=0.152586 YFRAC(45)=0.253109 ;YFRAC(56)=1. ************************************************************

Group 5. Z-Direction Grid Spacing PARAB = F NZ = 1 ZWLAST =0.1 ************************************************************ Group 6. Body-Fitted Coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T NAME(1)=P1 ;NAME(3)=U1 NAME(5)=V1 ;NAME(12)=KE NAME(13)=EP ;NAME(144)=PRPS NAME(145)=SHRY ;NAME(146)=SHRX NAME(147)=EPKE ;NAME(148)=DEN1 NAME(149)=EL1 ;NAME(150)=ENUT * Y in SOLUTN argument list denotes: * 1-stored 2-solved 3-whole-field * 4-point-by-point 5-explicit 6-harmonic averaging SOLUTN(P1,Y,Y,N,N,N,Y) SOLUTN(U1,Y,Y,N,N,N,Y) SOLUTN(V1,Y,Y,N,N,N,Y) SOLUTN(KE,Y,Y,N,N,N,N) SOLUTN(EP,Y,Y,N,N,N,N) SOLUTN(PRPS,Y,N,N,N,N,N) SOLUTN(SHRY,Y,N,N,N,N,Y) SOLUTN(SHRX,Y,N,N,N,N,Y) SOLUTN(EPKE,Y,N,N,N,N,Y) SOLUTN(DEN1,Y,N,N,N,N,Y) SOLUTN(EL1,Y,N,N,N,N,Y) SOLUTN(ENUT,Y,N,N,N,N,Y) DEN1 = 148 VIST = 150 LEN1 = 149 PRPS = 144 ************************************************************ 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(KE,N,Y,Y,N,Y,N) TERMS(EP,N,Y,Y,N,Y,N) DIFCUT =0.5 ;ZDIFAC =1. GALA = F ;ADDDIF = T NEWENT = T ISOLX = -1 ;ISOLY = -1 ;ISOLZ = 0 ************************************************************ Group 9. Properties used if PRPS is not stored, and where PRPS = -1.0 if it is! RHO1 =1.189 ;TMP1 =0. EL1 = GRND4 TSURR =0. ;TEMP0 =273. PRESS0 =1.01325E+05 DVO1DT =3.41E-03 ;DRH1DP =0. EMISS =0. ;SCATT =0.

RADIA =0. ;RADIB =0. EL1A =0. ;EL1B =0. EL1C =0. ENUL =1.544E-05 ;ENUT = GRND3 ENUTA =0. ;ENUTB =0. ENUTC =0. IENUTA = 0 PRNDTL(U1)=1. ;PRNDTL(V1)=1. PRNDTL(KE)=1. ;PRNDTL(EP)=1. PRT(U1)=1. ;PRT(V1)=1. PRT(KE)=1. ;PRT(EP)=1.314 CP1 =1005. ;CP2 =1. ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initial field variables (PHIs) FIINIT(P1)=0. ;FIINIT(U1)=1.0E-10 FIINIT(V1)=1.0E-10 ;FIINIT(KE)=0.694722 FIINIT(EP)=0.380602 ;FIINIT(PRPS)=-1. FIINIT(SHRY)=1.0E-10 ;FIINIT(SHRX)=1.0E-10 FIINIT(EPKE)=1.0E-10 ;FIINIT(DEN1)=1.189 FIINIT(EL1)=1.0E-10 ;FIINIT(ENUT)=1.0E-10 Parent VR object for this patch is: FORMULA PATCH(OB1 ,INIVAL, 30, 61, 1, 32, 1, 1, 1, 1) INIT(OB1 ,PRPS,0. ,198. ) 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 ,WEST , 1, 1, 1, 56, 1, 1, 1, 1) COVAL(OB2 ,P1 , FIXFLU ,19.820629 ) COVAL(OB2 ,U1 ,0. ,16.67 ) COVAL(OB2 ,V1 ,0. ,0. ) COVAL(OB2 ,KE ,0. ,0.694722 ) COVAL(OB2 ,EP ,0. ,0.380602 ) Parent VR object for this patch is: OUTLET PATCH(OB3 ,EAST , 110, 110, 1, 56, 1, 1, 1, 1) COVAL(OB3 ,P1 ,1000. ,0. ) COVAL(OB3 ,U1 ,0. ,0. ) COVAL(OB3 ,V1 ,0. ,0. ) COVAL(OB3 ,KE ,0. , SAME ) COVAL(OB3 ,EP ,0. , SAME ) Parent VR object for this patch is: CHAO PATCH(OB4 ,SWALL , 1, 110, 1, 1, 1, 1, 1, 1) COVAL(OB4 ,U1 , GRND2 ,16.67 ) COVAL(OB4 ,KE , GRND2 , GRND2 )

COVAL(OB4 ,EP , GRND2 , GRND2 ) PATCH(KESOURCE,PHASEM, 1, 110, 1, 56, 1, 1, 1, 1) COVAL(KESOURCE,KE , GRND4 , GRND4 ) COVAL(KESOURCE,EP , GRND4 , GRND4 ) XCYCLE = F EGWF = T WALLCO = GRND2 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 15000 ;ISWC1 = 1 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1 SELREF = T RESFAC =1.0E-03 ************************************************************ Group 16. Terminate Iterations LITER(P1)=200 ;LITER(U1)=10 LITER(V1)=10 ;LITER(KE)=20 LITER(EP)=20 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03 ENDIT(V1)=1.0E-03 ;ENDIT(KE)=1.0E-03 ENDIT(EP)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(P1,LINRLX,0.5) RELAX(U1,LINRLX,0.5) RELAX(V1,LINRLX,0.5) RELAX(KE,LINRLX,0.5) RELAX(EP,LINRLX,0.5) RELAX(PRPS,LINRLX,1.) RELAX(SHRY,LINRLX,1.) RELAX(SHRX,LINRLX,1.) RELAX(EPKE,LINRLX,1.) RELAX(DEN1,LINRLX,0.5) RELAX(EL1,LINRLX,1.) RELAX(ENUT,LINRLX,0.5) KELIN = 3 OVRRLX =0. EXPERT = F ;NNORSL = F ************************************************************ Group 18. Limits VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10 VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06 VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06 VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10 VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10 VARMAX(SHRY)=1.0E+10 ;VARMIN(SHRY)=-1.0E+10 VARMAX(SHRX)=1.0E+10 ;VARMIN(SHRX)=-1.0E+10 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=1.0E-10 VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=1.0E-06 VARMAX(EL1)=1.0E+10 ;VARMIN(EL1)=1.0E-10 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=1.0E-10 ************************************************************ Group 19. Data transmitted to GROUND

GENK = T CONWIZ = T CALFOR = T GEN1 = 681 ISG50 = 1 SPEDAT(SET,DOMAIN,PHASE_1_MAT,I,0) SPEDAT(SET,DRAG,MOMCENX,R,0.) SPEDAT(SET,DRAG,MOMCENY,R,0.) SPEDAT(SET,DRAG,MOMCENZ,R,0.) SPEDAT(SET,FORMULA,DATFILE,C,eletrico_side) SPEDAT(SET,OBJNAM,^OB1,C,FORMULA) SPEDAT(SET,OBJTYP,^OB1,C,BLOCKAGE) SPEDAT(SET,FORMULA,MATERIAL,R,198.) SPEDAT(SET,OBJNAM,!OB2,C,INLET) SPEDAT(SET,OBJTYP,!OB2,C,INLET) SPEDAT(SET,OBJNAM,!OB3,C,OUTLET) SPEDAT(SET,OBJTYP,!OB3,C,OUTLET) SPEDAT(SET,ARATIO,!OB3,R,1.) SPEDAT(SET,OBJNAM,!OB4,C,CHAO) SPEDAT(SET,OBJTYP,!OB4,C,PLATE) SPEDAT(SET,FACETDAT,NUMOBJ,I,4) SPEDAT(SET,MATERIAL,198,L,T) ************************************************************ Group 20. Preliminary Printout ************************************************************ 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(KE,Y,N,Y,Y,Y,Y) OUTPUT(EP,Y,N,Y,Y,Y,Y) OUTPUT(PRPS,Y,N,Y,N,N,N) OUTPUT(SHRY,Y,N,Y,N,N,N) OUTPUT(SHRX,Y,N,Y,N,N,N) OUTPUT(EPKE,Y,N,Y,N,N,N) OUTPUT(DEN1,Y,N,Y,N,N,N) OUTPUT(EL1,Y,N,Y,N,N,N) OUTPUT(ENUT,Y,N,Y,N,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 69 ;IYMON = 19 ;IZMON = 1 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = 10001 UWATCH = F ;USTEER = F HIGHLO = F ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 15000 ;NUMCLS = 5 NXPRIN = 22 ;IXPRF = 1 ;IXPRL = 110 NYPRIN = 11 ;IYPRF = 1 ;IYPRL = 56 IPLTF = 1 ;IPLTL = 15000 ;NPLT = 750 ISWPRF = 1 ;ISWPRL = 100000 ITABL = 3 ;IPROF = 1 ABSIZ =0.5 ;ORSIZ =0.4

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 6.667E-02 2.000E-01 3.333E-01 4.667E-01 6.000E-01 7.333E-01 8.667E-01 1.000E+00 1.133E+00 1.267E+00 1.400E+00 1.533E+00 1.667E+00 1.800E+00 1.933E+00 2.067E+00 2.200E+00 2.333E+00 2.467E+00 2.600E+00 2.733E+00 2.867E+00 3.000E+00 3.133E+00 3.267E+00 3.400E+00 3.533E+00 3.667E+00 3.800E+00 3.933E+00 4.026E+00 4.079E+00 4.131E+00 4.184E+00 4.236E+00 4.289E+00 4.341E+00 4.394E+00 4.446E+00 4.499E+00 4.551E+00 4.604E+00 4.656E+00 4.709E+00 4.761E+00 4.814E+00 4.866E+00 4.919E+00 4.971E+00 5.024E+00 5.076E+00 5.129E+00 5.181E+00 5.234E+00 5.286E+00 5.339E+00 5.391E+00 5.444E+00 5.496E+00 5.549E+00 5.719E+00 6.008E+00 6.296E+00 6.585E+00 6.873E+00 7.162E+00 7.450E+00 7.739E+00 8.027E+00 8.316E+00 8.604E+00 8.893E+00 9.181E+00 9.470E+00 9.758E+00 1.005E+01 1.034E+01 1.062E+01 1.091E+01 1.120E+01 1.149E+01 1.178E+01 1.207E+01 1.235E+01 1.264E+01 1.293E+01 1.322E+01 1.351E+01 1.380E+01 1.409E+01 1.437E+01 1.466E+01 1.495E+01 1.524E+01 1.553E+01 1.582E+01 1.611E+01 1.639E+01 1.668E+01 1.697E+01 1.726E+01 1.755E+01 1.784E+01 1.812E+01 1.841E+01 1.870E+01 1.899E+01 1.928E+01 1.957E+01 1.986E+01 Y-coordinates of the cell centres 1.500E-02 3.243E-02 3.777E-02 4.419E-02 5.188E-02 6.112E-02 7.220E-02 8.550E-02 1.015E-01 1.206E-01 1.436E-01 1.712E-01 2.043E-01 2.440E-01 2.916E-01 3.488E-01 4.112E-01 4.684E-01 5.160E-01 5.557E-01 5.888E-01 6.164E-01 6.394E-01 6.585E-01 6.745E-01 6.878E-01 6.989E-01 7.081E-01 7.158E-01 7.222E-01 7.276E-01 7.345E-01 7.445E-01 7.564E-01 7.707E-01 7.879E-01 8.086E-01 8.333E-01 8.631E-01 8.987E-01 9.415E-01 9.929E-01 1.054E+00 1.128E+00 1.217E+00 1.324E+00 1.451E+00 1.605E+00 1.789E+00 2.009E+00 2.274E+00 2.592E+00 2.974E+00 3.432E+00 3.981E+00 4.640E+00 Z-coordinates of the cell centres 5.000E-02 X-coordinates of the (higher) cell faces 1.333E-01 2.667E-01 4.000E-01 5.333E-01 6.667E-01 8.000E-01 9.333E-01 1.067E+00 1.200E+00 1.333E+00 1.467E+00 1.600E+00 1.733E+00 1.867E+00 2.000E+00 2.133E+00 2.267E+00 2.400E+00 2.533E+00 2.667E+00 2.800E+00 2.933E+00 3.067E+00 3.200E+00 3.333E+00 3.467E+00 3.600E+00 3.733E+00 3.867E+00 4.000E+00 4.053E+00 4.105E+00 4.157E+00 4.210E+00 4.263E+00 4.315E+00 4.367E+00 4.420E+00 4.472E+00 4.525E+00

4.577E+00 4.630E+00 4.682E+00 4.735E+00 4.787E+00 4.840E+00 4.892E+00 4.945E+00 4.997E+00 5.050E+00 5.103E+00 5.155E+00 5.207E+00 5.260E+00 5.312E+00 5.365E+00 5.418E+00 5.470E+00 5.523E+00 5.575E+00 5.864E+00 6.152E+00 6.441E+00 6.729E+00 7.017E+00 7.306E+00 7.595E+00 7.883E+00 8.172E+00 8.460E+00 8.749E+00 9.037E+00 9.326E+00 9.614E+00 9.903E+00 1.019E+01 1.048E+01 1.077E+01 1.106E+01 1.135E+01 1.163E+01 1.192E+01 1.221E+01 1.250E+01 1.279E+01 1.308E+01 1.336E+01 1.365E+01 1.394E+01 1.423E+01 1.452E+01 1.481E+01 1.510E+01 1.538E+01 1.567E+01 1.596E+01 1.625E+01 1.654E+01 1.683E+01 1.712E+01 1.740E+01 1.769E+01 1.798E+01 1.827E+01 1.856E+01 1.885E+01 1.913E+01 1.942E+01 1.971E+01 2.000E+01 Y-coordinates of the (higher) cell faces 3.000E-02 3.486E-02 4.069E-02 4.769E-02 5.608E-02 6.616E-02 7.825E-02 9.276E-02 1.102E-01 1.311E-01 1.561E-01 1.862E-01 2.223E-01 2.656E-01 3.176E-01 3.800E-01 4.424E-01 4.944E-01 5.377E-01 5.738E-01 6.039E-01 6.289E-01 6.498E-01 6.672E-01 6.818E-01 6.938E-01 7.039E-01 7.123E-01 7.193E-01 7.251E-01 7.300E-01 7.390E-01 7.499E-01 7.629E-01 7.786E-01 7.973E-01 8.198E-01 8.468E-01 8.793E-01 9.182E-01 9.648E-01 1.021E+00 1.088E+00 1.169E+00 1.266E+00 1.382E+00 1.521E+00 1.688E+00 1.889E+00 2.130E+00 2.419E+00 2.766E+00 3.182E+00 3.681E+00 4.281E+00 5.000E+00 Z-coordinates of the (higher) cell faces 1.000E-01 Total number of F-array elements used is 748890 --- INTEGRATION OF EQUATIONS BEGINS --- Flow field at ITHYD= 1, IZ= 1, ISWEEP= 15000, ISTEP= 1 Field Values of P1 IY= 56 3.384E+01 1.739E+01 -1.111E+01 -4.918E+01 -2.268E+01 IY= 45 5.105E+01 6.225E+01 -1.041E+01 -1.145E+02 -1.379E+01 IY= 34 5.309E+01 7.774E+01 2.885E+01 -1.191E+02 -1.257E+01 IY= 23 5.348E+01 8.167E+01 4.597E+01 -1.191E+02 -1.237E+01 IY= 12 5.442E+01 9.324E+01 none -1.184E+02 -1.159E+01 IY= 1 5.449E+01 9.403E+01 1.377E+02 -1.185E+02 -1.115E+01 IX= 1 23 45 67 89 Field Values of U1 IY= 56 1.667E+01 1.756E+01 1.885E+01 2.053E+01 1.926E+01 IY= 45 1.667E+01 1.598E+01 1.878E+01 1.640E+01 1.544E+01 IY= 34 1.667E+01 1.515E+01 1.492E+01 1.081E+00 9.371E+00 IY= 23 1.667E+01 1.493E+01 1.188E+01 -8.926E-01 8.108E+00 IY= 12 1.667E+01 1.423E+01 none -3.373E+00 4.826E+00 IY= 1 1.667E+01 1.418E+01 1.314E+01 1.912E+01 7.822E+00 IX= 1 23 45 67 89 Field Values of V1 IY= 55 1.102E-02 4.265E-01 5.949E-01 2.063E-01 -1.949E-01 IY= 44 3.411E-02 1.941E+00 6.493E+00 6.657E-03 -1.048E+00

IY= 33 2.394E-02 1.646E+00 7.888E+00 1.029E-01 -8.273E-01 IY= 22 2.045E-02 1.482E+00 6.006E+00 3.132E-01 -7.081E-01 IY= 11 5.239E-03 4.417E-01 none 5.339E-01 -1.081E-01 IX= 1 23 45 67 89 Field Values of KE IY= 56 6.917E-01 6.414E-01 6.383E-01 6.282E-01 5.601E-01 IY= 45 6.918E-01 6.794E-01 1.507E+00 1.210E+01 1.225E+01 IY= 34 6.917E-01 7.069E-01 2.303E+00 1.908E+01 1.828E+01 IY= 23 6.917E-01 7.129E-01 1.806E+00 1.540E+01 1.748E+01 IY= 12 6.898E-01 5.485E-01 none 5.531E+00 9.130E+00 IY= 1 1.472E-05 6.015E-02 1.672E-01 7.879E-02 6.192E-01 IX= 1 23 45 67 89 Field Values of EP IY= 56 3.775E-01 3.243E-01 3.148E-01 2.977E-01 2.360E-01 IY= 45 3.775E-01 3.516E-01 1.448E+00 1.100E+02 4.412E+01 IY= 34 3.775E-01 3.728E-01 1.747E+01 1.070E+02 6.284E+01 IY= 23 3.775E-01 3.786E-01 2.112E+01 6.905E+01 5.701E+01 IY= 12 3.770E-01 3.144E-01 none 3.820E+01 2.734E+01 IY= 1 1.508E-06 3.941E-01 1.827E+00 5.908E-01 1.302E+01 IX= 1 23 45 67 89 Field Values of PRPS IY= 56 pil prop pil prop pil prop pil prop pil prop IY= 45 pil prop pil prop pil prop pil prop pil prop IY= 34 pil prop pil prop pil prop pil prop pil prop IY= 23 pil prop pil prop pil prop pil prop pil prop IY= 12 pil prop pil prop blockage pil prop pil prop IY= 1 pil prop pil prop pil prop pil prop pil prop IX= 1 23 45 67 89 Field Values of SHRY IY= 56 1.000E-10 1.000E-10 1.000E-10 1.000E-10 1.000E-10 IY= 45 1.000E-10 1.000E-10 1.000E-10 1.000E-10 1.000E-10 IY= 34 1.000E-10 1.000E-10 1.000E-10 1.000E-10 1.000E-10 IY= 23 1.000E-10 1.000E-10 1.000E-10 1.000E-10 1.000E-10 IY= 12 1.000E-10 1.000E-10 none 1.000E-10 1.000E-10 IY= 1 1.000E-10 1.000E-10 1.000E-10 1.000E-10 1.000E-10 IX= 1 23 45 67 89 Field Values of SHRX IY= 56 1.000E-10 1.000E-10 1.000E-10 1.000E-10 1.000E-10 IY= 45 1.000E-10 1.000E-10 1.000E-10 1.000E-10 1.000E-10 IY= 34 1.000E-10 1.000E-10 1.000E-10 1.000E-10 1.000E-10 IY= 23 1.000E-10 1.000E-10 1.000E-10 1.000E-10 1.000E-10 IY= 12 1.000E-10 1.000E-10 none 1.000E-10 1.000E-10 IY= 1 1.000E-10 1.000E-10 1.000E-10 1.000E-10 1.000E-10 IX= 1 23 45 67 89 Field Values of EPKE IY= 56 5.457E-01 5.056E-01 4.932E-01 4.739E-01 4.213E-01 IY= 45 5.457E-01 5.175E-01 9.614E-01 9.088E+00 3.601E+00 IY= 34 5.457E-01 5.273E-01 7.584E+00 5.610E+00 3.438E+00 IY= 23 5.457E-01 5.310E-01 1.170E+01 4.483E+00 3.261E+00 IY= 12 5.465E-01 5.733E-01 none 6.905E+00 2.995E+00 IY= 1 1.025E-01 6.552E+00 1.093E+01 7.499E+00 2.102E+01

IX= 1 23 45 67 89 Field Values of DEN1 IY= 56 1.189E+00 1.189E+00 1.189E+00 1.189E+00 1.189E+00 IY= 45 1.189E+00 1.189E+00 1.189E+00 1.189E+00 1.189E+00 IY= 34 1.189E+00 1.189E+00 1.189E+00 1.189E+00 1.189E+00 IY= 23 1.189E+00 1.189E+00 1.189E+00 1.189E+00 1.189E+00 IY= 12 1.189E+00 1.189E+00 none 1.189E+00 1.189E+00 IY= 1 1.189E+00 1.189E+00 1.189E+00 1.189E+00 1.189E+00 IX= 1 23 45 67 89 Field Values of EL1 IY= 56 2.504E-01 2.603E-01 2.661E-01 2.748E-01 2.918E-01 IY= 45 2.504E-01 2.617E-01 2.098E-01 6.289E-02 1.597E-01 IY= 34 2.504E-01 2.620E-01 3.288E-02 1.279E-01 2.043E-01 IY= 23 2.504E-01 2.612E-01 1.888E-02 1.439E-01 2.106E-01 IY= 12 2.497E-01 2.123E-01 none 5.596E-02 1.658E-01 IY= 1 6.150E-03 6.150E-03 6.150E-03 6.150E-03 6.150E-03 IX= 1 23 45 67 89 Field Values of ENUT IY= 56 1.141E-01 1.142E-01 1.165E-01 1.193E-01 1.196E-01 IY= 45 1.141E-01 1.181E-01 1.411E-01 1.199E-01 3.062E-01 IY= 34 1.141E-01 1.207E-01 2.733E-02 3.061E-01 4.784E-01 IY= 23 1.141E-01 1.208E-01 1.390E-02 3.093E-01 4.824E-01 IY= 12 1.136E-01 8.611E-02 none 7.209E-02 2.744E-01 IY= 1 1.292E-05 8.262E-04 1.378E-03 9.456E-04 2.651E-03 IX= 1 23 45 67 89 ************************************************************ Whole-field residuals before solution with resref values determined by EARTH & resfac=1.0E-03 variable resref (res sum)/resref (res sum) P1 1.953E-04 2.315E-01 4.520E-05 U1 3.523E-03 8.907E-01 3.138E-03 V1 2.690E-04 1.979E+00 5.323E-04 KE 4.675E-04 2.594E+00 1.212E-03 EP 2.273E-03 1.605E+00 3.647E-03 ************************************************************ Sources and sinks !! Zero nett sources are not printed !!! Nett Sources have units of mass_per_unit_time * variable Average values have units of the variable Typically the units of the sources are: U1,V1,W1 - Force - Newtons R1 - Mass - kg/s TEM1 - Energy - Watts Nett source of U1 at patch named: OB2 (INLET ) = 1.652050E+02 Nett source of U1 at patch named: OB3 (OUTLET ) =-1.710171E+02 Nett source of U1 at patch named: OB4 (CHAO ) = 1.797176E-01 pos. sum=165.384674 neg. sum=-171.01712 nett sum=-5.632446

Nett source of V1 at patch named: OB3 (OUTLET ) = 2.056486E+00 pos. sum=2.056486 neg. sum=0. nett sum=2.056486 Nett source of R1 at patch named: OB2 (INLET ) = 9.910315E+00 Nett source of R1 at patch named: OB3 (OUTLET ) =-9.910315E+00 (Mass Out -9.910315E+00 In 0.000000E+00) pos. sum=9.910315 neg. sum=-9.910315 nett sum=0. Nett source of KE at patch named: OB2 (INLET ) = 6.884917E+00 (Average 6.947224E-01) Nett source of KE at patch named: OB3 (OUTLET ) =-2.307733E+01 (Ave Out 2.328618E+00 In 0.000000E+00) Nett source of KE at patch named: KESOURCE =-1.546713E+02 pos. sum=6.884917 neg. sum=-177.748657 nett sum=-170.863739 Nett source of EP at patch named: OB2 (INLET ) = 3.771889E+00 (Average 3.806024E-01) Nett source of EP at patch named: OB3 (OUTLET ) =-3.922397E+01 (Ave Out 3.957894E+00 In 0.000000E+00) Nett source of EP at patch named: KESOURCE =-3.063000E+03 pos. sum=3.771889 neg. sum=-3102.224121 nett sum=-3098.452148 Integrated force on object: FORMULA ----------------------------------- Fx = 1.572264E+01 (Pressure= 1.561001E+01, Friction= 1.126203E-01) Fy = 9.973325E+00 (Pressure= 9.961365E+00, Friction= 1.196007E-02) Force unit vector: 8.444388E-01 5.356520E-01 1.380849E-04 Total moment about Z axis = -3.892967E+01 Moment of Fx about Z axis = 5.323749E+00 at distance Y = 3.386041E-01 Moment of Fy about Z axis = -4.425343E+01 at distance X = 4.437179E+00 Integrated forces for all included objects ------------------------------------------ Fx = 1.572264E+01 Fy = 9.973325E+00 Ftot= 1.861904E+01 Force unit vector: 8.444388E-01 5.356520E-01 1.380849E-04 The total force acts at (centre of pressure): X = 4.437179E+00 Y = 3.386041E-01 Z = 5.000000E-02 Moments about origin: Total moment about Z axis = -3.892967E+01 Moment of Fx about Z axis = 5.323749E+00 at distance Y = 3.386041E-01 Moment of Fy about Z axis = -4.425343E+01 at distance X = 4.437179E+00

************************************************************ spot values vs sweep or iteration number IXMON= 69 IYMON= 19 IZMON= 1 TIMESTEP= 1 Tabulation of abscissa and ordinates... ISWP P1 U1 V1 KE EP 1 6.5276E+03 1.5053E+01 -2.6052E-01 6.9616E-01 5.3584E+00 751 -1.1890E+02 -1.2844E+00 1.6302E+00 1.7606E+01 7.2142E+01 1501 -1.1933E+02 -1.2388E+00 4.0167E-01 1.5878E+01 6.1833E+01 2251 -1.2051E+02 -1.3461E+00 2.6458E-01 1.5232E+01 5.8211E+01 3001 -1.1901E+02 -1.4083E+00 2.4665E-01 1.4913E+01 5.6399E+01 3751 -1.1681E+02 -1.4426E+00 2.4552E-01 1.4740E+01 5.5410E+01 4501 -1.1490E+02 -1.4616E+00 2.4581E-01 1.4644E+01 5.4867E+01 5251 -1.1353E+02 -1.4722E+00 2.4601E-01 1.4591E+01 5.4565E+01 6001 -1.1263E+02 -1.4782E+00 2.4611E-01 1.4561E+01 5.4396E+01 6751 -1.1206E+02 -1.4816E+00 2.4618E-01 1.4545E+01 5.4301E+01 7501 -1.1172E+02 -1.4835E+00 2.4623E-01 1.4535E+01 5.4247E+01 8251 -1.1152E+02 -1.4846E+00 2.4625E-01 1.4530E+01 5.4216E+01 9001 -1.1140E+02 -1.4853E+00 2.4627E-01 1.4527E+01 5.4198E+01 9751 -1.1133E+02 -1.4857E+00 2.4628E-01 1.4525E+01 5.4188E+01 10501 -1.1129E+02 -1.4859E+00 2.4628E-01 1.4524E+01 5.4182E+01 11251 -1.1127E+02 -1.4860E+00 2.4629E-01 1.4523E+01 5.4178E+01 12001 -1.1125E+02 -1.4861E+00 2.4629E-01 1.4523E+01 5.4176E+01 12751 -1.1124E+02 -1.4861E+00 2.4629E-01 1.4523E+01 5.4175E+01 13501 -1.1124E+02 -1.4861E+00 2.4629E-01 1.4522E+01 5.4175E+01 14251 -1.1124E+02 -1.4861E+00 2.4629E-01 1.4522E+01 5.4174E+01 Variable 1 = P1 2 = U1 3 = V1 4 = KE 5 = EP Minval= -1.205E+02 -1.486E+00 -2.605E-01 6.962E-01 5.358E+00 Maxval= 6.528E+03 1.505E+01 1.630E+00 1.761E+01 7.214E+01 Cellav= 2.183E+02 -6.208E-01 2.987E-01 1.413E+01 5.347E+01 1.00 2..5.+....+....+....+....+....+....+....+....+....+ . . 0.90 + 4 + . 5 4 4 4 . 0.80 + 5 4 4 4 4 4 4 4 4 4 4 4 4 4 4 . 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 0.70 + + . . 0.60 + + . . 0.50 + + . . 0.40 + + . 3 . 0.30 + 3 + . 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 0.20 + + . . 0.10 + + . . 0.00 5..2.2..2.+2.2.+2.2.+2..2+2..2+.2.2+.2.2+.2..2.2..2 0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 the abscissa is ISWP. min= 1.00E+00 max= 1.43E+04 ************************************************************

************************************************************ residuals vs sweep or iteration number Tabulation of abscissa and ordinates... ISWP P1 U1 V1 KE EP 1 5.1816E+06 6.1600E+06 4.3044E+06 4.8214E+06 6.0518E+06 751 1.2331E+01 5.4861E+03 2.1419E+04 3.3966E+03 3.5567E+03 1501 4.5686E+00 2.0047E+03 4.8611E+03 2.4208E+02 6.6243E+02 2251 1.8625E+00 1.0039E+03 2.0831E+03 1.1466E+02 2.8683E+02 3001 1.5777E+00 6.4650E+02 9.6973E+02 6.3663E+01 1.4516E+02 3751 1.3965E+00 4.5772E+02 5.1185E+02 3.9707E+01 8.1565E+01 4501 1.0876E+00 2.9859E+02 3.3425E+02 2.5135E+01 4.6761E+01 5251 8.4171E-01 1.8551E+02 2.2390E+02 1.6345E+01 2.7794E+01 6001 6.4076E-01 1.1213E+02 1.4081E+02 1.0818E+01 1.6695E+01 6751 5.1467E-01 6.6883E+01 8.5944E+01 7.4924E+00 1.0418E+01 7501 4.2697E-01 3.9525E+01 5.1734E+01 5.8101E+00 6.9552E+00 8251 3.6635E-01 2.3241E+01 3.0843E+01 4.8078E+00 4.6903E+00 9001 3.3551E-01 1.3698E+01 1.8412E+01 3.8426E+00 3.5792E+00 9751 3.0472E-01 8.1333E+00 1.1160E+01 3.7246E+00 2.8622E+00 10501 2.8642E-01 4.8850E+00 7.0331E+00 3.2368E+00 2.3307E+00 11251 2.8396E-01 3.0419E+00 4.6839E+00 2.8465E+00 1.9530E+00 12001 2.6021E-01 1.9621E+00 3.2648E+00 2.7101E+00 1.8882E+00 12751 2.4695E-01 1.3565E+00 2.5402E+00 2.6052E+00 1.7628E+00 13501 2.5308E-01 1.0157E+00 2.3528E+00 2.4805E+00 1.7033E+00 14251 2.2043E-01 9.0812E-01 2.0338E+00 2.5981E+00 1.6911E+00 Variable 1 = P1 2 = U1 3 = V1 4 = KE 5 = EP Minval= -1.512E+00 -9.638E-02 7.099E-01 9.085E-01 5.254E-01 Maxval= 1.546E+01 1.563E+01 1.528E+01 1.539E+01 1.562E+01 1.00 5....+....+....+....+....+....+....+....+....+....+ . . 0.90 + + . . 0.80 + + . . 0.70 + + . 3 . 0.60 + + . 2 3 . 0.50 + 5 2 3 + . 2 . 0.40 + 5 3 3 + . 5 3 2 . 0.30 + 4 5 3 3 + . 1 4 5 3 2 . 0.20 + 1 4 4 5 5 3 3 + . 1 4 4 5 3 2 . 0.10 + 1 1 1 1 4 5 5 3 3 2 + . 1 1 4 5 5 5 3 3 2 . 0.00 +....+....+....+....+....+....+.1.1+.5.5+.5..5.5..5 0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 the abscissa is ISWP. min= 1.00E+00 max= 1.43E+04 ************************************************************ ************************************************************ SATLIT RUN NUMBER = 1 ; LIBRARY REF.= 0 Run started at 17:51:05 on Wednesday, 30 May 2012

Run completed at 18:02:28 on Wednesday, 30 May 2012 CPU time of run 616 s This includes 610 seconds of user time and 5 seconds of system time. TIME/(VARIABLES*CELLS*TSTEPS*SWEEPS*ITS) = 1.333E-06 ************************************************************