Optimal control and numerical optimization applied to epidemiological models

A relação entre a epidemiologia, a modelação matemática e as ferramentas computacionais permite construir e testar teorias sobre o desenvolvimento e combate de uma doença. Esta tese tem como motivação o estudo de modelos epidemiológicos aplicados a doenças infeciosas numa perspetiva de Controlo Ótimo, dando particular relevância ao Dengue. Sendo uma doença tropical e subtropical transmitida por mosquitos, afecta cerca de 100 milhões de pessoas por ano, e é considerada pela Organização Mundial de Saúde como uma grande preocupação para a saúde pública. Os modelos matemáticos desenvolvidos e testados neste trabalho, baseiam-se em equações diferenciais ordinárias que descrevem a dinâmica subjacente à doença nomeadamente a interação entre humanos e mosquitos. É feito um estudo analítico dos mesmos relativamente aos pontos de equilíbrio, sua estabilidade e número básico de reprodução. A propagação do Dengue pode ser atenuada através de medidas de controlo do vetor transmissor, tais como o uso de inseticidas específicos e campanhas educacionais. Como o desenvolvimento de uma potencial vacina tem sido uma aposta mundial recente, são propostos modelos baseados na simulação de um hipotético processo de vacinação numa população. Tendo por base a teoria de Controlo Ótimo, são analisadas as estratégias ótimas para o uso destes controlos e respetivas repercussões na redução/erradicação da doença aquando de um surto na população, considerando uma abordagem bioeconómica. Os problemas formulados são resolvidos numericamente usando métodos diretos e indiretos. Os primeiros discretizam o problema reformulando-o num problema de optimização não linear. Os métodos indiretos usam o Princípio do Máximo de Pontryagin como condição necessária para encontrar a curva ótima para o respetivo controlo. Nestas duas estratégias utilizam-se vários pacotes de software numérico. Ao longo deste trabalho, houve sempre um compromisso entre o realismo dos modelos epidemiológicos e a sua tratabilidade em termos matemáticos.

Effects of the climatic factors and anthropogenic actions in the Ria de Aveiro

A compreensão dos impactes das alterações climáticas é fundamental para a gestão a longo do prazo dos ecossistemas estuarinos. Esta compreensão só poderá ser efectiva considerando a variabilidade climática natural e o papel relativo das intervenções antropogénicas nestes ecossistemas. Assim, a presente dissertação analisa a influência das alterações climáticas e pressões antropogénicas na qualidade da água e dinâmica ecológica da Ria de Aveiro com base numa abordagem integrada, que combinou a análise de séries temporais dos últimos 25 anos e a modelação numérica de elevada resolução de cenários futuros de alterações climáticas e intervenções antropogénicas. A componente de modelação de qualidade da água e ecológica foi melhorada a vários níveis. A análise de sensibilidade do modelo 3D hidrodinâmicoecológico ECO-SELFE aplicado à Ria de Aveiro e a revisão das constantes de semi-saturação para absorção de nutrientes pelo fitoplâncton contribuíram para a precisão e robustez das aplicações. A concentração do fitoplâncton foi significativamente influenciada pelas taxas de crescimento do fitoplâncton e de mortalidade e excreção do zooplâncton, e apresentou uma sensibilidade reduzida à variação das constantes de semi-saturação na gama identificada para as diatomáceas. O acoplamento do ECO-SELFE a um modelo de campo próximo e a integração do ciclo do oxigénio aumentaram a sua capacidade de representação dos processos e das escalas espaciais relevantes. A validação do ECO-SELFE foi realizada com base num conjunto de campanhas específicas realizadas no canal de Mira. Os padrões espaciais e temporais observados para as várias variáveis (clorofila a, nutrientes, oxigénio dissolvido, salinidade, temperatura da água, correntes e níveis) foram simulados com erros menores ou semelhantes aos obtidos neste tipo de aplicações. A análise dos padrões de variabilidade espacial e temporal da qualidade da água e ecológica na Ria de Aveiro a diferentes escalas, efectuada com base nos dados históricos de 1985 a 2010 complementados pelas campanhas realizadas, sugeriu uma influência combinada da variabilidade climática e das acções antropogénicas. Os cenários futuros de alterações climáticas e intervenções antropogénicas simulados evidenciaram uma influência mais significativa das alterações climáticas quando comparadas com os efeitos das acções antropogénicas analisadas. As variações mais significativas são previstas para os cenários de subida do nível do mar, seguidos dos cenários de alterações dos regimes hidrológicos, evidenciando o papel da circulação (maré e caudal fluvial) no estabelecimento da qualidade da água e dinâmica ecológica na laguna. Para os cenários de subida do nível do mar são previstos decréscimos significativos da clorofila a e dos nutrientes a jusante e nas zonas intermédias do canal, e um aumento significativo da salinidade a montante. Estas alterações poderão favorecer modificações da composição e distribuição das comunidades, afectando a cadeia alimentar e causando uma progressão para montante de espécies marinhas. Os resultados sugerem ainda que os efeitos poderão ser mais significativos em estuários pouco profundos.

