Este trabalho apresenta o resultado da aplicação do Método dos Volumes Finitos, adotado pelo software comercial Star-CCM+ na simulação para o regime permanente de escoamentos turbulentos incompressíveis e compressíveis sobre seções de aerofólios. Para o caso incompressível modelam-se seções do aerofólio NACA 0012 com ângulo de ataque zero. Para o caso compressível, uma seção do aerofólio supercrítico OAT15A em pequeno ângulo de ataque é modelada. Os domínios computacionais são discretizados por malhas não estruturadas de células poliédricas ou por malhas estruturadas de geração hiperbólica para diferentes topologias e parâmetros construtivos determinados pela estimativa de grandezas do fenômeno físico, como a altura da primeira camada de células quadrilaterais junto à parede. A qualidade e adequação dessas malhas para as simulações são verificadas por estudo de dependência quanto ao nível de refinamento e também quanto à posição do contorno onde o escoamento é livre no caso de escoamento compressível. Na metodologia de verificação, o índice de convergência de malha GCI e a ordem observada de convergência do método (dos Volumes Finitos) são obtidos para três níveis de refinamento com o propósito de selecionar uma malha de trabalho que concilie precisão e esforço computacional com os recursos disponíveis. As simulações são conduzidas para dois modelos de turbulência – o modelo Spalart-Allmaras e o modelo k-ω/SST. Os resultados obtidos pela aplicação desses modelos são interpretados sob o ponto de vista fenomenológico e comparados com os resultados experimentais disponíveis em literatura. / The Finite Volumes Method adopted by the commercial software Star-CCM+ is applied to the simulation of the steady state regime of incompressible and compressible turbulent flows over selected airfoil’s sections. The physical model used with the incompressible flow case is a NACA 0012 airfoil section at zero angle of attack. The ONERA’s OAT15A supercritical airfoil section at small angle of attack applies to the compressible flow case. The computational domains are discretized by structured and unstructured meshes for different topologies and far field configurations. The structured meshes are of the quadrilateral type with hyperbolic node distribution whilst the unstructured meshes use polyhedral cells. The grids are generated by applying a methodology where estimates of the flow variables are used as input for the grid’s constructive parameters like the near wall cell thickness. Grid dependency studies are carried out in order to verify the grid’s quality and suitability to represent the physical phenomena. The grid’s asymptotic convergence index GCI and its observed order of convergence are evaluated for three refinement levels and far field position for the compressible flow cases. The objective is to select the most suitable grid taking into account the accuracy requirements and the computational resources available. The one-equation Spalart-Allmaras turbulence model and the two-equation k-ω/SST turbulence models are used. The numerical results are discussed from the physical point of view and compared with the experimental ones available in literature.

Os conceitos de mecânica da fratura têm sido cada vez mais utilizados em projetos de engenharia, permitindo a determinação de tamanhos críticos de defeito ou limites de solicitação em função da presença de trincas. No entanto, para aplicar com sucesso esses conceitos em uma avaliação de integridade estrutural é de suma importância a determinação exata do comportamento a fratura do material. Materiais com comportamento elasto-plástico têm ocupado cada vez mais lugar de destaque em estruturas e componentes críticos de engenharia e a maneira mais utilizada para caracterização da tenacidade à fratura desse tipo de material é a curva J-R. Embora existam métodos normatizados para determinação dessas curvas, esses possuem muitas vezes inconvenientes e limitações que tem justificado o desenvolvimento de vários estudos com vista à determinação de métodos alternativos. Nesse trabalho foi investigada a hipótese de que, a partir do estudo do campo deformacional nas proximidades da ponta de uma trinca, é possível desenvolver uma metodologia para determinação da curva de resistência à fratura de um material com comportamento elasto-plástico. Para tanto, um sistema óptico de medição de deformação baseado na técnica de correlação de imagem digital foi empregado para monitorar a deformação na superfície de corpos de prova do tipo SE(B) durante o ensaio de tenacidade à fratura. O objeto de estudo foi uma junta do aço API X65 soldada por processo de soldagem a laser. Os resultados obtidos demonstram o sucesso da técnica utilizada para obtenção da evolução do campo deformacional durante o ensaio. Além disso, a metodologia empregada possibilitou a determinação de expressões matemáticas que permitem associar valores de deformação com valores de integral J e Δa. Essas expressões são válidas para qualquer uma das três distintas regiões da junta soldada e permitem determinar com sucesso a curva J-R a partir da evolução da deformação durante o ensaio de um único corpo de prova. / The use of fracture mechanics has become increasingly popular in damage tolerance analysis or during determination of load limits in engineering components and structures. Nevertheless, the success of applying the concepts of fracture mechanics to investigate structural integrity assessment is highly dependent upon the correct determination of fracture mechanics material properties. The accurate prediction of ductile fracture behavior plays an important role in structural integrity assessment of critical engineering structures, including nuclear reactors, petrochemical vessels and tanks, pipelines in oil and gas industries and aircraft structures. J-R curve has become the most popular material parameters in elastic-plastic fracture mechanics and has been applied widely in practical engineering. There are procedures to determine J-R curves and they are widely presented in books and international standards. However, these procedures sometimes present some practical limitations or restrictions which have prompted research effort to develop alternative methods. In this work the hypothesis that is possible to develop a methodology to determinate J-R curve to materials with elastic-plastic behavior from analysis of near crack tip strain field was verified. For this purpose, an optical strain measurement system (OSMS) was used in order to monitor strains on surface of SE(B) specimens during fracture toughness testing. The methodology was applied to a laser welded joint of API X65 steel. The measurements successfully evaluate the strain field during testing. Besides of this, the methodology proposed allowed the determination of mathematical expressions to associate strain values with values of J-integral and crack extension (Δa). These expressions are valid to characterize these fracture toughness parameters to any one of the three distinct regions of the weld joint and made possible to determine J-R curve by monitoring the surface strain field evolution during a testing of a single specimen.

