Formação congênita de satélites coorbitais: simulações numéricas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Jato transversal de gotas: simulações por ALE/FEM e efeitos interfaciais.

Um código computacional para escoamentos bifásicos incorporando metodologia híbrida entre oMétodo dos Elementos Finitos e a descrição Lagrangeana-Euleriana Arbitrária do movimento é usado para simular a dinâmica de um jato transversal de gotas na zona primária de quebra. Os corpos dispersos são descritos por meio de um método do tipo front-tracking que produz interfaces de espessura zero através de malhas formadas pela união de elementos adjacentes em ambas as fases e de técnicas de refinamento adaptativo. Condições de contorno periódicas são implementadas de modo variacionalmente consistente para todos os campos envolvidos nas simulações apresentadas e uma versão modificada do campo de pressão é adicionada à formulação do tipo um-fluido usada na equação da quantidade de movimento linear. Simulações numéricas diretas em três dimensões são executadas para diferentes configurações de líquidos imiscí veis compatíveis com resultados experimentais encontrados na literatura. Análises da hidrodinâmica do jato transversal de gotas nessas configurações considerando trajetórias, variação de formato de gota, espectro de pequenas perturbações, além de aspectos complementares relativos à qualidade de malha são apresentados e discutidos.

Identificação de danos estruturais utilizando dados no domínio do tempo provenientes de ensaios de vibração

O presente trabalho aborda o problema de identificação de danos em uma estrutura a partir de sua resposta impulsiva. No modelo adotado, a integridade estrutural é continuamente descrita por um parâmetro de coesão. Sendo assim, o Modelo de Elementos Finitos (MEF) é utilizado para discretizar tanto o campo de deslocamentos, quanto o campo de coesão. O problema de identificação de danos é, então, definido como um problema de otimização, cujo objetivo é minimizar, em relação a um vetor de parâmetros nodais de coesão, um funcional definido a partir da diferença entre a resposta impulsiva experimental e a correspondente resposta prevista por um MEF da estrutura. A identificação de danos estruturais baseadas no domínio do tempo apresenta como vantagens a aplicabilidade em sistemas lineares e/ou com elevados níveis de amortecimento, além de apresentar uma elevada sensibilidade à presença de pequenos danos. Estudos numéricos foram realizados considerando-se um modelo de viga de Euler-Bernoulli simplesmente apoiada. Para a determinação do posicionamento ótimo do sensor de deslocamento e do número de pontos da resposta impulsiva, a serem utilizados no processo de identificação de danos, foi considerado o Projeto Ótimo de Experimentos. A posição do sensor e o número de pontos foram determinados segundo o critério D-ótimo. Outros critérios complementares foram também analisados. Uma análise da sensibilidade foi realizada com o intuito de identificar as regiões da estrutura onde a resposta é mais sensível à presença de um dano em um estágio inicial. Para a resolução do problema inverso de identificação de danos foram considerados os métodos de otimização Evolução Diferencial e Levenberg-Marquardt. Simulações numéricas, considerando-se dados corrompidos com ruído aditivo, foram realizadas com o intuito de avaliar a potencialidade da metodologia de identificação de danos, assim como a influência da posição do sensor e do número de dados considerados no processo de identificação. Com os resultados obtidos, percebe-se que o Projeto Ótimo de Experimentos é de fundamental importância para a identificação de danos.

Simulador MATLAB para Sistemas Ópticos WDM Amplificados.

Neste trabalho é apresentado um simulador MATLAB para sistemas ópticos WDM amplificados baseado na solução das equações não lineares de Schrödinger acopladas, pelo método de Fourier de passo alternado. Este simulador permite o estudo da propagação de pulsos em fibras ópticas, considerando dispersão cromática, efeitos não lineares como automodulação de fase e modulação de fase cruzada e atenuação, prevendo também o emprego de amplificadores ópticos a fibra dopada com Érbio (EDFAs). Através de simulações numéricas, foi explorada a técnica de otimização do posicionamento de um EDFA ao longo de um enlace óptico, sem repetidores, que objetiva a redução dos custos de implantação de sistemas ópticos, seja pela diminuição da potência do transmissor ou pela relaxação da exigência de sensibilidade do receptor. Além disto, pode favorecer um aumento na capacidade do sistema, através do aumento do alcance ou da taxa de transmissão. A concordância dos resultados obtidos com os disponíveis na literatura confirmam a validade da técnica, bem como a versatilidade e robustez do simulador desenvolvido.

