Caracterização experimental e modelação numérica de espumas metálicas

As espumas de alumínio são materiais ultraleves, o que as torna atractivas para um largo espectro de aplicações comerciais na área da defesa, na indústria automóvel e aeroespacial, entre outras. Actualmente, há um grande interesse na utilização de espumas de alumínio em componentes estruturais ultraleves, incorporados em sistemas de absorção de energia para protecção contra o impacto. O recurso à simulação numérica para resolver problemas de engenharia em várias áreas é cada vez mais comum. A modelação numérica dos materiais assume vital importância quando o problema envolve a análise de processos tecnológicos como, por exemplo, a conformação plástica de materiais, ou a análise de estruturas. Deste modo, torna-se imprescindível garantir que a modelação dos materiais é de tal forma rigorosa que permite simular o melhor possível o seu comportamento real nas condições concretas da análise a realizar. A forma mais comum de garantir o rigor dos modelos utilizados é a validação dos modelos numéricos tendo por base resultados experimentais. Neste trabalho, fez-se a caracterização do comportamento mecânico das espumas de alumínio com nome comercial ALPORAS!, obtidas pelo processo de fabrico denominado expansão directa do metal fundido por adição de um agente expansor. Esta caracterização consistiu num conjunto de ensaios experimentais quer no regime quasi-estático, quer no regime dinâmico. No regime quasi-estático realizaram-se ensaios de compressão uniaxial e de compressão multiaxial. Para a caracterização no regime dinâmico foram realizados ensaios em barras de Hopkinson de polimetil-metacrilato (PMMA). Com base nos resultados experimentais obtidos determinaram-se os parâmetros dos dois modelos constitutivos para espumas metálicas implementados no programa comercial Abaqus™/Explicit. Estes modelos, e os respectivos parâmetros determinados, foram validados reproduzindo numericamente alguns ensaios experimentais quasi-estáticos e dinâmicos. Assim, verificou-se a adequabilidade dos modelos em diversas condições quer em termos de esforços quer em termos de regime de taxa de deformação Por último, desenvolveu-se uma estrutura inovadora para absorção de energia durante um impacto, constituída por componentes perfilados em liga de alumínio e por componentes em espumas de alumínio. Esta estrutura foi testada exclusivamente com recurso à simulação numérica, utilizando os modelos constitutivos validados anteriormente.

Análise experimental e numérica do comportamento de estruturas de solo reforçado com geossintéticos : solo granular versus fino

As estruturas de solo reforçado com geossintéticos são normalmente constituídas por solos granulares com boas propriedades físicas e mecânicas. O uso de apenas este tipo de solos pode proporcionar o aumento, por vezes insustentável, do custo da execução das estruturas e o aumento do seu impacto ambiental. Deste modo, as estruturas de solo reforçado perdem a sua vantagem competitiva em relação a outros tipos de estruturas (muros de betão, muros de gravidade, muros de gabiões, etc.). Para resolver este problema podem ser utilizados outros tipos de solos (solos locais, finos, com propriedades físicas e mecânicas piores mas, no entanto, mais baratos) para a execução deste tipo de estruturas. De forma geral, com este estudo pretendeu-se contribuir para o incremento do conhecimento sobre a utilização de solos finos para a construção de estruturas de solo reforçado (muros e taludes). Para tal avaliaram-se as diferenças no comportamento mecânico dos materiais compósitos (solo granular reforçado versus solo fino reforçado) e das estruturas de solo reforçado constituídas com os dois tipos de solos. Assim, os objetivos deste estudo foram avaliar: a influência de vários parâmetros nas propriedades mecânicas e na capacidade de carga dos solos reforçados com geossintéticos; a influência de vários parâmetros no dimensionamento das estruturas de solo reforçado; e o comportamento das estruturas dimensionadas (incluindo a estabilidade global e a influência do processo construtivo) recorrendo a uma ferramenta numérica (PLAXIS). Para cumprir os objetivos propostos foram realizadas análises experimentais em laboratório (análise do comportamento do solo reforçado através de ensaios triaxiais e de California Bearing Ratio) e análises numéricas (dimensionamento de estruturas de solo reforçado; modelação numérica do comportamento através de uma ferramenta numérica comercial com o método dos elementos finitos). Os resultados dos ensaios experimentais mostraram que o comportamento mecânico e a capacidade de carga do solo foram incrementados com a inclusão das camadas de geossintético. Este efeito variou com os diversos parâmetros analisados mas, de forma geral, foi mais importante no solo fino (solo com propriedades mecânicas piores). As análises numéricas mostraram que as estruturas de solo fino precisaram de maior densidade de reforços para serem estáveis. Além disso, as estruturas de solo fino foram mais deformáveis e o efeito do seu processo construtivo foi mais importante (principalmente para estruturas de solo fino saturado).

