Visando compreender os fluxos na fase líquida do crescimento de silício pelo método Czochralski, é feita a Simulação Numérica do silício fundido, resolvendo-se as equações que governam o fenômeno da convecção forçada no fluido: Balanço de Quantidade de Movimento e Balanço de Massa. A técnica numérica escolhida é a de Elementos Finitos, onde é utilizada a formulação de Galerkin, com aproximações quadráticas nas componentes da velocidade e linear na pressão. A partir de várias combinações de rotações cadinho-cristal, os perfis de velocidade obtidos são analisados com relação aos efeitos de incorporação de impurezas e/ou dopantes no cristal em crescimento. / In order to visualise the flow conditions during crystal growth of Silicon by the Czochralski technique, a numerical simulation is done. It is used the Finite Element Method with the Galerkin Formulation , and with quadratic approximations on the components of the velocity and linear approximations on the pressure. Many combinations of crystal and crucible rotations are analised and discussed considering optimal growth conditions.

O presente trabalho é focado no estudo dos violentos fenômenos de impacto hidrodinâmico que podem prejudicar a operação de plataformas offshore. São três os fenômenos abordados neste trabalho: o green water, o wave runup e o sloshing. O fenômeno de wave runup consiste na projeção vertical de uma coluna de água devido ao impacto de ondas em estruturas transversais. O fenômeno de green water consiste no embarque de água no convés, podendo danificar os equipamentos da planta de produção. Por fim, o sloshing consiste no movimento violento de fluído em tanques parcialmente preenchidos, resultando em perigosos carregamentos em suas paredes. Tais fenômenos possuem natureza altamente não linear e sua análise, considerando-se toda a sua complexidade, ainda constitui um desafio para a engenharia naval e oceânica. Os métodos de partículas têm se destacado no tratamento de tais fenômenos envolvendo interação fluído-estruturas, grandes deformações e fragmentação de superfície livre. Desta maneira, optou-se pelo emprego do método de partículas Moving Particles Semi-Implicit (MPS) neste trabalho para o estudo dos fenômenos de impacto hidrodinâmico. O MPS é um método totalmente lagrangeano para escoamentos incompressíveis. Para os três fenômenos abordados neste trabalho há uma primeira etapa de validação, na qual os resultados numéricos são comparados a resultados experimentais da literatura. Uma segunda etapa é baseada na aplicação do método numérico na análise de ferramentas para a mitigação dos esforços resultantes do impacto hidrodinâmico. Nesta etapa é investigada a influência do formato da proa no fenômeno de green water e a utilização de anteparas fixas e flutuantes para a mitigação de sloshing em tanques. / The present work is focused in the study of the violent hydrodynamic impact phenomenon which could jeopardize the offshore platforms operation. In this work three different phenomena involving hydrodynamic impact are studied: green water, wave runup and sloshing. The wave runup consists in the vertical projection of a water column due to wave impact on a transversal structure, such as submersible columns. The green water consists in the water boarding on the deck which could damage the equipment over the oil platform deck. Finally, the sloshing phenomenon is the violent movement of fluid in partially filled tanks, resulting in dangerous impact loads at its walls. The hydrodynamic impact phenomenon has strongly non linear nature and is still a challenge for the naval and offshore engineering its analysis considering all its complexity. The particle methods present advantages in the analysis of phenomena involving fluid structure interaction, large free surface deformation, fragmentation and merging. Therefore, in the present study the Moving Particles Semi-Implicit (MPS) method is used. The MPS is a fully lagrangian method for the simulation of incompressible flows. For the three phenomena studied in the present work a first step of validation is performed. In the validation step the numerical results obtained by the particle method are compared to experimental data presented in the literature. The second step consists in the application of the numerical method to investigate simple mechanisms to mitigate the hydrodynamic impact loads. For example, the effect of the bow shape in the green water phenomenon is studied. Also in this step the use of fixed and floating baffles in order to suppress the sloshing phenomenon are investigated.

