Cette thèse porte sur la résolution numérique du modèle multi-dimensionnel d'écoulement cisaillé en eau peu profonde. Dans le cas d'un mouvement unidimensionnel, ces équations coïncident avec les équations de la dynamique de gaz pour un choix particulier de l'équation d'état. Dans le cas multi-dimensionnel, le système est complètement différent du modèle de la dynamique de gaz. Il s'agit d'un système EDP hyperbolique 2D non-conservatif qui rappelle un modèle de turbulence barotrope. Le modèle comporte trois types d'ondes correspondant à la propagation des ondes de surface, des ondes de cisaillement et à celle de la discontinuité de contact. Nous présentons dans le cas 2D un schéma numérique basé sur une nouvelle approche de ``splitting" pour les systèmes d'équations non-conservatives. Chaque sous-système ne contient qu'une seule famille d'ondes: ondes de surface ou ondes de cisaillement, et discontinuité de contact. La précision d'une telle approche est testée sur des solutions exactes 2D décrivant l'écoulement lorsque la vitesse est linéaire par rapport aux variables spatiales, ainsi que sur des solutions décrivant des trains de rouleaux 1D. Finalement, nous modélisons un ressaut hydraulique circulaire formé dans un écoulement convergent radial d'eau. Les résultats numériques obtenus sont clairement similaires à ceux obtenus expérimentalement: oscillations du ressaut et son rotation avec formation du point singulier. L'ensemble des validations proposées dans ce manuscrit démontre les aptitudes du modèle et de la méthode numérique pour la résolution des problèmes complexes d'écoulements cisaillés en eau peu profonde multidimensionnels. / This thesis is devoted to the numerical modelling of multi-dimensional shear shallow water flows. In 1D case, the corresponding equations coincide with the equations describing non--isentropic gas flows with a special equation of state. However, in the multi-D case, the system differs significantly from the gas dynamics model. This is a 2D hyperbolic non-conservative system of equations which is reminiscent of a generic Reynolds averaged model of barotropic turbulent flows. The model has three families of characteristics corresponding to the propagation of surface waves, shear waves and average flow (contact characteristics). First, we show the ability of the one-dimensional conservative shear shallow water model to predict the formation of roll-waves from unstable initial data. The stability of roll waves is also studied.Second, we present in 2D case a new numerical scheme based on a splitting approach for non-conservative systems of equations. Each split subsystem contains only one family of waves (either surface or shear waves) and contact characteristics. The accuracy of such an approach is tested on exact 2D solutions describing the flow where the velocity is linear with respect to the space variables, and on the solutions describing 1D roll waves. Finally, we model a circular hydraulic jump formed in a convergent radial flow of water. Obtained numerical results are qualitatively similar to those observed experimentally: oscillation of the hydraulic jump and its rotation with formation of a singular point. These validations demonstrate the capability of the model and numerical method to solve challenging multi--dimensional problems of shear shallow water flows.

