Análise de guias de ondas pelos métodos vetorial magnético e dos elementos-finitos

Neste trabalho, é apresentada uma formulação apropriada à análise de guias de ondas eletromagnéticos, cobrindo do espectro de microondas até o da óptica. Nas regiões a partir do ultravioleta, os comprimentos de onda são equivalentes às dimensões atômicas e a formulação necessita de uma abordagem quântica, que não é considerada neste estudo. A formulação é fundamentada nos métodos vetorial magnético e dos elementos finitos (MEF), em meios não homogêneos, anisotrópicos e não dissipativos, embora a dissipação possa ser facilmente introduzida na análise. Deu-se preferência à formulação com o campo magnético em vez do elétrico, pelo fato do campo magnético ignorar descontinuidades elétricas. Ele é contínuo em regiões de permeabilidade homogênea, propriedade dos meios dielétricos em geral ( = 0), independente da permissividade dos respectivos meios, conquanto os campos elétricos sejam descontínuos entre regiões de permissividades diferentes.

Análise de transição entre cabos coaxiais constituída por corrugações e desníveis angulares.

Esta dissertação apresenta um formalismo baseado no Método dos Elementos Finitos (MEF), adequado a análise de descontinuidades coaxiais com simetria axial, entre duas linhas coaxiais quaisquer, incluindo corrugações nos tubos internos e externos. Pelo método de Galerkin-Budnov deduz-se um operador integral bilinear, aplicado ao campo magnético, expandido em todo o domínio da estrutura, coaxial-descontinuidade-coaxial. As portas de entrada e saída da estrutura são posicionadas distantes da descontinuidade, de forma que nelas só haja o modo Transversal Eletromagnético (TEM). O campo magnético procurado e obtido pelo MEF. Os resultados encontrados; perdas de retorno, comportamento do campo magnético e as equi-fases nas portas de entrada e saída da estrutura, foram calculados e confrontados com as do Método de Casamento de Modos (MCM), com um alto grau de concordância.

Um estudo de caso envolvendo a retroanálise de tensões residuais de cravação.

As tensões residuais de cravação ocorrem quando a estaca e o solo não são totalmente descomprimidos após dissipação completa da energia transferida à estaca durante sua instalação. Este assunto tem merecido a atenção de vários pesquisadores, uma vez que a presença das tensões residuais na estaca causa uma alteração no seu comportamento quando carregada. A presente dissertação reuniu estudos anteriores que estabeleceram os principais fatores que influenciam as tensões residuais de cravação. Procurou-se executar uma série de simulações objetivando melhor entender o desenvolvimento das tensões residuais com a porcentagem de resistência de ponta, em continuidade aos estudos de Costa (1994). Costa (1994) observou que a porcentagem de carga na ponta é um fator relevante na avaliação das tensões residuais de cravação. Uma análise paramétrica efetuada confirmou estudos anteriores, verificando que a razão entre a carga residual na ponta em relação à capacidade de carga global cresce, a medida que a porcentagem de carga na ponta aumenta, chegando a um valor máximo para, em seguida, diminuir. Este comportamento é similar ao da curva de compactação do solo, quando, ao aumentar a umidade, o peso específico seco aumenta, até um valor máximo, correspondente à umidade ótima, para em seguida reduzir. Verificouse, ainda, semelhantemente ao aumento da energia de compactação, que vai transladando a curva para cima e para à esquerda do gráfico umidade x peso específico seco, que o aumento do comprimento da estaca apresenta um comportamento similar. O aumento do comprimento leva a curva para cima e para a esquerda. Finalmente, selecionou-se um caso de obra com condições propícias para o desenvolvimento de altas tensões residuais de cravação. Trata-se de um caso de estacas metálicas longas, embutidas em solo residual jovem, de elevada resistência, bem documentado com extensa instrumentação, por ocasião da instalação. A retroanálise de cinco estacas do banco de dados mostrou que as cargas residuais previstas, num programa de simulação de cravação, se aproximaram muito dos valores experimentais. Foi observado também que a profundidade do ponto neutro previsto e medido apresentou uma excelente concordância. O resultado mais relevante desta pesquisa foi quando os valores de porcentagem de ponta foram introduzidos no eixo das abscissas e os valores da razão entre a carga residual na ponta e a capacidade de carga global, no eixo das ordenadas, e se observou o aspecto da curva semelhante à de compactação. A dissertação ilustra ser possível e simples a previsão das tensões residuais de cravação através de uma análise pela equação da onda.

