Análise computacional dos campos de velocidade e temperatura do ambiente interno da usina termelétrica Santana - Amapá

A simulação numérica do escoamento de ar em ambientes internos é na atualidade o método mais apropriado para análise de conforto térmico em ambientes internos. O escoamento de ar nesses ambientes configura-se como um escoamento complexo, pois, em regra geral, é uma combinação de escoamentos cisalhantes livres (jatos) e cisalhantes de parede, além disso, esses escoamentos são governados por forças de inércia e forças de empuxo, caracterizando-o como de convecção mista. A combinação desses mecanismos cria um escoamento com características complexas, como zonas de recirculação, vórtices, descolamento e recolamento de camada-limite dentre outras. Portanto, a precisão da solução estará diretamente ligada, principalmente, na habilidade do modelo de turbulência adotado de reproduzir as características turbulentas do escoamento de ar e da transferência térmica. O objetivo principal do presente trabalho foi a simulação computacional do ambiente térmico interno do galpão que abriga os geradores e motores Wärtzilä da Usina Termelétrica Santana no estado do Amapá. A formulação matemática baseada na solução das equações gerais de conservação inclui uma análise dos principais modelos de turbulência aplicados ao escoamento de ar em ambientes internos, assim como os processos de transferência de calor associados. Na modelagem numérica o método de volumes finitos é usado na discretização das equações de conservação, através do código comercial Fluent-Airpak, que foi usado nas simulações computacionais para a análise dos campos de velocidade e temperatura do ar. A utilização correta do programa computacional foi testada e validada para o problema através da simulação precisa de casos retirados da literatura. Os resultados numéricos foram comparados a dados obtidos de medições experimentais realizados no galpão e apresentou boa concordância, considerando a complexidade do problema simulado, o objetivo da simulação em face da diminuição da temperatura no interior do galpão e, também, em função das limitações encontradas quando da tomada das medições experimentais. Além disso, foram feitas simulações de estratégias de melhoria do ambiente térmico da Usina, baseadas na realidade levantada e nos resultados da simulação numérica. Finalmente, foram realizadas simulações do protótipo de solução proposto para a diminuição da temperatura interna do galpão o que possibilitará um aumento, na faixa de 20 a 30%, do tempo de permanência no interior do galpão.

Um modelo de turbulência baseado no conceito de vórtice

O fenômeno da turbulência está presente na maioria dos escoamentos observados na indústria e na natureza. Muitas são as considerações a respeito das dificuldades relacionadas à caracterização dos escoamentos turbulentos. Uma das muitas questões trata do procedimento de análise do problema através da descrição estatística dos campos por grandezas “médias”, o que leva ao problema de fechamento e à modelagem do tensor de Reynolds, normalmente com modelos baseados no conceito de viscosidade turbulenta. Os modelos de turbulência já existentes apresentam algumas deficiências na previsão do escoamento, além de outras limitações, o que justifica a busca por novas abordagens para o tratamento da turbulência. Neste trabalho, o problema de fechamento é tratado segundo a modelagem turbulenta baseada no conceito de viscosidade turbulenta. Um novo modelo de turbulência é proposto, que admite a existência de vórtices imersos no escoamento e aplica conceitos e definições relacionados à identificação de vórtices, com o uso do critério de identificação Q , que caracteriza a região do escoamento ocupada pelo vórtice. Propõe-se a investigação da aplicabilidade do critério Q em conjunto com o modelo k − ε , para o desenvolvimento de um novo modelo de turbulência chamado k − ε −Q . Validou-se a aplicabilidade do modelo através de um código numérico computacional para tratamento de escoamentos turbulentos. A solução numérica foi obtida através da discretização do domínio fluido, utilizando o método de volumes finitos e o método multigrid foi utilizado para resolver o sistema linear resultante. Como verificação, foi utilizado este modelo de turbulência para simular o escoamento em uma cavidade quadrada com tampa deslizante e o escoamento turbulento sobre um degrau. Os resultados obtidos foram confrontados com dados experimentais e demonstraram que o modelo aqui proposto se apresenta mais eficiente que o clássico modelo k − ε , no tratamento da turbulência nesses dois problemas clássicos.

Análise da resolução do modelo de Hughes em sondagens magnetotelúricas

