Avaliação temporal da poluição por antiincrustantes organoestânicos no litoral do Estado do Rio de Janeiro: antes e após o banimento nacional e internacional

Tintas antifoulings são utilizadas para evitar a incrustação de organismos em estruturas submersas, especialmente casco de embarcações. Os compostos organoestânicos (OTs), utilizados nessas tintas, entre eles o tributilestanho, são desreguladores endócrinos e causaram diversos danos aos ecossistemas marinhos. No caso dos moluscos gastrópodes, esse tipo de poluição faz com que as fêmeas adquiram características masculinas, como vaso deferente e pênis, fenômeno esse conhecido como imposex. A Organização Marítima Internacional (IMO) estabeleceu o banimento de tintas à base de COEs nas embarcações, em 2008. No Brasil, a NORMAM 23, proibiu o uso em 2007, contudo a Marinha já havia suspendido seu uso desde 2003. Entretanto, efeitos deletérios destes compostos ainda são detectados em vários países, inclusive ao longo do litoral brasileiro. Esse trabalho teve como objetivo principal fazer uma avaliação temporal (1997-2012) da poluição por organestânicos na costa do Estado do Rio de Janeiro utilizando como bioindicador a espécie Stramonita haemastoma. A área de estudo abrangeu cinco regiões: Paraty, Ilha Grande, Baia de Ilha Grande, Baia de Guanabara e Arraial do Cabo. Análises químicas de butilestânicos foram feitas em sedimentos superficiais de mangues de Paraty (Mangue do Estaleiro e mangue do Saco do Mamanguá) e Ilha Grande (Mangue do Aventureiro) sendo os dois últimos considerados áreas de referência. Em cada estação de biomonitoramento foram coletados 30 indivíduos sexualmente adultos da espécie S. haemastoma, através de mergulho livre em apnéia e analisados através do método não destrutivo proposto por nosso grupo de pesquisa. Em todas as cinco áreas analisadas ao longo da costa do Estado do Rio de Janeiro foram registradas estações com altos índices de imposex. Apesar de muitas estações apresentarem diminuição nos índices, na maioria ocorreu um aumento ou conservação alta na porcentagem de imposex depois do banimento. As concentrações médias de butilestânicos no mangue (S1), perto de fontes locais, foram 205,7 16,8 ng (Sn) g-1 de TBT, 16,4 1,3 ng (Sn) g-1 de DBT e 10,0 2,9 ng (Sn) g-1 de MBT. Nas áreas de referência: mangue do Saco do Mamanguá (S2) foram 16,0 0,8 ng (Sn) g-1 de TBT, 10,1 1,4 ng (Sn) g-1 de DBT e 10,1 2,2 ng (Sn) g-1 de MBT e mangue do Aventureiro (S3) com 18,1 4,2 ng (Sn) g-1 de TBT, 15,3 0,5 ng (Sn) g-1 de DBT e 10,2 1,5 ng (Sn). g-1 de MBT. As taxas de degradação foram de 01, 1,3 e 1,4 respectivamente indicando inputs recentes desses compostos. Os resultados deste estudo indicam que carbono orgânico dissolvido e particulado, bem como xenoestrógenos podem estar interferindo no desenvolvimento da síndrome, levando a subestimação de avaliação do imposex. As concentrações de butilestânicos, além do aumento ou continuidade alta na incidência de imposex após o banimento na maioria das estações indicam que, apesar da proibição do uso do TBT em tintas antiiincrustantes no Brasil, elas ainda estão sendo utilizadas de forma ilegal, especialmente em pequenos barcos. Esse estudo é fundamental para se propor medidas de mitigação e controle dos compostos organoestânicos, até mesmo dos novos antifoulings TBT- free, que também possuem efeitos prejudiciais ao ambiente. Além disso, as áreas onde se registrou altos índices de imposex serão essenciais para o monitoramento dos efeitos desses novos antifoulings.

Anatomia da madeira e dendrocronologia de espécies arbóreas do gênero Machaerium Pers. (Leguminosae-Papilionoideae)