Implementação de um algoritmo numérico para solução da equação de Christoffel generalizada em acustoelasticidade.

Extensos estudos realizados nas últimas décadas sobre a propagação de ondas ultrassônicas em sólidos levaram ao desenvolvimento de técnicas não destrutivas para a avaliação da segurança e integridade de estruturas e componentes industriais. O interesse na aplicação de técnicas ultrassônicas para medição de tensões aplicadas e residuais decorre da mudança mensurável da velocidade das ondas ultrassônicas na presença de um campo de tensões, fenômeno conhecido como efeito acustoelástico. Uma teoria de acustoelasticidade fornece um meio atrativo e não destrutivo de medir a tensão média ao longo do caminho percorrido pela onda. O estudo da propagação das ondas ultrassônicas em meios homogêneos anisotrópicos sob tensão conduz a um problema não linear de autovalores dado pela equação de Christoffel generalizada. A característica não linear deste problema decorre da interdependência entre as constantes elásticas efetivas do material e as tensões atuantes. A medição experimental de tensões por técnicas ultrassônicas é um problema inverso da acustoelasticidade. Esta dissertação apresenta a implementação de um algoritmo numérico, baseado no método proposto por Degtyar e Rokhlin, para solução do problema inverso da acustoelasticidade em sólidos ortotrópicos sujeitos a um estado plano de tensões. A solução da equação de Christoffel generalizada apresenta dificuldades de natureza numérica e prática. A estabilidade e a precisão do algoritmo desenvolvido, bem como a influência das incertezas na medição experimental das velocidades das ondas ultrassônicas, foram então investigadas. Dados sintéticos para as velocidades das ondas ultrassônicas de incidência oblíqua em uma placa sujeita a um estado plano de tensões foram gerados pela solução direta da equação de Christoffel generalizada para ilustrar a aplicação do algoritmo desenvolvido. O objetivo maior desta dissertação é a disponibilização de uma nova ferramenta de cálculo para suporte às atividades experimentais de medição de tensões por ultrassom no país.

A internacionalização privada na onda estatizante: as relações público-privadas na internacionalização das empreiteiras brasileiras entre 1974-1979

Durante a década de 1970, ocorreu, no Brasil, a fase pioneira de internacionalização das empresas privadas brasileiras do setor de prestação de serviços de engenharia de construção. Essa internacionalização deu-se, apesar da ausência de uma política estruturada do setor público. A incapacidade estatal de perceber esse processo como algo profícuo em seus planos, impediu o governo de criar uma política consolidada para a multiplicação dessas empresas no exterior. Essa miopia estatal é percebida tanto na esfera da política doméstica, voltada para manter uma política de substituição de importação dentro das fronteiras, como, na política externa, ao desenvolver uma internacionalização das empresas públicas ao invés das do setor privado. Essa situação pode ser percebida no caso analisado da empresa privada brasileira Mendes Júnior, que ingressou no mercado iraquiano mais por conta própria do que em decorrência de uma política estatal brasileira.

Sobre a influência da dimensionalidade do espaço no átomo de hidrogênio relativístico

Faz-se uma revisão do problema da dimensionalidade do espaço entendido como um problema de Física, enfatizando que algumas leis físicas dependem fortemente deste parâmetro topológico do espaço. Discute-se o que já foi feito tanto no caso da equação de Schrödinger quanto na de Dirac. A situação na literatura é bastante controversa e, no caso específico da equação de Dirac em D dimensões, não se encontra nenhum trabalho na literatura científica que leve em conta o potencial de intera coulombiana corretamente generalizado quando o número de dimensões espaciais é maior do que três. Discute-se, portanto, o átomo de hidrogênio relativístico em D dimensões. Novos resultados numéricos para os níveis de energia e para as funções de onda são apresentados e discutidos. Em particular, considera-se a possibilidade de existência de átomos estáveis em espaços com dimensionalidade 6= 3.

