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.

Networked tourism: from world destinations to inter-organizations

Na última década, a referência ao conceito de redes cresceu rapidamente entre a literatura sobre turismo, geralmente aplicado a tópicos como as interorganizações, estrutura de multi-destinos, espaços de Turismo online, entre outros. O conceito de rede difundiu-se na natureza e na sociedade, em áreas que vão desde a Biologia à Medicina, ou da Economia à Gestão, e o conhecimento sobre redes tem vindo a impulsionar uma teoria comum para facilitar a compreensão de diferentes sistemas complexos e a representação das ligações entre organizações, acções, bens, proteínas ou pessoas. A tese teve como propósito o encontro de um eixo comum entre dois campos férteis de investigação através de uma revisão teórica sistemática. A investigação sobre redes complexas é um campo recente na Física que tem vindo a desenvolver-se bastante na última década com fortes aplicações interdisciplinares. Por outro lado, a análise de redes sociais é uma área de investigação activa em Sociologia e Economia há bastante tempo. O estudo das implicações das redes complexas para a ciência das redes de turismo é uma área promissora já com resultados fascinantes. A tese tem três resultados principais. Primeiro, traz conhecimento das ricas áreas de conhecimento sobre redes complexas e redes sociais. Em segundo lugar, apresenta modelos evolutivos que melhor se adaptam às chegadas turísticas internacionais. Como se organizam as redes sociais? Como é que os indivíduos escolhem os seus destinos de viagem? Estes são exemplos de questões que serão abordadas na tese. Em terceiro lugar, discute resultados que fazem notar comportamentos comuns entre redes em turismo e outras redes reais. O que é comum a todas as redes na natureza? Adicionalmente, os padrões encontrados entre os destinos turísticos mostram um comportamento não social, com destinos mais característicos de redes económicas e sistemas tecnológicos que questionam a faceta social do sector do turismo. Por acréscimo, a rede de transportes aéreos e a rede de turismo mostram diferenças consideráveis que se podem dever a razões políticas ou outras que provavelmente explicam o aumento da utilização de voos charters.

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.

Estudo e aplicação de amplificadores ópticos para suporte a redes banda larga

