Esta tese apresenta um modelo matemático para simulação numérica do escoamento com turbilhonamento (swirl) em tochas de plasma térmico de arco não transferido que operam em corrente contínua, assim como os resultados obtidos com as simulações para estudo de casos. O modelo magneto-hidrodinâmico (modelo MHD) bidimensional permitiu simular a interação entre o escoamento e o arco elétrico usando uma configuração axissimétrica, que abrange as seguintes regiões: entrada do gás; interior da tocha; jato de plasma livre no ambiente. O modelo foi implementado num código numérico baseado no Método dos Volumes Finitos para a solução numérica das equações governantes. Para os estudos foram simulados casos com diferentes condições operacionais (vazão de gás; intensidade de corrente elétrica; gases plasmogênicos: ar e argônio; intensidade de turbilhonamento). A fim de verificar a qualidade do modelo, alguns resultados foram comparados com a literatura e apresentaram boa concordância: a maior diferença obtida entre valores de temperatura experimentais e valores calculados foi -10%, e a média das diferenças obtidas nas comparações foi de aproximadamente ±3,2%. Os perfis de temperatura e de velocidade obtidos para a região do arco e para o jato de plasma resultante permitiram o estudo do comportamento do escoamento na tocha de plasma em diferentes condições. Conclui-se que o modelo desenvolvido é apto à realização de investigações numéricas do escoamento em tochas de plasma e dos efeitos do turbilhonamento na interação arco/escoamento. / This thesis presents a mathematical model for numerical simulation of swirling flow in DC non-transferred arc thermal plasma torches, as well as the results obtained from simulations to case studies. The two-dimensional magnetohydrodynamic model (MHD model) allowed simulate the interaction between the flow and the electric arc using an axisymmetric configuration, covering the following areas: gas inlet; inside the torch; free jet of plasma in the environment. The model was implemented in a computer code based on the Finite Volume Method (FVM) to enable the numerical solution of the governing equations. For the study, cases were simulated with different operating conditions (gas flow rate; electric current intensity; plasmogenic gases: air and argon; swirl intensity). In order to verify the quality of the model, some results were compared with the literature and showed good agreement: the biggest difference between experimental temperature values and calculated values was 10%, and the average of the differences obtained in the comparisons was approximately ±3.2%. The resulting profiles of temperature and velocity obtained for the region of the arc and the plasma jet allowed the study of the flow behavior in the plasma torch in different conditions. It is concluded that the model developed is able to carry out numerical investigations of the flow in plasma torches and the effects of swirl in the interaction arc/flow.

