Simulação tridimensional adaptativa da separação das fases de uma mistura bifásica usando a equação de Cahn-Hilliard

Simulamos a separação dos componentes de uma mistura bifásica com a equação de Cahn-Hilliard. Esta equação contém intrincados termos não lineares e derivadas de alta ordem. Além disso, a delgada região de transição entre os componentes da mistura requer muita resolução. Assim, determinar a solução numérica da equação de Cahn-Hilliard não é uma tarefa fácil, principalmente em três dimensões. Conseguimos a resolução exigida no tempo usando uma discretização semi-implícita de segunda ordem. No espaço, obtemos a precisão requerida utilizando malhas refinadas localmente com a estratégia AMR. Essas malhas se adaptam dinamicamente para recobrir a região de transição. O sistema linear proveniente da discretização é solucionado por intermédio de técnicas multinível-multigrid.

Um estudo sobre a equação de Schrödinger biharmônica

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

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.

Implementação de um algoritmo numérico para solução da equação de Christoffel generalizada em acustoelasticidade.

Extensos estudos realizados nas últimas décadas sobre a propagação de ondas ultrassônicas em sólidos levaram ao desenvolvimento de técnicas não destrutivas para a avaliação da segurança e integridade de estruturas e componentes industriais. O interesse na aplicação de técnicas ultrassônicas para medição de tensões aplicadas e residuais decorre da mudança mensurável da velocidade das ondas ultrassônicas na presença de um campo de tensões, fenômeno conhecido como efeito acustoelástico. Uma teoria de acustoelasticidade fornece um meio atrativo e não destrutivo de medir a tensão média ao longo do caminho percorrido pela onda. O estudo da propagação das ondas ultrassônicas em meios homogêneos anisotrópicos sob tensão conduz a um problema não linear de autovalores dado pela equação de Christoffel generalizada. A característica não linear deste problema decorre da interdependência entre as constantes elásticas efetivas do material e as tensões atuantes. A medição experimental de tensões por técnicas ultrassônicas é um problema inverso da acustoelasticidade. Esta dissertação apresenta a implementação de um algoritmo numérico, baseado no método proposto por Degtyar e Rokhlin, para solução do problema inverso da acustoelasticidade em sólidos ortotrópicos sujeitos a um estado plano de tensões. A solução da equação de Christoffel generalizada apresenta dificuldades de natureza numérica e prática. A estabilidade e a precisão do algoritmo desenvolvido, bem como a influência das incertezas na medição experimental das velocidades das ondas ultrassônicas, foram então investigadas. Dados sintéticos para as velocidades das ondas ultrassônicas de incidência oblíqua em uma placa sujeita a um estado plano de tensões foram gerados pela solução direta da equação de Christoffel generalizada para ilustrar a aplicação do algoritmo desenvolvido. O objetivo maior desta dissertação é a disponibilização de uma nova ferramenta de cálculo para suporte às atividades experimentais de medição de tensões por ultrassom no país.

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.

Solução analítica da equação unidimensional de transporte de nêutrons monoenergéticos com espalhamento linearmente anisotrópico e aproximação sintética de difusão

Nesta dissertação, são apresentados os seguintes modelos matemáticos de transporte de nêutrons: a equação linearizada de Boltzmann e a equação da difusão de nêutrons monoenergéticos em meios não-multiplicativos. Com o objetivo de determinar o período fluxo escalar de nêutrons, é descrito um método espectronodal que gera soluções numéricas para o problema de difusão em geometria planar de fonte fixa, que são livres de erros de truncamento espacial, e que conjugado com uma técnica de reconstrução espacial intranodal gera o perfil detalhado da solução. A fim de obter o valor aproximado do fluxo angular de nêutrons em um determinado ponto do domínio e em uma determinada direção de migração, descreve-se também um método de reconstrução angular baseado na solução analítica da equação unidimensional de transporte de nêutrons monoenergéticos com espalhamento linearmente anisotrópico com aproximação sintética de difusão nos termos de fonte por espalhamento. O código computacional desenvolvido nesta dissertação foi implementado na plataforma livre Scilab, e para ilustrar a eficiência do código criado,resultados numéricos obtidos para três problemas-modelos são apresentados

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.

Posicionamento aproximado do estado final para sistemas térmicos descritos pela equação do calor.

