Modelagem matemática e simulação computacional de derrames de óleo em corpos d'água: um estudo de caso na região do Porto de Santos

The objective of this work is to accomplish studies of mathematical modeling and computational simulation of oil spills in water bodies. For this reason, a case study in the region of the Port of Santos was developed using the softwares SisBAHIA and ADIOS2 for the simulation of different hypothetical scenarios of oil spilling on the surface of water, aiming to obtain information that contribute to the reduction of the possible environmental impacts that can be caused by such accidents. The results generated in the different simulations had shown that the obtained data can be extremely useful to subsidize the elaboration of mitigation plans, the mapping of risk areas or even the proposal of emergencial strategies in cases of real accidents, configuring the modeling and the simulation as important and modern tools for the environmental planning and management.

Concepção, modelagem, simulação e usinagem de um aerofólio automotivo auxiliado por computador

With the advancement of computer technology and the availability of technology computer aided design (CAD) errors in the designs are getting smaller. To this end the project aims to assess the reliability of the machine (CNC), which was designed by students of mechanical engineering college engineering - UNESP Bauru, by designing, modeling, simulation and machining an airfoil automotive. The profile template selected for the study will be a NACA 0012 machined plates in medium density fiberboard (MDF) and will be performed with a structural analysis simulation using finite elements and a software CFD (Computational Fluid Dynamics), and test the real scale model in a wind tunnel. The results obtained in the wind tunnel and CFD software will be compared to see the error in the machining process.

Simulação de modelagem qualitativa para avaliação preliminar da qualidade da água na Bacia do Ribeirão das Perdizes em Campos do Jordão/SP, como subsídio ao enquadramento

The purpose of this research was to present the principles for the implementation of a future water body classification program at the UGRHI-1 (Management Unit of Water Resources-1/São Paulo, Brazil) using qualitative and quantitative water modeling. Our study area was in Campos do Jordão, a city in the state of São Paulo, specifically in the Perdizes river watershed, which is included in the UGRHI-1. The water quality assessment used a modeling mehtodology that can support a future study for the classification of water bodies in the region. The behavior of DO and BOD parameters was assessed in the water quality scenarios, considering the variations of flows, loads, reductions in withdrawals and different efficiencies in sewage treatment. The model used was QUAL2E developed by USEPA (United States Environmental Protection Agency) that received a new graphic interface named QUAL2R model. The reference flows Q7,10, Q95% and Qm used were obtained through the DAEE (Department of Water and Electric Energy) hydrologic regionalization method. It was noted that of the 9 proposed scenarios, only scenario 6 that was predicted in the flow regime Qm with the existence of a treatment system capable of removing 93% of the BOD, in the Perdizes River remained 67.5% of the time in class 2. Scenario 8 predicted a 50% reduction in In scenario 8 a reduction of 50% in the flow captured in CA-04, associated with the efficiency of 93% of BOD removal keeping the river in class 2. In scenario 9 the minimum required efficiency in sewage treatment in the flow regime Q95% for keeping the river in class 2 was also calculated and the value of 94.7% was obtained.

Simulação e modelagem computacional de dados de espalhamento à baixos ângulos - enfoque em estruturas de alta simetria

