Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / The models based on spatially distributed processes describe erosion and sediment yield using a physical approach of the mechanisms. In general, these models require high quality input data, description of the landscape, insertion of the physical processes involved, among others. Thus, it is not uncommon to find models that due to some operational difficulty do not fit well in the natural processes. Therefore, the objective of this study was to understand the hydrological and sedimentological dynamics of the Arvorezinha/RS experimental catchment in order to set the constraints of mathematical modeling that generate the overestimation in modeling the suspended sediment concentration (SSC) and sediment yield (SY). The study took place in the Arvorezinha experimental catchment, located in southern Brazil, which possesses a total drainage area of 1.23km2. Along with hidrossedimentological monitoring, there was an insertion of solute analysis through dissolved silicon (Sid) and particulate during and between precipitation events. This was carried out to understand the hydrological behavior of the catchment as well as assist in flow separation. The mathematical modelling of the erosive process was done with the LISEM model for the nine events monitored with Sid. In addition, we conducted an analysis of extreme events and variables responses, especially for SY, between 2010 and 2015. The years of 2014 and 2015 were the wettest since 2010, which had around 2.500mm of rain. The events monitored with Sid led us to conclude that the sub-surface flows have the biggest contribution during the event of precipitation, especially in events that take place in autumn-winter and low-average magnitude. Furthermore, it was possible to observe the effect of dilution that occurs during the events, that is, as liquid discharge increases, the concentration of Sid reduces. This is reflected in the suspended sediment concentration, because, as the event occurs, the SSC resulting from the surface flow will be diluted in the subsurface flow. Because, according to the present study, these water movements with lower speeds have a very important role in the hydrograph (it was higher than or close to 50% of the total volume). Therefore, the modeling of the erosive process is significantly influenced by the contribution of the slow flows (subsurface and groundwater flow), given that the model simulates the quick flow (superficial). In this way, LISEM model usage in areas with lower flow speeds can be limited; consequently, a greater level of care is required when concluding on the use of the model and their capacity to represent the hidrossedimentological processes. / Os modelos baseados em processos espacialmente distribuídos descrevem a erosão e a produção de sedimentos utilizando uma abordagem física dos mecanismos. Em geral, esses modelos exigem alta qualidade nos dados de entrada, descrição da paisagem, inserção dos processos físicos envolvidos, entre outros. Assim, não é raro encontrarmos modelos que, por alguma dificuldade operacional, não se ajustem bem aos processos naturais. Dessa forma, o objetivo do estudo foi compreender a dinâmica hidrológica e sedimentológica da bacia experimental de Arvorezinha/RS, para definir os condicionantes da modelagem matemática que geram a superestimava na modelagem da concentração de sedimentos em suspensão e da produção de sedimentos (PS). O estudo ocorreu na bacia experimental de Arvorezinha, sul do Brasil, com área de drenagem de 1,23km2. Juntamente ao monitoramento hidrossedimentológico, houve a inserção da análise de solutos por meio do silício dissolvido (Sid) e particulado durante eventos e entre eventos de precipitação. O monitoramento com Sid foi realizado para compreender o comportamento hidrológico da bacia hidrográfica, bem como auxiliar na separação dos escoamentos. A modelagem matemática do processo erosivo foi realizada com o modelo LISEM para os nove eventos monitorados com Sid. Além disso, foi realizada uma análise dos eventos extremos e das variáveis respostas, principalmente a PS, entre 2010 e 2015. Os anos de 2014 e 2015 foram os mais chuvosos desde 2010, com, aproximadamente, 2.500mm ao ano. Os eventos monitorados com Sid permitiram concluir que os escoamentos subsuperficiais possuem a maior contribuição de água no hidrograma durante o evento de precipitação, principalmente em eventos que ocorrem no outono-inverno e de baixa-média magnitude. Além disso, pôde-se verificar o efeito de diluição que ocorre durante os eventos, isto é, com o aumento da descarga líquida, a concentração de Sid reduz. Isso reflete-se na concentração de sedimentos em suspensão (Css), pois, havendo o evento, a Css de origem do escoamento superficial será diluída nos escoamentos subsuperficiais, uma vez que, conforme o estudo aqui realizado, esses movimentos de água mais lentos apresentam um papel muito importante no hidrograma (foi superior ou próximo a 50% do volume total). Portanto, a modelagem do processo erosivo na bacia experimental de Arvorezinha é significativamente influenciada pela contribuição dos fluxos lentos (subsuperficial e subterrâneo), uma vez que o modelo simula o escoamento rápido (superficial). Desse modo, o uso do modelo LISEM em áreas com a presença de escoamentos de menor velocidade pode ser limitado; por isso, é necessário cuidado ao obter conclusões sobre o uso do modelo e a capacidade do mesmo em representar os processos hidrossedimentológicos.

