Defeitos pontuais nos compostos intermetálicos ZrNi e Zr/sub 2/Ni estudados por dinâmica molecular

Empregamos a técnica de Dinâmica Molecular para estudar propriedades de defeitos pontuais nos compostos intermetálicos ZrNi e Zr2Ni. Descrevemos as configurações estáveis de defeitos e mecanismos de migração, assim como as energias envolvidas. Os potenciais interatômicos foram derivados do Embedded Atom Model. No intuito de levar em conta a variação de estequiometria causada pela presença de alguns tipos de defeitos em intermetálicos, apresentamos um método numérico que fornece a energia efetiva de formação de defeitos e aplicamos o método ao ZrNi e Zr2Ni. Os resultados mostraram que vacâncias são mais estáveis na sub—rede do Ni, com energia de formação de 0,83 e-0,61 eV em ZrNi e Zr2Ni, respectivamente. Vacâncias de Zr são instáveis em ambos compostos; elas decaem espontaneamente em pares anti—sítio e vacância de Ni. Configurações e energias de formação de intersticiais também foram calculadas e mostraram comportamentos similares. Em ZrNi, a migração de vacâncias ocorre preferencialmente nas direções [025] e [100], com as respectivas energias de migração 0,67 e 0,73 eV, e é um processo essencialmente bidimensional no plano (001). Em Zr,Ni, a migração de vacâncias é unidimensional, ocorrendo na direção [001], com energia de migração de 0,67 eV. Em ambos compostos a presença de defeitos de anti—sítio de Ni diminui a energia de migração da vacância de Ni em até 3 vezes e facilita a movimentação em três dimensões. Mecanismos de anel não são energeticamente eficientes em comparação com saltos diretos. As configurações estáveis de intersticiais em ambos compostos consistem em um átomo de Ni sobre o plano (001) entre dois vizinhos de Zr fora do plano. Intersticiais de Zr são instáveis e tendem a deslocar um átomo de Ni, ocupando seu sítio. Energias de deslocamento foram estudadas através de simulações de irradiação de ambos compostos. Durante o processo de colisão binária, um potencial universal ZBL foi usado para colisões a curta distância. Para distâncias intermediárias usamos um potencial de união arbitrário. Zr mostrou—se mais difícil de ser arrancado de seu sítio do que Ni. Encontramos valores de energia de deslocamento no intervalo de aproximadamente 29 eV até 546 eV. Alguns resultados experimentais são mostrados e apresentam boa concordância com os cálculos.

Respostas periódicas em sistemas lineares e fracamente não lineares não ressonantes e comportamento dinâmico de sistemas rotativos com o uso da base dinâmica

Este trabalho visa realizar o estudo do comportamento dinâmico de um eixo rotor flexível, modelado segundo a teoria de Euler-Bernoulli e caracterizar as respostas periódicas de sistemas LTI (sistemas lineares invariantes no tempo) e sistemas fracamente não lineares de ordem arbitrária. Para tanto, é utilizada a base dinâmica gerada pela resposta impulso ou solução fundamental. O comportamento dinâmico de um eixo rotor flexível foi discutido em termos da função de Green espacial e calculada de maneira não-modal. Foi realizado um estudo do problema de autovalor para o caso de um um eixo rotor biapoiado. As freqüências são obtidas e os modos escritos em termos da base dinâmica e da velocidade de rotação. As respostas periódicas de sistemas LTI, utilizadas nas aproximações com sistemas fracamente não lineares, são obtidas, independentemente da ordem do sistema, como um operador integral onde o núcleo é a função de Green T-periódica. Esta função é caracterizada em termos das propriedades de continuidade, periodicidade e salto da função de Green T-periódica, e da base dinâmica Simulações foram realizadas para sistemas concentrados, matriciais e escalares, com o objetivo de mostrar a validade da metodologia desenvolvida com as propriedades da função de Green T-periódica. Foi abordado um modelo não-linear para uma centrífuga utilizada na indústria textil [Starzinski, 1977].

Física de Pomeron a altas energias

