Aceleração em Cell/B.E. da animação de superfícies deform

A animação de superfícies deformáveis, nomeadamente a modelação de tecidos, atravessa hoje uma época de grande relevância na indústria do cinema e no mundo dos jogos. A grande dedicação a este tema, em termos de investigação e a evolução das capacidades das arquitecturas de computadores no que toca a poder de processamento, tornou hoje possível efectuar este tipo de simulações usando um vasto leque de técnicas com diferentes objectivos. Entre estas técnicas encontra-se a simulação através de modelos discretos. Geralmente, neste tipo de modelação, as características do tecido são discretizadas num sistema de partículas organizadas entre si segundo um esquema de forças ou energias internas. Assim, a simulação pode ser efectuada integrando o sistema de forma a calcular as novas posições das partículas ao longo do tempo. Este tipo de computação é normalmente caracterizado como sendo bastante intensivo. A aceleração da animação de superfícies deformáveis recorrendo ao poder de processamento para além do CPU convencional foi realizada em vários trabalhos. No entanto, apenas uma pequena parte desses artigos está relacionada com a arquitectura Cell/B.E. O Cell/B.E. foi desenvolvido por uma equipa de investigadores vindos da Toshiba, Sony e IBM. Esta equipa tinha como objectivo a criação de uma arquitectura que suportasse um elevado leque de aplicações, incluindo o suporte de uma consola de jogos, de forma eficaz e com baixo consumo de energia. Assim, o processador Cell/B.E. convencional pode ser descrito por um chip multicore heterogéneo composto por um processador PowerPC e oito processadores vectoriais (SIMD) de 128 bits, permitindo assim ao programador uma maior flexibilidade na forma de paralelização de um determinado processamento. O principal objectivo deste trabalho passou pelo estudo desta arquitectura e da forma de a explorar e avaliar as suas capacidades, aplicando-as na aceleração de um simulador de superfícies deformáveis com

Métodos analíticos e numéricos do tipo elementos finitos contínuos e contínuos/descontínuos para o estudo de modelos de Boussinesq melhorados para a propagação de ondas

Dissertação para obtenção do Grau de Doutor em Matemática na área de especialização de Análise Numérica

Desenvolvimento de formulações de elementos finitos para barras sujeitas a cargas transitórias

Com o presente trabalho pretende-se formular, implementar e validar duas classes de elementos finitos não-convencionais para problemas elastoestáticos e elastodinâmicos (harmónicos e transitórios) envolvendo barras solicitadas por cargas axiais. O desempenho numérico dos elementos não convencionais é estudado para uma larga gama de situações de interesse prático e comparado com o dos elementos finitos conformes de deslocamento (convencionais). A resolução de problemas transitórios envolve a integração no tempo e no espaço das equações diferenciais governativas, bem como a imposição das respetivas condições iniciais e de fronteira. A metodologia de integração no tempo adotada neste trabalho é baseada no método de Newmark. A resolução de problemas estáticos e harmónicos não carece de integração no tempo, ou a mesma é feita de forma trivial. Concluída a discretização no tempo, a segunda fase da resolução envolve a integração no espaço de cada uma das equações discretizadas, nomeadamente através do método dos elementos finitos. Para esse efeito, apresentam-se as formulações relativas aos elementos finitos convencionais, híbridos e híbridos-Trefftz. As três formulações têm como ponto de partida a forma fraca da equação diferencial de Navier, que é imposta utilizando o método de Galerkin. A principal diferença entre os elementos convencionais e não-convencionais prende-se com a maneira como são impostas as condições de fronteira de Dirichlet e as condições de compatibilidade nas fronteiras interiores. Os elementos não convencionais são implementados numa plataforma computacional desenvolvida de raiz no ambiente Matlab. A implementação é feita de maneira a permitir uma definição muito geral e flexível da estrutura e das respetivas ações, bem como das discretizações no tempo e no espaço e das bases de aproximação, que podem ser diferentes para cada elemento finito. Por fim, efetuam-se testes numéricos com o objetivo de analisar os resultados obtidos com os elementos não convencionais e de os comparar com as respetivas soluções analíticas (caso existam), ou com os resultados obtidos utilizando elementos convencionais. É especialmente focada a convergência das soluções aproximadas sob refinamentos da malha (h), no espaço e no tempo, e das funções de aproximação (p), sendo que o uso simultâneo dos dois refinamentos parece conduzir mais rapidamente a soluções próximas da solução exata. Analisam-se também problemas complexos, envolvendo propagação de ondas de choque, com o fim de se efetuar uma comparação entre os elementos convencionais, disponíveis no programa comercial SAP2000, e os elementos não convencionais fornecidos pela plataforma computacional desenvolvida neste trabalho.

