Cette thèse est consacrée à l'étude de méthodes numériques pour la simulation des écoulements polyphasiques en milieu poreux, en vue de leur application à des problèmes d'ingénierie pétrolière ou environnementale. Nous présentons une formulation générique du modèle d'écoulements à nombre quelconque de composants présents dans un nombre quelconque de phases. Dans notre approche l'approximation des flux diffusifs (issus, par exemple, de la loi de Darcy) s'appuie sur de nouveaux schémas, appelés schémas gradient, qui ont plusieurs avantages sur les schémas industriels standard : ces derniers, qui sont des schémas volumes finis multi-points centrés aux mailles, ne sont généralement pas symétriques et convergent difficilement sur des cas à forts rapports d'anisotropie. Nous montrons en revanche que les schémas gradient conduisent naturellement à des approximations symétriques et convergentes. Parmi cette classe de schémas, nous étudions plus particulièrement le schéma "VAG" qui fait intervenir des inconnues au centre des mailles et aux sommets du maillage. Ce schéma conduit à la définition de flux entre le centre d'une maille et ses sommets, qui sont utilisés pour généraliser la méthode "VAG" au contexte polyphasique. Des tests numériques montrent alors que ce schéma est robuste, et conduit à un très bon compromis précision/coût, ce qui en fait un candidat idoine pour les applications industrielles. Nous présentons notamment un cas test, basé sur des observations de terrains, d'injection et de dissolution de CO2 dans la région proche d'un puits foré dans un aquifère salin. Nous montrons alors que le schéma numérique permet de simuler l'assèchement et la précipitation de minéral observée en pratique. Un chapitre de la thèse est enfin consacré à l'étude pratique et théorique d'une méthode numérique générique pour contrôler l'effet d'axe lors de l'utilisation de schémas industriels / This thesis is focused on numerical methods dedicated to the simulation of multiphase flow in porous media, involved in petroleum or environmental engineering. We present a generic formulation of the flow model which is able to take into account any number of components within any number of phases. In our approach the approximation of the diffusive fluxes (mainly resulting from Darcy's law) is based on new schemes, called gradient schemes, which show several advantages over the standard industrial numerical schemes : these schemes, which belong to the class of the cell-centred MultiPoint Flux Approximation finite volume schemes, are not symmetric and may lead to difficulties of convergence in the case of high anisotropy ratios. We indeed show that gradient schemes are naturally providing symmetric and convergent approximations. We particularly study the Vertex Approximate Gradient scheme (called the VAG scheme), whose discrete unknowns are the values at the cell centres and at the vertices of the mesh. This scheme implies the definition of fluxes between the centre of a given cell and its vertices, used for the extension of the scheme to multiphase flow. Numerical tests show the robustness and the accuracy of the method for a low computational cost, which enables the use of the VAG scheme in an industrial framework. A test case, based on experimental data of injection and dissolution of CO2 in the near-well region within a saline aquifer, shows the aptitude of the scheme for reproducing drying and salt precipitation, which are practically observed. Finally, a chapter of the thesis is devoted to the theoretical and practical study of a general numerical method for controlling Grid Orientation Effect in industrial simulators

