En simulation dynamique des réservoirs, un des artéfacts les plus gênants pour la prédiction de production est l’effet de l’orientation du maillage. Bien que celui-ci soit « normal » pour tout schéma numérique, il se trouve amplifié par l’instabilité du modèle physique, ce qui a lieu lorsque le contraste de mobilités entre l’eau (fluide poussant, utilisé dans les procédés de récupération secondaires) et l’huile (fluide poussé, contenant les hydrocarbures) dépasse un certain seuil critique. On parle alors d’écoulements à rapport de mobilités défavorable. Connu depuis longtemps, ce problème a fait l'objet de nombreux travaux dans les années 1980 ayant abouti au schéma dit à neuf points. Actuellement implanté dans PumaFlow, logiciel développé et commercialisé par IFPEN, ce schéma fonctionne relativement bien en maillages carrés et dépend d’un paramètre scalaire dont le réglage varie selon les auteurs sur la base de considérations heuristiques. Dans cette thèse, nous proposons une nouvelle démarche méthodologique afin non seulement d’ajuster ce paramètre libre de manière optimale mais aussi de généraliser le schéma aux maillages rectangulaires. La stratégie que nous préconisons repose sur une analyse d’erreur du problème, à partir de laquelle il est possible de définir une notion d’erreur angulaire et de garantir que le comportement du schéma obtenu soit le « moins anisotrope » possible via une minimisation de son écart par rapport à un comportement idéal. Cette procédure de minimisation est ensuite appliquée à deux autres familles de schémas numériques~ : (1) un schéma multidimensionnel proposé par Kozdon, dans lequel le paramètre libre est une fonction~ ; (2) un autre schéma à neuf points faisant intervenir deux paramètres scalaires. C’est ce dernier qui réduit le mieux l’effet de l’orientation lorsque le rapport des pas de maillage s’éloigne de 1. Enfin, une extension de la méthode à des modèles physiques plus complets est envisagée. / In dynamic reservoir simulation, one of the most troublesome artifacts for the prediction of production is the grid orientation effect. Although this normally arises from any numerical scheme, it happens to be amplified by the instability of the physical model, which occurs when the mobility contrast between the water (pushing fluid, used in the processes of secondary recovery) and the oil (pushed fluid, containing the hydrocarbons) exceeds a some critical threshold. We then speak of flows with adverse mobility ratio. This GOE issue has received a lot of attention from the engineers. Numerous works dating back to the 1980s have resulted in the so-called nine-point scheme. Currently implemented in the IFPEN software PumaFlow, this scheme performs relatively well in square meshes and depends on a scalar parameter whose value varies from one author to another, on the grounds of heuristic considerations. In this thesis, we propose a new methodological approach in order not only to optimally adjust this free parameter, but also to extend the scheme to rectangular meshes. The strategy that we advocate is based on an error analysis of the problem, from which it is possible to define a notion of angular error and to guarantee that the behavior of the obtained scheme is the "least anisotropic" possible through a minimization of its deviation from some ideal behavior. This minimization procedure is then applied to two other families of numerical schemes: (1) a multidimensional scheme proposed by Kozdon, in which the free parameter is a function; (2) another nine-point scheme involving two scalar parameters. The latter provides the best results regarding GOE reduction when the ratio of the mesh steps is far away from 1. Finally, an extension of the method to more sophisticated physical models is envisaged.

Bezsíťové metody reprezentují alternativu ke standardním diskretizačním technikám, které pro svůj chod vyžadují síť. V posledních desetiletích bylo vynaloženo mnoho úsilí k ověření konkurenceschopnosti bezsíťových metod v různých inženýrských odvětvích. Diplomová práce je zaměřena na aplikaci vhodné bezsíťové metody ve výpočetní aeroakustice. Stěžejní část této práce se zabývá úlohami šíření zvuku, které lze modelovat pomocí linearizovaných Eulerových rovnic. Obecně lze tyto rovnice zařadit mezi lineární hyperbolické systémy. Pro úlohy aeroakustiky se jako vhodná bezsíťová metoda jeví Finite point method (FPM), která byla úspěšně použita pro řešení úloh dynamiky tekutin. Odvozením této metody a návrhy k dosažení vysoké přesnosti se věnuje další část práce. Úlohy šíření zvuku se známým řešením jsou testovány vlastním softwarem vyvinutým speciálně pro tyto účely.

