Avaliação da concentração de Cl, K e Ca na urina, hemolinfa e túbulos de Malpighi de Rhodnius prolixus usando a técnica de fluorescência de raios X por reflexão total por radiação síncrotron (SR-TXRF)

Neste trabalho utilizou-se a técnica fluorescência de raios X usando radiação síncrotron (SR-TXRF) para estudar, quantitativamente, o transporte de cloro, potássio e cálcio na hemolinfa, urina e túbulos de Malpighi (TM) em ninfas de quinto estágio do Rhodnius prolixus (R. prolixus), considerando a excreção destes elementos em diferentes dias após o repasto sanguíneo. R. prolixus é um dos principais vetores do Trypanosoma cruzi, agente causador da doença de Chagas. R. prolixus fornece um sistema modelo particularmente útil porque seus TMs tanto secretam quanto reabsorvem íons a taxas elevadas. Os TMs filtram a hemolinfa e secretam um líquido que é muitas vezes comparado com a urina primária em vertebrados. Os resultados obtidos mostram que a concentração de potássio na urina é substancialmente maior do que na hemolinfa. A concentração de cloro na hemolinfa é menor do que na urina, mas a diferença não é tão marcada como no caso do potássio. No caso do Rhodnius é razoável interpretar a elevada concentração de potássio na urina como adaptativo para o problema de excreção imediato do inseto. A concentração de cálcio nos TMs é substancialmente maior em comparação com os valores encontrados na hemolinfa e urina. Este resultado mostra que o cálcio é retido no corpo do R. prolixus e pouco eliminado. Os resultados obtidos estão coerentes com a literatura. Avaliou-se também o efeito no transporte de Cl, K e Ca após um repasto de sangue de coelho contaminado com HgCl2 de modo a avaliar o efeito da presença deste metal tóxico no balanço iônico nos fluidos de excreção urina e hemolinfa e também pelo principal órgão de transporte, os túbulos de Malpighi. As excreções de Cl e K pela urina são afetadas pela ingestão. Este resultado é esperado levando-se em consideração a ingestão de excesso de Cl através do HgCl2. O transporte de Cl, K e Ca na hemolinfa do Rhodnius prolixus não é afetada pela ingestão de HgCl2. Nos túbulos de Malpighi, as altas concentrações de Ca obtidas foram comparáveis àquelas encontradas nos insetos controle. Pode-se concluir que SR-TXRF é um método muito promissor para análises diretas, rápidas e confiáveis para a quantificação simultânea de elementos envolvidos na regulação do transporte e em todo o sistema excretor de insetos. Além disso, o estudo do transporte e a excreção de elementos no inseto Rhodnius prolixus abrem oportunidade para a maior compreensão de efeitos da poluição em espécies de invertebrados.

Simulação de escoamentos incompressíveis empregando o método Smoothed Particle Hydrodynamics utilizando algoritmos iterativos na determinação do campo de pressões

Nesse trabalho, foi desenvolvido um simulador numérico (C/C++) para a resolução de escoamentos de fluidos newtonianos incompressíveis, baseado no método de partículas Lagrangiano, livre de malhas, Smoothed Particle Hydrodynamics (SPH). Tradicionalmente, duas estratégias são utilizadas na determinação do campo de pressões de forma a garantir-se a condição de incompressibilidade do fluido. A primeira delas é a formulação chamada Weak Compressible Smoothed Particle Hydrodynamics (WCSPH), onde uma equação de estado para um fluido quase-incompressível é utilizada na determinação do campo de pressões. A segunda, emprega o Método da Projeção e o campo de pressões é obtido mediante a resolução de uma equação de Poisson. No estudo aqui desenvolvido, propõe-se três métodos iterativos, baseados noMétodo da Projeção, para o cálculo do campo de pressões, Incompressible Smoothed Particle Hydrodynamics (ISPH). A fim de validar os métodos iterativos e o código computacional, foram simulados dois problemas unidimensionais: os escoamentos de Couette entre duas placas planas paralelas infinitas e de Poiseuille em um duto infinito e foram usadas condições de contorno do tipo periódicas e partículas fantasmas. Um problema bidimensional, o escoamento no interior de uma cavidade com a parede superior posta em movimento, também foi considerado. Na resolução deste problema foi utilizado o reposicionamento periódico de partículas e partículas fantasmas.

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.