Esta tese tem como foco principal a análise dos principais tipos de amplificação óptica e algumas de suas aplicações em sistemas de comunicação óptica. Para cada uma das tecnologias abordadas, procurou-se definir o estado da arte bem como identificar as oportunidades de desenvolvimento científico relacionadas. Os amplificadores para os quais foi dirigido alguma atenção neste documento foram os amplificadores em fibra dopada com Érbio (EDFA), os amplificadores a semicondutor (SOA) e os amplificadores de Raman (RA). Este trabalho iniciou-se com o estudo e análise dos EDFA’s. Dado o interesse científico e económico que estes amplificadores têm merecido, apenas poucos nichos de investigação estão ainda em aberto. Dentro destes, focá-mo-nos na análise de diferentes perfis de fibra óptica dopada de forma a conseguir uma optimização do desempenho dessas fibras como sistemas de amplificação. Considerando a fase anterior do trabalho como uma base de modelização para sistemas de amplificação com base em fibra e dopantes, evoluiu-se para amplificadores dopados mas em guias de onda (EDWA). Este tipo de amplificador tenta reduzir o volume físico destes dispositivos, mantendo as suas características principais. Para se ter uma forma de comparação de desempenho deste tipo de amplificador com os amplificadores em fibra, foram desenvolvidos modelos de caixa preta (BBM) e os seus parâmetros afinados por forma a termos uma boa modelização e posterior uso deste tipo de amplificiadores em setups de simulação mais complexos. Depois de modelizados e compreendidos os processo em amplificadores dopados, e com vista a adquirir uma visão global comparativa, foi imperativo passar pelo estudo dos processos de amplificação paramétrica de Raman. Esse tipo de amplificação, sendo inerente, ocorre em todas as bandas de propagação em fibra e é bastante flexível. Estes amplificadores foram inicialmente modelizados, e algumas de suas aplicações em redes passivas de acesso foram estudadas. Em especial uma série de requisitos, como por exemplo, a gama de comprimentos de onda sobre os quais existem amplificação e os altos débitos de perdas de inserção, nos levaram à investigação de um processo de amplificação que se ajustasse a eles, especialmente para buscar maiores capacidades de amplificação (nomeadamente longos alcances – superiores a 100 km – e altas razões de divisão – 1:512). Outro processo investigado foi a possibilidade de flexibilização dos parâmetros de comprimento de onda de ganho sem ter que mudar as caractísticas da bomba e se possível, mantendo toda a referenciação no transmissor. Este processo baseou-se na técnica de clamping de ganho já bastante estudada, mas com algumas modificações importantes, nomeadamente a nível do esquema (reflexão apenas num dos extremos) e da modelização do processo. O processo resultante foi inovador pelo recurso a espalhamentos de Rayleigh e Raman e o uso de um reflector de apenas um dos lados para obtenção de laser. Este processo foi modelizado através das equações de propagação e optimizado, tendo sido demonstrado experimentalmente e validado para diferentes tipos de fibras. Nesta linha, e dada a versatilidade do modelo desenvolvido, foi apresentada uma aplicação mais avançada para este tipo de amplificadores. Fazendo uso da sua resposta ultra rápida, foi proposto e analisado um regenerador 2R e analisada por simulação a sua gama de aplicação tendo em vista a sua aplicação sistémica. A parte final deste trabalho concentrou-se nos amplificadores a semiconductor (SOA). Para este tipo de amplificador, os esforços foram postos mais a nível de aplicação do que a nível de sua modelização. As aplicações principais para estes amplificadores foram baseadas em clamping óptico do ganho, visando a combinação de funções lógicas essenciais para a concepção de um latch óptico com base em componentes discretos. Assim, com base num chip de ganho, foi obtido uma porta lógica NOT, a qual foi caracterizada e demonstrada experimentalmente. Esta foi ainda introduzida num esquema de latching de forma a produzir um bi-estável totalmente óptico, o qual também foi demonstrado e caracterizado. Este trabalho é finalizado com uma conclusão geral relatando os subsistemas de amplificação e suas aplicacações.

Diffusion coefficients in dense and supercritical fluids

