Computational analysis of quinoxalines N,N-Dioxide and ITS Derivates

A quinoxalina e seus derivativos são uma importante classe de compostos heterocíclicos, onde os elementos N, S e O substituem átomos de carbono no anel. A fórmula molecular da quinoxalina é C8H6N2, formada por dois anéis aromáticos, benzeno e pirazina. É rara em estado natural, mas a sua síntese é de fácil execução. Modificações na estrutura da quinoxalina proporcionam uma grande variedade de compostos e actividades, tais como actividades antimicrobiana, antiparasitária, antidiabética, antiproliferativa, anti-inflamatória, anticancerígena, antiglaucoma, antidepressiva apresentando antagonismo do receptor AMPA. Estes compostos também são importantes no campo industrial devido, por exemplo, ao seu poder na inibição da corrosão do metal. A química computacional, ramo natural da química teórica é um método bem desenvolvido, utilizado para representar estruturas moleculares, simulando o seu comportamento com as equações da física quântica e clássica. Existe no mercado uma grande variedade de ferramentas informaticas utilizadas na química computacional, que permitem o cálculo de energias, geometrias, frequências vibracionais, estados de transição, vias de reação, estados excitados e uma variedade de propriedades baseadas em várias funções de onda não correlacionadas e correlacionadas. Nesta medida, a sua aplicação ao estudo das quinoxalinas é importante para a determinação das suas características químicas, permitindo uma análise mais completa, em menos tempo, e com menos custos.

Efeito da razão de submersão sobre os parâmetros hidrodinâmicos de uma coluna de recirculação de líquido através de um tubo exterior