Development of numerical methodologies for parameter identification and shape optimization in metal forming simulations

Por parte da indústria de estampagem tem-se verificado um interesse crescente em simulações numéricas de processos de conformação de chapa, incluindo também métodos de engenharia inversa. Este facto ocorre principalmente porque as técnicas de tentativa-erro, muito usadas no passado, não são mais competitivas a nível económico. O uso de códigos de simulação é, atualmente, uma prática corrente em ambiente industrial, pois os resultados tipicamente obtidos através de códigos com base no Método dos Elementos Finitos (MEF) são bem aceites pelas comunidades industriais e científicas Na tentativa de obter campos de tensão e de deformação precisos, uma análise eficiente com o MEF necessita de dados de entrada corretos, como geometrias, malhas, leis de comportamento não-lineares, carregamentos, leis de atrito, etc.. Com o objetivo de ultrapassar estas dificuldades podem ser considerados os problemas inversos. No trabalho apresentado, os seguintes problemas inversos, em Mecânica computacional, são apresentados e analisados: (i) problemas de identificação de parâmetros, que se referem à determinação de parâmetros de entrada que serão posteriormente usados em modelos constitutivos nas simulações numéricas e (ii) problemas de definição geométrica inicial de chapas e ferramentas, nos quais o objetivo é determinar a forma inicial de uma chapa ou de uma ferramenta tendo em vista a obtenção de uma determinada geometria após um processo de conformação. São introduzidas e implementadas novas estratégias de otimização, as quais conduzem a parâmetros de modelos constitutivos mais precisos. O objetivo destas estratégias é tirar vantagem das potencialidades de cada algoritmo e melhorar a eficiência geral dos métodos clássicos de otimização, os quais são baseados em processos de apenas um estágio. Algoritmos determinísticos, algoritmos inspirados em processos evolucionários ou mesmo a combinação destes dois são usados nas estratégias propostas. Estratégias de cascata, paralelas e híbridas são apresentadas em detalhe, sendo que as estratégias híbridas consistem na combinação de estratégias em cascata e paralelas. São apresentados e analisados dois métodos distintos para a avaliação da função objetivo em processos de identificação de parâmetros. Os métodos considerados são uma análise com um ponto único ou uma análise com elementos finitos. A avaliação com base num único ponto caracteriza uma quantidade infinitesimal de material sujeito a uma determinada história de deformação. Por outro lado, na análise através de elementos finitos, o modelo constitutivo é implementado e considerado para cada ponto de integração. Problemas inversos são apresentados e descritos, como por exemplo, a definição geométrica de chapas e ferramentas. Considerando o caso da otimização da forma inicial de uma chapa metálica a definição da forma inicial de uma chapa para a conformação de um elemento de cárter é considerado como problema em estudo. Ainda neste âmbito, um estudo sobre a influência da definição geométrica inicial da chapa no processo de otimização é efetuado. Este estudo é realizado considerando a formulação de NURBS na definição da face superior da chapa metálica, face cuja geometria será alterada durante o processo de conformação plástica. No caso dos processos de otimização de ferramentas, um processo de forjamento a dois estágios é apresentado. Com o objetivo de obter um cilindro perfeito após o forjamento, dois métodos distintos são considerados. No primeiro, a forma inicial do cilindro é otimizada e no outro a forma da ferramenta do primeiro estágio de conformação é otimizada. Para parametrizar a superfície livre do cilindro são utilizados diferentes métodos. Para a definição da ferramenta são também utilizados diferentes parametrizações. As estratégias de otimização propostas neste trabalho resolvem eficientemente problemas de otimização para a indústria de conformação metálica.