Os coeficientes de difusão (D 12) são propriedades fundamentais na investigação e na indústria, mas a falta de dados experimentais e a inexistência de equações que os estimem com precisão e confiança em fases comprimidas ou condensadas constituem limitações importantes. Os objetivos principais deste trabalho compreendem: i) a compilação de uma grande base de dados para valores de D 12 de sistemas gasosos, líquidos e supercríticos; ii) o desenvolvimento e validação de novos modelos de coeficientes de difusão a diluição infinita, aplicáveis em amplas gamas de temperatura e densidade, para sistemas contendo componentes muito distintos em termos de polaridade, tamanho e simetria; iii) a montagem e teste de uma instalação experimental para medir coeficientes de difusão em líquidos e fluidos supercríticos. Relativamente à modelação, uma nova expressão para coeficientes de difusão a diluição infinita de esferas rígidas foi desenvolvida e validada usando dados de dinâmica molecular (desvio relativo absoluto médio, AARD = 4.44%) Foram também estudados os coeficientes de difusão binários de sistemas reais. Para tal, foi compilada uma extensa base de dados de difusividades de sistemas reais em gases e solventes densos (622 sistemas binários num total de 9407 pontos experimentais e 358 moléculas) e a mesma foi usada na validação dos novos modelos desenvolvidos nesta tese. Um conjunto de novos modelos foi proposto para o cálculo de coeficientes de difusão a diluição infinita usando diferentes abordagens: i) dois modelos de base molecular com um parâmetro específico para cada sistema, aplicáveis em sistemas gasosos, líquidos e supercríticos, em que natureza do solvente se encontra limitada a apolar ou fracamente polar (AARDs globais na gama 4.26-4.40%); ii) dois modelos de base molecular biparamétricos, aplicáveis em todos os estados físicos, para qualquer tipo de soluto diluído em qualquer solvente (apolar, fracamente polar e polar). Ambos os modelos dão origem a erros globais entre 2.74% e 3.65%; iii) uma correlação com um parâmetro, específica para coeficientes de difusão em dióxido de carbono supercrítico (SC-CO2) e água líquida (AARD = 3.56%); iv) nove correlações empíricas e semi-empíricas que envolvem dois parâmetros, dependentes apenas da temperatura e/ou densidade do solvente e/ou viscosidade do solvente. Estes últimos modelos são muito simples e exibem excelentes resultados (AARDs entre 2.78% e 4.44%) em sistemas líquidos e supercríticos; e v) duas equações preditivas para difusividades de solutos em SC-CO2, em que os erros globais de ambas são inferiores a 6.80%. No global, deve realçar-se o facto de os novos modelos abrangerem a grande variedade de sistemas e moléculas geralmente encontrados. Os resultados obtidos são consistentemente melhores do que os obtidos com os modelos e abordagens encontrados na literatura. No caso das correlações com um ou dois parâmetros, mostrou-se que estes mesmos parâmetros podem ser ajustados usando um conjunto muito pequeno de dados, e posteriormente serem utilizados na previsão de valores de D 12 longe do conjunto original de pontos. Uma nova instalação experimental para medir coeficientes de difusão binários por técnicas cromatográficas foi montada e testada. O equipamento, o procedimento experimental e os cálculos analíticos necessários à obtenção dos valores de D 12 pelo método de abertura do pico cromatográfico, foram avaliados através da medição de difusividades de tolueno e acetona em SC-CO2. Seguidamente, foram medidos coeficientes de difusão de eucaliptol em SC-CO2 nas gamas de 202 – 252 bar e 313.15 – 333.15 K. Os resultados experimentais foram analisados através de correlações e modelos preditivos para D12.

Novel critical phenomena in optimization driven processes on networks

The work presented in this Ph.D thesis was developed in the context of complex network theory, from a statistical physics standpoint. We examine two distinct problems in this research field, taking a special interest in their respective critical properties. In both cases, the emergence of criticality is driven by a local optimization dynamics. Firstly, a recently introduced class of percolation problems that attracted a significant amount of attention from the scientific community, and was quickly followed up by an abundance of other works. Percolation transitions were believed to be continuous, until, recently, an 'explosive' percolation problem was reported to undergo a discontinuous transition, in [93]. The system's evolution is driven by a metropolis-like algorithm, apparently producing a discontinuous jump on the giant component's size at the percolation threshold. This finding was subsequently supported by number of other experimental studies [96, 97, 98, 99, 100, 101]. However, in [1] we have proved that the explosive percolation transition is actually continuous. The discontinuity which was observed in the evolution of the giant component's relative size is explained by the unusual smallness of the corresponding critical exponent, combined with the finiteness of the systems considered in experiments. Therefore, the size of the jump vanishes as the system's size goes to infinity. Additionally, we provide the complete theoretical description of the critical properties for a generalized version of the explosive percolation model [2], as well as a method [3] for a precise calculation of percolation's critical properties from numerical data (useful when exact results are not available). Secondly, we study a network flow optimization model, where the dynamics consists of consecutive mergings and splittings of currents flowing in the network. The current conservation constraint does not impose any particular criterion for the split of current among channels outgoing nodes, allowing us to introduce an asymmetrical rule, observed in several real systems. We solved analytically the dynamic equations describing this model in the high and low current regimes. The solutions found are compared with numerical results, for the two regimes, showing an excellent agreement. Surprisingly, in the low current regime, this model exhibits some features usually associated with continuous phase transitions.