The numerics of hyperbolic conservation laws, e.g., the Euler equations of gas dynamics, shallow water equations and MHD equations, is non-trivial due to the convective terms being highly non-linear and equations being coupled. Many numerical methods have been developed to solve these equations, out of which central schemes and upwind schemes (such as Flux Vector Splitting methods, Riemann solvers, Kinetic Theory based Schemes, Relaxation Schemes etc.) are well known. The majority of the above mentioned schemes give rise to very dissipative solutions. In this thesis, we propose novel low dissipative numerical algorithms for some hyperbolic conservation laws representing fluid flows. Four different and independent numerical methods which give low diffusive solutions are developed and demonstrated. The first idea is to regulate the numerical diffusion in the existing dissipative schemes so that the smearing of solution is reduced. A diffusion regulator model is developed and used along with the existing methods, resulting in crisper shock solutions at almost no added computational cost. The diffusion regulator is a function of jump in Mach number across the interface of the finite volume and the average Mach number across the surface. The introduction of the diffusion regulator makes the diffusive parent schemes to be very accurate and the steady contact discontinuities are captured exactly. The model is demonstrated in improving the diffusive Local Lax-Friedrichs (LLF) (or Rusanov) method and a Kinetic Scheme. Even when employed together with accurate methods of Roe and Osher, improvement in solutions is demonstrated for multidimensional problems. The second method, a Central Upwind-Biased Scheme (CUBS), attempts to reorganize a central scheme such that information from irrelevant directions is largely reduced and the upwind biased information is retained. The diffusion co-efficient follows a new format unlike the use of maximum characteristic speed in the Local Lax-Friedrichs method and the scheme results in improved solutions of the flow features. The grid-aligned steady contacts are captured exactly with the reorganized format of diffusion co-efficient. The stability and positivity of the scheme are discussed and the procedure is demonstrated for its ability to capture all the features of solution for different flow problems. Another method proposed in this thesis, a Central Rankine-Hugoniot Solver, attempts to integrate more physics into the discretization procedure by enforcing a simplified Rankine-Hugoniot condition which describes the jumps and hence resolves steady discontinuities very accurately. Three different variants of the scheme, termed as the Method of Optimal Viscosity for Enhanced Resolution of Shocks (MOVERS), based on a single wave (MOVERS-1), multiple waves (MOVERS-n) and limiter based diffusion (MOVERS-L) are presented. The scheme is demonstrated for scalar Burgers equation and systems of conservation laws like Euler equations, ideal Magneto-hydrodynamics equations and shallow water equations. The new scheme uniformly improves the solutions of the Local Lax-Friedrichs scheme on which it is based and captures steady discontinuities either exactly or very accurately. A Grid-Free Central Solver, which does not require a grid structure but operates on any random distribution of points, is presented. The grid-free scheme is generic in discretization of spatial derivatives with the location of the mid-point between a point and its neighbor being used to define a relevant coefficient of numerical dissipation. A new central scheme based on convective-pressure splitting to solve for mid-point flux is proposed and many test problems are solved effectively. The Rankine-Hugoniot Solver, which is developed in this thesis, is also implemented in the grid-free framework and its utility is demonstrated. The numerical methods presented are solved in a finite volume framework, except for the Grid-Free Central Solver which is a generalized finite difference method. The algorithms developed are tested on problems represented by different systems of equations and for a wide variety of flow features. The methods presented in this thesis do not need any eigen-structure and complicated flux splittings, but can still capture discontinuities very accurately (sometimes exactly, when aligned with the grid lines), yielding low dissipative solutions. The thesis ends with a highlight on the importance of developing genuinely multidimensional schemes to obtain accurate solutions for multidimensional flows. The requirement of simpler discretization framework for such schemes is emphasized in order to match the efficacy of the popular dimensional splitting schemes.

Fundação Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro / O processo de recuperação secundária de petróleo é comumente realizado com a injeção de água no reservatório a fim de manter a pressão necessária para sua extração. Para que o investimento seja viável, os gastos com a extração têm de ser menores do que o retorno financeiro obtido com o petróleo. Para tanto, tornam-se extremamente importantes as simulações dos processos de extração. Neste trabalho são estudados os problemas de Burgers e de Buckley-Leverett visando o escoamento imiscível água-óleo em meios porosos, onde o escoamento é incompressível e os efeitos difusivos (devido à pressão capilar) são desprezados. Com o objetivo de incorporar conhecimento matemático mais avançado, para em seguida utilizá-lo no entendimento do problema estudado, abordou-se com razoável profundidade a teoria das leis de conservação. Foram consideradas soluções fracas que, fisicamente, podem ser interpretadas como ondas de choque ou rarefações, então, para que fossem distinguidas as fisicamente admissíveis, foi utilizado o princípio de entropia, nas suas diversas formas. Inicialmente consideramos alguns exemplos clássicos de métodos numéricos para uma lei de conservação escalar, os quais podem ser vistos como esquemas conservativos de três pontos. Entre eles, o método de Lax-Friedrichs (LF) e o método de Lax-Wendroff (LW). Em seguida, um esquema composto foi testado, o qual inclui na sua formulação os métodos LF e LW (chamado de LWLF-4). Respeitando a condição CFL, foram obtidas soluções numéricas de todos os problemas tratados aqui. Com o objetivo de validar tais soluções, foram utilizadas soluções analíticas oriundas dos problemas de Burgers e Buckley- Leverett. Também foi feita uma comparação com os métodos do tipo TVDs com limitadores de fluxo, obtendo resultado satisfatório. Vale à pena ressaltar que o esquema LWLF-4, pelo que nos consta, nunca foi antes utilizado nas resoluções das equações de Burgers e Buckley- Leverett. / The secondary recovery of petroleum is usually performed with injection of water through an oil reservoir to keep the oil pressure for the exploration. In order to make the exploration profitable, the extraction cost must be less than the financial return, which means that the simulation of the exploration process is extremely relevant. In this work, the Burgers- and- Buckley-Leverett problems are studied seeking a two-phase displacement in porous media. The flow is considered incompressible and capillary effects are ignored. In order to analyze the problem, it was necessary to use the theory of conservation law in a spatial variable. Weak solutions, which can be understood as shock or rarefaction waves, are studied with the entropy condition, so that only the physically correct solutions are considered. Some classical numerical methods, which can be seen as conservative schemes of three points, are studied, among them the Lax-Friedrichs (LF) and Lax-Wendroff (LW) methods. A composite scheme, called LWLF-k, is tested using LF and LW methods, being respected the CFL condition, with satisfactory results. In order to validate the numerical schemes, we consider analytical solutions of the Burgers-and-Buckley-Leverett equations. Was also made a comparison with TVDs methods with flux limiters, obtaining satisfactory results. We emphasize that to the best of our knowledge, the LWLF-4 scheme has never been used to solve the Buckley-Leverett equation.

