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.

Nous présentons une contribution à la modélisation tridimensionnelle des phénomènes électromagnétiques et thermiques couplés dans les matériaux supraconducteurs à haute température critique. La méthode des volumes finis est adoptée comme méthode de résolution des équations aux dérivées partielles caractéristiques aux phénomènes physiques traités. Le couplage électromagnétique thermique est assuré par un algorithme alterné. L’ensemble des modèles mathématico-numériques ainsi développés et implémentés sous Matlab, sont appliqués pour étudié le comportement des supraconducteurs dans le cadre des paliers magnétiques et pour étudier le comportement des supraconducteurs durant le processus d’aimantation. Les résultats à caractère magnétique et ceux à caractère thermique sont largement présentés. La validité du travail proposé est atteinte par comparaison des résultats ainsi obtenus à ceux donnés par l’expérimentation. / We present a contribution for three-dimensional modeling of coupled electromagnetic and thermal phenomena in high temperature superconductor. The control volume method is used for the resolution of the partial derivative equations characterising of the treated physical phenomena. The electromagnetic and thermal coupling is ensured by an alternate algorithm. All mathematical and numerical models thus developed and implemented in Matlab software, are used for the simulation. The results in magnetic term and those in thermal term are largely presented. The validity of the suggested work is reached by the comparison of the results so obtained to those given by the experiment.

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.

La corrosion des aciers dans les ouvrages en béton armé due aux chlorures est l’une des principales causes de dégradation des constructions notamment en façade maritime. La prédiction de la durée de vie des constructions vis-à-vis de cette pathologie, nécessite la compréhension des principaux mécanismes liés au transfert de chlorures surtout dans les matériaux insaturés (soumis au marnage ou au brouillard salin). Pour améliorer les travaux qui existent dans la littérature sur le sujet, les principaux objectifs de cette thèse sont : la caractérisation expérimentale du phénomène de double couche électrique (DCE) et l’étude de son influence sur le transfert ionique, l’étude des effets de la température et de l’hystérésis (sorption-désorption) sur le profil de chlorures, l’étude du couplage entre corrosion des aciers et transfert des chlorures. Pour atteindre ces objectifs, le travail est divisé en quatre parties. La première est consacrée à l’approche théorique du problème pour établir les équations permettant de décrire le couplage entre transferts ionique et hygrothermique dans les matériaux insaturés. La prise en compte de la DCE est également explicitée. La deuxième partie est expérimentale, elle est réalisée sur des éco-bétons fabriqués à base de laitier et de cendres volantes. Elle consiste à quantifier les paramètres d’entrée du modèle identifiés lors de la première partie. Il s’agit : de la porosité, de la perméabilité aux gaz, du coefficient de diffusion d’ions chlorures, de la distribution porale, des isothermes d’interactions chimiques et des isothermes de sorption-désorption et du potentiel Zeta caractérisant la DCE. Dans la troisième partie, les équations de transfert traduisant le modèle sont résolues numériquement. L’étude de sensibilité du modèle vis à vis de la DCE, de la température, de la perméabilité à l’eau est réalisée. Une confrontation des résultats du modèle (profil de chlorures) avec des résultats expérimentaux obtenus en situation de marnage a permis de montrer l’intérêt de considérer l’hystérésis des courbes de sorption-désorption. Quant à l’effet de la DCE il se manifeste pour potentiel Zeta dépassant 25mV. La dernière partie concerne l’étude du couplage entre transfert de chlorures et corrosion des aciers. Les estimations de la durée d’incubation, de la perte de section d’acier, les temps requis pour l’apparition des premières fissures et l’éclatement du béton d’enrobage ont été réalisées. Nous avons montré que l’utilisation de la loi de Fick peut sous-estimer ces grandeurs. / The corrosion of steels in the reinforced concrete structures caused by chlorides is one of the main causes of degradation of the constructions in particular in maritime facade. The prediction of the service life of the constructions towards this pathology requires the understanding of the main mechanisms relating to the chlorides transport especially in the unsaturated materials submitted to the tide or to the salt spray. To improve the works of the literature, the main objectives should be accomplished are : to study the experimental characterization of the electrical double layer (EDL) and to study its influence on the ionic transport ; to research on the temperature effect and the hysteresis effect on the chloride profile ; to examine chloride transfer coupled with corrosion of reinforcement. To reach these objectives, the work is divided into four parts. The first one is dedicated to the theoretical approach of the problem to establish the equations allowing describing the coupling between ionic and hydrothermal transfer to the unsaturated materials. The consideration of the EDL is also clarified. The second part is experiment. It is realized on eco-concretes made on base of slag and fly ash. It consists in quantifying the inputs of the model identified during the first part such as: porosity, gas permeability, chloride diffusion coefficients, pore size distribution, chemical interactions, isotherms of sorption-desorption and the Zeta potential characterizing the EDL. In the third part, the transfer equations translating to the model are numerically resolved. The study of model sensibility according to EDL, the temperature, the water permeability is performed. A comparison between the numerical results (chloride profile) obtained with the model with experimental ones obtained in the case of tidal zone allowed shows the interest to consider the hysteresis effects. Also, it is shown that the EDL effect is significant beyond Zeta potential of 25mV. The last part concerns the study of the coupling between chloride transfer and corrosion of steels. The estimations of the duration of incubation, the loss of section of steel, the times required to cover cracking were examined. The result showed that the use of the law of Fick can underestimate these parameters.