A esfera é utilizada como domínio computacional na modelagem de diversos fenômenos físicos, como em previsão numérica do tempo. Sua discretização pode ser feita de diversas formas, sendo comum o uso de malha regulares em latitude/longitude. Recentemente, também para melhor uso de computação paralela, há uma tendência ao uso de malhas mais isotrópicas, dentre as quais a icosaédrica. Apesar de já existirem modelos atmosféricos que usam malhas icosaédricas, não há consenso sobre as metodologias mais adequadas a esse tipo de malha. Nos propusemos, portanto, a estudar em detalhe diversos fatores envolvidos no desenvolvimento de modelos atmosféricos globais usando malhas geodésicas icosaédricas. A discretização usual por volumes finitos para divergente de um campo vetorial utiliza como base o Teorema da Divergência e a regra do ponto médio nas arestas das células computacionais. A distribuição do erro obtida com esse método apresenta uma forte relação com características geométricas da malha. Definimos o conceito geométrico de alinhamento de células computacionais e desenvolvemos uma teoria que serve de base para explicar interferências de malha na discretização usual do divergente. Destacamos os impactos de certas relações de alinhamento das células na ordem da discretização do método. A teoria desenvolvida se aplica a qualquer malha geodésica e também pode ser usada para os operadores rotacional e laplaciano. Investigamos diversos métodos de interpolação na esfera adequados a malhas icosaédricas, e abordamos o problema de interpolação e reconstrução vetorial na esfera em malhas deslocadas. Usamos métodos alternativos de reconstrução vetorial aos usados na literatura, em particular, desenvolvemos um método híbrido de baixo custo e boa precisão. Por fim, utilizamos as técnicas de discretização, interpolação e reconstrução vetorial analisadas em um método semi-lagrangiano para o transporte na esfera em malhas geodésicas icosaédricas. Realizamos experimentos computacionais de transporte, incluindo testes de deformações na distribuição do campo transportado, que mostraram a adequação da metodologia para uso em modelos atmosféricos globais. A plataforma computacional desenvolvida nesta tese, incluindo geração de malhas, interpolações, reconstruções vetoriais e um modelo de transporte, fornece uma base para o futuro desenvolvimento de um modelo atmosférico global em malhas icosaédricas. / Spherical domains are used to model many physical phenomena, as, for instance, global numerical weather prediction. The sphere can be discretized in several ways, as for example a regular latitude-longitude grid. Recently, also motivated by a better use of parallel computers, more isotropic grids have been adopted in atmospheric global circulation models. Among those, the icosahedral grids are promising. Which kind of discretization methods and interpolation schemes are the best to use on those grids are still a research subject. Discretization of the sphere may be done in many ways and, recently, to make better use of computational resources, researchers are adopting more isotropic grids, such as the icosahedral one. In this thesis, we investigate in detail the numerical methodology to be used in atmospheric models on icosahedral grids. The usual finite volume method of discretization of the divergence of a vector field is based on the divergence theorem and makes use of the midpoint rule for integration on the edges of computational cells. The error distribution obtained with this method usually presents a strong correlation with some characteristics of the icosahedral grid. We introduced the concept of cell alignment and developed a theory which explains the grid imprinting patterns observed with the usual divergence discretization. We show how grid alignment impacts in the order of the divergence discretization. The theory developed applies to any geodesic grid and can also be used for other operators such as curl and Laplacian. Several interpolation schemes suitable for icosahedral grids were analysed, including the vector interpolation and reconstruction problem on staggered grids. We considered alternative vector reconstruction methods, in particular, we developed a hybrid low cost and good precision method. Finally, employing the discretizations and interpolations previously analysed, we developed a semi-Lagrangian transport method for geodesic icosahedral grids. Several tests were carried out, including deformational test cases, which demonstrated that the methodology is suitable to use in global atmospheric models. The computational platform developed in this thesis, including mesh generation, interpolation, vector reconstruction and the transport model, provides a basis for future development of global atmospheric models on icosahedral grids.

Submitted by Renata Lopes (renatasil82@gmail.com) on 2017-03-02T18:09:12Z No. of bitstreams: 1 deniseschimitzdecarvalhotristao.pdf: 734334 bytes, checksum: 9fda9bda660d5bfec3204e328fe66d1c (MD5) / Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2017-03-06T19:58:58Z (GMT) No. of bitstreams: 1 deniseschimitzdecarvalhotristao.pdf: 734334 bytes, checksum: 9fda9bda660d5bfec3204e328fe66d1c (MD5) / Made available in DSpace on 2017-03-06T19:58:58Z (GMT). No. of bitstreams: 1 deniseschimitzdecarvalhotristao.pdf: 734334 bytes, checksum: 9fda9bda660d5bfec3204e328fe66d1c (MD5) Previous issue date: 2013-08-30 / O desenvolvimento de modelos matemáticos e métodos computacionais para a simulação de escoamentos em meios porosos é de grande interesse, devido à sua aplicação em diversas áreas da engenharia e ciências aplicadas. Em geral, na simulação numérica de um modelo de escoamento em meios porosos, são adotadas estratégias de desacoplamento dos sistemas de equações diferenciais parciais que o compõem. Este estudo recai sobre esquemas numéricos para leis de conservação hiperbólicas, cuja aproximação é não-trivial. Os esquemas de volumes finitos de alta resolução baseados no algoritmo REA (Reconstruct, Evolve, Average) têm sido empregados com considerável sucesso para a aproximação de leis de conservação. Recentemente, esquemas centrais de alta ordem, baseados nos métodos de Lax-Friedrichs e de Rusanov (Local Lax-Friedrichs) têm sido apresentados de forma a reduzir a excessiva difusão numérica característica destes esquemas de primeira ordem. Nesta dissertação apresentamos o estudo e a aplicação de esquemas de volumes finitos centrais de alta ordem para equações hiperbólicas que aparecem na modelagem de escoamentos em meios porosos. / The development of mathematical models and computational methods for the simulation of flow in porous media has a great interest because of its applications in engineering and other sciences. In general, in order to solve numerically the flow model in porous media the system of partial differential equations are decoupled. This study focus on the numerical schemes for the hyperbolic conservation laws, which solution is non-trivial. The finite volume schemes based on high order algorithm REA (Reconstruct, Evolve, Average) have been used with considerable success for the numerical solution of the conservation laws. Recently, high-order central schemes, based on the methods of Lax-Friedrichs and Rusanov (Local Lax-Friedrichs) have been presented, they reduce the excessive numerical diffusion presented in the first order schemes. In this dissertation we present the study and application of the high-order finite volume central schemes for hyperbolic equations as appear in the porous media flow modeling.