Análise de erros da equação de advecção unidimensional no Método de Volumes Finitos

Uma análise utilizando a série de Taylor é apresentada para se estimar a priori os erros envolvidos na solução numérica da equação de advecção unidimensional com termo fonte, através do Método dos Volumes Finitos em uma malha do tipo uniforme e uma malha não uniforme. Também faz-se um estudo a posteriori para verificar a magnitude do erro de discretização e corroborar os resultados obtidos através da análise a priori. Por meio da técnica de solução manufaturada tem-se uma solução analítica para o problema, a qual facilita a análise dos resultados numéricos encontrados, e estuda-se ainda a influência das funções de interpolação UDS e CDS e do parâmetro u na solução numérica.

Velocidade da onda de pulso em adultos jovens acompanhados por 18 anos: impacto de variáveis pressóricas, antropométricas, metabólicas, inflamatórias e de função endotelial. Estudo do Rio de Janeiro.

Dados sobre a avaliação não invasiva da rigidez vascular e suas relações com variáveis de risco cardiovascular são escassos em jovens. Objetiva avaliar a relação entre a velocidade de onda de pulso (VOP) e a pressão arterial (PA), variáveis antropométricas, metabólicas, inflamatórias e de disfunção endotelial em indivíduos adultos jovens. Foram estudados 96 indivíduos (51 homens) do Estudo do Rio de Janeiro, em duas avaliações, A1 e A2, com intervalo de 17,691,58 anos (16 a 21 anos). Em A1 foram avaliados em suas escolas (10-15 anos - média 12,421,47 anos) e em A2 foram novamente avaliados em nível ambulatorial (26-35 anos - média 30,091,92 anos). Em A1 foram obtidos pressão arterial (PA) e índice de massa corporal (IMC). Em A2 foram obtidos a velocidade da onda de pulso (VOP)-método Complior, PA, IMC, circunferência abdominal (CA), glicose, perfil lipídico, leptina, insulina, adiponectina, o índice de resistência à insulina HOMA-IR, proteína C-Reativa ultrassensível (PCRus) e as moléculas de adesão E-selectina, Vascular Cell Adhesion Molecule-1(VCAM-1) e Intercellular Adhesion Molecule-1 (ICAM-1). Foram obtidos, ainda, a variação da PA e do IMC entre as 2 avaliações. Em A2 os indivíduos foram estratificados segundo o tercil da VOP para cada sexo. Como resultados temos: 1) Os grupos foram constituídos da seguinte forma: Tercil 1:homens com VOP < 8,69 m/s e mulheres com VOP < 7,66 m/s; Tercil 2: homens com VOP ≥ 8,69 m/s e < 9,65m/s e mulheres com VOP ≥ 7,66 m/s e < 8,31m/s;Tercil 3:homens com VOP ≥ 9,65 m/s e mulheres com VOP ≥ 8,31 m/s. 2) O grupo com maior tercil de VOP mostrou maiores médias de PA sistólica (PAS) (p=0,005), PA diastólica (PAD) (p=0,007), PA média (PAM) (p=0,004), variação da PAD (p=0,032), variação da PAM (p=0,003), IMC (p=0,046), variação do IMC (p=0,020), insulina (p=0,019), HOMA-IR (p=0,021), E-selectina (p=0,032) e menores médias de adiponectina (p=0,016), além de maiores prevalências de diabetes mellitus/intolerância à glicose (p=0,022) e hiperinsulinemia (p=0,038); 3) Houve correlação significativa e positiva da VOP com PAS (p<0,001), PAD (p<0,001), PP (p=0,048) e PAM (p<0,001) de A2, com a variação da pressão arterial (PAS, PAD e PAM) (p<0,001) entre as duas avaliações, com o IMC de A2 (p=0,005) e com a variação do IMC (p<0,001) entre as duas avaliações, com CA (p=0,001), LDLcolesterol (p=0,049) e E-selectina (p<0,001) e correlação negativa com HDLcolesterol (p<0,001) e adiponectina (p<0,001); 4)Em modelo de regressão múltipla, após ajuste do HDL-colesterol, LDLcolesterol e adiponectina para sexo, idade, IMC e PAM, apenas o sexo masculino e a PAM mantiveram correlação significativa com a VOP. A VOP em adultos jovens mostrou relação significativa com variáveis de risco cardiovascular, destacando-se o sexo masculino e a PAM como importantes variáveis no seu determinismo. Os achados sugerem que a medida da VOP pode ser útil para a identificação do acometimento vascular nessa faixa etária.