Design and optimization next generation passive all-optical networks

Este trabalho investiga novas metodologias para as redes óticas de acesso de próxima geração (NG-OAN). O trabalho está dividido em quatro tópicos de investigação: projeto da rede, modelos numéricos para efeitos não lineares da fibra ótica, impacto dos efeitos não lineares da fibra ótica e otimização da rede. A rede ótica de acesso investigada nesse trabalho está projetado para suprir os requisitos de densidade de utilizadores e cobertura, isto é, suportar muitos utilizadores ( 1000) com altas velocidades de conexão dedicada ( 1 Gb/s) ocupando uma faixa estreita do espectro ( 25 nm) e comprimentos de fibra ótica até 100 km. Os cenários são baseados em redes óticas passivas com multiplexagem por divisão no comprimento de onda de alta densidade (UDWDM-PON) utilizando transmissores/receptores coerentes nos terminais da rede. A rede é avaliada para vários ritmos de transmissão usando formatos de modulação avançados, requisitos de largura de banda por utilizador e partilha de banda com tecnologias tradicionais de redes óticas passivas (PON). Modelos numéricos baseados em funções de transferência das séries de Volterra (VSTF) são demonstrados tanto para a análise dos efeitos não lineares da fibra ótica quanto para avaliação do desempenho total da rede. São apresentadas as faixas de potência e distância de transmissão nas quais as séries de Volterra apresentam resultados semelhantes ao modelo referência Split-Step Fourier (SSF) (validado experimentalmente) para o desempenho total da rede. Além disso, um algoritmo, que evita componentes espectrais com intensidade nulo, é proposto para realizar cálculos rápidos das séries. O modelo VSTF é estendido para identificar unicamente os efeitos não lineares da fibra ótica mais relevantes no cenário investigado: Self-Phase Modulation (SPM), Cross-Phase Modulation (XPM) e Four-Wave Mixing (FWM). Simulações numéricas são apresentadas para identificar o impacto isolado de cada efeito não linear da fibra ótica, SPM, XPM e FWM, no desempenho da rede com detecção coerente UDWDM-PON, transportando canais com modulação digital em fase (M-ária PSK) ou modulação digital em amplitude (M-ária QAM). A análise numérica é estendida para diferentes comprimentos de fibra ótica mono modo (SSMF), potência por canal e ritmo de transmissão por canal. Por conseguinte, expressões analíticas são extrapoladas para determinar a evolução do SPM, XPM e FWM em função da potência e distância de transmissão em cenários NG-OAN. O desempenho da rede é otimizada através da minimização parcial da interferência FWM (via espaçamento desigual dos canais), que nesse caso, é o efeito não linear da fibra ótica mais relevante. Direções para melhorias adicionas no desempenho da rede são apresentados para cenários em que o XPM é relevante, isto é, redes transportando formatos de modulação QAM. A solução, nesse caso, é baseada na utilização de técnicas de processamento digital do sinal.

Injection moulding process optimization for microstructured parts production

Ao longo das últimas décadas, a micromoldação (u-moldação) por injeção de termoplásticos ganhou um lugar de destaque no mercado de equipamentos eletrónicos e de uma ampla gama de componentes mecânicos. No entanto, quando o tamanho do componente diminui, os pressupostos geralmente aceites na moldação por injeção convencional deixam de ser válidos para descrever o comportamento reológico e termomecânico do polímero na microimpressão. Por isso, a compreensão do comportamento dinâmico do polímero à escala micro bem como da sua caraterização, análise e previsão das propriedades mecânicas exige uma investigação mais alargada. O objetivo principal deste programa doutoral passa por uma melhor compreensão do fenómeno físico intrínseco ao processo da μ-moldação por injeção. Para cumprir com o objetivo estabelecido, foi efetuado um estudo paramétrico do processo de μ-moldação por injeção, cujos resultados foram comparados com os resultados obtidos por simulação numérica. A caracterização dinâmica mecânica das μ-peças foi efetuada com o objetivo de recolher os dados necessários para a previsão do desempenho mecânico das mesmas, a longo prazo. Finalmente, depois da calibração do modelo matemático do polímero, foram realizadas análises estruturais com o intuito de prever o desempenho mecânico das μ-peças no longo prazo. Verificou-se que o desempenho mecânico das μ-peças pode ser significativamente afetado pelas tensões residuais de origem mecânica e térmica. Estas últimas, resultantes do processo de fabrico e das condições de processamento, por isso, devem ser consideradas na previsão do desempenho mecânico e do tempo de serviço das u-moldações.

