O presente trabalho tem como objetivo investigar o comportamento elasto-plástico em baixas deformações do fio máquina laminado em aço AISI 1045 destinado a trefilação e processos subsequentes de beneficiamento. Inicialmente, foi obtido material de quatro diferentes lotes de fabricação, sendo todas as amostras retiradas após a etapa de pré-endireitamento do processo de trefilação. O material foi caracterizado através de análise da composição química, metalográfica, determinação de tamanho de grão e segregação. Além disso, o comportamento elasto-plástico do material em ciclos de tração e compressão foi investigado, identificando-se o efeito Bauschinger para três diferentes temperaturas. Os modelos numéricos de Ramberg-Osgood e Ludwig-Hollomon foram utilizados para a definição de uma modelo numérico da curva para utilização em trabalhos futuros de simulação computacional. O modelo de Ludwig-Hollomon foi avaliado como mais satisfatório para baixas deformações e foi então utilizado para o cálculo da curva média de quatro lotes. A partir de curvas de escoamento de aço similar obtidas da literatura, foi possível obter-se uma extrapolação da curva para diferentes taxas de deformação, além da taxa de 0.05s-1, original dos ensaios: 1.6 s-1, 8 s-1, e 40 s-1. / The purpose of this work is to investigate the elastic-plastic behavior at low strains of steel rolled coil AISI 1045 intended for cold-drawing and subsequent processes. First, the material was withdrawn from four different manufacturing batches and all samples were taken from the pre-straightening process step. The material was characterized by means of chemical and metallography analysis, grain size and segregation. Furthermore, the elastic-plastic behavior in tension and compression cycles was investigated and the Bauschinger effect was identified for three different temperatures. The numerical models Ramberg-Osgood and Ludwig-Hollomon were the models used in this work to define the flow curve which can be used in future computer simulation works. The analysis showed that the Ludwig-Hollomon model is more adequate for low strains ranges. Ludwig-Hollomin was used to calculate the average flow curve based on the curves of the four batches. By using literature flow curves of similar steels, it was possible to make an extrapolation and to obtain flow curves at different strain rates beyond the experimental one: 1.6 s-1, 8 s-1, e 40 s-1.

Neste trabalho, relatam-se soluções analíticas para as equações da cinética da teoria de difusão de nêutrons. Para a solução das equações da cinética pontual consideram-se seis grupos de precursores de nêutrons atrasados e assume-se reatividade variável como uma função arbitrária do tempo. A ideia principal consiste inicialmente na determinação da solução das equações da cinética pontual com reatividade constante apenas usando os resultados bem conhecidos para a solução de sistemas de equações diferenciais matriciais lineares de primeira ordem com entradas constantes. Com a aplicação do método de Decomposição, é possível transformar as equações da cinética pontual com reatividade variável com o tempo em um conjunto de problemas recursivos semelhantes às equações da cinética pontual com reatividade constante, o que pode ser resolvido diretamente com a técnica mencionada anteriormente. Para ilustração, apresentam-se simulações para as funções com reatividade constante, linear e senoidal, bem como comparações com resultados na literatura. Já com relação às equações da cinética espacial, consideram-se um e seis grupos de precursores de nêutrons atrasados, modelo multigrupo de energia, meio homogêneo e dimensões espaciais bi e tridimensionais. O formalismo do procedimento da solução é geral em relação ao número de grupos de energia, famílias de precursores de nêutrons atrasados e regiões com diferentes composições químicas. O fluxo rápido e térmico e as concentrações de nêutrons atrasados são expandidos em uma série de termos de autofunções que, pela aplicação da técnica da GITT, resulta em uma equação diferencial matricial de primeira ordem semelhante às equações de cinética pontual. Por esse motivo, a solução deste problema transformado segue o formalismo do método da Decomposição aplicado às equações da cinética pontual. Por fim, apresentam-se simulações numéricas e comparações com resultados disponíveis na literatura. / In this work we report analytical solutions for the neutron kinetics diffusion equations. For the solution of the point kinetics equations we consider six groups of delayed neutron precursors and assume that the reactivity is an arbitrary function of time. The main idea initially consists in the determination of the solution of the point kinetics equations with constant reactivity by just using the well-known results of the solution of systems of first-order linear ordinary differential equations in matrix form with constant matrix entries. Applying the decomposition method, we are able to transform the point kinetics equations with time dependent reactivity into a set of recursive problems similar to the point kinetics equations with constant reactivity, which can be directly solved by the above mentioned technique. For illustration, we also report simulations for constant, linear and sinusoidal reactivity time functions of time as well as comparisons with results published in the literature. As for the space kinetics equations we consider six groups of delayed neutron precursors, energy multigroup model, homogeneous media and two and three-dimensional geometries. The solution procedure formalism is general with respect to the number of energy groups, neutron precursor families and regions with different chemical compositions. The fast and thermal flux and the delayed neutron precursors concentrations are expanded in a series in terms of eigenfunctions that, upon insertion into the kinetics equation and upon taking moments, result in a first order linear differential matrix equation with source terms similar to the point kinetics equations. The solution of this transformed problem follows the formalism of the decomposition method applied to the point kinetics equations. We present numerical simulations and comparisons with available results in the literature.

