Sistema de iluminação eficiente utilizando a tecnologia LED para espaços públicos interiores

Dissertação para obtenção do Grau de Mestre em Engenharia Electrotécnica e de Computadores

Modelação numérica de um dispositivo pontual de aproveitamento de energia das ondas

Dissertação para obtenção do Grau de Mestre em Engenharia Mecânica

Plano detalhado de chã de caldeiras na Ilha do Fogo – Cabo Verde

Dissertação para obtenção do Grau Mestre em Engenharia Civil – Perfil de Construção

Modelação automática 2D de mineralizações filonianas – aplicação às Minas da Panasqueira

Neste trabalho apresenta-se uma metodologia protótipo destinada a efectuar a modelação automática a 2D da morfologia de mineralizações filonianas. Em primeiro lugar procede-se à estimação do número de ocorrências de filões em cada bloco, de uma malha de blocos, em que se subdividiu o volume em estudo. Esta intensidade de ocorrências é quantificada por uma variável que representa o número de filões por metro linear(NFM) intersectados por uma sondagem vertical, e constitui o objectivo de intensidade a atingir. Seguidamente são gerados por simulação, linhas poligonais ou arcos correspondentes ao posicionamento dos filões no perfil. Estes arcos ligam alguns pares de intersecções de filões amostrados nas sondagens e são escolhidos ao acaso segundo regras de orientação e distância. Seguidamente faz-se a avaliação da intensidade local de filões do modelo e, para os locais onde exista défice de filões relativamente ao objectivo, são adicionadas intersecções virtuais, ou seja, que não foram reconhecidas por sondagens. Este procedimento continua até o modelo dos filões se aproximar do objectivo previamente definido. O conjunto dos arcos em cada perfil, associado às espessuras de intersecção observadas nas sondagens, constitui um modelo morfológico dos filões em estrutura vectorial. No final faz-se a avaliação quantitativa do modelo e da incerteza. Os dados do estudo prático que motivaram o desenvolvimento da metodologia foram recolhidos no depósito mineral das Minas da Panasqueira. Os resultados obtidos no estudo mostram que a introdução de uma metodologia automática para a modelação vectorial de filões mineralizados é uma mais-valia, porque permite gerar modelos mais realistas e com melhor controlo geológico e resolução do que as abordagens clássicas de possança e acumulação, constituindo um auxiliar precioso de avaliação de reservas de minério.

Simulação numérica da interação de ondas com um quebra-mar poroso utilizando o código OpenFoam

O principal objetivo desta dissertação é explorar o potencial do código OpenFoam, no estudo de problemas de engenharia costeira, mais especificamente a interação de ondas com quebra-mares porosos. Neste sentido, inicialmente é realizada a validação do código numérico porousInterFoam, onde é estudada a dependência de malha e de modelos de turbulência. Conclui-se que os resultados numéricos obtidos apresentam uma boa correlação com os experimentais. Em seguida é feita uma comparação com outros códigos numéricos, nomeadamente o porousWaveFoam e o IHFOAM onde se verificam resultados idênticos. Por fim, é realizada uma simulação através do porousWaveFoam, da interação de uma onda regular com um quebra-mar poroso, e conclui-se, comparando os resultados obtidos com os resultados dos ensaios experimentais no Laboratório Nacional de Engenharia Civil (LNEC), que a elevação da superfície livre é bem modelada.

Simulação de estruturas meandriformes por objectos e estatísticas multiponto e avaliação da porosidade – aplicação a reservatórios siliciclásticos