Nous étudions plusieurs aspects d'équations aux dérivées partielles multi-échelles. Pour trois exemples, la présence de multiples échelles, spatiales ou temporelles, motive un travail de modélisation mathématique ou constitue un enjeu de discrétisation. La première partie est consacrée à la construction et l'étude d'un système multicouche de type Saint-Venant pour décrire un fluide à surface libre (océan). Son obtention s'appuie sur l'analyse des échelles spatiales, précisément l'hypothèse « eau peu profonde ». Nous justifions nos équations à partir du modèle primitif et montrons un résultat d'existence locale de solution. Puis nous proposons un schéma volumes finis et des simulations numériques. Nous étudions ensuite un problème hyperbolique de relaxation, inspiré de la théorie cinétique des gaz. Nous construisons un schéma numérique via une stratégie préservant l'asymptotique : nous montrons sa convergence pour toute valeur du paramètre de relaxation, ainsi que sa consistance avec le problème à l'équilibre local. Des estimations d'erreurs sont établies et des simulations numériques sont présentées. Enfin, nous étudions un problème d'écoulement sanguin dans une artère avec stent, modélisé par un système de Stokes dans un domaine contenant une petite rugosité périodique (géométrie double échelle). Pour éviter une discrétisation coûteuse du domaine rugueux (l'artère stentée), nous formulons un ansatz de développement de la solution type Chapman-Enskog, et obtenons une loi de paroi implicite sur le bord du domaine lisse (artère seule). Nous montrons des estimations d'erreurs et des simulations numériques / This work is concerned with different aspects of multiscale partial differential equations. For three problems, we address questions of modelling and discretization thanks to the observation of the multiplicity of scales, time or space. We propose in the first part a model of approximation of a fluid with a free surface (ocean). The derivation of our multilayer shallow water type model is based on the analysis of the different space scales generally observed in geophysical flows, precisely the 'shallow water' assumption. We obtain an existence and uniqueness result of local in time solution and propose a finite volume scheme and numerical simulations. Next we study a hyperbolic relaxation problem, motivated by the kinetic theory of gaz. Adopting an Asymptotic Preserving strategy of discretization, we build and analyze a numerical scheme. The convergence is proved for any value of the relaxation parameter, as well as the consistency with the equilibrium problem, thanks to error estimates. We present some numerical simulations. The last part deals with a blood flow model in a stented artery. We consider a Stokes problem in a multiscale space domain, that is a macroscopic box (the artery) containing a microscopic roughness (the stent). To avoid expensive simulations when discretizing the whole rough domain, we perform a Chapman-Enskog type expansion of the solution and derive an implicit wall law on the boundary of the smooth domain. Error estimates are shown and numerical simulations are presented