Neste trabalho investigamos a física de Pomeron, tanto em sua formulação perturbativa dada pela Cromodinâmica Quântica como o aspecto não-perturbativo dado pela teoria de Regge. Estes formalismos teóricos levam em consideração o fato que o Pomeron porta os números quânticos do vácuo e deve dominar o espalhamento entre partículas em altas energias. Partindo do formalismo de dipolos de cor, estudamos a interação perturbativa usual corrigida por efeitos de unitariedade em altas energias para os processos DIS e Drell- Yan, determinando a importância destas correções para ambos processos. Referente ao Pomeron não-perturbativo, investigamos os aspectos do modelo de múltiplos espalhamentos, o qual também contribui também para a unitariedade no setor suave. Dado os resultados similares entre modelos QCD e de Regge na descrição de DIS difrativo, propomos um novo observável, o coeficiente angular difrativo da função de estrutura difrativa, que permitirá discriminar dinâmicas subjacentes ao processo. Finalmente estudamos a fenomenologia associada à truncagem da série assintótica advinda do formalismo BFKL, e a utilizamos para descrever os principais observáveis em colisões hadrônicas e elétron-próton.

Efeitos de parâmetros inerciais obtidos através de diferentes procedimentos na determinação de forças e torques articulares resultantes

O objetivo principal deste trabalho é avaliar os efeitos dos parâmetros inerciais (massa, localização do centro de massa e momento de inércia) obtidos através de diferentes procedimentos, no valor calculado para as forças e torques articulares resultantes, determinados através da dinâmica inversa. Para tal, pretende-se: (a) implementar um método para calcular a força nas articulações do tornozelo, joelho e quadril, em atividades motoras humanas consideradas bidimensionais, utilizando a técnica da dinâmica inversa, com os equipamentos disponíveis do Laboratório de Pesquisa do Exercício da Escola de Educação Física da Universidade Federal do Rio Grande do Sul; (b) desenvolver um protocolo de medição dos parâmetros inerciais dos segmentos corporais (massa, localização do centro de massa e momento de inércia) do membro inferior, baseado na pesagem hidrostática, e que apresente informações individualizadas; e (c) operacionalizar um protocolo de medição dos parâmetros inerciais dos segmentos corporais (massa, localização do centro de massa e momento de inércia) do membro inferior baseado no uso da tomografia computadorizada. Diversas situações consideradas bidemensionais foram avaliadas, como caminhada, corrida, agachamento, salto com amortecimento, e salto sem amortecimento. Para avaliar o modelo, os dados provenientes do cálculo foram confrontados com valores de força obtidos a partir da instrumentação de uma prótese de joelho (Cervieri, 2000). Os parâmetros inerciais obtidos através da tomografia computadorizada e da pesagem hidrostática apresentaram diferenças superiores a 100%, quando comparados com os valores fornecidos pelas tabelas antropométricas. Entretanto, no que se refere a fase de contato com o solo dos eventos realizados (instante de maiores forças envolvidas), os diversos métodos de obtenção dos parâmetros inerciais, não apresentam diferenças significativas no valor máximo calculado para as forças internas, através da dinâmica inversa.

Tópicos em programação dinâmica