Estudo de sistemas de protecção térmica húmidos sujeitos a condições de incêndio

O presente trabalho tem como principal objectivo o estudo e a interpretação dos mecanismos de transferência de calor e massa no interior de um meio poroso húmido actuando como barreira térmica de protecção em condições de incêndio, visando, sobretudo, melhorar o desempenho destes sistemas. A barreira térmica estudada é constituída por camadas de fibra de vidro, cuja face exposta à chama é protegida por uma folha de alumínio perfurada. A realização deste estudo envolveu a concepção, realização e operação de uma montagem experimental, com a qual se obteve a evolução temporal da temperatura e da taxa de evaporação de água no interior de uma amostra da barreira térmica, submetida à acção directa de uma chama. A análise do processo de protecção foi complementada com resultados obtidos por simulação numérica do processo, designadamente a evolução das propriedades térmicas, que determinam os mecanismos de transferência de calor e massa no interior da barreira. Para esse efeito foi desenvolvido um modelo matemático de simulação do processo de protecção para as barreiras em estudo, que permitiu obter resultados que acompanham bem as evoluções verificadas nos ensaios experimentais. A barreira térmica estudada mostrou claramente possuir um bom desempenho no processo de protecção, face às barreiras térmicas secas. Com a realização deste estudo identificaram-se e quantificaram-se os mecanismos de transferência de calor e massa dominantes nas diversas fases do processo de protecção, o que permitiu extrair conclusões valiosas quanto às características dos meios porosos que contribuem para o seu bom desempenho como barreiras térmicas húmidas.

Multi-scale dynamics on the Eastern Brazilian margin

Este trabalho combina esforços de simulação numérica e de análise de dados para investigar a dinâmica em diversos compartimentos (oceano aberto, plataforma continental e zona costeira-estuarina) e, em multiplas escalas, na Margem Continental Leste Brasileira (MCLB). A circulação de largo e mesoescala espacial e a propagação da maré barotrópica são investigadas através de uma configuração aninhada do modelo numérico ROMS. O estudo da dinâmica regional da Baía de Camamu (CMB) baseia-se na análise de dados locais. A MCLB, localizada a SW do Atlântico Sul entre 8±S e 20±S, possui plataforma estreita, batimetria complexa, e baixa produtividade primária. A sua dinâmica é influenciada pela divergência da Corrente Sul Equatorial (CSE). As simulações refletem as conexões sazonais e espaciais entre a Corrente do Brasil e a Contra Corrente Norte do Brasil , em conexão com a dinâmica da CSE. As simulações revelam atividades vorticais nas proximidades da costa e interações com a dinâmica costeira, cujos padrões são descritos. A validação do modelo em mesoescala é baseada em cálculos de energia cinética turbulenta e em dados históricos de transporte. A CMB, localizada a 13±400S, abriga uma comunidade piscatória tradicional e extenso de manguezal. Situa-se porém sobre uma bacia sedimentar com grande reservas de óleo e gás, estando em tensão permanente de impacto ambiental. Neste trabalho sumarizamos as condições físicas regionais e investigamos sua dinâmica interna, focando sua variabilidade em amostragens realizadas sob condições de seca (Setembro de 2004) e de chuva (Julho de 2005). Finalmente, o modelo numérico ROMS é forçado com o sinal de maré, empregando-se uma configuração simples (com coeficientes de atrito de fundo constantes e condições hidrográficas homogéneas), com o intuito de avaliar sua resposta e investigar a natureza da propagação da maré barotrópica na MCLB, convergindo na CMB. A análise da resposta do modelo à maré basea-se em séries históricas do nível do mar para a MCLB e dados recentes da CMB.