Caracterização de juntas soldadas por TIG autógeno manual com arco pulsado e arco não pulsado do aço inoxidável hiperduplex SAF 2707 HD.

O aço inoxidável hiperduplex, SAF 2707 HD, foi desenvolvido com o intuito de se obter novas ligas com maior resistência à corrosão do que a disponível nos aços inoxidáveis duplex / superduplex. Além da melhorada resistência à corrosão, este tb oferece propriedades mecânicas superiores. Os aços hiperduplex são aços de última geração que possuem elevados teores de elementos de liga, principalmente cromo, molibdênio e nitrogênio. Este tipo de aço caracteriza-se por apresentar estrutura bifásica, constituída de proporções praticamente iguais de ferrita e austenita devido à distribuição controlada dos elementos alfagênicos e gamagênicos. O interesse por esses aços cresce gradativamente com a necessidade de novos materiais para diversas aplicações na indústria petrolífera. Porém, quando expostos e mantidos a temperaturas elevadas, na faixa entre 600C e 1000C, algumas fases intermetálicas podem se formar, em que a fase sigma () é a mais proeminente. Possui uma estrutura cristalina tetragonal complexa rica em Cr e Mo, tendo efeito deletério no material afetando tanto a resistência à corrosão, quanto as propriedades mecânicas. Para este fim, faz-se necessário estudos da junta soldada para delinear as limitações desses aços e aperfeiçoar a aplicação. Essa pesquisa teve como objetivo caracterizar a junta soldada por TIG autógeno manual com arco pulsado e não pulsado do aço inoxidável hiperduplex SAF 2707 HD. As técnicas empregadas foram a metalografia por ataque eletrolítico (reagente NaOH) e color etching (reagente Behara), medidas de microdureza e quantificação microestrutural por Processamento Digital de Imagem. Os aspectos microestruturais foram observados por microscopia óptica (MO) e microscopia eletrônica de varredura (MEV), estes passaram por etapas de processamento digital de imagens (PDI) para quantificar a fração volumétrica da fase austenita. Realizou-se análise química semi-quatitativa por EDS. Os resultados foram analisados estatisticamente através do teste de hipóteses com distribuição t de Student. Pela técnica color etching observou-se que a fase austenita foi gerada com distribuição mais homogênea para o arco pulsado, que o não pulsado. O ataque eletrolítico não revelou uma terceira fase (fase ) na junta soldada, a análise química por EDS não identificou uma variação significativa nos elementos presentes ao longo da zona de transição do metal de base para a zona de fusão. Através do PDI foram obtidos os valores médios da fração volumétrica de austenita de 36,38% (desvio padrão 6,40%) e 32,41% (desvio padrão 6,67%) para os dois métodos, pulsado e não pulsado, respectivamente. Foram obtidos os valores de microdureza para o metal de base 355,10 HV (desvio padrão 28,60) e para a zona de fusão 343,60 HV (desvio padrão 20,51) da amostra soldada pelo modo pulsado, para o modo não pulsado foram apresentaram os valores médios de 370,30 HV (desvio padrão 34,51) para o metal de base e de 345,20 HV (desvio padrão 41,33) para a zona de fusão. A análise estatística indicou que não houve variação significativa da fração volumétrica da fase austenita no cordão de solda para as duas condições testadas e não houve variação da microdureza entre a zona de fusão e o metal de solda das amostras submetidas aos dois processos.

Modelagem matemática da transferência de calor numa placa plana sob o efeito de uma fonte pontual externa de radiação térmica.

Este trabalho apresenta uma modelagem matemática para o processo de aquecimento de um corpo exposto a uma fonte pontual de radiação térmica. O resultado original que permite a solução exata de uma equação diferencial parcial não linear a partir de uma seqüência de problemas lineares também é apresentado. Gráficos gerados com resultados obtidos pelo método de diferenças finitas ilustram a solução do problema proposto.

Estudo do comportamento em fluência de juntas metálicas coladas