Muitos problemas de Dinâmica em Economia se encaixam dentro de uma estrutura de modelos de decisão seqüencial, sendo resolvidos recursivamente. Programação Dinâmica uma técnica de otimização condicionada que se encarrega de solucionar problemas desse tipo. Esse trabalho tem como objetivo apresentar uma resenha dos principais resultados teóricos em Programação Dinâmica. Os métodos da Programação Dinâmica são válidos tanto para problemas determinísticos como para os que incorporam variável incerteza. esperada objetividade de uma dissertação de Mestrado, no entanto, nos impediu de extender análise, deixando assim de considerar explicitamente neste trabalho modelos estocásticos, que teria enriquecido bastante parte destinada aplicações Teor ia Econômica. No capítulo desenvolvemos instrumental matemático, introduzindo uma série de conceitos resultados sobre os quais se constrói análise nos capítulos subsequentes. Ilustramos tais conceitos com exemplos que seguem um certo encadeamento. Nas seções 1.1 1.2 apresentamos as idéias propriedades de espaços métricos espaços vetoriais. Na seção 1.3, prosseguimos com tópicos em análise funcional, introduzindo noção de norma de um vetor de espaços de Banach. seção 1.4 entra com idéia de contração, Teor ema do Ponto Fixo de Banach e o teor ema de Blackwell. O Teorema de Hahn-Banach, tanto na sua forma de extensão quanto na sua forma geométrica, preocupação na seção 1.5. Em particular, forma geométrica desse teorema seus corolários são importantes para análise conduzida no terceiro capítulo. Por fim, na seção 6, apresentamos Teorema do Máximo. Ao final deste capítulo, como também dos demais, procuramos sempre citar as fontes consultadas bem como extensões ou tratamentos alternativos ao contido no texto. No capítulo II apresentamos os resultados métodos da Programação Dinâmica em si seção 2.1 cuida da base da teoria, com Princípio da Otimal idade de Eellman e a derivação de um algoritmo de Programação Dinâmica. Na seção 2.2 mostramos que esse algoritmo converge para função valor ótima de um problema de horizonte infinito, sendo que esta última satisfaz chamada Equação de Bellman. seção seguinte se preocupa em fornecer caracterizaçBes para função valor mencionada acima, mostrando-se propriedades acerca de sua monotonicidade concavidade. seção 2.4 trata da questão da diferenciabi idade da função valor, que permite se obter alguns resultados de estática Cou dinâmica} comparativa partir da Equação de Bellman. Finalmente, na seção 2.5 apresentamos uma primeira aplicação Teoria Econômica, através de um modelo de crescimento econômico ótimo. No capítulo III introduzimos uma outra técnica de otimização Programação Convexa- mostramos dificuldade em se tentar estabelecer alguma relação de dominância entre Programação Dinâmica Programação Convexa. Na seção 3.2 "apresentamos os Teoremas de Separação, dos quais nos utilizamos na seção seguinte para demonstrar existência de Multiplicadores de Lagrange no problema geral da Programação Convexa. No final desta seção dizemos porque não podemos inferir que em espaços de dimensão infinita Programação Convexa não pode ser aplicada, ao contrário da Programação Dinâmica, que evidenciaria uma dominancia dessa última técnica nesses espaços. Finalmente, capítulo IV destinado uma aplicação imediata das técnicas desenvolvidas principalmente no segundo capítulo. Com auxílio dessas técnicas resolve-se um problema de maximização intertemporal, faz-se uma comparação dos resultados obtidos através de uma solução cooperativa de uma solução não-cooperativa.

Endividamento, custo do capital e valor da empresa: uma análise para o setor siderúrgico

Na utilização de variados tipos de financiamento, a empresa necessita, primariamente, da análise do custo do capital - custo básico para remuneração dos recursos próprios e de terceiros que compõem a estrutura de capital - o qual deve ser a taxa referencial para haver retorno do investimento e geração_ de valor, emergindo, por conseguinte, a discussão acerca da repercussão do nível de endividamento, e, conseqüentemente, do custo de capital sobre o valor da empresa. Nesse contexto, existem modelos que consideram o endividamento irrelevante no valor da empresa, como: análise pelo Lucro Operacional ou de Modiglianni e Miller sem Imposto de Renda. Em outro extremo, há abordagens que o têm como função direta para o valor, como a do Lucro Líquido e de Modigliani e Miller com imposto de renda. Entre essas posições antagônicas podem ser citadas a posição Tradicional e do Custo de Dificuldades Financeiras, que admitem a influência do endividamento após determinado nível de alavancagem. A Teoria da Informação Assimétrica, por sua vez, sustenta ser a sinalização enviada com mudanças na estrutura de capital fator de influência no valor da empresa. No modelo de avaliação por Opções, argumenta-se que, dentro de determinadas hipóteses, a riqueza dos acionistas é incrementada por uma alavancagem financeira mais elevada. A avaliação em situação de incerteza para a empresa é refletida no modelo que utiliza o Capital Asset Price Model, onde argumenta-se que um aumento no endividamento eleva o beta e, conseqüentemente, o retorno esperado. Por último, há a análise baseada em fluxos de caixa onde o valor presente da empresa é adicionado ao valor presente líquido dos efeitos intrínsecos ao financiamento. O cálculo do valor da empresa também pode ser feito através do fluxo de caixa do acionista descontado pela taxa de remuneração do capital próprio ou fluxo operacional descontado pelo custo médio ponderado do capital, subtraindo o valor do endividamento. Este trabalho se propôs a analisar, durante o período 1994-1998, a relação do endividamento com o custo médio ponderado do capital e com valor de mercado de oito empresas, componentes de amostra representativa do setor siderúrgico brasileiro, e também a estimar o valor do capital próprio, através do modelo de desconto do fluxo de caixa operacional pelo custo médio ponderado do capital, verificando a aderência dessa avaliação, através de comparação com o valor de mercado observado. Pela análise dos resultados observamos, uma relação inversa entre custo médio do capital e endividamento. No que se refere à avaliação da empresa, o modelo apresentado de desconto pelo custo médio ponderado do capital, revelou-se de modo geral, satisfatório, com o valor estimado explicando 68% do valor de mercado, na análise da regressão; sendo que as maiores distorções foram verificadas nas empresas que tiveram grau de alavancagem elevado.