Este trabalho apresenta um método numérico para simular escoamentos viscoelásticos bidimensionais governados pela equação constitutiva Giesekus [Schleiniger e Weinacht 1991]. As equações governantes são resolvidas pelo método de diferenças finitas numa malha deslocada. A superfície livre do fluido é modelada por partículas marcadoras possibilitando assim a sua visualização e localização. O cálculo da velocidade é efetuado por um método implícito enquanto a pressão é calculada por um método explícito. A equação constitutiva de Giesekus é resolvida pelo método de Euler modificado explícito. O método numérico desenvolvido nesse trabalho é verificado comparando-se a solução numérica com a solução analítica para o escoamento de um fluido Giesekus em um canal. Resultados de convergência são obtidos pelo uso de refinamento de malha. Os resultados alcançados incluem um estudo da aplicação do modelo Giesekus para simular o escoamento numa contração planar 4:1 e o problema de um jato incidindo sobre uma placa rígida, em que o fenômeno jet buckling é simulado. / This work presents a numerical method to simulate two-dimensional viscoelastic flows governed by the Giesekus constitutive equation [Schleiniger e Weinacht 1991]. The governing equations are solved by the finite difference method on a staggered grid. The free surface of the fluid is modeled by tracer particles thus enabling its visualization and location. The calculation of the velocity is performed by an implicit method while pressure is calculated by an explicit method. The Giesekus constitutive equation is resolved by the explicit modified Euler method. The numerical method developed in this work is verified by comparing the numerical solution with the analytical solution for the flow of a Giesekus fluid in a channel. Convergence results are obtained by the use of mesh refinement. Results obtained include a study of the application of the Giesekus model to simulate the flow through a 4:1 contraction and the problem of a jet flowing onto a rigid plate where the phenomenon of jet buckling is simulated.

Este trabalho apresenta um método numérico para simular escoamentos viscoelásticos tridimensionais com superfícies livres governados pela equação constitutiva de Maxwell. O método numérico é uma extensão da técnica newtoniana GENSMAC3D para escoamentos viscoelásticos. As equações governantes para escoamentos incompressíveis cartesianos isotérmicos são apresentadas em detalhes. O tratamento do tensor não-newtoniano em contornos rígidos tridimensionais é apresentado em detalhes, bem como o cálculo da condição de contorno na superfície livre. As equações governantes são resolvidas pelo método de diferenças finitas numa malha deslocada tridimensional. O fluido é modelado pela técnica das partículas marcadoras e tratado como uma superfície linear por partes. O método numérico desenvolvido foi implementado no sistema Freeflow3D, e resultados numéricos obtidos na simulação de escoamentos tridimensionais governados pela equação constitutiva de Maxwell são apresentados. Adicionalmente, apresentamos uma validação mostrando a convergência do método desenvolvido nesse trabalho / This work presents a numerical method for solving three-dimensional viscoelastic flows with free surfaces governed by the Maxwell constitutive equation. The numerical method is an extension of the Newtonian technique GENSMAC3D to viscoelastic flows. The governing equations for Cartesian incompressible isothermic flows are presented in details. The treatment of the non-Newtonian tensor on three-dimensional rigid boundaries is given in details as well as the calculation of the boundary conditions on the free surface. The governing equations are solved by a finite difference method using a three-dimensional staggered grid. The fluid is described by marker particles and is represented by a piecewise linear surface. The numerical method developed in this work was implemented into the Freeflow3D system, and numerical results obtained from the simulation of complex three-dimensional flows are presented. Additionally, we present validation results and demonstrate the convergence of the method by performing mesh refinement