O fenômeno da transferência de massa entre dois meios é altamente dependente da concentração do componente a ser transferido e da velocidade na camada limite de troca. A medida de velocidade na superfície da água em uma interface ar/água ainda não é realidade devido a falta de instrumental adequado. Procurando sanar este obstáculo, esta pesquisa propôs o desenvolvimento de uma sonda ótica para medida de turbulência na superfície da água em um tanque agitado por grade oscilante. O funcionamento desta sonda foi investigado e suas medidas comparadas a medidas realizadas pelo instrumento microADV, de forma a verificar se o parâmetro medido pelo novo instrumento era relacionado à componente vertical da velocidade próximo à superfície. O comportamento do escoamento no tanque foi modelado com o uso do software Open-FOAM, visando uma maior compreensão do fenômeno. A pesquisa resultou na modelação do escoamento e no desenvolvimento de um dispositivo auxiliar que possibilita o uso do microADV em medidas próximas à superfície livre. A comparação das medidas adquiridas pela sonda ótica e microADV mostram que, apesar de as medidas da sonda terem relação com a turbulência no tanque, este parâmetro turbulento apresenta fraca correlação com a medida da componente vertical da velocidade. / The mass transfer phenomena is highly dependent on the concentration of the component to be transferred and the velocity on the exchange layer. The velocity measurement at the surface of the water in an air / water interface is not yet a reality due to lack of appropriate instruments. In order to overcome this obstacle, this research proposed the development of an optical probe to measure turbulence on the surface of water in an oscillating grid tank. We investigated the operation of the proposed probe comparing its measurements to measurements made by microADV to check the relation between the parameter measured by the new instrument and the vertical velocity near to the surface. The flow behavior in the tank was MODELLED using the software OpenFOAM seeking a greater understanding of the phenomenon. The research resulted in the modeling of the flow in the tank and the development of an auxiliary device that allows the use of microADV in near-surface measurements. Although the probe measurement have relation with the turbulence in the tank, the comparison of measurements acquired by the optical probe and microADV shows the turbulent parameter presenting weak correlation to the measurement of the vertical velocity.

