Construção de quadrados mágicos pelo método do passo uniforme

Lehmer (1929) analisa matematicamente o método do passo uniforme para construção de quadrados mágicos de ordem impar. Ele divide sua análise em várias etapas. Na primeira delas, envolvendo a discussão de condições necessárias e suficientes para o preenchimento do quadrado pelo método, o autor afirma que se dois números guardarem entre si uma certa relação, eles serão designados a ocupar a mesma célula do quadrado causando seu não preenchimento. A análise do preenchimento pelo método do passo uniforme envolve a resolução de um sistema linear módulo n. Nesse trabalho, discutimos o comportamento das soluções desse sistema quando o método falha no preenchimento. Como consequência, concluímos que números que guardam a relação mencionada nunca ocupam a mesma célula. A análise das condições necessárias e suficientes para obter quadrados mágicos segundo a definição de Lehmer (1929) envolve a resolução de equações de congruências lineares a duas variáveis. Nesse trabalho, detalhamos os resultados de Lehmer (1929). A análise das condições necessárias e suficientes para obtenção de quadrados mágicos, como são reconhecidos usualmente, também envolve a resolução de equações de congruências lineares a duas variáveis. Discutimos o comportamento das soluções dessas equações para obter diagonais principais mágicas. Como consequência, mostramos que diagonais principais mágicas são obtidas se e somente se as coordenadas iniciais guardarem certas relações

Simulação numérica de escoamentos de fluídos pelo método de elementos finitos baseado em volumes de controle com a técnica de passo fracionado

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

Um método de linearização local com passo adaptativo para solução numérica de equações diferenciais estocásticas com ruído aditivo

Neste trabalho apresentamos um novo método numérico com passo adaptativo baseado na abordagem de linearização local, para a integração de equações diferenciais estocásticas com ruído aditivo. Propomos, também, um esquema computacional que permite a implementação eficiente deste método, adaptando adequadamente o algorítimo de Padé com a estratégia “scaling-squaring” para o cálculo das exponenciais de matrizes envolvidas. Antes de introduzirmos a construção deste método, apresentaremos de forma breve o que são equações diferenciais estocásticas, a matemática que as fundamenta, a sua relevância para a modelagem dos mais diversos fenômenos, e a importância da utilização de métodos numéricos para avaliar tais equações. Também é feito um breve estudo sobre estabilidade numérica. Com isto, pretendemos introduzir as bases necessárias para a construção do novo método/esquema. Ao final, vários experimentos numéricos são realizados para mostrar, de forma prática, a eficácia do método proposto, e compará-lo com outros métodos usualmente utilizados.

Esquemas alternativos para o passo de parametrização do método da continuação baseados em parâmetros físicos

O método de fluxo de carga convencional é considerado inadequado para se obter o ponto de máximo carregamento (PMC) de sistemas de potência, devido à singularidade da matriz Jacobiana neste ponto. Os métodos da continuação são ferramentas eficientes para a solução deste tipo de problema, visto que técnicas de parametrização podem ser utilizadas para evitar a singularidade da matriz Jacobiana. Neste trabalho, novas opções para a etapa de parametrização do método da continuação são apresentadas. Mostra-se que variáveis com claro significado físico podem ser utilizadas na etapa de parametrização. As seguintes variáveis foram testadas: perda total de potência ativa e reativa, potência ativa e reativa na barra de referência, potência reativa das barras de geração, e as perdas de potência ativa e reativa nas linhas de transmissão (LT). Além de facilitar a implementação computacional do método de continuação, as técnicas de parametrização apresentadas simplificam a definição matemática e o entendimento do método por parte de engenheiros de potência, visto que os métodos de continuação existentes na literatura sempre utilizam técnicas de parametrização complexas, e de interpretação puramente geométrica. Resultados obtidos com a nova metodologia para os sistemas testes do IEEE (14, 30, 57 e 118 barras) mostram que as características de convergência do método de fluxo de carga convencional são melhoradas na região do PMC. Além disso, durante o traçado das curvas PV, as diversas técnicas de parametrização podem ser comutadas entre si possibilitando o cálculo de todos os pontos da curva com um número reduzido de iterações. Diversos testes são realizados para proporcionar a comparação do desempenho dos esquemas de parametrização propostos.

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.

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.