Esta tese apresenta uma abordagem para a criação rápida de modelos em diferentes geometrias (complexas ou de alta simetria) com objetivo de calcular a correspondente intensidade espalhada, podendo esta ser utilizada na descrição de experimentos de es- palhamento à baixos ângulos. A modelagem pode ser realizada com mais de 100 geome- trias catalogadas em um Banco de Dados, além da possibilidade de construir estruturas a partir de posições aleatórias distribuídas na superfície de uma esfera. Em todos os casos os modelos são gerados por meio do método de elementos finitos compondo uma única geometria, ou ainda, compondo diferentes geometrias, combinadas entre si a partir de um número baixo de parâmetros. Para realizar essa tarefa foi desenvolvido um programa em Fortran, chamado de Polygen, que permite modelar geometrias convexas em diferentes formas, como sólidos, cascas, ou ainda com esferas ou estruturas do tipo DNA nas arestas, além de usar esses modelos para simular a curva de intensidade espalhada para sistemas orientados e aleatoriamente orientados. A curva de intensidade de espalhamento é calculada por meio da equação de Debye e os parâmetros que compõe cada um dos modelos, podem ser otimizados pelo ajuste contra dados experimentais, por meio de métodos de minimização baseados em simulated annealing, Levenberg-Marquardt e algorítmicos genéticos. A minimização permite ajustar os parâmetros do modelo (ou composição de modelos) como tamanho, densidade eletrônica, raio das subunidades, entre outros, contribuindo para fornecer uma nova ferramenta para modelagem e análise de dados de espalhamento. Em outra etapa desta tese, é apresentado o design de modelos atomísticos e a sua respectiva simulação por Dinâmica Molecular. A geometria de dois sistemas auto-organizado de DNA na forma de octaedro truncado, um com linkers de 7 Adeninas e outro com linkers de ATATATA, foram escolhidas para realizar a modelagem atomística e a simulação por Dinâmica Molecular. Para este sistema são apresentados os resultados de Root Mean Square Deviations (RMSD), Root Mean Square Fluctuations (RMSF), raio de giro, torção das hélices duplas de DNA além da avaliação das ligações de Hidrogênio, todos obtidos por meio da análise de uma trajetória de 50 ns.

Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico

LINS, Filipe C. A. et al. Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico. In: CONGRESSO BRASILEIRO DE PESQUISA E DESENVOLVIMENTO EM PETRÓLEO E GÁS, 5. 2009, Fortaleza, CE. Anais... Fortaleza: CBPDPetro, 2009.

Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico

LINS, Filipe C. A. et al. Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico. In: CONGRESSO BRASILEIRO DE PESQUISA E DESENVOLVIMENTO EM PETRÓLEO E GÁS, 5. 2009, Fortaleza, CE. Anais... Fortaleza: CBPDPetro, 2009.

Variação temporal de uma comunidade fitoplanctônica do reservatório de APM-Manso através de modelagem ecológica tridimensional.

A ecologia de reservatórios, que são ecossistemas complexos, dinâmicos e artificiais, vem assumindo destaque no Brasil. O objetivo deste trabalho foi avaliar a viabilidade da aplicação, no reservatório de APM-Manso, de um modelo ecológico tridimensional em estudos sobre a dinâmica fitoplanctônica, simulando a variação temporal do fitoplâncton para cenários distintos de carga de nutrientes. O modelo CAEDYM foi acoplado ao ELCOM e simulação foi realizada em duas etapas: uma hidrodinâmica e outra ecológica. Escolheu-se para as simulações o período de cinco meses, a partir de 1 de setembro de 2005. Foram construídos dois cenários de simulação, o primeiro contendo os valores reais de carga de nutrientes dos principais rios contribuintes medidos em campo, e o segundo com redução na carga nutricional destes rios, simulando um possível processo de substituição de áreas florestadas por áreas de pastagem na bacia do rio Manso. A comunidade fitoplanctônica simulada apresentou rápidas respostas à disponibilidade nutricional do ambiente, e os resultados obtidos corroboraram com diversas teorias sobre as estratégias adaptativas e sobre as dinâmicas algais. Dentre as classes simuladas, Bacillariophyceae e Cryptophyceae se mostraram mais sensíveis às reduções de carga, enquanto Chrolophyceae e Cyanophyceae, apesar de terem suas biomassas reduzidas, sofreram menos com o impacto, sugerindo estarem mais adaptadas à limitação de nutrientes. Os picos chuvosos influenciaram positivamente as taxas de crescimento das Bacillariophyceae apenas no Cenário 1, uma vez que a limitação por nutrientes foi mais decisiva para esta classe no Cenário 2. Observou-se em ambas as simulações uma tendência de substituição na dominância de Cyanophyceae por Chlorophyceae.