[pt] A motivação deste trabalho consiste no uso do four-roll mill para aumentar a separação de fases de emulsões água em óleo (A/O) presente no processamento primário da indústria de petróleo. A partir da conservação de massa e momento, a fase contínua foi modelada como escoamento incompressível, bi-dimensional e isotérmico. Simulações numéricas utilizando o método de elementos finitos foram implementadas para revelar a influência das diversas configurações de escoamento no comportamento mecânico do material. A partir dos resultados obtidos, a habitual forma de classificar o escoamento no four-roll mill de acordo com a literatura se demonstrou ineficiente. Este trabalho sugere classificações locais de escoamento a cada posição dependendo se a mesma está ocupada pela fase contínua ou dispersa da emulsão. O efeito da fase dispersa é descrito via pós-processamento. Microelementos no formato de vetores foram inseridos no domínio e investigou-se suas deformações e trajetórias. Consequentemente, analisou-se a deformação de gotas e a sua respectiva influência na instabilidade da emulsão. / [en] The motivation of this work consists in the use of four-roll mill in order to increase the phase separation of water-in-oil emulsions (W/O) present in the primary process of oil industry. With mass and momentum conservation, the continuous phase is modeled by an incompressible, bi-dimensional and isothermal flow. Numerical simulations employing the finite element method were implemented to reveal the influence of the several flow configurations in the material mechanical behavior. From the obtained results, the standard way of classifying the flow in the four-roll mill according to the literature was proved inefficient. This work suggests local flow classifications for each position depending if it is occupied by the continuous or dispersed phase. The effect of the dispersed phase was described by a post-processing scheme. Microelements in shape of vectors were inserted in the domain and their deformations and pathlines were investigated. Thus, the deformation of droplets and their respective influences in the emulsion instability were analyzed.