Orientador: Sônia Maria Gomes, Jorge Stolfi / Tese (doutorado) - Universidade Estadual de Campinas, Instituto de Matemática, Estatística e Computação Científica / Made available in DSpace on 2018-08-19T12:39:30Z (GMT). No. of bitstreams: 1 Castro_DouglasAzevedo_D.pdf: 8872633 bytes, checksum: a17b2761789c6a831631ac143fdf5ca7 (MD5) Previous issue date: 2011 / Resumo: O objetivo desta tese é desenvolver algoritmos baseados em malhas e bases funcionais inovadoras usando técnicas de multiescala para aproximação de funções e resolução de problemas de equações diferenciais. Para certas classes de problemas, é possível incrementar a eficiência dos algoritmos de multiescala usando bases adaptativas, associadas a malhas construídas de forma a se ajustarem com o fenômeno a ser modelado. Nesta abordagem, em cada nível da hierarquia, os detalhes entre a aproximação desse nível e a aproximação definida no próximo nível menos refinado pode ser usada como indicador de regiões que necessitam de mais ou menos refinamento. Desta forma, em regiões onde a solução é suave, basta utilizar os elementos dos níveis menos refinados da hierarquia, enquanto que o maior refinamento é feito apenas onde a solução tiver variações bruscas. Consideramos dois tipos de formulações para representações multiescala, dependendo das bases adotadas: splines diádicos e wavelets. A primeira abordagem considera espaços aproximantes por funções splines sobre uma hierarquia de malhas cuja resolução depende do nível. A outra abordagem considera ferramentas da analise wavelet para representações em multirresolução de médias celulares. O enfoque está no desenvolvimento de algoritmos baseados em dados amostrais d-dimensionais em malhas diádicas que são armazenados em uma estrutura de árvore binária. A adaptatividade ocorre quando o refinamento é interrompido em algumas regiões do domínio, onde os detalhes entre dois níveis consecutivos são suficientemente pequenos. Um importante aspecto deste tipo de representação é que a mesma estrutura de dados é usada em qualquer dimensão, além de facilitar o acesso aos dados nela armezenados. Utilizamos as técnicas desenvolvidas na construção de um método adaptativo de volumes finitos em malhas diádicas para a solução de problemas diferenciais. Analisamos o desempenho do método adaptativo em termos da compressão de memória e tempo de CPU em comparação com os resultados do esquema de referência em malha uniforme no nível mais refinado. Neste sentido, comprovamos a eficiência do método adaptativo, que foi avaliada levando-se em consideração os efeitos da escolha de diferentes tipos de fluxo numérico e dos parâmetros de truncamento / Abstract: The goal of this thesis is to develop algorithms based on innovative meshes and functional bases using multiscale techniques for function approximation and solution of differential equation problems. For certain classes of problems, one can increase the efficiency of multiscale algorithms using hierarchical adaptive bases, associated to meshes whose resolution varies according to the local features of the phenomenon to be modeled. In this approach, at each level of the hierarchy the details-differences between the approximation for that level and that of the next coarser level-can be used as indicators of regions that need more or less refinement. In this way, in regions where the solution is smooth, it suffices to use elements of the less refined levels of the hierarchy, while the maximum refinement is used only where the solution has sharp variations. We consider two classes of formulations for multiscale representations, depending on the bases used: dyadic splines and wavelets. The first approach uses approximation spaces consisting of spline functions defined over a mesh hierarchy whose resolution depends on the level. The other approach uses tools from wavelet analysis for multiresolu-tion representations of cell averages. The focus is on the development of algorithms based on sampled d-dimensional data on dyadic meshes which are stored in a binary tree structures. The adaptivity happens when the refinement is interrupted in certain regions of the domain, where the details between two consecutive levels are sufficiently small. This representation greatly simplifies the access to the data and it can be used in any dimension. We use these techniques to build an adaptive finite volume method on dyadic grids for the solution of differential problems. We analyze the performance of the method in terms of memory compression and CPU time, comparing it with the reference scheme (which uses a uniform mesh at the maximum refinement level). In these tests, we confirmed the efficiency of the adaptive method for various numeric flow formulas and various choices of the thresholding parameters / Doutorado / Matematica Aplicada / Doutor em Matemática Aplicada