Machaerium é um dos maiores gêneros arbóreos tropicais de leguminosas, com cerca de 130 espécies com distribuição predominantemente neotropical e centro de diversidade no Brasil, onde ocorrem cerca de 80 espécies. O gênero ocorre em todos os domínios fitogeográficos do país, porém a Mata Atlântica e a Floresta Amazônica possuem os maiores índices de riqueza e endemismo. As espécies do gênero estão classificadas em cinco seções infragenéricas, que se baseiam principalmente na forma e venação dos folíolos e na presença de estípulas espinescentes. Entretanto, esta classificação tem sido questionada por alguns autores, principalmente quando comparada com análises filogenéticas. Dessa forma, surge a necessidade de buscar outros caracteres que auxiliem na delimitação das espécies e que permitam uma reavaliação na classificação nfragenéricas, além de conhecer o potencial para estudos dendrocronológicos das espécies em um bioma tão rico e ameaçado como a Mata Atlântica. O presente trabalho visou estudar a anatomia do lenho de onze espécies arbóreas de Machaerium a fim de verificar a consistência das seções infragenéricas, fornecer caracteres diagnósticos para a delimitação das espécies e caracterizar as pontoações intervasculares ornamentadas, para verificar seu potencial diagnóstico no gênero em questão. Além disso, analisar a periodicidade de crescimento e a influência dos fatores climáticos no crescimento de Machaerium incorruptibile, espécie endêmica da Mata Atlântica. As amostras foram coletadas através de método não destrutivo e processadas seguindo os métodos usuais para anatomia do lenho, microscopia eletrônica de varredura e dendrocronologia. As espécies apresentaram as características anatômicas descritas para a família Leguminosae e para a subfamília Papilionoideae. A presença de faixas de parênquima não lignificado, fibras de paredes delgadas e raios irregularmente estratificados foram importantes na separação de Machaerium hirtum das outras dez espécies. As dez espécies restantes foram separadas entre si pelos dados quantitativos do lenho, principalmente diâmetro e frequência de vasos, e também pelos caracteres morfológicos das pontoações ornamentadas, como a projeção ou não das ornamentações na abertura da pontoação. A anatomia da madeira não correspondeu as seções infragenéricas tradicionalmente tratadas para o gênero. Para a dendrocronologia os dados foram analisados com uso do software ARSTAN e foi construída uma cronologia para M. incorruptibile. As cronologias foram analisadas juntamente com os dados de precipitação e temperatura, onde observou-se uma correlação significativa com a temperatura. A largura dos anéis de crescimento foi correlacionada positivamente com a temperatura média da primavera, época em que a temperatura se encontra mais amena e, negativamente com a temperatura média do verão, onde as temperaturas são mais altas. Também houve correlação negativa entre os eventos mais severos de El Niño e a largura dos anéis de crescimento, demonstrando o efeito deste fenômeno no crescimento da população arbórea, como encontrado em outras espécies tropicais.

Dinâmica de crescimento em Tabebuia rosea (Bertol.) Bertero ex A.DC. (Bignoniaceae): fenologia apical e cambial caulinar

O estudo da dinâmica de crescimento em espécies arbóreas permite melhor compreender a adaptação de uma espécie ao ambiente em que se desenvolve. Dados dessa natureza ainda são escassos no Brasil considerando-se a diversidade de espécies e fitofisionomias existentes. Esse tipo de estudo tem sido apontado como de grande relevância principalmente diante das situações atuais de desmatamento e mudanças climáticas, que exigem recuperação de áreas degradadas e prospectar como as espécies nativas sobreviveram diante deste ambiente de aceleradas mudanças. Nesta perspectiva, o objetivo deste trabalho foi investigar a dinâmica de crescimento de Tabebuia rosea, por meio do monitoramento mensal da fenologia apical e cambial caulinar e dos teores dos pigmentos fotossintetizantes: clorofilas a, b, totais e carotenoides. Os resultados obtidos foram correlacionados entre si e também com a sazonalidade das variáveis ambientais: temperatura, precipitação e fotoperíodo. No que se refere aos teores dos pigmentos fotossintetizantes, foi também avaliada as oscilações observadas nas folhas obtidas nas orientações geográficas: norte, sul, leste e oeste. O trabalho foi desenvolvido no Campus da Universidade do Estado do Rio de Janeiro, RJ, Brasil. A fenologia dos ápices caulinares foi acompanhada mensalmente a partir de observações qualitativas e quantitativas em 15 indivíduos. As amostras contendo a zona cambial foram obtidas por método não destrutivo e processadas segundo as técnicas usuais em histologia vegetal. A extração dos pigmentos fotossintéticos foi realizada em acetona 80%, com posterior centrifugação e análise ao espectrofotômetro em diferentes comprimentos de onda. Tabebuia rosea apresentou aumento na largura da zona cambial em número de camadas celulares e em micrômetros, diminuição do diâmetro radial da célula inicial fusiforme e aumento em micrômetros da camada de células em processo de alongamento e diferenciação do xilema secundário no mesmo período em que a copa das árvores estava ocupada por folhas adultas e com os teores mais elevados de clorofilas nas folhas, coincidindo com o período em que as variáveis ambientais apresentaram-se elevadas, corroborando os resultados já obtidos para outras espécies nativas da Mata Atlântica se desenvolvendo no estado do Rio de Janeiro. A maior concentração de carotenoides foi observada no período em que as variáveis ambientais apresentavam índices elevados e as plantas estavam mais expostas à radiação luminosa, possivelmente em função de um maior investimento em fotoproteção. Cabe destacar que o lado oeste da copa das árvores mostrou uma tendência em apresentar maiores concentrações dos pigmentos fotossintéticos estudados, resultado ainda não observado na literatura até o momento, o que pode ser uma estratégia da espécie em compensar níveis baixos de radiação solar com acúmulo de clorofilas.