Fundação Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro / O processo de recuperação secundária de petróleo é comumente realizado com a injeção de água no reservatório a fim de manter a pressão necessária para sua extração. Para que o investimento seja viável, os gastos com a extração têm de ser menores do que o retorno financeiro obtido com o petróleo. Para tanto, tornam-se extremamente importantes as simulações dos processos de extração. Neste trabalho são estudados os problemas de Burgers e de Buckley-Leverett visando o escoamento imiscível água-óleo em meios porosos, onde o escoamento é incompressível e os efeitos difusivos (devido à pressão capilar) são desprezados. Com o objetivo de incorporar conhecimento matemático mais avançado, para em seguida utilizá-lo no entendimento do problema estudado, abordou-se com razoável profundidade a teoria das leis de conservação. Foram consideradas soluções fracas que, fisicamente, podem ser interpretadas como ondas de choque ou rarefações, então, para que fossem distinguidas as fisicamente admissíveis, foi utilizado o princípio de entropia, nas suas diversas formas. Inicialmente consideramos alguns exemplos clássicos de métodos numéricos para uma lei de conservação escalar, os quais podem ser vistos como esquemas conservativos de três pontos. Entre eles, o método de Lax-Friedrichs (LF) e o método de Lax-Wendroff (LW). Em seguida, um esquema composto foi testado, o qual inclui na sua formulação os métodos LF e LW (chamado de LWLF-4). Respeitando a condição CFL, foram obtidas soluções numéricas de todos os problemas tratados aqui. Com o objetivo de validar tais soluções, foram utilizadas soluções analíticas oriundas dos problemas de Burgers e Buckley- Leverett. Também foi feita uma comparação com os métodos do tipo TVDs com limitadores de fluxo, obtendo resultado satisfatório. Vale à pena ressaltar que o esquema LWLF-4, pelo que nos consta, nunca foi antes utilizado nas resoluções das equações de Burgers e Buckley- Leverett. / The secondary recovery of petroleum is usually performed with injection of water through an oil reservoir to keep the oil pressure for the exploration. In order to make the exploration profitable, the extraction cost must be less than the financial return, which means that the simulation of the exploration process is extremely relevant. In this work, the Burgers- and- Buckley-Leverett problems are studied seeking a two-phase displacement in porous media. The flow is considered incompressible and capillary effects are ignored. In order to analyze the problem, it was necessary to use the theory of conservation law in a spatial variable. Weak solutions, which can be understood as shock or rarefaction waves, are studied with the entropy condition, so that only the physically correct solutions are considered. Some classical numerical methods, which can be seen as conservative schemes of three points, are studied, among them the Lax-Friedrichs (LF) and Lax-Wendroff (LW) methods. A composite scheme, called LWLF-k, is tested using LF and LW methods, being respected the CFL condition, with satisfactory results. In order to validate the numerical schemes, we consider analytical solutions of the Burgers-and-Buckley-Leverett equations. Was also made a comparison with TVDs methods with flux limiters, obtaining satisfactory results. We emphasize that to the best of our knowledge, the LWLF-4 scheme has never been used to solve the Buckley-Leverett equation.