O incêndio florestal é uma complexa combinação da energia liberada na forma de calor devido à combustão dos produtos oriundos da pirólise da vegetação e o transporte dessa energia para o ar e para a vegetação à sua volta. O primeiro é o domínio da química e ocorre na escala de moléculas e o segundo é o domínio da física e ocorre em escalas de até quilômetros. É a interação desses processos sobre uma ampla gama de escalas temporais e espaciais envolvidas no incêndio florestal que faz a modelagem do seu comportamento uma tarefa tão difícil. A propagação do incêndio através de vegetação rasteira e folhas mortas foi simulada numericamente usando a formulação física do WFDS. A abordagem utilizada foi tridimensional e transiente, e baseada em uma descrição dos fenômenos físicos que contribuem para a propagação de um incêndio de superfície através de uma camada de combustível. Neste cenário de incêndio, existem duas regiões: vegetação e ar, cada uma com suas propriedades físicas e químicas e, embora elas precisem ser integradas no mecanismo de solução, há diferentes fenômenos que ocorrem em cada uma. Na região de vegetação, a abordagem é representá-la como partículas submalha cercadas de ar. O caráter heterogêneo da vegetação, como sua natureza, folhagens, pequenos galhos, etc. foi levado em conta usando propriedades físicas médias características da floresta amazônica. Os fenômenos na região de vegetação são a evaporação da sua umidade, a pirólise e a transferência de calor por radiação e por convecção. Na região do ar, a combustão com chama ocorre em um ambiente turbulento, onde as transferências de calor por radiação e por convecção desempenham um papel significativo. Para incorporar a radiação dos gases de combustão, o modelo físico emprega o método de volumes finitos, que resolve a equação de transferência de calor por radiação como uma equação de transporte para um número finito de discretos ângulos sólidos, e que pode ser usado em uma ampla faixa de espessuras óticas e meios participantes. A combustão turbulenta para a fase gasosa é modelada com base no modelo Eddy Dissipation Concept (EDC). O modelo de combustão turbulenta adota a hipótese de reação química infinitamente rápida entre o combustível e o ar e é controlado apenas pela velocidade de mistura desses reagentes. Esse modelo representa bem a física de incêndios em ambientes ventilados, como é o caso dos incêndios florestais. Para incluir os efeitos do transporte turbulento é utilizado o método Large Eddy Simulation (LES), que calcula explicitamente as grandes estruturas turbulentas, mas trata a dissipação e a cascata inercial em escalas menores usando aproximações na escala submalha. As regiões de vegetação e ar trocam massa e energia. O comportamento da mistura gasosa resultante da degradação térmica da vegetação e das reações de combustão é regido pelas equações de Navier-Stokes. As equações que regem os modelos físicos são formuladas como equações diferenciais parciais que são resolvidas por métodos numéricos. O método utilizado para discretização das equações é o método de diferenças finitas em malha deslocada. O modelo numérico utilizado resolve as equações de Navier-Stokes para fluidos compressíveis usando o filtro de Favre. A dissipação de energia cinética é obtida através de um fechamento simples para a tensão turbulenta: o modelo de coeficiente constante de Deardorff. O transporte turbulento de energia e massa é contabilizado pelo uso, respectivamente, de números de Prandtl e de Schmidt turbulentos constantes. Os resultados das simulações do modelo físico descrito foram comparados aos dados experimentais obtidos em campo para a propagação do incêndio na floresta amazônica. Apesar da idealização das condições de combustível, vento e as incertezas dos dados experimentais, as previsões do modelo estão na mesma ordem de grandeza dos experimentos. As taxas de propagação do incêndio experimentais variam de 0,12 +/-0,06 a 0,35+/-0,07 m/min. Mesmo considerando-se o desvio padrão da taxa de propagação do incêndio experimental, os valores das taxas simuladas ficaram dentro do erro experimental somente em dois de sete casos. As simulações mostraram que os parâmetros importantes para o modelo são a área superficial por volume da vegetação, sua massa específica aparente e sua umidade. Como o coeficiente de absorção por radiação é função direta da massa específica aparente e da área superficial por volume da vegetação, esses parâmetros afetam o comportamento numérico do incêndio de superfície. De acordo com os resultados das simulações numéricas, a umidade da vegetação também tem importância no incêndio de superfície. A temperatura inicial da vegetação e a umidade do ar na faixa de variação analisada não influenciam a taxa de propagação do incêndio. As simulações também mostraram que o processo de radiação é muito importante, e afeta diretamente todos os demais processos e a taxa de propagação do incêndio. A convecção tem importância muito menor que a radiação na condição de ausência de vento externo. A coerência das taxas de propagação do incêndio experimental e numérica em função da massa específica aparente de material combustível e da umidade da vegetação foi investigada. O modelo numérico é coerente em todas as nove combinações de casos. Já o experimento é coerente em quatro combinações. Com base nas comparações entre cada dois casos experimentais e as respectivas simulações numéricas, nota-se que as taxas de propagação a partir das simulações numéricas foram mais coerentes que as experimentais. / Forest fire is a complex combination of energy released as heat due to the combustion of the products from the vegetation pyrolysis and the transport of this energy to the surrounding air and vegetation. The first is the domain of chemistry and occurs on the molecular scale, and the second is the domain of physics and occurs at scales up to kilometers. It is the interaction of these processes on a wide range of temporal and spatial scales involved in forest fires that makes modeling its behavior such a challenging task. The spread of fire through small plants and dead leaves was simulated numerically using WFDS physical formulation. The approach used was three-dimensional and transient, based on a description of the physical phenomena that contribute to the spread of a surface fire through a layer of fuel. In this fire scenario, there are two regions: vegetation and air, each one with its physical and chemical properties and, although they need to be integrated into the solution mechanism, there are different phenomena that occur in each one. In the vegetation region, the approach is to represent it as subgrid particles surrounded by air. The heterogeneity of the vegetation, such as its nature, leaves, twigs, etc. was taken into account by using average physical properties that are representative of the Amazon forest. The phenomena in the vegetation region are the evaporation of its moisture, pyrolysis, heat transfer by radiation and convection. In the air region, the flaming combustion occurs in a turbulent environment, and heat transfer by radiation and convection play a significant role. To incorporate the radiation from the combustion gases, the physical model employs the finite volumes method, solving the radiation transfer equation as a transport equation for a finite number of discrete solid angles, which can be used in a wide range of optical thicknesses and participating media. Turbulent combustion for the gaseous phase is modeled using the Eddy Dissipation Concept (EDC) model. The mixing controlled turbulent combustion model adopts the assumption of infinitely fast chemical reaction between the fuel and air. This model represents well the fire physics in ventilated areas, as is the case of forest fires. To include the turbulent flow effects, it is used the Large Eddy Simulation (LES) method, which explicitly calculates the large turbulent structures, but models the dissipation and inertial cascade using approximations in the sub-grid scale. The vegetation and air regions exchange mass and energy. The behavior of the gas mixture resulting from the vegetation thermal degradation and combustion reactions is governed by the Navier-Stokes equations. The equations governing the physical model are formulated as partial differential equations, which are solved by numerical methods. The method used for discretization of the equations is the finite difference method on a staggered grid. The numerical model solves the Navier-Stokes equations for compressible fluids using the Favre filter. Dissipation of kinetic energy is achieved through a simple closure for the turbulent stress: the constant coefficient Deardorff model. The turbulent transport of heat and mass is accounted for by use of constant turbulent Prandtl and Schmidt numbers, respectively. The physical model simulation results were compared to experimental data obtained in the field for the spread of fire in the Amazon forest. Despite of the idealized conditions of fuel, wind and the uncertainties of the experimental data, the model predictions and the experiments are in the same order of magnitude. Experimental rate of spread range from 0.12 +/- 0.06 to 0.35 +/- 0.07 m/min. Even considering rate of spread experimental standard deviation, simulated rate values were within experimental error only in two of seven cases. The simulations showed that the important parameters for the model are the vegetation surface area to volume ratio, its bulk density and moisture. As the radiation absorption coefficient is a direct function of vegetation bulk density and surface area to volume ratio, these parameters affect the numeric behavior of the surface fire. According to the numerical simulations results, vegetation moisture is also important in the surface fire scenario. Vegetation initial temperature and air humidity in the range analyzed does not influence the rate of spread. The simulations also showed that the radiation process is very important and directly affects all other processes and rate of spread. Convection heat transfer has much less significance than radiation heat transfer in the absence of external wind. The consistency of the experimental and numerical rate of spread, as a function of combustible material bulk density and vegetation moisture was investigated. The numerical model is consistent in all nine case combinations. The experiment is consistent in four cases. Based on comparisons between each two experiments and their numerical simulations, it is noted that the rate of spread variation from the numerical simulation is more consistent than the experimental one.