TROCAS DE MASSA E DE ENERGIA ENTRE SUPERFÍCIES NATURAIS E A ATMOSFERA

A modelização matemática das trocas de massa e de energia resultantes da interação entre vegetação e atmosfera é tratada pela combinação de análises das equações diferenciais representando a conservação de propriedades do ar (temperatura e umidade) com a equação do balanço de energia local do meio vegetal. Fluxos e fontes de calor e de massa, perfis de temperatura e umidade foram calculados através das relações matriciais das equações diferenciais discretizadas. As transferências radiativas são definidas no modelo, considerando a vegetação como um meio túrbido, onde as superfícies vegetais são distribuídas de forma aleatória no espaço. Os resultados produzidos pelo modelo para floresta são comparados com medidas micrometeorológicas realizadas na Reserva Florestal Ducke, Manaus - AM.

Da modelação matemática à simulação computacional: a experimentação matemática no ensino

Tese de Doutoramento em Ciências (área de especialização em Matemática).

Da modelação matemática à simulação computacional: a experimentação matemática no ensino

Tese de Doutoramento em Ciências (área de especialização em Matemática).

Estabilidade e Fenómenos de Localização em Elementos de Painel Sanduíche

É estudado o comportamento de pós-encurvadura de painéis sanduiche axialmente comprimidos. Aborda-se o problema com formulação alternativas, permitindo investigar a possibilidade de configurações deformadas deferentes daquelas obtidas por uma análise de Rayleigh-Ritz convencional. São usados alguns princípios gerais de elasticidade, deduzindo-se a função de emergia potencial total, dependente de um campo de deslocamentos. Admite-se um campo harmónico de amplitude constante, que é introduzido na função da energia potencial total, e discutem-se algumas diferenças observadas nos estudos, em relação a modelos estudados anteriormente. As técnicas utilizadas na pesquisa de soluções localizadas são introduzidas com o estudo do problema da viga em função elástica não linear. Analogia dinâmica e a técnica de dupla escala com perturbação são utilizadas para a introdução de conceitos de localização. Por último, aplica-se a um painel sanduíche a técnica de dupla escala com perturbação para deduzir as equações diferenciais nas amplitudes. Discute-se a solução destas equações com alguns exemplos, fazendo-se a resolução analítica e numérica das equações diferencias nas amplitudes.

Modelação cinética da conversão do furfural utilizando catalisadores multifuncionais

A utilização da biomassa como fonte de matéria-prima renovável para a obtenção de biocombustíveis e produtos de consumo pode reduzir os riscos geopolíticos e de segurança energética. As hemiceluloses – um dos principais componentes da biomassa vegetal - são convertidas industrialmente em furfural (Fur) que tem uma vasta gama de aplicações, sendo, por isso, considerado um dos compostos químicos renováveis de base mais importantes para biorefinarias. Este trabalho está relacionado com a valorização química do Fur, inserindo-se no conceito da Biorefinaria e nos objetivos de investigação do laboratório associado CICECO. Foi investigada a conversão catalítica do Fur em determinados bio-produtos de interesse para diversos setores da indústria química. Em particular, alquilfurfuril éteres, ésteres alquil-4-oxopentanoato, 5-metil-2-furanonas, ácido levulínico e gama-valerolactona são importantes para a produção de bio-aditivos para combustíveis, bio-polímeros, agroquímicos, aromas e fragrâncias, solventes, etc. A conversão do Fur nos bio-produtos envolve mecanismos reacionais complexos e requer o uso de catalisadores ácidos e de redução. A realização das múltiplas etapas num único reator é importante em termos de intensificação de processos. Neste trabalho, investigou-se a conversão do Fur em bio-produtos, num reator descontínuo fechado, a 120 ºC, utilizando catalisadores multifuncionais nanocristalinos do tipo zeólito beta contendo centros ativos de zircónio (Zr) e alumínio (Al). Com base numa revisão bibliográfica, propuseram-se dois modelos cinéticos pseudo-homogéneos para a conversão do Fur. Os modelos foram testados por ajuste das equações diferenciais ordinárias, representativas dos balanços às espécies reativas num reator descontínuo isotérmico, aos resultados experimentais de concentração versus tempo do Fur e cada um dos bio-produtos. Os estudos cinéticos foram combinados com os de caracterização dos catalisadores no estudo da influência das propriedades dos materiais nos diferentes passos do mecanismo global. Os centros ácidos de Al e Zr promoveram as reações ácidas e os de Zr eram essenciais para os passos de redução. No geral, as constantes cinéticas aparentes relativas às reações ácidas desejadas aumentaram com a razão Al/Zr dos catalisadores, e as dos passos de redução aumentaram com a diminuição da razão Si/Zr. O aumento da razão Al/Zr simultaneamente promoveu as reações indesejadas, diminuindo o rendimento total dos bio-produtos.

