Calamidades naturais são cada dia mais frequentes e a maior parte delas está relacionada aos danos ocasionados por tempestades severas. Desta forma, são necessários estudos que auxiliam tanto na compreensão desses fenômenos quanto na previsão dos mesmos para que sejam desenvolvidas estratégias mitigatórias e edificações cada vez mais preparadas para as adversidades meteorológicas. Intensas rajadas de vento provocam um dos principais danos provenientes de tempestades severas e podem ser originadas a partir de fenômenos tais como downbursts ou tornados. Downbursts são fenômenos meteorológicos que apresentam distinto perfil de velocidade de vento e turbulência que causam diversos acidentes no mundo e, por meio deste trabalho, verificou-se que são a causa de diversos acidentes também no Brasil, em especial na região da Bacia Amazônica, na Região Sul e Sudeste brasileira e na região próxima ao litoral do nordeste brasileiro. A previsibilidade de downbursts também é abordada nessa pesquisa; buscou-se verificar a relação existente entre altos valores do índice CAPE e baixos valores de CINE na ocorrência de ventos provenientes de downbursts, tendo sido constatado que esses índices não são eficientes se analisados de forma isolada na previsão de downbursts. Por fim, verificou-se que estudos mais aprofundados a respeito das ações de downbursts sobre edificações fazem-se necessários para o Brasil. Por fim desenvolveu-se uma revisão de literatura a respeito da simulação física e numérica do fenômeno, buscando gerar um embasamento teórico para futuros trabalhos nesta área. / Natural disasters are becoming more frequent every day and most of them are related to damages caused by severe storms, thus studies that help both in understanding these phenomena and in predicting are necessary to the development of mitigation strategies and to the development of building increasingly prepared to the weather adversity. Intense wind gusts are responsible for the major damage from severe thunderstorms and can arise from phenomena such as downbursts or tornadoes. Downbursts are meteorological phenomena that exhibit distinctive profiles of wind speed and turbulence causing many accidents in the world and through this work are found to be the cause of many accidents in Brazil too, especially in the Amazon Basin region in the South and Southeast Brazilian’s Regions and in the region next to the coast of Brazilian’s Northeastern. Predictability of downbursts is also addressed in this study, it was valued the relationship between high values of CAPE index and low values of CINE index for the occurrence of winds from downbursts, it was found that these index are not efficient when analyzed in isolation in predicting downbursts. Finally, it was concluded that more intensive studies on the effects of downbursts on structures are up needed in Brazil. Finally it was developed an investigation on literature about the physical and numerical simulation of this phenomena, seeking to generate a theoretical basis for future work in this area.

Desde a exploração até o planejamento de lavra, metodologias para quantificação do risco associado ao empreendimento mineiro podem melhorar drasticamente o processo de tomada de decisão. Pequenas variações nas condições de contorno de um projeto podem ter impacto significativo sobre o retorno final do mesmo. Dentro deste contexto, na mineração, a incerteza geológica é vista como um dos principais fatores que contribuem para o insucesso do projeto. A medição e gerenciamento do risco geológico na avaliação do projeto e tomada de decisão pode ser realizada em varias etapas do empreendimento mineiro, desde a quantificação dos recursos, sequenciamento de lavra até a exaustão das reservas. Além disso, a incerteza geológica pode ser traduzida no risco financeiro associado, proporcionando uma análise sobre a viabilidade econômica do empreendimento. A necessidade da quantificação da incerteza geológica já era reconhecida pela indústria de mineração desde a década de 1970, onde estimativas globais e locais eram consideradas insuficientes para a otimização do planejamento de produção, sequenciamento de lavra ou estratégias de homogeneização de pilhas. Desta forma, a simulação estocástica vem desempenhando um papel crucial na construção de modelos de incerteza de teores em depósitos minerais, proporcionando uma ferramenta para obter uma análise do risco associado. Este trabalho tem como objetivo o desenvolvimento de uma metodologia que contribua para a construção de cenários de planejamento de lavra e avaliação do impacto provocado pelas incertezas geológicas fornecidas a partir de modelos simulados. As respostas buscadas, usando a abordagem sugerida, compreendem a definição de uma metodologia para controle e planejamento de lavra incorporando a incerteza local associada ao minério, avaliando a sensibilidade do retorno financeiro e o impacto no planejamento de lavra, considerando variações e incertezas relacionadas ao modelo geológico e estimativa de teores. / From exploration to mine planning, methodologies for quantifying the risk associated with mining enterprise can dramatically improve decision making process. Small variations in the boundary conditions of a project may have significant impact on the final return of the same. Within this context, in mining, geological uncertainty is seen as a major factor contributing to the failure of the project. The need to quantify and manage the geological risk for project evaluation and decision making can be translated into the need to assess the geological risk of any parameter at all stages, from the resources quantification, mining sequencing to reserves exhaustion. In addition, the geological uncertainty can be translated into financial risk, providing an analysis of the economic feasibility of the project. However, the need for quantification of geological uncertainty was already recognized by the mining industry since the 1970, where local and global estimates were considered insufficient for the optimization of production planning, sequencing mining, strategies of homogenization. Thus, the stochastic simulation has played a crucial role in the construction of grade uncertainty models in mineral deposits, providing a tool to perform risk analysis. This study aims to develop a methodology able to contribute in building scenarios of mine planning and evaluation of the impact caused by geological uncertainties supplied from simulated models. Answers searched using the suggested approach include understand the definition of a methodology to control and mine planning incorporating the local uncertainty associated with the ore, assessing the sensitivity of financial return and mine planning considering variantions and uncertainty related to the geological model and grade estimation.