In this Thesis we explore, both theoretically and numerically, splitting strategies for a hyperbolic system of one-dimensional (1D) blood flow equations with a passive scalar transport equation. Our analysis involves a two-step framework that includes splitting at the level of partial differential equations (PDEs) and numerical methods for discretizing the ensuing problems. This study is inspired by the original flux splitting approach of Toro and Vázquez-Cendón (2012) originally developed for the conservative Euler equations of compressible gas dynamics. In this approach the flux vector in the conservative case, and the system matrix in the non-conservative one, are split into advection and pressure terms: in this way, two systems of partial differential equations are obtained, the advection system and the pressure system. From the mathematical as well as numerical point of view, a basic problem to be solved is the special Cauchy problem called the Riemann problem. This latter provides an analytical solution to evaluate the performance of the numerical methods and, in our approach, it is of primary importance to build the presented numerical schemes. In the first part of the Thesis a detailed theoretical analysis is presented, involving the exact solution of the Riemann problem for the 1D blood flow equations, depicted for a general constant momentum correction coefficient and a tube law that allows to describe both arteries and veins with continuous or discontinuous mechanical and geometrical properties and an advection equation for a passive scalar transport. In literature, this topic has been already studied only for a momentum correction coefficient equal to one, that is related to the prescribed velocity profile and in this case corresponds to a flat one, i.e. an inviscid fluid. In the case of discontinuous properties, only the subsonic regime is considered. In addition we propose a procedure to compute the obtained exact solution and finally we validate it numerically, by comparing exact solutions to those obtained with well-known, numerical schemes on a carefully designed set of test problems. Furthermore, an analogous theoretical analysis and resolution algorithm are presented for the advection system and the pressure system arising from the splitting at the level of PDEs of the complete system of 1D blood flow equations. It is worth noting that the pressure system, in case of veins, presents a loss of genuine non-linearity resulting in the formation of rarefactions, shocks and compound waves, these latter being a composition of rarefactions and shocks. In the second part of the Thesis we present novel finite volume-type, flux splitting-based, numerical schemes for the conservative 1D blood flow equations and splitting-based numerical schemes for the non-conservative 1D blood flow equations that incorporate an advection equation for a passive scalar transport, considering tube laws that allow to model blood flow in arteries and veins and take into account a general constant momentum correction coefficient. A detailed efficiency analysis is performed in order to showcase the advantages of the proposed methodologies in comparison to standard approaches.

