A intransitividade é uma propriedade dos grafos conectados e orientados que representam interações entre espécies. Essa propriedade está associada `a persistência da coexistência mesmo na presença de competição, sendo o jogo Pedra-Papel-Tesoura o exemplo padrão. Neste trabalho, consideramos uma generalização com quatro espécies, que é o número mínimo necessário para que haja outras interações além do ciclo único (um predador, uma presa). Além disso, introduzimos sítios vazios dinâmicos ao sistema, os quais afetam o nível de intransitividade nas interações. Verificamos então que, ao contrário do que o campo médio prevê, na rede quadrada o modelo apresenta duas transições, cujas localizações dependem das probabilidades de predação e reprodução. A primeira transição ocorre entre um estado absorvente, com apenas uma espécie, para um estado em que todas as espécies coexistem. A transição seguinte ocorre quando uma das quatro estratégias é extinta. Essa dependência com as probabilidades de predação e reprodução demonstra que a estrutura do grafo de interação sozinha não é suficiente para prever os resultados finais em modelos desse tipo. Adicionalmente, probabilidades diferentes de predação permitem ajustar o nível de transitividade do sistema, indicando que é necessária uma quantidade mínima de intransitividade para que a coexistência entre todas as espécies persista. / Intransitivity is a property of connected, oriented graphs representing species interactions that may drive their coexistence even in the presence of competition, the standard example being the three species Rock-Paper-Scissors game. We consider here a generalization with four species, the minimum number of species allowing other interaction beyond the single loop (one predator, one prey). We also consider dynamic empty sites, which affect the amount of intransitivity in the system. We show that contrary to the mean field prediction, in the square lattice the model presents two transitions as the parameters setting the predation and reproduction rates change, from an absorbing state with one species only, to a coexistent state first, and then to a state in which one species gets extinct. Such a dependence on the predation and reproduction rates shows that the interaction graph structure alone is not enough to predict the outcome of such models. In addition, different predation rates permit to tune the level of transitiveness, indicating that for the coexistence of all species to persist, there must be a minimum amount of intransitivity.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Um aspecto particular no dimensionamento de maciços reforçados com geossintéticos consiste na análise da estabilidade interna. A ruptura interna pode ocorrer quando as solicitações impostas ao elemento de reforço superam a resistência à tração, ou quando ocorre o arrancamento do reforço da massa de solo, por ancoragem insuficiente. A distribuição das deformações e das tensões ao longo do comprimento enterrado do reforço não é uniforme, e este aspecto não é considerado no dimensionamento. Desta forma, modelos analíticos e numéricos aparecem como alternativas capazes de reduzir incertezas no dimensionamento de maciços reforçados, permitindo a adoção de soluções menos conservativas. O presente trabalho propõe um modelo analítico para a reprodução do mecanismo de transferência de esforços e deslocamentos ao longo do comprimento de geotêxteis sob condição de arrancamento, e apresenta a modelagem numérica de ensaios de arrancamento, fazendo uso do programa Plaxis, de elementos finitos. A partir dos resultados de um extenso programa experimental de ensaios de arrancamento instrumentados, em geotêxteis (Espinoza,2000), os modelos analítico e numérico foram validados e discutidos. Posteriormente, apresenta-se a simulação de um ensaio de arrancamento hipotético fazendo uso de ambos os modelos. Os resultados sugerem que os modelos analítico e numérico foram adequados na previsão dos esforços, deformações e deslocamentos ao longo do comprimento de geotêxteis em solicitações de arrancamento. Observou-se um melhor ajuste entre as previsões do modelo analítico e os resultados experimentais, justificado pela adoção do modelo não linear para o elemento geotêxtil. A distribuição de esforços e deslocamentos ao longo de geotêxteis é complexa, e a boa concordância dos modelos com os resultados experimentais reforça a potencialidade dos modelos para uso futuro. / The internal stability analysis is a particular aspect in the design of reinforced soil with geotextiles. Internal failure may occur when the stresses transmitted to the reinforcing element exceeds the tensile strength or when it is pulled out, due to insufficient anchorage. The stress strain distribution along the embedded length of the reinforcement is not uniform, and this condition is not incorporated in the design. Thus, analytical and numerical models appear as alternatives to reduce uncertainties in the design of reinforced soil structures, allowing the adoption of less conservative solutions.The present research proposes an analytical model that reproduces load transfer mechanism and displacements along the length of geotextiles under pullout condiction, and also presents a numerical simulation of pullout tests, making use of Plaxis FEM program. Based on the results of an extensive program of pullout tests in instrumented samples of geotextiles (Espinoza, 2000), the analytical and numerical models were validated and discussed. The simulation of a hypothetical pullout test making use of both models is also presented.The results suggest that the analytical and numerical models are suitable to predict loads, strains and displacements along the geotextiles length, submitted to pullout. The analytical model provided a better fit for the experimental results, since it incorporates a non-linear behavior for the geotextile. The distribution of loads and displacements along the geotextiles is complex, and the good agreement between the models and the experimental results emphasizes the capability of the models for further use.

Com os grandes avanços ocorridos na disponibilização de computadores, existe uma tendência contínua para a utilização de técnicas computacionais auxiliando no projeto de equipamentos de engenharia. Cada vez mais estão se obtendo resultados bastante próximos às condições reais, incluindo a simulação de motores de combustão interna. Neste sentido o presente trabalho tem o objetivo de analisar o escoamento turbulento no processo de admissão de ar em um motor operando em ciclo Diesel. A investigação é focada na determinação da influência do passo de tempo no cálculo do coeficiente de descarga e razão de swirl. Adicionalmente, o campo de velocidades, pressão, energia cinética turbulenta e outros parâmetros são apresentados e analisados, com o objetivo de auxiliar no entendimento da dinâmica envolvida. Essencialmente, dois modelos de turbulência são empregados, juntamente com dois tratamentos de parede. Seus resultados também são confrontados e discutidos. A geometria considerada é de um motor Fiat 1.9 L quatro tempos com duas válvulas. A análise é concentrada em um único cilindro. O pacote computacional utilizado é o Star-cd, e seu aplicativo es-ice. A independência de malha foi obtida, chegando a 1.672.056 volumes. Os resultados são apresentados de duas formas. A primeira delas refere-se a resultados de simulações em regime permanente, realizadas em boa parte por outros autores, com ênfase na determinação do coeficiente de descarga e razão de swirl, estes confrontados com valores experimentais, visando à validação da metodologia. Fica evidente a importância da escolha do modelo de turbulência na simulação de motores de combustão interna, assim como das funções de interpolação utilizadas. Na segunda parte os resultados referem-se a uma análise transiente, considerando o movimento do pistão e válvulas, a 1500 RPM. Observa-se a grande exigência quanto ao passo de tempo requerido no transiente real, ficando demonstrado que para esta velocidade o menor passo de tempo utilizado, 0,05° (5.5555E-6 s), ainda é insuficiente para alguns momentos do ciclo. É possível notar maior influência no coeficiente de descarga do que na razão de swirl, em relação aos passos de tempo utilizados. A forte dependência do modelo de turbulência nos resultados obtidos é mais uma vez confirmada, conforme o esperado, já que as hipóteses sobre a física do fenômeno são diferentes em cada modelo. Os resultados quanto ao tratamento na parede não apresentaram significantes diferenças, quando aplicados junto ao modelo de turbulência k-ω SST. / Considering the increase in the availability of computers, there is a continuing trend toward the use of computational simulation aiding in the design of engineering equipments. Reasonable results, close to the real conditions, are obtained, including the simulation of internal combustion engines. In this way, the present work has the objective of analyzing the turbulent flow in the air intake process of an engine operating in Diesel cycle. The investigation focuses on the determination of the time step in the calculation of the air discharge coefficient and swirl ratio. Additionally, the turbulent kinetic energy, pressure and velocity fields, besides other parameters, are presented and analyzed, with the objective of aiding in the understanding of the involved dynamics. Essentially, two turbulence models are employed, together with two wall treatments. Their results are also confronted and discussed. The considered geometry is a four-stroke, 1.9-L FIAT engine, with two valves. The analysis is concentrated on a single cylinder. The software package used is the Star-cd, and its application es-ice. The mesh independence is carried out, arriving in 1.672.056 volumes. The results are presented in two ways. The first one refers to simulation results of the steady state, also accomplished by other authors, with emphasis in the determination of the discharge coefficient and swirl ratio. These data are confronted with experimental values, aiming to validate the applied methodology. The importance of the choice of the turbulent model becomes evident in the simulation of internal combustion engines, as well as the interpolation functions used. In the second part the results refer to a transient analysis, considering the valves and piston movement, at 1500 rpm. It is observed the great demand on time step required is observed for the real transient, demonstrating that, for this speed, the smallest time step used, 0.05º (5.5555E-6 s), is still insufficient for some moments of the cycle. Also regarding the time step, it is possible to notice a greater influence in the discharge coefficient than in the swirl ratio. The strong dependence of the turbulence model on the results is once again confirmed, as expected, since the hypotheses about the physics of the phenomenon are different in each model. The results, regarding the wall treatment, presented no significant differences, when applied together with the SST k-ω turbulence model.

Made available in DSpace on 2014-06-11T19:32:09Z (GMT). No. of bitstreams: 0 Previous issue date: 2011-06-30Bitstream added on 2014-06-13T19:02:30Z : No. of bitstreams: 1 mitmapillaca_ejd_dr_guara.pdf: 1950757 bytes, checksum: f8081b468d019f4dac4207cfc646916a (MD5) / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Implantação iônica por imersão em plasma com campo magnético (3IPCM) foi investigada usando técnicas numéricas e experimentais. O campo magnético considerado é essencialmente não uniforme e é produzido por duas bobinas magnéticas posicionadas ao redor da câmara de vácuo. O estudo é centrado na análise do efeito de dois dos parâmetros mais importantes: tensão e pressão no processo 3IPCM. Outro tema importante como a dinâmica dos elétrons secundários foi também abordado neste trabalho. Neste contexto, o processo 3IPCM foi pesquisado inicialmente usando o código computacional KARAT. Os resultados numéricos mostraram um incremento da densidade do plasma ao redor do alvo durante a variação dos parâmetros de tensão, pressão e campo magnético considerados. Como consequência deste aumento, um incremento da densidade de corrente iônica sobre o alvo foi observado. Os resultados numéricos mostraram que o sistema de campos cruzados E×B intensifica o processo 3IPCM. Posteriormente, 3IPCM foi realizado experimentalmente. Resultados experimentais mostraram que a densidade de corrente foi incrementada em aproximadamente 100 % em relação ao caso sem campo magnético quando os parâmetros externos foram variados. Todos estes resultados numéricos e experimentais são explicados através do mecanismo de ionização do gás por colisão com os elétrons magnetizados realizando movimento de deriva em campos E×B. Finalmente, para analisar os efeitos do processo 3IPCM no tratamento de materiais foram realizados implantações em amostras de silício. Os resultados mostraram que o processo 3IPCM promove mudanças nas propriedades superficiais das amostras, tornando-as hidrofóbicas. Esta técnica mostra ser atrativa posto que foi possível incrementar a dose e a profundidade de implantação em alta tensão. / Plasma immersion ion implantation (PIII) with magnetic field has been investigated using numerical and experimental methods. The magnetic field in consideration is essentially non-uniform and is generated by two magnetic coils installed outside the PIII vacuum chamber. The study is focused on analysis of the effect of two of the most important process parameters: voltage and gas pressure on the PIII with magnetic field. Another important subject such as the dynamics of secondary electrons has also been considered in this work. In this context, the PIII process with magnetic field has been initially analyzed numerically using the 2.5D computer code KARAT. The numeric results have shown an increase of the plasma density around of the target in the range of the considered parameters, voltage, pressure and magnetic field. As consequence of this an enhancement of the ion current density on the target was observed. The simulation results have demonstrated that the system of crossed E×B fields intensifies the PIII process with magnetic field. Later, the PIII process with magnetic field has been carried out experimentally. Experimental results have shown an increase of the current density in about 100 % in relation to the case without magnetic field when the external parameters have been varied. The numerical and experimental results are explained through the mechanism of gas ionization by collision with electrons drifting in crossed E×B field. Finally, to analyze the effect of the PIII process with magnetic field in material treatment implantation in Silicon samples has been carried out. The results indicate that the PIII process with magnetic field promotes changes of the samples surface properties, turning them hydrophobic. This PIII technique is attractive since it can increase the dose and the depth of implantation at high voltage.

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.

O objetivo deste trabalho é o cálculo modal da resposta impulso distribu ída de uma viga descrita pela equação de Timoshenko e das vibrações forçadas, devidas a influência de cargas externas. Os modos vibratórios foram obtidos com o uso da base dinâmica, gerada por uma resposta livre e suas derivadas. Esta resposta é caracterizada por condições iniciais impulsivas. Simulações foram realizadas para os modos, a resposta impulso distribuída e vibrações forçadas em vigas apoiadas em uma extremidade e na outra livre, fixa, deslizante ou apoiada, sujeitas a cargas oscilatórias espacialmente concentradas ou distribuídas através de pulsos.

Este trabalho tem como objetivo investigar quais são os valores dos principais parâmetros que a literatura apresenta para o processo de forjamento a morno de precisão. Os parâmetros escolhidos são: curva de escoamento do SAE-1050M a frio e a quente, fator de atrito e o coeficiente de transferência de calor entre peça e ferramenta. Estipulou-se uma faixa de variação destes parâmetros, simulou o processo através do método dos elementos finitos com o uso do pacote MSC.Superform 2005 para os diferentes parâmetros e cruzou os resultados por meio do conceito de sensitividade. Os resultados cruzados são: força, trabalho, deformação e temperatura. Simulações são realizadas com variações dos parâmetros de entrada baseadas nos ensaios e levantamento bibliográficos. As respostas das simulações são adimensionalizadas com o uso do conceito de sensitividade, podendo assim, serem comparadas entre si. O fator de atrito e o coeficiente de transferência de calor entre peça e ferramenta não influenciam tanto na resposta como a curva de escoamento. / This work has as objective to investigate which values of key parameters are presented by the literature for the warm precision forging process. The chosen parameters are: hot and cold SAE-1050M flow stress curve, friction factor and heat transfer coefficient between workpiece and die. Stipulate a range of variation of these parameters, simulate the process through the finite elements method with the use of the package MSC.Superform 2005 with different parameters and cross the results with concept of sensitivity. The crossed results are: force, work, strain and temperature. Simulations are carried out with variations of the input parameters based on the bibliographical research and tests. The simulations output are nondimensionalised with the use of the sensitivity concept, thus being able, to be comparative between themselves. The friction factor and the heat transfer coefficient between workpiece and die do not influence as much as the flow stress curve.

Neste trabalho é apresentado um estudo numérico sobre escoamentos incompressíveis, não isotérmicos, bi e tridimensionais nos regimes laminar e turbulento através da Simulação de Grandes Escalas e da utilização do Método de Elementos Finitos. Para tornar isso possível, é implementada a equação da energia e os termos de forças de campo (empuxo) em um algoritmo numérico desenvolvido em FORTRAN, já existente, que simula escoamentos incompressíveis, isotérmicos, tridimensionais, nos regimes laminar e turbulento. O código desenvolvido abrange escoamentos onde as formas básicas de troca térmica ocorrem por difusão e advecção. No que tange a natureza da convecção térmica é possível analisar escoamentos com convecção forçada, mista ou natural. O método numérico empregado é o de elementos finitos (FEM) e a discretização espacial das equações que governam o fenômeno (continuidade, conservação da quantidade de movimento e conservação da energia) é realizada através do método de Galerkin. Para a análise dos termos temporais nos escoamentos transientes aplica-se o esquema temporal explícito de Taylor-Galerkin. O elemento finito utilizado é o hexaedro isoparamétrico de oito nós. É empregado o método da pseudo-compressibilidade com o objetivo de manter os termos derivados da pressão na equação da continuidade, pois essa ausência gera uma dificuldade adicional na discretização das equações. Para a abordagem da turbulência é empregada a simulação de grandes escalas (LES) com modelagem sub-malha clássica de Smagorinsky para a viscosidade e a difusividade turbulenta. Visando a melhoria no tempo de processamento foi utilizada integração explícita das matrizes dos elementos e a técnica de processamento paralelo OpenMP. São apresentados resultados para escoamentos com vários números de Reynolds, Prandtl e de Grashoff dos campos de velocidade, pressão e temperatura para escoamentos em cavidade bidimensional, nos regimes laminar e turbulento, e para o degrau tri-dimensional no regime laminar. As simulações para escoamentos em cavidades nos regimes laminar e na região de transição são comparados com os resultados de outros autores, se mostrando bastante satisfatórios, tanto no regime transiente como no permanente. Além disso, a inserção das forças de campo no código melhorou os resultados obtidos com o mesmo. As outras simulações são apresentadas como novos casos e tiveram um comportamento qualitativamente satisfatório. / A numerical study about non-isothermal, bi and three-dimensional, incompressible, laminar and turbulent flows is done in this work using Large Eddy Simulation and Finite Element Method. To became this possible, is implemented the energy equation and buoyance forces (in the Navier-Stokes equations) in a numerical algorithm, developed in FORTRAN, already existent, that simulate isothermal, three-dimensional, incompressible, laminar and turbulent flows. The developed code includes flows where the basic forms of heat transfer are diffusion or advection. About the nature of thermal convection it is possible to analyze the forced, mixed or natural convection flows. The numerical method used is the finite element method (FEM) and the spatial discretization of governing equations of phenomena (mass, conservation of momentum and conservation of energy) is done through Galerkin method. To analyze the time-dependent terms in transient flows is employed a time-explicit Taylor-Galerkin scheme. The finite element used is the isoparametric hexahedral with eight nodes. It is used the pseudo-compressibility method to keep the pressure terms in continuity equation, because without these terms there are additional difficulties to obtain the discretizated equations. Regarding the turbulence approach, it is employed the large eddy simulation (LES) and for subgrid-scales is used the classical Smagorinsky model to turbulent viscosity and diffusivity. To minimize the processing time is used explicit integration of element matrix and the multiprocessing technical OpenMP. Results are presented to a wide range of Reynolds, Prandtl and Grashoff numbers for velocity, pressure and temperature fields to laminar and turbulent, bi-dimensional, lid-driven cavity flow and a laminar three-dimensional backward-facing step. Simulations of lid-driven cavity flows in laminar and transitional regimes is compared with others authors results, presenting good agreement, in both transient and permanent regimes. Besides that, the implementation of buoyance forces in the present code improved the results obtained by it. The others simulations are presented like new cases and had qualitatively good behavior.

A avaliação numérica das cargas limites (de colapso) de estruturas é obtida geralmente ou pela aplicação direta dos teoremas de análise limite junto com processos de otimização ou mediante análises incrementais, que levam em consideração o comportamento elasto-plástico do material. Entretanto, ambas as estratégias conduzem eventualmente a dificuldades numéricas, particularmente para cargas próximas à de colapso. Neste trabalho, emprega-se uma alternativa que consiste em simular assintoticamente o comportamento elasto-plástico mediante uma relação elástica não-linear. As vantagens deste tipo de formulação são a possibilidade de se lidar com leis de fluxo não-associadas e um custo computacional reduzido. A relação elástica não-linear é implementada no programa comercial de elementos finitos ABAQUS, através de uma sub-rotina externa ao programa escrita em linguagem FORTRAN. Diversos exemplos de estruturas cujos materiais são regidos pelos critérios de resistência de von Mises e Drucker-Prager são modelados, verificando-se que os resultados das cargas limites obtidas com essa formulação são muito próximos daqueles encontrados na literatura. Por fim, a relação elástica não-linear é empregada para a determinação do domínio de resistência de meios porosos, com diferentes níveis de porosidade. / The numerical assessment of limit loads of structures is generally achieved through the direct implementation of limit analysis theorems together with optimization processes, or through incremental analyses, which account for the elastic-plastic behavior of the material. However, both the strategies may lead to numerical difficulties, particularly when the load is close to its limit value. In this context, the alternative approach presented in this work consists in simulating asymptotically the regime of free plastic flow by means of a fictitious non-linear elastic material. One of the main advantages of this kind of formulation lies in its ability to deal with non-associated flow rules and a reduced computational cost. The non-linear elastic behavior is implemented into the finite element computational software ABAQUS, making use of an external subroutine written in FORTRAN language. Several examples of geotechnical and structural problems with materials ruled by von Mises and Drucker-Prager failure criteria are analyzed. The results obtained with this formulation prove to be very close to those obtained through analytical solutions. At last, the non-linear elastic relation is used in the determination of the resistance domain of porous media with different levels of porosity.

Com os grandes avanços ocorridos na disponibilização de computadores, existe uma tendência contínua para a utilização de técnicas computacionais auxiliando no projeto de equipamentos de engenharia. Cada vez mais estão se obtendo resultados bastante próximos às condições reais, incluindo a simulação de motores de combustão interna. Neste sentido o presente trabalho tem o objetivo de analisar o escoamento turbulento no processo de admissão de ar em um motor operando em ciclo Diesel. A investigação é focada na determinação da influência do passo de tempo no cálculo do coeficiente de descarga e razão de swirl. Adicionalmente, o campo de velocidades, pressão, energia cinética turbulenta e outros parâmetros são apresentados e analisados, com o objetivo de auxiliar no entendimento da dinâmica envolvida. Essencialmente, dois modelos de turbulência são empregados, juntamente com dois tratamentos de parede. Seus resultados também são confrontados e discutidos. A geometria considerada é de um motor Fiat 1.9 L quatro tempos com duas válvulas. A análise é concentrada em um único cilindro. O pacote computacional utilizado é o Star-cd, e seu aplicativo es-ice. A independência de malha foi obtida, chegando a 1.672.056 volumes. Os resultados são apresentados de duas formas. A primeira delas refere-se a resultados de simulações em regime permanente, realizadas em boa parte por outros autores, com ênfase na determinação do coeficiente de descarga e razão de swirl, estes confrontados com valores experimentais, visando à validação da metodologia. Fica evidente a importância da escolha do modelo de turbulência na simulação de motores de combustão interna, assim como das funções de interpolação utilizadas. Na segunda parte os resultados referem-se a uma análise transiente, considerando o movimento do pistão e válvulas, a 1500 RPM. Observa-se a grande exigência quanto ao passo de tempo requerido no transiente real, ficando demonstrado que para esta velocidade o menor passo de tempo utilizado, 0,05° (5.5555E-6 s), ainda é insuficiente para alguns momentos do ciclo. É possível notar maior influência no coeficiente de descarga do que na razão de swirl, em relação aos passos de tempo utilizados. A forte dependência do modelo de turbulência nos resultados obtidos é mais uma vez confirmada, conforme o esperado, já que as hipóteses sobre a física do fenômeno são diferentes em cada modelo. Os resultados quanto ao tratamento na parede não apresentaram significantes diferenças, quando aplicados junto ao modelo de turbulência k-ω SST. / Considering the increase in the availability of computers, there is a continuing trend toward the use of computational simulation aiding in the design of engineering equipments. Reasonable results, close to the real conditions, are obtained, including the simulation of internal combustion engines. In this way, the present work has the objective of analyzing the turbulent flow in the air intake process of an engine operating in Diesel cycle. The investigation focuses on the determination of the time step in the calculation of the air discharge coefficient and swirl ratio. Additionally, the turbulent kinetic energy, pressure and velocity fields, besides other parameters, are presented and analyzed, with the objective of aiding in the understanding of the involved dynamics. Essentially, two turbulence models are employed, together with two wall treatments. Their results are also confronted and discussed. The considered geometry is a four-stroke, 1.9-L FIAT engine, with two valves. The analysis is concentrated on a single cylinder. The software package used is the Star-cd, and its application es-ice. The mesh independence is carried out, arriving in 1.672.056 volumes. The results are presented in two ways. The first one refers to simulation results of the steady state, also accomplished by other authors, with emphasis in the determination of the discharge coefficient and swirl ratio. These data are confronted with experimental values, aiming to validate the applied methodology. The importance of the choice of the turbulent model becomes evident in the simulation of internal combustion engines, as well as the interpolation functions used. In the second part the results refer to a transient analysis, considering the valves and piston movement, at 1500 rpm. It is observed the great demand on time step required is observed for the real transient, demonstrating that, for this speed, the smallest time step used, 0.05º (5.5555E-6 s), is still insufficient for some moments of the cycle. Also regarding the time step, it is possible to notice a greater influence in the discharge coefficient than in the swirl ratio. The strong dependence of the turbulence model on the results is once again confirmed, as expected, since the hypotheses about the physics of the phenomenon are different in each model. The results, regarding the wall treatment, presented no significant differences, when applied together with the SST k-ω turbulence model.