Neste trabalho apresenta-se uma metodologia destinada a simular a morfologia e a porosidade de canais de areia em reservatórios fluviais. Como informação de partida utilizam-se um ou vários canais de areia de treino, em estrutura vectorial, que sejam representativos do reservatório, e leis de distribuição das dimensões largura e profundidade dos canais. Para a simulação da morfologia dos canais de areia, são calculados os ângulos azimutais de todos os segmentos de recta das linhas poligonais que constituem os canais de treino, a que se segue a codificação em classes de azimute e a determinação de estatísticas multiponto destas classes. Seguidamente, faz-se a simulação estocástica dos novos canais no volume do reservatório, ou seja, são geradas novas linhas poligonais condicionadas às estatísticas multiponto, a que se associam as dimensões largura e profundidade. Este algoritmo multiponto é muito eficiente, porque é a 1D, e contorna a problemática da variável azimute ser do tipo circular. A avaliação da porosidade é feita em simultâneo na conversão do modelo vectorial da morfologia dos canais de areia para o modelo matricial da malha de blocos do reservatório. Para cada bloco reservatório é contabilizada a proporção de não canal / canal e das fácies conforme um modelo conceptual de zonamento na secção do canal. A estimativa da porosidade média de cada bloco do reservatório é obtida pela soma do produto da proporção de cada fácies pela respectiva porosidade de referência. Para ilustrar as potencialidades da metodologia proposta apresenta-se um caso de estudo com dados sintéticos. De acordo com critérios geológicos e estratigráficos, consideraram-se quatro fácies na secção do canal que são discriminadas lateralmente e em profundidade por valores de referência de porosidade. O algoritmo de simulação morfológica capturou adequadamente a morfologia das imagens de treino e gerou várias imagens equiprováveis que depois foram convertidas em porosidade.

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.

Sintonização de controladores difusos através de técnicas de optimização não-linear com restrições

A presente dissertação tem como objetivo principal a implementação de uma arquitetura baseada em algoritmos evolutivos para a sintonização dos parâmetros do controlador PID (Proporcional-Integral-Derivativo) difuso, sendo o conceito de desempenho em malha fechada explicitamente tido em conta. A sintonização dos parâmetros do controlador difuso é realizada tendo em conta um problema de otimização com restrições, em que a função de custo a ser minimizada é descrita em termos do desempenho em malha fechada, com a dinâmica do sistema a ser aproximada por um modelo não linear. Como nas metodologias de otimização existentes, a incorporação de mecanismos de adaptação referentes às funções de pertença não é comum, na presente dissertação é tido em conta, para além da usual sintonização dos fatores de escala, a sintonização dos fatores de escala e funções de pertença em simultâneo. Os resultados experimentais realizados num sistema de referência, visam demonstrar os benefícios de incorporar as funções de pertença no processo de otimização em diferido. É também utilizado um método analítico de segunda ordem como referência, por forma a comparar o desempenho de uma abordagem de otimização global contra uma de otimização local. Finalmente é implementada uma abordagem em-linha, usando o método analítico de segunda ordem, na otimização dos fatores de escala e funções de pertença.

Modelação numérica do carregamento de uma estrutura confinada sujeita a uma onda de choque explosiva

A ocorrência de explosões em indústrias e em infraestruturas civis quer acidentais quer devido a ataques terroristas, vêm colocando à engenharia cada vez maiores desafios relativos à conceção de soluções que aumentem a proteção dessas infraestruturas contra os efeitos das explosões. Neste trabalho, utilizou-se um hidrocódigo comercial, o LS-DYNA, no tratamento e simulação de problemas envolvendo ondas de choque explosivas tendo-se realizado um estudo paramétrico do efeito da detonação de um alto explosivo no interior de um compartimento fechado, de geometria prismática retangular, em que as dimensões foram sendo variadas. Foi necessário calibrar previamente o modelo geométrico escolhido bem como os parâmetros utilizados para gerar a onda de choque explosiva, de forma a assegurar a consistência dos resultados, dado que não foi possível a realização de ensaios experimentais. Foi utilizada uma malha multimaterial Arbitrária Lagrangiana-Euleriana (MMALE), com recurso a elementos sólidos com um único ponto de integração, recorrendo ao método das diferenças finitas centrais para as integrações em relação ao tempo, tendo-se utilizado como base de calibração a literatura consultada. Os dados obtidos no estudo paramétrico foram tratados recorrendo a métodos estatísticos de regressão de forma a se obter uma relação fenomenológica para prever as pressões máximas dentro de um compartimento.