Metrologia óptica com fibra polimérica

O presente trabalho teve por objectivo global o estudo e desenvolvimento de sensores baseados em fibra óptica polimérica. O crescimento da tecnologia polimérica nos últimos anos permitiu a introdução deste tipo de fibras ópticas na área das telecomunicações e no desenvolvimento de sensores. As vantagens associadas à metrologia óptica com fibra polimérica têm vindo a atrair as atenções da comunidade científica dado que permitem o desenvolvimento de sistemas de baixo-custo ou custo competitivo face às tecnologias convencionais. Dada a actualidade do tema proposto, descreve-se, numa primeira fase, a tecnologia em fibra óptica polimérica existente no mercado e o estado de arte de sensores em fibra óptica polimérica. Segue-se a descrição de dois tipos de sensores baseados em modulação de intensidade. Projectou-se um sensor extrínseco capaz de avaliar a quantidade de luz dispersa e absorvida por partículas suspensas num líquido. Foi efectuada a caracterização do sensor quanto à concentração de partículas suspensas, tamanho e reflectividade. O sensor foi testado no âmbito da monitorização ambiental, designadamente, na análise de turbidez em amostras de sedimentos recolhidos em áreas ardidas. O sistema desenvolvido foi comparado com um sistema comercial. Um sensor intrínseco, baseado no polimento lateral de fibra óptica polimérica, foi analisado analiticamente. O modelo teórico avalia o sensor em diferentes condições de macroencurvamento e de índice de refracção do meio envolvente. O modelo teórico foi validado positivamente através de resultados experimentais. Foi avaliada a sensibilidade à temperatura e os conhecimentos adquiridos foram aplicados no desenvolvimento de um sistema capaz de monitorizar a cura de diferentes materiais. É ainda apresentada uma técnica para melhorar a sensibilidade do sensor de curvatura através da aplicação de um revestimento na zona sensível. A dependência na curvatura da potência transmitida por uma fibra óptica polida lateralmente serviu de base ao desenvolvimento de uma joelheira e de uma cotoveleira instrumentada, capazes de avaliar quantitativamente o movimento articular. A necessidade de portabilidade levou ao desenvolvimento de um sistema sem fios para aquisição e transmissão de dados. Espera-se que os protótipos desenvolvidos venham a ter um impacto significativo em sistemas futuros aplicados à medicina física e reabilitação.

Hydraulics of stepped chutes: experimental-numerical-theoretical study

