Abordagens do tipo livre de jacobiana na simulação do escoamento de fluidos compressíveis em meios porosos

O desenvolvimento de software livre de Jacobiana para a resolução de problemas formulados por equações diferenciais parciais não-lineares é de interesse crescente para simular processos práticos de engenharia. Este trabalho utiliza o chamado algoritmo espectral livre de derivada para equações não-lineares na simulação de fluxos em meios porosos. O modelo aqui considerado é aquele empregado para descrever o deslocamento do fluido compressível miscível em meios porosos com fontes e sumidouros, onde a densidade da mistura de fluidos varia exponencialmente com a pressão. O algoritmo espectral utilizado é um método moderno para a solução de sistemas não-lineares de grande porte, o que não resolve sistemas lineares, nem usa qualquer informação explícita associados com a matriz Jacobiana, sendo uma abordagem livre de Jacobiana. Problemas bidimensionais são apresentados, juntamente com os resultados numéricos comparando o algoritmo espectral com um método de Newton inexato livre de Jacobiana. Os resultados deste trabalho mostram que este algoritmo espectral moderno é um método confiável e eficiente para a simulação de escoamentos compressíveis em meios porosos.

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 espectral das equações de transporte de nêutrons na formulação de ordenadas discretas em meios multiplicativos

É presentada nesta dissertação uma análise espectral das equações de transporte de nêutrons, independente do tempo, em geometria unidimensional e bidimensional, na formulação de ordenadas discretas (SN), utilizando o modelo de uma velocidade e multigrupo, considerando meios onde ocorrem o fenômeno da fissão nuclear. Esta análise espectral constitui-se na resolução de problemas de autovalores e respectivos autovetores, e reproduz a expressão para a solução geral analítica local das equações SN (para geometria unidimensional) ou das equações nodais integradas transversalmente (geometria retangular bidimensional) dentro de cada região homogeneizada do domínio espacial. Com a solução geral local determinada, métodos numéricos, tais como os métodos de matriz de resposta SN, podem ser derivados. Os resultados numéricos são gerados por programas de computadores implementados em MatLab, versão 2012, a fim de verificar a natureza dos autovalores e autovetores correspondentes no espaço real ou complexo.

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

Levantamento e estudo das ocorrências de grafita do Distrito Grafitífero Aracoiába-Baturité, CE

O Distrito Grafitífero Aracoiába-Baturité apresenta depósitos do tipo gnaisse grafitoso (minério disseminado) e veio (minério maciço) com diferentes origens genéticas e com características físicas e ambientes geológicos de formação próprios. O minério tipo gnaisse grafitoso é de origem sedimentar, singenético, com teores de 1,5 a 8% de C, que se distribuem ao longo de duas extensas faixas paralelas, hospedadas na Subunidade Baturité, que constitui um importante metalotecto regional. A associação de grafita metamórfica disseminada em metassedimentos da Sequência Acarápe constitui um geoindicador de antiga bacia sedimentar neoproterozóica e, também, pode ser considerado como zona de geosutura resultante do subsequente fechamento de um oceano primitivo. As rochas desta subunidade correspondem na paleogeografia da Sequência Acarápe aos fácies de sopé de talude e de planície abissal. O minério tipo veio (fluido depositado) é epigenético e, com teores entre 20% e 70% de C, forma corpos tabulares e bolsões, controlados em escala local por estruturas de alívio (falhas, fraturas, zonas de contato, eixos de dobras etc.) que permitiram a percolação de soluções penumatolíticas relacionadas ao corpo plutônico de Pedra Aguda. As variações dos valores das relações entre isótopos estáveis de carbono (δ13C) na grafita do minério disseminado são de -26,72 a -23,52 e do minério maciço de -27,03 a -20,83, revelando sinal de atividades biológicas (bioassinaturas) e permitem afirmar que a grafita das amostras acima são derivadas de matéria orgânica. Foram apresentados os principais guias de prospecção para grafita e testados os seguintes métodos geofísicos: Eletro-Resistividade; GPR - Ground Penetrating Radar; Magnetometria; VLF (Very Low Frequency); e Polarização Induzida Espectral (IPS) / Resistividade (ER). A conjugação dos métodos de Polarização Induzida Espectral (IPS) e Eletro Resistividade (ER) foi o que demonstrou a melhor eficiência. Com relação à determinação do teor de carbono por termogravimetria (ATG), que é o método mais utilizado para este elemento. Verificou-se, que as faixas de queima atribuídas ao carbono no minério do Distrito de Aracoiába-Baturité (340 a 570C e de 570 a 1050C) eram diferentes das faixas do minério de Minas Gerais (350C a 650C e 650C a 1.050C). Esta constatação indica a necessidade de se determinar previamente as faixas de temperatura para cada região pesquisada.

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.

Métodos espectrais de agrupamento

Os métodos espectrais são ferramentas úteis na análise de dados, sendo capazes de fornecer informações sobre a estrutura organizacional de dados. O agrupamento de dados utilizando métodos espectrais é comumente baseado em relações de similaridade definida entre os dados. O objetivo deste trabalho é estudar a capacidade de agrupamento de métodos espectrais e seu comportamento, em casos limites. Considera-se um conjunto de pontos no plano e usa-se a similaridade entre os nós como sendo o inverso da distância Euclidiana. Analisa-se a qual distância mínima, entre dois pontos centrais, o agrupamento espectral é capaz de reagrupar os dados em dois grupos distintos. Acessoriamente, estuda-se a capacidade de reagrupamento caso a dispersão entre os dados seja aumentada. Inicialmente foram realizados experimentos considerando uma distância fixa entre dois pontos, a partir dos quais os dados são gerados e, então, reduziu-se a distância entre estes pontos até que o método se tornasse incapaz de efetuar a separação dos pontos em dois grupos distintos. Em seguida, retomada a distância inicial, os dados foram gerados a partir da adição de uma perturbação normal, com variância crescente, e observou-se até que valor de variância o método fez a separação dos dados em dois grupos distintos de forma correta. A partir de um conjunto de pontos obtidos com a execução do algoritmo de evolução diferencial, para resolver um problema multimodal, testa-se a capacidade do método em separar os indivíduos em grupos diferentes.

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.