Este trabalho tem por objetivo simular e analisar uma usina termelétrica a carvão em várias condições de funcionamento. A usina simulada neste trabalho é a AVV 1 localizada em Copenhague, Dinamarca. A AVV 1 é uma usina de geração de potência e aquecimento distrital que pode funcionar em várias condições de carga. A simulação da usina supracitada foi tema de um concurso de simuladores proposto no congresso ECOS 2003 realizado em Copenhague, Dinamarca. Para a realização deste trabalho foi construído um programa na linguagem FORTRAN 90. Cada componente da usina é modelado através de equações de balanço de massa e energia, e o sistema completo tem sua solução obtida pelo método de substituição sucessiva. Para viabilizar essa solução é necessário também implementar uma rotina de cálculo de propriedades do fluido de trabalho. No caso estudado, o fluido de trabalho da usina é a água e a formulação utilizada para o cálculo de suas propriedades nos diversos estados é a IAPWS IF-97. A usina é simulada em dois modos de operação: modo de condensação, onde ocorre apenas geração de eletricidade, e em modo de contrapressão, onde há geração de eletricidade e aquecimento distrital, conforme nomenclatura sugerida pela organização do concurso No modo de condensação, são feitas quatro séries de simulações variando a carga de 100% a 40%. Cada série contém um conjunto de hipóteses quanto à variação das eficiências isoentrópicas e pressões das turbinas em função da vazão mássica. No modo de contrapressão, a usina é simulada funcionando com 100% da carga. O programa desenvolvido calcula as propriedades para qualquer ponto de trabalho ao longo da planta, assim como a eficiência da mesma, a potência gerada, e todas as vazões mássicas pertinentes. Além disso, é feita também uma análise exergética da planta. A simulação demonstrou que a planta possui uma eficiência global de 42,02% com uma geração de 250,2 MW em 100% de carga no modo de condensação. Nessas mesmas condições, do ponto de vista exergético, a eficiência encontrada é de 37,21%. No modo de contrapressão, a usina apresenta uma eficiência exergética de 40,19% com um aproveitamento energético de 90,55%. Por fim, é possível também verificar a comportamento da eficiência da planta e a variação de água de resfriamento do condensador com a carga. Os resultados gerados são próximos àqueles encontrados pelos diversos pesquisadores que abordaram o problema.