Este trabalho tem como propósito, investigar a validade da lei de Hess-Murray, através da experimentação numérica, aplicando o método do Design Construtal associado ao método de otimização da busca exaustiva, no escoamento interno em estruturas em forma de T com seção circular. Variação do número de Reynolds, escoamento de fluidos newtonianos e não newtonianos, estrutura em forma de T com paredes impermeáveis e permeáveis, foram as principais características avaliadas para confrontar a lei de Hess-Murray. O estudo proposto assume escoamento tridimensional, laminar, incompressível, regime permanente e propriedades fluidodinâmicas constantes, sendo o regime de escoamento governado pelo número de Reynolds (Re). O objetivo principal consiste em determinar as configurações ótimas que facilitem o escoamento de fluido ou minimizem as resistências ao escoamento, quando a área ocupada pelos dutos (A) e o volume ocupado pelos dutos (V) são mantidos constantes, variando-se as razões de diâmetros (aD) e comprimentos (aL). As equações de conservação de massa e quantidade de movimento, foram resolvidas através do método de volumes finitos. A geometria foi discretizada através de uma malha tridimensional composta por aproximadamente 1.950.000 elementos. Como resultados, obteve-se as geometrias ótimas que apresentaram resistências ao escoamento até 30 vezes menor do que as demais configurações. Além disso, foi possível verificar que a lei de Hess-Murray nem sempre é válida, visto que o sistema adapta sua geometria ótima para cada condição de escoamento, a fim de proporcionar a melhor arquitetura de escoamento para atender ao objetivo de minimizar as resistências ao escoamento em acordo com a Lei Construtal. Esta dissertação avançou no presente estado da arte, pois desenvolveu um modelo tridimensional sem simplificações, aplicado ao sistema de escoamento de fluidos em estrutura em forma de T utilizando o método do Design Construtal, validando os resultados analíticos apresentados na bibliografia e apresentando novas referências que permitem ampliar a complexidade dos sistemas de escoamento bem como a implementação de métodos de otimização mais avançados. / This work investigates, through the numerical experimentation together with the Construtal Design method, the Hess-Murray Law in the internal flow in T-shaped structures with a circular section for the laminar flow of Newtonian and non-Newtonian fluids with impermeable and permeable walls, determining the optimal configurations that facilitate fluid flow or minimize flow resistance. The geometric global constants, the volume occupied by the ducts (V) and the area occupied by the ducts (A), delimit the space occupied by the T-shaped structure and the degrees of freedom, the ratio between the diameter of the parent duct and daughter (aD) and the ratio between parent duct length and daughter (aL), are the main geometric parameters to be evaluated. The proposed study is assumed three-dimensional, laminar, incompressible, permanent and constant fluidodynamic properties being the flow regime governed by Reynolds number (Re). Construtal Design method, associated with the exhaustive search, was used to determine the global geometric constants, degrees of freedom and objective function in the geometric evaluation of the system. The numerical solution of the mass conservation and momentum equations is solved based on the finite volume method. The geometries and mesh of the computational domain was discretized through a three-dimensional composed of approximately 1,950,000 elements. The results show that the optimal geometries that presented resistance to the flow up to 30 times smaller than the other configurations. In addition, it was possible to verify that the Hess-Murray Law is not always valid, since the system adapts its optimal geometry to each flow condition, in order to provide a better flow architecture to meet the objective of minimizing resistance to flow in agreement with the Constructal Law. This work advanced in the present state of the art, since it developed a three-dimensional model without simplifications, applied to the fluid flow system in T-shaped structure using the Construtal Design method, validating the analytical results presented in the bibliography and presenting new references that allow increase the complexity of flow systems as well as the implementation of more advanced optimization methods.