Ce travail porte sur la résolution numérique des transferts couplés par rayonnement et conduction au sein d'un milieu semi-transparent pour une géométrie 3D de forme complexe. Le rayonnement thermique est simulé par un code de calcul développé durant cette thèse. Ce code résout l'équation du transfert radiatif (ETR) par une méthode aux volumes finis (MVF) avec une formulation " cell-vertex " s'appliquant à des maillages tétraédriques non structurés. Il utilise un schéma de fermeture de type exponentiel, un ordre de parcours ainsi qu'une résolution matricielle innovante pour la MVF appliquée à l'ETR. Le modèle mis en place traite des milieux absorbants, émettants, gris ou non-gris bordés par des surfaces noires ou opaques à réflexion diffuse. Le couplage rayonnement-conduction s'effectue sur le même maillage avec un code d'éléments finis pour la conduction. La validation du code de rayonnement et du couplage passe par de nombreux cas tests issus de la littérature. Il aborde les milieux gris, isotherme avec différentes géométries où les effets de la discrétisation spatiale et angulaire sont observés au travers d'une étude de sensibilité. Trois schémas de fermeture ont été étudiés sur un milieu transparent pour montrer leurs influences sur la précision et la diffusion numérique. Les études des transferts de chaleur couplés traitent le problème de l'équilibre radiatif et du couplage conduction-rayonnement en régime stationnaire ou instationnaire avec les équations adimensionnées. La dernière étude porte sur un milieu non-gris tel que le verre en considérant la conduction et le rayonnement en régime stationnaire avec une méthode spectrale par bande pour la partie radiative / This work deals with the numerical solution of coupled radiative and conductive heat transfer in participating media in complex 3D geometries. Thermal radiation is simulated by a numerical code developed during this thesis. This code solves the radiative transfer equation (RTE) by a modified finite volume method (FVM) with a cell-vertex formulation applied to unstructured tetrahedral meshes. It uses a closure relation based on an exponential scheme, a marching order map and an innovative matrix solution for the FVM applied to the RTE. The model is applied to absorbing-emitting, grey or non-grey media bounded by black or opaque walls with diffuse reflection. The mesh used for the radiation-conduction coupling is the one used by the finite element code for the conduction. The validation of the radiative code and the coupling are carried out through several test cases taken from the literature. Grey and isothermal media with different geometries are considered, and the effects of the spatial and angular discretizations are observed through a sensitivity study. Three closure schemes have been studied on a transparent medium in order to show their influence on the accuracy and false scattering. Studies of coupled heat transfer are carried out on radiative equilibrium problems and coupled radiation-conduction problems in steady or transient states with the dimensionless equations. Finally a non-grey medium such glass is also studied, considering conduction and radiation in steady state with a spectral band model for radiation

Notre travail de thèse présente une contribution à l'étude du comportement d'un jet flottant horizontal, représentatif de la dispersion de rejets dans un milieu marin. Il consiste à modéliser ce type d'écoulement par une approche mathématique basée sur la résolution numérique moyennant la méthode des volumes finis, à valider le modèle numérique par des mesures à échelle réduite sur des maquettes expérimentales, et, enfin, à simuler la dispersion de polluants à grande échelle sur un cas réel. Trois cas d’étude ont été abordés:Le premier cas est relatif à l’étude numérique et expérimentale d'un jet flottant turbulent rond et "non-Boussinesq", injecté horizontalement dans un milieu statique et homogène. Les résultats ont permis de décrire la nature du jet et son comportement en fonction des conditions initiales d’éjection. La deuxième étude a concerné un autre cas représentatif des rejets des émissaires de stations d'épuration. Il s’agit d’un jet pariétal admettant l’approximation de Boussinesq en régime de convection mixte. L’objectif ici est d’étudier l'influence de l'effet combiné de la turbulence et de la présence de la paroi sur le comportement du jet. Dans le dernier cas d’étude, on a modélisé, en 2D et en 3D, à grande échelle un processus côtier de dispersion de rejets en surface libre appliqué sur le cas de la pollution de la baie de Tanger. Les résultats ont permis de visualiser le mécanisme de la dispersion et d’avoir des informations précieuses sur l’écoulement généré au voisinage des plages par l’interaction des rejets et des mouvements de flux et de reflux de la marée. / This work is a contribution to the study of horizontal buoyant jet behavior that presents the dispersion of discharges into the marine environment. It consists in the modeling of this flow by a mathematical approach based on numerical simulation by means of the finite volume method, the validation of a numerical model by measurements on experimental model at a small, and, finally, the simulation of pollutant dispersion on a large scale on a real case. Three cases of study were broached: The first case relates to the experimental and numerical study of horizontal round turbulent non-Boussinesq buoyant jet in a static homogeneous environment. The results were used to describe the nature and the behavior of the jet as a function of the initial conditions of ejection. The second study involved another case of discharges from outfalls. We investigated a numerical and experimental study about a horizontal buoyant wall turbulent jet in a static homogeneous environment. The aim was to analyze the influence the effect of both turbulence and wall boundary on the behavior of the jet. The latter case of study focused on numerical simulation in 2D and 3D of the coastal process of discharges dispersion on a free surface. This modeling dealt more precisely with the dispersion of discharges into a marine environment in the presence of cross flows. The bay of Tangier in Morocco was chosen as an application site. The results made it possible to visualize the dispersion mechanism and to gain valuable information on the flow generated by the interaction of discharges and high/low tide movements near the beaches of the bay.