Modelagem hidrológica da bacia hidrográfica do rio Bengalas, Nova Friburgo, RJ, utilizando o potencial de geotecnologias na definição de áreas de risco de inundação.

A presente pesquisa concentra-se no estudo hidrológico utilizando o potencial das geotecnologias na modelagem do escoamento na bacia do rio Bengalas, cujo rio principal de mesmo nome, corta o município de Nova Friburgo, RJ, no sentido Sul-Norte. Esse município, um dos mais importantes da região serrana, sofre frequentemente com inundações e deslizamentos, onde, dados históricos e acontecimentos recentes mostram que a ocupação inadequada de encostas e calhas dos rios são as áreas destacada e negativamente afetadas. A metodologia tem suporte no uso de um SIG, extraindo informações, que por sua vez, serão entrada de dados na fase de modelagem, e reforçando a apresentação dos resultados das simulações através de mapas. Ela está divida basicamente em três etapas: "SIG", "Modelagem" e "Suporte à Tomada de Decisão/Simulação". Esse primeiro estudo permitiu compor um banco de dados geográfico com as características fisiográficas da bacia; a seleção criteriosa de uma modelagem matemática e encadeamento de seus parâmetros com os componentes do ciclo hidrológico; realizar a calibração do modelo de transformação chuva-vazão, Soil Conservation Service (CN); e simular a passagem dos volumes gerados pela precipitação efetiva na calha do rio Bengalas, com o objetivo de identificar e analisar as áreas suscetíveis a inundações na porção central da cidade de Nova Friburgo. Modelagem dessa natureza vem sendo empregada, principalmente, no gerenciamento de recursos hídricos, onde a tomada de decisões embasada nos resultados de simulações computacionais, contribuem para evitar prejuízos materiais e financeiros, e ainda, perdas de vidas humanas em áreas de risco, neste caso, aquelas suscetíveis a inundações. Analisando os resultados encontrados temos que a área suscetível à inundação para uma chuva com tempo de recorrência de 50 anos, o mais crítico estudado, seja de aproximadamente 1,0 km, distribuídos nos seus 8,5km na região central de Nova Friburgo-RJ, sendo está, ora delimitada, prioritariamente edificada.

Sim-Colmeia: ambiente de simulação da dinâmica de uma colmeia para o ensino de Biologia

As ferramentas computacionais estão apoiando, de maneira crescente, o processo de ensino e aprendizagem em diversas áreas. Elas aumentam as possibilidades do docente para ministrar um conteúdo e interagir com seus alunos. Neste grupo de ferramentas estão as simulações baseadas em sistemas multiagentes. Neste contexto, este trabalho tem como objetivo apresentar um ambiente de simulação do crescimento populacional de uma colmeia para o ensino de Biologia. As variáveis do sistema podem ser alteradas visando analisar os diferentes resultados obtidos. Aspectos como duração e tempo da florada das plantações, conhecidos como campos de flores, podem ser manipulados pelo aluno. A abordagem multiagentes em Inteligência Artificial Distribuída foi a solução escolhida, para que o controle das atividades do aplicativo fosse feito de maneira automatizada. A Realidade Virtual foi utilizada para acrescentar aspectos importantes do processo que não podem ser visualizados pela simulação matemática. Uma síntese da utilização de tecnologias na educação, em especial da Informática, é discutida no trabalho. Aspectos da aplicação no ensino de Biologia são apresentados, assim como resultados iniciais de sua utilização.

Simulação numérica do processo de alteamento de áreas de deposição de resíduos pelo método a montante.