Na maioria dos métodos de exploração geofísica, a interpretação é feita assumindo-se um modelo da Terra uniformemente estratificado. Todos os métodos de inversão, inclusive o de dados eletromagnéticos, exigem técnica de modelamento teórico de modo a auxiliar a interpretação. Na literatura os dados são geralmente interpretados em termos de uma estrutura condutiva unidimensional; comumente a Terra é assumida ser horizontalmente uniforme de modo que a condutividade é função somente da profundidade. Neste trabalho uma técnica semi-analítica de modelagem desenvolvida por Hughes (1973) foi usada para modelar a resposta magnética de duas camadas na qual a interface separando as camadas pode ser representada por uma expansão em série de Fourier. A técnica envolve um método de perturbação para encontrar o efeito de um contorno senoidal com pequenas ondulações. Como a perturbação é de primeira ordem a solução obtida é linear, podemos então usar o princípio da superposição e combinar soluções para várias senoides de forma a obter a solução para qualquer dupla camada expandida em série de Fourier. Da comparação com a técnica de elementos finitos, as seguintes conclusões podem ser tiradas: • Para um modelo de dupla camada da Terra, as camadas separadas por uma interface cuja profundidade varia senoidalmente em uma direção, as respostas eletromagnética são muito mais fortes quando a espessura da primeira camada é da ordem do skin depth da onda eletromagnética no meio, e será tanto maior quanto maior for o contraste de condutividade entre as camadas; • Por outro lado, a resistividade aparente para este modelo não é afetada pela mudança na frequência espacial (v) do contorno; • Em caso do uso da solução geral para qualquer dupla camada na Terra cuja interface possa ser desenvolvida em série de Fourier, esta técnica produziu bons resultados quando comparado com a técnica de elementos finitos. A linerização restringe a aplicação da técnica para pequenas estruturas, apesar disso, uma grande quantidade de estruturas pode ser modelada de modo simples e com tempo computacional bastante rápido; • Quando a dimensão da primeira camada possui a mesma ordem de grandeza da estrutura, esta técnica não é recomendada, porque para algumas posições de sondagem, as curvas de resistividade aparente obtidas mostram um pequeno deslocamento quando comparadas com as curvas obtidas por elementos finitos.

Simulação computacional da adsorção dos poluentes benzeno, tolueno e p-xileno sobre carvão ativado

Os maiores problemas de contaminação de aquíferos e solos são atribuídos aos hidrocarbonetos monoaromáticos, que são os constituintes mais solúveis e mais móveis da fração de algumas substâncias, como por exemplo, da gasolina. Para a remoção destes contaminantes, a adsorção por carvão ativado é o método mais utilizado, pois o carvão apresenta uma habilidade significativa para adsorver componentes orgânicos de baixo peso molecular, como o benzeno, o tolueno e o p-xileno. Neste trabalho, verificou-se a adsorção dos mesmos sobre carvão ativado via simulação computacional. Como base, utilizou-se o modelo postulado de carvão preparado por Bourke et al. (2007). Várias etapas foram concluídas desde os modelos das estruturas do carvão e dos poluentes até as simulações de dinâmica molecular. Para a análise conformacional da estrutura do carvão, foi utilizado o método semi-empírico PM3 e para o processo de dinâmica, o campo de força AMBER FF99SB. A estrutura passou por um aquecimento, à pressão constante, até alcançar uma temperatura final de 298K (25ºC), sendo suas informações coletadas a cada 50ps. Posteriormente, a estrutura foi submetida a equilíbrio de sistema, à temperatura constante de 298K (25ºC), por 500ps para então suas informações serem analisadas. Por fim, o sistema foi então submetido à dinâmica molecular durante 15 ns. Após análise dos resultados, constatou-se que os grupos éter, lactona e carbonila (cetona) presentes na estrutura de carvão ativado conferem caráter ácido à mesma e devido a isto e à sua consequente carga superficial negativa, a adsorção tornou-se viável uma vez que os poluentes apresentavam carga superficial positiva, o que corrobora o entendimento que já se tem a respeito desse tipo de fenômeno.

Modelo dinâmico de um túnel de vento de sopro

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Análise do escoamento turbulento por um queimador industrial a gás utilizando dinâmica dos fluidos computacional

Pós-graduação em Engenharia Mecânica - FEG

Momentos fracionários e localização dinâmica para o modelo de Anderson discreto multidimensional

Pós-graduação em Matematica Aplicada e Computacional - FCT

Um modelo para a dengue com influência sazonal

In this study, a mathematical model with temporal dependence for dengue transmission was developed, considering coupling between human population and the vector mosquito, and a sorotype circulating on population. This model was analysed with the goal to explain disease's periodicity. Finally, a genetic algorithm was set up to study model's sensibility.

Análises, aplicações e validações–numérico/experimentais do modelo SWAN em áreas restritas e ao largo

Pós-graduação em Engenharia Elétrica - FEIS