Pico de força das articulações do membro inferior e ativação muscular da coluna dorso-lombar durante o manuseio de carga com estilo livre

o objetivo do presente estudo foi analisar o pico de força articular (FA) do tornozelo, joelho e quadril do membro inferior direito, e a ativação eletromiográfia (AE) da musculatura da coluna dorso-lombar no manuseio de carga, no movimento de colocar e retirar uma caixa situada no chão e na altura dos olhos. Foram utilizadas cargas correspondentes a I% e 10% da massa corporal do sujeito, totalizando-se oito tarefas. A amostra foi composta por 4 mulheres e 4 homens, na faixa etária entre 23 a 36 anos. Nenhuma instrução foi dada em relação à forma de execução da tarefa caracterizando-a como estilo livre. Para a análise do pico de FA utilizou-se a dinâmica inversa 2D, com auxílio de rotinas desenvolvidas no software MATLAB. Para análise da AE foi utilizada a eletromiografia de superfície. Adotou-se um Índice Postural (IP) nas tarefas de altura baixa para quantificar a postura no início da retirada da carga, e no final da colocação da mesma. Para o tratamento estatístico foi realizado o teste de normalidade de Shapiro-Wlik, e de homogeneidade de Levene após a ANOVA e post-hoc Tukey-b. O nível de significância foi p< 0.05. Em todas as situações o pico de FA apresentou um padrão em que os valores foram maiores no tornozelo, seguidos pelo joelho e depois pelo quadril. O fator que influenciou significativamente o pico de FA, o valor RMS e pico de AE foi a altura, sendo os maiores valores encontrados na altura baixa com peso pesado no movimento de retirar. Em uma relação temporal, os picos de FA e picos de AE na tarefa da altura baixa no movimento de retirar, em sua maioria, ocorreram na primeira metade do movimento, e no movimento de colocar na segunda metade. Já em uma relação temporal na tarefa da altura alta, não houve um padrão. O IP identificou que a amostra realizou uma postura predominantemente de agachamento, não havendo diferença intra e entre os sujeitos. Durante o manuseio de carga na altura baixa, o instante do pico de FA ocorreu quando o joelho e o quadril estavam em flexão máxima, e o tornozelo em desiflexão. Infere-se, então, que as articulações estudadas podem estar mais suscetíveis a lesões nestes espaços temporais ao se manusearem cargas que se encontram ao nível do chão.

Modelos cinéticos da equação linearizada de Boltzmann e um problema de transferência de calor em microescala

Neste trabalho, um problema de transferência de calor da dinâmica de gases rarefeitos, causado pela diferença de temperaturas nas superfícies de um canal, é abordado. O problema é formulado através dos modelos cinéticos BGK, S e Gross-Jackson da equação linearizada de Boltzmann e resolvido, de forma unificada, pelo método analítico de ordenadas discretas (método ADO). Resultados numéricos para as perturbações de densidade e temperatura e também para o fluxo de calor são apresentados e comparados, mostrando que não se pode dizer que algum dos três modelos seja uma melhor aproximação da solução aos resultados da equação linearizada de Boltzmann.

Estudo da passividade da Liga 600 em meio ácido