O foco principal do presente trabalho é estender uma metodologia numérica embasada no uso de uma técnica de refinamento adaptativo de malha (AMR - Adaptive Mesh Refinement) e no uso de esquemas temporais multipasso implícitos-explícitos (IMEX) a aplicações envolvendo escoamentos reativos com baixo número de Mach. Originalmente desenvolvida para escoamentos incompressíveis, a formulação euleriana daquela metodologia emprega as equações de Navier-Stokes como modelo matemático para descrever a dinâmica do escoamento e o Método da Projeção, baseado no divergente nulo da velocidade do escoamento, para tratar o acoplamento pressão-velocidade presente na formulação com variáveis primitivas. Tal formulação euleriana original é estendida para acomodar novas equações agregadas ao modelo matemático da fase contínua: conservação de massa, fração de mistura (para representar as concentrações de combustível e oxidante), e energia. Além disso, uma equação termodinâmica de estado é integrada ao modelo matemático estendido e é empregada juntamente com a equação de conservação de massa para produzir uma nova restrição (não nula desta vez) ao divergente do campo de velocidade. Assume-se que o escoamento ocorre a baixo número de Mach (hipótese principal). O Método de Diferença Finita é empregado na discretização espacial das variáveis eulerianas de estado, empregando-se uma malha AMR. As vantagens e dificuldades desta extensão são cuidadosamente investigadas e reportadas. Pela importância, do ponto de vista de aplicações práticas, alguns estudos numéricos preliminares envolvendo escoamentos incompressíveis turbulentos com sprays são realizados (as gotículas compõem a fase dispersa). Num primeiro momento, apenas sprays com gotículas inertes são considerados. Embora ainda apenas iniciais, tais estudos já se mostram importantes pois identificam com clareza, em primeira instância, algumas das dificuldades inerentes a serem enfrentadas ao se tratar dentro desta nova metodologia um conjunto relativamente grande de gotículas lagrangianas. No caso de escoamentos incompressíveis turbulentos com sprays, a integração temporal se dá com métodos IMEX para a fase contínua e com o Método de Euler Modificado para a fase dispersa. A turbulência, em todos os casos que a envolvem, é tratada pelo modelo de Simulação das Grandes Escalas (LES - Large Eddy Simulation). As simulações computacionais se dão em um domínio tridimensional, um parelelepípedo, e empregam uma extensão (resultante do presente trabalho) do código AMR3D, um programa de computador sequencial implementado em Fortran90, oriundo de uma colaboração de longa data entre o IME-USP e o MFLab/FEMEC-UFU (Laboratório de Dinâmica de Fluidos da Universidade Federal de Uberlândia). O processamento foi efetuado no LabMAP (Laboratório da Matemática Aplicada do IME-USP). / It is the main goal of the present work to extend a numerical methodology based on both the use of an adaptive mesh refinement technique (AMR) and the use of a multistep, implicit-explicit time-step strategy (IMEX) to applications involving low Mach number reactive flows. Originally developed for incompressible flows, the Eulerian formulation of that methodology employs the Navier-Stokes equations to model the flow dynamics and the Projection Method, based on the vanishing divergence of the velocity field, to tackle the pressure-velocity coupling present when using primitive variables. That Eulerian formulation is extended by adding a new set of equations to the original mathematical model, describing the various properties of the continuous phase: mass conservation, mixture fraction (to represent concentrations of fuel and oxidizer) and energy. Also, a thermodynamic equation of state is included into the extended mathematical model which is employed, along with the equation for the conservation of mass, to derive a new restriction (this time, different from zero) to the divergence of the velocity field. It is assumed that one is dealing with a low Mach number flow (the main hipothesis). The discretization in space employs the Finite Difference Method for the Eulerian variables on a AMR mesh. Advantages and difficulties of such an extension of the previous methodology are carefully investigated and reported. For its importance in the real-world applications, few preliminary numerical studies involving incompressible turbulent flows with sprays are performed (the droplets form what it is called the dispersed phase). Only sprays formed by inert droplets are considered. Even though initial yet, such studies are most important because they clearly identify, first hand, certain difficulties in handling relatively large sets of Lagrangian droplets in the context of this new AMR methodology. In the context of turbulent incompressible flows with sprays, the overall time-step scheme is given by IMEX methods for the continuous phase and by the Improved Euler Method for the dispersed phase. In all the cases in which it is considered, turbulence is modeled by the Large Eddy Simulation (LES) model. The computational simulations are held in a tridimensional domain given by a paralellepiped and all of them employ the extention (resulting of the present work) of the AMR3D code, a sequencial computer program implemented in Fortran90, whose origin is the collaborative work between IMEUSP and MFLab/FEMEC-UFU (Fluid Dynamics Laboratory, Federal University of Uberlândia). Computations were performed at LabMAP (Applied Mathematics Laboratory at IME-USP).