Este trabalho tem como propósito investigar como a variação geométrica de determinados parâmetros envolvidos na construção de uma geometria bifurcada de seção circular, em forma de Y, afeta a resistência ao escoamento, tanto de fluidos newtonianos como não newtonianos. As geometrias estudadas foram construídas utilizando-se o princípio do Design Construtal. Os parâmetros variados foram a relação entre os comprimentos dos dutos pais e filhos, a relação entre os diâmetros dos mesmos dutos, e o ângulo central da estrutura em forma de Y. Para as relações geométricas lineares foram utilizados os valores de 0,5; 0,6; 0,7; 0,8; 0,9 e 1, enquanto para os ângulos, foram utilizados os valores de 155°, 135°, 115°, 95°, 75°, 45°, 25° e 10°. Os fluidos utilizados foram do tipo newtoniano e não newtoniano, dentre estes últimos, foram estudados fluidos pseudoplásticos e dilatantes. O trabalho foi realizado através de simulações numéricas, implementadas com a utilização do software comercial Ansys Fluent, o qual resolve as equações governantes através do método dos volumes finitos. As malhas utilizadas foram do tipo poliédrica. Os resultados indicam que há uma diferença em relação ao que se espera da literatura para as relações entre os diâmetros e os comprimentos. A Lei Hess-Murray indica que estas relações ótimas seriam de 2-1/3 para as relações entre os diâmetros e comprimentos. No presente trabalho, foram determinadas relações entre os diâmetros próximas de 0,6, e entre os comprimentos, iguais a 1. Os ângulos ótimos ficaram localizados no intervalo entre 100° e 135°. / This work aims to investigate how the geometric variation of certain parameters involved in the construction of a bifurcated Y-shaped circular cross-section geometry affects the flow resistance of both Newtonian and non-Newtonian fluids. The geometries studied were constructed using the Constructal Design principle. The parameters were the relationship between the lengths of the daughter and parent ducts, the relationship between the diameters of the same ducts, and the central angle of the Y-shaped structure. For the linear geometric relations, values of 0.5; 0.6; 0.7; 0.8; 0.9 and 1 where used, for the angles, the values of 155 °, 135 °, 115 °, 95°, 75 °, 45 °, 25 ° and 10 ° were used. The fluids used were of the Newtonian and non-Newtonian type, among the latter, pseudo plastic and dilatant fluids were studied. The work was carried out through numerical simulations, implemented with the commercial software Ansys Fluent, which solves the governing equations through the finite volume method. The meshes used were of the polyhedral type. The results indicate that there is a difference in relation to what is expected from the literature for the relationships between diameters and lengths. The Hess-Murray Law indicates that these optimal relations would be 2-1/3 for the relationships between diameters and lengths. In the present work, relationships between the diameters close to 0,6 were found and s equal to 1 between the lengths. The optimum angles were located in the range between 100 ° and 135 °.