A passividade da Liga 600 (76Ni 16Cr 8Fe), em Na2SO4 0,5 M, pH=2,0, em atmosfera desarejada e à temperatura ambiente, foi estudada empregando-se diferentes métodos eletroquímicos e não eletroquímicos. A voltametria cíclica, com eletrodo rotatório de disco, revelou um comportamento típico ativo-passivo, com valores para a densidade de corrente anódica bastante baixos, na ordem de alguns poucos mA/cm2. Dois picos de corrente anódica pouco resolvidos foram observados e atribuídos à provável dissolução ativa de níquel, cromo e ferro. A ausência de picos catódicos e a existência de uma histerese na região de potenciais negativos sugerem que o filme passivo formado na varredura direta não é totalmente reduzido na varredura inversa, permanecendo sempre algum tipo de filme residual sobre a superfície da liga. A região passiva se estende de aproximadamente 100 a 700 mV e corresponde à região onde níquel e cromo puros também encontram-se passivos nas condições experimentais empregadas. Na região de potenciais mais positivos do que 700 mV tem início o processo de dissolução transpassiva da liga. Constatou-se, também, que o comportamento ativo-passivo da liga é essencialmente influenciado pelo comportamento do cromo, o qual é conhecido ser bastante complexo. Através das medidas de impedância eletroquímica foi possível sugerir três circuitos equivalentes para o sistema liga/filme/solução, um para cada região de potencial (de dissolução ativa, passiva e transpassiva). Através dos mesmos pôdese caracterizar a composição química e transformações mais importantes apresentadas pelos filmes passivos formados sobre a Liga 600. As espectroscopias eletrônicas (Auger e XPS) revelaram que os filmes passivos formados são extremamente finos, na faixa de 1,2 a 1,8 nm, e que apresentam uma estrutura duplex, com uma região interna (em contato com a liga) enriquecida em cromo e uma região externa (em contato com a solução) rica em níquel e ferro. Além disso, com base nos resultados obtidos e no modelo previamente proposto para filmes passivos formados sobre o aço inoxidável 304 em solução de borato, é sugerida uma representação esquemática das prováveis estruturas dos óxidos e dos possíveis processos de transporte, para os filmes passivos formados sobre a liga. O comportamento capacitivo dos filmes passivos foi estudado empregando-se a equação de Mott-Schottky. Os resultados obtidos mostram que os filmes formados se comportam como semicondutores degenerados do tipo n e do tipo p, na região de potenciais situada maiores e menores do que o potencial de banda plana, respectivamente. Esse comportamento é considerado conseqüência das propriedades semicondutoras dos óxidos de ferro (tipo n) e cromo (tipo p) presentes nos filmes passivos. Essa interpretação é fortalecida pelos resultados obtidos através das espectroscopias eletrônicas, as quais possibilitam o estabelecimento de uma relação direta entre a composição química das duas regiões de óxidos e a análise de Mott-Schottky. O comportamento dos filmes formados na região de potenciais próximos ao potencial de banda plana é essencialmente controlado pelo óxido de níquel, cujo comportamento pode ser comparado ao de um dielétrico, sem alterar a semicondutividade do óxido de ferro, quando ambos encontram-se misturados. O alto grau de degenerescência se deve ao valor elevado da densidade de doadores e aceptores (~ 1021 cm-3). Baseado nos resultados obtidos, o modelo da estrutura eletrônica previamente proposto para explicar a semicondutividade de filmes passivos e térmicos crescidos sobre o aço inoxidável 304, pode também ser aplicado no presente estudo. Segundo tal modelo, a estrutura eletrônica dos filmes formados pode ser comparada a de uma heterojunção do tipo p–n, onde as regiões de carga espacial encontram-se localizadas nas interfaces liga-filme e filme-solução.

Esterificação direta de butenos com ácido acético em fase líquida

Este trabalho investiga o comportamento da reação de esterificação direta de butenos com ácido acético em fase líquida a partir de uma carga contendo butenos e butano. O estudo faz parte de um projeto que visa desenvolver uma rota alternativa para a produção do acetato de sec-butila, um solvente amplamente utilizado na indústria de tintas e vernizes e de elevado valor agregado. Foram testados três tipos de catalisadores para a reação de interesse: zeólitas, resinas de troca iônica e um catalisador a base de óxido de nióbio. As resinas mostraram-se mais eficientes sendo que, entre as resinas, aquela que apresentou maior conversão a acetato foi a resina Amberlyst 36 Dry. O estudo avaliou também a desativação do catalisador num tempo de aproximadamente 300h de operação a 90 ºC. Observou-se que, além do acetato de sec-butila, os principais sub-produtos formados foram o álcool sec-butílico, o éter butílico e oligômeros provenientes do buteno. Foram realizados experimentos a 50 ºC, 70 ºC, 90 ºC e 110 ºC com a razão molar inicial de ácido/buteno variando entre 0,6 e 2,6. Os testes foram conduzidos em um reator batelada com volume útil de um litro acoplado a um cromatógrafo, na ausência de efeitos difusivos. Foi proposto um modelo pseudo-homogêneo para descrever o comportamento da reação e foram estimadas as energias de ativação e os coeficientes pré-exponenciais das reações direta e inversa. A constante de equilíbrio termodinâmico também foi calculada, utilizando o modelo UNIFAC para a determinação dos coeficientes de atividade. Finalmente, os dados de equilíbrio foram comparados com os dados cinéticos, observando-se boa correlação principalmente para as temperaturas de trabalho mais altas.