Juntas coladas têm se mostrado como a principal alternativa de união entre componentes fabricados em materiais distintos ou não, tendo como vantagem principal a ausência de aporte de calor e melhor distribuição de tensões. A literatura mostra que estudos anteriores foram realizados com o objetivo de caracterizar diversas propriedades associados às juntas, porém, pouco ainda se conhece sobre o comportamento dessas uniões relacionado à fluência. Os ensaios de fluência têm como principal desvantagem a baixa produtividade de resultados visto que, na maioria dos casos, os equipamentos aptos a realizarem este tipo de teste ensaiam apenas um corpo de prova por vez. Neste caso, o tratamento estatístico é baseado em uma pequena quantidade de resultados, reduzindo a confiabilidade na predição do comportamento em fluência. O Laboratório de Adesão e Aderência (LAA/UERJ/IPRJ) desenvolveu o equipamento pneumático de fluência (EPF) capaz de realizar dez ensaios simultâneos com parâmetros distintos. Este trabalho investiga o comportamento de juntas de cisalhamento simples (SLJ) compostas por substratos metálicos colados com adesivo epóxi e poliuretano. As juntas foram fabricadas em conformidade a norma ASTM D 1002, seguindo procedimento de preparo RT-003/08 TMC/CENPES/PETROBRAS. Os resultados mostram que, para projetos onde esteja prevista utilização das juntas durante longos períodos de tempo sob carregamento, conhecer apenas suas tensões médias de ruptura não é o suficiente para garantir a segurança do empreendimento. É proposto um modelo inicial do comportamento sob fluência de juntas coladas.

O fluxo de massa do Vieira, Teresópolis-RJ, megadesastre da serra em 2011: análise das feições sedimentológicas e de sua fenomenalogia

Em janeiro de 2011, o grande número de escorregamentos naturais e induzidos que afetaram a Região Serrana do Estado do Rio de Janeiro, durante o Megadesastre, são processos destruidores de grande impacto, mas com tempo de recorrência maior. Em particular a bacia do rio Vieira (noroeste de Teresópolis), afetada pelo evento, apresentou magnitude espetacular do desastre registrado, com os materiais mobilizados em fluxos gravitacionais de massa atingindo quilômetros de distâncias das suas escarpas principais, o que causou, horas depois, prejuízos significativos. A pesquisa da presente dissertação discutiu a caracterização destes processos e dos ambientes deposicionais, a partir da análise das feições sedimentológicas e apreciação das características mecânicas, do vale do rio Vieira para estabelecer base de conhecimento da evolução da paisagem local. O método de trabalho utilizou o reconhecimento dos condicionantes geológicos e a leitura do registro sedimentar. A corrida de massa do Vieira foi zoneada em diferentes estágios de comportamento baseado na classificação de fluxos gravitacionais de massa de Middleton e Hampton (1973). Esta classificação é adequada para linha de pesquisa, pois os autores definem fluxos gravitacionais de sedimentos como um termo geral para fluxos de sedimentos ou misturas sedimentos/fluidos, que fluem sob a ação da gravidade. Os resultados obtidos caracterizam o fenômeno do ponto de vista geológico, estratigráfico e geomorfológico. Uma corrida de massa afetou o canal do rio Vieira desde a sua nascente, na cota 1750m, até o fim da sua planície de inundação, na cota 900m; o movimento de massa com extensão de 7,5km, 30-100m de largura máxima e 4m de espessura média, causou diretamente 86 mortes e a destruição de casas e lavouras da zona rural do bairro Vieira, Teresópolis-RJ. O fluxo gravitacional do Vieira aparentemente contou com diferentes condicionantes e envolveu diferentes fases, passando pelo debrisflow e pelo mudflow, provavelmente em diferentes pulsos, e no seu estágio final para o waterflow.

Captação de recursos por meio de debêntures conversíveis em ações: uma opção para os investimentos no Pré-Sal

