Estimação robusta em modelos lineares de equações simultâneas

Os Modelos de Equações Simultâneas (SEM) são modelos estatísticos com muita tradição em estudos de Econometria, uma vez que permitem representar e estudar uma vasta gama de processos económicos. Os estimadores mais usados em SEM resultam da aplicação do Método dos Mínimos Quadrados ou do Método da Máxima Verosimilhança, os quais não são robustos. Em Maronna e Yohai (1997), os autores propõem formas de “robustificar” esses estimadores. Um outro método de estimação com interesse nestes modelos é o Método dos Momentos Generalizado (GMM), o qual também conduz a estimadores não robustos. Estimadores que sofrem de falta de robustez são muito inconvenientes uma vez que podem conduzir a resultados enganadores quando são violadas as hipóteses subjacentes ao modelo assumido. Os estimadores robustos são de grande valor, em particular quando os modelos em estudo são complexos, como é o caso dos SEM. O principal objectivo desta investigação foi o de procurar tais estimadores tendo-se construído um estimador robusto a que se deu o nome de GMMOGK. Trata-se de uma versão robusta do estimador GMM. Para avaliar o desempenho do novo estimador foi feito um adequado estudo de simulação e foi também feita a aplicação do estimador a um conjunto de dados reais. O estimador robusto tem um bom desempenho nos modelos heterocedásticos considerados e, nessas condições, comporta-se melhor do que os estimadores não robustos usados no estudo. Contudo, quando a análise é feita em cada equação separadamente, a especificidade de cada equação individual e a estrutura de dependência do sistema são dois aspectos que influenciam o desempenho do estimador, tal como acontece com os estimadores usuais. Para enquadrar a investigação, o texto inclui uma revisão de aspectos essenciais dos SEM, o seu papel em Econometria, os principais métodos de estimação, com particular ênfase no GMM, e uma curta introdução à estimação robusta.

Exact and numerical solutions for diffusio-reaction equations with convection term

Esta dissertação estuda essencialmente dois problemas: (A) uma classe de equações unidimensionais de reacção-difusão-convecção em meios não uniformes (dependentes do espaço), e (B) um problema elíptico não-linear e paramétrico ligado a fenómenos de capilaridade. A Análise de Perturbação Singular e a dinâmica de Hamilton-Jacobi são utilizadas na obtenção de expressões assimptóticas para a solução (com comportamento de frente) e para a sua velocidade de propagação. Os seguintes três métodos de decomposição, Adomian Decomposition Method (ADM), Decomposition Method based on Infinite Products (DIP), e New Iterative Method (NIM), são apresentados e brevemente comparados. Adicionalmente, condições suficientes para a convergência da solução em série, obtida pelo ADM, e uma aplicação a um problema da Telecomunicações por Fibras Ópticas, envolvendo EDOs não-lineares designadas equações de Raman, são discutidas. Um ponto de vista mais abrangente que unifica os métodos de decomposição referidos é também apresentado. Para subclasses desta EDP são obtidas soluções numa forma explícita, para diferentes tipos de dados e usando uma variante do método de simetrias de Bluman-Cole. Usando Teoria de Pontos Críticos (o teorema usualmente designado mountain pass) e técnicas de truncatura, prova-se a existência de duas soluções não triviais (uma positiva e uma negativa) para o problema elíptico não-linear e paramétrico (B). A existência de uma terceira solução não trivial é demonstrada usando Grupos Críticos e Teoria de Morse.

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.

Photonic devices for optical and RF signal processing