Este trabalho apresenta uma sistemática para realizar a otimização numérica de pré-formas e de matrizes em problemas de forjamento axissimétricos e em estado plano de deformações. Para este fim, desenvolveu-se um código computacional composto basicamente de três módulos: módulo de pré-processamento, módulo de análise e módulo de otimização. Cada um destes foi elaborado acrescentando rotinas em programas comerciais ou acadêmicos disponíveis no GMAp e no CEMACOM. Um programa gerenciador foi desenvolvido para controlar os módulos citados no processo de otimização. A abordagem proposta apresenta uma nova função objetivo a minimizar, a qual está baseada em uma operação booleana XOR (exclusive or) sobre os dois polígonos planos que representam a geometria desejada para o componente e a obtida na simulação, respectivamente. Esta abordagem visa eliminar possíveis problemas geométricos associados com as funções objetivo comumente utilizadas em pesquisas correlatas. O trabalho emprega análise de sensibilidade numérica, via método das diferenças finitas. As dificuldades associadas a esta técnica são estudadas e dois pontos são identificados como limitadores da abordagem para problemas de conformação mecânica (grandes deformações elastoplásticas com contato friccional): baixa eficiência e contaminação dos gradientes na presença de remalhamentos. Um novo procedimento de diferenças finitas é desenvolvido, o qual elimina as dificuldades citadas, possibilitando a sua aplicação em problemas quaisquer, com características competitivas com as da abordagem analítica Malhas não estruturadas são tratadas mediante suavizações Laplacianas, mantendo as suas topologias. No caso de otimização de pré-formas, o contorno do componente a otimizar é parametrizado por B-Splines cujos pontos de controle são adotados como variáveis de projeto. Por outro lado, no caso de otimização de matrizes, a parametrização é realizada em termos de segmentos de reta e arcos de circunferências. As variáveis de projeto adotadas são, então, as coordenadas das extremidades das retas, os raios e centros dos arcos, etc. A sistemática é fechada pela aplicação dos algoritmos de programação matemática de Krister Svanberg (Método das Assíntotas Móveis Globalmente Convergente) e de Klaus Schittkowski (Programação Quadrática Sequencial – NLPQLP). Resultados numéricos são apresentados mostrando a evolução das implementações adotadas e o ganho de eficiência obtido.

O presente trabalho apresenta uma modelagem detalhada de processos de combustao turbulentos para um jato concentrico de combustıvel e ar. A modelagem é fundamentada nas equacões de conservacão de massa, de quantidade de movimento, de energia e de espécies quımicas. A turbulencia é resolvida pela utilizacão do modelo k- padrão. Dois modelos de reacões quımicas são apresentados. O modelo SCRS – Simple Chemically-Reacting Systems, que assume taxas instantâneas de reacões quımicas. Também é abordado o modelo E-A – Eddy Breakup - Arrhenius, que assume taxas finitas de reacões quımicas. A radiacão térmica, fenômeno de grande importância devido as altas temperaturas alcancadas em processos de combustão,é modelada através do Método das Zonas. O modelo da soma ponderada de gases cinzas – WSGGM,é usado para determinar o espectro de emissão e absorcão dos gases no processo. Para a solucão destas equacões diferenciais, juntamente com os modelos de turbulência, de reaçõoes químicas e radiação térmica, faz-se o uso do Método dos Volumes Finitos. Para validar a modelagem apresentada resolve-se o processo de combustão em uma câmara cilíndrica. A câmara de combustão usada áa mesma abordada no First Workshop on Aerodynamics of Steady State Combustion Chambers and Furnaces, organizado pela ERCORTAC - European Research Community On Flow Turbulence And Combustion, em outubro de 1994, que apresenta dados experimentais de temperatura e concentração das espécies químicas para várias posições de interesse no interior da câmara. Utiliza-se o gás natural como combustível e o ar atmosférico como oxidante. O processo de combustão sem pré-mistura é resolvido para a condição de excesso de combustível de 5 % para ambos os modelos, onde o gás natural é injetado por um duto circular central, e o ar atmosférico por um orifício anular externo a esse duto, no mesmo plano Uma reação química não estagiada é assumida para o modelo SCRS. Para o modelo E-A duas situações são resolvidas: combustão não estagiada, com uma etapa global de reação química; e reação quımica estagiada, com duas etapas globais. Os resultados obtidos com o modelo SCRS para a distribuição de temperaturas, em termos de tendências gerais, são razoáveis. Já as concentrações de espécies químicas não apresentam dados satisfatórios para este modelo. Para o modelo E-A os resultados apresentam boa concordância com os dados experimentais, principalmente para a situação em que o processo de combustão é assumido em duas etapas globais. ´E analisado em detalhe o papel desempenhado pela transferencia de calor por radiacao, com meio participante. Para melhor verificar as trocas de calor, assume-se uma camara de combustao cilındrica com paredes d’agua. A injecao do combustıvel e do oxidante e feita atraves de um queimador central, semelhante ao usado para validar a modelagem, porem com dois orifıcios concentricos para injecao de combustıvel. Nesta situação o efeito do turbilhonamento (swril), assumido como 20 % da velocidade axial de entrada, sobre a injecao de ar e computado atraves da condicao contorno da equacao de conservacao da quantidade de movimento angular. Nesta fase apenas o modelo E-A, com duas etapas globais de reacoes quımicas, e considerado, ja que o mesmo apresenta os melhores resultados. O processo de combustao e simulado com e sem a presenca da radiacao termica. Verifica-se que a presenca da radiacao termica homogeneiza a temperatura dos gases no interior da camara. Com isso verifica-se tambem alterações nas taxas de reacoes quımicas, modificando a magnitude das fracoes das especies quımicas Quando a radiacao termica e considerada efeitos de extinção local da chama sao verificados nas regioes de temperaturas mais altas, diminuindo o consumo de oxigenio e aumentando a producao de monoxido de carbono, caracterizando assim uma combustao incompleta. Em algumas situacoes tem-se uma variacao de temperatura de ate 500 K, a montante da chama. A radiacao termica tambem aumenta a taxa de transferencia de calor dos gases quentes para as paredes da camara, e desta para o seu exterior. Com os resultados obtidos a partir desta modelagem e possıvel determinar o perfil da zona de combustao, a distribuicao de concentracoes de especies quımicas, o campo de velocidades e as taxas de transferencia de calor para as paredes da camara de combustao, total, por conveccao superficial e por radiacao. Estes resultados sao de extrema importancia para prever a performance de camaras de combustao, assim como auxiliar na sua otimizacao.