Neste trabalho, será considerado um problema de controle ótimo quadrático para a equação do calor em domínios retangulares com condição de fronteira do tipo Dirichlet é nos quais, a função de controle (dependente apenas no tempo) constitui um termo de fonte. Uma caracterização da solução ótima é obtida na forma de uma equação linear em um espaço de funções reais definidas no intervalo de tempo considerado. Em seguida, utiliza-se uma sequência de projeções em subespaços de dimensão finita para obter aproximações para o controle ótimo, o cada uma das quais pode ser gerada por um sistema linear de dimensão finita. A sequência de soluções aproximadas assim obtidas converge para a solução ótima do problema original. Finalmente, são apresentados resultados numéricos para domínios espaciais de dimensão 1.

Formação de nanopadrões em superfícies por sputtering iônico: Estudo numérico da equação anisotrópica amortecida de Kuramoto-Sivashinsky.

Apresenta-se uma abordagemnumérica para ummodelo que descreve a formação de padrões por sputtering iônico na superfície de ummaterial. Esse processo é responsável pela formação de padrões inesperadamente organizados, como ondulações, nanopontos e filas hexagonais de nanoburacos. Uma análise numérica de padrões preexistentes é proposta para investigar a dinâmica na superfície, baseada em ummodelo resumido em uma equação anisotrópica amortecida de Kuramoto-Sivashinsky, em uma superfície bidimensional com condições de contorno periódicas. Apesar de determinística, seu caráter altamente não-linear fornece uma rica gama de resultados, sendo possível descrever acuradamente diferentes padrões. Umesquema semi implícito de diferenças finitas com fatoração no tempo é aplicado na discretização da equação governante. Simulações foram realizadas com coeficientes realísticos relacionados aos parâmetros físicos (anisotropias, orientação do feixe, difusão). A estabilidade do esquema numérico foi analisada por testes de passo de tempo e espaçamento de malha, enquanto a verificação do mesmo foi realizada pelo Método das Soluções Manufaturadas. Ondulações e padrões hexagonais foram obtidos a partir de condições iniciais monomodais para determinados valores do coeficiente de amortecimento, enquanto caos espaço-temporal apareceu para valores inferiores. Os efeitos anisotrópicos na formação de padrões foramestudados, variando o ângulo de incidência.

Aplicação da equação de Poisson-Boltzmann modificada em sistemas biológicos: análise da partição iônica em um eritrócito

Neste trabalho, a partição iônica e o potencial de membrana em um eritrócito são analisados via equação de Poisson-Boltzmann modificada, considerando as interações não eletrostáticas presentes entre os íons e macromoléculas, assim como, o potencial β. Este potencial é atribuído à diferença de potencial químico de referência entre os meios intracelular e extracelular e ao transporte ativo de íons. O potencial de Gibbs-Donnan via equação de Poisson-Boltzmann na presença de carga fixa em um sistema contendo uma membrana semipermeável também é estudado. O método de aproximação paraboloide em elementos finitos em um sistema estacionário e unidimensionalé aplicado para resolver a equação de Poisson-Boltzmann em coordenadas cartesianas e esféricas. O parâmetro de dispersão relativo às interações não eletrostáticas écalculado via teoria de Lifshitz. Os resultados em relação ao potencial de Gibbs-Donnan mostram-se adequados, podendo ser calculado pela equação de Poisson-Boltzmann. No sistema contendo um eritrócito, quando o potencial β é considerado igual a zero, não se verifica a diferença iônica observada experimentalmente entre os meios intracelular e extracelular. Dessa forma, os potenciais não eletrostáticos calculados via teoria de Lifshitz têm apenas uma pequena influência no que se refere à alta concentração de íon K+ no meio intracelular em relação ao íon Na+

Equação de referência do Teste do Degrau medida em distância para indivíduos saudáveis e sedentários