Modelagem e simulação do escoamento imiscível em meios porosos fractais descritos pela equação de Kozeny-Carman Generalizada

O presente trabalho trata do escoamento bifásico em meios porosos heterogêneos de natureza fractal, onde os fluidos são considerados imiscíveis. Os meios porosos são modelados pela equação de Kozeny-Carman Generalizada (KCG), a qual relaciona a porosidade com a permeabilidade do meio através de uma nova lei de potência. Esta equação proposta por nós é capaz de generalizar diferentes modelos existentes na literatura e, portanto, é de uso mais geral. O simulador numérico desenvolvido aqui emprega métodos de diferenças finitas. A evolução temporal é baseada em um esquema de separação de operadores que segue a estratégia clássica chamada de IMPES. Assim, o campo de pressão é calculado implicitamente, enquanto que a equação da saturação da fase molhante é resolvida explicitamente em cada nível de tempo. O método de otimização denominado de DFSANE é utilizado para resolver a equação da pressão. Enfatizamos que o DFSANE nunca foi usado antes no contexto de simulação de reservatórios. Portanto, o seu uso aqui é sem precedentes. Para minimizar difusões numéricas, a equação da saturação é discretizada por um esquema do tipo "upwind", comumente empregado em simuladores numéricos para a recuperação de petróleo, o qual é resolvido explicitamente pelo método Runge-Kutta de quarta ordem. Os resultados das simulações são bastante satisfatórios. De fato, tais resultados mostram que o modelo KCG é capaz de gerar meios porosos heterogêneos, cujas características permitem a captura de fenômenos físicos que, geralmente, são de difícil acesso para muitos simuladores em diferenças finitas clássicas, como o chamado fenômeno de dedilhamento, que ocorre quando a razão de mobilidade (entre as fases fluidas) assume valores adversos. Em todas as simulações apresentadas aqui, consideramos que o problema imiscível é bidimensional, sendo, portanto, o meio poroso caracterizado por campos de permeabilidade e de porosidade definidos em regiões Euclideanas. No entanto, a teoria abordada neste trabalho não impõe restrições para sua aplicação aos problemas tridimensionais.

Solução analítica da equação unidimensional de transporte de nêutrons monoenergéticos com espalhamento linearmente anisotrópico e aproximação sintética de difusão

Nesta dissertação, são apresentados os seguintes modelos matemáticos de transporte de nêutrons: a equação linearizada de Boltzmann e a equação da difusão de nêutrons monoenergéticos em meios não-multiplicativos. Com o objetivo de determinar o período fluxo escalar de nêutrons, é descrito um método espectronodal que gera soluções numéricas para o problema de difusão em geometria planar de fonte fixa, que são livres de erros de truncamento espacial, e que conjugado com uma técnica de reconstrução espacial intranodal gera o perfil detalhado da solução. A fim de obter o valor aproximado do fluxo angular de nêutrons em um determinado ponto do domínio e em uma determinada direção de migração, descreve-se também um método de reconstrução angular baseado na solução analítica da equação unidimensional de transporte de nêutrons monoenergéticos com espalhamento linearmente anisotrópico com aproximação sintética de difusão nos termos de fonte por espalhamento. O código computacional desenvolvido nesta dissertação foi implementado na plataforma livre Scilab, e para ilustrar a eficiência do código criado,resultados numéricos obtidos para três problemas-modelos são apresentados

Avaliação dos algoritmos de Picard-Krylov e Newton-Krylov na solução da equação de Richards