O tema a respeito do financiamento da exploração e produção de petróleo, gás natural e outros hidrocarbonetos fluidos nas áreas do Pré-Sal tem sido motivo de debate entre estudiosos e profissionais da área. No novo regime de exploração e produção previsto na Lei n 12.351/2010, a Petróleo Brasileiro S.A. (Petrobras) será a operadora de todos os blocos contratados, ou seja, responsável pela condução e execução, direta ou indireta, de todas as atividades de exploração, avaliação, desenvolvimento, produção e desativação das instalações de exploração e produção. Tal incumbência legal, por si só, denota o volume de investimentos que a Petrobras deverá realizar nos próximos anos, seja com recursos próprios ou de terceiros, para exploração dessas novas áreas descobertas. Ademais, as sociedades contratadas para empreender as operações exploratórias também necessitarão de recursos, uma vez que, junto com a Petrobras, assumirão os custos e os investimentos referentes às atividades de exploração, avaliação, desenvolvimento, fabricação e desativação das instalações de exploração e produção decorrentes dos contratos de partilha. As debêntures conversíveis em ações apresentam-se como uma alternativa viável para a captação de recursos financeiros, além de proporcionarem vantagens, se comparadas com outras formas de financiamento, como empréstimos junto a instituições financeiras e aumento de capital, uma vez que é a companhia emissora responsável por decidir, por exemplo, a data e prazo de vencimento das debêntures e a forma de remuneração dos subscritores de tais títulos. O novo regime legal das debêntures, instituído pela Lei n 12.431/2011, também veio a facilitar o uso desse instrumento para financiar as atividades de produção e exploração de petróleo, gás natural e outros hidrocarbonetos fluidos nas áreas do Pré-Sal. A Petrobras, mesmo enquanto sociedade de economia mista, poderá emitir debêntures conversíveis em ações, observadas algumas regras para a manutenção do seu controle pela União.

Uma dedução dos critérios de multicriticalidade e o cálculo de pontos tricríticos via otimização global

Uma dedução dos critérios de multicriticalidade para o cálculo de pontos críticos de qualquer ordem representa a formalização de ideias utilizadas para calcular pontos críticos e tricríticos e ainda amplia tais ideias. De posse desta dedução pode-se compreender os critérios de tricriticalidade e, com isso, através de uma abordagem via problema de otimização global pode-se fazer o cálculo de pontos tricríticos utilizando um método numérico adequado de otimização global. Para evitar um excesso de custo computacional com rotinas numéricas utilizou-se aproximações na forma de diferenças finitas dos termos que compõem a função objetivo. Para simular a relação P v - T optou-se pela equação de estado cúbica de Peng-Robinson e pela regra clássica de fluidos de van der Vaals, para modelagem do problema também se calculou os tensores de ordem 2, 3, 4 e 5 da função do teste de estabilidade. Os resultados obtidos foram comparados com dados experimentais e por resultados obtidos com outros autores que utilizaram métodos numéricos, equação de estado ou abordagem diferente das utilizadas neste trabalho.

Estudo teórico e experimental da velocidade do som em sistemas contendo hidrocarbonetos

A descoberta de reservatórios de petróleo na camada do pré-sal abriu novas perspectivas para a economia brasileira. A grandes profundidades, com elevado teor de CO2 e submetidos a altas pressões e temperaturas, a exploração de petróleo das camadas de pré-sal constitui um considerável desafio tecnológico. Neste contexto, o conhecimento dos coeficientes de compressibilidade são fundamentais para o estudo do comportamento termodinâmico e escoamento dos fluidos destes campos de produção. Com este objetivo, foram determinadas experimentalmente nesta dissertação a massa específica e a velocidade do som de sistemas envolvendo quatro hidrocarbonetos puros: n-decano, n-hexadecano, decalina e tetralina, assim como para sistemas multicompostos envolvendo estes hidrocarbonetos e para uma amostra de condensado de petróleo. Os resultados obtidos permitiram o cálculo dos coeficientes de compressibilidade isentrópica, isotérmica e isentálpica. É apresentada uma avaliação sobre a capacidade de se correlacionar dados experimentais de velocidade do som aos dados calculados pelas equações de estado de Peng-Robinson e Lee-Kesler com regras de mistura: (i) Spencer e Danner e (ii) Lee-Kesler

Simulação numérica do escoamento de água em áreas alagáveis da floresta amazônica com a utilização da estrutura de dados Autonomous Leaves Graph

A Amazônia exibe uma variedade de cenários que se complementam. Parte desse ecossistema sofre anualmente severas alterações em seu ciclo hidrológico, fazendo com que vastos trechos de floresta sejam inundados. Esse fenômeno, entretanto, é extremamente importante para a manutenção de ciclos naturais. Neste contexto, compreender a dinâmica das áreas alagáveis amazônicas é importante para antecipar o efeito de ações não sustentáveis. Sob esta motivação, este trabalho estuda um modelo de escoamento em áreas alagáveis amazônicas, baseado nas equações de Navier-Stokes, além de ferramentas que possam ser aplicadas ao modelo, favorecendo uma nova abordagem do problema. Para a discretização das equações é utilizado o Método dos Volumes Finitos, sendo o Método do Gradiente Conjugado a técnica escolhida para resolver os sistemas lineares associados. Como técnica de resolução numérica das equações, empregou-se o Método Marker and Cell, procedimento explícito para solução das equações de Navier-Stokes. Por fim, as técnicas são aplicadas a simulações preliminares utilizando a estrutura de dados Autonomous Leaves Graph, que tem recursos adaptativos para manipulação da malha que representa o domínio do problema