A maioria dos líquidos encontrados na natureza são não newtonianos e o estudo do seu comportamento reológico tem uma importância significante em diferentes áreas da engenharia. Entre eles, existe uma classe de fluidos que exibem pequena deformação aparente quando sujeitos a um nível de tensões inferior a uma tensão limite de escoamento, referido como comportamento viscoplástico. Nesta classe de materiais, alguns apresentam também comportamento elástico quando submetidos a baixas taxas de cisalhamento. A presente Tese tem como objetivo o estudo numérico de escoamentos bidimensionais em regime permanente de fluidos elasto-viscoplásticos através de uma expansão-contração planar. O modelo mecânico é definido pelas equações de conservação de massa e de balanço de momentum acopladas ao modelo elasto-viscoplástico proposto nesta Tese. Esta modelagem é aproximada por um método de elementos finitos multi-campos estabilizado baseado na metodologia de Galerkin mínimos-quadrados que possui como variáveis primais os campos de tensão extra polimérica, velocidade e pressão. As condições de compatibilidade entre os sub-espaços de elementos finitos para tensão extra-velocidade e velocidade-pressão são violadas, permitindo assim a utilização de interpolações de igual ordem. O método estabilizado foi implementado no código de elementos finitos para fluidos não newtonianos em desenvolvimento no Laboratório de Mecânica dos Fluidos Aplicada e Computacional (LAMAC) da Universidade Federal do Rio Grande do Sul. Nesta Tese é adotada uma metodologia alternativa para a definição das zonas rígidas do escoamento como sendo a posição onde a taxa de cisalhamento é igual a um valor dado pela relação de parâmetros reológicos do fluido, especificamente a tensão limite de escoamento e a viscosidade newtoniana para baixas taxas de cisalhamento. Nas simulações numéricas realizadas, o tempo de relaxação adimensional, o número de salto, o coeficiente power-law, a vazão adimensional e a massa específica adimensional são variados de forma a avaliar de que modo influenciam na dinâmica de escoamentos elastoviscoplásticos. Os resultados obtidos estão qualitativamente de acordo com a literatura, atestando a estabilidade da formulação empregada. / Non-Newtonian fluids are the majority of liquids found in nature and the study of their rheological behavior has a significant importance on different areas of engineering. Among them, there is a class of materials that exhibits little apparent deformation when subjected to a stress level behind an yield stress, referenced as viscoplastic material. In this class of materials, some fluids also exhibit elastic behavior at low shear rates. The present work aimed to a numerical study of two-dimensional steady state laminar flows of elasto-viscoplastic fluids through a planar expansion-contraction cavity. The mechanical model was defined by the mass conservation and momentum balance equations coupled to the elasto-viscoplastic model porposed in this work. This modeling has been approximated by a stabilized multi-field finite element method based on the Galerkin least-squares methodology, having as primal variables the elastic extra-stress component, velocity and pressure fields. In this way, the compatibility conditions between the extra-stress-velocity and pressure-velocity (Babuška- Brezzi condition) finite element subspaces are violated, allowing to use equal-order finite element interpolations. The stabilized method has been implemented in the finite element code for non-Newtonian fluids under development at the Laboratory of Applied and Computational Fluid Mechanics (LAMAC) of the Federal University of Rio Grande do Sul. An alternative methodology is adopted to define the yield surface as the position where the strain rate is equal to a value given by the relation of the rheological parameters of the fluid, namely the yield stress and the viscosity at low shear rates. In the performed numerical simulations, the non-dimensional relaxation time, the jump number, the power-law coefficient, the non-dimensional flow rate, and the non-dimensional density are varied in order to evaluate their influence on the elasto-viscoplastic fluid dynamics. All results found are in qualitatively accordance with the affine literature, and attesting the good stability features of the formulation.