O presente trabalho tem por objectivo o estudo de novos dispositivos fotónicos aplicados a sistemas de comunicações por fibra óptica e a sistemas de processamento de sinais RF. Os dispositivos apresentados baseiam-se em processamento de sinal linear e não linear. Dispositivos lineares ópticos tais como o interferómetro de Mach-Zehnder permitem adicionar sinais ópticos com pesos fixos ou sintonizáveis. Desta forma, este dispositivo pode ser usado respectivamente como um filtro óptico em amplitude com duas saídas complementares, ou, como um filtro óptico de resposta de fase sintonizável. O primeiro princípio de operação serve como base para um novo sistema fotónico de medição em tempo real da frequência de um sinal RF. O segundo princípio de operação é explorado num novo sistema fotónico de direccionamento do campo eléctrico radiado por um agregado de antenas, e também num novo compensador sintonizável de dispersão cromática. O processamento de sinal é não linear quando sinais ópticos são atrasados e posteriormente misturados entre si, em vez de serem linearmente adicionados. Este princípio de operação está por detrás da mistura de um sinal eléctrico com um sinal óptico, que por sua vez é a base de um novo sistema fotónico de medição em tempo real da frequência de um sinal RF. A mistura de sinais ópticos em meios não lineares permite uma operação eficiente numa grande largura espectral. Tal operação é usada para realizar conversão de comprimento de onda sintonizável. Um sinal óptico com multiplexagem no domínio temporal de elevada largura de banda é misturado com duas bombas ópticas não moduladas com base em processos não lineares paramétricos num guia de ondas de niobato de lítio com inversão periódica da polarização dos domínios ferroeléctricos. Noutro trabalho, uma bomba pulsada em que cada pulso tem um comprimento de onda sintonizável serve como base a um novo conversor de sinal óptico com multiplexagem no domínio temporal para um sinal óptico com multiplexagem no comprimento de onda. A bomba é misturada com o sinal óptico de entrada através de um processo não linear paramétrico numa fibra óptica com parâmetro não linear elevado. Todos os dispositivos fotónicos de processamento de sinal linear ou não linear propostos são experimentalmente validados. São também modelados teoricamente ou através de simulação, com a excepção dos que envolvem mistura de sinais ópticos. Uma análise qualitativa é suficiente nestes últimos dispositivos.

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.

Multiplicity results for some classes of Schrödinger-Poisson systems