Este trabalho tem como principal objetivo o desenvolvimento de um método numérico para simulação das grandes escalas de escoamentos turbulentos tridimensionais utilizando uma modelagem de turbulência baseada em filtragem temporal (denominada TLES - Temporal Large Eddy Simulation). O método desenvolvido combina discretizações temporais com ordem de mínima precisão 2 (Adams-Bashforth, QUICK, Runge-Kutta), um método de projeção de ordem 2, com discretizações espaciais também de ordem 2 obtidas pelo método de volumes finitos. Esta metodologia foi empregada na simulação de problemas teste turbulentos como o canal e a cavidade impulsionada, sendo este último resultado simulado pela primeira vez com modelagem TLES. Os resultados mostram uma excelente concordância quando comparado com resultados de simulações diretas (DNS) e dados experimentais, superando resultados clássicos obtidos com formulação LES com filtragem espacial. / The main objective of this work is to develop a numerical method for large eddy simulation of tridimensional turbulent flows using a model based on temporal filtering (TLES - Temporal Large Eddy Simulation). The developed method combines at least 2nd order temporal discretizations (Adams-Bashforth, QUICK, Runge-Kutta), a 2nd order projection method, and 2nd order spatial discretizations obtained by the finite volume method. This methodology was employed to the simulation of turbulent benchmark problems such as channel and lid-driven cavity flows. The latter is simulated for the first time using a TLES turbulence modelling. Results show excellent agreement when compared to Direct Numerical Simulations (DNS) and experimental data, with better results than classical results produced by standard LES formulation with spatial filtering.