Neste trabalho é apresentada uma análise da antena planar F-invertido (PIFA) com o objetivo de identificar os parâmetros geométricos que influenciam nas características de desempenho: frequência de ressonância, perda de retorno e largura de banda. Para realizar este estudo é desenvolvido um algoritmo em linguagem C baseado no método das Diferenças Finitas no Domínio do Tempo (FDTD). Duas PIFAs de geometria distintas foram investigadas. O primeiro modelo trata da geometria tradicional de elemento irradiador e plano de terra, ambos retangulares. Já o segundo modelo apresenta a fenda em formato de L no elemento irradiador, para operação em frequência dupla, e o plano de terra em formato de T, visando a aumentar a largura de banda. Os resultados da simulação do algoritmo identificam os parâmetros que influenciam no aumento nas duas larguras de banda, e são aplicados na prototipação de uma PIFA GSM-900 e GSM-1800. As medidas de laboratório do protótipo evidenciam os resultados previstos pelo algoritmo. / A Planar Inverted-F Antenna (PIFA) analysis is presented in this work, aiming to identify the geometric parameters that influence the antenna performance, such as resonant frequency, return loss and bandwidth. To accomplish this goal, an algorithm in C language based on the Finite-Difference Time-Domain (FDTD) method is developed. Two PIFAs with distinct geometries are investigated. The first one corresponds to a typicall rectangular patch and ground plane, The second model presents an L-shaped slot for dual frequency band operation, combined with a T-shaped ground plane, in order to increase the bandwidth. The simulated results obtained by algorithm identify the parameters that respond for both bandwidths enhancement, and are applied in a PIFA GSM-900/GSM-1800 practical project. The prototype’s measurements confirm the simulated results.

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.

Sistemas de Distribuição de Energia Elétrica (SDEE) estão, pelas próprias condições de funcionamento, sujeitos a perdas de energia e quedas nos níveis de tensão. Uma das formas de reduzir as perdas e ao mesmo tempo melhorar o perfil de tensão do SDEE consiste na utilização de bancos de capacitores e reguladores de tensão ao longo dos alimentadores. Para isto a alocação ótima destes equipamentos torna-se de fundamental importância, logo este trabalho tem por objetivo a utilização de um método quadrático de programação por metas multiobjetivo para a ótima alocação desses equipamentos visando a redução de perdas (meta 1) com o menor custo de investimento (meta 2) sujeito a restrições do fluxo de potência, limites técnicos de operação, mudança no controle dos capacitores chaveados e dos taps dos reguladores de tensão. O constante aumento na penetração da geração distribuída traz novos desafios para o controle e a operação dos SDEE. A metodologia proposta considera os efeitos destas fontes de geração e os seus impactos na alocações ótima dos bancos de capacitores e reguladores de tensão. Com o intuito de resolver o problema de otimização em programas comerciais como Gurobi e Cplex uma metodologia linear de cálculo de fluxo de potência é utilizada. Os algoritmos foram desenvolvidos no software Matlab, e posteriormente foi realizada a conversão para linguagem GAMS (Generic Algebraic Modeling System) para utilização dos solvers na plataforma NEOS (Network Enabled Optimization System). / Power distribution systems (DS) are susceptible to energy loss and voltage drops, mainly because of the way that they are operated and built. One way to minimize system loss and improve the voltage profile of DS is with the use of capacitor banks and voltage regulators (VRs) along the distribution feeders. In this scenario, the optimal placement of these devices becomes crucial. In this paper, a quadratic multi objective goal programming method is proposed to obtain the optimal placement of these devices, where the goals are energy loss (goal 1) and the total investment cost (goal 2). The model is subject to system operation constraints, such as power flow demands, limits of voltages and currents, number of stages of switched capacitors and the tap control of VRs. The increase of distributed generation penetration brings new challenges for control and operation of DS. The proposed methodology considers the effects of these power sources and their impact on optimal allocations of capacitors and VRs. With the aim of solving the optimization problem in commercial programs, a linear power flow is used. The algorithms were developed with Matlab software and the conversion was made to a General Algebraic Modeling System (GAMS) language and solved in the NEOS (Network Enabled Optimization System) platform.