Submitted by Eduarda Figueiredo (eduarda.ffigueiredo@ufpe.br) on 2015-03-10T13:05:04Z No. of bitstreams: 2 Tese__EngMec2012_Lira Júnior, José Claudino.pdf: 8591000 bytes, checksum: f9939b3d89bca78a83d391b60adecbc1 (MD5) license_rdf: 1232 bytes, checksum: 66e71c371cc565284e70f40736c94386 (MD5) / Made available in DSpace on 2015-03-10T13:05:05Z (GMT). No. of bitstreams: 2 Tese__EngMec2012_Lira Júnior, José Claudino.pdf: 8591000 bytes, checksum: f9939b3d89bca78a83d391b60adecbc1 (MD5) license_rdf: 1232 bytes, checksum: 66e71c371cc565284e70f40736c94386 (MD5) Previous issue date: 2011-07-20 / Vórtice Engenharia / No processo de laminação de aço, os cilindros que conformam o material sofrem desgastes constantes, devido ao aquecimento continuo provocado pelas barras em processamento. Esse desgaste interfere na vida útil do cilindro. Para manter o cilindro em operação são feitos câmbios para reposicionar o ponto de passagem da barra, porém, esta operação leva à parada total do laminador causando perdas financeiras para a produção. Este trabalho tem como objetivo apresentar um modelo computacional para simular o comportamento térmico dos cilindros de laminação. O modelo foi elaborado utilizando o método dos volumes finitos para um sistema bidimensional, transiente e móvel. O modelo calcula a distribuição de temperatura superficial do cilindro em diversas condições de serviço. Os resultados apresentados auxiliaram na validação do modelo e realização de estudos paramétricos a fim de investigar a influência das variáveis na temperatura circunferencial do cilindro.