Avaliação do método de imunofenotipagem por citometria de fluxo no diagnóstico das doenças displásicas (SMD, LMA-relacionada à mielodisplasia, LMMC e LMMJ) em adultos e crianças

As síndromes mielodisplásicas (SMD) se caracterizam por terem uma hematopoese displásica, citopenias e pelo risco de progressão para leucemia mielóide aguda. O diagnóstico baseia-se na clínica e nos achados citomorfológicos da medula óssea (MO) e citogenéticos. Na fase inicial ou quando a MO é hipocelular o diagnóstico é difícil e a citogenética frequentemente é normal. A imunofenotipagem (IMF) tem sido cada vez mais utilizada nos casos de SMD em adultos e pouco explorada na SMD pediátrica. Os nossos objetivos foram: estudar os casos de SMD e doenças correlatas (LMA relacionada à SMD: LMA-rMD; leucemia mielomonocítica crônica: LMMC e leucemia mielomonocítica juvenil: LMMJ) em adultos e crianças, associando os dados clínicos e laboratoriais aos obtidos pela IMF, que utilizou um painel de anticorpos monoclonais para as várias linhagens medulares. No período compreendido entre 2000 e 2010 foram estudados 87 pacientes (64 adultos e 23crianças) oriundos do HUPE/UERJ e IPPMG/UFRJ e 46 controles (23 adultos e 20 crianças). Todos os doentes realizaram mielograma, biópsia óssea, citogenética, citoquímica e estudo imunofenotípico. Segundo os critérios da OMS 50 adultos foram classificados como SMD, 11 como LMA-rMD e 3 LMMC. Entre as crianças 18 eram SMD, 2 LMA e 3 LMMJ. Os pacientes adultos com SMD foram divididos em alto risco (n = 9; AREB-1 e AREB-2) e baixo risco (n=41; CRDU, CRDM, CRDM-SA, SMD-N e SMD-5q-). As crianças com SMD em CR (n=16) e AREB (n = 2). Anormalidades clonais recorrentes foram encontradas em 22 pacientes adultos e em 7 crianças. Na análise da IMF foi utilizada a metodologia da curva ROC para a determinação dos valores de ponto de corte a fim de identificar os resultados anormais dos anticorpos monoclonais nos pacientes e nos controles, permitindo determinar a sensibilidade e especificidade desses em cada linhagem. A IMF foi adequada para a análise em todos os pacientes e 3 ou mais anormalidades foram encontradas. A associação da IMF aumentou a sensibilidade da análise morfológica na linhagem eritróide de 70 para 97% nos adultos e de 59 para 86% nas crianças; na linhagem granulocítica de 53 para 98% nos adultos e de 50 para 100% nas crianças. Nos monócitos, onde a morfologia não foi informativa, mostrou uma sensibilidade de 86% nos adultos e 91% nas crianças. Enquanto que na linhagem megacariocítica, não analisada pela IMF, a morfologia mostrou uma sensibilidade de 95% nos adultos e 91% nas crianças. Na população de blastos foi expressiva a ausência de precursores linfóide B (em 92% dos adultos e em 61% das crianças). Os resultados observados nas crianças com SMD foram semelhantes aos encontrados nos adultos. Em conclusão, nossos resultados mostraram que a IMF é um método complementar ao diagnóstico da SMD e doenças correlatas tanto em adultos quanto em crianças podendo contribuir para o reconhecimento rápido e precoce dessas enfermidades, devendo ser incorporado aos procedimentos de rotina diagnóstica.