Modeling the Western Iberian margin circulation: present and future

The purpose of this work is to carry out a comprehensive study on the Western Iberian Margin (WIM) circulation my means of numerical modeling, and to postulate what this circulation will be in the future. The adopted approach was the development of a regional ocean model configuration with high resolution, capable of reproducing the largeand small-scale dynamics of the coastal transition zone. Four numerical experiences were carried out according to these objectives: (1) a climatological run, in order to study the system’s seasonal behavior and its mean state; (2) a run forced with real winds and fluxes for period 2001-2011 in order to study the interannual variability of the system; (3) a run forced with mean fields from Global Climate Models (GCMs) for the present, in order to validate GCMs as adequate forcing for regional ocean modeling; (4) a similar run (3) for period 2071-2100, in order to assess possible consequences of a future climate scenario on the hydrography and dynamics of the WIM. Furthermore, two Lagrangian particle studies were carried out: one in order to trace the origin of the upwelled waters along the WIM; the other in order to portrait the patterns of larval dispersal, accumulation and connectivity. The numerical configuration proved to be adequate in the reproduction of the system’s mean state, seasonal characterization and an interannual variability study. There is prevalence of poleward flow at the slope, which coexists with the upwelling jet during summer, although there is evidence of its shifting offshore, and which is associated with the Mediterranean Water flow at deeper levels, suggesting a barotropic character. From the future climate scenario essay, the following conclusions were drawn: there is general warming and freshening of upper level waters; there is still poleward tendency, and despite the upwellingfavorable winds strengthening in summer the respective coastal band becomes more restricted in width and depth. In what concerns larval connectivity and dispersion along the WIM, diel vertical migration was observed to increase recruitment throughout the domain, and while smooth coastlines are better suppliers, there is higher accumulation where the topography is rougher.

Modelling the Minho river plume intrusion into the rias baixas