Es wird gezeigt, dass ein Satz über die Abbildung spektraler Lücken, welcher exponentielle Dichotomie charakterisiert, für eine allgemeine Klasse (SH) von semilinearen hyperbolischen Systemen von partiellen Differentialgleichungen in einem Banach-Raum X von stetigen Funktionen gilt. Dies beantwortet ein Schlüsselproblem für die Existenz und Glattheit invarianter Mannigfaltigkeiten semilinearer hyperbolischer Systeme. Unter natürlichen Annahmen an die Nichtlinearitäten wird gezeigt, dass schwache Lösungen von (SH) einen glatten Halbfluß im Raum X bilden. Für Linearisierungen werden hochfrequente Abschätzungen für Spektren sowie Resolventen unter Verwendung von reduzierten (block)diagonal Systemen hergestellt. Darauf aufbauend wird der Abbildungssatz für spektrale Lücken im kleinen Raum X bewiesen: Eine offene spektrale Lücke des Generators wird exponentiell auf eine offene spektrale Lücke der Halbruppe abgebildet und umgekehrt. Es folgt, dass ein Phänomen wie im Gegenbeispiel von Renardy nicht auftreten kann. Unter Verwendung der allgemeinen Theorie implizieren die Ergebnisse die Existenz von glatten Zentrumsmannigfaltigkeiten für (SH). Die Ergebnisse werden auf traveling wave Modelle für die Dynamik von Halbleiter Lasern angewandt. Für diese werden Moden Approximationen (Systeme von gewöhnlichen Differentialgleichungen, welche die Dynamik auf gewissen Zentrumsmannigfaltigkeiten approximativ beschreiben) hergeleitet und gerechtfertigt, die generische Bifurkation von modulierten Wellen aus rotierenden Wellen wird gezeigt. Globale Existenz und glatte Abhängigkeit von nichtautonomen traveling wave Modellen werden betrachtet, außerdem werden Moden Approximationen für solche nichtautonomen Modelle rigoros hergeleitet. Insbesondere arbeitet die Theorie für die Stabilitäts- und Bifurkationsanalyse von Turing Modellen mit korellierter Zufallsbewegung. Ferner beinhaltet die Klasse (SH) neutrale und retardierte funktionale Differentialgleichungen. / A spectral gap mapping theorem, which characterizes exponential dichotomy, is proven for a general class of semilinear hyperbolic systems of PDEs in a Banach space X of continuous functions. This resolves a key problem on existence and smoothness of invariant manifolds for semilinear hyperbolic systems. It is shown that weak solutions to (SH) form a smooth semiflow in X under natural conditions on the nonlinearities. For linearizations high frequency estimates of spectra and resolvents in terms of reduced diagonal and blockdiagonal systems are given. Using these estimates a spectral gap mapping theorem in the small Banach space X is proven: An open spectral gap of the generator is mapped exponentially to an open spectral gap of the semigroup and vice versa. Hence, a phenomenon like in Renardy''s counterexample cannot appear for linearizations of (SH). By the general theory the results imply existence of smooth center manifolds for (SH). Moreoever, the results are applied to traveling wave models of semiconductor laser dynamics. For such models mode approximations (ODE systems which approximately describe the dynamics on center manifolds) are derived and justified, and generic bifurcations of modulated waves from rotating waves are shown. Global existence and smooth dependence of nonautonomous traveling wave models with more general solutions, which possess jumps, are considered, and mode approximations are derived for such nonautonomous models. In particular the theory applies to stability and bifurcation analysis for Turing models with correlated random walk. Moreover, the class (SH) includes neutral and retarded functional differential equations.

La Commande prédictive ou Commande Optimale à Horizon Glissant (COHG) devient de plus en plus populaire dans de nombreuses applications pratiques en raison de ses avantages importants tels que la stabilisation et la prise en compte des contraintes. Elle a été bien étudiée pour des systèmes en dimension finie même dans le cas non linéaire. Cependant, son extension aux systèmes en dimension infinie n'a pas retenu beaucoup d'attention de la part des chercheurs. Ce travail de thèse apporte des contributions à l'application de cette approche aux systèmes de lois de conservation. Nous présentons tout d'abord une preuve de stabilité complète de la COHG pour certaines classes de systèmes en dimension infinie. Ce résultat est ensuite utilisé pour les systèmes hyperboliques 2x2 commandés aux frontières et appliqué à un problème de contrôle de canal d'irrigation. Nous proposons aussi l'extension de cette stratégie au cas de réseaux de systèmes hyperboliques 2x2 en cascade avec une application à un ensemble de canaux d'irrigation connectés. Nous étudions également les avantages de la COHG dans le contexte des systèmes non linéaires et semi-linéaires notamment vis-à-vis des chocs. Toutes les analyses théoriques sont validées par simulation afin d'illustrer l'efficacité de l'approche proposée. / The predictive control or Receding Horizon Optimal Control (RHOC) is becoming increasingly popular in many practical applications due to its significant advantages such as the stabilization and constraints handling. It has been well studied for finite dimensional systems even in the nonlinear case. However, its extension to infinite dimensional systems has not received much attention from researchers. This thesis proposes contributions on the application of this approach to systems of conservation laws. We present a complete proof of stability of RHOC for some classes of infinite dimensional systems. This result is then used for 2x2 hyperbolic systems with boundary control, and applied to an irrigation canal. We also propose the extension of this strategy to networks of cascaded 2x2 hyperbolic systems with an application to a set of connected irrigation canals. Furthermore, we study the benefits of RHOC in the context of nonlinear and semi-linear systems in particular with respect to the problem of shocks. All theoretical analyzes are validated by simulation in order to illustrate the effectiveness of the proposed approach.