Estudio de las ecuaciones de Boussinesq para la modelización de tsunamis

[es]Podemos encontrar las ecuaciones de Boussinesq en la descripción de playas, rios y lagos. Estas ecuaciones estudian la dinámica de las aguas poco profundas como las ecuaciones “ Korteweg-deVries (KdV)". Sin embargo, a pesar de ser más conocidas, las ecuaciones de KdV, no son capaces de modelar olas solitarias propagándose en distintas direcciones. Entre muchas otras aplicaciones de las ecuaciones de Boussinesq destaca la de modelar olas de tsunamis. Estos tipos de olas ya son perfectamente descritos por las ecuaciones de Navier Stokes, pero todavía no existen técnicas que permitan resolverlas en un dominio tridimensional. Para ello se usan las ecuaciones de Boussinesq, pensadas como una simplificación de las ecuaciones de Navier Stokes. Los años 1871 y 1872 fueron muy importantes para el desarrollo de las ecuaciones de Boussinesq. Fue en 1871 cuando Valentin Joseph Boussinesq recibió el premio de la “Academy of Sciences”, por su trabajo dedicado a las aguas poco profundas. Ahí fue donde Boussinesq introdujo por primera vez los efectos dispersivos en las ecuaciones de Saint-Venant. Por ello, se puede decir que las ecuaciones de Boussinesq son más completas físicamente que las ecuaciones de Saint-Venant. Las ecuaciones de Boussinesq contienen una estructura hiperbólica (al igual que las ecuaciones no lineales de aguas poco profundas) combinada con derivadas de orden elevado para modelar la dispersión de la ola. Las ecuaciones de Boussinesq pueden aparecer de muchas formas distintas. Dependiendo de como hayamos escogido la variable de la velocidad podemos obtener un modelo u otro. El caso más usual es escoger la variable velocidad en un nivel del agua arbitrario. La efectividad de la ecuación de Boussinesq seleccionada variará dependiendo de la dispersión. Una buena elección de la variable velocidad puede mejorar significativamente la modelización de la propagación de ondas largas. Formalmente, como veremos en el capítulo 1, podemos transformar términos de orden elevado en términos de menor orden usando las relaciones asintóticas. Esto nos proporciona una forma elegante de mejorar las relaciones de dispersi\'on. Las ecuaciones de Boussinesq más conocidas son las que resolveremos en el capítulo 2. En dicho capítulo veremos la ecuación cúbica de Boussinesq, que sirve para describir el movimiento de ondas largas en aguas poco profundas; las ecuaciones de Boussinesq acopladas, que describen el movimiento de dos fluidos distintos en aguas poco profundas (como puede ser el caso de un barco que desprende accidentalmente aceite, el aceite va creando una capa que flota encima de la superficie del agua); la ecuación de Boussinesq estándar, que describe un gran número de fenómenos de olas dispersivas no lineales como la propagaci\ón en ambas direcciones de olas largas en la superficie de aguas poco profundas. Pero en olas de longitud de onda corta presenta una inestabilidad y la ecuación es incorrecta para el problema de Cauchy, por ello Bogolubsky propuso la ecuación de Boussinesq mejorada. Esta ecuación es la última que estudiaremos en el capítulo 2 y es una ecuación físicamente estable, correcta para el problema de Cauchy y además como veremos en el capítulo 3, apropiada para las simulaciones numéricas. Como ya indicado, en el capi tulo 1 deduciremos las ecuaciones de Boussinesq a partir de las ecuaciones físicas del flujo potencial. El objetivo principal es deducir dos modelos de ecuaciones de Boussinesq acopladas y obtener su relación de dispersión. Para llegar a ello, se usa un método de la expansión asintótica de la velocidad potencial en términos de un pequeño parámetro. De esta manera conseguimos dos modelos distintos, cada uno asociado a uno de los dos modelo de disipación que hemos establecido. Por último dado que las ecuaciones siempre vienen dadas en variables dimensionales, volveremos a la notación dimensional para analizar la relación de dispersión de las ecuaciones de Boussinesq disipativas. En el capí tulo 2 pasaremos a su resolución analítica, buscando soluciones de tipo solitón. Introduciremos el método de la tangente hiperbólica, muy útil para encontrar soluciones exactas de ecuaciones no lineales. Usaremos este método para resolver la ecuación cúbica de Boussinesq, un sistema de ecuaciones acopladas de Boussinesq, la ecuación estandar de Boussinesq y la mejorada. Los sistemas que aparecen en la aplicación del método de la tangente hiperbólica estan resueltos usando el software Mathematica y uno de ellos irá incluido en el apéndice A. En el capíulo 3 se introduce un esquema en diferencias finitas, que sirve para convertir problemas de ecuaciones diferenciales en problemas algebraicos fácilmente resolubles numéricamente. Este método nos ayudaráa estudiar la estabilidad y a resolver la ecuación mejorada de Boussinesq numéricamente en dos ejemplos distintos. En el apéndice B incluiremos el programa para la resolución numérica del primer ejemplo con el Mathematica.

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.