Método numérico de Matriz Resposta acoplado a um esquema de reconstrução espacial analítica para cálculos unidimensionais de transporte de nêutrons na formulação de ordenadas discretas multigrupo de energia com fonte fixa

Um método de Matriz Resposta (MR) é descrito para gerar soluções numéricas livres de erros de truncamento espacial para problemas multigrupo de transporte de nêutrons com fonte fixa e em geometria unidimensional na formulação de ordenadas discretas (SN). Portanto, o método multigrupo MR com esquema iterativo de inversão nodal parcial (NBI) 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 multigrupo SN, afora os erros de arredondamento da aritmética finita computacional. É também desenvolvido um esquema numérico de reconstrução espacial, que fornece a saída para os fluxos escalares de nêutrons em cada grupo de energia em um intervalo qualquer 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.

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

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

Elementos finitos em fluidos dominados pelo fenômeno de advecção: um método semi-Lagrangeano.

Os escoamentos altamente convectivos representam um desafio na simulação pelo método de elementos finitos. Com a solução de elementos finitos de Galerkin para escoamentos incompressíveis, a matriz associada ao termo convectivo é não simétrica, e portanto, a propiedade de aproximação ótima é perdida. Na prática as soluções apresentam oscilações espúrias. Muitos métodos foram desenvolvidos com o fim de resolver esse problema. Neste trabalho apresentamos um método semi- Lagrangeano, o qual é implicitamente um método do tipo upwind, que portanto resolve o problema anterior, e comparamos o desempenho do método na solução das equações de convecção-difusão e Navier-Stokes incompressível com o Streamline Upwind Petrov Galerkin (SUPG), um método estabilizador de reconhecido desempenho. No SUPG, as funções de forma e de teste são tomadas em espaços diferentes, criando um efeito tal que as oscilações espúrias são drasticamente atenuadas. O método semi-Lagrangeano é um método de fator de integração, no qual o fator é um operador de convecção que se desloca para um sistema de coordenadas móveis no fluido, mas restabelece o sistema de coordenadas Lagrangeanas depois de cada passo de tempo. Isto prevê estabilidade e a possibilidade de utilizar passos de tempo maiores.Existem muitos trabalhos na literatura analisando métodos estabilizadores, mas não assim com o método semi-Lagrangeano, o que representa a contribuição principal deste trabalho: reconhecer as virtudes e as fraquezas do método semi-Lagrangeano em escoamentos dominados pelo fenômeno de convecção.

Método de avaliação de tolerância à toxidez de alumínio em trigo, em condições de hidroponia, na Embrapa Trigo.

2006

L'ensenyament a Catalunya durant el trienni liberal (1820-1823) : el mètode dels escolapis. 'La enseñanza en Cataluña durante el trienio liberal (1820-1823) : el método de los escolapios'.

Estudiar la realidad educativa durante el trienio liberal en Cataluña y las posibles relaciones entre la escuela o la enseñanza y el momento de cambio que vive el país.. Enseñanza durante el trienio liberal en Cataluña.. Hace una exposición del marco histórico en el que se mueve el trienio y examina aspectos como la historia, la sociedad, la economía, la cultura, la iglesia y los escolapios. Presenta la legislación relacionada con la enseñanza. Expone el mapa escolar de Cataluña: escuelas de las comarcas, oferta escolar en la ciudad de Barcelona y los colegios escolapios. Explica características de los agentes educativos y del proceso educativo. Examina las diferentes posturas ideológicas: los ilustrados, los afrancesados, los liberales, la enseñanza tradicional y los escolapios. Aporta documentos como programaciones, memoriales, reglamentos, láminas con planos, etc.. Archivos históricos. Bibliografía. Era necesaria una legislación sobre la enseñanza, la Constitución sentaba los puntos que eran necesarios desarrollar. En 1821 se publica un Reglamento y en 1822 se da a conocer el Proyecto del Reglamento General de Primera Enseñanza. Los maestros no tenían ninguna formación específica, el maestro escolapio destacaba porque recibía una preparación en conocimientos y didáctica de las materias que había de enseñar. Loa alumnos pobres no podían asistir a la escuela. Los escolapios seguían el método uniforme..

Análise da interação solo-estrutura em fundações de concreto armado pelo método dos elementos finitos

