Modelando evolução por endossimbiose

Nesta dissertação é apresentada uma modelagem analítica para o processo evolucionário formulado pela Teoria da Evolução por Endossimbiose representado através de uma sucessão de estágios envolvendo diferentes interações ecológicas e metábolicas entre populações de bactérias considerando tanto a dinâmica populacional como os processos produtivos dessas populações. Para tal abordagem é feito uso do sistema de equações diferenciais conhecido como sistema de Volterra-Hamilton bem como de determinados conceitos geométricos envolvendo a Teoria KCC e a Geometria Projetiva. Os principais cálculos foram realizados pelo pacote de programação algébrica FINSLER, aplicado sobre o MAPLE.

Aplicaciones de las tecnologías de fabricación aditiva en la industria actual

[ES]En el ámbito industrial actual, la fabricación aditiva supone una alternativa para dar respuesta a las nuevas necesidades del mercado, que no pueden verse resueltas mediante las técnicas de fabricación clásicas. La fabricación aditiva (AM Adittive Manufacturing) representa un cambio en métodos y procesos industriales, permitiendo llevar a cabo piezas antes imposibles de realizar. Tiene sus principales aplicaciones en el sector médico, aeronáutico o el sector de la automoción. Se dan numerosas técnicas clasificadas por diferentes criterios, como el medio a través del cual se aporta la energía. Presenta numerosas ventajas como la complejidad en la geometría de las piezas que se pueden fabricar, a pesar de que aún no está totalmente desarrollada y presenta limitaciones como el coste de material. Todo el diseño parte de las tecnologías CAD/CAM 3D, y se fabrica ya en una amplia variedad de materiales, tanto metálicos como plásticos. Aún bajo poca normativa, se espera que en las próximas décadas esta tecnología desarrolle todo su potencial y sea de uso habitual en los sectores propicios para su explotación.

Estimativa do erro de discretização analítico na solução de equações diferenciais utilizando o Método de Volumes Finitos

As análises de erros são conduzidas antes de qualquer projeto a ser desenvolvido. A necessidade do conhecimento do comportamento do erro numérico em malhas estruturadas e não-estruturadas surge com o aumento do uso destas malhas nos métodos de discretização. Desta forma, o objetivo deste trabalho foi criar uma metodologia para analisar os erros de discretização gerados através do truncamento na Série de Taylor, aplicados às equações de Poisson e de Advecção-Difusão estacionárias uni e bidimensionais, utilizando-se o Método de Volumes Finitos em malhas do tipo Voronoi. A escolha dessas equações se dá devido a sua grande utilização em testes de novos modelos matemáticos e função de interpolação. Foram usados os esquemas Central Difference Scheme (CDS) e Upwind Difference Scheme(UDS) nos termos advectivos. Verificou-se a influência do tipo de condição de contorno e a posição do ponto gerador do volume na solução numérica. Os resultados analíticos foram confrontados com resultados experimentais para dois tipos de malhas de Voronoi, uma malha cartesiana e outra triangular comprovando a influência da forma do volume finito na solução numérica obtida. Foi percebido no estudo que a discretização usando o esquema CDS tem erros menores do que a discretização usando o esquema UDS conforme literatura. Também se percebe a diferença nos erros em volumes vizinhos nas malhas triangulares o que faz com que não se tenha uma uniformidade nos gráficos dos erros estudados. Percebeu-se que as malhas cartesianas com nó no centróide do volume tem menor erro de discretização do que malhas triangulares. Mas o uso deste tipo de malha depende da geometria do problema estudado

Evolução estrutural das zonas de cisalhamento dúcteis na porção centro-leste do domínio da Zona Transversal na Província Borborema