A engenharia geotécnica é uma das grandes áreas da engenharia civil que estuda a interação entre as construções realizadas pelo homem ou de fenômenos naturais com o ambiente geológico, que na grande maioria das vezes trata-se de solos parcialmente saturados. Neste sentido, o desempenho de obras como estabilização, contenção de barragens, muros de contenção, fundações e estradas estão condicionados a uma correta predição do fluxo de água no interior dos solos. Porém, como a área das regiões a serem estudas com relação à predição do fluxo de água são comumente da ordem de quilômetros quadrados, as soluções dos modelos matemáticos exigem malhas computacionais de grandes proporções, ocasionando sérias limitações associadas aos requisitos de memória computacional e tempo de processamento. A fim de contornar estas limitações, métodos numéricos eficientes devem ser empregados na solução do problema em análise. Portanto, métodos iterativos para solução de sistemas não lineares e lineares esparsos de grande porte devem ser utilizados neste tipo de aplicação. Em suma, visto a relevância do tema, esta pesquisa aproximou uma solução para a equação diferencial parcial de Richards pelo método dos volumes finitos em duas dimensões, empregando o método de Picard e Newton com maior eficiência computacional. Para tanto, foram utilizadas técnicas iterativas de resolução de sistemas lineares baseados no espaço de Krylov com matrizes pré-condicionadoras com a biblioteca numérica Portable, Extensible Toolkit for Scientific Computation (PETSc). Os resultados indicam que quando se resolve a equação de Richards considerando-se o método de PICARD-KRYLOV, não importando o modelo de avaliação do solo, a melhor combinação para resolução dos sistemas lineares é o método dos gradientes biconjugados estabilizado mais o pré-condicionador SOR. Por outro lado, quando se utiliza as equações de van Genuchten deve ser optar pela combinação do método dos gradientes conjugados em conjunto com pré-condicionador SOR. Quando se adota o método de NEWTON-KRYLOV, o método gradientes biconjugados estabilizado é o mais eficiente na resolução do sistema linear do passo de Newton, com relação ao pré-condicionador deve-se dar preferência ao bloco Jacobi. Por fim, há evidências que apontam que o método PICARD-KRYLOV pode ser mais vantajoso que o método de NEWTON-KRYLOV, quando empregados na resolução da equação diferencial parcial de Richards.

Roteamento e alocação de comprimento de onda em redes WDM segundo algoritmo baseado em regras da natureza.

O surgimento de novos serviços de telecomunicações tem provocado um enorme aumento no tráfego de dados nas redes de transmissão. Para atender a essa demanda crescente, novas tecnologias foram desenvolvidas e implementadas ao longo dos anos, sendo que um dos principais avanços está na área de transmissão óptica, devido à grande capacidade de transporte de informação da fibra óptica. A tecnologia que melhor explora a capacidade desse meio de transmissão atualmente é a multiplexação por divisão de comprimento de onda ou Wavelength Division Multiplexing (WDM) que permite a transmissão de diversos sinais utilizando apenas uma fibra óptica. Redes ópticas WDM se tornaram muito complexas, com enorme capacidade de transmissão de informação (terabits por segundo), para atender à explosão de necessidade por largura de banda. Nesse contexto, é de extrema importância que os recursos dessas redes sejam utilizados de forma inteligente e otimizada. Um dos maiores desafios em uma rede óptica é a escolha de uma rota e a seleção de um comprimento de onda disponível na rede para atender uma solicitação de conexão utilizando o menor número de recursos possível. Esse problema é bastante complexo e ficou conhecido como problema de roteamento e alocação de comprimento de onda ou, simplesmente, problema RWA (Routing and Wavelentgh Assignment problem). Muitos estudos foram realizados com o objetivo de encontrar uma solução eficiente para esse problema, mas nem sempre é possível aliar bom desempenho com baixo tempo de execução, requisito fundamental em redes de telecomunicações. A técnica de algoritmo genético (AG) tem sido utilizada para encontrar soluções de problemas de otimização, como é o caso do problema RWA, e tem obtido resultados superiores quando comparada com soluções heurísticas tradicionais encontradas na literatura. Esta dissertação apresenta, resumidamente, os conceitos de redes ópticas e de algoritmos genéticos, e descreve uma formulação do problema RWA adequada à solução por algoritmo genético.

Posicionamento aproximado do estado final para sistemas térmicos descritos pela equação do calor.