Em problemas de viscoelasticidade computacional, a discretização espacial para a solução global das equações de equilíbrio é acoplada à discretização temporal para a solução de um problema de valor inicial local do fluxo viscoelástico. É demonstrado que este acoplamento espacial-temporal (ou global-local) éconsistente se o tensor de deformação total, agindo como elemento acoplador, tem uma aproximação de ordem p ao longo do tempo igual à ordem de convergência do método de integração de Runge-Kutta (RK). Para a interpolação da deformação foram utilizados polinômios baseados em soluções obtidas nos tempos tn+1, tn, . . ., tn+2−p, p ≥ 2, fornecendo dados consistentes de deformação nos estágios do RK. Em uma situação onde tal regra para a interpolação da deformação não é satisfeita, a integração no tempo apresentará, consequentemente, redução de ordem, baixa precisão e, por conseguinte, eficiência inferior. Em termos gerais, o propósito é generalizar esta condição de consistência proposta pela literatura, formalizando-a matematicamente e o demonstrando através da utilização de métodos de Runge-Kutta diagonalmente implícitos (DIRK) até ordem p = 4, aplicados a modelos viscoelásticos não-lineares sujeitos a deformações finitas. Através de exemplos numéricos, os algoritmos de integração temporal adaptados apresentaram ordem de convergência nominal e, portanto, comprovam a validade da formalização do conceito de interpolação consistente da deformação. Comparado com o método de integração de Euler implícito, é demonstrado que os métodos DIRK aqui aplicados apresentam um ganho considerável em eficiência, comprovado através dos fatores de aceleração atingidos. / In computational viscoelasticity, spatial discretization for the solution of the weak form of the balance of linear momentum is coupled to the temporal discretization for solving a local initial value problem (IVP) of the viscoelastic flow. It is shown that this spatial- temporal (or global-local) coupling is consistent if the total strain tensor, acting as the coupling agent, exhibits the same approximation of order p in time as the convergence order of the Runge-Kutta (RK) integration algorithm. To this end we construct interpolation polynomials based on data at tn+1, tn, . . ., tn+2−p, p ≥ 2, which provide consistent strain data at the RK stages. If this novel rule for strain interpolation is not satisfied, time integration shows order reduction, poor accuracy and therefore less efficiency. Generally, the objective is to propose a generalization of this consistency idea proposed in the literature, formalizing it mathematically and testing it using diagonally implicit Runge-Kutta methods (DIRK) up to order p = 4 applied to a nonlinear viscoelasticity model subjected to finite strain. In a set of numerical examples, the adapted time integrators obtain full convergence order and thus approve the novel concept of consistency. Substantially high speed-up factors confirm the improvement in the efficiency compared with Backward Euler algorithm.

O comportamento estrutural de veículos de transporte coletivo, quando submetido a eventos de impacto, é de grande relevância na engenharia automobilística, pois principalmente o ônibus, vem se tornando um dos meios de transporte de elevada importância, sendo que seu uso vem aumentando a cada ano. No entanto, as estruturas de ônibus que circulam nas estradas brasileiras não são capazes de resistir, sem causar danos aos passageiros, a eventos de impacto em acidentes, o que pode ser verificado em pesquisas divulgadas pela Agência Nacional de Transportes Terrestres (ANTT) em 2007, que mostra a evolução de acidentes e vítimas envolvendo ônibus, que cada ano aumentam consideravelmente. Neste contexto, pretende-se desenvolver um absorvedor de impacto que aumente a capacidade de absorção da energia de impacto e evite a invasão da região de sobrevivência dos passageiros, e desta forma diminuindo o número de vítimas em acidentes. Serão abordados neste trabalho eventos de impacto semi-frontal, justificado pela gravidade dos acidentes desse tipo nas estradas, sendo este o tipo de acidente que deixa um maior número de vítimas, pois geralmente parte da estrutura lateral do ônibus é removida, expondo os passageiros no período de impacto. Esse efeito é conhecido como “Efeito Abridor de Latas”. O modelo numérico deste trabalho será formado por elementos unifilares, sendo a estrutura do ônibus formada por vigas flexíveis e/ ou rígidas, unidas através de juntas não lineares rotacionais e translacionais. A rigidez de cada junta não linear é obtida a partir de métodos analíticos que descrevem o comportamento do tubo de parede fina quando submetido a impacto. / The structural behavior of transportation vehicles when subjected to events of impact is of great relevance in automotive engineering, mainly because the bus has become a highly important alternative means of transportation, and its use is increasing every year. However, the structures of the buses circulating on Brazilian roads are not able to withstand strong impact, without causing passengers’ injury in events where accidents occur, which can be shown in surveys released by the National Transport Agency (ANTT) in 2007, presenting that the great incidence of accidents with victims involving buses is increasing considerably every year. In this context, the development of an impact absorber that had the capability to increase the impact energy absorption and prevent the invasion of part of bus on the region where passengers are located (reducing the number of fatalities in accidents) are the objectives of this work. In this work, a semi-frontal impact will be presented and studied due to the severity of such accident, as the side of the bus structure is removed, exposing the passengers during the impact. This is the effect known as “Can opener effect”. The numerical model of this study will consist of a single-wire element, with a bus structure formed by flexible and/or rigid beams joined by nonlinear rotational and translational joints. The rigidity of each nonlinear joint is obtained by analytical methods describing the behavior of a thin-walled tube subjected to an impact force.