Interação heparina-antitrombina : reconhecimento molecular caracterizado por ferramentas de modelagem molecular

A heparina foi isolada no início do século XX e permanece, até os dias atuais, como um dos mais importantes e eficientes agentes terapêuticos de ação antitrombótica. Sua atividade anticoagulante deve-se à ativação da antitrombina (AT), uma serpina responsável pela inibição fisiológica de serino-proteinases plasmáticas, tais como fIIa e fXa. Os esforços no sentido da elucidação dos aspectos estruturais e dinâmicos associados ao reconhecimento molecular da heparina pela AT, em nível atômico, vêm encontrando diversas dificuldades, principalmente associadas aos compostos sacarídicos. Em decorrência de sua elevada polaridade e flexibilidade, glicosaminoglicanos como a heparina são difíceis de estudar e modelar. Soma-se a isto o fato de que os resíduos de iduronato presentes na heparina (IdoA) apresentam um incomum equilíbrio conformacional entre estados de cadeira (1C4) e bote-torcido (2SO), sendo esta última estrutura postulada como a possível conformação bioativa. Sendo assim, este trabalho apresenta um estudo de modelagem molecular do perfil conformacional da heparina, tanto em solução quanto complexada à AT, utilizando cálculos ab initio e simulações de dinâmica molecular (DM) Em decorrência da ausência de parâmetros capazes de descrever polissacarídeos nos campos de força atualmente disponíveis, cargas atômicas foram geradas, utilizando-se os esquemas de Mulliken, Löwdin e cargas ajustadas ao Potencial Eletrostático, através de cálculos quantum-mecânicos na base 6-31G** e testadas na simulação de DM da heparina. Diversas condições de simulação, tais como modelos de água, concentrações de sais e as conformações 1C4 e 2SO do IdoA, foram avaliadas de forma a identificar as melhores condições para a descrição conformacional da heparina em solução. O protocolo de DM obtido foi então utilizado no estudo do complexo AT-heparina. Os resultados obtidos estão de acordo com os dados experimentais atualmente disponíveis acerca da conformação da heparina e da contribuição energética de cada resíduo de aminoácido da AT para a formação do complexo com a heparina. A partir dos cálculos ab initio realizados, foi proposto um refinamento na estrutura da heparina determinada por RMN, enquanto que as simulações de DM permitiram reinterpretar aspectos da estrutura tridimensional da AT determinada por cristalografia de raios-X. Adicionalmente, os dados obtidos sugerem que não há requerimento conformacional para a interação do IdoA com a AT Globalmente, os dados indicam que simulações de DM podem ser utilizadas para representar adequadamente a conformação da heparina, assim como para caracterizar e quantificar suas interações com a AT. Assim sendo, propomos o uso de simulações de DM do complexo entre a AT e a heparina, ou entre a AT e compostos derivados da heparina, como uma ferramenta útil no processo de desenvolvimento de novos agentes antitrombóticos e anticoagulantes.

Aprimoramento estatístico da regionalização de vazões máximas e médias. Aplicação a bacias hidrográficas do Rio Grande do Sul e Santa Catarina