A área estudada está inserida no Domínio Transversal da Província Borborema. As unidades litoestratigráficas que compõem o embasamento paleoproterozócio (riaciano) são representadas por rochas ortoderivadas dos Complexos Salgadinho e Cabaceiras. Esses complexos foram individualizados de acordo com as suas diferenças composicionais, texturais e/ou geocronológicas. As rochas metassedimentares de idade paleoproterozóica (Orosiriano) foram interpretadas como constituintes do Complexo Sertânia. O magmatismo no estateriano é caracterizado por ortognaisses sienogranítcos da Suíte Carnoió-Caturité e por metanortositos do Complexo Metanortosítico Boqueirão. As unidades litoestratigráficas do Neoproterozóico são representadas por sucessões metassedimentares Criogenianas do Complexo Surubim e ortognaisses granodioríticos e sienograníticos do início do Ediacarano, denominados de Complexo Sumé e Ortognaisse Riacho de Santo Antônio, respectivamente. O magmatismo granítico do Ediacarano foi caracterizado pelo alojamento dos Plutons Inácio Pereira e Marinho. Os dados geocronológicos (U-Pb em zircão) obtidos indicam, no mínimo, o desenvolvimento de três eventos tectono-magmáticos. As idades de 2042 + 11Ma e 1996 + 13Ma obtidas nos ortoanfibolitos do Complexo Cabaceiras foram interpretadas como a idade de cristalização do protólito e metamorfismo, respectivamente. A idade de 1638 + 13Ma proveniente de hornblenda ortognaisse sienogranítico da Suíte Carnoió-Caturité foi interpretada como a idade de cristalização do protólito, marcando um evento magmático Estateriano de afinidade anorogênica. A idade de 550 + 3.1Ma encontrada em monzogranito porfirítico do Pluton Marinho é um registro do último evento magmático no final do Ediacarano, associado ao estágio tardio de desenvolvimento da Zona de Cisalhamento Coxixola. Os dados estruturais permitiram a individualização de três fases de deformação dúcteis, individualizadas como D1, D2 e D3. A fase D1 foi responsável pela geração de uma foliação S1, observada somente na charneira de dobras F2. O evento D2 é assinalado por uma tectônica contracional com transporte para NNW, observado a partir de bandas de cisalhamento assimétricas e dobras de arrasto em cortes paralelos a lineação de estiramento (L2x). Zonas de cisalhamento dúcteis de geometria e cinemática distintas desenvolveram-se durante a fase D3. As zonas de Cisalhamento Boa Vista, Carnoió e Congo estão orientadas na direção NE-SW e exibem cinemática sinistral em cortes paralelos à lineação de estiramento (L3x). As terminações meridionais dessas zonas de cisalhamento estão conectadas com a Zona de Cisalhamento Coxixola. Essa zona de cisalhamento, de direção WSW-ENE e cinemática destral, atravessa toda a área de estudo, com uma espessura média de rochas miloníticas de 300m. A Zona de Cisalhamento Inácio Pereira ocorre na porção leste da área de estudo, orientada na direção WNW-ESE. A análise geométrica e cinemática dessa zona de cisalhamento sugere uma evolução deformacional através de regime transpressivo oblíquo sinistral. O padrão anastomosado final resultante do desenvolvimento de todas as zonas de cisalhamento da área é relacionado à evolução estrutural de um sistema de zonas de cisalhamento dúcteis conjugadas.

Análise global da estabilidade termodinâmica de misturas: um estudo com o método do conjunto gerador

O cálculo do equilíbrio de fases é um problema de grande importância em processos da engenharia, como, por exemplo, na separação por destilação, em processos de extração e simulação da recuperação terciária de petróleo, entre outros. Mas para resolvê-lo é aconselhável que se estude a priori a estabilidade termodinâmica do sistema, a qual consiste em determinar se uma dada mistura se apresenta em uma ou mais fases. Tal problema pode ser abordado como um problema de otimização, conhecido como a minimização da função distância do plano tangente à energia livre de Gibbs molar, onde modelos termodinâmicos, de natureza não convexa e não linear, são utilizados para descrevê-lo. Esse fato tem motivado um grande interesse em técnicas de otimização robustas e eficientes para a resolução de problemas relacionados com a termodinâmica do equilíbrio de fases. Como tem sido ressaltado na literatura, para proporcionar uma completa predição do equilíbrio de fases, faz-se necessário não apenas a determinação do minimizador global da função objetivo do teste de estabilidade, mas também a obtenção de todos os seus pontos estacionários. Assim, o desenvolvimento de metodologias para essa tarefa desafiadora tem se tornado uma nova área de pesquisa da otimização global aplicada à termodinâmica do equilíbrio, com interesses comuns na engenharia química e na engenharia do petróleo. O foco do presente trabalho é uma nova metodologia para resolver o problema do teste de estabilidade. Para isso, usa-se o chamado método do conjunto gerador para realizar buscas do tipo local em uma rede de pontos previamente gerada por buscas globais efetuadas com uma metaheurística populacional, no caso o método do enxame de partículas.Para se obter mais de um ponto estacionário, minimizam-se funções de mérito polarizadas, cujos pólos são os pontos previamente encontrados. A metodologia proposta foi testada na análise de quatorze misturas polares previamente consideradas na literatura. Os resultados mostraram que o método proposto é robusto e eficiente a ponto de encontrar, além do minimizador global, todos os pontos estacionários apontados previamente na literatura, sendo também capaz de detectar, em duas misturas ternárias estudadas, pontos estacionários não obtidos pelo chamado método de análise intervalar, uma técnica confiável e muito difundida na literatura. A análise do teste de estabilidade pela simples utilização do método do enxame de partículas associado à técnica de polarização mencionada acima, para a obtenção de mais de um ponto estacionário (sem a busca local feita pelo método do conjunto gerador em uma dada rede de pontos), constitui outra metodologia para a resolução do problema de interesse. Essa utilização é uma novidade secundária deste trabalho. Tal metodologia simplificada exibiu também uma grande robustez, sendo capaz de encontrar todos os pontos estacionários pesquisados. No entanto, quando comparada com a abordagem mais geral proposta aqui, observou-se que tal simplificação pode, em alguns casos onde a função de mérito apresenta uma geometria mais complexa, consumir um tempo de máquina relativamente grande, dessa forma é menos eficiente.