Malgré l'efficacité informatique et la faible dissipation numérique de la méthode de Boltzmann sur réseau (LBM) classique reposant sur un algorithme de propagation-collision, cette méthode est limitée aux maillages cartésiens uniformes. L'adaptation de l'étape de discrétisation à différentes échelles de la mécanique des fluides est généralement réalisée par des schémas LBM à échelles multiples, dans lesquels le domaine de calcul est décomposé en plusieurs sous-domaines uniformes avec différentes résolutions spatiales et temporelles. Pour des raisons de connectivité, le facteur de résolution des sous-domaines adjacents doit être un multiple de deux, introduisant un changement abrupt des échelles spatio-temporelles aux interfaces. Cette spécificité peut déclencher des instabilités numériques et produire des sources de bruit parasite rendant l'exploitation de simulations à finalités aéroacoustiques impossible. Dans la présente thèse, nous avons d'abord élucidé le sujet du raffinement de maillage dans la LBM classique en soulignant les défis et les sources potentielles d'erreur. Par la suite, une méthode de Boltzmann sur réseau hybride (HLBM) est proposée, combinant l'algorithme de propagation-collision avec un algorithme de flux au sens eulérien obtenu à partir d'une discrétisation en volumes finis des équations de Boltzmann à vitesse discrète. La HLBM combine à la fois les avantages de la LBM classique et une flexibilité géométrique accrue. La HLBM permet d'utiliser des maillages cartésiens non-uniformes. La validation de la méthode hybride sur des cas tests 2D à finalité aéroacoustique montre qu'une telle approche constitue une alternative viable aux schémas Boltzmann sur réseau à échelles multiples, permettant de réaliser des raffinements locaux en H. Enfin, un couplage original, basé sur l'algorithme de propagation-collision et une formulation isotherme des équations de Navier-Stokes en volumes finis, est proposé. Une telle tentative présente l'avantage de réduire le nombre d'équations du solveur volumes finis tout en augmentant la stabilité numérique de celui-ci, en raison d'une condition CFL plus favorable. Les deux solveurs sont couplés dans l'espace des moments, où la solution macroscopique du solveur Navier-Stokes est injectée dans l'algorithme de propagation-collision à l'aide de la collision des moments centrés. La faisabilité d'un tel couplage est démontrée sur des cas tests 2D, et les résultas obtenus sont comparés avec la HLBM. / Despite the inherent efficiency and low dissipative behaviour of the standard lattice Boltzmann method (LBM) relying on a two step stream and collide algorithm, a major drawback of this approach is the restriction to uniform Cartesian grids. The adaptation of the discretization step to varying fluid dynamic scales is usually achieved by multi-scale lattice Boltzmann schemes, in which the computational domain is decomposed into multiple uniform subdomains with different spatial resolutions. For the sake of connectivity, the resolution factor of adjacent subdomains has to be a multiple of two, introducing an abrupt change of the space-time discretization step at the interface that is prone to trigger instabilites and generate spurious noise sources that contaminate the expected physical pressure signal. In the present PhD thesis, we first elucidate the subject of mesh refinement in the standard lattice Boltzmann method and point out challenges and potential sources of error. Subsequently, we propose a novel hybrid lattice Boltzmann method (HLBM) that combines the stream and collide algorithm with an Eulerian flux-balance algorithm that is obtained from a finite-volume discretization of the discrete velocity Boltzmann equations. The interest of a hybrid lattice Boltzmann method is the pairing of efficiency and low numerical dissipation with an increase in geometrical flexibility. The HLBM allows for non-uniform grids. In the scope of 2D periodic test cases, it is shown that such an approach constitutes a valuable alternative to multi-scale lattice Boltzmann schemes by allowing local mesh refinement of type H. The HLBM properly resolves aerodynamics and aeroacoustics in the interface regions. A further part of the presented work examines the coupling of the stream and collide algorithm with a finite-volume formulation of the isothermal Navier-Stokes equations. Such an attempt bears the advantages that the number of equations of the finite-volume solver is reduced. In addition, the stability is increased due to a more favorable CFL condition. A major difference to the pairing of two kinetic schemes is the coupling in moment space. Here, a novel technique is presented to inject the macroscopic solution of the Navier-Stokes solver into the stream and collide algorithm using a central moment collision. First results on 2D tests cases show that such an algorithm is stable and feasible. Numerical results are compared with those of the previous HLBM.