Os registros de vazões líquidas obtidos das estações fluviométricas são válidos, unicamente, para o local de onde foram coletados ou muito próximo a ele. Na maioria das vezes, a região de influência deste não inclui o local de interesse para o desenvolvimento de projetos de recursos hídricos. Este inconveniente, geralmente, pode ser resolvido através do uso de métodos de regionalização hidrológica. Para determinar os coeficientes da equação de regionalização, o procedimento mais usado consiste em: i) estabelecer uma relação do tipo exponencial multivariada entre a variável dependente (vazão média de longo prazo ou média anual de cheia) e as covariáveis (variáveis climáticas e fisiográficas da bacia hidrográfica); ii) linearizar a equação anterior mediante a transformação logarítmica de ambos os membros; iii) utilizar modelos lineares de regressão para estimar os coeficientes, geralmente, o método dos mínimos quadrados ordinários; e iv) aplicar a transformação inversa para definir a equação. A aplicação deste procedimento implica assumir certas propriedades dos dados (assimetria positiva, registros da mesma extensão e que os mesmos possuem o mesmo período de início e fim, entre outros) que dificilmente podem ser atendidas, prejudicando a verificação das hipóteses nas quais estão baseados os métodos lineares de regressão e, em conseqüência, seu adequado uso, bem como a confiabilidade dos resultados obtidos. Esta pesquisa apresenta um aprimoramento dos métodos de regionalização de vazões geralmente empregados, incluindo-se técnicas que levam em consideração as limitações anteriores. Estas técnicas foram: i) uso da transformada de Box-Cox na linearização da equação exponencial multivariada; ii) determinação dos coeficientes da equação de regionalização usando mínimos quadrados ponderados; e iii) verificação se os resíduos da regressão estão correlacionados ou não. Para o desenvolvimento e verificação da metodologia proposta foram usados somente registros fluviométricos de Bacias Hidrográficas Brasileiras, que drenam suas águas para o Rio Grande do Sul e/ou que estejam localizadas dentro dele. Geograficamente, a área de estudo inclui a totalidade do estado do Rio Grande do Sul e parte de Santa Catarina. As equações de regionalização foram definidas usando dados de vazões médias de longo prazo e média de cheia, para tempo de retorno de 2,33 e 50 anos. Neste último caso, as freqüências foram estimadas através do método dos momentos-L.Comparando os resultados obtidos utilizando o modelo de regionalização hidrológica proposto neste trabalho (transformada de Box-Cox / mínimos quadrados ponderados) junto a seus similares gerados usando a metodologia convencional (transformada logarítmica / mínimos quadrados ordinários) e de modelos intermediários (transformada logarítmica / mínimos quadrados ponderados e transformada de Box-Cox / mínimos quadrados ordinários), os mesmos podem ser considerados satisfatórios, visto que em todas as simulações realizadas o modelo proposto forneceu melhores resultados que aqueles obtidos com os outros modelos, sendo utilizado como padrão de comparação: 1) a qualidade do ajuste, 2) o grau de verificação das hipóteses dos métodos lineares de regressão e 3) os erros na estimativa das descargas, em termos de vazão específica. Nas simulações realizadas usando os modelos intermediários, observou-se que: i) na regionalização de vazões médias, o ganho de considerar a heterogeneidade temporal dos dados é maior do que corrigir a assimetria dos mesmos; ii) quando são usadas séries de descargas máximas, ocorre o efeito contrário, visto que o ganho de corrigir a assimetria das séries é maior do que o efeito da heterogeneidade temporal dos dados. Com relação aos resíduos da regressão, contrariamente ao esperado, os mesmos não sugerem estar correlacionados; isto pode ser conseqüência de utilizar como variável dependente um único registro por estação (vazão média de longo prazo ou média anual de cheia).

Solução da equação de difusão unidimensional transiente para o estudo da dispersão de poluentes na camada limite planetária

Neste trabalho apresenta-se uma solução analítica para a dispersão vertical turbulenta em uma Camada Limite Convectiva e em uma Camada Limite Estável. A equação analisada considera a difusão com velocidades finitas, o que representa o transporte turbulento fisicamente correto. Considerando o caráter não-local, adicionam-se na equação que representa uma fonte área instantânea, termos como: o tempo de relaxação, a assimetria, a escala de tempo Lagrangeana e a velocidade turbulenta vertical. A solução é obtida utilizando-se a técnica da Transformada de Laplace. Os parâmetros que encerram a turbulência são derivados da teoria de difusão estatística de Taylor combinada com a teoria de similaridade. Foram utilizados coeficientes de difusão especáficos para cada uma das camadas. A transformada inversa é obtida através do esquema numérico de quadratura Gaussiana. São apresentadas várias simulações para diferentes alturas de fonte área e obtém-se o valor da concentração para alturas próximas ao solo e próximas ao topo da Camada Limite Planetária. A inserção do termo de contra-gradiente na equação resultou em uma pequena influência na concentração de poluentes, observada de forma mais expressiva na Camada Limite Convectiva.

Avaliação do ensino superior : uma visão geral e uma análise comparativa entre os contextos brasileiro e português