A velocidade de queima laminar adiabática é um importante parâmetro da combustão que dita o comportamento de chamas pré-misturadas. Dos métodos disponíveis para a medição desse parâmetro, o método do fluxo de calor destaca-se pela simplicidade e precisão. No presente trabalho, esse método é utilizado para medir a velocidade de queima de biogás (modelado como CH4 com diferentes níveis de diluição com CO2) e de gás de síntese (modelado como uma mistura de CH4, H2, CO, CO2 e N2) em ar a 298 K e 1 atm. Tais gases são de crescente interesse para a sociedade em função de aspectos ambientais, porém, suas velocidades de queima não foram amplamente estudadas ainda. Os resultados obtidos são comparados com as previsões de cinco mecanismos cinéticos (GRI-Mech 3.0, Davis et al., Konnov, San Diego e USC Mech II) a fim de avaliar a sua capacidade preditiva. Os resultados experimentais e numéricos das velocidades de queima de biogás e ar apresentam uma boa concordância e as incertezas encontradas foram condizentes com as relatadas na literatura. Os resultados experimentais desse gás foram parametrizados em uma correlação empírica de fácil utilização em modelos numéricos. As medições da velocidade de queima de gás de síntese e ar, por outro lado, apresentaram valores inferiores às previsões numéricas de todos os mecanismos estudados. Os dados experimentais da literatura, para a mesma mistura, diferem tanto em valores quanto em comportamento dos resultados do presente trabalho. Tal comportamento está provavelmente relacionado a alguma contaminação no CO utilizado, já que quando esse gás está presente observa-se uma chemi-luminescência não relatada na literatura. / The adiabatic laminar burning velocity is an important combustion parameter that dictates premixed flames characteristics. Among the measuring methods available in literature, the heat flux method stands out for its simplicity and accuracy. In the present work, this method is used to measure the adiabatic laminar burning velocity of biogas (modeled as CH4 with different dilution levels with CO2) and syngas (modeled as a CH4, H2, CO, CO2 and N2 mixture) in air at 298 K and 1 atm. Such gases are of growing society interest due to environmental aspects, however, their adiabatic laminar burning velocity have not been widely studied yet. The experimental results are compared to predictions of five kinetic mechanisms (GRI-Mech 3.0, Davis et al., Konnov, San Diego e USC Mech II) to evaluate their predictive capacity. Experimental and numerical results of biogas/air mixtures adiabatic laminar burning velocity show good agreement and the found uncertainties are in agreement with literature. Experimental results of this gas were fitted in an empiric correlation of simple numerical application. Experimental results of the laminar burning velocity of syngas/air, on the other hand, show lower values than the numerical predictions of all studied kinetic mechanisms. Literature available data for the same mixture differ both in values and behavior of the present work results. Such behavior is probably related to some contamination on the CO used since a chemi-luminescence not reported in literature can be noted when this gas is present.