Ce travail s'inscrit dans l'étude mathématique d'écoulements de fluides environnementaux. Nous en abordons deux aspects, à travers deux contextes distincts d'application.En lien avec la simulation des écoulements en milieux poreux, on s'intéresse dans une première partie à la discrétisation d'opérateurs de diffusion anisotropes hétérogènes par des méthodes de volumes finis sur des maillages généraux. Dans le but d'obtenir des solutions approchées qui respectent les bornes physiques des modèles, notre attention se porte sur la conservation du principe du maximum pour les opérateurs elliptiques. Nous présentons des mécanismes généraux permettant de corriger tout schéma volumes finis afin de garantir un principe du maximum discret tout en préservant certaines de ses propriétés principales. On étudie en particulier les propriétés de coercivité et de convergence des schémas corrigés.La deuxième partie est consacrée à la construction de modèles approchés pour la propagation des vagues en eaux peu profondes et sur des topographies irrégulières. A cet effet, nous proposons tout d'abord une adaptation de la démarche d'étude classique à des écoulements bidimensionnels sur des topographies polygonales. Dans un cadre plus général, nous développons ensuite une démarche formelle qui débouche sur des alternatives non locales à quelques modèles classiques (équations de Saint-Venant, équations de Serre, système de Boussinesq). Ces nouveaux modèles contiennent des termes régularisants pour les contributions du fond. / This work investigates two research questions associated with environmental flows and their mathematical modeling.The first part is devoted to the development of finite volume methods for anisotropic and heterogeneous diffusion operators arising in models of porous media flows. To ensure that the approximate solutions lie within physical bounds, we aim at maintaining a discrete analogous of the maximum principle for elliptic operators. Starting from any given cell-centered finite volume scheme, we present a general approach to devise non-linear corrections providing a discrete maximum principle while retaining some main properties of the scheme. In particular, we study the coercivity and convergence properties of the modified schemes.The second part of this work focuses on the derivation of approximate models for shallow water wave propagation over rough topographies. In the particular case of one-dimensional polygonal bottom profiles, we first propose an adaptation of the usual derivation method using complex analysis tools. We then develop a formal approach to account for more general topographies. We propose nonlocal alternatives to some classical models (namely Saint-Venant equations, Serre equations and Boussinesq system). All these alternative models only involve smoothing contributions of the bottom.

Lors de la rentrée atmosphérique nous sommes amenés à modéliser trois phénomènes physiques différents. Tout d'abord, l'écoulement autour du véhicule entrant dans l'atmosphère est hypersonique,il est caractérisé par la présence d'un choc fort et provoque un fort échauffement du véhicule. Nous modélisons l'écoulement par les équations de Navier-Stokes compressibles et l'échauffement du véhicule au moyen de la thermique anisotrope. De plus le véhicule est protégé par un bouclier thermique siège de réactions chimiques que l'on nomme communément ablation.Dans le premier chapitre de cette thèse nous présentons le schéma numérique de diffusion CCLAD (Cell-Centered LAgrangian Diffusion) que nous utilisons pour résoudre la thermique anisotrope. Nous présentons l'extension en trois dimensions de ce schéma ainsi que sa parallélisation.Nous continuons le manuscrit en abordant l'extension de ce schéma à une équation de diffusion tensorielle. Cette équation est obtenue en supprimant les termes convectifs de l'équation de quantité de mouvement des équations de Navier-Stokes. Nous verrons qu'une pénalisation doit être introduite afin de pouvoir inverser la loi constitutive et ainsi appliquer la méthodologie CCLAD. Nous présentons les propriétés numériques du schéma ainsi obtenu et effectuons des validations numériques.Dans le dernier chapitre, nous présentons un schéma numérique de type Volumes Finis permettant de résoudre les équations de Navier-Stokes sur des maillages non-structurés obtenu en réutilisant les deux schémas de diffusion présentés précédemment. / When studying the problem of atmospheric reentry we need to model three different physical phenomenons. First, the ow around the atmospheric reentry vehicle is hypersonic, it is characterized by the presence of a strong shock which leads to a rapid heating of the vehicle. We model the ow using the compressible Navier-Stokes equations and the heating of the vehicle is modeled with the anisotropic heat transfer equation. Furthermore the vehicle is protected by an heat shield, where thermochemical reactions, commonly named ablation, occurs.In the first chapter of this thesis we introduce the numerical diffusion scheme CCLAD (Cell-Centered LAgrangian Diffusion) that we use to solve the anisotropic heat diffusion. We develop its non trivial extension to three-dimensional geometries and present its parallelization. We continue this thesis by the presentation of the extension of this scheme to tensorial diffusion. This equation is obtained by suppressing the convective terms of the momentum equation of the Navier-Stokes equations. We show that we need to introduce a penalization term in order to be able to invert the constitutive law. The invertibility of the constitutive law allows us to apply the CCLAD methodology to this equation straightforwardly. We present the numerical properties of this scheme and show numerical validations.In the last chapter, we present a Finite Volume scheme on unstructured grids that solves the compressible Navier-Stokes equations. This numerical scheme is mainly obtained by gathering the contributions of the two diffusion schemes we developed in the previous chapters.