Um método de matriz resposta com esquema iterativo de inversão parcial por região para problemas unidimensionais de transporte de nêutrons monoenergéticos na formulação de ordenadas discretas

Um método de matriz resposta (RM) é descrito para gerar soluções numéricas livres de erros de truncamento espacial para problemas de transporte de nêutrons monoenergéticos e com fonte fixa, em geometria unidimensional na formulação de ordenadas discretas (SN). O método RM com esquema iterativo de inversão parcial por região (RBI) converge valores numéricos para os fluxos angulares nas fronteiras das regiões que coincidem com os valores da solução analítica das equações SN, afora os erros de arredondamento da aritmética finita computacional. Desenvolvemos um esquema numérico de reconstrução espacial, que fornece a saída para os fluxos escalares de nêutrons em qualquer ponto do domínio definido pelo usuário, com um passo de avanço também escolhido pelo usuário. Resultados numéricos são apresentados para ilustrar a precisão do presente método em cálculos de malha grossa.

Avaliação númérica de ensaios push-out com conectores de cisalhamento tipo pino e perfobond.

Recentemente a utilização de estruturas mistas, aço concreto, vem ganhando espaço nas construções. Isso deve-se principalmente a economia de tempo em sua execução. Além disto, nesta solução a utilização do aço e do concreto otimiza as características estruturais de seus elementos: a resistência a tração do aço e a compressão do concreto. A transferência dos esforços entre os dois materiais possui grande influência no desempenho de uma estrutura mista, sendo comum a utilização de conectores de cisalhamento na região da interface entre estes dois materiais. O Eurocode 4 define um ensaio experimental denominado push-out de modo a determinar a resistência e ductilidade de conectores de cisalhamento. Seu desempenho é influenciado pelas resistências do concreto a compressão, as dimensões e taxa de armadura da laje de concreto, dimensões do perfil de aço, a disposição e a geometria dos conectores e pelas características dos aços utilizados no conector, no perfil e nas barras de reforço. Nota-se com isso uma grande quantidade de variáveis que influenciam o ensaio. Assim, o presente trabalho apresenta o desenvolvimento de um modelo em elementos finitos com base no programa ANSYS para simulação de ensaios push-out. Os resultados numéricos apresentados neste trabalho foram calibrados com resultados obtidos em ensaios experimentais existentes na literatura de ensaios push-out para conectores do tipo pino com cabeça (stud) e conectores tipo perfobond. Estes últimos apresentam elevada resistência sendo influenciados por inúmeros fatores como: número e diâmetro dos furos no conector e a inclusão ou não de barras de reforço extras nestes furos.

Intensidade acústica útil: um novo método para identificação de regiões radiantes em superfícies com geometrias arbitrárias

Neste trabalho é descrita a teoria necessária para a obtenção da grandeza denominada intensidade supersônica, a qual tem por objetivo identificar as regiões de uma fonte de ruído que efetivamente contribuem para a potência sonora, filtrando, consequentemente, a parcela referente às ondas sonoras recirculantes e evanescentes. É apresentada a abordagem de Fourier para a obtenção da intensidade supersônica em fontes com geometrias separáveis e a formulação numérica existente para a obtenção de um equivalente à intensidade supersônica em fontes sonoras com geometrias arbitrárias. Este trabalho apresenta como principal contribuição original, uma técnica para o cálculo de um equivalente à intensidade supersônica, denominado aqui de intensidade acústica útil, capaz de identificar as regiões de uma superfície vibrante de geometria arbitrária que efetivamente contribuem para a potência sonora que será radiada. Ao contrário da formulação numérica existente, o modelo proposto é mais direto, totalmente formulado na superfície vibrante, onde a potência sonora é obtida através de um operador (uma matriz) que relaciona a potência sonora radiada com a distribuição de velocidade normal à superfície vibrante, obtida com o uso do método de elementos finitos. Tal operador, chamado aqui de operador de potência, é Hermitiano, fato crucial para a obtenção da intensidade acússtica útil, após a aplicação da decomposição em autovalores e autovetores no operador de potência, e do critério de truncamento proposto. Exemplos de aplicações da intensidade acústica útil em superfícies vibrantes com a geometria de uma placa, de um cilindro com tampas e de um silenciador automotivo são apresentados, e os resultados são comparados com os obtidos via intensidade supersônica (placa) e via técnica numérica existente (cilindro), evidenciando que a intensidade acústica útil traz, como benefício adicional, uma redução em relação ao tempo computacional quando comparada com a técnica numérica existente.