Neste trabalho, é estudada a descrição do fluxo de um gás rarefeito em um duto cilíndrico de comprimento infinito. A formulação matemática do problema está baseada na forma integral de equações cinéticas derivadas da Equação de Boltzmann. Particularmente são estudados os modelos cinéticos conhecidos como BGK e S. Métodos espectrais são propostos para obtenção de soluções, em forma fechada, para quantidades de interesse como o perfil de velocidade do gás, bem como taxas de fluxo. As formulações espectrais são baseadas em duas abordagens: expansão clássica em termos de Polinômios de Legendre e expansão em termos de splines cúbicas de Hermite, neste caso, associada a um esquema de colocação. A implementação das propostas produz resultados computacionais satisfatórios do ponto de vista prático. Para obtenção de resultados com maior precisão, técnicas de tratamento da singularidade do núcleo da equação integral foram introduzidas, resultando em ganho computacional significativo. Finalmente, a proposta de solução espectral para problemas em geometria cilíndrica se mostrou adequada para problemas em que se admite reflexão especular na superfície do cilindro, situação onde outras abordagens clássicas disponíveis na literatura não podem ser utilizadas. / In this work, rarefied gas flows in cylindrical ducts are studied. The mathematical formulation of the problems are based on the integral form of kinetic equations derived from the Boltzmann equation. Particularly, the BGK and S models are studied. Spectral methods are proposed to obtain closed form solutions for quantities of interest as velocity profile of the gas as well as flow rates. The spectral formulations are based on two approaches: classical expansions in terms of Legendre Polynomials and Hermite cubic splines expansions. In this case, associated with a collocation scheme. The approaches provide good computational results, from the practical point of view. On the other hand, for obtaining higher accuracy, some techniques were introduced to deal with the inherent singularity of the integral kernel. In this context, a significant computational gain is achieved. Finally, this spectral approach has shown to be adequate to solve problems where specular reflection is assumed at the surface, in which cases, classical approaches available in the literature can not be used.

Neste trabalho, são inicialmente discutidas as dificuldades referentes à obtenção de resultados numéricos para a operação de uma turbina Savonius independentes do grau de discretização, tamanho de domínio de cálculo e de máximo tempo físico simulado. Também são relatadas as divergências entre as metodologias numéricas e experimentais adotadas por diversos autores, que dificultam análises e comparações dos resultados obtidos por meio dessas metodologias com os resultados próprios obtidos. Devido a esses fatos, no presente trabalho, uma série de procedimentos experimentais e numéricos são realizados para efetuar análises do escoamento sobre uma turbina eólica Savonius. Nos experimentos em canal aerodinâmico, perfis de velocidade e parâmetros da turbulência são obtidos pela técnica de anemometria de fio quente. Medições com o uso de tubos de Pitot e manômetros eletrônicos são efetuadas para avaliar a variação da pressão e os perfis de velocidade média em posições selecionadas. Além de dados para análise, informações úteis para uso como condições de contorno nas simulações numéricas também são obtidas. Os fenômenos são reproduzidos através de simulações numéricas pelo Método de Volumes Finitos, que solucionam as equações da continuidade, de Navier-Stokes com médias de Reynolds e do modelo de turbulência k-ω SST. Análises experimentais e numéricas considerando o escoamento sobre um cilindro, que mantém semelhanças com o escoamento sobre o rotor, também são realizadas. Simulações numéricas do escoamento sobre o cilindro são efetuadas, fornecendo resultados representativos do escoamento real, quando geometrias tridimensionais são aplicadas na modelagem numérica. Nas simulações do escoamento sobre o rotor Savonius em condição estática, resultados representativos do escoamento real são obtidos com o uso de uma modelagem que leva em consideração a rugosidade das pás do rotor, estacionado na posição angular de 90°. Para posições angulares menores, não se obteve uma boa concordância entre os resultados experimentais e numéricos. A realização deste trabalho fornece informações úteis para a análise do fenômeno e tem potencial para contribuir com futuros trabalhos desse tema. / This research work initially presents a discussion about the difficulties related to obtaining numerical results for the operation of a turbine Savonius independent of the degree of discretization, calculation domain size and maximum physical time of the simulation. The differences between the numerical and experimental methodologies adopted by various authors difficult the analysis and comparisons of the results obtained through these methods with the results obtained by the methodology. Due to these facts, in this research work, a series of experimental and numerical procedures are performed to conduct analyzes of flow over a Savonius wind turbine. In the experiments on aerodynamic channel, velocity and turbulence profiles parameters are obtained by the technique of hot wire anemometry. Measurements using Pitot tubes and electronic manometers are made to evaluate the variation of pressure and mean velocity profiles at selected positions. In addition to data analysis, useful information for use as boundary conditions in the numerical simulations are also obtained. The phenomena are reproduced through numerical simulations by Finite Volume Method, that solve the equations of continuity, Reynolds-averaged Navier–Stokes equations and the equation of the turbulence model k-ω SST. Experimental and numerical analyzes considering the flow over a cylinder, which holds similarities with the flow over the rotor, are also performed. Numerical simulations of the flow over the cylinder are made, providing results representative of the actual flow when three-dimensional geometries are applied in numerical modeling. In flow simulations over the Savonius rotor in static condition, at 90°, representative results of the actual flow are obtained using a model that takes into account the roughness of the rotor blades. For smaller angular positions a good agreement between experimental and numerical results was not obtained. This work provides useful information for the analysis of the phenomenon and has the potential to contribute to future work on this theme.

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.