In this thesis, we study the existence and multiplicity of solutions of the following class of Schr odinger-Poisson systems: u + u + l(x) u = (x; u) in R3; = l(x)u2 in R3; where l 2 L2(R3) or l 2 L1(R3). And we consider that the nonlinearity satis es the following three kinds of cases: (i) a subcritical exponent with (x; u) = k(x)jujp 2u + h(x)u (4 p < 2 ) under an inde nite case; (ii) a general inde nite nonlinearity with (x; u) = k(x)g(u) + h(x)u; (iii) a critical growth exponent with (x; u) = k(x)juj2 2u + h(x)jujq 2u (2 q < 2 ). It is worth mentioning that the thesis contains three main innovations except overcoming several di culties, which are generated by the systems themselves. First, as an unknown referee said in his report, we are the rst authors concerning the existence of multiple positive solutions for Schr odinger- Poisson systems with an inde nite nonlinearity. Second, we nd an interesting phenomenon in Chapter 2 and Chapter 3 that we do not need the condition R R3 k(x)ep 1dx < 0 with an inde nite noncoercive case, where e1 is the rst eigenfunction of +id in H1(R3) with weight function h. A similar condition has been shown to be a su cient and necessary condition to the existence of positive solutions for semilinear elliptic equations with inde nite nonlinearity for a bounded domain (see e.g. Alama-Tarantello, Calc. Var. PDE 1 (1993), 439{475), or to be a su cient condition to the existence of positive solutions for semilinear elliptic equations with inde nite nonlinearity in RN (see e.g. Costa-Tehrani, Calc. Var. PDE 13 (2001), 159{189). Moreover, the process used in this case can be applied to study other aspects of the Schr odinger-Poisson systems and it gives a way to study the Kirchho system and quasilinear Schr odinger system. Finally, to get sign changing solutions in Chapter 5, we follow the spirit of Hirano-Shioji, Proc. Roy. Soc. Edinburgh Sect. A 137 (2007), 333, but the procedure is simpler than that they have proposed in their paper.

A influência do terrorismo nos comportamentos de viagem internacionais

A partir de uma amostra de 600 turistas internacionais que circulam em Portugal, Espanha e Itália, este estudo identifica as principais os conceitos chave relacionados com o terrorismo, a percepção de risco, envolvimento e motivação para a segurança dos turistas internacionais. Diferentes níveis de preocupação relativamente à segurança pode influenciar as decisões dos turistas. No seu processo de decisão, os turistas avaliam vários factores, nomeadamente, o nível de risco ou de segurança que consideram nos destinos (Sonmez, 1998). Os turistas adoptam uma atitude protectora alterando os seus comportamentos durante os processos de decisão, substituindo os destinos que consideram inseguros por outros associados a uma maior segurança (Gu & Martin, 1992; Mansfeld, 1996). O terrorismo exacerbado pelos media tem efeitos graves nas receitas dos destinos turísticos (Taylor, 2006). Através da publicidade negativa, um destino turístico que experiencia um incidente terrorista pode ver a sua reputação danificada e a actividade turística severamente comprometida (Sonmez, 1998). Inclusivamente, a imageme negativa de um destino pode ser generalizada e pode também afectar outros países ou regiões por períodos de tempo indeterminados (Taylor, 2006). Um modelo de equações estruturais revela que os turistas são motivados para adquirir informação sobre o terrorismo nos media, nomeadamente mostram atenção e interesse sobre essas notícias e esse facto influencia directamente o seu risco percebido. A percepção de risco influencia directamente o envolvimento dos turistas no planeamento da viagem, especificamente a procura de informação antes da viagem e quando estão no destino. A percepção de risco e o envolvimento dos turistas influencia a percepção da importância da segurança.A discussão foca as implicações deste modelo para a teoria e para as instituições e organizações turísticas. São igualmente apresentadas recomendações para os gestores e promotores dos destinos e para os gestores das organizações turísticas. Direcções futuras de investigação são igualmente apresentadas.

Analysis of multirate behavior in electronic systems

Esta tese insere-se na área da simulação de circuitos de RF e microondas, e visa o estudo de ferramentas computacionais inovadoras que consigam simular, de forma eficiente, circuitos não lineares e muito heterogéneos, contendo uma estrutura combinada de blocos analógicos de RF e de banda base e blocos digitais, a operar em múltiplas escalas de tempo. Os métodos numéricos propostos nesta tese baseiam-se em estratégias multi-dimensionais, as quais usam múltiplas variáveis temporais definidas em domínios de tempo deformados e não deformados, para lidar, de forma eficaz, com as disparidades existentes entre as diversas escalas de tempo. De modo a poder tirar proveito dos diferentes ritmos de evolução temporal existentes entre correntes e tensões com variação muito rápida (variáveis de estado activas) e correntes e tensões com variação lenta (variáveis de estado latentes), são utilizadas algumas técnicas numéricas avançadas para operar dentro dos espaços multi-dimensionais, como, por exemplo, os algoritmos multi-ritmo de Runge-Kutta, ou o método das linhas. São também apresentadas algumas estratégias de partição dos circuitos, as quais permitem dividir um circuito em sub-circuitos de uma forma completamente automática, em função dos ritmos de evolução das suas variáveis de estado. Para problemas acentuadamente não lineares, são propostos vários métodos inovadores de simulação a operar estritamente no domínio do tempo. Para problemas com não linearidades moderadas é proposto um novo método híbrido frequência-tempo, baseado numa combinação entre a integração passo a passo unidimensional e o método seguidor de envolvente com balanço harmónico. O desempenho dos métodos é testado na simulação de alguns exemplos ilustrativos, com resultados bastante promissores. Uma análise comparativa entre os métodos agora propostos e os métodos actualmente existentes para simulação RF, revela ganhos consideráveis em termos de rapidez de computação.

Arsénio em águas subterrâneas em Portugal

Da totalidade da água existente na Terra, somente 1% diz respeito à água subterrânea, recurso finito e essencial para a sobrevivência do Homem. A composição química das águas subterrâneas é influência não só pela actividade humana, mas também pelas características químicas do substrato geológico. As rochas e os sedimentos funcionam como meio de transporte e armazém desse tipo de águas e a presença de certos minerais pode originar a contaminação natural das águas subterrâneas em elementos prejudiciais para o Homem. Um desses elementos é o arsénio. Uma vez que a concentração máxima admitida de arsénio para água consumo humano diminuiu dos anteriores 50 mg.l-1 para os actuais 10 mg.l-1,a probabilidade de encontrar concentrações naturais em arsénio em águas subterrâneas acima desse limite aumentou. Este trabalho visa estudar e compreender, através de ensaios laboratoriais, a interacção água/rocha em termos químicos e a mobilidade/comportamento do arsénio em seis zonas de Portugal Continental com contextos geológicos diferentes (Vila Flor – Formação Filito-Quartzítica, Baião – granitos, Cacia – sedimentos do Cretácico, Mamodeiro – sedimentos do Neogénico, Escusa – calcários dolomíticos e Beja - gabros). Para o efeito, nas zonas de estudo mencionadas foram colhidas amostras de rocha ou sedimento e de água subterrânea que foram analisadas para 36 e 76 elementos, respectivamente. A componente mineralógica e as fases suporte do As foram estudadas recorrendo à análise por difracção de raios X, microssonda electrónica e à aplicação da extracção química selectiva sequencial. A mobilização do As foi avaliada através da realização de ensaios de coluna e de agitação com hidrocarbonetos derivados do petróleo. No primeiro ensaio, a circulação de diferentes soluções com diferentes velocidades lineares no interior de colunas de vidro cheias de rocha moída, avaliou a presença de As associado a fases hidrossolúveis e biodisponíveis. No segundo, a agitação de rocha moída com água e hidrocarbonetos derivados do petróleo pretendeu avaliar o comportamento do As quando um aquífero apresenta uma contaminação com aquele tipo de hidrocarbonetos.

Estabilidade de sistemas comutados com reset parcial

A presente tese é dedicada ao estudo da estabilidade de sistemas definidos por famílias finitas de sistemas lineares invariantes e por regras de comutação que coordenam a comutação entre eles. Assumimos que, em cada instante de tempo onde ocorre comutação, a trajectória do estado do sistema possa sofrer um "salto" desencadeado pela aplicação de um reset. Estes sistemas são, neste trabalho, designados por sistemas comutados com reset. Os resets podem ser de dois tipos - totais ou parciais, dependendo se a totalidade ou apenas uma parte das componentes do estado está disponível para reset. Neste sentido, distinguimos sistemas comutados com reset (total) e sistemas comutados com reset parcial. Analisamos a estabilidade dos dois tipos de sistemas comutados referidos à luz da teoria de Lyapunov e sob duas perspectivas; por um lado determinamos sob que condições um sistema comutado com reset é estável e por outro, identificamos resets que, quando aplicados, asseguram a estabilidade do sistema. Neste último ponto, a escolha dos resets adequados a aplicar pode por si só revelar-se insuficiente para obter estabilidade, especialmente se apenas parte das componentes do estado estiver disponível para reset (caso de reset parcial).

Measurements of thermophysical properties of perfluorocarbons and ionic liquids

No presente trabalho propõe-se estudar a tensão superficial de vários perfluorocarbonetos lineares, cíclicos, aromáticos e [alfa]-substituídos bem como líquidos iónicos com o catião imidazolium em comum. Apesar do seu interesse inerente, informação sobre esta propriedade para os compostos seleccionados é escassa e quando disponível apresenta discrepâncias consideráveis entre si. As medições foram realizadas no intervalo de temperaturas (283 to 353) K usando o método do anel de Du Noüy. Para os fluorocarbonetos, os dados experimentais demonstram que a estrutura molecular é o factor primordial no comportamento da superfície uma vez que os fluorocarbonetos aromáticos apresentam a tensão superficial mais elevada, seguida pelos fluorocarbonetos cíclicos e substituídos. Os perfluorocarbonetos lineares apresentam os menores valores de tensão superficial, aumentando ligeiramente com o aumento do número de carbonos. Os líquidos iónicos estudados foram seleccionados com o objectivo de fornecerem um estudo compreensivo sobre a influência do tamanho da cadeia alquílica do catião, o número de substituições no catião e a influência do anião. A influência do conteúdo de água na tensão superficial foi estudada em função da temperatura e da fracção molar de água para o liquido iónico mais hidrofóbico, [omim][PF6], e para o mais higroscópico, [bmim][PF6]. As funções termodinâmicas de superfície, como a entropia e entalpia de superfície, foram derivadas a partir da dependência da tensão superficial com a temperatura. Os dados obtidos para o fluorocarbonetos foram comparados com a correlação proposta por Faizullin, apresentando um desvio inferior a 4 % e demonstrando a sua aplicabilidade para com esta classe de compostos. A metodologia adoptada neste trabalho requer o conhecimento das densidades dos compostos de modo a aplicar a necessária correcção hidrostática. Contudo, para os líquidos iónicos esta informação é limitada ou mesmo inexistente. Por este motivo realizaram-se medições de densidade em função da pressão (0.10 < p/MPa < 10.0) e da temperatura (293.15 < T/K < 393.15). Desta dependência, as propriedades termodinâmicas, tais como compressibilidade isotérmica, expansividade isobárica, coeficiente térmico da pressão e dependência da capacidade calorífica com a pressão foram investigadas. A influência do teor de água na densidade foi também estudada para o líquido iónico mais hidrofóbico, [omim][PF6]. Um modelo simples de volume-ideal foi aplicado de forma preditiva para os volumes molares dos líquidos iónicos, em condições ambientais, descrevendo bem os dados experimentais. ABSTRACT: This work aims at studying the surface tension of some linear, cyclic, aromatic, [alfa]-substituted perfluorocarbons and imidazolium based ionic liquids. Despite its fundamental interest, information about this property for these compounds is scarce and the available data present strong discrepancies among each other. The measurements were carried out in the temperature range (283 to 353) K with the Du Noüy ring method. For the fluorocarbons, the analysis of the experimental data shows that the molecular structure is the main factor in the surface since the aromatic fluorocompounds present the highest surface tensions, followed by the cyclic and substituted fluorocompounds. The linear n-perfluoroalkanes exhibit the lowest surface tension values, slightly increasing with the carbon number. The set of selected ionic liquids was chosen to provide a comprehensive study of the influence of the cation alkyl chain length, the number of cation substitutions and the anion on the properties under study. The influence of water content in the surface tension was studied for several ILs as a function of the temperature as well as a function of water mole fraction, for the most hydrophobic IL investigated, [omim][PF6], and one hygroscopic IL, [bmim][PF6]. The surface thermodynamic functions such as surface entropy and enthalpy were derived from the temperature dependence of the surface tension values. The perfluorocarbons experimental data were compared against the Faizullin correlation, and it is shown that this correlation describes the measured surface tensions with deviations inferior to 4 %. The methodology adopted in this work requires the knowledge of the densities of the compounds under study in order to apply an hydrostatic correction. However, for ionic liquids these information is scarse and in some cases unavailable. Therefore, experimental measurements of the pressure (0.10 < p/MPa < 10.0) and temperature (293.15 < T/K < 393.15) dependence of the density and derived thermodynamic properties, such as the isothermal compressibility, the isobaric expansivity, the thermal pressure coefficient, and the pressure dependence of the heat capacity of several imidazolium-based ionic were determined. The influence of water content in the density was also studied for the most hydrophobic IL used, [omim][PF6]. A simple ideal-volume model was employed for the prediction of the imidazolium molar volumes at ambient conditions, which proved to agree well with the experimental results.

Caracterização de novas membranas de titanossilicatos microporosos por ensaios de permeabilidade

O interesse crescente das membranas inorgânicas deve-se à potencial aplicação em novas áreas de investigação e da indústria, e em alternativa a operações mais convencionais. Em particular, as membranas de titanossilicatos oferecem vantagens importantes sobre as de zeólitos, pois podem ser sintetizadas sem agentes estruturantes orgânicos, para evitar a calcinação subsequente usualmente responsável por defeitos irreversíveis, exibem novas possibilidades de substituição isomórfica da matriz, permitindo um ajuste mais fino das propriedades catalíticas e de adsorção, e são capazes de separar misturas com base em diferenças de afinidade e tamanho molecular (efeito de peneiro). Os objectivos principais deste trabalho foram: i) a caracterização dinâmica de membranas do tipo zeolítico sintetizadas no Laboratório Associado CICECO, realizando-se experiências de permeação com gases puros e misturas; ii) o desenvolvimento e validação de novos modelos para a transferência de massa multicomponente através de membranas porosas pela abordagem de Maxwell-Stefan, tendo em conta os mecanismos específicos encontrados, particularmente a contribuição por difusão superficial; e iii) a modelação dos pontos experimentais medidos, bem como dados compilados da literatura. De forma a realizar os ensaios de permeação, desenhou-se, montou-se e testou-se uma instalação experimental. Para gases puros, os objectivos principais foram a medição de permeâncias a temperatura constante, por variação da pressão transmembranar r ( ΔP ), e de permeâncias a temperatura programada, conduzidas a ΔP constante. Seguidamente, calcularam-se as selectividades ideais. Em relação a misturas, a determinação de selectividades reais requer as fracções molares no permeado e no retido. Na globalidade, estudaram-se três suportes diferentes (aço inoxidável e α − alumina) e dezanove membranas de AM-3, ETS-10, ZSM-5 e zeólito 4A, utilizando-se H2, He, N2, CO2, e O2. A primeira avaliação exploratória da qualidade das membranas foi feita permeando azoto à temperatura ambiente. Assim, permeâncias superiores a 10−6 mol/m2s.Pa evidenciavam defeitos grosseiros, levando-nos a efectuar cristalizações adicionais sobre as primeiras camadas. Este procedimento foi implementado com oito membranas. Um trabalho experimental mais detalhado foi conduzido com cinco membranas. Membranas com curvas permeância-temperatura ( Π −T ) decrescentes indicam tipicamente transporte viscoso e de Knudsen, i.e. meso e macrodefeitos. Por exemplo, a membrana nº 3 de AM-3 exibiu este comportamento com H2, He, N2 e CO2 puros. A contribuição de Knudsen foi confirmada pela relação linear encontrada entre as permeâncias e o inverso da raiz quadrada da massa molar. O mecanismo viscoso foi também identificado, pois as permeâncias eram inversamente proporcionais à viscosidade do gás ou, atendendo a equações do tipo de Chapman-Enskog, directamente proporcionais a 2 0.5 k d M (onde k d é o diâmetro cinético e M a massa molar). Um comportamento de permeação distinto observou-se com a membrana nº 5 de AM-3. As permeâncias registadas a temperatura programada eram aproximadamente constantes para o N2, CO2 e O2, enquanto com o H2 cresciam significativamente. Conjuntamente elas evidenciam a ocorrência de macro, meso e microdefeitos intercristalinos. O transporte gasoso activado através dos microporos compensa o impacto diminuidor dos meso e macroporos. Ao contrário do N2, CO2 e O2, o pequeno diâmetro do hidrogénio torna-lhe possível permear através dos microporos intracristalinos, o que lhe adiciona um mecanismo de transferência responsável por esse crescimento. No que respeita à difusão superficial, o sistema CO2/ZSM-5 pode ser tomado como um exemplo paradigmático. Uma vez que este zeólito adsorve o CO2, as permeâncias diminuem com o crescimento de ΔP , em virtude de as concentrações no sólido aumentarem de forma não linear e tenderem para a saturação. Os resultados contrastantes obtidos com azoto realçam ainda mais o mecanismo superficial, pois o N2 não é adsorvido e as permeâncias medidas são constantes. Globalmente, as selectividades ideais calculadas ( α* ) variam de cerca de 1 a 4.2. Este parâmetro foi também utilizado para discriminar as melhores membranas, uma vez que baixos valores de α* denotam o escoamento viscoso não-selectivo típico de macrodefeitos. Por exemplo, o H2/CO2 na membrana nº 3 de AM-3 apresentou α* = 3.6 − 4.2 para 40–120ºC, enquanto que na membrana nº 5 de AM-3 originou α* = 2.6 − 3.1. Estes resultados corroboraram as observações anteriores, segundo as quais a membrana nº 5 era melhor do que a nº 3. Alguns ensaios foram realizados com membranas saturadas com água para aumentar a selectividade: as medições mostraram claramente uma melhoria inicial seguida de uma redução consistente de α* com o aumento da temperatura, devido à remoção das moléculas de água responsáveis pela obstrução de alguns poros. Em relação às selectividades reais de misturas contendo hidrogénio, devem ser realizadas mais experiências e a quantificação do hidrogénio deve ser melhorada. No que concerne à modelação, novos factores termodinâmicos de Maxwell- Stefan foram derivados para as isotérmicas mono e multicomponente de Nitta, Langmuir-Freundlich e Toth, tendo sido testadas com dados de equilíbrio e de permeação da literatura. (É importante realçar que só estão publicadas equações para Langmuir e Dual-Site Langmuir de componentes puros e misturas). O procedimento de validação adoptado foi exigente: i) as isotérmicas multicomponente foram previstas a partir das de gás puro; ii) os parâmetros de difusão dos componentes puros foram ajustados a dados de permeação de cada gás; iii) depois, as difusividades cruzadas de Maxwell- Stefan foram estimadas pela relação de Vignes; finalmente, v) as novas equações foram testadas usando-se estes parâmetros, tendo sido capazes de estimar com sucesso fluxos binários. Paralelamente ao enfoque principal do trabalho, derivou-se um novo modelo para permuta iónica em materiais microporosos baseado nas equações de Maxwell-Stefan. Este foi validado com dados experimentais de remoção de Hg2+ e Cd2+ de soluções aquosas usando ETS-4. A sua capacidade preditiva foi também avaliada, sendo possível concluir que se comporta muito bem. Com efeito, conseguiram-se boas previsões com parâmetros optimizados a partir de conjuntos de dados independentes. Este comportamento pode ser atribuído aos princípios físicos sólidos da teoria de Maxwell-Stefan.