Les méthodes hybrides Moments/PDF sont bien adaptées pour la description des écoulements diphasiques turbulents, polydispersés, hors équilibre thermodynamique. Ces méthodes permettent d'avoir une description assez fine de la polydispersion, de la convection et des termes sources non-linéaires. Cependant, les approximations issues de telles simulations sont bruitées ce qui, dans certaines situations, occasionne un biais. L'approche alternative étudiée dans ce travail consiste à coupler une description Eulerienne des moments avec une description stochastique Lagrangienne à l'intérieur de la phase dispersée, permettant de réduire l'erreur statistique et d'éliminer le biais. La mise en oeuvre de cette méthode nécessite le développement de schémas numériques robustes. Les approches proposées sont basées sur l'utilisation simultanée des techniques de relaxation et de décentrement, et permettent d'obtenir des approximations stables des solutions instationnaires du système d'équations aux dérivées partielles, avec des données peu régulières provenant du couplage avec le modèle stochastique. Une comparaison des résultats de la méthode hybride Moments-Moments/PDF avec ceux issus de la méthode hybride "classique'' est présentée en termes d'analyse des erreurs numériques sur un cas de jet co-courant gaz-particules. / Hybrid Moments/PDF methods have shown to be well suitable for the description of polydispersed turbulent two-phase flows in non-equilibrium which are encountered in some industrial situations involving chemical reactions, combustion or sprays. hey allow to obtain a fine enough physical description of the polydispersity, non-linear source terms and convection phenomena. However, their approximations are noised with the statistical error, which in several situations may be a source of a bias. An alternative hybrid Moments-Moments/PDF approach examined in this work consists in coupling the Moments and the PDF descriptions, within the description of the dispersed phase itself. This hybrid method could reduce the statistical error and remove the bias. However, such a coupling is not straightforward in practice and requires the development of accurate and stable numerical schemes. The approaches introduced in this work rely on the combined use of the upwinding and relaxation-type techniques. They allow to obtain stable unsteady approximations for a system of partial differential equations containing non-smooth external data which are provided by the PDF part of the model. A comparison of the results obtained using the present method with those of the ``classical'' hybrid approach is presented in terms of the numerical errors for a case of a co-current gas-particle wall jet.

On étudie dans ce travail des systèmes de lois de conservation hyperboliques. La première partie étudie le temps d'existence des solutions régulières et régulières par morceaux de la dynamique des fluides compressibles. Après avoir présenté l'état de l'art en matière de solutions régulières, on montre une extension d'un théorème de Grassin à des gaz de Van der Waals. On étudie ensuite les solutions ondes de chocs : on poursuit l'approche de T. T. Li pour estimer leur temps d'existence dans le cas isentropique à symétrie sphérique, et l'approche de Whitham afin d'obtenir une équation approchée vérifiée par la surface de discontinuité. Dans une deuxième partie, motivée par la modélisation d'un rond-point en trafic routier, on étudie une extension multi-classe du modèle macroscopique de Lighthill-Whitham-Richards sur une route infinie avec des jonctions. On différencie les véhicules selon leur origine et leur destination et on introduit des conditions aux bords adaptées au niveau des jonctions. On obtient existence et unicité d'une solution au problème de Riemann pour ce modèle. Des simulations numériques attestent que les solutions obtenues existent en temps long. On aborde enfin le problème de Cauchy par la méthode de front tracking. La dernière partie concerne les lois de conservation scalaires. La première question abordée est le contrôle de la variation totale de la solution et la stabilité des solutions faibles entropiques par rapport au flux et à la source. Ce résultat nous permet d'étudier des équations avec flux non-local. Une fois établi leur caractère bien posé, on montre la Gâteaux-différentiabilité du semi-groupe obtenu par rapport aux conditions initiales. / In this work, we study hyperbolic systems of balance laws. The first part is devoted to compressible fluid dynamics, and particularly to the lifespan of smooth or piecewise smooth solutions. After presenting the state of art, we show an extension to more general gases of a theorem by Grassin.We also study shock waves solutions: first, we extend T. T. Li's approach to estimate the time of existence in the isentropic spherical case; second, we develop Whitham's ideas to obtain an approximated equation satisfied by the discontinuity surface. In the second part, we set up a new model for a roundabout. This leads us to study a multi-class extension of the macroscopic Lighthill-Whitham-Richards' model. We study the traffic on an infinite road, with some points of junction. We distinguish vehicles according to their origin and destination and add some boundary conditions at the junctions. We obtain existence and uniqueness of a weak entropy solution for the Riemann problem. As a complement, we provide numerical simulations that exhibit solutions with a long time of existence. Finally, the Cauchy problem is tackled by the front tracking method. In the last part, we are interested in scalar hyperbolic balance laws. The first question addressed is the control of the total variation and the stability of entropy solutions with respect to flow and source. With this result, we can study equations with non-local flow, which do not fit into the framework of classical theorems. We show here that these kinds of equations are well posed and we show the Gâteaux-differentiability with respect to initial conditions, which is important to characterize maxima or minima of a given cost functional.

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