Estudo numérico do fenómeno de vibração induzida por vórtices num cilindro com 1 grau de liberdade em OpenFoam

Esta dissertação tem o intuito de validar o código numérico OpenFoam para problemas de fenómeno de vibração induzida por desprendimento de vórtices num cilindro circular com um grau de liberdade transversal ao escoamento. Na primeira parte é simulado o escoamento em torno de um cilindro circular fixo, e na segunda parte o escoamento em torno do cilindro oscilatório. O número de Reynolds é 200 para todas as simulações e o movimento do cilindro é descrito segundo um sistema massa-mola e massa-mola-amortecedor. Para o movimento do cilindro utilizou-se o solver sixDoFRigidDisplacement e para a resolução de problemas de malha móvel o solver displacementLaplacian. O código numérico utilizado para o caso do cilindro fixo foi pimpleFoam e para o caso do cilindro oscilatório o pimpleDyMFoam. As simulações foram feitas através da resolução das equações de Navier-Stokes num domínio computacional bidimensional. Na simulação do escoamento em redor do cilindro fixo, os coeficientes de arrasto e sustentação e o número de Strouhal foram calculados e comparados com resultados da literatura. Os resultados estão próximos da realidade, no entanto observa-se que a dimensão da malha é muito importante para a precisão dos mesmos. Na simulação do escoamento em torno do cilindro oscilatório foram calculados os coeficientes de arrasto e sustentação e deslocamento do cilindro, para os sistemas massa-mola e massa-mola-amortecedor. Para os dois sistemas foi apresentada e comparada a evolução destes resultados em função da velocidade reduzida e por fim foram comparados. Observa-se que o fator de amortecimento afeta as respostas do cilindro, dependo da velocidade reduzida. Os resultados obtidos em ambos os estudos foram satisfatórios e conclui-se que o código OpenFoam é uma boa ferramenta para resolver problemas com o fenómeno de vibração induzida por vórtices.

Reforço de placas de betão armado com argamassas armadas para acções de explosão

Quando uma explosão, com origem terrorista ou acidental, ocorre perto de estruturas vulneráveis, pode causar grandes estragos e perdas de vida. Existem diversos estudos experimentais e numéricos que analisaram reforços de estruturas de betão com Polímeros Reforçados com Fibras (FRP) contra o efeito de explosivos. Apesar dos bons resultados, estas soluções são pouco económicas. Existe, portanto, a necessidade de continuar a desenvolver novos métodos mais económicos, para aumentar a resistência dos vários elementos que compõem a estrutura de um edifício contra os efeitos de explosões. Nesta dissertação estudou-se a utilização de 3 diferentes tipos de reforço em placas de betão armado. Estas placas, com 2,6x2,17 m e 12 cm de espessura, originalmente foram projectadas como painéis para fachadas de edifícios, foram ensaiadas para acções de explosões normais ao seu plano. Os reforços estudados nesta dissertação podem assim ser utilizados para reforçar lajes e paredes de betão armado de fachadas de edifícios. Os 3 diferentes tipos de reforço, com 2 cm de espessura em toda a área da face que representa o interior do edifício, têm como base uma calda de cimento. A sua diferença está no tipo de armadura utilizada. Na solução A, é utilizada uma malha distendida galvanizada em aço. Nas soluções B e C são utilizadas fibras de aço contínuas. Estudos recentes desenvolveram e caracterizaram esta Calda Reforçada com Fibras Unidirecionais (CRFU) [1] [2]. Este compósito diverge face aos outros, desenvolvidos até então, por ser reforçado com uma manta não-tecida de fibras de aço contínuas e unidireccionais. Nesta dissertação, esta Calda Reforçada com Fibras Contínuas de aço é utilizada como reforço contra os efeitos de uma explosão nas soluções B e C. Na solução B, é utilizado 1% de fibras na direcção do vão, e na solução C, é utilizado 0,5 % de fibras em direcções ortogonais. Para ensaio experimental foram preparadas 4 placas de betão armado: três com os reforços já mencionados e uma sem reforço, como referência. Estas placas foram ensaiadas, simplesmente apoiadas em dois bordos paralelos, com um vão livre de 2,3 m. A placa A não foi ensaiada por rotura do modelo aquando da sua montagem. A placa C foi utilizada num ensaio preliminar com o objectivo de testar o sistema de ensaio. No primeiro ensaio desta placa, foram utilizados 2n kg de TNT a 3 m de distância, não provocando qualquer dano na placa. No segundo ensaio foram utilizados 8 kg de TNT a 2 m de distância, levando a placa à rotura. A placa B e a placa de referência foram ensaiadas com 8 kg de TNT a 3 m de distância. Os resultados mostram que os danos na placa reforçada diminuíram cerca de 40%.