Simulação numérica da dispersão advecção de pesticidas no solo sob efeito da temperatura.

É apresentado um modelo de dispersão-advecção de evolução unidimensional que simula a lixiviação de pesticidas em lisímetros ou colunas de solo sob efeito da temperatura média diária do perfil do solo. O modelo matemático e todo o seu conjunto de hipóteses será denominado DAPESTE. Nas simulações numéricas do modelo DAPESTE serão utilizados o método dos elementos finitos para a semi-discretização da variável espacial e o método de Eüler regressivo para a discretização da variável temporal. Serão utilizados métodos de elementos finitos apropriados para problemas de dispersão-advecção nos quais o transporte advectivo predomina sobre o transporte dispersivo. O modelo DAPESTE supõe que as difusividades do pesticida nas fases gasosa e aquosa do solo dependem da temperatura média diária do solo a qual varia periodicamente com a profundidade e com o tempo. O coeficiente de dispersão hidrodinâmico do modelo DAPESTE dependerá da temperatura do solo. O coeficiente de partição água-ar, variando com a temperatura, será determinado pela equação de Clausius-Clapeyron. A equação de van?t Ho® será usada para determinar o coeficiente de sorção do pesticida no solo como função da temperatura. Com a equação de Arrhenius será estimado o efeito da temperatura na taxa de degradação do pesticida. Estas relações de dependência entre os parâmetros do modelo e a temperatura do solo auxiliam sobremaneira na compreensão do destino de pesticidas no solo sob diferentes cenários de temperaturas médias diárias, especialmente a meia vida do pesticida e a concentração lixiviada no perfil do solo.

Numerical modelling of radionuclide migration in near-surface systems

Em todo o mundo são usados, hoje em dia, modelos numéricos hidrogeoquímicos para simular fenómenos naturais e fenómenos decorrentes de actividades antrópicas. Estes modelos ajudam-nos a compreender o ambiente envolvente, a sua variabilidade espacial e evolução temporal. No presente trabalho apresenta-se o desenvolvimento de modelos numéricos hidrogeoquímicos aplicados no contexto do repositório geológico profundo para resíduos nucleares de elevada actividade. A avaliação da performance de um repositório geológico profundo inclui o estudo da evolução geoquímica do repositório, bem como a análise dos cenários de mau funcionamento do repositório, e respectivas consequências ambientais. Se se escaparem acidentalmente radionuclídeos de um repositório, estes poderão atravessar as barreiras de engenharia e barreiras naturais que constituem o repositório, atingindo eventualmente, os ecosistemas superficiais. Neste caso, os sedimentos subsuperficiais constituem a última barreira natural antes dos ecosistemas superficiais. No presente trabalho foram desenvolvidos modelos numéricos que integram processos biogeoquímicos, geoquímicos, hidrodinâmicos e de transporte de solutos, para entender e quantificar a influência destes processos na mobilidade de radionuclídeos em sistemas subsuperficiais. Os resultados alcançados reflectem a robustez dos instrumentos numéricos utilizados para desenvolver simulações descritivas e predictivas de processos hidrogeoquímicos que influenciam a mobilidade de radionuclídeos. A simulação (descritiva) de uma experiência laboratorial revela que a actividade microbiana induz a diminuição do potencial redox da água subterrânea que, por sua vez, favorece a retenção de radionuclídeos sensíveis ao potencial redox, como o urânio. As simulações predictivas indicam que processos de co-precipitação com minerais de elementos maioritários, precipitação de fases puras, intercâmbio catiónico e adsorção à superfície de minerais favorecem a retenção de U, Cs, Sr e Ra na fase sólida de uma argila glaciar e uma moreia rica em calcite. A etiquetagem dos radionuclídeos nas simulações numéricas permitiu concluir que a diluição isotópica joga um papel importante no potencial impacte dos radionuclídeos nos sistemas subsuperficiais. A partir dos resultados das simulações numéricas é possivel calcular coeficientes de distribuição efectivos. Esta metodologia proporciona a simulação de ensaios de traçadores de longa duração que não seriam exequíveis à escala da vida humana. A partir destas simulações podem ser obtidos coeficientes de retardamento que são úteis no contexto da avaliação da performance de repositórios geológicos profundos.

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.

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.