Dans divers domaines de la physique, certains phénomènes sont modélisés par des systèmes hyperboliques non-conservatifs. En particulier, dans le domaine de la physique des plasmas, dont l'un des champs d'application majeur est la Fusion par Confinement Inertiel, le système d'Euler bi-température, modélisant les phénomènes de transport de particules chargées, en est un exemple. La difficulté de l'étude de ces systèmes réside dans la présence de termes non-conservatifs, qui empêchent la définition classique des solutions faibles. Pour parvenir à une définition de ce type de solutions, on a recours à l'emploi de systèmes cinétiques sous-jacents. Dans ce manuscrit, on s'intéresse à l'étude numérique de ces systèmes cinétiques pour la résolution du système d'Euler bi-température.Ce manuscrit se divise en deux parties. La première partie contient l'étude numérique du système d'Euler bi-température. Dans un premier chapitre, on résout numériquement les équations en dimension 1 d'espace par le biais d'un système sous-jacent issu de la physique des plasmas: le système de Vlasov-BGK-Ampère. On présente une méthode numérique préservant l'asymptotique pour ce système sous-jacent et on montre, par des simulations numériques, que le schéma limite obtenu donne des résultats consistants avec Euler bi-température. Dans un second chapitre, on résout le même modèle en dimension 2 d'espace par un système sous-jacent de type BGK discret. On démontre une inégalité d'entropie pour les solutions issues du modèle sous-jacent, ainsi qu'une inégalité discrète de dissipation d'entropie pour le schéma.Dans la deuxième partie de ce manuscrit, on s'intéresse au développement de méthodes numériques pour quelques modèles cinétiques. On considère ici le cas des écoulements raréfiés de mélanges de gaz, dans l'optique d'une application aux cas des plasmas. Premièrement, on présente un schéma cinétique adaptatif et dynamique en vitesse pour les gaz inertes. Par l'emploi de lois de conservation discrètes, la solution est approchée sur un ensemble de vitesses discrètes local et dynamique. Dans un second temps, on propose une extension de cette méthode visant à améliorer les performances de celle-ci. Puis, ces deux versions de la méthode sont comparées à la méthode classique sur grille fixe uniforme sur une série de cas tests.Enfin, dans le dernier chapitre, on propose une méthode numérique pour la résolution d'une extension de ces équations, prenant en compte la présence de réactions chimiques au sein du mélange. Le contexte considéré est celui des réactions chimiques bi-moléculaires réversibles lentes. La méthode proposée, de type implicite-explicite, est linéaire, stable et conservative. / In various domains of physics, several phenomena can be modeled via the use of nonconservative hyperbolic systems. In particular, in plasma physics, in the process of developping and understanding the phenomena leading to Inertial Confinement Fusion, the bi-temperature Euler sytem can be used to model particle transport phenomena in a plasma. The difficulty of the mathematical study of such systems dwells in the presence of so-called non-conservative products, which prevent the classical definition of weak solutions via distribution theory. To attempt to define these quantities, it is useful to supplement the hyperbolic system with an underlying kinetic model. In this work, the objective is the numerical study of such kinetic systems in order to solve the bi-temperature Euler system.This manuscript is split in two parts. The first one contains the study of the bi-temperature Euler system. In the first chapter, this system in dimension 1 is solved by the use of an underlying kinetic model sprung from plasma physics: the Vlasov-BGK-Ampère system. An asymptotic-preserving numerical method is introduced, and it is shown that the scheme obtained in the limit is consistant with a scheme for teh bi-temperature Euler system. In the following chapter, the same hyperbolic model in dimension 2 is studied, this time via a discrete-BGK type underlying model. An entropy inequality is proved for solutions coming from the kinetic model, as well as a discrete entropy dissipation inequality.In the second part of the manuscript, we are interested in the development of numerical schemes for gas mixture rarefied flows. Firstly, an adaptive kinetic scheme is introduced for inert gas mixtures. By the use of discrete conservation laws, the solution is approximated on a set of discrete velocities that depends on space, time and species. Secondly, an extension of the method is proposed in order to improve the efficiency of the first method. Finally, the two methods are compared to the classical fixed grid method on a series of test cases.In the last chapter, a numerical method is proposed for rarefied flows of reacting mixtures. The setting considered is the case of slow bimolecular reversible chemical reactions. The method introduced is an explicit-implicit treatment of the relaxation operator, which is shown to be stable, linear and conservative.