The Minho River, situated 30 km south of the Rias Baixas is the most important freshwater source flowing into the Western Galician Coast (NW of the Iberian Peninsula). This discharge is important to determine the hydrological patterns adjacent to its mouth, particularly close to the Galician coastal region. The buoyancy generated by the Minho plume can flood the Rias Baixas for long periods, reversing the normal estuarine density gradients. Thus, it becomes important to analyse its dynamics as well as the thermohaline patterns of the areas affected by the freshwater spreading. Thus, the main aim of this work was to study the propagation of the Minho estuarine plume to the Rias Baixas, establishing the conditions in which this plume affects the circulation and hydrographic features of these coastal systems, through the development and application of the numerical model MOHID. For this purpose, the hydrographic features of the Rias Baixas mouths were studied. It was observed that at the northern mouths, due to their shallowness, the heat fluxes between the atmosphere and ocean are the major forcing, influencing the water temperature, while at the southern mouths the influence of the upwelling events and the Minho River discharge were more frequent. The salinity increases from south to north, revealing that the observed low values may be caused by the Minho River freshwater discharge. An assessment of wind data along the Galician coast was carried out, in order to evaluate the applicability of the study to the dispersal of the Minho estuarine plume. Firstly, a comparative analysis between winds obtained from land meteorological stations and offshore QuikSCAT satellite were performed. This comparison revealed that satellite data constitute a good approach to study wind induced coastal phenomena. However, since the numerical model MOHID requires wind data with high spatial and temporal resolution close to the coast, results of the forecasted model WRF were added to the previous study. The analyses revealed that the WRF model data is a consistent tool to obtain representative wind data near the coast, showing good results when comparing with in situ wind observations from oceanographic buoys. To study the influence of the Minho buoyant discharge influence on the Rias Baixas, a set of three one-way nested models was developed and implemented, using the numerical model MOHID. The first model domain is a barotropic model and includes the whole Iberian Peninsula coast. The second and third domains are baroclinic models, where the second domain is a coarse representation of the Rias Baixas and adjacent coastal area, while the third includes the same area with a higher resolution. A bi-dimensional model was also implemented in the Minho estuary, in order to quantify the flow (and its properties) that the estuary injects into the ocean. The chosen period for the Minho estuarine plume propagation validation was the spring of 1998, since a high Minho River discharge was reported, as well as favourable wind patterns to advect the estuarine plume towards the Rias Baixas, and there was field data available to compare with the model predictions. The obtained results show that the adopted nesting methodology was successful implemented. Model predictions reproduce accurately the hydrodynamics and thermohaline patterns on the Minho estuary and Rias Baixas. The importance of the Minho river discharge and the wind forcing in the event of May 1998 was also studied. The model results showed that a continuous moderate Minho River discharge combined with southerly winds is enough to reverse the Rias Baixas circulation pattern, reducing the importance of the occurrence of specific events of high runoff values. The conditions in which the estuarine plume Minho affects circulation and hydrography of the Rias Baixas were evaluated. The numerical results revealed that the Minho estuarine plume responds rapidly to wind variations and is also influenced by the bathymetry and morphology of the coastline. Without wind forcing, the plume expands offshore, creating a bulge in front of the river mouth. When the wind blows southwards, the main feature is the offshore extension of the plume. Otherwise, northward wind spreads the river plume towards the Rias Baixas. The plume is confined close to the coast, reaching the Rias Baixas after 1.5 days. However, for Minho River discharges higher than 800 m3 s-1, the Minho estuarine plume reverses the circulation patterns in the Rias Baixas. It was also observed that the wind stress and Minho River discharge are the most important factors influencing the size and shape of the Minho estuarine plume. Under the same conditions, the water exchange between Rias Baixas was analysed following the trajectories particles released close to the Minho River mouth. Over 5 days, under Minho River discharges higher than 2100 m3 s-1 combined with southerly winds of 6 m s-1, an intense water exchange between Rias was observed. However, only 20% of the particles found in Ria de Pontevedra come directly from the Minho River. In summary, the model application developed in this study contributed to the characterization and understanding of the influence of the Minho River on the Rias Baixas circulation and hydrography, highlighting that this methodology can be replicated to other coastal systems.

Optimisation of viscoelastic treatments using genetic algorithms

Viscoelastic treatments are one of the most efficient treatments, as far as passive damping is concerned, particularly in the case of thin and light structures. In this type of treatment, part of the strain energy generated in the viscoelastic material is dissipated to the surroundings, in the form of heat. A layer of viscoelastic material is applied to a structure in an unconstrained or constrained configuration, the latter proving to be the most efficient arrangement. This is due to the fact that the relative movement of both the host and constraining layers cause the viscoelastic material to be subjected to a relatively high strain energy. There are studies, however, that claim that the partial application of the viscoelastic material is just as efficient, in terms of economic costs or any other form of treatment application costs. The application of patches of material in specific and selected areas of the structure, thus minimising the extension of damping material, results in an equally efficient treatment. Since the damping mechanism of a viscoelastic material is based on the dissipation of part of the strain energy, the efficiency of the partial treatment can be correlated to the modal strain energy of the structure. Even though the results obtained with this approach in various studies are considered very satisfactory, an optimisation procedure is deemed necessary. In order to obtain optimum solutions, however, time consuming numerical simulations are required. The optimisation process to use the minimum amount of viscoelastic material is based on an evolutionary geometry re-design and calculation of the modal damping, making this procedure computationally costly. To avert this disadvantage, this study uses adaptive layerwise finite elements and applies Genetic Algorithms in the optimisation process.

O desenho da comunicação como conhecimento