Análise da dissipação das tensões em dentes humanos restaurados com facetas laminadas de cerâmica, com três tipos de preparos, através do método dos elementos finitos

O objetivo deste trabalho é analisar in vitro a dissipação de tensões em incisivos centrais superiores humanos restaurados com facetas de cerâmica feldspática, através da análise do método dos elementos finitos, considerando cargas funcionais de mastigação e corte dos alimentos, em função de três tipos de preparos utilizados: sem proteção incisal; com proteção incisal em ângulo e com proteção incisal em degrau palatino. Foram utilizadas modelagens bidimensionais de um incisivo central superior e suas estruturas de suporte, simulando três situações: (Primeira modelagem) incisivo central superior com desgaste vestibular (em forma de janela); (Segunda modelagem) incisivo central superior com desgaste vestibular e proteção incisal em plano inclinado; (Terceira modelagem) incisivo central superior com desgaste vestibular, e proteção incisal com degrau palatino. Foi considerada uma carga (P=100N) com uma inclinação de 45 concentrada, simulando a região de contato do incisivo central inferior com o superior durante a mastigação e uma na região de contato topo a topo dos incisivos superior e inferior, simulando o corte dos alimentos. Após a análise dos dados obtidos pela distribuição de tensões, pode-se concluir que quanto à dissipação das tensões em todo o sistema proposto, com a aplicação de carga em 45, não foram observadas mudanças no estado tensional nos três diferentes preparos. Quando foi aplicada carga vertical, simulando o contato de topo, houve variação no estado tensional no sistema do dente com preparo em janela. Nas facetas, com a aplicação de carga em 45, nos preparos em janela e com proteção incisal em plano inclinado o resultado foi semelhante nos valores tensionais enquanto, nas facetas em dentes preparados com proteção incisal com degrau palatino, a distribuição foi mais homogênea tendo valores superiores, mostrando que o abraçamento do dente diminuiu a flexão.

Estimativa do erro de discretização analítico na solução de equações diferenciais utilizando o Método de Volumes Finitos

As análises de erros são conduzidas antes de qualquer projeto a ser desenvolvido. A necessidade do conhecimento do comportamento do erro numérico em malhas estruturadas e não-estruturadas surge com o aumento do uso destas malhas nos métodos de discretização. Desta forma, o objetivo deste trabalho foi criar uma metodologia para analisar os erros de discretização gerados através do truncamento na Série de Taylor, aplicados às equações de Poisson e de Advecção-Difusão estacionárias uni e bidimensionais, utilizando-se o Método de Volumes Finitos em malhas do tipo Voronoi. A escolha dessas equações se dá devido a sua grande utilização em testes de novos modelos matemáticos e função de interpolação. Foram usados os esquemas Central Difference Scheme (CDS) e Upwind Difference Scheme(UDS) nos termos advectivos. Verificou-se a influência do tipo de condição de contorno e a posição do ponto gerador do volume na solução numérica. Os resultados analíticos foram confrontados com resultados experimentais para dois tipos de malhas de Voronoi, uma malha cartesiana e outra triangular comprovando a influência da forma do volume finito na solução numérica obtida. Foi percebido no estudo que a discretização usando o esquema CDS tem erros menores do que a discretização usando o esquema UDS conforme literatura. Também se percebe a diferença nos erros em volumes vizinhos nas malhas triangulares o que faz com que não se tenha uma uniformidade nos gráficos dos erros estudados. Percebeu-se que as malhas cartesianas com nó no centróide do volume tem menor erro de discretização do que malhas triangulares. Mas o uso deste tipo de malha depende da geometria do problema estudado

Um método SN híbrido direto para cálculos de sistemas combustível-moderador em geometria unidimensional