O estudo de escoamentos turbulentos em descarregadores em degraus tem sido um desafio para os investigadores. A macro-rugosidade do leito, a ondulação da superfície livre, a intermitência da localização da secção inicial de entrada de ar e o escoamento bi-fásico a jusante da secção inicial de entrada de ar fazem com que a caracterização do escoamento deslizante sobre turbilhões em descarregadores em degraus não seja simples. Actualmente, é possível combinar técnicas de medição fiáveis com simulações numéricas e análise teórica. Nesta dissertação, o estudo experimental baseia-se em resultados experimentais obtidos em duas instalações experimentais: a instalação A, do Laboratório Nacional de Engenharia Civil (LNEC), dotada de um descarregador em degraus com declive de 1V:0.75H e a instalação B, do Instituto Superior Técnico (IST), dotada de um descarregador em degraus com declive de 1V:2H. Uma sonda de ar, um tubo de Pitot modificado e vários hidrómetros permitiram o estudo do escoamento deslizante sobre turbilhões mono- e bifásico ao longo do descarregador em degraus da instalação A. Na bacia de dissipação de energia, a observação visual do escoamento foi auxiliada pelas leituras efectuadas em tomadas de pressão localizadas na soleira da bacia e pelas leituras da altura do escoamento efectuadas em réguas graduadas localizadas nas paredes da bacia. Na instalação B foram utilizados um tubo de Pitot e vários hidrómetros para estudar a região não arejada do escoamento deslizante sobre turbilhões no descarregador em degraus. O estudo numérico é baseado em simulações numéricas da região não arejada do escoamento deslizante sobre turbilhões sobre descarregadores em degraus com o código comercial de CFD FLOW-3D® de modo a reproduzir as condições ensaiadas experimentalmente. As simulações beneficiaram da técnica de blocos múltiplos (multi-block) num sistema de coordenadas cartesianas, da determinação da superfície livre pelo método TruVOF e da utilização de dois modelos de turbulência: os modelos k- e RNG k-. Por último, o estudo teórico consistiu em desenvolver um modelo simplificado 1D para determinar as características hidráulicas principais do trecho não arejado do escoamento deslizante sobre turbilhões em descarregadores em degraus. O modelo foi desenvolvido a partir das equações de Navier-Stokes, conjuntamente com resultados experimentais e numéricos. Os resultados apresentados nesta dissertação contribuem para o conhecimento do escoamento deslizante sobre turbilhões em descarregadores em degraus, nomeadamente na região não arejada, na secção inicial de entrada de ar e na região arejada. A hidráulica dos dissipadores de energia a jusante de descarregadores em degraus, em particular das bacias tipo III do USBR, é também objecto de estudo. Em relação à região não arejada do escoamento em descarregadores com declive acentuado, e com base em resultados experimentais e numéricos, são propostas expressões para estimar o desenvolvimento da altura equivalente de água, da espessura da camada limite, da concentração média de ar, do coeficiente de energia cinética, da dissipação de energia, do factor de resistência e do coeficiente n da fórmula de Manning. São ainda propostas expressões adimensionais para a energia cinética turbulenta e sua dissipação. Para declives moderados, são propostas expressões para estimar o desenvolvimento da altura equivalente de água, do coeficiente de energia cinética e da energia específica residual. São ainda apresentados valores do expoente 1/N da expressão adimensional da distribuição de velocidades, quer para descarregadores com declive acentuado quer com declive moderado. Em conformidade com outros estudos centrados em escoamentos de parede e com derivações teóricas, para a região não arejada do escoamento em descarregadores em degraus de acentuado declive, observa-se que o factor de resistência depende da macro-rugosidade criada pelos degraus e da geometria da secção transversal e que o coeficiente n da fórmula de Manning aumenta com a rugosidade. A descrição estatística da turbulência do escoamento é igualmente explorada, contribuindo para o conhecimento da estrutura do escoamento. Observou-se que para números de Reynolds rugoso não superiores a 6.8x104 a energia cinética turbulência e a sua dissipação cumprem leis de semelhança. Estas expressões adimensionais estão de acordo com os resultados obtidos por outros autores para escoamentos completamente desenvolvidos em canais abertos e no escoamento em rios com leito de gravilha. Em acréscimo, a taxa de dissipação de energia, quer para descarregadores de declive acentuado quer de moderado declive, é baixa. Por último, observa-se que os valores da média temporal da concentração de ar entre 0 e 1 medidos na região não arejada do escoamento dizem respeito não só ao ar capturado entre ondas de água, na zona de ondulação da superfície livre, mas também ao ar emulsionado no escoamento, i.e., sob a forma de bolhas de ar, quando perto da secção média inicial de entrada de ar, devido à diferença entre localizações instantânea e média temporal. Foram revistas metodologias e fórmulas para estimar a localização da secção inicial de entrada de ar e apresentadas expressões para estimar a concentração média de ar e a altura equivalente de água nessa secção. Relativamente à região de escoamento arejado em descarregadores em degraus com declive acentuado, os resultados experimentais apresentados nesta dissertação permitiram estimar a influência da definição da superfície livre nos parâmetros hidráulicos da região do escoamento arejado e estimar a máxima elevação do escoamento nesta região do escoamento. Com base nos resultados experimentais obtidos na bacia de dissipação de energia do tipo III do USBR localizada a jusante do descarregador em degraus da instalação A, observou-se que os perfis da altura piezométrica e da altura do escoamento tendem a seguir o perfil recomendado pelo USBR para bacias tipo III. A excepção ocorre à entrada da bacia, onde as alturas piezométricas apresentadas nesta dissertação excedem largamente as apresentadas pelo USBR. É ainda observado que, tal como entre as bacias tipo I e tipo III do USBR, o ressalto hidráulico estabiliza muito mais rapidamente numa bacia tipo III a jusante de um descarregador em degraus do que uma bacia tipo I a jusante do mesmo descarregador em degraus. Finalmente, observa-se que os blocos de amortecimento a colocar no descarregador não têm influência visível nos resultados da altura piezométrica nem da altura do escoamento ao longo da bacia. Relativamente às simulações numéricas do escoamento não arejado, a proximidade entre resultados experimentais e numéricos permite validar o modelo teórico e a integração numérica usados no FLOW-3D®. As simulações desenvolvidas também mostraram que o modelo de turbulência k- permite representar as características do escoamento não arejado em descarregadores em degraus, uma vez que não foram observadas diferenças significativas entre as simulações com este modelo e com o modelo RNG k-. Finalmente, observou-se que o modelo de entrada de ar usado no FLOW-3D® é válido para estimar a localização da secção inicial de entrada de ar. Por último, a proximidade entre os resultados obtidos da aplicação do modelo teórico desenvolvido no âmbito desta dissertação e os resultados experimentais indica que as hipóteses e simplificações consideradas no desenvolvimento do modelo são adequadas.