A previsão do comportamento de resíduos constituiu-se em um desafio geotécnico, uma vez que estes materiais apresentam uma resposta distinta dos materiais usualmente encontrados em depósitos naturais. A análise dos recalques da fundação, decorrentes da sobrecarga imposta pelo alteamento, é complexa, tendo em vista que o adensamento de resíduos pressupõe grandes deformações, invalidando o uso de teorias clássicas de adensamento. Atualmente, no Brasil, a técnica de disposição de resíduos de bauxita prevê uma operação inicial de lançamento no interior de lagos artificiais, em forma de polpa. Após o esgotamento do lago e ressecamento do resíduo, inicia-se o lançamento pelo método a montante. Neste método, a polpa é lançada sobre o resíduo pré-existente, que se encontra em processo de adensamento. O presente trabalho tem como objetivo reproduzir numericamente o comportamento de áreas de resíduos durante a etapa de alteamento a montante. A pesquisa tem como enfoque 2 áreas de resíduos de bauxita. Uma delas encontra-se em fase de reabilitação e dispõe de instrumentação de campo (recalques e deslocamentos horizontais). A outra se encontra em fase de operação do alteamento e dispõe de dados experimentais. Desta forma, a metodologia consistiu na reprodução numérica do processo de alteamento da área instrumentada e comparação dos resultados com a instrumentação de campo, com objetivo de avaliar o modelo numérico e os parâmetros do resíduo. Posteriormente, realizou-se a previsão do comportamento do resíduo de fundação da área em fase de alteamento. Os parâmetros geotécnicos foram definidos a partir de um extenso programa de ensaios de campo e laboratório, executado no local em estudo, fazendo-se uso de um tratamento estatístico dos dados experimentais. Os resultados numéricos mostraram a potencialidade do programa na previsão do comportamento de áreas de resíduos durante o alteamento a montante, com previsões de recalques e deslocamentos horizontais coerentes com a instrumentação de campo.

Modelagem matemática da transferência de calor numa placa plana sob o efeito de uma fonte pontual externa de radiação térmica.

Este trabalho apresenta uma modelagem matemática para o processo de aquecimento de um corpo exposto a uma fonte pontual de radiação térmica. O resultado original que permite a solução exata de uma equação diferencial parcial não linear a partir de uma seqüência de problemas lineares também é apresentado. Gráficos gerados com resultados obtidos pelo método de diferenças finitas ilustram a solução do problema proposto.

Desenvolvimento e avaliação de novas abordagens de modelagem de processos de separação em leito móvel simulado