Submitted by Johnny Rodrigues (johnnyrodrigues@ufcg.edu.br) on 2018-03-23T01:44:32Z No. of bitstreams: 1 GUSTAVO GOMES SAMPAIO CURSINO - TESE PPGEQ 2016.pdf: 5939681 bytes, checksum: 518baf4150ea2fb7085706252276b9fa (MD5) / Made available in DSpace on 2018-03-23T01:44:32Z (GMT). No. of bitstreams: 1 GUSTAVO GOMES SAMPAIO CURSINO - TESE PPGEQ 2016.pdf: 5939681 bytes, checksum: 518baf4150ea2fb7085706252276b9fa (MD5) Previous issue date: 2016-04-18 / Este trabalho teve o propósito de determinar o comportamento dos gases na seção de radiação de um combustor de ar que pertence a uma planta industrial. O corpo metálico do equipamento rompeu em seu primeiro ano de operação, devido a um problema conceitual em sua geometria. A fluidodinâmica computacional (CFD), por meio do método dos volumes finitos, foi utilizada para desenvolver um modelo tridimensional que pudesse reproduzir o perfil de temperatura e o comportamento do fluxo do ar de combustão no equipamento. Na simulação, através do uso do software ANSYS CFX, foram utilizados: (i) o modelo de turbulência Reynolds Stress Model (RSM); (ii) as malhas hexaédrica, tetraédrica e prismática; (iii) o modelo de radiação P-1; e (iv) o modelo de combustão Eddy Dissipation Concept (EDC). Como resultado, foram apresentadas quatro possíveis mudanças na geometria do combustor de ar que, caso adotadas, eliminariam os riscos de novas falhas e garantiriam a continuidade operacional da unidade de processo. / This paper has the objective to describe the behavior of the flow and temperature of the flue gas in the radiation section of the vessel used to preheat air in a combustor. The equipment failed in its first operational year, due to a conceptual problem in its geometry. The CFD code based on finite volume method was applied to simulate the physical model of combustor using the ANSYS CFX software, reproducing the main features of the preheater. The simulation had considered: (i) Reynolds Stress Model (RSM) as turbulence model, (ii) The meshes applied were the hexahedral, tetrahedral and prismatic, (iii) P-1 was used as the radiation model and (iv) Eddy Dissipation Concept (EDC) as combustion model. Through the simulation was possible to propose four different kind of combustor geometry modification, that the application of anyone of them would eliminate the risk of new failures, ensuring the unit production availability.

The advance of the materials science has motivated the advent of composite materials with different characteristics that assure high performance thermo-mechanical, such as the advanced fiber reinforced composites and those that present a gradual variation in its microstructure. Nowadays, many computational models and analytical methods are being employed for evaluation of the behavior of such materials. An alternative technique, applied to the steady-state thermo-mechanical analysis, which considers the coupling between microstructure and macrostructure behaviors, is that originally denominated of Higher-Order Theory. In this work is used the same theoretical base of the Higher-Order Theory, with a simplification in the discretization and assembly of the system of equations. In this way, this theory presents some similarities in relation to the finite-volume technique used in fluid dynamics problems, reason for which is enough reasonable to adopt the denomination finitevolume theory for this method. Besides, as a contribution of this study, it is presented a parametric formulation that allows a larger flexibility in the mesh generation and a reduction of the problem in relation to the number of variables, particularly appropriate for analysis of structures with curved contour. The formulation was also extended for the accomplishment of transient thermo-mechanical analysis. In the present study, a three-dimensional formulation of the method is also used for the determination of the effective properties of fiber reinforced composites and particulate materials, where comparisons were accomplished with micromechanics simplified models and with those based on the mean field theory (selfconsistent, Mori-Tanaka and differential scheme). Besides, there is a series of numerical applications in bi-dimensional thermo-elastic and elastic problems, where are accomplished verifications using analytical solutions and comparisons with the finite element method. / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / O avanço da ciência de materiais tem proporcionado o advento de materiais compósitos com características peculiares que asseguram elevado desempenho termo-mecânico, tais como os compósitos avançados reforçados por fibras e aqueles dotados de microestrutura com variação gradual. Atualmente, muitos modelos computacionais, assim como métodos analíticos, vêm sendo empregados para avaliação do comportamento de tais materiais. Uma técnica alternativa, voltada à análise termo-mecânica em regime estacionário, na qual o comportamento do material é analisado considerando-se o acoplamento entre microestrutura e macroestrutura, é aquela originalmente denominada de Higher-Order Theory. Neste trabalho, utiliza-se a mesma base teórica da Higher-Order Theory, com uma simplificação em termos de discretização e montagem do sistema de equações. Neste sentido, esta teoria apresenta algumas semelhanças em relação à técnica de volumes finitos usada em problemas de dinâmica dos fluidos, razão pela qual é bastante razoável adotar a denominação teoria de volumes finitos para o método. Além disso, como uma contribuição deste estudo, apresenta-se uma formulação paramétrica que permite uma maior flexibilidade na geração da malha e uma diminuição do problema em relação ao número de incógnitas, particularmente apropriada para análise de estruturas com contorno curvo. A formulação também foi ampliada para possibilitar a execução de análises termo-mecânicas transientes. No presente estudo, também é utilizada uma formulação tridimensional do método para a determinação das propriedades efetivas de materiais compósitos reforçados por fibras e particulados, onde foram realizadas comparações com modelos simplificados da micromecânica e com aqueles baseados na teoria de campos médios (Auto-consistente, Mori-Tanaka e Esquema Diferencial). Além disso, há uma série de aplicações numéricas em problemas termo-elásticos e elásticos bidimensionais, onde são realizadas verificações a partir de soluções analíticas e comparações com o método dos elementos finitos.