Influência de plumas de água doce na hidrodinâmica do estuário da Ria de Arousa (NW Espanha), uma abordagem de modelação matemática

Localizada na parte sul da costa da Galiza, a ria de Arousa, tem um interesse particular devido à sua importante indústria de exploração de marisco, especialmente culturas de bivalves. O crescimento, qualidade e mortalidade destes bivalves não só são afetados por possíveis incidentes de poluição, mas também por episódios de grande run-off, que podem produzir decréscimos sazonais significativos de salinidade. Neste contexto, a presente dissertação pretende contribuir para uma melhor compreensão de como as plumas de água doce provenientes dos rios vão influenciar os gradientes de salinidade e temperatura no estuário da ria de Arousa. O objetivo é proceder a uma melhoria da previsão hidrodinâmica do modelo operacional oceanográfico desta ria, que é executado diariamente pela MeteoGalicia. Para atingir esse objectivo executaram-se simulações numéricas para diferentes cenários de forçamento. Foram escolhidos dois períodos com diferentes condições meteorológicas: um período de inverno, coincidindo com chuvas fortes e um período de verão, coincidindo com a ausência de precipitação. Foi aplicado o modelo hidrodinâmico MOHID (Hydrodynamic Model) para simular os processos físicos da ria de Arousa, utilizando a metodologia de modelos encaixados. O forçamento para a superfície foi imposto através do modelo atmosférico WRF (Weather Regional Forecast), que fornece as variáveis meteorológicas necessárias. Para fronteira terra, uma ferramenta hidrológica foi implementada, o modelo SWAT (Soil water Assessment Tool), este providencia as descargas de água doce para os principais rios da região. É feita neste trabalho uma descrição da implementação deste sistema, assim como a caracterização da área de estudo. Os resultados do modelo foram comparados com medidas provenientes de bóias oceanográficas e perfis de CTDs (Conductivity-Temperature-Depth instrument) em diferentes localizações.

Dynamics and variability of the alongshore flows on the northwestern iberian margin

Tese de doutoramento, Ciências Geofísicas e da Geoinformação (Oceanografia), Universidade de Lisboa, Faculdade de Ciências, 2014

Controlo directo de binário de uma máquina assíncrona trifásica

Dissertação para obtenção do grau de Mestre em Engenharia Electrotécnica Ramo de Automação e Electrónica Industrial

Estudo do escoamento dos produtos da combustão no interior de uma conduta de evacuação de turbinas a gás utilizadas para propulsão naval

Dissertação apresentada na Faculdade de Ciências e Tecnologia da Universidade Nova de Lisboa para obtenção do grau de Mestre em Engenharia Mecânica

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

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

Avaliação do desempenho de ligações tubulares em estruturas de veículos pesados de passageiros por simulação numérica

Nesta Dissertação ir-se-á avaliar o desempenho de algumas ligações existentes nos veículos pesados de passageiros, através do Regulamento Eurocódigo 3. Nos últimos anos ocorreram diversos acidentes envolvendo este tipos de veículos, em que os mesmos causaram vítimas mortais e feridos graves. Serão testadas por simulação numérica algumas ligações pertencentes a elementos constituintes da superestrutura, em que esta é normalmente afectada com a ocorrência de acidentes. Assim sendo, o estudo de nós de ligação tem uma importância fulcral para que uma superestrutura suporte situações extremas e que resista a solicitações externas aplicadas. Iniciou-se esta Dissertação com o estudo da sinistralidade e de acidentes que envolvem veículos pesados de passageiros. No capítulo 2 abordou-se um programa que promove simulações numéricas de acidentes e estudo do comportamento dos passageiros em caso de acidente, sendo referido o Regulamento que homologa os veículos pesados de passageiros e os seus principais métodos. Abordaram-se os principais constituintes da estrutura de um veículo pesado de passageiros. No capítulo 3, é referido o Eurocódigo 3 em termos do estudo das ligações tubulares usadas neste tipo de veículos. No capítulo 4, fez-se o estudo e selecção de elementos a utilizar para a simulação numérica de casos preconizados pelo Eurocódigo 3 e estudaram-se três tipos de ligações que são usadas na construção da superestrutura deste tipo de veículos, tendo-se retirado conclusões deste estudo.