Dans cette thèse, nous nous intéressons à l'étude d'un modèle de transport de sédiments en nous plaçant sous deux angles d'approche différents. L'un concerne la modélisation numérique du problème et propose une méthode de résolution numérique basée sur un solveur de Riemann approché pour le système de Saint-Venant-Exner qui reste un des systèmes les plus répandus pour traiter le transport sédimentaire par charriage. Ce dernier repose sur un couplage du modèle hydraulique de Saint-Venant et du modèle morphodynamique d'Exner. Le point essentiel de la méthode proposée se situe au niveau du traitement du couplage de ce système. En effet, il existe deux stratégies; la première consiste à découpler la résolution de la partie fluide de la partie solide et les faire interagir à des instants donnés alors que la seconde considère une résolution couplée du système en mettant à jour conjointement les grandeurs hydrauliques et solides aux mêmes instants. Se posera alors la question du choix de la stratégie de résolution pour laquelle nous apporterons des éléments de réponses en comparant les deux approches. L'autre se concentre sur la mise en place d'une méthodologie pour l'étude des incertitudes liées au même modèle. Pour ce faire, nous proposons une formulation stochastique du système de Saint-Venant-Exner et nous cherchons à caractériser la variabilité des sorties par rapport aux paramètres d'entrée naturellement aléatoires. Cette première étude révélera la nécessité de revenir à un système de Saint-Venant avec un fond bruité pour étudier la sensibilité des grandeurs hydrauliques par rapport aux perturbations topographiques. / In this thesis, we are interested on the study of a sediment transport model through two different approaches. One of them concerns the numerical modelling of the problem and proposes a numerical problem-solving method based on an approximate Riemann solver for the Saint-Venant-Exner system which is one of the most common model to deal with sedimentary bed-load transport. This last one is based on a coupling between the hydraulic model of Saint-Venant and the morphodynamic model of Exner. The key point of the proposed method is the treatment of the coupling issue. Indeed, there exists two strategies; the first one consists on decoupling the resolution of the fluid part from the solid part and making them interact at fixed times whereas the second one considers a coupled approach to solve the system by jointly updating the hydraulic and solid quantities at same times. We then raise the issue of the choice of the strategy for which we suggest answers by comparing both approaches. The other one focuses on the development of a methodology to study the uncertainties related to the model previously mentioned. To this end, we propose a stochastic formulation of the Saint-Venant-Exner system and we look for characterizing the variabilities of the outputs in relation to the naturally random input parameters. This first study reveals the need for a return to the Saint-Venant system with a perturbed bed to understand the sensitivity of the hydraulic quantities on the topographical perturbations.