Fluidos Não-newtonianos estão presentes em grande quantidade na natureza e em muitas aplicações industriais. Tais fluidos apresentam uma variedade de comportamentos específicos. Dentre eles podemos citar: elasticidade, plasticidade, viscosidade variando com a taxa de deformação, tixotropia. Esta tese de doutorado apresenta um recente modelo para Fluidos Estruturados Tixotrópicos (de Souza Mendes, 2009), no qual baseasse nas tradicionais equações de balanço de momentum e massa, em uma equação constitutiva de Maxwell-b modificada e uma equação de evolução para o nível de estrutura do material. Na equação constitutiva de Maxwell-b modificada, tempo de relaxação e viscosidade estrutural são funções da estrutura do material. A equação evolutiva para o parâmetro de estruturação prevê efeitos de tixotropia, dependentes do tempo de equilíbrio do material. As equações que constituem o modelo mecânico para Fluidos Estruturados Tixotrópicos são aproximadas numericamente através do Método de Elementos Finitos de Galerkin Mínimos-quadrados (GLS). Tal método adicionou termos dependentes de malha ao método clássico de Galerkin, de forma a aumentar a estabilidade da formulação original, proporcionando compatibilizar os sub-espaços de velocidade e pressão. São investigados escoamentos creeping flow de Fluidos Estruturados Tixotrópicos em dois problemas clássicos de engenharia: em torno de um cilindro confinado entre duas placas e em uma expansão planar abrupta. Parâmetros relevantes ao modelo são variados em faixas pertinentes ao problema físico simulado e os resultados demonstraram uma boa capacidade de previsão. / Non-Newtonian fluids are present in large quantities in nature and in many industrial applications. Such fluids exhibit a variety of specific behaviors. Among them we can mention: elasticity, plasticity, viscosity varying with the strain rate, thixotropy, etc. This thesis presents a new model for thixotropic Structured Fluids (de Souza Mendes, 2009), in which were a basis in traditional balance equations of mass and momentum, the constitutive equation of Maxwell-type and an evolution equation for the level structure of the material. In the modified constitutive equation of Maxwell-b, relaxation time and structural viscosity are functions of material structure. The structure parameter evolutive equation provides effects of thixotropy, dependent on the equilibration time of the material. The equations which constitute the mechanical model for Structured Thixotropic fluids are approximated numerically using the Galerkin Least-squares Finite Element Method (GLS). Such a method added mesh dependent terms to the classical method of Galerkin in order to increase the stability of the original formulation, providing compatibility between sub-spaces of velocity and pressure. Are investigated Structured Thixotropic Fluids creeping flow in two classical problems in engineering: around a cylinder confined between two plates and in a planar abrupt expansion. Relevant parameters to the model are varied in pertinent ranges to physical problem simulated and the results demonstrated a good predictive capability.