O presente trabalho tem como foco principal os aspectos físicos da camada limite de um corpo rombudo. Os dados experimentais foram obtidos através da técnica de anemometria à Laser e os dados numéricos foram obtidos através de Simulação Numérica Direta (DNS). O modelo é composto de uma semi-elipse seguida de uma base reta, este modelo foi extensivamente estudado pelo Prof. P. Bearman nos anos 60, por esta razão será aqui chamado de corpo de Bearman. A razão de esbeltez (AR) é definida por c/h, onde c é a corda do modelo e h é a altura da base. Por se tratar de uma geometria com ponto de separação fixo, torna-se possível o estudo da influência da camada limite nos parâmetros médios da esteira. Os experimentos foram realizados entre os números de Reynolds de 600 e 3000 e as razões de esbeltez estudadas são AR=3, AR=5 e AR=7. Com o auxílio de simulações DNS, desenvolveu-se uma metodologia para definir o momento exato do desprendimento em função da variação dos parâmetros médios, tais como vazão instantânea e pressão na base. Na parte experimental, mostrou-se que o aumento da espessura da camada limite com o aumento da razão de esbeltez acarreta na diminuição da frequência de desprendimento. Em geral, os resultados numéricos e experimental mostraram boa concordância exceto o comprimento de formação dos vórtices.

Todos os cursos de água carream sedimentos com maior ou menor intensidade, sendo responsáveis pela redução da calha de rios e volume de reservatórios. O assoreamento dos reservatórios reduz sua vida útil, comprometendo os sistemas de geração de energia, de navegação, de irrigação, gerando problemas de ordem ambiental. Neste contexto, subsídios ao processo de gerenciamento dos reservatórios pela predição do assoreamento, tornam-se necessários. O reservatório Represa Velha, situado na bacia hidrográfica do Mogi Guaçu, foi escolhido para estudo de caso pela disponibilidade de dados batimétricos, de vazão e de sedimentos. Com os recursos da dinâmica dos fluidos computacional (DFC), é utilizado o programa MIKE 21C. A solução das equações hidrodinâmicas (Saint-Venant 2D) integradas na vertical, tiveram como condição de contorno a vazão constante em cada período da medição, na entrada do reservatório, e a altura de água constante, para cada dado de vazão, no canal de saída. A calibração do modelo numérico foi feita para os dados temporais de vazão, concentração de sedimentos e batimetrias do reservatório. Com dados da primeira batimetria e campanhas batimétricas sucessivas, foi realizada a comparação entre os resultados obtidos pela simulação numérica e aqueles medidos em campo. As dificuldades e limitações encontradas neste processo são comentadas. / All rivers transport sediments with bigger or smaller intensity, being responsible for the reduction of their channel and volume of reservoirs. The sedimentation of reservoirs reduces their useful life, compromising energy generation systems, navigation, irrigation, and generating problems of environmental order. In this context, subsidies to the process of reservoirs management for the prediction of sedimentation become necessary. The reservoir Represa Velha, located in the Mogi Guaçu\'s watershed, was chosen as a case study because of its available data of bathymetry, flow and sediment discharge. With the resources of computational fluids dynamics (CFD), the program MIKE 21C was used. In the solution of the problem are used the hydrodynamic equations (Saint-Venant 2D) integrated in the vertical. The boundary conditions are the constant flow, in each period of the measurement, at the entrance of the reservoir, and the constant water level, for each flow measured, at the exit channel. The calibration of the numeric model was developed with data from previews work, of flow, sediment concentration and reservoir bathymetry. The numerical data were compared with the measured sedimentation. The difficulties and limitations found in this process are commented on.