Desenvolvemos nesta dissertação um método híbrido direto para o cálculo do fator de desvantagem e descrição da distribuição do fluxo de nêutrons em sistemas combustível-moderador. Na modelagem matemática, utilizamos a equação de transporte de Boltzmann independente do tempo, considerando espalhamento linearmente anisotrópico no modelo monoenergético e espalhamento isotrópico no modelo multigrupo, na formulação de ordenadas discretas (SN), em geometria unidimensional. Desenvolvemos nesta dissertação um método híbrido direto para o cálculo do fator de desvantagem e descrição da distribuição do fluxo de nêutrons em sistemas combustível-moderador. Na modelagem matemática, utilizamos a equação de transporte de Boltzmann independente do tempo, considerando espalhamento linearmente anisotrópico no modelo monoenergético e espalhamento isotrópico no modelo multigrupo, na formulação de ordenadas discretas (SN), em geometria unidimensional. Descrevemos uma análise espectral das equações de ordenadas discretas (SN)a um grupo e a dois grupos de energia, onde seguimos uma analogia com o método de Case. Utilizamos, neste método, quadraturas angulares diferentes no combustível (NC) e no moderador (NM), onde em geral assumimos que NC > NM . Condições de continuidade especiais que acoplam os fluxos angulares que emergem do combustível (moderador) e incidem no moderador (combustível), foram utilizadas com base na equivalência entre as equações SN e PN-1, o que caracteriza a propriedade híbrida do modelo proposto. Sendo um método híbrido direto, utilizamos as NC + NM equações lineares e algébricas constituídas pelas (NC + NM)/2 condições de contorno reflexivas e (NC + NM)/2 condições de continuidade para determinarmos as NC + NM constantes. Com essas constantes podemos calcular os valores dos fluxos angulares e dos fluxos escalares em qualquer ponto do domínio. Apresentamos resultados numéricos para ilustrar a eficiência e a precisão do método proposto.

Análise global da estabilidade termodinâmica de misturas: um estudo com o método do conjunto gerador

O cálculo do equilíbrio de fases é um problema de grande importância em processos da engenharia, como, por exemplo, na separação por destilação, em processos de extração e simulação da recuperação terciária de petróleo, entre outros. Mas para resolvê-lo é aconselhável que se estude a priori a estabilidade termodinâmica do sistema, a qual consiste em determinar se uma dada mistura se apresenta em uma ou mais fases. Tal problema pode ser abordado como um problema de otimização, conhecido como a minimização da função distância do plano tangente à energia livre de Gibbs molar, onde modelos termodinâmicos, de natureza não convexa e não linear, são utilizados para descrevê-lo. Esse fato tem motivado um grande interesse em técnicas de otimização robustas e eficientes para a resolução de problemas relacionados com a termodinâmica do equilíbrio de fases. Como tem sido ressaltado na literatura, para proporcionar uma completa predição do equilíbrio de fases, faz-se necessário não apenas a determinação do minimizador global da função objetivo do teste de estabilidade, mas também a obtenção de todos os seus pontos estacionários. Assim, o desenvolvimento de metodologias para essa tarefa desafiadora tem se tornado uma nova área de pesquisa da otimização global aplicada à termodinâmica do equilíbrio, com interesses comuns na engenharia química e na engenharia do petróleo. O foco do presente trabalho é uma nova metodologia para resolver o problema do teste de estabilidade. Para isso, usa-se o chamado método do conjunto gerador para realizar buscas do tipo local em uma rede de pontos previamente gerada por buscas globais efetuadas com uma metaheurística populacional, no caso o método do enxame de partículas.Para se obter mais de um ponto estacionário, minimizam-se funções de mérito polarizadas, cujos pólos são os pontos previamente encontrados. A metodologia proposta foi testada na análise de quatorze misturas polares previamente consideradas na literatura. Os resultados mostraram que o método proposto é robusto e eficiente a ponto de encontrar, além do minimizador global, todos os pontos estacionários apontados previamente na literatura, sendo também capaz de detectar, em duas misturas ternárias estudadas, pontos estacionários não obtidos pelo chamado método de análise intervalar, uma técnica confiável e muito difundida na literatura. A análise do teste de estabilidade pela simples utilização do método do enxame de partículas associado à técnica de polarização mencionada acima, para a obtenção de mais de um ponto estacionário (sem a busca local feita pelo método do conjunto gerador em uma dada rede de pontos), constitui outra metodologia para a resolução do problema de interesse. Essa utilização é uma novidade secundária deste trabalho. Tal metodologia simplificada exibiu também uma grande robustez, sendo capaz de encontrar todos os pontos estacionários pesquisados. No entanto, quando comparada com a abordagem mais geral proposta aqui, observou-se que tal simplificação pode, em alguns casos onde a função de mérito apresenta uma geometria mais complexa, consumir um tempo de máquina relativamente grande, dessa forma é menos eficiente.