Simulações de problemas inversos com aplicações em engenharia nuclear usando técnicas de transporte de partículas neutras monoenergéticas na formulação unidimensional de ordenadas discretas

Neste trabalho, três técnicas para resolver numericamente problemas inversos de transporte de partículas neutras a uma velocidade para aplicações em engenharia nuclear são desenvolvidas. É fato conhecido que problemas diretos estacionários e monoenergéticos de transporte são caracterizados por estimar o fluxo de partículas como uma função-distribuição das variáveis independentes de espaço e de direção de movimento, quando os parâmetros materiais (seções de choque macroscópicas), a geometria, e o fluxo incidente nos contornos do domínio (condições de contorno), bem como a distribuição de fonte interior são conhecidos. Por outro lado, problemas inversos, neste trabalho, buscam estimativas para o fluxo incidente no contorno, ou a fonte interior, ou frações vazio em barras homogêneas. O modelo matemático usado tanto para os problemas diretos como para os problemas inversos é a equação de transporte independente do tempo, a uma velocidade, em geometria unidimensional e com o espalhamento linearmente anisotrópico na formulação de ordenadas discretas (SN). Nos problemas inversos de valor de contorno, dado o fluxo emergente em um extremo da barra, medido por um detector de nêutrons, por exemplo, buscamos uma estimativa precisa para o fluxo incidente no extremo oposto. Por outro lado, nos problemas inversos SN de fonte interior, buscamos uma estimativa precisa para a fonte armazenada no interior do domínio para fins de blindagem, sendo dado o fluxo emergente no contorno da barra. Além disso, nos problemas inversos SN de fração de vazio, dado o fluxo emergente em uma fronteira da barra devido ao fluxo incidente prescrito no extremo oposto, procuramos por uma estimativa precisa da fração de vazio no interior da barra, no contexto de ensaios não-destrutivos para aplicações na indústria. O código computacional desenvolvido neste trabalho apresenta o método espectronodal de malha grossa spectral Greens function (SGF) para os problemas diretos SN em geometria unidimensional para gerar soluções numéricas precisas para os três problemas inversos SN descritos acima. Para os problemas inversos SN de valor de contorno e de fonte interior, usamos a propriedade da proporcionalidade da fuga de partículas; ademais, para os problemas inversos SN de fração de vazio, oferecemos a técnica a qual nos referimos como o método físico da bissecção. Apresentamos resultados numéricos para ilustrar a precisão das três técnicas, conforme descrito nesta tese.

Condições de contorno albedo para cálculos globais de reatores nucleares térmicos com o modelo de ordenadas discretas a dois grupos de energia

Como eventos de fissão induzida por nêutrons não ocorrem nas regiões nãomultiplicativas de reatores nucleares, e.g., moderador, refletor, e meios estruturais, essas regiões não geram potência e a eficiência computacional dos cálculos globais de reatores nucleares pode portanto ser aumentada eliminando os cálculos numéricos explícitos no interior das regiões não-multiplicativas em torno do núcleo ativo. É discutida nesta dissertação a eficiência computacional de condições de contorno aproximadas tipo albedo na formulação de ordenadas discretas (SN) para problemas de autovalor a dois grupos de energia em geometria bidimensional cartesiana. Albedo, palavra de origem latina para alvura, foi originalmente definido como a fração da luz incidente que é refletida difusamente por uma superfície. Esta palavra latina permaneceu como o termo científico usual em astronomia e nesta dissertação este conceito é estendido para reflexão de nêutrons. Este albedo SN nãoconvencional substitui aproximadamente a região refletora em torno do núcleo ativo do reator, pois os termos de fuga transversal são desprezados no interior do refletor. Se o problema, em particular, não possui termos de fuga transversal, i.e., trata-se de um problema unidimensional, então as condições de contorno albedo, como propostas nesta dissertação, são exatas. Por eficiência computacional entende-se analisar a precisão dos resultados numéricos em comparação com o tempo de execução computacional de cada simulação de um dado problema-modelo. Resultados numéricos para dois problemas-modelo com de simetria são considerados para ilustrar esta análise de eficiência.