Phase equilibria modeling for biofuels production

Em consequência de uma série de problemas ambientais, económicos e políticos relacionados com o uso de combustíveis convencionais, vários países estão agora a focar as suas atenções em combustíveis alternativos. O biodiesel está na linha da frente das alternativas ao petróleo no sector dos transportes, sendo considerado uma opção a curto prazo visto que o seu preço é competitivo e não são necessárias mudanças nos motores para implementar o seu uso. De entres os possíveis processos de produzir biodiesel, a reacção de transesterificação com catálise básica é o método preferido. Depois da reacção são sempre necessários processos de purificação de modo ao biodiesel produzido cumprir os standards definidos para os combustíveis alternativos, reduzindo problemas de motor e consequentemente aumentando a sua aceitação por parte dos consumidores. De entre as especificações encontram-se o conteúdo em água, em álcool e em glicerol. Ser-se capaz de descrever correctamente o equilíbrio de fases de sistemas que são de interesse para os processos de purificação de biodieseis numa gama alargada de condições termodinâmicas é uma condição necessária para uma correcta simulação do processo industrial, de modo a se atingir uma elevada produtividade a baixos custos de operação. O uso de moléculas oxigenadas como combustíveis representa uma alteração significativa em termos da termodinâmica de soluções. Para combustíveis baseados em petróleo as equações de estado cúbicas e os modelos clássicos de coeficientes de actividade mostraram ser apropriados, no entanto para combustíveis novos como o biodiesel, sendo mais complexos do ponto de vista das interacções intermoleculares com formação de dipolos e pontes de hidrogénio, são necessários modelos termodinâmicos mais complexos para descrever essas interacções. Neste trabalho a CPA EoS (Cubic-Plus-Association Equation of State) será desenvolvida de modo a permitir uma descrição adequada dos equilíbrios líquido-vapor e líquido-líquido para uma serie de sistemas binários e multicomponentes contendo água, ácidos gordos, ésteres de ácidos gordos, glicerol e álcoois.