GUERRERO, Jorge Recarte Henriquez, também é conhecido em citações bibliográficas por: HENRIQUEZ, Jorge Recarte / Submitted by Pedro Barros (pedro.silvabarros@ufpe.br) on 2018-08-06T18:42:44Z No. of bitstreams: 2 license_rdf: 811 bytes, checksum: e39d27027a6cc9cb039ad269a5db8e34 (MD5) DISSERTAÇÃO João Felipe Vital Gonçalves de Souza.pdf: 1824224 bytes, checksum: ef97947078137d9ff4ee1630a019a15d (MD5) / Rejected by Alice Araujo (alice.caraujo@ufpe.br), reason: Por favor, inserir lattes. on 2018-08-09T21:41:57Z (GMT) / Submitted by Pedro Barros (pedro.silvabarros@ufpe.br) on 2018-08-13T21:08:29Z No. of bitstreams: 2 license_rdf: 811 bytes, checksum: e39d27027a6cc9cb039ad269a5db8e34 (MD5) DISSERTAÇÃO João Felipe Vital Gonçalves de Souza.pdf: 1824224 bytes, checksum: ef97947078137d9ff4ee1630a019a15d (MD5) / Approved for entry into archive by Alice Araujo (alice.caraujo@ufpe.br) on 2018-08-13T21:17:56Z (GMT) No. of bitstreams: 2 license_rdf: 811 bytes, checksum: e39d27027a6cc9cb039ad269a5db8e34 (MD5) DISSERTAÇÃO João Felipe Vital Gonçalves de Souza.pdf: 1824224 bytes, checksum: ef97947078137d9ff4ee1630a019a15d (MD5) / Made available in DSpace on 2018-08-13T21:17:56Z (GMT). No. of bitstreams: 2 license_rdf: 811 bytes, checksum: e39d27027a6cc9cb039ad269a5db8e34 (MD5) DISSERTAÇÃO João Felipe Vital Gonçalves de Souza.pdf: 1824224 bytes, checksum: ef97947078137d9ff4ee1630a019a15d (MD5) Previous issue date: 2016-09-02 / CAPES / O fenômeno de fusão e solidificação está presente na natureza e em diversos processos industriais. Entre outros podem ser citados o crescimento de cristais; o congelamento do solo e lagos em ambientes de baixa temperatura; a fundição e solidificação de metais e ligas; o congelamento e descongelamento de alimentos; o controle térmico em edificações através do uso de paredes, janelas e tetos preenchidos com material de mudança de fase e o armazenamento de energia térmica por calor latente. É nesta última aplicação que o presente estudo tem a sua maior utilidade. A busca por novos materiais que possam ter propriedades vantajosas para o processo de termoacumulação tem voltado à atenção para os nanofluidos que são substâncias formadas por nanopartículas dispersadas num fluido. O presente trabalho pretende contribuir com informações técnicas para o estado da arte no processo de solidificação e nanofluidos por meio da integração destes dois temas num estudo que analisará o processo de transferência de calor com mudança de fase no interior de uma cápsula esférica que contem um nanofluido. O estudo foi realizado por meio de um modelo numérico construído com base no método de diferenças finitas com malha móvel e utilizando a plataforma Matlab para programação do algoritmo de simulação. Foram analisados nanofluidos a base de água e parafina C14 com nanopartículas de oxido de alumínio (Al2O3), óxido de cobre (CuO), nanotubos de carbono e grafeno, e para cada uma destas combinações foram geradas curvas de resfriamento e posição da interface sólido-líquido em função do tempo, assim como os tempos de solidificação completa para cada condição operacional em função da temperatura e vazão do fluido de resfriamento, diâmetro da capsula esférica e também em função do tipo de nanofluido. Os resultados indicam que a redução do tempo de solidificação é maior, dentre as simulações realizadas, quando a combinação é entre água e nanotubos de carbono e parafina C14 e grafeno. / The melting and solidification phenomenon is present in the Nature and some industrial processes. Growing crystals; freezing of lakes at low temperatures; melting and freezing of alloys and metals; freezing and thaw of food; thermal control in buildings through filling of windows, walls and roof with phase change materials; latent heat thermal energy storage. This paper focuses on thislast mentioned application. The research on materials with improved properties for thermal energy storage aims nanofluids, that is, substances composed of nanoparticles dispersed in a base fluid. A priori, the same fluids used like phase change materials can play the role of nanofluids, dispersing nanoparticles with size and proportions propers. The present paper aims contribute with technical informations about art state of solidification process and nanofluids through the integration of these subjects in a study that will analyze the heat transfer process with phase change in a spherical capsule containing nanofluids. The study used a numerical model based on the finite differences method with mobile mesh in the MatLab. Water and paraffine C14 were the base fluids, and alumine oxide, cupper oxide, carbon nanotubes and graphene were the nanoparticles, and all combinantions were analyzed through curves of solid-liquid interface position versus time and the curves of complete solidification time considering temperature and flow rate of cooling fluid and diameter of spherical capsule being constant. The results indicate that the solidification time are reduced when the combinations are water and carbon nanotubes and paraffin C14 and graphene.

Made available in DSpace on 2018-08-02T00:03:15Z (GMT). No. of bitstreams: 1 tese_9683_ALINE CRISTINA SILVA DOS ANJOS.pdf: 3641799 bytes, checksum: f90e5717dc480d7f0db0e2a66b2bfcb0 (MD5) Previous issue date: 2015-11-30 / Um material suscetível a alguma transformação química é considerado um meio poroso reativo. Combustíveis são materiais capazes de liberar energia ao reagirem com o oxigênio e podem ser considerados como um meio poroso reativo. Este trabalho teve o objetivo de modelar a combustão em um meio poroso reativo. O modelo foi implementado em um programa computacional para simular a propagação da frente de combustão em um meio poroso reativo em leito fixo, levando-se em conta mecanismos de reação que podem ocorrer durante o processo. As sub-rotinas que resolvem as equações de conservação da modelagem foram validadas por meio de análises comparativas entre as soluções numéricas e soluções analíticas ou experimentais para garantir que o modelo utilizado fosse capaz de gerar resultados condizentes com a realidade. Além disso, uma investigação foi realizada sobre a influência que o mecanismo de reação, adotado na modelagem, tem sobre o perfil de temperatura na ignição do leito de combustão. Os testes mostraram que o modelo é sensível às mudanças dos mecanismos. Por fim, o modelo aponta a possibilidade de prever, na evolução da temperatura, o fenômeno da ignição. Palavras chave: combustível sólido, meio poroso reativo, simulação numérica.