Simulação com o código GEANT4 de medida de espessura de revestimento metálico em metal por XRF

Nesta dissertação são apresentados resultados de simulações Monte Carlo de fluorescência de raios X (XRF), utilizando o programa GEANT4, para medidas de espessura de revestimento metálico (Ni e Zn) em base metálica (Fe). As simulações foram feitas para dois tamanhos de espessura para cada metal de revestimento, (5μm e 10μm), com passos de 0,1 μm e 0,001 μm e com 106 histórias. No cálculo da espessura do revestimento foram feitas as aproximações de feixe de raios X monoenegético, com a análise da transmissão apenas da energia do K-alfa e para uma geometria compatível com um sistema real de medição (ARTAX-200). Os resultados mostraram a eficiência da metodologia de simulação e do cálculo da espessura do revestimento, o que permitirá futuros cálculos, inclusive para multirevestimentos metálicos em base metálica.

Simulación de la refrigeración en la estampación en caliente para una pieza de chapa.

[ES]El objetivo principal de este proyecto se centra en conseguir las características mecánicas requeridas por los componentes del chasis de los vehículos con una reducción de peso y mejora de productividad, para ellos es necesario simular el proceso de templado que se da durante la estampación en caliente de una pieza de chapa mediante utillajes con conductos de refrigeración. Para ello se ha utilizado el programa de elementos finitos ANSYS hasta obtener un proceso de embutición simulado que represente con una desviación dentro de los límites aceptables el comportamiento real de la chapa en la etapa de refrigeración. Como objetivo secundario se encuentra el afianzamiento de las bases teóricas de ciencia de materiales y la adquisición de más conocimientos relacionados con la transmisión de calor entre cuerpos sólidos, centrándose sobre todo en la distribución de temperaturas sobre la superficie de éstos. En una primera parte se tratarán los conceptos generales de la estampación en caliente y sus posibles variantes. También se explicará la necesidad del uso de nuevos materiales para la industria automovilística, así como la razón por la cual se utilizan conductos de refrigeración. A continuación, se definirá la geometría de la chapa a analizar, tanto las consecuentes geometrías de los utillajes, que tendrán diferentes distribuciones de conductos de refrigeración. Además se establecerán los criterios para realizar el análisis térmico transitorio del conjunto troquel – chapa – matriz. Una vez seleccionado el tipo de análisis se profundizará en su estudio, aplicándolo a los diferentes utillajes ya citados. Se analizarán los resultados obtenidos y los errores y se buscarán posibles alternativas. Finalmente, se procederá a sacar las conclusiones de la simulación realizada y se procederá a comparar los resultados obtenidos con las diferentes distribuciones de conductos de refrigeración en los utillajes.

Diseño y fabricación de una boquilla coaxial de aporte de material para un láser.

[ES]Este trabajo de fin de grado tiene como objetivo principal el diseño y la posterior fabricación de una boquilla coaxial que aporte material de cobertura para “laser cladding”. Para ello se utilizarán programas de diseño gráfico para hacer la geometría y programas de elementos finitos para simular el comportamiento del polvo y el gas dentro de la boquilla.

Estampación con Pam-Stamp de la parte interna de un capó.

[ES]El objetivo principal de este documento es diseñar todo el proceso de fabricación de la parte interna de capó. Esta acción se realiza con la ayuda de un software informático llamado Pam-Stamp. Es un programa de simulación virtual de procesos de conformado de chapa. El proceso elegido es la estampación en frío. Como resultados final se conocen la geometría tras la estampación, así como todas las propiedades físicas del producto elegido.

Simulación por elementos finitos del conformado de una pieza para carrocería de un automóvil.

[ES]El presente proyecto tiene como objetivo optimizar los parámetros de conformado a partir de la geometría de los troqueles de una pieza de un automóvil, con el fin de asegurar los parámetros de calidad que se exigen a este tipo de procesos se cumplan. Para ello, se dispone de un programa de simulación en elementos finitos cuyo nombre es Pam-Stamp 2G, en el cual se va a llevar a cabo el diseño. Además, se dispone de la matriz de la pieza requerida para una mayor facilidad de resolución. El proyecto constará de una descripción detallada del proceso de diseño del conformado durante todo el curso académico, complementado con un diagrama Gantt, así como las decisiones adoptadas durante el trabajo para lograr un desarrollo óptimo del mismo. Por último, contendrá también un presupuesto en el que se detallará los gastos que ha supuesto el desarrollo completo del trabajo.