Dans ce travail, nous nous intéressons à deux problèmes de la théorie des systèmes hyperboliques faisant intervenir des interfaces. Le premier concerne des modèles de couplages entre un fluide compressible et une particule ponctuelle et le second concerne la capture numérique précise des chocs, ces discontinuités qui apparaissent dans les solutions des systèmes hyperboliques.Sur la première thématique, nous commençons par introduire les différents modèles, dans lesquels la particule et le fluide interagissent à travers une force de frottement qui tend à rapprocher leurs vitesses. Le couplage est singulier car il fait intervenir le produit d’une fonction discontinue par une mesure de Dirac. On peut toutefois définir précisément le système en voyant la particule comme une interface à travers laquelle des relations liant les propriétés du fluide et celle de la particule sont imposées. Lorsque le fluide suit une équation de Burgers, nous démontrons la convergence d’une classe de schéma numérique, et nous obtenons l’existence d’une solution au problème de Cauchy pour une donnée initiale à variation totale bornée. Dans le cas plus complexe où le fluide est décrit par les équa- tions d’Euler isothermes, on prouve l’existence et l’unicité d’une solution autosemblable au problème de Riemann lorsque la particule est immobile. Des simulations numériques sont également présentées.La dernière partie de la thèse est consacrée à la construction de schémas non diffusifs pour les systèmes hyperboliques. Ces schémas, de type volumes finis, sont construits pour être exact lorsque la donnée initiale est un choc isolé. Ils sont basé sur une reconstruction discontinue de la solution au début de chaque itération en temps, dans le but de reconstituer des chocs à l’intérieur de certaines cellules du maillage. Cette stratégie mène à des schémas très peu diffusifs qui, lorsque l’opérateur de reconstruction est bien choisi, approchent correctement les solutions de cas tests problématiques (chocs lents, chocs forts, réflexions pour la dynamique des gaz, chocs non classiques pour les systèmes qui ne sont pas vraiment non linéaires). / In this work, we study two problems concerning hyperbolic systems involving interfaces. The first one concerns the study of models of coupling between a compressible fluid and a pointwise particle. The second one deals with the sharp numerical approximation of shocks, which are discontinuities that appear in the solutions of hyperbolic systems.In the first two parts of the manuscript, we introduce different models of fluid-particle couplings. The fluid and the particle interact on each other through a drag force, which brings their velocities closer to one another. The coupling is singular because it can be written as the product of a discontinuous function by a Dirac measure. However, the system can be precisely defined as follows. The particle is seen as an interface through which interface conditions linking the properties of the fluid with those of the particle are imposed. When the fluid follows the compressible Burgers equations, we prove the convergence of a family of finite volume schemes and obtain the existence of a solution when the initial data has total bounded variation. In the more difficult case where the fluid is described by the isothermal Euler equations, we prove the existence and uniqueness of a selfsimilar solution to the Riemann problem, when the particle is motionless. Numerical experiments are also presented.In the last part of this work, we build non diffusive numerical schemes for different hyperbolic systems. These finite volume schemes are built to be exact when the initial data is an isolated shock. They are based on a discontinuous reconstruction of the solution at the beginning of each time step, in order to reconstruct shocks inside some specific cells of the mesh. The schemes we present have a very low numerical diffusion and, when the reconstruction operator is well chosen, they are able to correctly approximate the solution on various problematic test cases. These cases include slowly moving shocks, strong shocks and shock reflections for gas dynamics, as well as the apparition of nonclassical shocks for systems that are not truely nonlinear.

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.

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.