Técnica para dedução do modelo algébrico PTT, aplicações e análises numéricas em escoamentos bidimensionais

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Modelo dinâmico bayesiano multivariado para análise espaço-temporal de dados de área

Pós-graduação em Matematica Aplicada e Computacional - FCT

Instrumentação eletrônica para experimentação do modelo Bouncer Ball dissipativo

In our experiments, we studied the classical dynamic of a particle in vertical motion subject to a constant gravitational field and the partial shock with an elastic wall with oscillatory motion located below the particle. The motion of the particle was confined to the vertical direction by a glass guide tube, which was initially evacuated, allowing viscous drag forces to be disregarded. The time between impacts, as well as the time between the reference phase of the movement of the base and the moment of impact and the period of oscillation of the base will be acquired by a hardware and software to obtain the phases space to be compared later with the description of the system through mapping discrete variables via the 'particle velocity immediately after shock 'and' phase of the movement of the base at the instant of shock 'obtained by computer simulation. This requires developing an electronic analog system followed by Digital implemented in reconfigurable logic, more specifically a sequential machine able to discriminate the impact with the metal base from the impacts of the glass guide tube, based on the frequency spectrum of the response of the microphone to these different impacts

Sistema de cálculo computacional da unidade monitora para feixes de elétrons em radioterapia

In radiation theraphy with electron beam, the electrons are produced in linear accelerators, and energy the most used have between 4MeV and 20MeV. Generally, the treatments are done for superficial injuries, because the low penetration of these particles. In this work a system for calculation of monitor units (U.M.) for cases of treatments with electron beam was developed. The Excel program of Microsoft was used and is easily found in the operational system of the personal microcomputers. In the Excel has been inserted the pertinent data of the linear accelerator of Varian, model 2100C, used in the Service of radiation theraphy of the Hospital of the Clinics of the College of Medicine of the UNESP of Botucatu. For some values of the physical parameters, such as: factors field and factors calibration, not supplied in the tests of acceptance of the machine, still proceeded calculations from interpolation and extrapolation. The mathematical formulas for automatic search of these and others factors used in the calculations of the determination of the U.M had been developed in agreement available routines in Excel. For this the functions had been used the function IF (that it imposes search condition) and the PROCH (that looks a value in a column from determined line), beyond the basic functions of addition, multiplication and division. It is intended to optimize the routine of the Services of radiation theraphy that perform through eletrontheraphy procedures, speeding the calculations and minimizing the occurrence of errors and uncertainties deriving of the maken a mistake manipulation of the parameters gotten in tables of data of electron beams

Modelagem computacional das túnicas cardíacas para processamento e controle de qualidade

This paper considers a study of the anatomical features of the cardiac system and a three-dimensional model of the different tunics that comprise the heart wall, for processing and quality control of radiological images. The structures are built by the layer overlapping method, where a layer can be understood as a slice of the three-dimensional object. The pericardium, myocardium and endocardium were represented with three-dimensional cylinders and hexagons. The spatial arrangement of the cardiac system is determined by an background image of a real model, which values are defined according to the shape of the region and on the anatomical patients characteristics. The results are significant, considering the anatomical structures details, as well as the representation of the thicknesses of the regions of the heart wall. The validation of the anatomical model was accomplished through comparisons with dimensions obtained from a real model and allows verifying that the model is appropriate. The degree of representation will allow the verification of the influence of radiological parameters, morphometric peculiarities and stage of the diseases on the quality of the images, as well as on the performance of the Computer-Aided Diagnosis (CAD).

Considerações e avaliações das solicitações em estruturas de edificações de paredes de concreto pré-moldado utilizando o modelo de pórtico tridimensional: estudo de caso

The present study is in reference to precast concrete (panels) walls as the main structural system. The diversification of the structural systems is a reality due to specific characteristics', as a result there is a lack of study and analysis. Some systems are already established in their applications, for instance: conventional reinforced concrete, structural masonry, metal structures, and wood structure. Apart from this precast concrete elements has had a growth in building sites in Brazil, therefore professionals should be more aware. Affordable house funding initiatives to address the housing shortage has been recurrent throughout the country, additionally precast concrete walls has been an alternative that meets the three basic concepts of engineering, which are: the technique, economy, and security. The objective of this study is to gather concepts from the literature and TOMO (2012) about precast concrete walls and make a didactic synthesis of how they are viewed. The modeling structure method of the system was performed following the concept of assigning bars, using the computer program of calculations SAP2000. Initially will be presented a theoretical part, furthermore a case study to illustrate the sizing of the walls using EXCEL tables programmed with calculation routine. The results will be presented in calculated efforts by the program for further analysis. Finally, will be highlighted important items of modeling and interpretation of the results