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.

Metrologia de instrumentos endodônticos por microscopia eletrônica de varredura e análise digital de imagens: uma proposta metodológica

A padronização para a fabricação de instrumentos endodônticos em aço inoxidável contribuiu para o desenvolvimento de novos aspectos geométricos. Surgiram propostas de alterações no desenho da haste helicoidal, da seção reta transversal, da ponta, da conicidade e do diâmetro na extremidade (D0). Concomitantemente, o emprego de ligas em Níquel-Titânio possibilitou a produção de instrumentos acionados a motor, largamente empregados hoje. A cada ano a indústria lança instrumentos com diversas modificações, sem, contudo, disponibilizar informações suficientes quanto às implicações clínicas destas modificações. Existe um crescente interesse no estudo dos diferentes aspectos geométricos e sua precisa metrologia. Tradicionalmente, a aferição de aspectos geométricos de instrumentos endodônticos é realizada visualmente através de microscopia ótica. Entretanto, esse procedimento visual é lento e subjetivo. Este trabalho propõe um novo método para a metrologia de instrumentos endodônticos baseado no microscópio eletrônico de varredura e na análise digital das imagens. A profundidade de campo do MEV permite obter a imagem de todo o relevo do instrumento endodôntico a uma distância de trabalho constante. Além disso, as imagens obtidas pelo detector de elétrons retro-espalhados possuem menos artefatos e sombras, tornando a obtenção e análise das imagens mais fáceis. Adicionalmente a análise das imagens permite formas de mensuração mais eficientes, com maior velocidade e qualidade. Um porta-amostras específico foi adaptado para obtenção das imagens dos instrumentos endodônticos. Ele é composto de um conector elétrico múltiplo com terminais parafusados de 12 pólos com 4 mm de diâmetro, numa base de alumínio coberta por discos de ouro. Os nichos do conector (terminais fêmeas) têm diâmetro apropriado (2,5 mm) para o encaixe dos instrumentos endodônticos. Outrossim, o posicionamento ordenado dos referidos instrumentos no conector elétrico permite a aquisição automatizada das imagens no MEV. Os alvos de ouro produzem, nas imagens de elétrons retro-espalhados, melhor contraste de número atômico entre o fundo em ouro e os instrumentos. No porta-amostras desenvolvido, os discos que compõem o fundo em ouro são na verdade, alvos do aparelho metalizador, comumente encontrados em laboratórios de MEV. Para cada instrumento, imagens de quatro a seis campos adjacentes de 100X de aumento são automaticamente obtidas para cobrir todo o comprimento do instrumento com a magnificação e resolução requeridas (3,12 m/pixel). As imagens obtidas são processadas e analisadas pelos programas Axiovision e KS400. Primeiro elas são dispostas num campo único estendido de cada instrumento por um procedimento de alinhamento semi-automático baseado na inter-relação com o Axiovision. Então a imagem de cada instrumento passa por uma rotina automatizada de análise de imagens no KS400. A rotina segue uma sequência padrão: pré-processamento, segmentação, pós-processamento e mensuração dos aspectos geométricos.

Segmentação e classificação de imagens destinadas ao mapeamento digital do Parque Nacional da Serra dos Órgãos (PARNASO)

Na perspectiva ambiental, o Parque Nacional da Serra dos Órgãos (PARNASO) é o Parque Nacional mais pesquisado no Brasil e configura-se como uma importante Unidade de Conservação inserida no estado do Rio de Janeiro, devido à sua importância ambiental para o estado. Localizado em quatro municípios da região serrana: Teresópolis, Petrópolis, Magé e Guapimirim foi constatado que essa área tem passado por alguns problemas, relativamente recentes, de ocupação desordenada devido à expansão urbana em sua vizinhança caracterizada por pressão antrópica. Através do processamento de imagens digitais, mais especificamente as etapas de segmentação e classificação, foi possível ilustrar o processo de ocupação humana por meio de documentos cartográficos. Além de estes processos possibilitarem a geração de mapas de uso da Terra e cobertura vegetal, com o intuito de auxiliar e dar fomento à execução de atividades, o mapeamento digital configura-se numa importante ferramenta para a análise ambiental, contribuindo para o posterior zoneamento da área de estudo. Adotaram-se classes temáticas de uso e ocupação da Terra com o propósito de permitir a classificação das imagens digitais trabalhadas. São elas: afloramento rochoso, área urbana, agricultura e vegetação. Estudos foram feitos no sentido de indicar e explorar as funcionalidades das ferramentas SPRING e DEFINIENS e resultados foram comparados a partir do uso de imagens LANDSAT, CBERS, SPOT e IKONOS chegando-se a resultados de que no sistema SPRING, os melhores parâmetros a serem escolhidos foram similaridade 10 e área 400. Já para o sistema DEFINIENS, constatou-se que o processo de segmentação multinível permitiu o alcance de resultados mais rápidos, do ponto de vista computacional, do que o processo de segmentação único utilizado normalmente entre os sistemas de processamento de imagens digitais como o SPRING. Já sob a ótica do processo de classificação de imagens, a pesquisa constituiu em avaliar este mecanismo por meio de dois indicadores: o de exatidão/acurácia e o índice Kappa. Neste sentido, observaram-se tendências de melhores resultados no sistema SPRING.

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.

Avaliação não supervisionada do processo de segmentação de imagens utilizadas em geociências

Com a necessidade de extrair as informações contidas nas imagens de satélite de forma rápida, eficiente e econômica, são utilizadas cada vez mais as técnicas computacionais de processamento de imagens como a de segmentação automática. Segmentar uma imagem consiste em dividí-la em regiões através de um critério de similaridade, onde os pixels que estão contidos nestas possuem características semelhantes, como por exemplo, nível de cinza, textura, ou seja, a que melhor represente os objetos presentes na imagem. Existem vários exemplos de algoritmos segmentadores, como o de crescimento de regiões onde os pixels crescem e são aglutinados formando regiões. Para determinar quais os melhores parâmetros utilizados nestes algoritmos segmentadores é necessário que se avalie os resultados a partir dos métodos mais utilizados, que são os supervisionados onde há necessidade de uma imagem de referência, considerada ideal fazendo com que se tenha um conhecimento a priori da região de estudo. Os não supervisionados, onde não há a necessidade de uma imagem de referência, fazendo com que o usuário economize tempo. Devido à dificuldade de se obter avaliadores para diferentes tipos de imagem, é proposta a metodologia que permite avaliar imagens que possuam áreas com vegetação, onde serão formadas grandes regiões (Crianass) e o que avaliará as imagens com áreas urbanas onde será necessário mais detalhamento (Cranassir).

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).