Reconhecendo que a busca da qualidade no ensino superior se tornou hoje uma preocupação maior dos sistemas educativos de muitos paises tentou-se, neste estudo, compreender este fenômeno, particularmente na forma em que se desenvolveu como componente dos sistemas de avaliação do ensino superior de Brasil e Portugal. Tomando como ponto de partida realidades que manifestam diferentes aproximações e interpretações do tema em discussão, formularam-se hipóteses de trabalho para considerar o que poderá ter influenciado as discussões e negociações que ocorreram para implantação dos sistemas de avaliação do ensino superior de ambos os países: por exemplo, o efeito de novas formas de regulação dos sistemas de ensino superior; o grau de autonomia das instituições públicas de ensino superior; o comportamento e o relacionamento dos Conselhos de Reitores com os governos respectivos; e o funcionamento do mercado e das zonas de influência econômica (no caso do Brasil, o Mercosul e, no caso de Portugal, a União Européia). Tendo em conta estas hipóteses, foram analisadas as atas dos Conselhos de Reitores em ambos os países e feitas entrevistas com diversos atores com o objetivo de compreender de forma comparada os processos que decorreram em ambos os países em estudo Com base nesta análise foi possível concluir que, em ambos os países, os sistemas de regulamentação do ensino superior influenciaram a implantação dos sistemas de avaliação. Contudo, no Brasil, esta influência foi caracterizada pelo modelo de controle do Estado, enquanto que em Portugal foi o modelo de supervisão estatal que predominou. No que respeita à influência do grau de autonomia das instituições de ensino superior, existem diferenças importantes entre os dois países: no Brasil, por exemplo, há falta de autonomia enquanto em Portugal a autonomia tem uma expressão importante. Em relação aos Conselhos de Reitores o seu papel foi importante: no Brasil pela sua ausência e em Portugal por meio de uma participação ativa e decisiva. Em relação ao mercado, foi considerado que tinha alguma importância na universidade mas que a sua influência não se fez sentir diretamente no estabelecimento dos sistemas de avaliação do ensino superior. Quanto à influência das zonas econômicas ficou claro que o Mercosul não teve participação na implementação do sistema de avaliação do ensino superior brasileiro. Contudo, a União Européia foi uma das razões que levaram o Conselho de Reitores Português a iniciar o processo de avaliação nas universidades. Além de tudo isto, verificou-se que a lógica por trás do desenvolvimento da autonomia e da avaliação foi inversa nos dois países: no Brasil o ponto de partida foi a avaliação vista como conduzindo à autonomia; em Portugal, a autonomia foi o ponto de partida para o desenvolvimento da avaliação.

Solução LTSn para problemas de transferência radiativa com polarização em geometria plana

O método LTSN tem sido utilizado na resolução de uma classe abrangente de problemas de transporte de partículas neutras que são reduzidos a um sistema linear algébrico depois da aplicação da transformada de Laplace. Na maioria dos casos estudados os autovalores associados são reais e simétricos. Para o problema de criticalidade os autovalores associados são reais ou imaginários puros e simétricos, e para o o problema de multigrupo podem aparecer autovalores complexos. O objetivo deste trabalho consiste na generalização da formulação LTSN para problemas de transporte com autovalores complexos. Por esse motivo é focada a solução de um problema radiativo de transporte com polarização em uma placa plana. A solução apresentada fundamenta-se na aplicação da transformada de Laplace ao conjunto de equações SN dos problemas resultantes da decomposição da equação de transferência radiativa com polarização em série de Fourier, seguindo o procedimento de Chandrasekhar. Esse procedimento gera 2L + 2 sistemas lineares de ordem 4N dependentes do parâmetro complexo "s". Aqui, L é o grau de anisotropia e N a ordem de quadratura. A solução desse sistema simbólico é obtida através da aplicação da transformada inversa de Laplace depois da inversão da matriz simbólica pelo método da diagonalização. Para a obtenção das constantes de integração é assumido que os componentes do vetor de Stokes são reais e as matrizes dos autovalores e autovetores são separadas em suas partes real e imaginária. A solução LTSN para autovalores complexos é validada através da comparação da solução para uma placa com espessura unitária, grau de anisotropia L = 13, albedo de espalhamento simples $ = 0:99, coe ciente de re exão de Lambert ¸0 = 0:1 e N = 150, segundo dados da literatura consultada.