Neste trabalho apresentamos um esquema numérico para a simulação computacional de escoamentos bifásicos, água-óleo, em reservatórios de petróleo. O modelo matemático consiste em um sistema de equações diferenciais parciais não-linear nas incógnitas velocidade, pressão e saturação. Uma quebra de operadores a dois níveis possibilita uma maior eficiência ao método permitindo que a velocidade, fornecida pelo problema de velocidade-pressão, seja atualizada somente para determinados intervalos de tempo associados ao problema de transporte advectivo-difusivo em termos da saturação. O método dos volumes finitos é empregado na resolução numérica do problema de velocidade-pressão e do transporte de massa por advecção e difusão. Na solução do problema de transporte de massa utilizamos limitadores de fluxo na aproximação dos termos advectivos e diferenças centradas para os termos difusivos. O nosso simulador foi validado a partir de confrontações dos seus resultados com as soluções teóricas conhecidas para os problemas unidimensionais, equações de Burgers e de Buckley-Leverett, e com outros resultados numéricos em se tratando do escoamento bifásico água-óleo bidimensional em meios porosos heterogêneos. / A new numerical method is proposed for the solution of two-phase flow problem in petroleum reservoirs. The two-phase (water and oil) flow problem is governed by a pressure-velocity equation coupled to a saturation equation. For computational eficiency an operator spliting technique is used; distinct time steps can be used for the computation of transport and pressure-velocity problems. The finite volume method is used in the numerical solution of the velocity-pressure and mass transport problems. A flux limiter is used for the numerical discretization of the advective terms while centered schemes are employed for the diffusion terms in the mass transport problem. In the validation of our numerical method we compared numerical and theoretical solutions for one dimensional problems, Burgers and Buckley-Leverett equations, and compared our numerical results to others, in the case of oil-water flows in two dimensions for an heterogeneous porous media.

Neste trabalho apresentamos um esquema numérico para a simulação computacional de escoamentos bifásicos, água-óleo, em reservatórios de petróleo. O modelo matemático consiste em um sistema de equações diferenciais parciais não-linear nas incógnitas velocidade, pressão e saturação. Uma quebra de operadores a dois níveis possibilita uma maior eficiência ao método permitindo que a velocidade, fornecida pelo problema de velocidade-pressão, seja atualizada somente para determinados intervalos de tempo associados ao problema de transporte advectivo-difusivo em termos da saturação. O método dos volumes finitos é empregado na resolução numérica do problema de velocidade-pressão e do transporte de massa por advecção e difusão. Na solução do problema de transporte de massa utilizamos limitadores de fluxo na aproximação dos termos advectivos e diferenças centradas para os termos difusivos. O nosso simulador foi validado a partir de confrontações dos seus resultados com as soluções teóricas conhecidas para os problemas unidimensionais, equações de Burgers e de Buckley-Leverett, e com outros resultados numéricos em se tratando do escoamento bifásico água-óleo bidimensional em meios porosos heterogêneos. / A new numerical method is proposed for the solution of two-phase flow problem in petroleum reservoirs. The two-phase (water and oil) flow problem is governed by a pressure-velocity equation coupled to a saturation equation. For computational eficiency an operator spliting technique is used; distinct time steps can be used for the computation of transport and pressure-velocity problems. The finite volume method is used in the numerical solution of the velocity-pressure and mass transport problems. A flux limiter is used for the numerical discretization of the advective terms while centered schemes are employed for the diffusion terms in the mass transport problem. In the validation of our numerical method we compared numerical and theoretical solutions for one dimensional problems, Burgers and Buckley-Leverett equations, and compared our numerical results to others, in the case of oil-water flows in two dimensions for an heterogeneous porous media.