A study of singular integral operators with shift

Nesta tese, consideram-se operadores integrais singulares com a acção extra de um operador de deslocacamento de Carleman e com coeficientes em diferentes classes de funções essencialmente limitadas. Nomeadamente, funções contínuas por troços, funções quase-periódicas e funções possuíndo factorização generalizada. Nos casos dos operadores integrais singulares com deslocamento dado pelo operador de reflexão ou pelo operador de salto no círculo unitário complexo, obtêm-se critérios para a propriedade de Fredholm. Para os coeficientes contínuos, uma fórmula do índice de Fredholm é apresentada. Estes resultados são consequência das relações de equivalência explícitas entre aqueles operadores e alguns operadores adicionais, tais como o operador integral singular, operadores de Toeplitz e operadores de Toeplitz mais Hankel. Além disso, as relações de equivalência permitem-nos obter um critério de invertibilidade e fórmulas para os inversos laterais dos operadores iniciais com coeficientes factorizáveis. Adicionalmente, aplicamos técnicas de análise numérica, tais como métodos de colocação de polinómios, para o estudo da dimensão do núcleo dos dois tipos de operadores integrais singulares com coeficientes contínuos por troços. Esta abordagem permite também a computação do inverso no sentido Moore-Penrose dos operadores principais. Para operadores integrais singulares com operadores de deslocamento do tipo Carleman preservando a orientação e com funções contínuas como coeficientes, são obtidos limites superiores da dimensão do núcleo. Tal é implementado utilizando algumas estimativas e com a ajuda de relações (explícitas) de equivalência entre operadores. Focamos ainda a nossa atenção na resolução e nas soluções de uma classe de equações integrais singulares com deslocamento que não pode ser reduzida a um problema de valor de fronteira binomial. De forma a atingir os objectivos propostos, foram utilizadas projecções complementares e identidades entre operadores. Desta forma, as equações em estudo são associadas a sistemas de equações integrais singulares. Estes sistemas são depois analisados utilizando um problema de valor de fronteira de Riemann. Este procedimento tem como consequência a construção das soluções das equações iniciais a partir das soluções de problemas de valor de fronteira de Riemann. Motivados por uma grande diversidade de aplicações, estendemos a definição de operador integral de Cauchy para espaços de Lebesgue sobre grupos topológicos. Assim, são investigadas as condições de invertibilidade dos operadores integrais neste contexto.

Contribution to the study of the Ria de Aveiro inlet morphodynamics