As últimas décadas do séc. XX assistiram a um crescente protagonismo do Design de Informação que, desde então, tem sofrido inúmeras formas e designações, num processo de afirmação e auto descoberta. A proliferação de dados disponíveis deu ao design e em particular a este ramo do desenho para a compreensão, uma visibilidade crescente e a responsabilidade de encontrar, a partir da informação, novos meios para a construção de sentido. Do design à engenharia informática, são várias as disciplinas que convergem hoje nesse desígnio ainda que sob diferentes modelos e ferramentas. Esta convergência promove uma comparação entre modelos, que tendem a ser tanto mais valorizados quanto mais objectivas as representações. De Playfair a Bertin, as representações gráficas dos últimos duzentos anos têm-se situado no âmbito de disciplinas como a Economia, Sociologia ou a Gestão, explorando metáforas funcionais com vista à evidência da tradução numérica. O Design, enquanto mediador cultural e através do Desenho, tende a acrescentar ao mesmo exercício uma dimensão narrativa ou ilustrativa, convocando a própria existência do autor na interpretação dos mesmos dados numéricos. Com esta investigação, novos processos de semiose se oferecem, associando à objectividade dos dados quantitativos, a subjectividade da cultura formulada a partir do indivíduo enquanto intérprete. Na procura do conhecimento, reconhece-se assim que o desenho da informação ganha competências pela mediação da experiência, religando ética, técnica e estética.

Ensino superior como princípio hegemónico de cidadania na sociedade civil: um diálogo com António Gramsci no sul da Amazónia brasileira

A tese tem como foco quatro eixos centrais: o ensino superior, a sociedade civil, a cidadania e a hegemonia. Na primeira parte da tese, estes eixos teórico-conceptuais são explanados numa perspectiva da promoção de uma compreensão mais ampla da sociedade civil e do estado, nomeadamente o contratualismo, o liberalismo, o materialismo e o neoliberalismo. Um protagonismo fundamental é atribuído às concepções de Antonio Gramsci de sociedade civil, cidadania e estado aplicadas no ensino superior no Sul da Amazónia brasileira. A segunda parte da tese concentra-se num estudo de caso com três premissas de análise: a fragmentação do ensino superior brasileiro; a reconfiguração da cidadania e hegemonia; e o ensino superior no contexto do Amazonas. O ensino superior no Brasil teve início com as escolas jesuíticas, que, depois de encerradas pelo Marquês de Pombal não tiveram sucessoras em solo brasileiro, ao nível do que poderia ter sido o embrião de um sistema de ensino superior. A chegada da Corte Imperial Portuguesa, em 1808, permite reinstalar novos cursos. A primeira unidade de ensino superior no Amazonas surge apenas no princípio do século XX. Um século depois, tem início o processo de expansão com a implantação de unidades no interior do estado. O ensino superior no Brasil, nas últimas décadas, assumiu um caráter híbrido e de massificação. A massificação, no entanto, ainda é uma realidade a ser alcançada. Apesar dos avanços realizados na última década, ainda subsiste ainda uma forte exclusão de estudantes. O enfoque sobre o ensino superior e a cidadania, no contexto do Amazonas, surge na articulação de um conjunto de informação empírica, extraída de entrevistas realizadas com atores chave locais, com as categorias de pensamento de Antonio Gramsci, que sustentam teoricamente o estudo. Esta articulação tem no ensino superior a possibilidade de desenvolver a cidadania como o princípio organizador e fim último deste nível de ensino. Neste sentido, a hegemonia ganha um estatuto de orientação e direção que permite aos sujeitos envolvidos no ensino superior maior possibilidade de autonomia, liberdade, justiça social, empregabilidade e desenvolvimento social. O modelo de universidade para a cidadania apresenta-se como uma possibilidade de mudanças no horizonte social, económico e também político, no interior da própria universidade. O ensino superior, perspectivado como um instrumento essencial para a cidadania, tem como objetivo primordial a qualificação de professores para a educação de base, Reflexivamente, esta qualificação não deixará potencialmente de produzir retornos positivos na própria expansão e abrangência numérica e educativa do próprio ensino superior. Desta forma, a cidadania no ensino superior é perspectivada neste estudo como um deslocamento de lógicas de compreensão individual e de individualização elitista dos benefícios para lógicas assentes em construções mais coletivas, portadoras de benefícios sociais. A interiorização do ensino superior no Amazonas é assumida no estudo como uma possibilidade de formação académica para o desenvolvimento de práticas pedagógicas e científicas críticas e mais conscientes, constituindo um horizonte determinante para a ativação de processos de integração regional e nacional. Em suma, os pontos de confluência entre o ensino superior, a sociedade civil, a cidadania e as propostas de Gramsci, estão relacionados com a formação, a conscientização política, e o bem-estar económico e social.