Orientadores: Rubens Maciel Filho, Liliane Maria Ferrareso Lona / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-04T02:31:19Z (GMT). No. of bitstreams: 1 MendozaMarin_FlorentinoLazaro_D.pdf: 8929110 bytes, checksum: 9041b9e7f02f9a324fa10cde9a962f54 (MD5) Previous issue date: 2004 / Resumo: O modelo determinístico e processo homopolimerização na emulsão do estireno são aplicados em reator tubular contínuo sem e com deflectores angulares internos sob condição isotérmica e não isotérmica. Os resultados de modelagem e simulação foram realizados a estado estacionário, modelo unidimensional, coordenada cilíndrica, fluxo pistão laminar completamente desenvolvido, modelo Smith-Ewart para estimar a conversão do monômero, cinética química de Arrhenius corno modelo de velocidade finita laminar para computar a geração química. O objetivo é modelar, simular e analisar o comportamento do reator de homopolimerização na emulsão do estireno com deflectores angulares inclinados internos, e comparar com reator tubular. Os métodos experimental e matemático-dedutivo foram aplicados para obter resultados, por meio de programação computacional, usando Dinâmica de Fluido Computacional através do método de volumes finitos. As seguintes variáveis como temperatura de reação constante e variável, reator tubular sem e com deflectores, temperatura de alimentação, diâmetro de reator, processo adiabático e exotérmico, calor de reação constante e velocidade axial completamente desenvolvida foram investigados. Os efeitos de conversão de monômero, área transversal interna, temperatura axial, concentração do polímero, radicais e iniciador, outros corno densidade de polímero e monômero, perda de carga e queda de pressão foram determinados e simulados. Os produtos foram caracterizados com Número de Partículas (nucleação homogênea e heterogênea), distribuição de peso molecular, tamanho de partículas de polímero e distribuição de viscosidade. Estes resultados foram validados com resultados da literatura sob condição igualou aproximada. Os resultados sob condições não isotérmicas foram melhores que os resultados isotérmicos em termos de caracterização do polímero. Isso mostra que o desenho alternativo proposto (com deflectores) permite obter o polímero com propriedades melhores em termos de número de partículas, distribuição de peso molecular, distribuição do tamanho de partículas e viscosidade / Abstract: Deterministic model and emulsion homopolymerization process of styrene are applied in continuous tubular reactor without and with internal angular baffles under isothermic and no isothermic conditions. The modeling and simulation results were approximate to steady state, one-dimensional model, cylindrical coordinate, fully developed laminar plug flow, Smith-Ewart model to estimate the monomer conversion, Arrhenius chemical kinetics as laminar finite-rate model to compute chemical source. The objective is to model, simulate and to analyze the emulsion homopolymerization reactor performance of styrene with internal-inc1ined angular baffles, and to compare with continuous tubular reactor. The experimental and mathematical-deductive methods were applied to obtain results, by means of computational programming, using Computational Fluid Dynamics (program code), finite volume method. The following variables such as constant and variable reaction temperature, tubular reactor without and with baffles, feed temperature, reactor diameter, adiabatic and exothermic process, constant reaction heat and fully developed axial velocity were investigated. The monomer conversion, internal transversal are a, axial temperature, concentration of polymer, radicals and initiator, others as density of polymer and monomer, head loss and pressure drop effects were determined and simulated. The products were characterized by partic1es number (homogeneous and heterogeneous nuc1eation), molecular weight distribution, polymer partic1es size and polymer viscosity distribution. These results were validated with literature results under same or approximate condition. The results under no isothermic conditions were better than isothermic results in terms of polymer characterization. It is shown that the proposed alternative design (with baffles) allow to obtain the polymer with better properties in terms of number of partic1es, molecular weight distribution, particle size distribution and viscosity / Doutorado / Desenvolvimento de Processos Químicos / Mestre em Engenharia Química