Resolução numérica de escoamentos compressíveis empregando um método de partículas livre de malhas e o processamento em paralelo (CUDA)

Os métodos numéricos convencionais, baseados em malhas, têm sido amplamente aplicados na resolução de problemas da Dinâmica dos Fluidos Computacional. Entretanto, em problemas de escoamento de fluidos que envolvem superfícies livres, grandes explosões, grandes deformações, descontinuidades, ondas de choque etc., estes métodos podem apresentar algumas dificuldades práticas quando da resolução destes problemas. Como uma alternativa viável, existem os métodos de partículas livre de malhas. Neste trabalho é feita uma introdução ao método Lagrangeano de partículas, livre de malhas, Smoothed Particle Hydrodynamics (SPH) voltado para a simulação numérica de escoamentos de fluidos newtonianos compressíveis e quase-incompressíveis. Dois códigos numéricos foram desenvolvidos, uma versão serial e outra em paralelo, empregando a linguagem de programação C/C++ e a Compute Unified Device Architecture (CUDA), que possibilita o processamento em paralelo empregando os núcleos das Graphics Processing Units (GPUs) das placas de vídeo da NVIDIA Corporation. Os resultados numéricos foram validados e a eficiência computacional avaliada considerandose a resolução dos problemas unidimensionais Shock Tube e Blast Wave e bidimensional da Cavidade (Shear Driven Cavity Problem).

Um método sintético de difusão para aceleração do esquema de fonte de espalhamento em cálculos SN unidimensionais de fonte fixa

O esquema iterativo de fonte de espalhamento (SI) é tradicionalmente aplicado para a convergência da solução numérica de malha fina para problemas de transporte de nêutrons monoenergéticos na formulação de ordenadas discretas com fonte fixa. O esquema SI é muito simples de se implementar sob o ponto de vista computacional; porém, o esquema SI pode apresentar taxa de convergência muito lenta, principalmente para meios difusivos (baixa absorção) com vários livres caminhos médios de extensão. Nesta dissertação descrevemos uma técnica de aceleração baseada na melhoria da estimativa inicial para a distribuição da fonte de espalhamento no interior do domínio de solução. Em outras palavras, usamos como estimativa inicial para o fluxo escalar médio na grade de discretização de malha fina, presentes nos termos da fonte de espalhamento das equações discretizadas SN usadas nas varreduras de transporte, a solução numérica da equação da difusão de nêutrons em grade espacial de malha grossa com condições de contorno especiais, que aproximam as condições de contorno prescritas que são clássicas em cálculos SN, incluindo condições de contorno do tipo vácuo. Para aplicarmos esta solução gerada pela equação da difusão em grade de discretização de malha grossa nas equações discretizadas SN de transporte na grade de discretização de malha fina, primeiro implementamos uma reconstrução espacial dentro de cada nodo de discretização, e então determinamos o fluxo escalar médio em grade de discretização de malha fina para usá-lo nos termos da fonte de espalhamento. Consideramos um número de experimentos numéricos para ilustrar a eficiência oferecida pela presente técnica (DSA) de aceleração sintética de difusão.

Um método de matriz resposta com esquema iterativo de inversão parcial por região para problemas unidimensionais de transporte de nêutrons monoenergéticos na formulação de ordenadas discretas

Um método de matriz resposta (RM) é descrito para gerar soluções numéricas livres de erros de truncamento espacial para problemas de transporte de nêutrons monoenergéticos e com fonte fixa, em geometria unidimensional na formulação de ordenadas discretas (SN). O método RM com esquema iterativo de inversão parcial por região (RBI) converge valores numéricos para os fluxos angulares nas fronteiras das regiões que coincidem com os valores da solução analítica das equações SN, afora os erros de arredondamento da aritmética finita computacional. Desenvolvemos um esquema numérico de reconstrução espacial, que fornece a saída para os fluxos escalares de nêutrons em qualquer ponto do domínio definido pelo usuário, com um passo de avanço também escolhido pelo usuário. Resultados numéricos são apresentados para ilustrar a precisão do presente método em cálculos de malha grossa.