Estabelecimento de rotação econômica para uma floresta regulada

O objetivo deste trabalho foi determinar as rotações econômicas para um povoamento florestal, considerando-se um único corte e infinitos cortes, e para uma floresta regulada, por meio de modelos matemáticos, bem como compará-las através da análise marginal das condições de otimalidade. Para validar o modelo e verificar a magnitude dos efeitos, foi utilizado um estudo de caso de um projeto florestal com valores reais de produção, custos, receitas e taxa de desconto. Os resultados indicaram que os modelos foram eficientes para determinar a rotação econômica e que a rotação do povoamento para um único corte e infinitos cortes e para a floresta regulada foi, respectivamente, de 6,5; 5,5; e 5 anos. O Valor Esperado do Solo (VES) para a floresta regulada foi inferior ao VES do povoamento, em virtude de se impor ao manejo a condição de regulação.

Biomatemática: métodos e limitações

Este trabalho trata de uma visão panorâmica dos principais métodos utilizados, atualmente, na modelagem matemática de processos biológicos. Esses métodos são muito variados, com enfoques que, às vezes, são determinísticos, utilizando a teoria dos sistemas dinâmicos e, às vezes, são probabilísticos, lançando mão das equações diferenciais estocásticas. Em todos eles, entretanto, três fatores inerentes aos processos biológicos não podem deixar de ser considerados: a complexidade, a imprevisibilidade e a diversidade de escalas temporais.

A importância da matemática e da modelação matemática no mestrado integrado em ciências farmacêuticas

A Matemática e as Ciências Farmacêuticas encontram-se relacionadas desde há muito, no entanto, foi a partir do séc. XVII, período de notável agitação cultural e científico que os métodos experimentais foram sustentados com cálculos matemáticos. Esta ciência e as técnicas de modelagem matemática tornaram-se numa ferramenta amplamente utilizada, de tal modo, que nos dias de hoje são consideradas como fundamentais na generalidade das profissões e em especial nas Ciências Farmacêuticas. Contudo, para muitos ainda não é vista como fundamental e essencial para a formação de futuros farmacêuticos. Deste modo, pretende-se demonstrar como a Matemática e as técnicas de modelagem se tornaram ao longo dos anos nesta poderosa ferramenta. Quer pelos instrumentos, quer pelas competências que nos proporcionam. Pretende-se também, com recurso aos conteúdos programáticos desta unidade curricular, avaliar se os conhecimentos, sistemas de avaliação e distribuição da carga horária são efetuados de forma homogénea pelas diferentes instituições portuguesas, públicas ou privadas que lecionam o Mestrado Integrado em Ciências Farmacêuticas. Verificou-se que a Matemática é uma ciência plena de capacidades e recursos e que estabelece uma relação interdisciplinar com as Ciências Farmacêuticas. Quer pela componente utilitária, quer pela componente formativa que proporciona. A análise dos conteúdos programáticos demonstra que apesar de serem transversais, as Universidades que não lecionam Sistemas de Equações Lineares e Equações diferenciais deveriam faze-lo e também realizarem um melhor controlo da carga horária por temática.

Simulação numérica de um sistema de armazenamento de calor em estufas de plasticultura