O estudo de métodos de modelagem de escoamentos aerodinâmicos em regime transônico é de grande importância para a engenharia aeronáutica. O maior desafio no tratamento desses escoamentos está na sua característica não linear devido aos efeitos de compressibilidade e formação de ondas de choque. Tais efeitos não lineares influenciam no desempenho de superfícies aerodinâmicas em geral, bem como são responsáveis pelo aparecimento de fenômenos danosos para a resposta aeroelástica de aeronaves. O equacionamento para esses tipos de escoamentos pode ser obtido via as equações básicas da mecânica dos fluidos. No entanto, apenas soluções numéricas de tais equações são possíveis de ser obtidas de forma prática no presente momento. Para o caso específico do tratamento de problemas transônicos, as equações de Euler formam um conjunto de equações diferenciais a derivadas parciais capazes de capturar os efeitos não lineares de escoamentos compressíveis, porém os efeitos da viscosidade não são levados em consideração. O objetivo desse trabalho é implementar uma rotina computacional capaz de resolver numericamente escoamentos em regime transônico em torno de aerofólios. Para isso as equações de Euler não lineares são utilizadas e o campo de fluido ao redor de um perfil aerodinâmico é discretizado pelo método das diferenças finitas. Uma malha estruturada do tipo C discretizando o fluido ao redor de um aerofólio NACA0012 é considerada. A metodologia para solução numérica é baseada no método explícito de MacCormack de segunda ordem de precisão no tempo e espaço. Baseados na aproximação upwind, termos de dissipação artificial com coeficientes não lineares também são adicionados ao método. A solução do escoamento transônico estacionário ao redor do aerofólio NACA0012 é obtida e as principais propriedades do escoamento são apresentadas. Observa-se a formação de ondas de choque através de contornos de número de Mach ao redor do aerofólio. Gráficos das distribuições de pressão no intra e extradorso do aerofólio são mostrados, onde se identificam aos efeitos da brusca variação de pressão devido as ondas de choque. Os resultados são validados com valores de distribuição de pressão para o mesmo aerofólio encontradas na literatura técnica. Os resultados obtidos combinam bem com os fornecidos em códigos computacionais para solução do mesmo problema aerodinâmico / The study of aerodynamic modeling methods for the transonic flow regime is of great importance in aeronautical engineering. Major challenge on the treatment of those flows is on their nonlinear features due to compressibility effects and shock waves (appearance). Such nonlinear effects present a strong influence on aerodynamic performance, as well as they are responsible for harmful aeroelastic response phenomena in aircraft. Equations for transonic flows can be obtained from the basic fluid mechanic equations. However, only numerical methods are able to attain practical solutions for those set of differential equations at the present moment. For the specific case of treating transonic flow problems, the nonlinear Euler equations provide a set of partial differential equations with features to capture nonlinear effects of typical compressible flows, despite of not accounting for viscous flows effects. The aim of this work is to implement a computational routine for the numerical solution of transonic flows around airfoils. The Euler equations are used and the flow field around a aerodynamic profile is discretized by finite difference method. A C-type structured mesh is used to discretize the flow around a NACA0012 airfoil. The methodology for numerical solution is based on the explicit MacCormack method which has second order accuracy in time and space. Based on the upwind approximation, artificial dissipation with nonlinear coefficients is incorporated to the method. The steady transonic flow around the NACA0012 airfoil numerical solution is assessed and the main flow properties are presented. Shock wave structure can also be observed by means of the Mach number contours around the airfoil. Pressure distributions on upper and lower surfaces for different flow conditions are also shown, thereby allowing the observation of the effects of the abrupt pressure change due to shock waves. The results are validated using data presented in the technical literature. The present code solutions agree well with the solution obtained in other computational codes used for the same problem

Este trabalho apresenta uma técnica de diferenças finitas para resolver a equação constitutiva Fluido de Segunda Ordem para escoamentos tridimensionais com superfície livre. As equações governantes são resolvidas pelo método de diferenças finitas em uma malha deslocada 3D. A superfície livre é modelada por células marcadoras (Marker-and-Cell) e as condições de contorno a superfície livre são empregadas. O método numérico apresentado neste trabalho foi validado pela comparação entre as soluções numéricas obtidas para o escoamento em um tubo com a solução analítica correspondente para Fluidos de Segunda Ordem. Ao fazer refinamento de malha, a convergência do método numérico foi verificada. Resultados numéricos da simulação do problema do inchamento do extrudado para números de Deborah De \'< OU =\' 0:3 são apresentados / This work presents a finite difference method to simulate three-dimensional viscoelastic flow with free surfaces governed by the constitutive equation Second Order Fluid. The governing equations are solved by the finite difference method in a three-dimensional shifted mesh. The free surface of fluid is modeled by the Marker-and-Cell method which allows for the visualization and the location of the free surface of fluid. The full free surface stress conditions are employed. The numerical method developed in this work is validated by comparing the numerical and analytic solutions for the steady state flow of a Second Order Fluid in a pipe. By using mesh refinement convergence results are given. Numerical results of the simulation of the transient extrudate swell of a Second Order Fluid of the Deborah number De \'< OR =\' 0:3 are presented