A medição de vazão de escoamentos multifásicos é uma necessidade constante em diversas atividades industriais como exploração de óleo e gás, controle de linhas de transporte de vapor e monitoramento de sistemas de resfriamento de usinas nucleares. Dentre os meios disponíveis para a realização da medição de vazão mássica, os dispositivos de medição de pressão diferencial constituem um dos métodos mais simples, sendo sua construção, aplicação e operação em escoamentos monofásicos bem conhecidas e definidas por normas técnicas. No entanto, sua aplicação tem sido estendida a escoamentos multifásicos, geralmente estando aliada a uma técnica adicional de medição de fração de vazio ou fração volumétrica das fases. Este trabalho descreve o estudo numérico de escoamentos multifásicos através de medidores de vazão baseados em pressão diferencial como placas de orifício e bocais de vazão de raio longo. Para tal, primeiramente foram conduzidas simulações de escoamentos monofásicos através de placas de orifício e bocais de vazão de raio longo na faixa de número de Reynolds 15.000 500.000. Os resultados de coeficiente de descarga obtidos foram quantitativamente comparados com os valores preditos por norma ISO, apresentando desvio máximo de aproximadamente 4, 9% para as placas e de 1,0% para os bocais. Em uma segunda etapa, escoamentos do tipo gás úmido (wet gas) através de placas de orifício foram simulados através de três abordagens diferentes. Os resultados de vazão mássica total obtidos foram comparados com dados experimentais fornecidos pela PETROBRAS. As abordagens que consideram o escorregamento entre as fases apresentaram previsões mais próximas dos experimentos, com desvio relativo médio de 3,9%, enquanto a modelagem homogênea apresentou um desvio médio de 6, 6%. Nestes estudos, foram também avaliadas as estruturas desenvolvidas no escoamento através de visualizações da distribuição de fases. São também apresentadas duas sugestões para complementação da caracterização de um escoamento multifásico: (1) a introdução da informação de fração de vazio na formulação apresentada por Paz (2011) e (2) a análise estatística do sinal de pressão diferencial em placas de orifício. Com relação ao primeiro item, comparações quantitativas com dados experimentais sugeriram que a alternativa apresentada é viável para operações de monitoramento da produção. Já o último estudo mostrou qualitativamente a influência da quantidade de líquido na flutuação da pressão diferencial / The flowrate measurement of multiphase flows is a constant need at many industrial activities such as oil and gas exploration, steam transport lines control and monitoring of nuclear plants cooling systems. Within the available means for performing flowrate measurement, the differential pressure devices constitute one of the simplest methods, with their construction, application and operation in single phase flows being well known and defined by technical standards. However, their application has been extended to multiphase flows, usually being allied to a void fraction or phase volume fraction measurement technique. This work describes a numerical study of multiphase flows through differential pressure-based flowrate meters such as orifice plates and long radius nozzles. Firstly simulations of single phase flows through orifice plates and long radius nozzles were conducted in the Reynolds number range 15.000500.000. The obtained results of discharge coefficients were quantitatively compared to ISO Standard predicted values, showing a maximum deviation of approximately 4,9% for the orifice plates and of 1,0% for the nozzles. In a second stage, wet gas flows through orifice plates were simulated by means of three approaches. The calculated results of total mass flowrate were compared to experimental data provided by PETROBRAS. The approaches that considered the slip between phases provided the closest results to the experiments, with a mean relative error of 3, 9%, while the homogeneous modeling presented an error of 6, 6%. In these studies, the structures developed within the domain were also evaluated through the visualization of the phases distribution. Two suggestions for complementing the characterization of a multiphase flow are presented: (1) the introduction of void fraction information into the formulation presented by Paz (2011) and (2) the statistical analysis of the orifice plate pressure drop signal. Regarding the first item, quantitative comparison with experimental data suggested that the presented alternative is viable for production monitoring operations. The last study qualitatively revealed the influence of the liquid loading in the pressure drop fluctuation.

O objetivo desse trabalho consiste no estudo de métodos de identificação de parâmetros em problemas envolvendo a equação de advecção-difusão 2D. Essa equação é resolvida utilizando o método dos volumes finitos, sendo empregada métodos de reconstrução de alta ordem em malhas não-estruturadas de triângulos para calcular os fluxos nas faces dos volumes de controle. Como ferramenta de busca dos parâmetros é empregada a técnica baseadas em gradientes, sendo os mesmos calculados utilizando processos baseados em métodos adjuntos. / The aim of this work concern to study parameter identification methods on problems involving the advection-diffusion equation in two dimensions. This equation is solved employing the finite volume methods, and high-order reconstruction methods, on triangle unstructured meshes to solve the fluxes across the faces of control volumes. As parameter searching tool is employed technicals based on gradients. The gradients are solved using processes based on adjoint methods.

O objetivo desse trabalho consiste no estudo de métodos de identificação de parâmetros em problemas envolvendo a equação de advecção-difusão 2D. Essa equação é resolvida utilizando o método dos volumes finitos, sendo empregada métodos de reconstrução de alta ordem em malhas não-estruturadas de triângulos para calcular os fluxos nas faces dos volumes de controle. Como ferramenta de busca dos parâmetros é empregada a técnica baseadas em gradientes, sendo os mesmos calculados utilizando processos baseados em métodos adjuntos. / The aim of this work concern to study parameter identification methods on problems involving the advection-diffusion equation in two dimensions. This equation is solved employing the finite volume methods, and high-order reconstruction methods, on triangle unstructured meshes to solve the fluxes across the faces of control volumes. As parameter searching tool is employed technicals based on gradients. The gradients are solved using processes based on adjoint methods.