A análise da iteração solo-estrutura em fundações é um importante campo de pesquisa que ainda tem um grande progresso a ser feito. No presente trabalho foi desenvolvido um programa computacional para a análise da interação solo-estrutura de fundações de concreto armado. Este tema abrange duas áreas da engenharia civil: estruturas e geotecnia. O método dos elementos finitos foi usado no trabalho na seqüência para resolver o problema considerando estado plano de defonnação e comportamento elastoplásti.co dos materiais estudados (solo, concreto e aço). A linguagem de programação MATLAB foi usada em toda esta pesquisa como alternativa ao FORTRAN. O MATLAB foi escolhido uma vez que é uma linguagem de programação que permite facilmente construir uma interfàce de pré e pósprocessamento amigável. Os passos para a solução completa do problema foram os seguintes: Primeiramente um programa foi desenvolvido considerando o comportamento elastoplástico com critérios de plastificação e ruptura específicos para o concreto e solo. Soluções analíticas fechadas foram usadas para checar a precisão do programa. O segundo passo foi a introdução do reforço de aço no concreto por meio de um modelo para armaduras. Logo após, um modelo de fissuras para o concreto 1racionado foi in1roduzido no programa. Na seqüência o programa de pré e pós-processamento foi desenvolvido para gerar a malha de elementos finitos (pré-processamento), distribuição tensões e deformações, mapa de fissuras, etc (pósprocessamento). Finalmente, os parâme1ros constitutivos do solo, concreto e aço foram calibrados e várias situações reais de interação do solo-concreto de fundações de concreto armado foram simuladas. Nesta dissertação são encontrados resultados para as pressões de contato sapata-solo. Diferentes diagramas de tensões de interfàce foram obtidos em função rigidez relativa do elemento estrutural de concreto armado-solo. Na análise mnnérica, rigidez relativa desempenhou uma relevante função no comportamento mecânico do elemento estrutural de concreto armado (sapata) e da base de assentamento (solo), uma vez ruptura em ambos os casos esteve diretamente relacionada a esta grandeza. São encon1rados, em função da rigidez relativa, resultados indicativos dos modos de falha da fundação, excessiva plastificação do solo em fundações com rigidez relativa alta, e a plastificação armaduras, esmagamento do concreto, formação de fissuras, bielas e confinamento concreto para fundações de rigidez relativa baixa. Na análise numérica, obteve-se resultados importantes com relação ao projeto de fundações. Estes resultados foram cOnITontadoscom normas, destacando-se as discordâncias com relação às recomendações da nonna brasileira Projeto e Execução de Fundações" NBR-6122 (1996) para os diagramas de tensões interface sapata-solo usados no dimensionamento de fundações de concreto armado.

Simulação numérica de esteiras em transição utilizando o método dos contornos virtuais

Simulações Numéricas são executadas em um código numérico de alta precisão resolvendo as equações de Navier-Stokes e da continuidade para regimes de escoamento incompressíveis num contexto da turbulência bidimensional. Este código utiliza um esquema compacto de diferenças finitas de sexta ordem na aproximação das derivadas espaciais. As derivadas temporais são calculadas usando o esquema de Runge-Kuta de terceeira ordem com baixo armazenamento. Tal código numérico fornece uma representação melhorada para uma grande faixa de escalas de comprimento e de tempo. As técnicas dos contornos imersos acopladas ao método dos contornos virtuais permitem modelar escoamentos não-estacionários sobre geometrrias complexas, usando simplesmente uma malha Cartesiana uniforme. Por meio de procedimentos de aproximação/interpolação, as técnicas dos contornos imersos (aproximação Gaussiana, interpolação bilinear e redistribuição Gaussiana), permitem a representação do corpo sólido no interior do campo de escoamento, com a superfície não coincidindo com a malha computacional. O método dos contornos virtuais, proposto originalmente por Peskin, consiste, basicamente, na imposição na superfície e/ou no interior do corpo, de um termo de força temporal acrescentando às equações do momento. A aplicação deste campo de força local leva o fluido ao repouso na superfície do corpo, permitindo obter as condições de contorno de não-deslizamento e de não penetração de fluido na parede. A análise das oscilações induzidas no escoamento-contorno pelo processo de desprendimento de vórtices na esteira do cilindro circular e de geometria retangulares na incidência, para números de Reybolds variando de 40 a 400, confirma a eficiência computacional e a aplicabilidade das técncias implementadas.