Diversas aplicações industriais utilizam escoamentos de fluidos viscoelásticos, e em muitos casos é necessário saber se os escoamentos propagam-se no estado laminar ou no turbulento. Embora a hidrodinâmica de fluidos viscoelásticos seja fortemente afetada pelo balanço entre forças inerciais e elásticas no escoamento, o efeito da elasticidade sobre a estabilidade de escoamentos inerciais não foi completamente estabelecido. Neste trabalho, estuda-se o que ocorre durante a transição laminar-turbulenta, investigando a convecção de ondas de Tollmien-Schlichting para o escoamento incompressível, para um fluido viscoelástico, entre placas paralelas, utilizando a equação constitutiva Giesekus. Para isto, adotou-se a simulação numérica direta para verificar a estabilidade dos escoamentos à perturbações não estacionárias deste fluido. Experimentos computacionais para verificação do código foram realizados. Com os resultados numéricos obtidos, foi possível verificar e analizar a estabilidade de escoamentos utilizando-se o modelo não newtoniano Giesekus. / Several industrial applications use viscoelastic fluid flows, and it is necessary to know if the flows propagate in the laminar or turbulent state. Although the hydrodynamics of viscoelastic fluids is strongly affected by the balance between inertial and elastic forces in the flow, the effect of elasticity on the stability of inertial flows has not been completely established. In this work we study what happens during the laminar-turbulent transition, investigating the convection of Tollmien-Schlichting waves for the incompressible flow, for a viscoelastic fluid, between parallel plates, using the constitutive equation Giesekus. For this, the direct numerical simulation was used to verify the stability of the flows to the non-stationary perturbations of this fluid. Computational experiments to verify the code were performed. With the numerical results obtained, it was possible to verify and analyze the stability of flows modelled by Giesekus non-newtonian model.