Desenvolvimento de um método espectronodal livre de erros de truncamento espacial para problemas adjuntos de transporte de partículas neutras monoenergéticas na formulação de ordenadas discretas em geometria unidimensional

Um método numérico nodal livre de erros de truncamento espacial é desenvolvido para problemas adjuntos de transporte de partículas neutras monoenergéticas em geometria unidimensional com fonte fixa na formulação de ordenadas discretas (SN). As incógnitas no método são os fluxos angulares adjuntos médios nos nodos e os fluxos angulares adjuntos nas fronteiras dos nodos, e os valores numéricos gerados para essas quantidades são os obtidos a partir da solução analítica das equações SN adjuntas. O método é fundamentado no uso da convencional equação adjunta SN discretizada de balanço espacial, que é válida para cada nodo de discretização espacial e para cada direção discreta da quadratura angular, e de uma equação auxiliar adjunta não convencional, que contém uma função de Green para os fluxos angulares adjuntos médios nos nodos em termos dos fluxos angulares adjuntos emergentes das fronteiras dos nodos e da fonte adjunta interior. Resultados numéricos são fornecidos para ilustrarem a precisão do método proposto.

Intensidade acústica útil: um novo método para identificação de regiões radiantes em superfícies com geometrias arbitrárias

Neste trabalho é descrita a teoria necessária para a obtenção da grandeza denominada intensidade supersônica, a qual tem por objetivo identificar as regiões de uma fonte de ruído que efetivamente contribuem para a potência sonora, filtrando, consequentemente, a parcela referente às ondas sonoras recirculantes e evanescentes. É apresentada a abordagem de Fourier para a obtenção da intensidade supersônica em fontes com geometrias separáveis e a formulação numérica existente para a obtenção de um equivalente à intensidade supersônica em fontes sonoras com geometrias arbitrárias. Este trabalho apresenta como principal contribuição original, uma técnica para o cálculo de um equivalente à intensidade supersônica, denominado aqui de intensidade acústica útil, capaz de identificar as regiões de uma superfície vibrante de geometria arbitrária que efetivamente contribuem para a potência sonora que será radiada. Ao contrário da formulação numérica existente, o modelo proposto é mais direto, totalmente formulado na superfície vibrante, onde a potência sonora é obtida através de um operador (uma matriz) que relaciona a potência sonora radiada com a distribuição de velocidade normal à superfície vibrante, obtida com o uso do método de elementos finitos. Tal operador, chamado aqui de operador de potência, é Hermitiano, fato crucial para a obtenção da intensidade acússtica útil, após a aplicação da decomposição em autovalores e autovetores no operador de potência, e do critério de truncamento proposto. Exemplos de aplicações da intensidade acústica útil em superfícies vibrantes com a geometria de uma placa, de um cilindro com tampas e de um silenciador automotivo são apresentados, e os resultados são comparados com os obtidos via intensidade supersônica (placa) e via técnica numérica existente (cilindro), evidenciando que a intensidade acústica útil traz, como benefício adicional, uma redução em relação ao tempo computacional quando comparada com a técnica numérica existente.

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.

Desenvolvimento de um simulador numérico empregando o método Smoothed Particle Hydrodynamics para a resolução de escoamentos incompressíveis. Implementação computacional em paralelo (CUDA)

Neste trabalho, foi desenvolvido um simulador numérico baseado no método livre de malhas Smoothed Particle Hydrodynamics (SPH) para a resolução de escoamentos de fluidos newtonianos incompressíveis. Diferentemente da maioria das versões existentes deste método, o código numérico faz uso de uma técnica iterativa na determinação do campo de pressões. Este procedimento emprega a forma diferencial de uma equação de estado para um fluido compressível e a equação da continuidade a fim de que a correção da pressão seja determinada. Uma versão paralelizada do simulador numérico foi implementada usando a linguagem de programação C/C++ e a Compute Unified Device Architecture (CUDA) da NVIDIA Corporation. Foram simulados três problemas, o problema unidimensional do escoamento de Couette e os problemas bidimensionais do escoamento no interior de uma Cavidade (Shear Driven Cavity Problem) e da Quebra de Barragem (Dambreak).