No trabalho realizado nesta tese, procedeu-se ao estudo hidrodinâmico de uma coluna de borbulhamento de recirculação externa do líquido (CREL), que permitiu ampliar o conhecimento já existente sobre este tipo de colunas. Para realizar este estudo utilizaram-se líquidos viscosos Newtonianos, nomeadamente soluções aquosas de glicerina com viscosidades entre 0,007 e 0,522 Pa.s. Para a gama de caudais de ar injectados, 1,5x10-5 até 1,35x10-4 m3/s, o ar ascendia ao longo da coluna de borbulhamento sob a forma de bolhas tubulares. Após a realização dos ensaios, verificou-se que o regime de escoamento do líquido entre bolhas tubulares variou desde o tipo laminar até transição (1,9 de escoamento na esteira das referidas bolhas variou do regime laminar até ao regime turbulento (12,1 < ReV <3297). Nos ensaios experimentais foram utilizadas duas colunas de borbulhamento com diferentes diâmetros internos de 22 e 32 mm. A razão de submersão, , seleccionada foi de 0,85. Este valor foi escolhido por ser significativamente diferente dos usados em estudos anteriores (0,935 por Santos (2004), Santos (2005) e 0,91 por Rangel (2008)) de forma a poder analisar a influência da razão de submersão sobre os parâmetros hidrodinâmicos. Nos ensaios realizados foram determinadas a retenção do gás na coluna e a velocidade superficial do líquido em escoamento no interior da coluna de borbulhamento. Os valores da retenção do gás foram comparados com valores teóricos previstos por equações empíricas de Nicklin et al (1962). Estes valores também foram comparados com valores obtidos experimentalmente em trabalhos anteriores, onde se utilizaram razões de submersão de 0,935 e 0,91. Após análise dos valores da retenção verificou-se que estes eram mais elevados para as soluções de maior viscosidade e menor diâmetro da coluna. Para as experiências realizadas com as soluções com viscosidades entre 0,52 Pa.s e 0,053 Pa.s, os valores experimentais eram inferiores aos previstos por Nicklin et al (1962) para regime laminar tanto na coluna de 22 mm como na coluna de 32 mm. Para a solução com 0,007 Pa.s de viscosidade os valores experimentais situavam-se na zona de transição entre o regime laminar e o regime turbulento em ambas as colunas. Verificou-se que a diminuição da razão de submersão de 0,935 para 0,85 aumentou significativamente a retenção em todas as experiências realizadas. A velocidade superficial média de circulação do líquido variava na razão inversa do diâmetro e da viscosidade. A diminuição da retenção provocou uma diminuição significativa na velocidade do líquido.

Um novo gerador central de padrões de locomoção para geração de trajectórias em tempo real de robôs quadrúpedes

Neste trabalho estuda-se a geração de trajectórias em tempo real de um robô quadrúpede. As trajectórias podem dividir-se em duas componentes: rítmica e discreta. A componente rítmica das trajectórias é modelada por uma rede de oito osciladores acoplados, com simetria 4 2 Z  Z . Cada oscilador é modelado matematicamente por um sistema de Equações Diferenciais Ordinárias. A referida rede foi proposta por Golubitsky, Stewart, Buono e Collins (1999, 2000), para gerar os passos locomotores de animais quadrúpedes. O trabalho constitui a primeira aplicação desta rede à geração de trajectórias de robôs quadrúpedes. A derivação deste modelo baseia-se na biologia, onde se crê que Geradores Centrais de Padrões de locomoção (CPGs), constituídos por redes neuronais, geram os ritmos associados aos passos locomotores dos animais. O modelo proposto gera soluções periódicas identificadas com os padrões locomotores quadrúpedes, como o andar, o saltar, o galopar, entre outros. A componente discreta das trajectórias dos robôs usa-se para ajustar a parte rítmica das trajectórias. Este tipo de abordagem é útil no controlo da locomoção em terrenos irregulares, em locomoção guiada (por exemplo, mover as pernas enquanto desempenha tarefas discretas para colocar as pernas em localizações específicas) e em percussão. Simulou-se numericamente o modelo de CPG usando o oscilador de Hopf para modelar a parte rítmica do movimento e um modelo inspirado no modelo VITE para modelar a parte discreta do movimento. Variou-se o parâmetro g e mediram-se a amplitude e a frequência das soluções periódicas identificadas com o passo locomotor quadrúpede Trot, para variação deste parâmetro. A parte discreta foi inserida na parte rítmica de duas formas distintas: (a) como um offset, (b) somada às equações que geram a parte rítmica. Os resultados obtidos para o caso (a), revelam que a amplitude e a frequência se mantêm constantes em função de g. Os resultados obtidos para o caso (b) revelam que a amplitude e a frequência aumentam até um determinado valor de g e depois diminuem à medida que o g aumenta, numa curva quase sinusoidal. A variação da amplitude das soluções periódicas traduz-se numa variação directamente proporcional na extensão do movimento do robô. A velocidade da locomoção do robô varia com a frequência das soluções periódicas, que são identificadas com passos locomotores quadrúpedes.

Robôs quadrúpedes: geração de trajectórias em tempo real usando sistemas dinâmicos não-lineares

A geração de trajectórias de robôs em tempo real é uma tarefa muito complexa, não existindo ainda um algoritmo que a permita resolver de forma eficaz. De facto, há controladores eficientes para trajectórias previamente definidas, todavia, a adaptação a variações imprevisíveis, como sendo terrenos irregulares ou obstáculos, constitui ainda um problema em aberto na geração de trajectórias em tempo real de robôs. Neste trabalho apresentam-se modelos de geradores centrais de padrões de locomoção (CPGs), inspirados na biologia, que geram os ritmos locomotores num robô quadrúpede. Os CPGs são modelados matematicamente por sistemas acoplados de células (ou neurónios), sendo a dinâmica de cada célula dada por um sistema de equações diferenciais ordinárias não lineares. Assume-se que as trajectórias dos robôs são constituídas por esta parte rítmica e por uma parte discreta. A parte discreta pode ser embebida na parte rítmica, (a.1) como um offset ou (a.2) adicionada às expressões rítmicas, ou (b) pode ser calculada independentemente e adicionada exactamente antes do envio dos sinais para as articulações do robô. A parte discreta permite inserir no passo locomotor uma perturbação, que poderá estar associada à locomoção em terrenos irregulares ou à existência de obstáculos na trajectória do robô. Para se proceder á análise do sistema com parte discreta, será variado o parâmetro g. O parâmetro g, presente nas equações da parte discreta, representa o offset do sinal após a inclusão da parte discreta. Revê-se a teoria de bifurcação e simetria que permite a classificação das soluções periódicas produzidas pelos modelos de CPGs com passos locomotores quadrúpedes. Nas simulações numéricas, usam-se as equações de Morris-Lecar e o oscilador de Hopf como modelos da dinâmica interna de cada célula para a parte rítmica. A parte discreta é modelada por um sistema inspirado no modelo VITE. Medem-se a amplitude e a frequência de dois passos locomotores para variação do parâmetro g, no intervalo [-5;5]. Consideram-se duas formas distintas de incluir a parte discreta na parte rítmica: (a) como um (a.1) offset ou (a.2) somada nas expressões que modelam a parte rítmica, e (b) somada ao sinal da parte rítmica antes de ser enviado às articulações do robô. No caso (a.1), considerando o oscilador de Hopf como dinâmica interna das células, verifica-se que a amplitude e frequência se mantêm constantes para -5de Hopf, a amplitude e a frequência têm o mesmo comportamento, crescendo e diminuindo nos intervalos de g [-0.5,0.34] e [0.4,1.83], sendo nos restantes valores de g nulas. Isto traduz-se em variações na extensão do movimento e na velocidade do robô, proporcionais à amplitude e à frequência, respectivamente. Ainda com o oscilador Hopf, no caso (b), a frequência mantêm-se constante enquanto a amplitude diminui para g<0.2 e aumenta para g>0.2. A extensão do movimento varia de forma directamente proporcional à amplitude. No caso das equações de Morris-Lecar, quando a componente discreta é embebida (a.2), a amplitude e a frequência aumentam e depois diminuem para - 0.17de diminui para g<0.5 e aumenta para g>0.5 Pode concluir-se que: (1) a melhor forma de inserção da parte discreta que menos perturbação insere no robô é a inserção como offset; (2) a inserção da parte discreta parece ser independente do sistema de equações diferenciais ordinárias que modelam a dinâmica interna de cada célula. Como trabalho futuro, é importante prosseguir o estudo das diferentes formas de inserção da parte discreta na parte rítmica do movimento, para que se possa gerar uma locomoção quadrúpede, robusta, flexível, com objectivos, em terrenos irregulares, modelada por correcções discretas aos padrões rítmicos.

Simulação e controlo de veículos autónomos terrestres

A implementação e venda de robôs autónomos tem sido um sector que nos últimos anos tem adquirido cada vez mais quota no mercado, nomeadamente no sector militar, agrícola e da vigilância. Como tal, tem sido também de grande importância a capacidade de implementar e testar robôs por parte das entidades que os fabricam. Uma das formas que tem garantido o sucesso do desenvolvimento de robôs é a simulação prévia dos mesmos antes que estes passem a fase de produção. Sendo assim, o LSA como entidade de desenvolvimento de robôs autónomos, tem necessidade de adquirir um sistema que simule os robôs em desenvolvimento. O trabalho desta tese consiste na realização de um sistema que simule robôs autónomos terrestres de forma que se possa observar o comportamento da cinemática, dinânica e hardware dos robôs em ambiente 3D. Esta aplicação de simulação pode mais tarde ser utilizada pelo laboratório para testar missões, validar alterações de estrutura, sensores, etc. Para além disso, com recurso ao simulador Player/Stage/Gazebo testar o robô LINCE e implementar algoritmos de controlo para o mesmo. Os algoritmos de controlo implementados baseiam-se em primitivas de controlo básico para serem utilizadas pelo sistema de navegação e gerar trajectórias complexas. Os algoritmos desenvolvidos nesta tese baseiam-se nas equações cinemáticas do veículo estudado. Estes algoritmos depois de testados no simulador, poderão ser colocados no Hardware do robô. Desta forma consegue-se desenvolver algoritmos para determinado robô sem que este esteja operacional.

Abordagem prática para execução de taludes em obras rodoviárias com recurso a pré‐corte

Este trabalho é realizado no domínio das obras de engenharia, área onde o desmonte de rocha com recurso a explosivos em obras rodoviárias é uma actividade específica e consistiu no acompanhamento e execução de três obras rodoviárias de média e grande dimensão. A necessidade de executar escavações, recorrendo a técnicas de desmonte cuidadoso de contorno, onde o plano de corte do talude final deve obedecer a requisitos de localização, alinhamento, inclinação, estabilidade e também estéticos, acrescendo a isto a necessidade de optimizar os meios envolvidos, obriga a que esta actividade seja encarada de uma forma sistematizada, visando o racional aproveitamento de recursos. A execução desta actividade requer conhecimentos no domínio das técnicas de desmonte de contorno, dos explosivos, do mecanismo de rotura de rochas, da operação de perfuração e da geomecânica dos maciços. A abordagem deste trabalho incide sobre a técnica denominada de pré‐corte e tem como objectivo encontrar uma equação característica que permita relacionar diferentes parâmetros envolvidos nesta actividade. Este objectivo é alcançado recorrendo à correlação entre equações relativas à pressão de detonação, à pressão no furo e ao espaçamento entre furos consecutivos, desenvolvidas por outros autores. Desta forma obteve‐se uma equação que relaciona parâmetros relativos ao maciço rochoso (resistência à tracção), ao explosivo (velocidade de detonação e densidade) e ao diagrama de fogo (concentração de carga – volume de explosivo e comprimento do furo – volume do furo). A comparação entre os valores destes parâmetros obtidos na produção e os obtidos com recurso à equação característica permite concluir que a sua aplicação para execução de futuras obras possibilita uma optimização dos meios envolvidos.

Desenvolvimento de uma ferramenta para cálculo de estruturas em material compósito

Os desafios à engenharia moderna são cada vez maiores, pretendendo-se quase sempre obter estruturas mais leves, com propriedades mecânicas atrativas e muitas vezes com geometrias complexas. Com tais requisitos, um dos materiais que tem vindo a ter uma crescente aplicação é o material compósito. Contudo, no que toca ao cálculo estrutural destes materiais, tudo se torna mais complexo, já que são materiais que geralmente são formados por empilhamento de várias camadas de material heterogéneo, podendo estas encontrarem-se dispostas segundo diferentes orientações. Assim, a utilização de um software que permita a previsão das propriedades mecânicas de uma estrutura em material compósito através da micromecânica, a aplicação da Teoria Clássica dos Laminados e de um critério de rotura, como por exemplo o de Tsai-Hill, é fundamental para agilizar o processo de estudo da estrutura a fabricar. Para dar uma resposta a tal necessidade foi desenvolvida uma aplicação, em MATLAB® GUI, denominada CAFE – Composite Analysis For Engineers, com ambiente gráfico apelativo, que permite determinar todas as variáveis importantes no estudo de estruturas em material compósito. Esta aplicação visa suportar e agilizar a aprendizagem desta área do conhecimento, permitindo também o acesso ao código de cálculo por parte do utilizador, de modo a conhecerem-se as equações utilizadas e, eventualmente, ser alvo de futuros desenvolvimentos. O programa desenvolvido foi alvo de validação, recorrendo-se para tal, a uma comparação dos resultados obtidos entre o respetivo programa e por um outro programa de grande fiabilidade. Assim sendo, concluiu-se que o software CAFE apresenta resultados válidos, encontrando-se apto a ser utilizado.

Dimensionamento de um sistema UPS flywheel

Nesta dissertação, numa primeira fase, é efetuado um estudo ao estado da arte da qualidade de energia elétrica e do volante de inércia mais conhecido por flywheel. Como primeiro tópico da primeira fase é efetuado um estudo das perturbações possíveis, baseado na sua forma de onda, são representadas graficamente as diversas perturbações, é efetuada uma análise das possíveis causas para cada perturbação e por fim, também são apresentadas as possíveis soluções para as mesmas perturbações. Como possíveis soluções, são apresentados os vários sistemas de armazenamento nomeadamente, UPS (fonte de alimentação ininterrupta), volantes de inércia, entre outros sistemas. No segundo tópico é efetuada uma introdução ao sistema flywheel, apresentando alguns aspetos importantes, tais como as várias vantagens do seu uso como sistema de armazenamento de energia, diferentes tipos de aplicações, sua constituição (formas, tipos de matérias, máquina elétrica, entre outros) e são apresentados também os diferentes tipos de aplicações, utilizando‐a como sistema de armazenamento. A segunda fase é dedicada às equações fundamentais para o desenvolvimento duma flywheel. Neste caso são apresentadas várias equações para o cálculo da energia armazenada, momento de inércia, tensões, entre muitas outras equações essenciais. Na terceira fase é efetuada uma comparação entre diversos sistemas de armazenamento de energia, usando tanto a tecnologia flywheel como baterias tradicionais. Esta comparação foi realizada em três categorias, 150 kVA, 600 kVA e 1200 kVA. Em cada categoria é realizada uma comparação entre diversos fatores tais como: peso, custos, dimensões, limitações, regulações possíveis e questões ambientais. Como quarta fase, é apresentado um modelo de uma draga utilizando o sistema flywheel como armazenamento de energia. Neste modelo é apresentada a flywheel como solução, uma vez que armazena a energia recuperada dos processos eletromecânicos que necessitam de absorver grandes quantidades de energia. Desta forma, consegue‐se uma redução no consumo de energia. Os resultados da simulação mostram uma redução de 25% de energia e uma redução entre 80% a 90% em todas as categorias de emissões.

Modelagem e simulação da dinâmica do arrastamento por ar

Os compostos orgânicos voláteis constituem uma fonte vulgar de contaminação da água subterrânea, a qual pode ser eliminada pela tecnologia do arrastamento por ar (air stripping) em colunas com enchimento desordenado e utilizando fluxos das fases em contra-corrente. Propõe-se neste trabalho uma nova metodologia de dimensionamento destas colunas, para qualquer tipo de enchimento e de contaminante, onde não há necessidade de se arbitrar nenhum diâmetro, onde se evita o recurso a ábacos experimentais e onde o regime hidráulico conveniente é seleccionado à partida. O procedimento proposto foi algoritmizado e convertido num programa em linguagem C++. Para verificar e testar não só o dimensionamento mas também o comportamento teórico estacionário e dinâmico construiu-se de raiz uma coluna experimental. Seleccionou-se como contaminante uma solução de clorofórmio em água destilada. A experimentação permite, ainda, corrigir o coeficiente de transferência de massa global teórico estimado pelas correlações de Onda e que depende de inúmeros parâmetros nem sempre controláveis experimentalmente. Apresenta-se, em seguida, um modelo original de simulação dinâmica do comportamento da coluna e que é constituído por um sistema de equações diferenciais não lineares (parâmetros distribuidos). No entanto, se os débitos forem arbitrados como constantes, o sistema passa a ser linear apesar de não possuir solução analítica evidente (p.e. por transformações integrais). A discretização por diferenças finitas permitiu superar estas dificuldades. Existe uma notável concordância entre os valores experimentais e os previstos no modelo.

A modelação nas tecnologias de tratamento de água – o caso particular do arrastamento por ar

Os recursos hídricos subterrâneos constituem a maior fonte disponível de água doce e são, cada vez mais, utilizados no abastecimento de água para consumo humano e em múltiplas actividades industriais e agrícolas. A poluição destes recursos, fruto de acidente ou negligência, é uma realidade que se verifica com uma frequência alarmante e que atinge níveis preocupantes no nosso país. Nestas Jornadas pretende-se apresentar, numa primeira fase, uma breve alusão à operação de arrastamento por ar em colunas com enchimento desordenado, e suas vantagens face a outros processos de despoluição e, numa segunda fase, a utilidade de desenvolvimento de modelos de simulação e optimização numérica desta operação unitária. Para tal, é exposta a formulação do modelo matemático fenomenológico concebido por Castro e Fiúza (Castro, 1997) desenvolvido sob uma fundamentação teórica com carácter sistémico, em termos de equações de balanço de massa para uma determinada secção de volume infinitesimal, traduzido num sistema de EDOs e num sistema de EDPs, para situações de estado estacionário e transiente respectivamente e são apresentados os resultados entretanto obtidos.

Aplicação do Método de TAU a um modelo de descontaminação de água subterrânea

Neste trabalho são apresentados alguns resultados obtidos com o método de TAU implementado sobre um modelo dinâmico de transferência de massa em colunas com enchimento desordenado. A solução do sistema de equações em derivadas parciais que traduz o modelo em estado transiente revelou-se ajustada com resultados neméricos muito interessantes.

Transferência de massa através de arrastamento por ar: resultados de simulação laboratorial e numérica

O tolueno tem uma vasta aplicação industrial e apresenta alguns riscos para a saúde humana. Legislação comunitária, recentemente publicada, exige que se cumpra a médio prazo padrões mínimos de qualidade apontando como limite máximo de concentração de tolueno na água destinada a consumo humano o valor de 1 ppm. O processo de arrastamento por ar em colunas com enchimento desordenado é uma das técnicas de remoção disponíveis, embora ainda pouco implementada no nosso país. Neste trabalho simulou-se, laboratorial e numericamente, a remoção do tolueno da água através desta tecnologia sob condições operatórias pré-estabelecidas. A simulação à escala laboratorial foi executada numa instalação piloto constituída por uma coluna de vidro acrílico com 0.137metros de diâmetro interno preenchida desordenadamente com anéis Raflux de polietileno com 15milímetros de calibre numa altura de 1.10 metros. A experimentação efectuada abrangeu ensaios realizados quer em regime permanente quer em regime transitório. As simulações numéricas do comportamento de uma coluna com idênticas características, em estado estacionário e em estado transiente, foram efectuadas com “software” específico concebido para MatLab. Este software permite a descrição do regime hidrodinâmico e da transferência de massa deste processo de acordo com os modelos de Robbins [Robbins, 1991] e Onda [Onda, 1968], respectivamente, inseridos num modelo global fenomenológico traduzido num sistema de equações às derivadas parciais de 1ª ordem, do tipo hiperbólico, com uma variável espacial axial e uma variável temporal. A resolução deste modelo é efectuada através de algoritmos e estratégias específicas nomeadamente a implementação numérica dos métodos explícito de Euler para resolução do sistema de ED0s que traduz o estado estacionário primitivo e de diferenças finitas com esquema explícito associado a manipulação de matrizes esparsas para resolução do sistema de EDPs que traduz o estado dinâmico. A concordância entre os resultados obtidos na teoria e na prática, pese embora as dificuldades encontradas no controlo analítico das amostras, permite concluir que o modelo proposto possui um rigor fenomenológico suficiente e que o modelo numérico associado possibilita um elevado grau de exactidão e de fiabilidade. A aplicação desenvolvida pode constituir uma importante ferramenta de apoio à análise dinâmica desta operação unitária à escala industrial, reduzindo o esforço experimental e os respectivos custos associados. Constitui igualmente a base para uma implementação de um potencial sistema de controlo por actuação antecipada com base num modelo que reveste a forma de um sistema de parâmetros distribuídos.

Álgebra matricial: Conceitos, Exercícios e Aplicações

Este livro pretende ser um documento onde a ligação entre a abordagem clássica da Álgebra Linear habitualmente encontrada na literatura e a Teoria de Matrizes seja apresentada de forma simples e rigorosa em simultâneo com a exposição de aplicações. Conscientes da vastidão de possíveis caminhos a seguir na apresentação das matérias inerentes à Álgebra e ao Cálculo Matricial, os autores optaram por seguir uma orientação que tivesse em linha de conta a atual tendência para a diminuição dos tempos letivos e incentivo à utilização de software MATLAB®, principalmente nos cursos de Engenharia. Neste sentido, este livro está organizado em cinco capítulos – Revisão de conceitos elementares, Cálculo matricial e determinantes, Sistemas de equações lineares, Espaços vetoriais e transformações lineares e Geometria analítica – ao longo dos quais se procurou obedecer a uma estrutura evolutiva em torno do rigor e da formalidade, mas sem excessos de nomenclatura. No final de cada capítulo, é proporcionado um conjunto de exercícios variados e não repetitivos, em número suficiente e equilibrado, apresentando-se alguns deles já resolvidos, propondo-se outros para resolução e ilustrando algumas aplicações práticas de integração de conhecimentos, recorrendo ao software MATLAB®.

Cartografia e avaliação geotécnica de blocos rochosos em plataformas costeiras (Gaia, NW Portugal)

Mestrado em Engenharia Geotécnica e Geoambiente

Modelação numérica de uma caldeira doméstica a pellets

A biomassa é uma das fontes de energia renovável com maior potencial em Portugal, sendo a capacidade de produção de pellets de biomassa atualmente instalada superior a 1 milhão de toneladas/ano. Contudo, a maioria desta produção destina-se à exportação ou à utilização em centrais térmicas a biomassa, cujo crescimento tem sido significativo nos últimos anos, prevendo-se que a capacidade instalada em 2020 seja de aproximadamente 250 MW. O mercado português de caldeiras a pellets é bastante diversificado. O estudo que realizamos permitiu concluir que cerca de 90% das caldeiras existentes no mercado português têm potências inferiores a 60 kW, possuindo na sua maioria grelha fixa (81%), com sistema de ignição eléctrica (92%) e alimentação superior do biocombustível sólido (94%). O objetivo do presente trabalho foi o desenvolvimento de um modelo para simulação de uma caldeira a pellets de biomassa, que para além de permitir otimizar o projeto e operação deste tipo de equipamento, permitisse avaliar as inovações tecnológicas nesta área. Para tal recorreu-se o BiomassGasificationFoam, um código recentemente publicado, e escrito para utilização com o OpenFOAM, uma ferramenta computacional de acesso livre, que permite a simulação dos processos de pirólise, gasificação e combustão de biomassa. Este código, que foi inicialmente desenvolvido para descrever o processo de gasificação na análise termogravimétrica de biomassa, foi por nós adaptado para considerar as reações de combustão em fase gasosa dos gases libertados durante a pirólise da biomassa (recorrendo para tal ao solver reactingFoam), e ter a possibilidade de realizar a ignição da biomassa, o que foi conseguido através de uma adaptação do código de ignição do XiFoam. O esquema de ignição da biomassa não se revelou adequado, pois verificou-se que a combustão parava sempre que a ignição era inativada, independentemente do tempo que ela estivesse ativa. Como alternativa, usaram-se outros dois esquemas para a combustão da biomassa: uma corrente de ar quente, e uma resistência de aquecimento. Ambos os esquemas funcionaram, mas nunca foi possível fazer com que a combustão fosse autossustentável. A análise dos resultados obtidos permitiu concluir que a extensão das reações de pirólise e de gasificação, que são ambas endotérmicas, é muito pequena, pelo que a quantidade de gases libertados é igualmente muito pequena, não sendo suficiente para libertar a energia necessária à combustão completa da biomassa de uma maneira sustentável. Para tentar ultrapassar esta dificuldade foram testadas várias alternativas, , que incluíram o uso de diferentes composições de biomassa, diferentes cinéticas, calores de reação, parâmetros de transferência de calor, velocidades do ar de alimentação, esquemas de resolução numérica do sistema de equações diferenciais, e diferentes parâmetros dos esquemas de resolução utilizados. Todas estas tentativas se revelaram infrutíferas. Este estudo permitiu concluir que o solver BiomassGasificationFoam, que foi desenvolvido para descrever o processo de gasificação de biomassa em meio inerte, e em que a biomassa é aquecida através de calor fornecido pelas paredes do reator, aparentemente não é adequado à descrição do processo de combustão da biomassa, em que a combustão deve ser autossustentável, e em que as reações de combustão em fase gasosa são importantes. Assim, é necessário um estudo mais aprofundado que permita adaptar este código à simulação do processo de combustão de sólidos porosos em leito fixo.