Over the years it was observed at the Ria de Aveiro lagoon inlet, near the head of the north breakwater, a depth increase that might threaten the stability of this structure. A trend of accretion in the navigation channel of this lagoon is observed, endangering the navigation in this region. In order to understand the origin of these and other trends observed, the knowledge of the sediment transport in the study area is imperative. The main aim of this work is understanding the dominant physical processes in the sediment transport of sediment at the Ria de Aveiro lagoon inlet and adjacent area, improving knowledge of this region morphodynamics. The methodology followed in this study consisted in the analyzes of the topohydrographic surveys performed by the Administration of the Aveiro Harbor, and in the numerical simulations results performed with the morphodynamic modeling system MORSYS2D. The analysis of the surveys was performed by studying the temporal evolution of the bathymetry. The numerical analysis was based on the implementation of the model at the study area, sensitivity analysis of the formulations used to compute the sediment transport to the variation of input parameters (e.g. depth, sediment size, tidal currents) and analysis of the sediment uxes and bathymetric changes predicted. The simulations considered as sediment transport forcing the tidal currents only and the coupled forcing of tides and waves. Considering the wave e ect as sediment transport forcing, both monochromatic waves and a wave regime were simulated. The results revealed that the observed residual sediment transport patterns are generated due to the channel con guration. Inside the lagoon the uxes are mainly induced by the tidal currents action, restricting the action of waves to the inlet and adjacent coast. In the navigation channel the residual sediment uxes predicted are directed o - shore with values between 7 and 40 m3=day generating accretions of approximately 10 m3=day for the shallower region and 35 m3=day for the region between the tidal gauge and the tri^angulo das mar es. At the inlet, the residual uxes are approximately 30 m3=day inducing trends of erosion of approximately 20 m3=day. At the North side of the nearshore accretion is predicted, while at the South side is predicted erosion, at the rates of 250 and 1500 m3=day, respectively. It was also concluded that the waves with higher contribution to the residual sediment uxes are those with heights between 4 and 5 m. However, the storm waves with heights bigger than 5 m, despite their 10% of frequency of occurrence are responsible for 25% of the observed sediment transport.

Cyclic behaviour of RC elements with plain reinforcing bars

O comportamento cíclico das estruturas de betão armado é fortemente condicionado pelo mecanismo de aderência entre o betão e o aço. O escorregamento relativo entre os dois materiais, resultante da degradação progressiva da aderência em elementos solicitados por ações cíclicas, é uma causa frequente de danos graves e até do colapso de estruturas devido à ocorrência de sismos. Entre as estruturas existentes de betão armado que foram dimensionadas e construídas antes da entrada em vigor dos regulamentos sísmicos atuais, muitas foram construídas com armadura lisa, e portanto, possuem fracas propriedades de aderência. A informação disponível na literatura sobre o comportamento cíclico de elementos estruturais de betão armado com armadura lisa é reduzida e a influência das propriedades da aderência associadas a este tipo de armadura no comportamento cíclico das estruturas existentes não se encontra ainda devidamente estudada. O objectivo principal desta tese foi estudar a influência do escorregamento na resposta cíclica de elementos estruturais de betão armado com armadura lisa. Foram realizados ensaios cíclicos em elementos do tipo nó viga-pilar, construídos à escala real, representativos de ligações interiores em edifícios existentes sem pormenorização específica para resistir às ações sísmicas. Para comparação, foi realizado o ensaio de um nó construído com armadura nervurada. Foi ainda realizado o ensaio cíclico de uma viga de betão armado recolhida de uma estrutura antiga. Foram elaborados modelos numéricos não-lineares para simular a resposta dos elementos ensaiados, concentrando especial atenção no mecanismo do escorregamento. Os resultados obtidos no âmbito desta tese contribuem para o avanço do conhecimento sobre o comportamento cíclico de elementos estruturais de betão armado com armadura lisa. As análises numéricas realizadas comprovam a necessidade de incluir os efeitos do escorregamento na modelação numérica deste tipo de estruturas de forma a representar com rigor a sua resposta às ações cíclicas.

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.