Subspace-based techniques and applications

Este trabalho focou-se no estudo de técnicas de sub-espaço tendo em vista as aplicações seguintes: eliminação de ruído em séries temporais e extracção de características para problemas de classificação supervisionada. Foram estudadas as vertentes lineares e não-lineares das referidas técnicas tendo como ponto de partida os algoritmos SSA e KPCA. No trabalho apresentam-se propostas para optimizar os algoritmos, bem como uma descrição dos mesmos numa abordagem diferente daquela que é feita na literatura. Em qualquer das vertentes, linear ou não-linear, os métodos são apresentados utilizando uma formulação algébrica consistente. O modelo de subespaço é obtido calculando a decomposição em valores e vectores próprios das matrizes de kernel ou de correlação/covariância calculadas com um conjunto de dados multidimensional. A complexidade das técnicas não lineares de subespaço é discutida, nomeadamente, o problema da pre-imagem e a decomposição em valores e vectores próprios de matrizes de dimensão elevada. Diferentes algoritmos de préimagem são apresentados bem como propostas alternativas para a sua optimização. A decomposição em vectores próprios da matriz de kernel baseada em aproximações low-rank da matriz conduz a um algoritmo mais eficiente- o Greedy KPCA. Os algoritmos são aplicados a sinais artificiais de modo a estudar a influência dos vários parâmetros na sua performance. Para além disso, a exploração destas técnicas é extendida à eliminação de artefactos em séries temporais biomédicas univariáveis, nomeadamente, sinais EEG.