Uma formulação implícita para o método Smoothed Particle Hydrodynamics

Em uma grande gama de problemas físicos, governados por equações diferenciais, muitas vezes é de interesse obter-se soluções para o regime transiente e, portanto, deve-se empregar técnicas de integração temporal. Uma primeira possibilidade seria a de aplicar-se métodos explícitos, devido à sua simplicidade e eficiência computacional. Entretanto, esses métodos frequentemente são somente condicionalmente estáveis e estão sujeitos a severas restrições na escolha do passo no tempo. Para problemas advectivos, governados por equações hiperbólicas, esta restrição é conhecida como a condição de Courant-Friedrichs-Lewy (CFL). Quando temse a necessidade de obter soluções numéricas para grandes períodos de tempo, ou quando o custo computacional a cada passo é elevado, esta condição torna-se um empecilho. A fim de contornar esta restrição, métodos implícitos, que são geralmente incondicionalmente estáveis, são utilizados. Neste trabalho, foram aplicadas algumas formulações implícitas para a integração temporal no método Smoothed Particle Hydrodynamics (SPH) de modo a possibilitar o uso de maiores incrementos de tempo e uma forte estabilidade no processo de marcha temporal. Devido ao alto custo computacional exigido pela busca das partículas a cada passo no tempo, esta implementação só será viável se forem aplicados algoritmos eficientes para o tipo de estrutura matricial considerada, tais como os métodos do subespaço de Krylov. Portanto, fez-se um estudo para a escolha apropriada dos métodos que mais se adequavam a este problema, sendo os escolhidos os métodos Bi-Conjugate Gradient (BiCG), o Bi-Conjugate Gradient Stabilized (BiCGSTAB) e o Quasi-Minimal Residual (QMR). Alguns problemas testes foram utilizados a fim de validar as soluções numéricas obtidas com a versão implícita do método SPH.

Aperfeiçoamento do algoritmo algébrico sequencial para a identificação de variações abruptas de impedância acústica via otimização

Neste trabalho são utilizados a técnica baseada na propagação de ondas acústicas e o método de otimização estocástica Luus-Jaakola (LJ) para solucionar o problema inverso relacionado à identificação de danos em barras. São apresentados o algoritmo algébrico sequencial (AAS) e o algoritmo algébrico sequencial aperfeiçoado (AASA) que modelam o problema direto de propagação de ondas acústicas em uma barra. O AASA consiste nas modificações introduzidas no AAS. O uso do AASA resolve com vantagens o problema de identificação de danos com variações abruptas de impedância. Neste trabalho são obtidos, usando-se o AAS-LJ e o AASA-LJ, os resultados de identificação de cinco cenários de danos. Três deles com perfil suave de impedância acústica generalizada e os outros dois abruptos. Além disso, com o objetivo de simular sinais reais de um experimento, foram introduzidos variados níveis de ruído. Os resultados alcançados mostram que o uso do AASA-LJ na resolução de problemas de identificação de danos em barras é bastante promissor, superando o AAS-LJ para perfis abruptos de impedância.