Apresentar um modelo para simular um sistema de armazenamento de calor no solo em estufas para plasticultura é o objetivo do presente trabalho. O sistema consiste num feixe de tubos enterrados no solo. A convecção forçada de ar no seu interior realiza a troca térmica necessária para manter as estufas sob faixas desejadas de temperatura. O objetivo do modelo é investigar os efeitos no calor armazenado e a influência das variáveis, tais como diâmetro, comprimento, espaçamento entre os tubos e a velocidade de ar no canal provocam no sistema. O solo é tratado como um meio difusivo e avalia-se a contribuição do termo de condensação e evaporação da água contida no ar em escoamento nos tubos. A equação da energia é resolvida para o solo e para o ar. Os tubos de seção transversal circular são modelados como tubos de seção transversal quadrada com o objetivo de que as simulações possam ser processadas em coordenadas cartesianas. O programa resolve situações tridimensionais, transientes e emprega o Método dos Volumes Finitos para integrar as equações diferenciais governantes. O modelo original é baseado no modelo de Gauthier et al., 1997, tendo sido os resultados do mesmo foram usados para a validação do presente estudo. Um circuito de água quente é também projetado e apresentado para o aquecimento das estufas. A água circula através de mangueiras sobre o solo e é aquecida por um sistema de queimadores a gás liqüefeito de petróleo ou óleo combustível, transferindo assim calor para o interior da mesma. O projeto de aquecimento foi realizado através de um programa de parceria entre a Ufrgs, Sebrae, Fapergs e a Agropecuária Clarice.

Otimização de geometrias aerodinâmicas utilizando métodos inversos

O objetivo deste trabalho é a obtenção de uma técnica para a modelagem otimizada de corpos submetidos a fluxos de alta velocidade, como aerofólios em escoamentos transônicos e outras geometrias aerodinâmicas. A técnica é desenvolvida através de expansões em séries de Fourier para um conjunto de equações diferenciais com interrelação com as condições de contorno, sendo uma equação para a parte superior e outra para a parte inferior do aerofólio. O método de integração temporal empregado baseia-se no esquema explícito de Runge-Kutta de 5 estágios para as equações da quantidade de movimento e na relação de estado para a pressão. Para a aproximação espacial adota-se um esquema em volumes finitos no arranjo co-localizado em diferenças centrais. Utiliza-se dissipação artificial para amortecer as frequências de alta ordem do erro na solução das equações linearizadas. A obra apresenta a solução de escoamentos bi e tridimensionais de fluidos compressíveis transônicos em torno de perfis aerodinâmicos. Os testes num´ericos são realizados para as geometrias do NACA 0012 e 0009 e asas tridimensionais usando as equações de Euler, para número de Mach igual a 0.8 e ® = 0o. Os resultados encontrados comparam favoravelmente com os dados experimentais e numéricos disponíveis na literatura.

Aplicações do método dos elementos discretos em estruturas de concreto

A ciência moderna apresentou significativo avanço a partir do desenvolvimento da análise diferencial. A transformação de equações diferenciais de alta ordem em sistemas de equações algébricas foi possível através do desenvolvimento de métodos numéricos, constituindo este, outro grande avanço. Dentro desses pode-se destacar os métodos de diferenças finitas, dos elementos finitos, dos elementos discretos e mais recentemente, os elementos de contorno. Neste trabalho, faz-se uma contribuição ao desenvolvimento do Método dos Elementos Discretos para aplicações na Mecânica do Contínuo, na Mecânica da Fratura, assim como na determinação do dano em elementos estruturais submetidos a cargas. Neste método, a discretização espacial no modelo se realiza mediante um conjunto de massas ligadas entre se por forças materializadas como um arranjo de barras de treliça com rigidez equivalente ao contínuo que se quer representar, e mediante um esquema de integração explícita, se realiza a integração das equações de movimento no tempo. Verifica-se a validade e a capacidade do método em predizer o efeito de tamanho em elementos de concreto e concreto armado, obtendo-se uma excelente correlação com ensaios encontrados na literatura técnica, além de importantes conclusões a respeito da aplicação de cargas estáticas e dinâmicas, tanto em padrões de fissuração ou ruptura, quanto aos valores limites de resistência dos materiais ou cargas aplicadas, dando-se importância na geração aleatória das propriedades dos materiais mediante o uso do Método de Representação Espectral.