Control of uncertain compartmental systems

Os sistemas compartimentais são frequentemente usados na modelação de diversos processos em várias áreas, tais como a biomedicina, ecologia, farmacocinética, entre outras. Na maioria das aplicações práticas, nomeadamente, aquelas que dizem respeito à administração de drogas a pacientes sujeitos a cirurgia, por exemplo, a presença de incertezas nos parâmetros do sistema ou no estado do sistema é muito comum. Ao longo dos últimos anos, a análise de sistemas compartimentais tem sido bastante desenvolvida na literatura. No entanto, a análise da sensibilidade da estabilidade destes sistemas na presença de incertezas tem recebido muito menos atenção. Nesta tese, consideramos uma lei de controlo por realimentação do estado com restrições de positividade e analisamos a sua robustez quando aplicada a sistemas compartimentais lineares e invariantes no tempo com incertezas nos parâmetros. Além disso, para sistemas lineares e invariantes no tempo com estado inicial desconhecido, combinamos esta lei de controlo com um observador do estado e a robustez da lei de controlo resultante também é analisada. O controlo do bloqueio neuromuscular por meio da infusão contínua de um relaxante muscular pode ser modelado como um sistema compartimental de três compartimentos e tem sido objecto de estudo por diversos grupos de investigação. Nesta tese, os nossos resultados são aplicados a este problema de controlo e são fornecidas estratégias para melhorar os resultados obtidos.