Wind energy resource modelling in Portugal and its future large-scale alteration due to anthropogenic induced climate changes

The high dependence of Portugal from foreign energy sources (mainly fossil fuels), together with the international commitments assumed by Portugal and the national strategy in terms of energy policy, as well as resources sustainability and climate change issues, inevitably force Portugal to invest in its energetic self-sufficiency. The 20/20/20 Strategy defined by the European Union defines that in 2020 60% of the total electricity consumption must come from renewable energy sources. Wind energy is currently a major source of electricity generation in Portugal, producing about 23% of the national total electricity consumption in 2013. The National Energy Strategy 2020 (ENE2020), which aims to ensure the national compliance of the European Strategy 20/20/20, states that about half of this 60% target will be provided by wind energy. This work aims to implement and optimise a numerical weather prediction model in the simulation and modelling of the wind energy resource in Portugal, both in offshore and onshore areas. The numerical model optimisation consisted in the determination of which initial and boundary conditions and planetary boundary layer physical parameterizations options provide wind power flux (or energy density), wind speed and direction simulations closest to in situ measured wind data. Specifically for offshore areas, it is also intended to evaluate if the numerical model, once optimised, is able to produce power flux, wind speed and direction simulations more consistent with in situ measured data than wind measurements collected by satellites. This work also aims to study and analyse possible impacts that anthropogenic climate changes may have on the future wind energetic resource in Europe. The results show that the ECMWF reanalysis ERA-Interim are those that, among all the forcing databases currently available to drive numerical weather prediction models, allow wind power flux, wind speed and direction simulations more consistent with in situ wind measurements. It was also found that the Pleim-Xiu and ACM2 planetary boundary layer parameterizations are the ones that showed the best performance in terms of wind power flux, wind speed and direction simulations. This model optimisation allowed a significant reduction of the wind power flux, wind speed and direction simulations errors and, specifically for offshore areas, wind power flux, wind speed and direction simulations more consistent with in situ wind measurements than data obtained from satellites, which is a very valuable and interesting achievement. This work also revealed that future anthropogenic climate changes can negatively impact future European wind energy resource, due to tendencies towards a reduction in future wind speeds especially by the end of the current century and under stronger radiative forcing conditions.

The impact of biofuels for road traffic on air quality : a modelling approach

The selection of the energy source to power the transport sector is one of the main current concerns, not only relative with the energy paradigm but also due to the strong influence of road traffic in urban areas, which highly affects human exposure to air pollutants and human health and quality of life. Due to current important technical limitations of advanced energy sources for transportation purposes, biofuels are seen as an alternative way to power the world’s motor vehicles in a near-future, helping to reduce GHG emissions while at the same time stimulating rural development. Motivated by European strategies, Portugal, has been betting on biofuels to meet the Directive 2009/28/CE goals for road transports using biofuels, especially biodiesel, even though, there is unawareness regarding its impacts on air quality. In this sense, this work intends to clarify this issue by trying to answer the following question: can biodiesel use contribute to a better air quality over Portugal, particularly over urban areas? The first step of this work consisted on the characterization of the national biodiesel supply chain, which allows verifying that the biodiesel chain has problems of sustainability as it depends on raw materials importation, therefore not contributing to reduce the external energy dependence. Next, atmospheric pollutant emissions and air quality impacts associated to the biodiesel use on road transports were assessed, over Portugal and in particular over the Porto urban area, making use of the WRF-EURAD mesoscale numerical modelling system. For that, two emission scenarios were defined: a reference situation without biodiesel use and a scenario reflecting the use of a B20 fuel. Through the comparison of both scenarios, it was verified that the use of B20 fuels helps in controlling air pollution, promoting reductions on PM10, PM2.5, CO and total NMVOC concentrations. It was also verified that NO2 concentrations decrease over the mainland Portugal, but increase in the Porto urban area, as well as formaldehyde, acetaldehyde and acrolein emissions in the both case studies. However, the use of pure diesel is more injurious for human health due to its dominant VOC which have higher chronic hazard quotients and hazard indices when compared to B20.

Characterisation of the cyclic behaviour of reinforced concrete elements with plain bars