O Leito Móvel Simulado (LMS) é um processo de separação de compostos por adsorção muito eficiente, por trabalhar em um regime contínuo e também possuir fluxo contracorrente da fase sólida. Dentre as diversas aplicações, este processo tem se destacado na resolução de petroquímicos e principalmente na atualidade na separação de misturas racêmicas que são separações de um grau elevado de dificuldade. Neste trabalho foram propostas duas novas abordagens na modelagem do LMS, a abordagem Stepwise e a abordagem Front Velocity. Na modelagem Stepwise as colunas cromatográficas do LMS foram modeladas com uma abordagem discreta, onde cada uma delas teve seu domínio dividido em N células de mistura interligadas em série, e as concentrações dos compostos nas fases líquida e sólida foram simuladas usando duas cinéticas de transferência de massa distintas. Essa abordagem pressupõe que as interações decorrentes da transferência de massa entre as moléculas do composto nas suas fases líquida e sólida ocorram somente na superfície, de forma que com essa suposição pode-se admitir que o volume ocupado por cada molécula nas fases sólida e líquida é o mesmo, o que implica que o fator de residência pode ser considerado igual a constante de equilíbrio. Para descrever a transferência de massa que ocorre no processo cromatográfico a abordagem Front Velocity estabelece que a convecção é a fase dominante no transporte de soluto ao longo da coluna cromatográfica. O Front Velocity é um modelo discreto (etapas) em que a vazão determina o avanço da fase líquida ao longo da coluna. As etapas são: avanço da fase líquida e posterior transporte de massa entre as fases líquida e sólida, este último no mesmo intervalo de tempo. Desta forma, o fluxo volumétrico experimental é utilizado para a discretização dos volumes de controle que se deslocam ao longo da coluna porosa com a mesma velocidade da fase líquida. A transferência de massa foi representada por dois mecanismos cinéticos distintos, sem (tipo linear) e com capacidade máxima de adsorção (tipo Langmuir). Ambas as abordagens propostas foram estudadas e avaliadas mediante a comparação com dados experimentais de separação em LMS do anestésico cetamina e, posteriormente, com o fármaco Verapamil. Também foram comparados com as simulações do modelo de equilíbrio dispersivo para o caso da Cetamina, usado por Santos (2004), e para o caso do Verapamil (Perna 2013). Na etapa de caracterização da coluna cromatográfica as novas abordagens foram associadas à ferramenta inversa R2W de forma a determinar os parâmetros globais de transferência de massa apenas usando os tempos experimentais de residência de cada enantiômero na coluna de cromatografia líquida de alta eficiência (CLAE). Na segunda etapa os modelos cinéticos desenvolvidos nas abordagens foram aplicados nas colunas do LMS com os valores determinados na caracterização da coluna cromatográfica, para a simulação do processo de separação contínua. Os resultados das simulações mostram boa concordância entre as duas abordagens propostas e os experimentos de pulso para a caracterização da coluna na separação enantiomérica da cetamina ao longo do tempo. As simulações da separação em LMS, tanto do Verapamil quando da Cetamina apresentam uma discrepância com os dados experimentais nos primeiros ciclos, entretanto após esses ciclos iniciais a correlação entre os dados experimentais e as simulações. Para o caso da separação da cetamina (Santos, 2004), a qual a concentração da alimentação era relativamente baixa, os modelos foram capazes de predizer o processo de separação com as cinéticas Linear e Langmuir. No caso da separação do Verapamil (Perna, 2013), onde a concentração da alimentação é relativamente alta, somente a cinética de Langmuir representou o processo, devido a cinética Linear não representar a saturação das colunas cromatográficas. De acordo como o estudo conduzido ambas as abordagens propostas mostraram-se ferramentas com potencial na predição do comportamento cromatográfico de uma amostra em um experimento de pulso, assim como na simulação da separação de um composto no LMS, apesar das pequenas discrepâncias apresentadas nos primeiros ciclos de trabalho do LMS. Além disso, podem ser facilmente implementadas e aplicadas na análise do processo, pois requer um baixo número de parâmetros e são constituídas de equações diferenciais ordinárias.

Modelagem computacional da estabilização de sistemas subcríticos segundo o modelo unidimensional de difusão de nêutrons monoenergéticos

Um método numérico espectronodal (END) livre de erros de truncamento espacial é desenvolvido para problemas unidimensionais de difusão de nêutrons monoenergéticos em duas versões. Na versão de problemas de autovalor, o método gera soluções numéricas para o perfil do fluxo escalar e para o fator de multiplicação efetivo (k), que coincidem com a solução analítica dominante, afora os erros da aritmética finita computacional. Na versão de fonte fxa, o método também gera soluções numéricas analíticas para o problema de fonte fixa correspondente, onde a fonte de fissão, com dependência espacial, é obtida analiticamente, a partir da reconstrução espacial do fluxo escalar gerado pelo método END para problemas de autovalor. Alguns experimentos numéricos são apresentados para dois problemas modelos a fim de ilustrar a precisão do método.

Equilíbrio líquido-líquido entre ácido acético, água e fenólicos: experimental e modelagem