Neste trabalho, será considerado um problema de controle ótimo quadrático para a equação do calor em domínios retangulares com condição de fronteira do tipo Dirichlet é nos quais, a função de controle (dependente apenas no tempo) constitui um termo de fonte. Uma caracterização da solução ótima é obtida na forma de uma equação linear em um espaço de funções reais definidas no intervalo de tempo considerado. Em seguida, utiliza-se uma sequência de projeções em subespaços de dimensão finita para obter aproximações para o controle ótimo, o cada uma das quais pode ser gerada por um sistema linear de dimensão finita. A sequência de soluções aproximadas assim obtidas converge para a solução ótima do problema original. Finalmente, são apresentados resultados numéricos para domínios espaciais de dimensão 1.

Apocalipse não: o estado do bem-estar social sobrevive à onda neoliberal.

Esta Tese aborda a questão do Estado e dos direitos de cidadania, particularmente dos direitos sociais e políticas a eles correspondentes, que constituem o que se convencionou tratar por Estado do bem-estar social, ou welfare state. Nela é analisado o processo recente de contradição entre as idéias e políticas neoliberais e o sistemas de welfare state, procurando verificar o impacto daquelas políticas sobre os sistemas. O objetivo central é mostrar que a ofensiva neoliberal vem sendo mal sucedida em sua tentativa de desarticular as instituições dos direitos e políticas sociais no que diz respeito às sociedades mais complexas e desenvolvidas. A principal razão disto deriva da própria história do welfare state, que corresponde a um longo processo de transformação da relação política entre Estado e sociedade. O estudo do impacto efetivo das políticas neoliberais sobre os sistemas de bem-estar ou proteção social, a que se convencionou tratar por welfare state está centrada no que vem ocorrendo nos países desenvolvidos, fazendo um contraponto com o que se passa no Brasil. Este tema tem sido objeto de preocupação de um grande número de analistas no exterior e no Brasil em função tanto do conteúdo das propostas neoliberais em relação às políticas sociais fortalecimento das soluções de mercado, restrição e das políticas públicas de proteção social nos grupos mais pobres da população quanto do sucesso avassalador que o neoliberalismo alcançou nos campos da economia e do trabalho. A Tese defende a idéia de que a sociedade brasileira e seu welfare state chegaram a um nível tal de complexidade, que não pode mais ser compreendida com base num ferramental analítico que tende a reduzir suas singularidades à configurações muito amplas ora denominadas por periferia do sistema capitalista, ora por realidade latino-americana, as contribuem para simplificar a análise e dificultar o aprofundamento da discussão em relação ao conteúdo, ao rumo e ao alcance do processo político e das políticas sociais.

Formação de nanopadrões em superfícies por sputtering iônico: Estudo numérico da equação anisotrópica amortecida de Kuramoto-Sivashinsky.

Apresenta-se uma abordagemnumérica para ummodelo que descreve a formação de padrões por sputtering iônico na superfície de ummaterial. Esse processo é responsável pela formação de padrões inesperadamente organizados, como ondulações, nanopontos e filas hexagonais de nanoburacos. Uma análise numérica de padrões preexistentes é proposta para investigar a dinâmica na superfície, baseada em ummodelo resumido em uma equação anisotrópica amortecida de Kuramoto-Sivashinsky, em uma superfície bidimensional com condições de contorno periódicas. Apesar de determinística, seu caráter altamente não-linear fornece uma rica gama de resultados, sendo possível descrever acuradamente diferentes padrões. Umesquema semi implícito de diferenças finitas com fatoração no tempo é aplicado na discretização da equação governante. Simulações foram realizadas com coeficientes realísticos relacionados aos parâmetros físicos (anisotropias, orientação do feixe, difusão). A estabilidade do esquema numérico foi analisada por testes de passo de tempo e espaçamento de malha, enquanto a verificação do mesmo foi realizada pelo Método das Soluções Manufaturadas. Ondulações e padrões hexagonais foram obtidos a partir de condições iniciais monomodais para determinados valores do coeficiente de amortecimento, enquanto caos espaço-temporal apareceu para valores inferiores. Os efeitos anisotrópicos na formação de padrões foramestudados, variando o ângulo de incidência.