Estudo numérico da interação de ondas com um quebra-mar poroso utilizando o código IHFOAM

A presente dissertação tem como principal objetivo a validação da ferramenta IHFOAM do código OpenFOAM para problemas de interação de ondas com estruturas costeiras porosas, concretamente a interação de ondas com o quebra-mar poente do porto de pesca de Albufeira. Numa primeira parte deste trabalho analisa-se o comportamento do código com malhas de diferentes tipos, mais especificamente, estuda-se os efeitos sobre o escoamento de uma malha não estruturada no meio poroso em relação a uma malha estruturada. Isto é de extrema importância, devido ao tempo dispensado na construção de uma malha estruturada com elementos rectangulares. É ainda estudada a dependência da simulação numérica com o modelo de turbulência. O estudo mostra que os efeitos dos diferentes tipos de malha na zona porosa influenciam pouco o escoamento, sendo aceitável o recurso a malhas não estruturadas na zona porosa, em contraste com a escolha do modelo de turbulência que se mostra muito importante para a simulação. Numa segunda parte, são calibrados os parâmetros do meio poroso usados pelo modelo numérico. Apresenta-se o método utilizado para a calibração dos parâmetros e são estudados os seus efeitos sobre o comportamento hidráulico da estrutura. O estudo é realizado com duas alturas de ondas, uma altura de onda que não gera galgamento e uma altura de onda que gera galgamento. Foram colocadas várias sondas numéricas ao longo de um canal de ondas que medem a elevação de superfície livre simulada. Para os ensaios com galgamento, mediu-se a altura da água galgada e o volume galgado, podendo-se assim estudar os efeitos dos parâmetros do meio poroso do quebra-mar sobre o galgamento. A calibração sequencial dos parâmetros mostrou-se muito eficaz, gerando resultados muito bons. Ao longo do estudo entende-se a forte dependência do comportamento hidrodinâmico da estrutura, dos parâmetros e a necessidade de os calibrar, para cada caso estudado. Conclui-se que o IHFOAM é muito interessante para o estudo numérico de interação de onda-estrutura porosa, obtendo resultados satisfatórios para a superfície livre dentro e fora do quebra-mar. O galgamento, sendo um processo muito rápido, foi razoavelmente simulado e captado pelo código.

Sensores “em movimento"