O presente trabalho utiliza o método Constructal Design para desenvolver o estudo numérico da configuração de materiais de alta condutividade térmica em forma de “Y” que minimiza a resistência ao fluxo de calor, quando áreas ocupadas pelos materiais de alta e baixa condutividades são mantidas constantes. Para a solução numérica da equação diferencial da difusão do calor e suas respectivas condições de contorno, foi utilizado o software MATLAB ®, mais especificamente a ferramenta PDETOOL, Partial Differential Equations Tool. O objetivo deste trabalho é a minimização da resistência térmica do sistema gerador de calor com baixa condutividade térmica com a utilização de vias em formato de Y com material de alta condutividade térmica e volume constante, sendo variáveis os comprimentos e espessuras do material dos ramos simples e bifurcados. Todas as possibilidades geométricas foram avaliadas e a geometria ótima foi aquela que conduziu a menor resistência térmica. Duas condições são apresentadas, a primeira tem os ramos e a base da geometria “Y” com igual condutividade térmica. Os resultados para esta configuração mostram que existem valores específicos para os graus de liberdade que minimizam a resistência térmica. Nesse caso, os ramos se degeneraram e a configuração ótima tem a forma de um “V”. A segunda configuração apresenta combinações de condutividade térmica diferentes, para os ramos e a bases. Para estes casos obteve-se um valor otimizado próximo de 1 para a razão entre os comprimentos dos ramos simples e bifurcados, indicando que a configuração otimizada tem realmente a forma de um “Y” o que demonstra a dependência entre a geometria e as condições impostas pelo meio. Embora o design inicial do Y possa assumir diversas configurações, quando comparado o primeiro design com o design final, no caso do Y com iguais condutividades térmicas se conseguiu uma melhora superior a 28% e no caso do Y com condutividades diferentes mais de 30 %. Finalmente, este trabalho mostra que a geometria otimizada é aquela que melhor distribui as imperfeições, isto é, os pontos quentes (pontos de temperatura máxima), o que está de acordo com o princípio da ótima distribuição das imperfeições. / The present work used the method Constructal Design to develop numerical analyses of pathways of high thermal conductivity in "Y" shape which minimizes the thermal resistance when areas occupied by the materials of high and low conductivities are kept constant. For the numerical solution of the differential equations of heat diffusion and their boundary conditions, we used the MATLAB ® software, specifically the PDETOOL tool. The aim was to minimize the thermal resistance of the heat generator system with low thermal conductivity with the use of Y-shaped pathways with high thermal conductivity and constant volume, with variable lengths and thicknesses of material from stem and forked branches. All geometric possibilities were evaluated and the optimal geometry was that which resulted in lower thermal resistance. Two conditions were studied. In the first one the stem and branches of the "Y" have equal thermal conductivity. The results for this configuration show that there are specific values for the degrees of freedom to minimize the thermal resistance. In this case, the branches have degenerated and the optimum configuration has the shape of a "V". The second configuration offers different combinations of thermal conductivity, for branches and bases. For these cases we obtained a optimized value close to 1 for the ratio between the lengths of stem and bifurcated branches, indicating that the optimized configuration actually has the shape of a "Y" which shows the dependency of geometry and condition imposed by the environment. Although the initial design of Y can take various configurations, when compared the first design to the final design, in the case of Y with equal thermal conductivity it this improvement was achieved an improvement greater than 28% and in the case of Y with different conductivities over 30%. Finally, this study showed that the optimized geometry is the one that better distributes imperfections, this is, hot spots (points of maximum temperature), which is in accordance with the principle of the optimal distribution of imperfections.

As imagens são sinais que possuem muita informação. Os objetos representados em ima- gens podem sofrer deformações fazendo com que suas características mudem, o que difi- culta o reconhecimento do objeto. A chave para o sucesso de um sistema de reconheci- mento de padrões em imagens é escolher adequadamente a sua abordagem e um modelo para as feições presentes nas imagens. Uma das dificuldades é extrair e selecionar as feições que são mais discriminantes entre as diferentes classes usadas para modelar um objeto. Os modelos ativos de formas (Active Shape Models ASM) adaptam-se às defor- mações de um objeto. O objeto a ser modelado, pode ser representado com um modelo de pontos distribuídos (Point Distribution Models PDM). O PDM consiste de pontos de interesse ou marcadores, que permitem a extração de feições em localizações específicas do objeto. Após tratar a incerteza da localização e oclusão dos marcadores é possível ex- trair as feições mais representativas, obtendo-se um desempenho alto em termos da taxa de reconhecimento. Nesta tese são introduzidas novas formas para extrair e selecionar feições com modelos ativos de formas, que melhoram a taxa de classificação onde há objetos deformáveis. Esta tese é inovadora no sentido de aperfeiçoar o uso de ASMs na classificação de faces humanas, e na sua aplicação no monitoramento visual de outros tipos de objetos deformáveis. / Images are signals that have a lot of information. The objects depicted in images may suf- fer deformations causing changes in their characteristics, which hinders the recognition of the object. The key to the success of a system of pattern recognition in images is to choose properly your approach and a model for the features in the images. One difficulty is to extract and select the features that are most discriminating between different classes used to model an object. The Active Shape Models (ASM) adapt to deformations of an object. The object to be modeled can be represented with a Points Distribution Model (PDM). The PDM consists of points of interest or landmarks that allow the extraction of features in specific locations of the object. After treating the uncertainty of the location and occlu- sion of the landmarks it is possible to extract the most representative features, obtaining a high performance in terms of recognition rate. This thesis introduces new ways to extract and select features with ASMs, which improve the classification rate where deformable objects are present. This thesis is innovative in the sense that improves the use of ASMs in the classification of human faces, and that can be applied in visual monitoring of other types of deformable objects.