Made available in DSpace on 2018-08-24T22:53:32Z (GMT). No. of bitstreams: 1 tese_7561_L_C_Jesus_2014_Dissertacao.pdf: 15142064 bytes, checksum: 5d3f671f6ffba1ee71659099f67a5c5e (MD5) Previous issue date: 2014-02-25 / A maré é um dos principais processos físicos que provocam mistura no oceano e também pela variância de uma grande quantidade de variáveis oceanográficas, como salinidade, temperatura, clorofila, velocidades, etc. A previsão de mapas acurados de correntes e elevações e do comportamento energético da maré é indispensável para muitas aplicações práticas e científicas como transporte marítimo, dispersão de poluentes, operações náuticas, a determinação da circulação termohalina e a mistura no oceano entre outros. A maré astronômica barotrópica e baroclínica é estudada neste trabalho. O modelo numérico Delft3D, baseado nas equações de conservação, é utilizado para reconstruir as elevações e velocidades geradas pela maré para a região na região da cadeia Vitória-Trindade, ao largo da costa do Espírito Santo. Os resultados numéricos foram comparados a dados pré-existentes para estações de maré e boia fundeada. As séries temporais de velocidades e elevação foram transformadas para o espectro da frequência através de análise harmônica de Fourier e as equações de conservação de energia foram solucionadas para estimar as taxas de conversão de energia, do fluxo de energia baroclínico e barotrópico e dissipação. Os resultados mostram que as maiores elipses de maré estão localizadas na região do Banco de Abrolhos e Cadeia Vitória-Trindade e que as componentes harmônicas da maré possuem propagação essencialmente progressiva. Na Cadeia Vitória-Trindade foram encontrados as maiores taxas de conversão e dissipação de energia e também de focalização de energia. O modelo se mostrou satisfatório para a simulação das elevações e fases da maré astronômica e das flutuações de massa específica. Os montes submarinos que formam a Cadeia Vitória-Trindade são uma importante feição geomorfológica responsável pela maior parte das transformações de energia da maré.

Durante a realização de projetos de sistemas de refrigeração assistidos por sistemas de armazenagem de energia, na forma latente, é necessário que se saiba de antemão as taxas de transferência de energia dos mesmos. Na literatura soluções para este problema são escassos e não abordam todas as geometrias utilizadas para o armazenamento de energia latente, inexistindo correlaçõesmatemáticas que permitam aos projetistas prever o comportamento da fusão e da solidificação do material de mudança de fase (MMF) e as taxas de transferência de energia térmica. Na solução destes problemas são usualmente adotadas hipóteses simplificativas para amenizar o grau de complexidade da descrição matemática do problema Uma destas hipóteses constitui-se em desconsiderar o processo advectivo na transferência de calor no MMF, em sua fase líquida. Outra hipótese adotada, quando o MMF é a água pura, consiste em utilizar um comportamento linear para a massa específica ou um comportamento médio, consistindo em um afastamento da realidade. Este trabalho descreve esforços no sentido de aproximar estas soluções cada vez mais da realidade fisica do problema. São apresentadas simulações para a mudança de fase da água pura levando em conta a inversão da massa específica, em geometria polar, na forma transiente, utilizando o método de Volumes Finitos para solução das equações governantes. As soluções apresentadas, otimizadas quanto às malhas espacial e temporal, envolvem condições de contorno de primeira e terceira espécies Como resultado das soluções são apresentados o Nusselt local, ao longo da parede do cilindro além do Nusselt médio. Correlações matemáticas para o número de Nusselt médio durante a fusão e para o volume de MMF fundido, são propostos e apresentados. A evolução geométrica da fronteira de mudança de fase também é avaliada para as duas condições de contorno e diversos números de Rayleigh. Comparações entre o volume fundido considerando apenas o processo difusivo e o volume fundido utilizando a convecção são realizadas. Para a verificação da independ~encia de malha aplica-se o índice de convergência de malha (GCI) . Algumas comparações qualitativas com soluções experimentais disponíveis na literatura são apresentadas.