O presente trabalho apresenta resultados de medidas de campo, laboratório e simulações numéricas para o diagnóstico de um sítio de descarte de resíduos de re-refino de óleo lubrificante, localizado no Município de Ribeirão Preto, no Estado de São Paulo. A metodologia laboratorial utilizada mostrou-se viável, sendo descritos também casos nos quais a primeira proposta de trabalho não foi adequada. Adicionalmente, o programa numérico comercial escolhido para simular a situação do sítio mostrou-se suficientemente manuseável e adequado para estudos expeditos. Os dados de contaminação da água subterrânea por metais mostram que os principais contaminantes são o chumbo e o alumínio. Em termos de contaminação do solo por HPAs, verificou-se, de forma geral, baixa concentração (em alguns casos não detectáveis), ficando o fenantreno como o composto com a concentração mais elevada, porém sem padrão de comparação. O naftaleno, que possui esse padrão, mostrou concentração muito inferior àquela considerada como referência. Os resultados da simulação numérica são apresentados em sete camadas ao longo da profundidade, sendo sua resolução, no plano, efetuada com malhas de 50,0 m por 50,0 m. Simulou-se uma pluma poluente de chumbo, com evolução para um tempo de 10 anos, com concentração inicial e de contorno compatível com as medidas efetuadas nos poços de monitoramento (25 &#956g/L). Os resultados apontam para contaminação efetiva da área utilizada para disposição do resíduo, bem como alteração nas circunvizinhanças. O presente trabalho representa uma contribuição para a sistematização da metodologia de estudo de áreas similares. / This research introduces results of field and laboratory applications and numerical simulation to assess a site of oil waste disposal, in Ribeirão Preto, São Paulo, Brazil. The laboratorial methodology applied was useful, but there were also cases where the first try was not adequate. Additionally, the commercial numeric software selected to simulate site conditions was properly handled and applicable to expedited site assessment. The results of metals concentration in ground water highlight lead and aluminum as the main contaminants. In terms of PAHs in soil, in general, low concentrations (in some cases not detectible) were observed, where phenanthrene was the compound with the highest concentration, however without reference level. The naphthalene, has reference level, present concentrations much lower than the reference level. The numerical simulation results are presented in seven subsurface layers, which in a plan view the grid was divided in 50,0 m per 50,0 m. A lead contamination was simulated in the transport model, with 10 years of temporal discretization, and with initial and boundary concentrations compatible with measures in monitoring wells (25 &#956g/L). The results show contamination in the disposal area, and alteration in the neighboring areas. The present research represents a contribution to the methodology for similar site assessment.

O objetivo desse trabalho foi avaliar o processo de formação do cavaco durante o torneamento utilizando simulação numérica pelo método dos elementos finitos. Para realizar o estudo foram definidos dois tipos de aços inoxidáveis austeníticos, um com matriz metálica sem a presença significativa de inclusões, do tipo ABNT 304, e outro com a presença de inclusões não metálicas, do tipo ABNT 303. O estudo foi focado nos mecanismos de formação e ruptura do cavaco, na determinação das forças de usinagem, no campo de tensões, deformações, e temperaturas durante o processo, que foram relacionados com aspectos e características da microestrutura do material. Os resultados obtidos foram comparados com as forças de usinagem experimentais, com a espessura e morfologia do cavaco. O desenvolvimento do trabalho, de acordo com a metodologia adotada, foi realizado em diferentes etapas. Inicialmente foi elaborado e aplicado um modelo de simulação da usinagem considerando o material homogêneo. Em outra etapa, foi realizada a modelagem de uma microestrutura submetida a um estado de tensão e deformação semelhante ao encontrado na simulação da usinagem realizada com material homogêneo. Os resultados mostraram que as partículas das inclusões maiores, alongadas, e em maior quantidade aumentam a tensão e a deformação na microestrutura. As elevadas temperaturas obtidas na usinagem dos aços inoxidáveis austeníticos aumentam a ductilidade dos sulfetos, esses se deformam em compressão junto com a matriz, e têm um efeito limitado como agente de redução dos esforços de usinagem. Por outro lado, os sulfetos facilitam a etapa de ruptura do cavaco em tensões trativas, e tendem a se romper facilitando o processo de quebra. / The aim of this study was to evaluate the chip formation process during turning using a numerical simulation by finite element method. Two different austenitic stainless steel were evaluated, one without a significant presence of inclusions, ABNT 304, and another one with the presence of non-metallic inclusions, ABNT 303. The study was focused on the chip formation and rupture, in determining the machining forces, stresses, strains and temperatures during the process, which were related to aspects and characteristics of the microstructure. The results obtained were compared with the experimental machining forces, thickness and morphology of the chip. The development of the work, according to the methodology, was carried out at different steps. The first one implemented a simulation model considering a homogeneous material. In the next step, the microstructure was subjected to a state of stress similar to that found in the simulation of the machining performed with homogeneous material. The results showed that the largest inclusion particles, elongated and in greater quantity, increase stresses and strains in the microstructure. The high temperatures obtained when machining austenitic stainless steel increase the ductility of sulfides. These are deformed in compression along with the matrix, and have a limited effect in reducing machining forces. On the other hand, sulfides facilitate the chip fracture at tensile stresses.