A domótica e a Internet das Coisas (Internet of Things – IoT) são dois conceitos que estão a ter bastante impacto na sociedade e a mudar o seu quotidiano, pelo que o utilizador é mais exigente na utilização dos seus equipamentos, pretendendo, que sejam mais simples e eficazes. A IoT é também muito utilizada na vertente da qualidade de vida (verificação da qualidade do ar, temperatura etc.) pelo que é necessário a utilização de uma rede de sensores de forma que toda a informação recolhida seja disponibilizada ao utilizador. Os sensores estão ligados a nós da rede de forma a recolherem informação. Essa rede pode ser constituída por diferentes topologias, sendo que a mais usual é a rede em Malha, pois no caso de um dos nós deixar de funcionar, a mensagem é enviada através de outros nós até chegar ao nó principal. O passo seguinte consiste em enviar toda essa informação para a Cloud, sendo que, nas soluções existentes, é necessário que o utilizador tenha um ponto fixo (gateway) com acesso à internet. Quando não é possível internet por cabo, a solução é a utilização de redes sem fios ou a utilização de, por exemplo, um cartão 3G/4G que implica o pagamento de taxas às operadoras móveis. Estas soluções implicam, na sua grande maioria uma instalação elétrica para alimentação dos nós dos sensores. O grande problema surge quando o utilizador não possui nenhum acesso à Internet (no local onde são instalados os sensores), ou no caso de não existência de nenhuma instalação elétrica para alimentação dos nós dos sensores. A solução proposta nesta dissertação consiste na utilização de um telemóvel, de um utilizador aleatório, como gateway. Assim, um utilizador comum pode ligar os seus sensores onde for necessário. Esses sensores são configurados Over-The-Air através de uma camada de aplicação do dispositivo de comunicação, como por exemplo, o Synapse. Após configuração, os sensores estão prontos a recolher toda a informação e enviá-la através da rede até ao nó principal. O nó principal é constituído pelo dispositivo de comunicação referido anteriormente, ligado a um microcontrolador (Arduino, por exemplo) que tem agregado um leitor de cartões SD e um dispositivo Bluetooth (BLE). Toda a informação recolhida pelos sensores é guardada no cartão SD até que um utilizador, com um telemóvel (Smartphone Android com a aplicação desenvolvida instalada) com o Bluetooth ligado, se aproxime do nó principal. Assim que a ligação é aceite e estabelecida, a aplicação envia a data mais recente de um sensor específico presente na sua base de dados na Cloud para o Arduino, permitindo que este apague os dados mais antigos presentes no cartão e envie, como resposta para o telemóvel, a informação mais antiga recolhida pelo sensor, atualizando assim a informação. A aplicação deste conceito pode ser útil quando não existe nenhuma ligação à internet ou quando um utilizador, por exemplo uma entidade responsável pelo meio ambiente e que seja necessário inserir sensores numa floresta, para prevenção de fogos. Assim os nós vão enviar toda a informação recolhida através da sua rede. Posteriormente, cabe ao utilizador escolher um sítio estratégico, onde saiba que irão passar indivíduos com alguma frequência, de modo a que estes recebam essa informação para o seu telemóvel.

Reservatórios turbidíticos – simulação da morfologia por estatísticas multiponto e avaliação das propriedades petrofísicas

Esta dissertação apresenta uma metodologia original para simular a morfologia e as propriedades petrofísicas de reservatórios de hidrocarbonetos em sistemas de canais turbidíticos. Estes sistemas são constituídos por complexos, ou seja, conjuntos de canais de arquitetura meandriforme, e são considerados importantes alvos para a indústria petrolífera. A simulação da morfologia divide-se em duas partes, primeiramente é simulada a trajetória do complexo e depois são simuladas as trajetórias dos canais propriamente ditos condicionadas à trajetória do complexo. O algoritmo de simulação utiliza as classes de ângulos azimutais de linhas poligonais de treino como uma variável aleatória. As trajetórias são simuladas também como linhas poligonais, condicionais a estatísticas multiponto das trajetórias de treino e a pontos de controlo. As estatísticas multiponto são organizadas em árvore, que guarda sequências de classes de orientação que ocorrem na trajetória de treino e as respetivas probabilidades de ocorrência. Para avaliar as propriedades petrofísicas, o modelo morfológico das trajetórias é convertido para uma malha de blocos de alta resolução, identificando-se, em cada bloco, a fácies preponderante de acordo com um modelo conceptual de zonamento da secção dos canais. A conversão prioriza os canais mais recentes (do topo) sobre os mais antigos (da base). A cada fácies é associada uma lei de distribuição da porosidade e permeabilidade, assim são geradas imagens destas propriedades petrofísicas por Simulação Sequencial Direta com histogramas locais. Finalmente, o número de blocos da malha é reduzido por upscaling e as simulações são ordenadas para poderem ser utilizadas nos simuladores de fluxo. Para ilustrar a metodologia, utilizaram-se imagens de sísmica 3D de um reservatório turbidítico na Bacia do Baixo Congo para extrair leis de distribuição das dimensões dos canais e trajetórias de treino. Os resultados representam corretamente a arquitetura complexa destes sistemas.