O estudo de métodos de modelagem de escoamentos aerodinâmicos em regime transônico é de grande importância para a engenharia aeronáutica. O maior desafio no tratamento desses escoamentos está na sua característica não linear devido aos efeitos de compressibilidade e formação de ondas de choque. Tais efeitos não lineares influenciam no desempenho de superfícies aerodinâmicas em geral, bem como são responsáveis pelo aparecimento de fenômenos danosos para a resposta aeroelástica de aeronaves. O equacionamento para esses tipos de escoamentos pode ser obtido via as equações básicas da mecânica dos fluidos. No entanto, apenas soluções numéricas de tais equações são possíveis de ser obtidas de forma prática no presente momento. Para o caso específico do tratamento de problemas transônicos, as equações de Euler formam um conjunto de equações diferenciais a derivadas parciais capazes de capturar os efeitos não lineares de escoamentos compressíveis, porém os efeitos da viscosidade não são levados em consideração. O objetivo desse trabalho é implementar uma rotina computacional capaz de resolver numericamente escoamentos em regime transônico em torno de aerofólios. Para isso as equações de Euler não lineares são utilizadas e o campo de fluido ao redor de um perfil aerodinâmico é discretizado pelo método das diferenças finitas. Uma malha estruturada do tipo C discretizando o fluido ao redor de um aerofólio NACA0012 é considerada. A metodologia para solução numérica é baseada no método explícito de MacCormack de segunda ordem de precisão no tempo e espaço. Baseados na aproximação upwind, termos de dissipação artificial com coeficientes não lineares também são adicionados ao método. A solução do escoamento transônico estacionário ao redor do aerofólio NACA0012 é obtida e as principais propriedades do escoamento são apresentadas. Observa-se a formação de ondas de choque através de contornos de número de Mach ao redor do aerofólio. Gráficos das distribuições de pressão no intra e extradorso do aerofólio são mostrados, onde se identificam aos efeitos da brusca variação de pressão devido as ondas de choque. Os resultados são validados com valores de distribuição de pressão para o mesmo aerofólio encontradas na literatura técnica. Os resultados obtidos combinam bem com os fornecidos em códigos computacionais para solução do mesmo problema aerodinâmico / The study of aerodynamic modeling methods for the transonic flow regime is of great importance in aeronautical engineering. Major challenge on the treatment of those flows is on their nonlinear features due to compressibility effects and shock waves (appearance). Such nonlinear effects present a strong influence on aerodynamic performance, as well as they are responsible for harmful aeroelastic response phenomena in aircraft. Equations for transonic flows can be obtained from the basic fluid mechanic equations. However, only numerical methods are able to attain practical solutions for those set of differential equations at the present moment. For the specific case of treating transonic flow problems, the nonlinear Euler equations provide a set of partial differential equations with features to capture nonlinear effects of typical compressible flows, despite of not accounting for viscous flows effects. The aim of this work is to implement a computational routine for the numerical solution of transonic flows around airfoils. The Euler equations are used and the flow field around a aerodynamic profile is discretized by finite difference method. A C-type structured mesh is used to discretize the flow around a NACA0012 airfoil. The methodology for numerical solution is based on the explicit MacCormack method which has second order accuracy in time and space. Based on the upwind approximation, artificial dissipation with nonlinear coefficients is incorporated to the method. The steady transonic flow around the NACA0012 airfoil numerical solution is assessed and the main flow properties are presented. Shock wave structure can also be observed by means of the Mach number contours around the airfoil. Pressure distributions on upper and lower surfaces for different flow conditions are also shown, thereby allowing the observation of the effects of the abrupt pressure change due to shock waves. The results are validated using data presented in the technical literature. The present code solutions agree well with the solution obtained in other computational codes used for the same problem

O estudo do comportamento de fluidos é um antigo domínio das ciências da natureza. Ultimamente, fenômenos de engenharia que eram estudados empiricamente passaram a ser estudados com auxílio computacional. A Dinâmica de Fluidos Computacional (DFC) é a área da ciência da computação que estuda métodos computacionais para simulação de escoamento de fluidos, e muitas vezes é a forma mais prática, ou a única, de se observar fenômenos de interesse no escoamento. Este projeto de Mestrado procurou investigar, no âmbito da simulação de um escoamento bifásico, métodos computacionais para representar a interface entre dois fluidos imiscíveis. A separação dos fluidos por meio de uma interface é necessária para assegurar que, propriedades como viscosidade e densidade, específicas de cada fluido, sejam utilizadas corretamente para o cálculo do movimento de seus respectivos fluidos. Desenvolvemos um método lagrangeano sem a utilização de malhas com o objetivo de suprir algumas restrições de trabalhos prévios. Para representar a interface entre os dois fluidos, este método utiliza uma técnica de reconstrução de superfícies baseada em aproximações de superfícies algébricas de alta ordem. Os resultados numéricos reportados neste documento evidenciam o potencial da nossa abordagem / The study of the behaviour of fluids is an ancient field in natural sciences. Recently, engineering phenomena that were empirically studied started to be done with computacional aid. The Computational Fluid Dynamics (CFD) is the area of science that studies computational methods for computer simulation of fluid flow, and often is the most practical way, or the only, to observe phenomena of interest in flow. This Masters degree project sought to investigate, in the context of the simulation of biphasic flows, computational methods to represent the interface between two immiscible fluids. The separation of fluids by the means of an interface is required to ensure that, during the simulation, the physical properties of a fluid, like density and viscosity (specific of each fluid) are properly used in the calculus of the respective fluid motion. We developed a lagrangean method without the use of mesh with the goal of alleviating some of the previous works restrictions. To represent the interface between the two fluids, this method uses a surface reconstruction technique based on approximations of high order algebraic surfaces. The numerical results reported herein show the potential of our approach