FRIDA: um método para elicitação e modelagem de RNFs

Os requisitos direcionam o desenvolvimento de software porque são cruciais para a sua qualidade. Como conseqüência tanto requisitos funcionais quanto não funcionais devem ser identificados o mais cedo possível e sua elicitação deve ser precisa e completa. Os requisitos funcionais exercem um papel importante uma vez que expressam os serviços esperados pela aplicação. Por outro lado, os requisitos não funcionais estão relacionados com as restrições e propriedades aplicadas ao software. Este trabalho descreve como identificar requisitos não funcionais e seu mapeamento para aspectos. O desenvolvimento de software orientado a aspectos é apontado como a solução para os problemas envolvidos na elicitação e modelagem dos requisitos não funcionais. No modelo orientado a aspectos, o aspecto é considerado o elemento de primeira ordem, onde o software pode ser modelado com classes e aspectos. As classes são comumente usadas para modelar e implementar os requisitos funcionais, já os aspectos são adotados para a modelagem e implementação dos requisitos não funcionais. Desse modo, é proposta a modelagem dos requisitos não funcionais através das fases do ciclo de vida do software, desde as primeiras etapas do processo de desenvolvimento. Este trabalho apresenta o método chamado FRIDA – From RequIrements to Design using Aspects, cujo objetivo é determinar uma forma sistemática para elicitar e modelar tanto os requisitos funcionais quanto os não funcionais, desde as fases iniciais do ciclo de desenvolvimento. Em FRIDA, a elicitação dos requisitos não funcionais é realizada usando-se checklists e léxicos, os quais auxiliam o desenvolvedor a descobrir os aspectos globais – utilizados por toda a aplicação – bem como, os aspectos parciais que podem ser empregados somente a algumas partes da aplicação. O próximo passo consiste na identificação dos possíveis conflitos gerados entre aspectos e como resolvê-los. No método FRIDA, a identificação e resolução de conflitos é tão importante quanto a elicitação de requisitos não funcionais, nas primeiras fases do ciclo de vida do software. Além disso, é descrito como usar a matriz de conflitos para automatizar esse processo sempre que possível. A extração dos aspectos e sua modelagem visual são características muito importantes, suportadas pelo método, porque elas possibilitam a criação de modelos que podem ser reutilizados no futuro. Em FRIDA, é demonstrado como transformar os requisitos em elementos da fase de projeto (classes e aspectos) e como traduzir esses elementos em código. Outra característica do método FRIDA é que a conexão entre diagramas, que pertencem a diferentes fases do processo de desenvolvimento do software, permite um alto nível de rastreabilidade. Em resumo, FRIDA requer que o desenvolvedor migre de uma visão puramente funcional para outra que contemple também os requisitos não funcionais.

A base dinâmica na existência de autovalores duplos no modelo de Timoshenko para uma viga uniforme livre-livre

Considerando uma viga uniforme do tipo Timoshenko com condições de contorno livre-livre, Geist e McLaughlin em [8]apresentam uma condição necessária e suficiente que garante a existência de freqüências naturais duplas. Esta condição foi obtida usando a formulação espectral, método clássico encontrado na literatura, para as equações de quarta ordem desacopladas do modelo de Timoshenko. O método clássico requer a obtenção de um vetor constante com oito componentes para que a solução deste modelo seja conhecida. Segundo Claeyssen [2], [3], [4], [5], [6], a solução do modelo de Timoshenko pode ser obtida usando a base dinâmica gerada por uma resposta impulso-matricial fundamental. Este método permite encontrar a solução do modelo de Timoshenko usando as equações de segunda ordem acopladas. Além disso, para que a solução seja conhecida é necessário obter um vetor constante com quatro componentes. O objetivo deste trabalho é estudar a condição necessária e suficiente que garante a existência de freqüências naturais duplas, apresentada por Geist e McLaughlin, para uma viga uniforme do tipo timoshenko com condições de contorno livre-livre e verificar se é possível obter esta mesma condição quando é utilizada a base dinâmica para obter a solução deste modelo.