Estudo do galgamento em estruturas marítimas verticais utilizando o modelo numérico AMAZON

O galgamento oceânico é um fenómeno associado à interação de ondas com estruturas marítimas. Nesta dissertação estudou-se o galgamento em estruturas marítimas verticais através da aplicação do modelo numérico AMAZON, que é um modelo bidimensional, baseado na resolução das equações não lineares de águas pouco profundas. Para validar a sua aplicação a estruturas marítimas verticais impermeáveis, os resultados de caudal médio galgado obtidos com o modelo AMAZON foram comparados com os resultados de ensaios em modelo físico bidimensional realizados anteriormente e disponíveis na base de dados do projeto europeu CLASH. Os resultados foram ainda comparados com os das fórmulas empíricas e semi-empíricas de Van der Meer & Bruce de 2014, de Goda de 2008 e de Hedges & Reis de 2015. Antes da sua validação, o modelo foi calibrado utilizando um dos ensaios do CLASH. Para tal, realizou-se uma análise de sensibilidade a alguns parâmetros do modelo, concretamente da dimensão mínima da malha, da escala de aplicação do AMAZOM (à escala reduzida ou à do protótipo), da forma de modelação da parede vertical, já que, dado este modelo ser integrado na vertical, não é possível considerar paredes verticais sem utilizar um artifício de cálculo, dos valores de alguns dos parâmetros de entrada do modelo (profundidade mínima molhada e a profundidade mínima de fricção) e da dimensão do domínio de cálculo. Os resultados da calibração do AMAZON mostram que, por causa da precisão do modelo para os menores caudais médios de galgamento, se deve aplicar o modelo à escala do protótipo. Para estruturas verticais, a aproximação que o modelo faz, por defeito, para a secção vertical da estrutura através de um talude fictício a 10:1 apresenta melhores resultados que a utilização de um talude a 20:1. Contudo, ambos requerem a adoção de uma dimensão da malha bastante reduzida nesta zona e na sua vizinhança (0.03 m). No caso em estudo, os dois parâmetros de entrada do modelo não têm impacto significativo nos resultados, ao contrário do domínio de cálculo que, tal como recomendado na literatura de referência, se deve iniciar (fronteira de entrada) a aproximadamente um comprimento de onda do pé da estrutura. Os resultados obtidos mostram que, em geral, o AMAZON subestima o valor medido do caudal médio galgado com erros relativos variáveis entre -91% e -49%; apenas para um caso o modelo sobrestima o caudal médio, com um erro relativo de 43%. Quando comparado com as formulações empíricas e semi-empíricas, a fórmula de Van der Meer & Bruce apresenta melhores aproximações aos caudais médios medidos que constam na base de dados do CLASH, seguida da fórmula de Goda e dos resultados do modelo AMAZON, embora estes dois últimos apresentem valores bastante próximos. Assim, uma vez que o modelo está em constante desenvolvimento, espera-se que no futuro seja uma ferramenta mais poderosa. Contudo, enquanto não se fizerem melhorias no modelo e/ou não se estender a outros ensaios a análise feita no âmbito desta dissertação, os resultados sugerem que não há vantagem na sua utilização face às formulações empíricas e semi-empíricas disponíveis, tendo também em conta o significativo tempo computacional exigido para executar o AMAZON com a dimensão da malha escolhida.