objetivo principal do trabalho é determinar a capacidade última de tração (arrancamento) de uma fundação com solo tratado no topo a partir de uma modelagem numérica com o programa computacional Abaqus a partir dos ensaios feitos em solos tratados. Os modelos numéricos foram feitas a partir dos ensaios ao arrancamento em solos tratados por Ruver, 2011 e Girardello, 2014 para um aterro de areia (proveniente da cidade de Osório/RS) e areia com fibra (monofilamentos de polipropileno). Foi feito um modelo numérico inicial com as propriedades geotécnicas de ensaios feitos em laboratório e campo para depois fazer um ajuste do modelo com os ensaios ao arrancamento. Foi avaliado o modelo numérico de areia e areia-fibra para as tensões iniciais, deslocamentos inciais e a etapa geostatic porque não só é ter bons resultados nas figuras do ajuste, mas se ter um modelo numérico bem feito que apresente um comportamento próximo da realidade numericamente. Foram avaliados os benefícios gerados pelo aumento da capacidade de carga de fundações escavadas e reaterradas com areia fina e homogênea, areia-fibra, avaliando sua influência geotécnica (coesão, ângulo de atrito e módulo de elasticidade) a partir de uma análise paramétrica. Foi feito um método de dimensionamento e estimativa da capacidade ao arrancamento para um aterro com areia e areia-fibra a partir das metodologias internacionais e os modelos numéricos. / The main objective of the work is to determine the ultimate traction (pullout) capacity of a foundation with treated soil at the top from a numerical modeling with the Abaqus computational program from the tests on treated soils. The numerical models were made from the pullout tests on soils treated by Ruver, 2011 and Girardello, 2014 for a sand dump (from the city of Osório / RS) and sand with fiber (polypropylene monofilaments). An initial numerical model was made with the geotechnical properties of laboratory and field tests to later calibrate the model with the pullout tests. It was evaluated the numerical model of sand and sand-fiber for the initial stresses, displacements and the geostatic step because not only is it to have good results in the figures of the calibration, but it is necessary to have a well-designed numerical model that presents behavior close to reality numerically. It was evaluated the benefits generated by the increase of the load capacity of excavated foundations and recovered with fine sand and homogeneous, sand-fiber, evaluating its geotechnical influence (cohesion, friction angle and modulus of elasticity) from a parametric analysis. A method of dimensioning and estimation of the pullout capacity for a backfill with sand and sand-fiber was done from the international methodologies and the numerical models.

Este trabalho de pesquisa tem por objetivo apresentar e investigar a viabilidade de um método numérico que contempla o paralelismo no tempo. Este método numérico está associado a problemas de condição inicial e de contorno para equações diferenciais parciais (evolutivas). Diferentemente do método proposto neste trabalho, a maioria dos métodos numéricos associados a equações diferencias parciais evolutivas e tradicionalmente encontrados, contemplam apenas o paralelismo no espaço. Daí, a motivação em realizar o presente trabalho de pesquisa, buscando não somente um método com paralelismo no tempo mas, sobretudo, um método viável do ponto de vista computacional. Para isso, a implementação do esquema numérico proposto está por conta de um algoritmo paralelo escrito na linguagem C e que utiliza a biblioteca MPI. A análise dos resultados obtidos com os testes de desempenho revelam um método numérico escalável e que exige pouco nível de comunicação entre processadores. / This research aims to present and investigate the feasibility of a numerical method that considers the parallelism in time. This numerical method is associated with problems of initial and boundary conditions for (evolutionary) partial differential equations. Unlike the method proposed in this work, most of the numerical methods associated with evolutionary partial diferential equations and traditionally found include only parallelism in space. Hence, the motivation for carrying out the present research work,seeking not only a method with parallelism in time but,above all, a viable method. The implementation of this proposed computational parallel algorithm was written with the language C and uses the MPI library. The results obtained from performance tests reveal a scalable and numerical method that requires little level of communication amount processors