O teste de caminhada em seis minutos (TC6M) avalia a capacidade respiratória durante o exercício. Recentemente, o teste do degrau em seis minutos (TD6M) está sendo estudado como uma proposta para essa mesma avaliação. Diante do exposto, o desenvolvimento de uma equação de referência se torna importante. O objetivo desse estudo foi desenvolver uma equação de referência padrão para o Teste do Degrau em Seis minutos. Esse estudo foi do tipo transversal, em que foram selecionados 452 indivíduos. Após a aplicação dos critérios de inclusão/exclusão, foram selecionados 326 sujeitos saudáveis e sedentários com idade entre 20 e 80 anos. Para serem considerados saudáveis, os participantes não podiam ter história de doenças (exceto hipertensão arterial sistêmica ou diabetes mellitus em tratamento) e foram submetidos à radiografia de tórax, espirometria e eletrocardiograma, que deveriam ser normais. O nível de sedentarismo foi avaliado através do International Physical Activity Questionnaire (IPAQ). Foram coletados os seguintes dados demográficos: idade, peso e altura. Todos os indivíduos realizaram TD6M na sua própria cadência (autocadenciado) em um degrau de 16,5 cm de altura, 65 de largura e 30 cm de comprimento. O número de subidas e descidas foi contado por um pedômetro digital. Foram mensuradas a pressão arterial, a frequência cardíaca e a saturação de oxigênio, antes e depois do TD6M. A análise estatística foi realizada pelo software STATA 12.0, e as equações foram desenvolvidas pelo modelo estatístico de regressão linear múltipla. Como resultado, observou-se que os participantes apresentaram exames normais, sendo 135 homens e 191 mulheres. O IPAQ demonstrou 157 ativos, 114 irregularmente ativos e 40 inativos, sendo que 14 indivíduos não responderam. A análise da distância alcançada em relação à idade e à diferença da frequência cardíaca tanto para homens quanto para mulheres mostrou significância estatística, demonstrando a sua importância para o desenvolvimento das equações para cada gênero. Para mulheres: Distância (m) = 88,83 - [(FC final- FC inicial)* 0,23] (Idade *0,37); para homens: Distância (m)= 110,20 - [(FC final- FC inicial)*0,18] (Idade * 0,59). O estudo concluiu que as equações de referência desenvolvidas nesse estudo foram realizadas em uma população de indivíduos saudáveis e sedentários e pode ser usada como padrão de referência para o TD6M.

Equação de volume para araucárias centenárias da Reserva Florestal Embrapa/Epagri.

2010

PASSADO E FUTURO NO PRESENTE DOS PERCURSOS DE QUALIFICAÇÃO: A EQUAÇÃO DA APRENDIZAGEM

É ponto assente que o desenvolvimento social e económico de um país dependedo conhecimento e do nível educativo dos seus cidadãos. Num momento de crise económica profunda em que a sustentabilidade do ensino superior se encontra ameaçada e a sua organização com rumo incerto, importa, mais do que nunca, criar um espaço de reflexão em que a ciência se pronuncie, aproximando e colocando em evidência os contributos que em vários domínios se vão produzindo.

Reformulação da definição: conceito de "equação" para os alunos do 1º ciclo do ensino secundário em Angola

A presente dissertação intitulada Reformulação da definição: conceito de «equação» para os alunos do 1.º Ciclo do Ensino Secundário em Angola tem como objectivo analisar as definições que se encontram nos manuais do aluno de Matemática do ciclo em referência. Para materialização deste desiderato organizamos um corpus constituído por textos dos três manuais que constituem o 1.º Ciclo, designadamente, manual do aluno da 7.ª, 8.ª e 9.ª classe. A análise da definição foi feita com base nas linhas orientadoras para a elaboração de uma definição terminológica. Estas linhas orientadoras correspondem aos procedimentos que adoptámos nesta dissertação, para a análise e para a proposta de reformulação da definição terminológica. O interesse pelo tema surgiu da necessidade de contribuir para a melhoria do processo de ensino-aprendizagem da disciplina de Matemática. Tendo em conta que a definição terminológica visa em primeiro lugar estabilizar o conceito por meio da identificação de características que lhes são peculiares permitindo a diferenciação entre conceitos no seio de outros conceitos. É a partir das definições que os alunos terão acesso aos conceitos, podendo até relacionar os vários conceitos do domínio. Desta feita, tanto o professor como o aluno beneficiarão da existência de uma boa definição. O aluno, por intermédio da definição, chegará ao conceito; o professor, por sua vez, terá melhores ferramentas – as definições – para transmitir conhecimentos.

Solução da equação de transporte multidimensional em geometria cartesiana e meio infinito usando derivada fracionária

Neste trabalho, foi construída uma forma integral para a solução das equações de transporte em uma, duas e três dimensões, considerando o núcleo de espalhamento de Klein-Nishina, espalhamento isotrópico e o núcleo de espalhamento de Rutherford, respectivamente, seguindo a mesma idéia proposta em trabalhos recentes, nos quais foi construída uma solução para a equação de transporte de nêutrons em geometria cartesiana, usando derivada fracionária. A metodologia consiste em igualar a derivada fracionária do fluxo angular à equação integral, determinar a ordem da derivada fracionária comparando o núcleo da equação integral com o da definição de Riemann-Liouville. Essa formulação foi aplicada ao cálculo de dose absorvida. São apresentadas soluções geradas a partir do emprego do método da derivada fracionária e comparadas a resultados disponíveis na literatura.