Com o passar do tempo, a população mundial vem se conscientizando mais sobre problemas ambientais. Isso fez surgir uma demanda por tecnologias novas que possam se encaixar no cenário de sustentabilidade. Instabilidades frequentes no cenário político-econômico global acabam por elevar o preço do barril do petróleo. Assim a indústria química começa a buscar por alternativas que tenham a mesma versatilidade do petróleo. Dentre as opções de combustíveis renováveis, destaca-se o bio-óleo de pirólise. Seu interesse em estudos científicos vem do fato de poder-se utilizar do rejeito de processos como matéria prima, não necessitando competir por espaço com a plantação de alimentos. Sua composição pode ser representada por ácidos e fenóis. Em especial destacamos o ácido acético e fenóis oxigenados como m-cresol, o-cresol, p-cresol e guaiacol por estarem presentes em grandes quantidades. Sua separação das frações menos polares do bio-óleo pode ser realizada por meio de extração com água que é um reagente abundante e de baixo custo. O conhecimento das propriedades desses componentes puros é amplamente disponível na literatura, porém dados de composições das fases corrosivas, como misturas ternárias de água-ácido acético-m-cresol/o-cresol/p-cresol/guaiacol nas temperaturas de (298,15, 310,65 e 323,15) K são escassos. Devido a isso, o uso de modelos termodinâmicos para simulação do comportamento desses sistemas torna-se interessante. Todavia, quando são testadas as capacidades dos modelos clássicos, como o UNIFAC e NRTL, percebe-se que os mesmos não conseguem reproduzir o comportamento da binodal dos componentes corrosivos. Sendo assim, essa dissertação investigou soluções para melhorar a descrição desses sistemas, assim como obteve dados experimentais para tais sistemas de misturas ternárias de água-ácido acético-m-cresol/o-cresol/p-cresol/guaiacol nas temperaturas de (298,15, 310,65 e 323,15) K; desenvolveu-se uma metodologia para estimar parâmetros do modelo NRTL a partir de dados de composição da binodal e verificou-se a possibilidade de utilizar o modelo UNIFAC para prever o comportamento de equilíbrio de fases. Como resultado foram obtidos 314 novos dados experimentais, obtiveram-se parâmetros para o modelo NRTL que conseguem reproduzir com acurácia a forma da binodal com a metodologia proposta e verificou-se a necessidade de aperfeiçoamento no estudo do modelo UNIFAC para otimizar sua utilização na previsão do comportamento dos sistemas estudados

Impacto do uso de biocombustíveis na qualidade do ar da cidade do Rio de Janeiro um estudo de simulação

A frota de veículos movidos por biocombustíveis tem aumentado nos últimos dez anos e com isso as emissões para a atmosfera nas cidades têm sofrido alterações. Devido a esta problemática, estudos com modelos de previsão de emissões são o foco deste trabalho. Alguns modelos, homologados por agências de regulamentação de alguns países ou comunidades, servem de base para as análises de risco com auxílio de simulação.O objetivo geral deste trabalho é avaliar o impacto do uso de biocombustíveis na qualidade do ar na cidade do Rio de Janeiro e estudar cenários de qualidade do ar, em função do aumento da quantidade de biodiesel e álcool adicionada ao diesel e gasolina. Nesta dissertação, monóxido de carbono, óxidos de nitrogênio, compostos orgânicos voláteis e ozônio são os principais poluentes estudados para a cidade do Rio de Janeiro. A avaliação de cenários foi realizada empregando um modelo de qualidade do ar com base no modelo de trajetórias OZIPR e no modelo químico SAPRC. Os resultados demonstram que o aumento do uso de biodiesel diminui a concentração de ozônio na atmosfera em relação ao caso base estudado, em 10,23% utilizando a mistura BE Diesel, em 5,28% utilizando a mistura B20 e apenas 0,33% utilizando a mistura B10. Isso de fato acontece pois o biodiesel possuiu mais oxigenados na sua estrutura. O estudo revelou o aumento nas emissões e concentrações de NOx na troposfera da cidade do Rio de Janeiro, para as mistura de BE- Diesel, B20 e o B10 em relação ao caso base, foram de 5,66%, 2,83% e 0% respectivamente. Em especial na cidade do Rio de janeiro existe um consumodo O3 pelo NO presente na atmosfera urbana. A redução na concentração de CO para as misturas de BE- Diesel, B20 e o B10 em relação ao caso base, foram de 13,11%, 6,56% e 4,12% respectivamente. O aumento da disponibilidade de oxigênio presente no biodiesel produz uma melhor queima do combustível, contribuindo para uma menor concentração de CO na troposfera