Sismos recentes comprovam a elevada vulnerabilidade dos edifícios existentes de betão armado. A resposta das estruturas aos sismos é fortemente condicionada pelas características da aderência aço-betão, que exibe degradação das propriedades iniciais quando sujeitas a carregamentos cíclicos e alternados. Este fenómeno é ainda mais gravoso para elementos com armadura lisa, predominantes na maioria das estruturas construídas até à década de 70 nos países do sul da Europa. A prática corrente de conceção, dimensionamento e pormenorização das estruturas antigas leva a que tenham características de comportamento e níveis de segurança associados não compatíveis com as exigências atuais. Os estudos realizados sobre o comportamento cíclico de elementos estruturais de betão armado com armadura lisa são ainda insuficientes para a completa caracterização deste tipo de elementos. Esta tese visou a caraterização da relação tensão de aderência versus escorregamento para elementos estruturais com armadura lisa e o estudo da resposta cíclica de pilares e nós viga-pilar de betão armado com armadura lisa. Foram realizados dez séries de ensaios de arrancamento (nove monotónicos e um cíclico) em provetes com varões lisos. Os resultados destes ensaios permitiram propor novas expressões empíricas para a estimativa dos parâmetros usados num modelo disponível na literatura para representação da relação tensão de aderência versus escorregamento. É ainda proposto um novo modelo monotónico para a relação tensão de aderência versus escorregamento que representa melhor a resposta após a resistência máxima de aderência. Uma campanha de ensaios unidirecionais em pilares e nós viga-pilar foi também realizada com o objetivo principal de caracterizar o comportamento cíclico deste tipo de elementos. No total foram realizados oito ensaios em pilares, sete ensaios em nós viga-pilar interiores e seis ensaios em nós viga-pilar exteriores representativos de estruturas antigas de betão armado com armadura lisa. Os resultados experimentais permitiram avaliar a influência do escorregamento e estudar o mecanismo de corte em nós e a evolução dos danos para elementos com armadura lisa. Com base nos resultados experimentais foi proposta uma adaptação na expressão do Eurocódigo 8-3 para o cálculo da capacidade última de rotação de elementos com armadura lisa. Foi também desenvolvido um estudo paramétrico, com diferentes estratégias de modelação não linear, para a simulação da resposta de pilares considerando o escorregamento da armadura lisa. Por último, foi proposto um novo modelo simplificado trilinear para o aço que contempla o efeito do escorregamento da armadura lisa.

Modelação lagrangeana dos processos de transporte na Ria Formosa

Apresenta-se um modelo de advecção-difusão do tipo Lagrangeano para a simulação da dispersão de massas de água provenientes de fontes pontuais e descreve-se o modelo hidrodinâmico usado na simulação do escoamento. Estes modelos são utilizados para a caracterização hidrodinâmica da Ria Formosa do ponto de vista da qualidade da água. Tentam-se encontrar razões para os problemas associados aos viveiros de bivalves, sugerindo-se formas de melhorar a qualidade da água.

Aplicação de técnicas de optimização na utilização de acelerogramas registados

Actualmente, existe uma grande difusão de aparelhos digitais do tipo “strong-motion”, por todo mundo, que permitem registar movimentos sísmicos intensos. Também começa a ser frequente esses registos estarem disponíveis, na Internet, para acesso livre. O Eurocódigo 8 (EN 1998-1:2004) permite a representação da acção sísmica, no domínio do tempo, por intermédio de um número, não inferior a três, de acelerogramas registados. No entanto, a escolha desses acelerogramas, de forma a cumprir as regras impostas pelo EC8, reveste-se de alguma complexidade. Neste trabalho, apresentamos uma metodologia para selecção dos acelerogramas a utilizar nas análises sísmicas de estruturas. O que propomos é a escolha de um conjunto de factores multiplicadores dos acelerogramas, que ajustam os valores dos respectivos espectros de resposta às regras impostas pelo EC8. Recorre-se a técnicas de optimização numérica, por forma a quantificar os valores desses factores. Eles são determinados de modo a minimizarem a soma dos desvios em relação ao espectro de resposta objectivo, para o tipo de terreno em causa. A mesma técnica pode ser utilizada na determinação dos valores de TB, TC e TD dos espectros do EC8, depois de normalizados os espectros de resposta dos acelerogramas registados. São apresentados exemplos de aplicação das metodologias propostas a alguns casos de estudo.